Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

Résumé Technique : Interaction de Configuration Complète par Méthode Monte Carlo pour des Calculs de Structure Nucléaire Ab Initio Précis

Énoncé du Problème
Le développement de méthodes de corps multiples ab initio pour les noyaux atomiques fortement interagissants est un défi central de la théorie nucléaire moderne. Bien que la théorie du champ effectif chirale (EFT) fournisse un cadre systématique pour dériver les interactions nucléaires à partir de la chromodynamique quantique, la résolution de l'équation de Schrödinger à corps multiples pour les noyaux de masse moyenne et lourde reste informatiquement prohibitive.

• Modèle de Coquille Sans Noyau (NCSM) : Bien que fructueux pour les noyaux légers, le NCSM est limité par la croissance exponentielle de la dimension de l'espace de modèle, qui dépasse rapidement la capacité de calcul des supercalculateurs actuels (atteignant $\sim 10^{12}$ dimensions). Par conséquent, le NCSM repose lourdement sur des extrapolations pour les systèmes plus importants.
• Expansions de Corps Multiples Tronquées : Des méthodes telles que la théorie de la perturbation de corps multiples (MBPT), le groupe de renormalisation de similitude en milieu intérieur (IMSRG) et la théorie de Coupled-Cluster (CC) offrent une mise à l'échelle polynomiale, permettant des calculs pour des noyaux plus lourds. Cependant, elles souffrent d'erreurs de troncature. L'obtention d'une grande précision nécessite de pousser ces expansions à des ordres plus élevés, ce qui augmente la complexité mathématique et le coût computationnel. De plus, fournir des estimations fiables des incertitudes de corps multiples dans ces schémas tronqués est difficile, particulièrement pour les systèmes présentant de fortes corrélations de regroupement (clustering) ou collectives, où la convergence peut être lente.

Il existe un besoin critique pour une méthode capable de conserver la précision d'une Interaction de Configuration Complète (FCI) sans le stockage explicite du Hamiltonien, servant de référence pour valider et guider le développement des méthodes d'expansion tronquées.

Méthodologie
Les auteurs introduisent la méthode de l'Interaction de Configuration Complète par Méthode Monte Carlo (FCIQMC) dans le domaine de la structure nucléaire ab initio. La FCIQMC est une méthode de projection stochastique qui résout l'équation de Schrödinger en temps imaginaire :
$$-\frac{d}{d\tau} |\Psi(\tau)\rangle = (\hat{H} - E_0) |\Psi(\tau)\rangle$$
Les caractéristiques clés de l'implémentation incluent :

• Échantillonnage Stochastique : La fonction d'onde est représentée par un ensemble de « marcheurs » (particules signées) distribués à travers des déterminants de Slater dans l'espace de configuration complet. Le coefficient d'un déterminant est déterminé par la population nette signée de marcheurs sur ce déterminat.
• Étapes de Propagation : L'ensemble des marcheurs évolue à travers trois étapes stochastiques par pas de temps imaginaire ($\Delta\tau$) :
  1. Spawn (Éclosion) : Échantillonne les couplages de Hamiltonien hors-diagonaux ($H_{ij}$) pour créer de nouveaux marcheurs sur les déterminants connectés.
  2. Mort/Clonage : Échantillonne le terme diagonal ($H_{ii}$) décalé par un paramètre $S$, ajustant la population de marcheurs sur le déterminant actuel.
  3. Annihilation : Annule les marcheurs de signes opposés résidant sur le même déterminant, ce qui est crucial pour surmonter le problème de signe de Fermi sans nécessiter d'approximation de nœud fixe ou de fonction d'onde de guidage.
• Approximation de l'Initiateur : Pour contrôler le problème de signe dans les grands espaces de modèles, les auteurs emploient l'approximation de l'initiateur combinée à une correction de décalage adaptatif. Cela permet d'obtenir des estimations contrôlables et améliorables de la limite FCI à mesure que la population de marcheurs ($N_w$) augmente.
• Configuration du Système : Les calculs ont été réalisés pour $^4$He, $^8$Be, $^{12}$C et $^{16}$O en utilisant une interaction deux-nucléons chirale optimisée (N2LOopt) dans une base d'oscillateur harmonique sphérique ($\hbar\omega = 20$ MeV). Aucune troncature de configuration n'a été appliquée dans les calculs FCIQMC.

Résultats Clés
L'étude présente les énergies de l'état fondamental ($E$) et les rayons de charge ($R_{ch}$) pour les noyaux sélectionnés avec des incertitudes de corps multiples de l'ordre du sous-percentile.

• Validation dans de Petits Espaces : Dans un petit espace de modèle ($e_{max}=2$), les résultats FCIQMC concordent avec les résultats FCI exacts obtenus par diagonalisation directe, validant la précision de l'échantillonnage stochastique.
• Performance dans de Grands Espaces de Modèles : Dans un grand espace de modèle ($e_{max}=10$), où la dimension de l'espace de Hilbert pour $^{16}$O atteint $\sim 10^{33}$ (rendant la diagonalisation directe impossible), la FCIQMC a fourni des résultats précis.
  • Pour $^4$He, les résultats FCIQMC concordent parfaitement avec les extrapolations quasi-exactes du NCSM.
  • Pour les noyaux plus lourds ($^8$Be, $^{12}$C, $^{16}$O), où le NCSM ne parvient pas à converger, la FCIQMC a délivré des résultats précis.
• Évaluation des Méthodes Tronquées :
  • MBPT : A montré les plus grands écarts par rapport aux résultats FCIQMC.
  • IMSRG : Généralement raisonnable mais a présenté des écarts significatifs pour $^8$Be, où les corrections d'ordre supérieur (IMSRG(3f2)+T) s'écartaient encore davantage des valeurs exactes, indiquant l'absence de corrélations collectives.
  • Coupled-Cluster (CC) : Montre une amélioration systématique de CCSD vers CCSDT-3. Cependant, pour le rayon de charge de $^{16}$O, l'inclusion des triples itératifs dominants a conduit à un écart plus important par rapport à la FCIQMC, soulignant la difficulté de convergence pour certains observables.
• Effets de Regroupement (Clustering) : Les résultats pour $^8$Be, qui possède une structure moléculaire/de regroupement prononcée, ont démontré que les expansions basées sur une référence (comme l'IMSRG à référence unique) peinent à converger en raison de l'omission des corrélations collectives d'ordre élevé. La FCIQMC, en incluant les corrélations à tous les rangs d'excitation au sein de la base, maintient sa précision pour ces états regroupés.

Signification et Revendications
L'article affirme que la FCIQMC établit un nouveau standard pour les calculs de structure nucléaire ab initio :

1. Solveur de Plein Espace Stochastique : Elle est capable de traiter des systèmes dépassant la portée du NCSM conventionnel en échantillonnant directement la fonction d'onde exacte dans l'espace de configuration complet.
2. Incertitudes Quantifiées : Contrairement aux méthodes d'expansion tronquées où l'estimation de l'incertitude est difficile, la FCIQMC fournit des résultats avec des incertitudes de corps multiples quantifiées et inférieures à un pour cent.
3. Référence pour la Théorie : Les résultats servent de références exactes pour évaluer la précision des interactions chirales et la convergence des méthodes de corps multiples à mise à l'échelle polynomiale (MBPT, IMSRG, CC). L'étude démontre explicitement le rôle crucial des corrélations d'ordre supérieur qui font souvent défaut dans les schémas tronqués.
4. Impact Futur : Les auteurs postulent que cette méthode facilitera le développement de méthodes de corps multiples de prochaine génération avec une plus grande précision et aidera à résoudre des énigmes de longue date en théorie nucléaire, telles que l'évolution des rayons de charge dans la chaîne du calcium et les corrections de rupture d'isospin dans les désintégrations Fermi superpermises, où les effets de corps multiples d'ordre élevé sont centraux.

Ce travail représente la première application de la FCIQMC aux calculs ab initio précis de noyaux finis en interactions fortes, marquant une étape importante vers la compréhension des forces nucléaires et des symétries fondamentales avec une précision sans précédent.