Comme déjà souligné dans un commentaire de John Custer, Dirac était 100% correct dans
Les lois physiques sous-jacentes nécessaires à la théorie mathématique d'une grande partie de la physique et de l'ensemble de la chimie sont donc parfaitement connues, et la difficulté est seulement que l'application exacte de ces lois conduit à des équations beaucoup trop compliquées pour être solubles. Il devient donc souhaitable que des méthodes pratiques approximatives d'application de la mécanique quantique soient développées, ce qui peut conduire à une explication des principales caractéristiques de systèmes atomiques complexes sans trop de calculs.
Pour les systèmes légers, généralement considérée comme H-Kr, la mécanique quantique non relativiste, c'est-à-dire que l'équation de Schrödinger est suffisante pour décrire la chimie; pour les noyaux plus lourds, vous devez utiliser la mécanique quantique relativiste, c'est-à-dire l'équation de Dirac qui est un peu plus compliquée. Dans de nombreux cas, nous pouvons invoquer la méthode de Born-Oppenheimer et supposer que les noyaux se déplacent dans le champ instantané des électrons; maintenant tout ce que nous devons résoudre est le problème électronique.
Nous savons que la solution exacte au problème électronique dans ce cas est réalisable avec l'interaction de configuration complète (FCI) aka méthode de diagonalisation exacte, dans laquelle vous décrivez la fonction d'onde électronique à plusieurs particules par une somme pondérée de configurations électroniques aka les déterminants de Slater comme $ | \ Psi \ rangle = \ sum_k c_k | \ Phi_k \ rangle $ . Ces configurations électroniques sont construites en distribuant des électrons $ N $ dans des états de particules uniques $ K $ , alias orbitales. Pour que la méthode soit précise, $ K \ gg N $ , et pour obtenir la solution exacte dont vous avez besoin $ K \ à \ infty $ .
Maintenant, pour trouver l'état fondamental (ainsi que tous les états excités), il vous suffit de diagonaliser l'hamiltonien à plusieurs électrons sur la base des configurations électroniques. Mais, le problème est que le nombre de configurations d'électrons augmente extrêmement rapidement.
Si nous supposons que nous examinons un état singulet de spin, alors vous avez $ N / 2 $ spin-up et $ N / 2 $ électrons spin-down. Pour chaque spin, il y a $ {K \ choose N / 2} = \ frac {K!} {\ Frac N 2! (K- \ frac N 2)!} $ façons de peupler les orbitales. Le nombre total de configurations d'électrons pour les électrons $ N $ dans les orbitales $ K $ , généralement désignés par ( $ N $ e, $ K $ o), est alors $ {K \ choose N / 2} ^ 2 = \ left [\ frac {K!} {\ Frac N 2! (K- \ frac N 2)!} \ Right] ^ 2 $ .
Même pour le cas d'un très petit nombre d'orbitales, $ K = N $ , le nombre de configurations devient rapidement énorme. (8e, 8o) a 4900 configurations, (10e, 10o) a 63504, (12e, 12o) a 853776, (14e, 14o) a 11778624, (16e, 16o) a 165 636900, (18e, 18o) a 2 363 904 400, (20e, 20o) a 34 134 779 536, et (22e, 22o) a 497 634 306 624.
Bien que vous puissiez toujours résoudre le problème (8e, 8o) avec l'algèbre matricielle dense sur les ordinateurs modernes, vous voyez que très rapidement, vous devez devenir très intelligent dans la façon de diagonaliser la matrice. Parce que l'hamiltonien est un opérateur à deux particules, il est extrêmement rare dans la base des configurations électroniques: si deux configurations diffèrent par l'occupation de plus de deux orbitales, l'élément de la matrice hamiltonien est nul selon les règles de Slater et Condon. De plus, pour les problèmes de grande taille, vous souhaitez également éviter de stocker la matrice, c'est pourquoi vous souhaitez utiliser une méthode itérative. (La célèbre méthode Davidson pour la diagonalisation itérative a en fait été développée exactement dans ce but!)
Avec des algorithmes intelligents, les calculs de milliards de configuration, c'est-à-dire le problème (18e, 18o), sont possibles depuis le début des années 1990, voir par ex. Chem. Phys. Lett. 169, 463 (1990). Cependant, malgré l'énorme augmentation de la puissance de calcul au cours des 30 dernières années, la barrière a à peine bougé: à ma connaissance, le plus gros problème de FCI qui ait été résolu est le calcul (22e, 22o) dans J. Chem. Phys. 147, 184111 (2017).
La chose à noter ici est que même le calcul (22e, 22o) n'est pas assez grand pour résoudre un seul atome de manière exacte: vous avez besoin d'un beaucoup plus d'orbitales pour obtenir une précision quantitative avec l'expérience. Bien qu'une orbitale haute ne contribue que très peu à l'énergie de corrélation, il y en a BEAUCOUP.
Exactement comme l'a écrit Dirac, des approximations sont nécessaires. Les approximations fonctionnelles de densité sont extrêmement populaires dans les applications, mais elles sont loin d'être exactes. D'autre part, les études de haute précision utilisent souvent la méthode des clusters couplés, qui est une reparamétrisation de la méthode FCI; cependant, il présenterait également une mise à l'échelle exponentielle s'il n'était pas tronqué - c'est-à-dire approximé.