Full Uncertainty Quantification of Sign-Problem-Free Quantum Monte Carlo Methods and Nuclear Lattice Effective Field Theory Benchmarks

Résumé Technique : Quantification complète de l'incertitude des méthodes de Monte Carlo quantique sans problème de signe et des bancs d'essai de la Théorie Effective de Noyau sur Réseau

Énoncé du Problème
La théorie nucléaire ab initio moderne a réussi de manière significative à reproduire les observables pour les noyaux légers et la matière nucléaire infinie en combinant la théorie effective de champ (EFT) chirale avec des solveurs de systèmes à plusieurs corps. Cependant, l'extension de ces méthodes aux noyaux de masse moyenne et lourde (jusqu'à $^{100}$Sn) fait face à un goulot d'étranglement critique : l'absence de quantification rigoureuse des incertitudes liées aux erreurs algorithmiques systématiques. Sans un tel contrôle, les écarts entre la théorie et l'expérience sont souvent attribués aux forces nucléaires sous-jacentes, ce qui conduit au « réajustement » des constantes de couplage faible (LEC) pour compenser les déficiences du solveur. Cela masque l'origine réelle des erreurs théoriques et compromet le pouvoir prédictif des approches ab initio.

Plus précisément, des critiques récentes (Réf. [1]) ont contesté la validité des bancs d'essai de la Théorie Effective de Noyau sur Réseau (NLEFT), citant des problèmes de bornes variationnelles, d'énergies de corrélation et de mécanismes de saturation nucléaire. Ces critiques découlent souvent de comparaisons entre des théories de régularisation différentes (par exemple, le mélange des formalismes de matrice de transfert et de Hamiltonien) ou de traitements de volume fini incohérents, plutôt que de véritables échecs des méthodes à plusieurs corps.

Méthodologie
Les auteurs utilisent des méthodes de Monte Carlo quantique (QMC) sans problème de signe, spécifiquement le Monte Carlo quantique à champ auxiliaire (AFQMC), au sein du cadre de la NLEFT. L'étude se concentre sur l'action récemment développée LAT-OPT1, un Hamiltonien de réseau minimal sans problème de signe qui incorpore une interaction spin-orbite non perturbative tout en préservant la symétrie anti-unitaire requise pour des poids de Monte Carlo positifs dans les noyaux pairs-pairs.

La méthodologie implique une analyse complète du budget d'erreurs et un réexamen diagnostique des comparaisons de bancs d'essai :

1. Quantification de l'incertitude : Les auteurs quantifient systématiquement les incertitudes statistiques, de temps imaginaire, de volume fini, de discrétisation temporelle, induites par le champ auxiliaire et de signe résiduel pour des noyaux allant de $^4$He à $^{100}$Sn.
2. Traitement non perturbatif de l'interaction spin-orbite : Contrairement aux approches précédentes où les forces spin-orbite étaient traitées comme des insertions perturbatives, LAT-OPT1 intègre la densité spin-orbite directement dans l'Hamiltonien de projection dans le temps euclidien.
3. Vérifications de la cohérence du régulateur : Le document compare rigoureusement différents formalismes de calcul :
   • Matrice de transfert vs Hamiltonien : Analyse des différences entre les calculs de matrice de transfert à pas de temps sur réseau fini ($a_t$) et la limite Hamiltonienne $a_t \to 0$.
   • Conditions aux limites : Étude de l'impact des conditions aux limites périodiques (PBC) par rapport aux conditions aux limites de bord moyen (ATBC) sur les effets de couche de volume fini dans la matière nucléaire.
   • Renormalisation : Examen de la limite de continuum et de la renormalisation des interactions à courte portée.

Contributions Clés et Résultats

• Contrôle Quantitatif des Erreurs Systématiques :
  L'étude démontre que l'incertitude computationnelle globale de nombreux corps pour les énergies de l'état fondamental des noyaux doubles-magiques jusqu'à $^{100}$Sn est bien inférieure au niveau du pourcentage.

  • Temps Imaginaire et Volume Fini : Les extrapolations vers un temps de projection infini ($\tau \to \infty$) et un volume infini ($L \to \infty$) sont stables. Une taille de boîte $L=11$ (avec $a=1.32$ fm) est jugée suffisamment convergente pour les noyaux étudiés.
  • Discrétisation Temporelle : L'utilisation d'un formalisme de matrice de transfert avec un $a_t$ fini induit un Hamiltonien effectif ($H_{eff}$) qui diffère de l'Hamiltonien cible ($H$). Les auteurs montrent que bien que $H_{eff}$ dépende de $a_t$, l'erreur d'énergie suit une loi en $O(a_t^2)$, permettant un contrôle précis.
  • Termes Induits par le Champ Auxiliaire : La complétion par carré requise pour maintenir l'échantillonnage sans problème de signe pour le terme spin-orbite génère des termes induits à deux corps ($\rho_s^2$) et à trois corps. Ils sont quantifiés et s'avèrent être de petites corrections (quelques MeV), confirmant la validité de la construction.
  • Problème de Signe Résiduel : Pour les noyaux à nombre de nucléons impair ($A$) et les noyaux pairs-impairs, un léger problème de signe résiduel persiste. Cependant, la phase moyenne reste proche de l'unité (0,94–0,98), et les incertitudes statistiques restent comparables à celles des systèmes pairs-pairs.

• Rôle Non Perturbatif du Couplage Spin-Orbite :
  Une conclusion critique est que l'interaction spin-orbite doit être incluse de manière non perturbative dans la projection du temps euclidien. Traiter l'interaction uniquement comme une insertion perturbative sur une fonction d'onde à symétrie SU(4) occulte des contributions de structure de couche de plusieurs à dix MeV dans les noyaux de masse moyenne. L'intégration dans la propagation génère dynamiquement la structure de couche correcte.

• Résolution des Écarts de Banc d'Essai :
  Le document répond aux critiques récentes (Réf. [1]) en identifiant des erreurs conceptuelles et techniques dans la comparaison de différentes théories de réseau :

  • Comparaisons Inéquivalentes : Comparer les bornes de Hartree-Fock (HF) d'un Hamiltonien $a_t \to 0$ avec les résultats AFQMC d'une matrice de transfert à $a_t$ fini est invalide. La matrice de transfert à $a_t$ fini définit une théorie régulée distincte avec des interactions induites. Lorsque la même théorie régulée est utilisée de manière cohérente, le principe variationnel est respecté (les énergies HF bornent les énergies projetées).
  • Énergies de Corrélation : Les faibles énergies de corrélation dans la matière neutronique pour l'interaction « Éléments Essentiels » (EE) n'impliquent pas des fonctions d'onde triviales. Une analyse détaillée des densités irréductibles révèle une superfluidité multimodale significative (corrélations s-wave, p-wave et quartet) malgré des différences d'énergie minimes.
  • Saturation Nucléaire : L'aplatissement de l'énergie d'interaction à haute densité n'est pas un artefact de l'empilement sur le réseau. Des dérivations analytiques dans l'espace continu avec des régulateurs gaulien-invariants montrent que l'aplatissement de l'énergie d'interaction est une conséquence naturelle de la régulation dans l'espace des moments. Les propriétés de saturation des interactions NLEFT suivent l'échelle du régulateur induit par l'étalement, et non l'échelle d'occupation des sites du réseau.
  • Limite de Continuum : L'affirmation selon laquelle les interactions sur réseau mènent à un effondrement du système dans la limite de continuum est réfutée. L'argument d'échelle utilisé dans la critique repose sur une renormalisation incorrecte au niveau de Born ; une renormalisation appropriée par rapport aux observables de diffusion physique produit un comportement de mise à l'échelle différent ($C_c \propto a$) qui empêche l'effondrement.

Signification et Revendications
Le document affirme établir une base solide pour une analyse complète des incertitudes dans les simulations NLEFT avec des forces chirales complètes. En quantifiant rigoureusement les incertitudes et en démontrant que les méthodes sans problème de signe peuvent atteindre une précision sub-percent pour les noyaux lourds, ce travail valide la NLEFT comme un banc d'essai fiable pour les méthodes à plusieurs corps.

Les auteurs soutiennent que les divergences apparentes dans la littérature récente sont résolues lorsque les comparaisons sont faites entre des théories régulées équivalentes avec des paramètres, des volumes et des conditions aux limites cohérents. L'action LAT-OPT1 est présentée non pas simplement comme un outil phénoménologique, mais comme un modèle de référence non perturbatif et exactement soluble qui capture la dynamique essentielle de la structure de couche. Ce cadre permet de séparer les déficiences réelles des forces nucléaires des erreurs de solveurs à plusieurs corps, offrant une voie pour des études prédictives de l'évolution des couches, de la spectroscopie et de la matière dense dans les noyaux de masse moyenne et lourde.