¿Cuál es el estado actual de la técnica en la resolución de PDE parabólicas de dimensiones superiores (ecuación de Schrödinger de múltiples electrones)

¿Cuál es el estado actual de la técnica para resolver PDE parabólicas de dimensiones superiores (3-10) en el dominio complejo con polos simples (de la forma ) y absorbe las condiciones de contorno?$\frac{1}{|\vec{r}_1 - \vec{r}_2|}$

Específicamente, estoy interesado en resolver la ecuación de Schrödinger de múltiples electrones:

$\left( \sum_i \sum_{j\neq i}\left[ -\frac{\nabla_i^2}{2 m} - \frac{Z_i Z_j}{|\vec{r}_i - \vec{r}_j|} + V(\vec{r}_i, t) \right]\right)\psi = -i\partial_t \psi$

Para una molécula diatómica con más de 1 electrón.

Las soluciones para la ecuación están en Si la cantidad de electrones es losuficientemente pequeña, puede usar cualquiermétodotradicional. Como un método de discretización de dominio (Diferencia finita, Elemento finito, Elemento límite) o un método pseudoespectral. Dado que resolver esta ecuación no es más difícil que resolver una ecuación de onda multidimensional.

En el caso de sistemas más grandes, es necesario algún truco para obtener la solución. Reemplazamos la interacción electrón-electrón por la interacción de un electrón con una nube de electrones (una aproximación de campo medio del resto de ellos), y luego resolvemos de manera coherente (debido a la no linealidad que proviene del campo medio término). Esto se hace en Hartree-Fock y la teoría funcional de la densidad (DFT). Donde la ecuación diferencial original se transforma en una formulación variacional.

DFT es el método más común hoy en día, y la ventaja es que todas las ecuaciones están formuladas en términos de densidad de electrones y no en términos de ecuaciones de onda. Entonces, las ecuaciones se encuentran en un espacio tridimensional. Un libro que describe ambos métodos es

• Thijssen, Jos. Física computacional. Cambridge University Press, 2007. Enlace de Amazon .

¿Desea resolver de 3 a 10 sistemas de partículas (3D por partícula)? Hasta donde yo sé, las teorías de campo medio no funcionan especialmente bien para tan pocas partículas, pero parece que ha habido trabajo de DFT en moléculas diatómicas.

¿Es este un sistema en el que Born-Oppenheimer es válido? Si es así, podría inclinarme a expandir la función de onda electrónica utilizando una combinación lineal de determinantes Slater, posiblemente utilizando una cuadrícula dispersa o cuadrículas espectrales dispersas. Este artículo quizás podría ayudar .

Otra opción es intentar utilizar un enfoque de unión estrecha, aunque el hecho de que haya mencionado la absorción de condiciones límite sugiere que puede estar pensando en problemas relacionados con la ionización / disociación. La TB sería principalmente útil si intentara aproximarse a estados de bajo nivel.

Posiblemente, algo como el método Hartree-Fock de configuración múltiple dependiente del tiempo podría funcionar aquí MCTDHF .

Finalmente, podrías mirar los métodos cuánticos de Monte Carlo. Estos son los métodos por los cuales se obtienen modelos funcionales de intercambio y correlación para átomos individuales para hacer cálculos DFT. Parece que hay extensiones poliatómicas. (Estoy fuera de los privilegios de enlace).

$M$$3M$$N$$N$$N^{3M}$$M$$N=10$$10^9$

De esta consideración se deduce que no es posible considerar el problema con todos los electrones al mismo tiempo; debe restringirse a uno o dos electrones a la vez. Naturalmente, esto lo lleva a métodos como el método Hartree Fock que itera sobre los electrones mientras mantiene el resto del sistema fijo.

No conozco el campo lo suficientemente bien, pero imagino que hay una serie de artículos de revisión muy bien citados y bien escritos sobre el tema.