Efficient classical simulation of Fermi-Hubbard models with Kramers symmetry

Énoncé du problème

Les méthodes de Monte Carlo quantique (QMC) sont l'approche classique standard pour simuler les systèmes fermioniques en interaction, pourtant les garanties théoriques rigoureuses concernant leur efficacité sont rares. Bien que l'absence de « problème de signe » soit une condition nécessaire pour une simulation QMC efficace, elle n'est pas suffisante pour garantir la convergence en temps polynomial des chaînes de Markov sous-jacentes. De plus, les résultats rigoureux existants pour la simulation classique efficace des systèmes fermioniques sont largement confinés aux régimes de haute température, aux limites perturbatives ou à des correspondances spécifiques avec des systèmes de spins. Il existe un fossé important dans la compréhension de la complexité computationnelle des systèmes fermioniques fortement corrélés, même lorsqu'ils sont connus pour être exempts de problème de signe.

Ce travail aborde l'intersection spécifique de deux propriétés structurelles : la symétrie de renversement du temps de Kramers (qui assure une expansion exempte de problème de signe) et le théorème de l'absence de zéros de Lee-Yang (qui concerne l'emplacement des zéros de la fonction de partition). Les auteurs étudient si cette intersection permet des algorithmes classiques prouvablement efficaces pour calculer les fonctions de partition et les énergies du fondamental pour de larges familles de modèles de Hubbard de Fermi.

Méthodologie

Les auteurs développent un nouvel algorithme classique randomisé basé sur une chaîne de Markov de type « worm » de déterminant (determinantal worm Markov chain). La méthodologie se déroule en trois étapes :

1. Trotterisation et expansion d'interaction :
   La fonction de partition $Z = \text{Tr}(e^{-\beta H})$ est approximée par une décomposition de Trotter en temps discret. Le hamiltonien $H$ est divisé en une partie libre $H_0$ et des termes d'interaction $V_i$. L'expansion de $\text{Tr}(L^M)$ (où $L$ est l'étape de Trotter) est réécrite comme une somme sur des sous-ensembles d'« insertions » (termes d'interaction). En utilisant le théorème de Wick, chaque terme de cette somme est exprimé comme le déterminant d'une sous-matrice principale d'une matrice de fonction de Green globale $G$.

   • Positivité du signe : Les auteurs démontrent que pour les hamiltoniens appartenant à la « classe de Kramers attractive », la matrice $G$ est quaternioniènne. Une propriété clé des matrices quaternioniènnes est que leurs valeurs propres vont par paires complexes conjuguées, garantissant que le déterminant de toute sous-matrice carrée formée par des blocs de lignes et de colonnes entiers est non négatif. Cela établit que les poids dans l'expansion sont non négatifs, évitant ainsi le problème de signe.

2. Espace de configuration élargi (Le Worm) :
   L'échantillonnage direct des mineurs principaux (où les indices de lignes et de colonnes sont identiques) est difficile. Pour faciliter un échantillonnage efficace, les auteurs introduisent un espace de configuration élargi. Ils ajoutent un petit « champ de couplage » au hamiltonien, permettant des configurations où l'ensemble des blocs de lignes sélectionnés ($A$) et des blocs de colonnes ($B$) diffèrent, à condition que $|A| = |B|$.

   • Le poids d'une configuration $(A, B)$ est proportionnel à $\lambda^{|A \setminus B|} \det G[A, B]$, où $\lambda$ est un petit paramètre.
   • Les configurations valides (où $A=B$) correspondent aux termes originaux de la fonction de partition. La masse de probabilité de ces configurations valides reste constante (bornée inférieurement par une constante) même lorsque $\lambda \to 0$, si $\lambda$ est choisi de manière appropriée.

3. Mixage rapide via la log-concavité :
   La contribution technique centrale est de prouver que la chaîne de Markov définie sur cet espace élargi se mélange en temps polynomial.

   • Les auteurs construisent un polynôme générateur pour les déterminants $\det G[A, B]$. Ils prouvent que pour les matrices quaternioniènnes, ce polynôme est fortement log-concave.
   • Ils exploitent les résultats récents d'Anari, Liu, Oveis Gharan et Vinzant [ALOGV24] ainsi que de Chen et Liu [CL26], qui montrent que les chaînes d'échange de bases (et leurs variantes « inclinées » ou tilted) se mélangent rapidement pour les mesures possédant des polynômes générateurs fortement log-concaves.
   • Le champ de couplage introduit un facteur d'inclinaison $\lambda^{|A \setminus B|}$. Les auteurs vérifient que cette inclinaison satisfait les conditions de courbure spectrale locale requises par Chen et Liu [CL26] pour préserver un mixage rapide.

Contributions clés

• Algorithme en temps polynomial : L'article fournit la première preuve rigoureuse d'un algorithme classique en temps polynomial pour estimer la fonction de partition de larges familles de systèmes fermioniques en interaction (spécifiquement, les hamiltoniens invariants de Kramers attractifs).
• Chaîne de Worm de déterminant : Les auteurs introduisent une chaîne de Markov de type « worm de déterminant » qui étend le concept de l'algorithme worm (précédemment utilisé pour les spins et les bosons) aux systèmes fermioniques possédant la symétrie de Kramers. Contrairement aux tentatives précédentes de worm fermionique, cette technique est systématiquement applicable à l'ensemble de la classe de Kramers avec des garanties de convergence prouvables.
• Connexion des propriétés structurelles : Le travail relie rigoureusement l'absence de problème de signe (via la symétrie de Kramers) avec la propriété de non-existence de zéros de Lee-Yang (via la log-concavité du polynôme générateur) pour établir une simulabilité classique efficace.
• Étendue de l'applicabilité : L'algorithme s'applique aux :
  • Modèles de Hubbard de Fermi attractifs avec des amplitudes de saut réelles sur n'importe quel graphe.
  • Modèles de Hubbard de Fermi répulsifs à demi-remplissage sur des graphes bipartites (via une transformation particule-trou).
  • Crucialement, il n'y a aucune restriction sur la force ou l'extensivité des interactions.

Résultats

Le théorème principal stipule que pour tout hamiltonien $H$ dans la classe de Kramers attractive, il existe un algorithme classique randomisé qui estime la fonction de partition $Z$ avec une erreur relative $\epsilon$ en un temps $\text{poly}(N, r, \beta, \Lambda, 1/\epsilon)$, où $N$ est le nombre de modes, $r$ est le nombre de termes d'interaction, $\beta$ est l'inverse de la température et $\Lambda$ borne la force de l'interaction.

Par conséquent, l'algorithme peut estimer efficacement :

• L'énergie libre $F_\beta$ avec une erreur additive $\epsilon$.
• L'énergie du fondamental $E_0$ avec une erreur additive $\epsilon$ (en prenant la limite $\beta \to \infty$).

Le gap spectral de la chaîne de Markov est prouvé être au moins égal à $\lambda / (2|K|)$, assurant un temps de mixage polynomial.

Signification et affirmations

Les auteurs affirment que ce travail constitue les premières garanties de temps polynomial pour le Monte Carlo quantique fermionique. Bien que l'algorithme ne soit pas une variante du Monte Carlo quantique de déterminant (DQMC) lui-même, il valide rigoureusement des décennies de succès empiriques du DQMC et de ses variantes sur des systèmes possédant la symétrie de Kramers.

L'article postule que l'intersection de la propriété de Lee-Yang et de la stoquasticité (ou l'absence de problème de signe) peut constituer un principe directeur plus large pour identifier les systèmes quantiques qui admettent une simulation classique efficace. Les auteurs notent que bien qu'ils aient résolu le problème pour la classe de Kramers attractive, d'autres systèmes fermioniques exempts de problème de signe (par exemple, ceux possédant une positivité de réflexion de Majorana) qui possèdent également la propriété de Lee-Yang restent ouverts à l'exploration future. Ce travail suggère que les propriétés structurelles permettant une simulation efficace ne sont pas limitées aux systèmes de spins mais s'étendent à la matière fermionique fortement corrélée.