Wie bereits in einem Kommentar von John Custer ausgeführt, war Dirac in
den zugrunde liegenden physikalischen Gesetzen, die für die mathematische Theorie eines großen Teils der Physik und der gesamten Chemie erforderlich sind, zu 100% korrekt sind somit vollständig bekannt, und die Schwierigkeit besteht nur darin, dass die genaue Anwendung dieser Gesetze zu Gleichungen führt, die viel zu kompliziert sind, um löslich zu sein. Es wird daher wünschenswert, ungefähre praktische Methoden zur Anwendung der Quantenmechanik zu entwickeln, die zu einer Erklärung der Hauptmerkmale komplexer Atomsysteme ohne zu viel Berechnung führen können.
Für Lichtsysteme, typischerweise als H-Kr gedacht, reicht eine nicht-relativistische Quantenmechanik, dh die Schrödinger-Gleichung, aus, um die Chemie zu beschreiben; Für schwerere Kerne müssen Sie die relativistische Quantenmechanik verwenden, d. h. die Dirac-Gleichung, die etwas komplizierter ist. In vielen Fällen können wir die Born-Oppenheimer-Methode aufrufen und annehmen, dass sich die Kerne im Momentanfeld der Elektronen bewegen; Jetzt müssen wir nur noch das elektronische Problem lösen.
Wir wissen, dass die genaue Lösung des elektronischen Problems in diesem Fall mit der von Ihnen beschriebenen Methode der vollständigen Konfigurationsinteraktion (FCI) (Exact Diagonalization) erreicht werden kann Die elektronische Vielteilchenwellenfunktion durch eine gewichtete Summe elektronischer Konfigurationen, auch Slater-Determinanten genannt, als $ | \ Psi \ rangle = \ sum_k c_k | \ Phi_k \ rangle $ span>. Diese elektronischen Konfigurationen werden erstellt, indem $ N $ span> -Elektronen in $ K $ span> Einzelteilchenzustände a.k.a. Orbitale verteilt werden. Damit die Methode genau ist, $ K \ gg N $ span> und um genau die Lösung zu erhalten, benötigen Sie $ K \ zu \ infty $ span>.
Um nun den Grundzustand (sowie alle angeregten Zustände) zu ermitteln, müssen Sie nur den Vielelektronen-Hamilton-Operator anhand der Elektronenkonfigurationen diagonalisieren. Das Problem ist jedoch, dass die Anzahl der Elektronenkonfigurationen extrem schnell zunimmt.
Wenn wir davon ausgehen, dass es sich um einen Spin-Singulett-Zustand handelt, haben Sie $ N / 2 $ span> Spin-up- und $ N / 2 $ span> Spin-down-Elektronen. Für jeden Spin gibt es $ {K \ wähle N / 2} = \ frac {K!} {\ Frac N 2! (K- \ frac N 2)!} $ Span> Möglichkeiten, die Orbitale zu füllen. Die Gesamtzahl der Elektronenkonfigurationen für $ N $ span> -Elektronen in $ K $ span> -Orbitalen, typischerweise bezeichnet als ( $ N $ span> e, $ K $ span> o) ist dann $ {K \ wähle N / 2} ^ 2 = \ left [\ frac {K!} {\ Frac N 2! (K- \ frac N 2)!} \ Right] ^ 2 $ span>.
Selbst für den Fall einer sehr kleinen Anzahl von Orbitalen $ K = N $ span>, die Anzahl der Konfigurationen wird schnell riesig. (8e, 8o) hat 4900 Konfigurationen, (10e, 10o) hat 63 504, (12e, 12o) hat 853 776, (14e, 14o) hat 11 778 624, (16e, 16o) hat 165 636 900, (18e, 18o) hat 2 363 904 400, (20e, 20o) hat 34 134 779 536 und (22e, 22o) hat 497 634 306 624.
Obwohl Sie das Problem (8e, 8o) mit dichter Matrixalgebra auf modernen Computern immer noch lösen können, sehen Sie, dass Sie sehr schnell sehr klug werden müssen, um die Matrix zu diagonalisieren. Da der Hamilton-Operator ein Zwei-Teilchen-Operator ist, ist er aufgrund der Elektronenkonfigurationen äußerst spärlich: Wenn sich zwei Konfigurationen durch die Besetzung von mehr als zwei Orbitalen unterscheiden, ist das Hamilton-Matrixelement nach den Regeln von Slater und Condon Null. Darüber hinaus möchten Sie bei großen Problemgrößen auch das Speichern der Matrix vermeiden, weshalb Sie eine iterative Methode verwenden möchten. (Die berühmte Davidson-Methode zur iterativen Diagonalisierung wurde tatsächlich genau für diesen Zweck entwickelt!)
Mit intelligenten Algorithmen waren seit Anfang der neunziger Jahre Berechnungen der Milliardenkonfiguration möglich, d. H. Das (18e, 18o) -Problem, siehe z. Chem. Phys. Lette. 169, 463 (1990). Trotz des enormen Anstiegs der Rechenleistung in den letzten 30 Jahren hat sich die Barriere kaum bewegt: Soweit mir bekannt ist, ist das größte FCI-Problem, das gelöst wurde, die (22e, 22o) -Berechnung in J. Chem. Phys. 147, 184111 (2017).
Hier ist zu beachten, dass selbst die Berechnung (22e, 22o) nicht groß genug ist, um ein einzelnes Atom genau zu lösen: Sie benötigen a viel mehr Orbitale, um mit Experimenten eine quantitative Genauigkeit zu erreichen. Obwohl ein hoch liegendes Orbital nur sehr wenig zur Korrelationsenergie beiträgt, gibt es VIELE davon.
Genau wie Dirac schrieb, sind Annäherungen erforderlich. Dichtefunktionale Approximationen sind in Anwendungen äußerst beliebt, aber alles andere als exakt. Andererseits verwenden hochgenaue Studien häufig die Coupled-Cluster-Methode, bei der es sich um eine Reparametrisierung der FCI-Methode handelt. es würde jedoch auch eine exponentielle Skalierung zeigen, wenn es nicht abgeschnitten würde - d. h. angenähert