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

Résumé technique : Amélioration de la théorie de la perturbation avec la somme de carrés II : Termes de densité-densité importants

Énoncé du problème
L'article traite d'une limitation significative de l'application des décompositions auto-cohérentes de sommes de carrés (SOS) aux hamiltoniens de la chimie quantique. Des travaux antérieurs (Réf. [1]) ont démontré qu'une méthode SOS de degré six, utilisant un choix perturbatif d'opérateurs cubiques ($\tau_i$), pouvait surpasser la méthode SOS standard de degré quatre (2RDM) en termes de vitesse et de précision pour des hamiltoniens de modèle présentant de faibles interactions quartiques aléatoires. Cependant, cette méthode repose sur l'hypothèse que les termes d'interaction sont faibles (d'ordre $\epsilon$). Les hamiltoniens de la chimie quantique réelle contiennent souvent des termes d'interaction « densité-densité » forts ($U_{ij} n_i n_j$), ainsi que des sauts dépendants de la densité et des interactions spin-spin. Lorsque ces termes sont importants (d'ordre unité), le choix perturbatif des opérateurs cubiques échoue à traiter ces termes de manière précise, rendant la méthode précédente inefficace pour ces applications. Le défi central consiste à étendre le cadre SOS auto-cohérent pour prendre en compte ces termes d'interaction forts tout en maintenant l'efficacité computationnelle de la méthode et sa capacité à fournir des bornes inférieures strictes sur l'énergie du fondamental.

Méthodologie
Les auteurs proposent une modification de la définition des opérateurs cubiques ($\tau_i$) et de la structure de la décomposition SOS pour accommoder les termes de densité-densité forts.

1. Structure de l'Hamiltonien Modifiée : L'hamiltonien est redéfini comme $H = H_0 + \epsilon V$, où $H_0$ inclut désormais explicitement la partie quadratique ($H_{quad}$) et les termes d'interaction densité-densité forts ($\sum U_{ij} n_i n_j$). Le terme $\epsilon V$ contient les autres interactions quartiques plus faibles.
2. Renormalisation des dénominateurs d'énergie : Dans l'approche perturbative précédente, les coefficients des opérateurs cubiques étaient déterminés par des dénominateurs d'énergie de la forme $1/(e_i + e_j + e_k + e_l)$. Pour gérer les termes de densité-densité forts, les auteurs redéfinissent ces dénominateurs pour inclure l'énergie d'interaction. Spécifiquement, le dénominateur devient l'espérance mathématique de $H_0$ dans un état où les modes $i, j, k, l$ sont occupés et les autres vides :
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{paires}} U_{pair}$$
   Cela garantit que l'opérateur $\psi_i + \tau_i$ annule toujours l'état fondamental jusqu'à une erreur de $O(\epsilon^2)$, préservant la « propriété d'erreur » requise pour que la méthode produise une limite d'énergie précise à $O(\epsilon^4)$.
3. Gestion de la clôture et des termes dépendants de la densité : Pour maintenir la « propriété de clôture » (garantissant que l'hamiltonien effectif reste dans la même classe d'hamiloniens quartiques) et pour gérer les termes de saut dépendants de la densité ($n_i \psi^\dagger_j \psi_k$), les auteurs introduisent une structure de décomposition plus générale. Ils étendent l'hamiltonien d'essai pour inclure des termes de la forme $\sum \epsilon n_i V_i$.
   La nouvelle décomposition prend la forme $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$, où :
   $$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$$
   Ici, $\mu$ et $\nu$ sont des opérateurs cubiques additionnels choisis pour corriger la propriété d'erreur et générer les termes dépendants de la densité nécessaires, respectivement. Cette structure permet de reproduire l'hamiltonien cible $H$ via une procédure itérative auto-cohérente.

Contributions clés

• Extension aux interactions fortes : L'article fournit une méthode pour générer des décompositions SOS pour des hamiltoniens contenant des termes d'interaction densité-densité forts, ce qui constituait auparavant une barrière à l'application de la méthode SOS auto-cohérente à la chimie quantique.
• Définitions d'opérateurs généralisées : Les auteurs introduisent une définition modifiée des opérateurs cubiques $\tau_i$ utilisant des dénominateurs d'énergie renormalisés qui tiennent compte des forces d'interaction, garantissant que la « propriété d'erreur » est respectée même lorsque les interactions ne sont pas perturbativement petites.
• Préservation de la propriété de clôture : En introduisant des termes auxiliaires impliquant des opérateurs de nombre ($n_i$) et des corrections cubiques spécifiques ($\mu, \nu$), les auteurs démontrent comment maintenir la propriété de clôture, assurant que l'hamiltonien effectif reste un opérateur quartique adapté à la solution itérative.
• Formulation de flux réseau : L'article caractérise la gamme d'hamiloniens pouvant être décomposés à l'aide de cette méthode. Il réduit le problème de construction d'une décomposition valide pour un hamiltonien possédant des termes d'interaction négatifs ($U_{ij} < 0$) à un problème de flux réseau. Cela fournit un critère concret pour déterminer si un ensemble de paramètres d'interaction donnés peut être traité par la décomposition proposée.

Résultats et affirmations
L'article affirme que la méthode proposée :

• Maintient la précision : Elle produit une borne inférieure sur l'énergie du fondamental qui est précise à $O(\epsilon^4)$ par rapport à l'énergie du fondamental réel, même lorsque les termes de densité-densité sont d'ordre unité.
• Préserve l'efficacité : L'effort de calcul requis n'est pas significativement plus élevé que lorsque tous les termes sont petits ($O(\epsilon)$), car la méthode évite de résoudre des programmes semi-définis (SDP) complets de haut degré, s'appuyant plutôt sur la diagonalisation de matrices et des mises à jour itératives.
• Gère des termes spécifiques : La méthode traite avec succès les interactions densité-densité et, par effet secondaire, les termes de saut dépendants de la densité (de la forme $n_i(\psi^\dagger_j \psi_k + h.c.)$) ainsi que les interactions spin-spin (qui peuvent être mappées sur des sommes de termes densité-densité).
• Limitations : Les auteurs notent que la méthode ne traite actuellement les sauts dépendants de la densité que sous une forme spécifique (après conjugaison particule-trou). Les termes impliquant la création de paires ($n_i \psi^\dagger_j \psi^\dagger_k$) ne sont pas inclus, bien que les auteurs suggèrent que ceux-ci sont typiquement plus petits et peuvent être traités de manière perturbative. De plus, les termes de saut de singulet ne sont pas abordés par cette extension spécifique, car ils ne partagent pas la propriété d'annihiler l'état fondamental quadratique.

Signification
La signification de ce travail réside dans le fait de combler l'écart entre le succès théorique des méthodes SOS auto-cohérentes sur des hamiltoniens de modèles aléatoires et leur application pratique à de réels problèmes de chimie quantique. En surmontant la barrière des interactions densité-densité fortes, l'article étend la portée d'application d'une méthode qui offre des bornes inférieures strictes sur l'énergie sans le coût computationnel des solveurs SDP complets. Les auteurs positionnent ce travail comme une suite nécessaire d'études antérieures, permettant le traitement d'hamiloniens physiquement pertinents qui étaient auparavant hors de portée de cette approche spécifique guidée par la théorie des perturbations non-perturbative.