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

Technische Samenvatting: Verbetering van Perturbatietheorie met de Som-van-kwadraten II: Grote Dichtheids-dichtheids-termen

Probleemstelling
Het artikel behandelt een significante beperking bij het toepassen van zelfconsistente som-van-kwadraten (SOS) decomposities op kwantumchemische Hamiltonians. Voorgaand werk (Ref. [1]) toonde aan dat een graad-zes SOS-methode, die gebruikmaakt van een perturbatieve keuze van cubische operatoren ($\tau_i$), een standaard graad-vier SOS (2RDM) kon overtreffen in zowel snelheid als nauwkeurigheid voor model-Hamiltonians met zwakke, willekeurige quartische interacties. Deze methode rust echter op de aanname dat interactietermen klein zijn (orde $\epsilon$). Realistische kwantumchemische Hamiltonians bevatten vaak sterke "dichtheids-dichtheids" interactietermen ($U_{ij} n_i n_j$), evenals dichtheidsafhankelijke hopping en spin-spin interacties. Wanneer deze termen groot zijn (orde eenheid), faalt de perturbatieve keuze van cubische operatoren om deze accuraat te behandelen, wat de vorige methode ineffectief maakt voor dergelijke toepassingen. De kernuitdaging is het uitbreiden van het zelfconsistente SOS-raamwerk om deze sterke dichtheids-dichtheids termen te kunnen hanteren, terwijl de computationele efficiëntie van de methode en het vermogen om strikte ondergrenzen op de grondstoestandsenergie te bieden, behouden blijven.

Methodologie
De auteurs stellen een modificatie voor van de definitie van de cubische operatoren ($\tau_i$) en de structuur van de SOS-decompositie om sterke dichtheids-dichtheids termen te accommoderen.

1. Gemodificeerde Hamiltonian-structuur: De Hamiltonian wordt hergedefinieerd als $H = H_0 + \epsilon V$, waarbij $H_0$ nu expliciet de kwadratische component ($H_{quad}$) en de sterke dichtheids-dichtheids interactietermen ($\sum U_{ij} n_i n_j$) bevat. De term $\epsilon V$ bevat de resterende zwakkere quartische interacties.
2. Gerenormaliseerde Energiedenominatoren: In de vorige perturbatieve benadering werden de coëfficiënten van de cubische operatoren bepaald door energiedenominatoren van de vorm $1/(e_i + e_j + e_k + e_l)$. Om sterke dichtheids-dichtheids termen te hanteren, herdefiniëren de auteurs deze denominatoren om de interactie-energie te bevatten. Specifiek wordt de denominator de verwachtingswaarde van $H_0$ in een toestand waarin modi $i, j, k, l$ bezet zijn en anderen leeg:
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{paren}} U_{pair}$$
   Dit zorgt ervoor dat de operator $\psi_i + \tau_i$ de grondtoestand nog steeds annihileert tot een fout van $O(\epsilon^2)$, waardoor de "fout-eigenschap" behouden blijft die vereist is voor de methode om een $O(\epsilon^4)$ nauwkeurige energie-ondergrens te leveren.
3. Afhandeling van Closure en Dichtheidsafhankelijke Termen: Om de "closure-eigenschap" te handhaven (het waarborgen dat de effectieve Hamiltonian binnen dezelfde klasse van quartische Hamiltonians blijft) en om dichtheidsafhankelijke hopping-termen ($n_i \psi^\dagger_j \psi_k$) te behandelen, introduceren de auteurs een meer algemene decompositiestructuur. Ze breiden de trial-Hamiltonian uit met termen van de vorm $\sum \epsilon n_i V_i$.
   De nieuwe decompositie neemt de vorm aan $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$, waarbij:
   $$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$$
   Hier zijn $\mu$ en $\nu$ aanvullende cubische operatoren die gekozen zijn om de fout-eigenschap te corrigeren en de noodzakelijke dichtheidsafhankelijke termen te genereren, respectievelijk. Deze structuur maakt het mogelijk om de doel-Hamiltonian $H$ te reproduceren via een iteratief zelfconsistent proces.

Belangrijkste Bijdragen

• Extensie naar Sterke Interacties: Het artikel biedt een methode om SOS-decomposities te genereren voor Hamiltonians die sterke dichtheids-dichtheids interactietermen bevatten, wat voorheen een barrière vormde voor de toepassing van de zelfconsistente SOS-methode op kwantumchemie.
• Gegeneraliseerde Operatordefinities: De auteurs introduceren een gemodificeerde definitie voor de cubische operatoren $\tau_i$ met behulp van gerenormaliseerde energiedenominatoren die rekening houden met interactiekrachten, waardoor de "fout-eigenschap" behouden blijft, zelfs wanneer interacties niet perturbatief klein zijn.
• Behoud van de Closure-eigenschap: Door aanvullende termen te introduceren die betrokken zijn bij getaloperatoren ($n_i$) en specifieke cubische correcties ($\mu, \nu$), demonstreren de auteurs hoe de closure-eigenschap behouden kan blijven, wat garandeert dat de effectieve Hamiltonian een quartische operator blijft die geschikt is voor de iteratieve oplossing.
• Netwerkstroomformulering: Het artikel karakteriseert het bereik van Hamiltonians die met deze methode gedecomposeerd kunnen worden. Het reduceert het probleem van het construeren van een geldige decompositie voor een Hamiltonian met negatieve interactietermen ($U_{ij} < 0$) tot een netwerkstroomprobleem. Dit biedt een concreet criterium om te bepalen of een gegeven set interactieparameters door de voorgestelde decompositie afgehandeld kan worden.

Resultaten en Claims
Het artikel claimt dat de voorgestelde methode:

• Nauwkeurigheid Handhaaft: Het levert een ondergrens op de grondstoestandsenergie die nauwkeurig is tot $O(\epsilon^4)$ relatief aan de ware grondstoestandsenergie, zelfs wanneer dichtheids-dichtheids termen van de orde eenheid zijn.
• Efficiëntie Behoudt: De computationele inspanning die vereist is, is niet significant hoger dan wanneer alle termen klein zijn ($O(\epsilon)$), aangezien de methode het oplossen van volledige semi-definite programma's (SDP's) van hoge graad vermijdt en vertrouwt op matrixdiagonalisatie en iteratieve updates.
• Specifieke Termen Afhandelt: De methode handelt succesvol dichtheids-dichtheids interacties af en, als bijproduct, ook dichtheidsafhankelijke hopping-termen (van de vorm $n_i(\psi^\dagger_j \psi_k + h.c.)$) en spin-spin interacties (die gemapt kunnen worden naar sommen van dichtheids-dichtheids termen).
• Beperkingen: De auteurs merken op dat de methode momenteel dichtheidsafhankelijke hopping alleen in een specifieke vorm afhandelt (na deeltjes-gat conjugatie). Termen die betrokken zijn bij paarcreatie ($n_i \psi^\dagger_j \psi^\dagger_k$) zijn niet inbegrepen, hoewel de auteurs suggereren dat deze doorgaans kleiner zijn en als perturbatief behandeld kunnen worden. Daarnaast worden singlet hopping-termen niet geadresseerd door deze specifieke uitbreiding, omdat zij niet de eigenschap delen van het annihileren van de kwadratische grondtoestand.

Significantie
De significantie van dit werk ligt in het overbruggen van de kloof tussen het theoretische succes van zelfconsistente SOS-methoden op willekeurige model-Hamiltonians en hun praktische toepassing op echte kwantumchemische problemen. Door de barrière van sterke dichtheids-dichtheids interacties te overwinnen, breidt het artikel het toepasbare bereik van een methode uit die strikte ondergrenzen op energie biedt zonder de computationele kosten van volledige SDP-solvers. De auteurs positioneren dit werk als een noodzakelijk vervolg op eerdere studies, wat de behandeling mogelijk maakt van fysiek relevante Hamiltonians die voorheen buiten bereik lagen voor deze specifieke niet-perturbatieve, perturbatief-gestuurde benadering.