Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

Technische Samenvatting: Full Configuration Interaction Quantum Monte Carlo voor Nauwkeurige Ab Initio Kernstructuurberekeningen

Probleemstelling
Ab initio kernwetenschap streeft ernaar om eindige kernen te beschrijven vanuit microscopische Hamiltoniaanse functies afgeleid van chirale effectieve veldentheorie. Hoewel veel-deeltjes expansiemethoden (bijv. Coupled Cluster, IMSRG, MBPT) berekeningen naar zwaardere kernen uitbreiden, vertrouwen zij op afkappingen die residuele veel-deeltjes onzekerheden introduceren die moeilijk te kwantificeren zijn. Daarentegen biedt de No-Core Shell Model (NCSM) directe toegang tot de volledige configuration interaction (FCI) oplossing, maar is deze beperkt door de exponentiële groei van de Hilbertruimte. Er is een kritieke behoefte aan stochastische solvers die in staat zijn de FCI-limiet te benaderen in grote modelruimtes zonder de veel-deeltjes Hamiltonian matrix expliciet op te slaan, met name voor systemen met sterke hogere-orde correlaties voortkomend uit vervorming, clustering of open-shell structuren.

Methodologie
Dit werk past het Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algoritme, oorspronkelijk ontwikkeld voor kwantumchemie, aan voor kernstructuurberekeningen. De methode maakt gebruik van een stochastische projectie-aanpak in imaginaire tijd om de Schrödinger-vergelijking op te lossen binnen een gekozen symmetriesector (deeltjesaantal, pariteit en impulsmomentprojectie).

Belangrijke algoritmische componenten beschreven zijn:

• Voortplanting in Imaginaire Tijd: De golffunctie wordt gerepresenteerd door een populatie van getekende 'walkers' verdeeld over Slater-determinanten. De voortplanting omvat stappen van spawning (geboorte), dood/klonen en annihilatie om de fermionische tekenstructuur stochastisch vast te stellen.
• Generatie van Excitaties: Om de efficiëntie te verbeteren ten opzichte van uniforme generatie, maken de auteurs gebruik van een vooraf berekende heat-bath generator die bij voorkeur excitaties sampleert met grote Hamiltonian matrix-elementen.
• Initiator Approximatie (i-FCIQMC): Om het tekenprobleem in grote Hilbertruimtes te overwinnen, wordt spawning vanuit zwak bezette determinanten beperkt. Dit introduceert een eindige-walker bias die systematisch wordt verminderd naarmate de walker-populatie toeneemt.
• Adaptive Shift (AS) Correctie: Een lokale verschuiving wordt toegewezen aan niet-initiator determinanten om de convergentie naar de FCI-limiet te versnellen, terwijl de exacte limiet behouden blijft wanneer het aantal walkers naar oneindig gaat.
• Observabele Estimators:
  • Projected/Trial Estimators: Gebruikt voor energieën, gebruikmakend van referentie-determinanten of compacte trial-golffuncties.
  • Pure Estimators: Essentieel voor niet-commuterende observabelen zoals radii. Deze worden berekend via Reduced Density Matrix (RDM) sampling met behulp van twee onafhankelijke walker-replica's om de bias te elimineren die inherent is aan kwadratische estimators.
• Geëxciteerde Toestanden: De methode wordt uitgebreid naar geëxciteerde toestanden via orthogonalized multi-state propagatie, waarbij meerdere walker-ensembles parallel worden voortgeplant en georthogonaliseerd bij elke stap.

Belangrijkste Bijdragen

1. Formalisme Adaptatie: Het artikel beschrijft de specifieke implementatie van FCIQMC voor kern-Hamiltonia, inclus inclusief de behandeling van drie-deeltjes krachten en de specifieke excitatie-generatieschema's die vereist zijn voor kerninteracties.
2. Validatie tegenover Deterministische FCI: In kleine modelruimtes ($e_{\text{max}}=2$) voor $^4$He, $^8$Be, $^{12}$C en $^{16}$O, laten de stochastische FCIQMC resultaten voor energieën en gekwadrateerde punt-proton radii ($R_p^2$) zien dat zij de exacte deterministische FCI-waarden reproduceren binnen statistische onzekerheden.
3. Systematische Foutcontrole: De auteurs laten zien dat de residuele eindige-walker bias in initiator-approximaties kan worden gecontroleerd en gekwantificeerd. Zij introduceren een strategie met behulp van machtswet-extrapolaties ($O(N_w) = O_\infty + \lambda N_w^{-\alpha}$) in grote modelruimtes ($e_{\text{max}}=10$). De exponent $\alpha$ wordt geconstraint door consistentie te eisen tussen de snel convergerende geprojecteerde energie en de langzamer convergerende pure energie.
4. Extensie naar Geëxciteerde Toestanden: Het artikel demonstreert de capaciteit om laagliggende spectra (specifiek voor $^6$Li) te berekenen door meerdere toestanden binnen dezelfde symmetriesector te resolveren zonder expliciete diagonalisatie.

Resultaten

• Benchmarking: In de $e_{\text{max}}=2$ ruimte komen FCIQMC energieën en radii overeen met FCI-benchmarks. De adaptive shift correctie vermindert de eindige-walker bias aanzienlijk en verzwakt de afhankelijkheid van de initiator drempel ($n_{\text{init}}$).
• Grote Modelruimtes ($e_{\text{max}}=10$):
  • Voor $^4$He convergeert de berekening direct met toenemend aantal walkers ($N_w$), overeenkomend met NCSM-resultaten.
  • Voor $^8$Be, $^{12}$C en $^{16}$O convergeert de geprojecteerde energie snel, terwijl RDM-gebaseerde pure estimators (energie en radii) een tragere convergentie vertonen. Door de energie-geconstreinde machtswet-extrapolatie toe te passen, verkrijgen de auteurs resultaten die correlaties vangen die verder gaan dan getrunceerde veel-deeltjes expansiemethoden (MBPT, CC, IMSRG).
  • De extrapolatiestrategie wordt gevalideerd door aan te tonen dat resultaten onafhankelijk zijn van de adaptive-shift offset parameter ($\tilde{\Delta}$) en consistent zijn over verschillende initiator drempels.
• Geëxciteerde Toestanden: Voor $^6$Li slaagt de methode erin de vijf laagste toestanden ($1^+, 3^+, 0^+, 2^+_1, 2^+_2$) te resolveren. Het excitatie-spectrum is correct geordend, zelfs bij lagere aantallen walkers, en convergeert sneller dan NCSM-resultaten met betrekking tot de veel-deeltjes afkapping ($N_{\text{max}}$).

Betekenis
Het artikel positioneert FCIQMC als een systematisch verbeterbare stochastische solver die in staat is de FCI-limiet te benaderen voor kernsystemen waar deterministische methoden falen door de grootte van de Hilbertruimte. Het dient twee primaire rollen:

1. Nauwkeurige Solver: Het biedt hoog-precieze benchmarks voor grondtoestandsenergieën en radii in grote modelruimtes, waarbij hoge-orde correlaties worden gevangen die getrunceerde expansies missen.
2. Benchmark Tool: Het biedt een middel om de residuele veel-deeltjes onzekerheden van andere getrunceerde veel-deeltjes expansiemethoden te beoordelen door een referentie te bieden die dicht bij de exacte oplossing ligt.

Het werk legt de algoritmische en numerieke fundering voor de natuurkundige resultaten gerapporteerd in een begeleidende Letter [80], en bevestigt dat FCIQMC zowel gesloten als open-shell kernen binnen hetzelfde kader kan behandelen en kan worden uitgebreid naar excited-state spectroscopie. De auteurs merken op dat toekomstig werk zal bestaan uit expliciete drie-deeltjes interacties, grotere enkelvoudige deeltjes bases, en toepassingen op collectieve en elektrowerkelijke observabelen.