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

Sintesi Tecnica: Full Configuration Interaction Quantum Monte Carlo per calcoli accurati della struttura nucleare ab initio

Enunciato del Problema
La teoria nucleare ab initio mira a descrivere i nuclei finiti a partire da Hamiltoniane microscopiche derivate dalla teoria del campo efficace chirale. Mentre i metodi di espansione many-body (ad es. Coupled Cluster, IMSRG, MBPT) estendono i calcoli a nuclei più pesanti, essi si basano su troncamenti che introducono incertezze many-body residue difficili da quantificare. Al contrario, il No-Core Shell Model (NCSM) fornisce un accesso diretto alla soluzione di Full Configuration Interaction (FCI), ma è limitato dalla crescita esponenziale dello spazio di Hilbert. Esiste una necessità critica di solver stocastici capaci di avvicinarsi al limite FCI in grandi spazi di modello senza dover memorizzare esplicitamente la matrice dell'Hamiltoniana many-body, particolarmente per sistemi con forti correlazioni di ordine superiore derivanti da deformazione, clustering o strutture a guscio aperto.

Metodologia
Questo lavoro adatta l'algoritmo Full Configuration Interaction Quantum Monte Carlo (FCIQMC), originariamente sviluppato per la chimica quantistica, ai calcoli della struttura nucleare. Il metodo impiega un approccio di proiezione stocastica nel tempo immaginario per risolvere l'equazione di Schrödinger all'interno di un settore di simmetria scelto (numero di particelle, parità e momento angolare).

Componenti algoritmiche chiave descritte includono:

• Propagazione nel Tempo Immaginario: La funzione d'onda è rappresentata da una popolazione di walker con segno distribuiti su determinanti di Slater. La propagazione comporta passaggi di generazione (spawning), morte/clonazione e annichilazione per stabilire stocasticamente la struttura del segno fermionico.
• Generazione delle Eccitazioni: Per migliorare l'efficienza rispetto alla generazione uniforme, gli autori utilizzano un generatore heat-bath pre-calcolato che campiona preferenzialmente le eccitazioni con grandi elementi di matrice dell'Hamiltoniana.
• Approssimazione Initiator (i-FCIQMC): Per superare il problema del segno in grandi spazi di Hilbert, la generazione di walker da determinanti debolmente occupati è limitata. Ciò introduce un bias di popolazione finita che viene sistematicamente ridotto all'aumentare del numero di walker.
• Correzione Adaptive Shift (AS): Uno spostamento locale viene assegnato ai determinanti non-initiator per accelerare la convergenza verso il limite FCI, preservando al contempo il limite esatto quando il numero di walker tende all'infinito.
• Estimatori di Osservabili:
  • Estimatori Proiettati/Trial: Utilizzati per le energie, impiegano determinanti di riferimento o funzioni d'onda trial compatte.
  • Estimatori Puri: Essenziali per osservabili non commutanti come i raggi, vengono calcolati tramite campionamento della Matrice di Densità Ridotta (RDM) utilizzando due repliche indipendenti di walker per eliminare il bias inerente agli estimatori quadratici.
• Stati Eccitati: Il metodo viene esteso agli stati eccitati tramite propagazione multi-stato ortogonalizzata, dove molteplici ensemble di walker vengono propagati in parallelo e ortogonalizzati ad ogni passo.

Contributi Chiave

1. Adattamento del Formalismo: Il documento dettaglia l'implementazione specifica di FCIQMC per Hamiltoniane nucleari, includendo la gestione delle forze a tre nucleoni e gli schemi specifici di generazione delle eccitazioni richiesti dalle interazioni nucleari.
2. Validazione contro l'FCI Deterministico: In piccoli spazi di modello ($e_{\text{max}}=2$) per $^4$He, $^8$Be, $^{12}$C e $^{16}$O, i risultati stocastici di FCIQMC per energie e raggi protonici quadrati ($R_p^2$) mostrano di riprodurre i valori esatti dell'FCI deterministico entro le incertezze statistiche.
3. Controllo Sistematico dell'Errore: Gli autori dimostrano che il bias residuo del numero finito di walker nell'approssimazione degli initiator può essere controllato e quantificato. Introducono una strategia utilizzando estrapolazioni di legge di potenza ($O(N_w) = O_\infty + \lambda N_w^{-\alpha}$) in grandi spazi di modello ($e_{\text{max}}=10$). L'esponente $\alpha$ è vincolato dal requisito di coerenza tra l'energia proiettata, a rapida convergenza, e l'energia pura, a convergenza più lenta.
4. Estensione agli Stati Eccitati: Il lavoro dimostra la capacità di calcolare spettri a bassa energia (specificamente per $^6$Li) risolvendo molteplici stati all'interno dello stesso settore di simmetria senza diagonalizzazione esplicita.

Risultati

• Benchmarking: Nello spazio $e_{\text{max}}=2$, le energie e i raggi di FCIQMC corrispondono ai benchmark FCI. La correzione Adaptive Shift riduce significativamente il bias del numero finito di walker e indebolisce la dipendenza dalla soglia dell'initiator ($n_{\text{init}}$).
• Grandi Spazi di Modello ($e_{\text{max}}=10$):
  • Per $^4$He, il calcolo converge direttamente con l'aumento del numero di walker ($N_w$), corrispondendo ai risultati NCSM.
  • Per $^8$Be, $^{12}$C e $^{16}$O, l'energia proiettata converge rapidamente, mentre gli estimatori puri basati su RDM (energia e raggi) mostrano una convergenza più lenta. Applicando la strategia di estrapolazione a legge di potenza vincolata dall'energia, gli autori ottengono risultati che catturano le correlazioni oltre i metodi di espansione many-body troncati (MBPT, CC, IMSRG).
  • La strategia di estrapolazione è validata mostrando che i risultati sono indipendenti dal parametro di offset dell'adaptive-shift ($\tilde{\Delta}$) e coerenti tra diversi livelli di soglia dell'initiator.
• Stati Eccitati: Per $^6$Li, il metodo risolve con successo i cinque stati più bassi ($1^+, 3^+, 0^+, 2^+_1, 2^+_2$). Lo spettro di eccitazione è correttamente ordinato anche a numeri di walker inferiori, convergendo più velocemente dei risultati NCSM rispetto al troncamento many-body ($N_{\text{max}}$).

Significatività
Il documento posiziona l'FCIQMC come un solver stocastico sistematicamente migliorabile capace di avvicinarsi al limite FCI per sistemi nucleari dove i metodi deterministici falliscono a causa della dimensione dello spazio di Hilbert. Esso svolge due ruoli primari:

1. Solver Accurato: Fornisce benchmark ad alta precisione per energie e raggi dello stato fondamentale in grandi spazi di modello, catturando le correlazioni di ordine superiore che i metodi di espansione troncati trascurano.
2. Strumento di Benchmark: Offre un mezzo per valutare le incertezze many-body residue di altri metodi di espansione many-body troncati fornendo un riferimento vicino alla soluzione esatta.

Il lavoro stabilisce la base algoritmica e numerica per i risultati fisici riportati in una lettera gemella [80], confermando che l'FCIQMC può trattare nuclei a guscio chiuso e a guscio aperto all'interno dello stesso framework ed estendersi alla spettroscopia degli stati eccitati. Gli autori osservano che il lavoro futuro comporterà l'uso esplicito di interazioni a tre nucleoni, basi a singola particella più ampie e applicazioni a osservabili collettive ed elettrodeboli.