Improving Perturbation Theory with the Sum-of-squares II: Large Density-Density Terms

Technische Zusammenfassung: Verbesserung der Störungstheorie mit der Summe-der-Quadrate II: Große Dichte-Dichte-Terme

Problemstellung
Die Arbeit adressiert eine signifikante Einschränkung bei der Anwendung selbstkonsistenter Summe-der-Quadrate-Zerlegungen (Sum-of-Squares, SOS) auf Quantenchemie-Hamilton-Operatoren. Vorherige Arbeiten (Ref. [1]) zeigten, dass eine Grad-sechs-SOS-Methode, die eine perturbativ gewählte Struktur kubischer Operatoren ($\tau_i$) nutzt, die Standard-Grad-vier-SOS (2RDM) sowohl in der Geschwindigkeit als als auch in der Genauigkeit für Modell-Hamilton-Operatoren mit schwachen, zufälligen Quartik-Wechselwirkungen übertreffen kann. Dieses Verfahren beruht jedoch auf der Annahme, dass Wechselwirkungsterme klein sind (Ordnung $\epsilon$). Reale Quantenchemie-Hamilton-Operatoren enthalten jedoch oft starke „Dichte-Dichte“-Wechselwirkungsterme ($U_{ij} n_i n_j$) sowie dichteabhängige Hopping- und Spin-Spin-Wechselwirkungen. Wenn diese Terme groß sind (Ordnung eins), versagt die perturbative Wahl der kubischen Operatoren, was das vorherige Verfahren für diese Anwendungen unwirksam macht. Die zentrale Herausforderung besteht darin, den selbstkonsistenten SOS-Rahmen so zu erweitern, dass er diese starken Dichte-Dichte-Wechselwirkungen handhaben kann, während die rechnerische Effizienz der Methode und ihre Fähigkeit, strikte untere Schranken für die Grundzustandsenergie zu liefern, erhalten bleiben.

Methodik
Die Autoren schlagen eine Modifikation der Definition der kubischen Operatoren ($\tau_i$) und der Struktur der SOS-Zerlegung vor, um die starken Dichte-Dichte-Wechselwirkungen zu berücksichtigen.

1. Modifizierte Hamiltonian-Struktur: Der Hamilton-Operator wird als $H = H_0 + \epsilon V$ neu definiert, wobei $H_0$ nun explizit den quadratischen Teil ($H_{quad}$) und die starken Dichte-Dichte-Wechselwirkungsterme ($\sum U_{ij} n_i n_j$) enthält. Der Term $\epsilon V$ enthält die verbleibenden, schwächeren Quartik-Wechselwirkungen.
2. Renormierte Energiedenominatoren: In dem vorherigen perturbativen Ansatz wurden die Koeffizienten der kubischen Operatoren durch Energiedenominatoren der Form $1/(e_i + e_j + e_k + e_l)$ bestimmt. Um starke Dichte-Dichte-Wechselwirkungen zu handhaben, definieren die Autoren diese Denominatoren neu, indem sie die Wechselwirkungsenergie einbeziehen. Konkret wird der Nenner zum Erwartungswert von $H_0$ in einem Zustand, in dem die Modi $i, j, k, l$ besetzt sind und alle anderen leer sind:
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{Paare}} U_{\text{paar}}$$
   Dies stellt sicher, dass der Operator $\psi_i + \tau_i$ den Grundzustand bis zu einem Fehler von $O(\epsilon^2)$ annihiliert, wodurch die „Fehler-Eigenschaft“ gewahrt bleibt, die für die Genauigkeit der Methode von $O(\epsilon^4)$ erforderlich ist.
3. Handhabung von Closure und dichteabhängigen Termen: Um die „Closure-Eigenschaft“ (die Sicherstellung, dass der effektive Hamilton-Operator innerhalb derselben Klasse von Quartik-Hamilton-Operatoren bleibt) aufrechtzuerhalten und dichteabhängige Hopping-Terme ($n_i \psi^\dagger_j \psi_k$) zu behandeln, führen die Autoren eine allgemeinere Zerlegungsstruktur ein. Sie erweitern den Trial-Hamilton-Operator um Terme der Form $\sum \epsilon n_i V_i$.
   Die neue Zerlegung nimmt die Form $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$ an, wobei:
   $$F_a = \psi_i + n_j \sum_{k \neq j} d_k \psi_k + \theta_a$$
   $$\theta_a = \tau_i + \mu + \epsilon n_j \nu$$
   Hierbei sind $\mu$ und $\nu$ zusätzliche kubische Operatoren, die gewählt werden, um die Fehler-Eigenschaft zu korrigieren und die notwendigen dichteabhängigen Terme zu erzeugen. Diese Struktur ermöglicht es, den Ziel-Hamilton-Operator $H$ über ein iteratives selbstkonsistentes Verfahren zu reproduzieren.

Wesentliche Beiträge

• Erweiterung auf starke Wechselwirkungen: Die Arbeit liefert eine Methode zur Erzeugung von SOS-Zerlegungen für Hamilton-Operatoren, die starke Dichte-Dichte-Wechselwirkungsterme enthalten, welche zuvor ein Hindernis für die Anwendung der selbstkonsistenten SOS-Methode in der Quantenchemie darstellten.
• Generalisierte Operatorendefinitionen: Die Autoren führen eine modifizierte Definition für die kubischen Operatoren $\tau_i$ unter Verwendung renormierter Energiedenominatoren ein, die die Wechselwirkungsstärken berücksichtigen, wodurch die „Fehler-Eigenschaft“ auch dann gewahrt bleibt, wenn die Wechselwirkungen nicht perturbativ klein sind.
• Erhaltung der Closure-Eigenschaft: Durch die Einführung zusätzlicher Terme unter Verwendung von Zahlungsoperatoren ($n_i$) und spezifischen kubischen Korrekturen ($\mu, \nu$) zeigen die Autoren, wie die Closure-Eigenschaft aufrechterhalten werden kann, was sicherstellt, dass der effektive Hamilton-Operator ein Quartik-Operator bleibt.
• Netzwerkfluss-Formulierung: Die Arbeit charakterisiert den Bereich der Hamilton-Operatoren, die mit dieser Methode zerlegt werden können. Sie reduziert das Problem, eine gültige Zerlegung für einen Hamilton-Operator mit negativen Wechselwirkungstermen ($U_{ij} < 0$) zu konstruieren, auf ein Netzwerkfluss-Problem. Dies liefert ein konkretes Kriterium dafür, ob ein gegebener Satz von Wechselwirkungsparametern durch die vorgeschlagene Zerlegung handhabbar ist.

Ergebnisse und Ansprüche
Die Arbeit behauptet, dass die vorgeschlagene Methode:

• Genauigkeit beibehält: Sie liefert eine untere Schranke für die Grundzustandsenergie, die im Verhältnis zur wahren Grundzustandsenergie genau $O(\epsilon^4)$ ist, selbst wenn die Dichte-Dichte-Wechselwirkungen Ordnung eins erreichen.
• Effizienz bewahrt: Der Rechenaufwand ist nicht signifikant höher als wenn alle Terme klein ($O(\epsilon)$) sind, da die Methode vermeidet, vollständige semidefiniten Programme (SDP) hohen Grades zu lösen, und stattdessen auf Matrixdiagonalisierung und iterativen Updates basiert.
• Spezifische Terme handhabt: Die Methode handhabt erfolgreich Dichte-Dichte-Wechselwirkungen und, als Nebeneffekt, dichteabhängige Hopping-Terme (der Form $n_i(\psi^\dagger_j \psi_k + h.c.)$) sowie Spin-Spin-Wechselwirkungen (die auf Summen von Dichte-Dichte-Wechselwirkungen abgebildet werden können).
• Limitierungen: Die Autoren merken an, dass die Methode dichteabhängiges Hopping derzeit nur in einer spezifischen Form handhabt (nach Teilchen-Loch-Konjugation). Terme, die Paarerzeugung beinhalten ($n_i \psi^\dagger_j \psi^\dagger_k$), sind nicht enthalten, obwohl die Autoren darauf hinweisen, dass diese typischerweise kleiner sind und als perturbativ behandelt werden können. Zudem werden Singlett-Hopping-Terme durch diese spezifische Erweiterung nicht adressiert, da sie nicht die Eigenschaft besitzen, den quadratischen Grundzustand zu annihilieren.

Bedeutung
Die Bedeutung dieser Arbeit liegt in der Überbrückung der Lücke zwischen dem theoretischen Erfolg selbstkonsistenter SOS-Methoden an zufälligen Modell-Hamilton-Operatoren und ihrer praktischen Anwendung auf reale Quantenchemie-Probleme. Durch die Überwindung der Barriere starker Dichte-Dichte-Wechselwirkungen erweitert das Paper den Anwendungsbereich einer Methode, die strikte untere Energien-Schranken bietet, ohne die hohen Kosten vollständiger SDP-Solver zu verursachen. Die Autoren positionieren diese Arbeit als notwendige Fortsetzung vorangegangener Studien, um die Behandlung physikalisch relevanter Hamilton-Operatoren zu ermöglichen, die mit diesem spezifischen nicht-perturbativen, aber perturbativ geleiteten Ansatz zuvor unerreichbar waren.