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

Technische Samenvatting: Volledige Onzekerheidskwantificering van Sign-Problem-Free Quantum Monte Carlo-methoden en Nucleaire Lattice Effective Field Theory Benchmarks

Probleemstelling
Moderne ab initio kernfysica heeft aanzienlijk succes geboekt bij het reproduceren van grootheden voor lichte kernen en oneindige nucleaire materie door chirale effectieve veldentheorie (EFT) te combineren met many-body solvers. Het uitbreiden van deze methoden naar middelzware en zware kernen (tot $^{100}$Sn) staat echter voor een kritieke flessenhals: het gebrek aan rigoureuze onzekerheidskwantificering voor systematische algoritmische fouten. Zonder dergelijke controle worden discrepanties tussen theorie en experiment vaak toegeschreven aan de onderliggende nucleaire krachten, wat leidt tot de "heraanpassing" van lage-energieconstanten (LEC's) om tekortkomingen in de solver te compenseren. Dit maskeert de werkelijke oorsprong van theoretische fouten en ondermijnt de voorspellende kracht van ab initio benaderingen.

Specifiek hebben recente kritieken (Ref. [1]) de geldigheid van Nuclear Lattice Effective Field Theory (NLEFT) benchmarks uitgedaagd, waarbij zij wijzen op problemen met variationele grenzen, correlatie-energieën en nucleaire verzadigingsmechanismen. Deze kritieken komen vaak voort uit het vergelijken van resultaten die afkomstig zijn van verschillende gereguleerde lattice-theorieën (bijv. het mengen van transfer-matrix en Hamiltonian-formalismen) of inconsistente eindige-volume behandelingen, in plaats van werkelijke mislukkingen van de many-body methoden.

Methodologie
De auteurs maken gebruik van sign-problem-free Quantum Monte Carlo (QMC) methoden, specifiek Auxiliary-Field Quantum Monte Carlo (AFQMC), binnen het NLEFT-kader. De studie richt zich op de recent ontwikkelde LAT-OPT1 actie, een minimale sign-problem-free lattice Hamiltonian die een niet-perturbatieve spin-orbitaal interactie incorporeert terwijl de anti-unitaire symmetrie behouden blijft die vereist is voor positieve Monte Carlo-gewichten in even-oneven kernen.

De methodologie omvat een uitgebreide foutenbudgetanalyse en een diagnostische herbeoordeling van benchmark-vergelijkingen:

1. Onzekerheidskwantificering: De auteurs kwantificeren systematisch statistische, imaginaire-tijd, eindige-volume, temporele-discretisatie, auxiliary-field-geïnduceerde en residuele sign-onzekerheden voor kernen variërend van $^4$He tot $^{100}$Sn.
2. Niet-perturbatieve Spin-Orbitaal Behandeling: In tegen tegenstelling tot eerdere benaderingen waar spin-orbitaal krachten werden behandeld als perturbatieve inserties, implementeert LAT-OPT1 de spin-orbitaal dichtheid direct in de Euclidean-tijd projectie Hamiltonian.
3. Regulator Consistentiecontroles: Het artikel vergelijkt rigoureus verschillende computationele formalismen:
   • Transfer-Matrix versus Hamiltonian: Het analyseren van de verschillen tussen eindige temporele lattice-afstand ($a_t$) transfer-matrix berekeningen en de $a_t \to 0$ Hamiltonian limiet.
   • Randvoorwaarden: Het onderzoeken van de impact van Periodieke Randvoorwaarden (PBC) versus Average Twisted Boundary Conditions (ATBC) op eindige-volume shell-effecten in nucleaire materie.
   • Renormalisatie: Het onderzoeken van de continuüm-limiet en de renormalisatie van kort-bereik interacties.

Belangrijkste Bijdragen en Resultaten

• Kwantitatieve Controle van Systematische Fouten:
  De studie toont aan dat de cumulatieve many-body computationele onzekerheid in grondtoestandsenergieën voor dubbel magische kernen tot $^{100}$Sn ruim onder het procentuele niveau ligt.

  • Imaginaire Tijd & Eindige Volume: Extrapolaties naar oneindige projectietijd ($\tau \to \infty$) en oneindig volume ($L \to \infty$) zijn stabiel. Een box-grootte van $L=11$ (met $a=1.32$ fm) blijkt voldoende geconvergeerd te zijn voor de bestudeerde kernen.
  • Temporele Discretisatie: Het gebruik van een transfer-matrix formalisme met een eindige $a_t$ induceert een effectieve Hamiltonian ($H_{eff}$) die verschilt van de doel-Hamiltonian ($H$). De auteurs tonen aan dat hoewel $H_{eff}$ afhankelijk is van $a_t$, de energie-fout schaalt als $O(a_t^2)$, wat precieze controle mogelijk maakt.
  • Auxiliary-Field Geïnduceerde Termen: De square-completion die nodig is om sign-problem-free sampling voor de spin-orbitaal term te handhaven, genereert geïnduceerde twee-body ($\rho_s^2$) en drie-body termen. Deze zijn gekwantificeerd en bevonden als kleine correcties (enkele MeV), wat de validiteit van de constructie bevestigt.
  • Residueel Sign-Probleem: Voor odd-$A$ en odd-odd kernen blijft een mild residueel sign-probleem bestaan. Echter, de gemiddelde fase blijft dicht bij één (0.94–0.98), en de statistische onzekerheden blijven vergelijkbaar met even-even systemen.

• Niet-perturbatieve Rol van Spin-Orbitaal Koppeling:
  Een cruciale bevinding is dat de spin-orbitaal interactie niet-perturbatief moet worden opgenomen in de Euclidean-tijd projectie. Het behandelen ervan enkel als een perturbatieve insertie op een SU(4)-symmetrische golffunctie mist shell-structuur bijdragen van enkele tot tien MeV in middelzware kernen. Het inbedden ervan in de propagatie genereert dynamisch de correcte shell-structuur.

• Resolutie van Benchmark Discrepanties:
  Het artikel adresseert recente kritieken (Ref. [1]) door conceptuele en technische fouten in de vergelijking van verschillende lattice-theorieën te identificeren:

  • Inequivalente Vergelijkingen: Het vergelijken van Hartree-Fock (HF) grenzen van een $a_t \to 0$ Hamiltonian met AFQMC resultaten van een eindige-$a_t$ transfer-matrix is ongeldig. De eindige-$a_t$ transfer-matrix definieert een distinct gereguleerde theorie met geïnduceerde interacties. Wanneer dezelfde gereguleerde theorie consistent wordt gebruikt, houdt het variationele principe stand (HF-energieën begrenzen geprojecteerde energieën).
  • Correlatie-energieën: Kleine correlatie-energieën in neutronenmaterie voor de "Essential Elements" (EE) interactie impliceren geen triviale golffuncties. Gedetailleerde analyse van irreducibele dichtheden onthult significante multimodale superfluiditeit (s-golf, p-golf, en quartet correlaties) ondanks kleine energieverschillen.
  • Nucleaire Verzadiging: Het afvlakken van de interactie-energie bij hoge dichtheden is geen lattice-packing artefact. Analytische afleidingen in continue ruimte met Galileïsche-invariante regulatoren tonen aan dat het afvlakken van de interactie-energie een natuurlijk gevolg is van momentum-ruimte regulatie. De verzadigingseigenschappen van NLEFT-interacties volgen de smearing-geïnduceerde regulator-schaal, niet de lattice-site-bezetting-schaal.
  • Continuüm-limiet: De claim dat lattice-interacties leiden tot systeeminstorting in de continuüm-limiet wordt weerlegd. Het schaalargument dat in de kritiek wordt gebruikt, berust op een incorrecte Born-niveau renormalisatie; correcte renormalisatie naar fysieke verstrooiings-observabelen levert een andere schaalgedrag op ($C_c \propto a$) die instorting voorkomt.

Betekenis en Claims
Het artikel beweert een solide fundament te hebben gelegd voor volledige foutenanalyse in NLEFT-simulaties met volledige chirale krachten. Door onzekerheden rigoureus te kwantificeren en aan te tonen dat sign-problem-free methoden sub-procent precisie kunnen bereiken voor zware kernen, valideert dit werk NLEFT als een betrouwbare benchmark voor many-body methoden.

De auteurs stellen dat de schijnbare discrepanties in recente literatuur (Ref. [1]) worden opgelost wanneer vergelijkingen worden gemaakt tussen equivalente gereguleerde theorieën met consistente parameters, volumes en randvoorwaarden. De LAT-OPT1 actie wordt gepresenteerd niet enkel als een fenomenologisch instrument, maar als een exact oplosbaar, niet-perturbatief benchmark-model dat essentiële shell-structuur dynamica vangt. Dit kader maakt de scheiding mogelijk tussen genuanteerde tekortkomingen in nucleaire krachten en many-body solver fouten, wat een pad biedt voor voorspellende studies van shell-evolutie, spectroscopie en dichte materie in middelzware en zware kernen.