Quel est l'état actuel de la technique dans la résolution des PDE paraboliques de dimension supérieure (équation de Schrödinger multi-électrons)

11

Quel est l'état actuel de la technique pour résoudre les PDE paraboliques de dimension supérieure (3-10) dans le domaine complexe avec des pôles simples (de la forme ) et absorber les conditions aux limites?1|r1-r2|

Plus précisément, je suis intéressé par la résolution de l'équation de Schrödinger multi-électrons:

(jejje[-je22m-ZjeZj|rje-rj|+V(rje,t)])ψ=-jetψ

Pour une molécule diatomique avec plus d'un électron.

Andrew Spott
la source

Réponses:

6

Les solutions pour l'équation sont dans Si le nombre d'électrons est suffisamment petit, vous pouvez simplement utiliser n'importe quelle méthode traditionnelle . Comme une méthode de discrétisation de domaine (différence finie, élément fini, élément frontière) ou une méthode pseudospectrale. Puisque résoudre cette équation n'est pas plus difficile que de résoudre une équation d'onde multidimensionnelle.

ψC3M×R+.

Dans le cas de systèmes plus gros, une astuce est nécessaire pour obtenir la solution. Nous remplaçons l'interaction électron-électron pour l'interaction d'un électron avec un nuage d'électrons (une approximation de champ moyen du reste d'entre eux), puis résolvons de manière auto-cohérente (en raison de la non-linéarité qui provient du champ moyen terme). Cela se fait dans Hartree-Fock et la théorie fonctionnelle de la densité (DFT). Où l'équation différentielle d'origine est transformée en une formulation variationnelle.

La DFT est la méthode la plus courante de nos jours, et l'avantage est que toutes les équations sont formulées en termes de densité électronique et non en termes d'équations d'onde. Ainsi, les équations se situent dans un espace tridimensionnel. Un livre qui décrit ces deux méthodes est

  • Thijssen, Jos. Physique numérique. Cambridge University Press, 2007. Lien Amazon .
nicoguaro
la source
3

Vous souhaitez résoudre pour 3 à 10 systèmes de particules (3D par particule)? Pour autant que je sache, les théories de champ moyen ne fonctionnent pas particulièrement bien pour si peu de particules, mais il semble qu'il y ait eu des travaux DFT sur les molécules diatomiques.

Est-ce un système où Born-Oppenheimer est valable? Si tel est le cas, je pourrais être enclin à étendre la fonction d'onde électronique en utilisant une combinaison linéaire de déterminants de Slater, éventuellement en utilisant une grille clairsemée ou des grilles spectrales clairsemées. Cet article pourrait peut-être aider .

Une autre option consiste à essayer d'utiliser une approche à liaison étroite, bien que le fait que vous ayez mentionné l'absorption des conditions aux limites suggère que vous pensez peut-être à des problèmes impliquant l'ionisation / la dissociation. La TB serait surtout utile si vous tentiez d'approcher des états de bas niveau.

Peut-être quelque chose comme la méthode Hartree-Fock multi-configurationnelle en fonction du temps pourrait fonctionner ici MCTDHF .

Enfin, vous pouvez regarder les méthodes quantiques de Monte Carlo. Ce sont les méthodes par lesquelles des modèles fonctionnels d'échange et de corrélation pour des atomes uniques sont obtenus pour effectuer des calculs DFT. Il semble qu'il existe des extensions poly-atomiques. (Je n'ai plus de privilèges de lien).

Greg
la source
3-10 dimensions, pas des particules: spécifiquement 1 à 3 électrons, 2 noyaux (1d pour les noyaux, 6d pour les particules), sans approximation de Born-Oppenheimer. Et je fais des trucs de type ionisation.
Andrew Spott
1

M3MNNN3MMN=dixdix9

De cette considération découle qu'il n'est pas possible de considérer le problème avec tous les électrons en même temps - vous devez vous limiter à un ou deux électrons à la fois. Cela vous conduit naturellement à des méthodes telles que la méthode Hartree Fock qui itère sur les électrons tout en gardant le reste du système fixe.

Je ne connais pas assez bien le domaine, mais imaginez qu'il existe un certain nombre d'articles de synthèse très cités et bien écrits sur le sujet.

Wolfgang Bangerth
la source
dix9
Eh bien, les systèmes fermioniques ont pas mal de (anti) symétries en raison du principe de Pauli que vous pouvez exploiter pour réduire considérablement le nombre de degrés de liberté (au lieu de l'hypercube à 3 dimensions, vous n'avez qu'à considérer le simplexe correspondant, dont le cube contient (3M)! copies). Vous n'avez donc besoin que de fonctions de base binom (N, 3M) - toujours exponentielles, mais de plus en plus lentes. Cela pourrait mettre l'extrémité inférieure de la gamme à la portée d'un poste de travail costaud.
Christian Clason
Pour un système à 3 électrons, peut-être. Mais vous ne pourrez toujours rien faire au-delà de cela. Cela ne laisse pas un grand nombre de molécules :-)
Wolfgang Bangerth
Mais la question ne demandait que 3 à 10 variables :) (Mais votre point est valable: pour tout ce qui a plus d'un petit nombre d'électrons, vous devez considérer des modèles de champ approximatifs tels que DFT; mon point était qu'entre "peut être résolus avec des approches standard "et" ne peuvent être résolus qu'environ ", il y a une gamme non triviale de problèmes qui" peuvent (seulement) être résolus en utilisant des approches non standard ".)
Christian Clason