Realization of all-to-all fermion propagator for the first principle high accuracy strong interaction prediction

Problematica

La Cromodinamica Quantistica su Reticolo (LQCD) è uno strumento primario per calcolare la struttura degli adroni da primi principi, come l'accoppiamento assiale vettoriale del nucleone ($g_A$). Tuttavia, raggiungere una precisione sub-percentuale richiede di affrontare due grandi sfide:

1. Costo Computazionale dei Propagatori Fermionici All-to-All: Il calcolo di funzioni di correlazione $N$-punto (dove $N \ge 4$), essenziali per le correzioni QED, la rottura della simmetria di isospin (ISB) e i loop di quark disconnessi, richiede propagatori fermionici all-to-all. Il calcolo esatto scala come $O(V^2)$ con il volume del reticolo $V$, rendendolo proibitivamente costoso.
2. Contaminazione da Stati Eccitati (ESC): L'estrazione dei matrici di stato fondamentale è complicata dalla contaminazione da stati eccitati (ad esempio, stati di tipo Roper, stati multi-particellari come $\pi N$). Sebbene il metodo della distillazione (proiezione sui modi a basso momento) sopprima efficacemente l'ESC per gli spettri degli adroni, esso scarta i modi ad alto momento essenziali per le inserzioni di correnti locali, limitando la sua utilità per le funzioni $N$-punto e gli elementi di matrice.

Le tecniche esistenti, come le sorgenti sequenziali, gli approcci stocastici con sostituzione dei modi bassi (LMS) e il diradamento del campo (field sparsening), offrono soluzioni parziali ma faticano a gestire simultaneamente funzioni $N$-punto ad alta precisione, diagrammi disconnessi e la soppressione dell'ESC senza incorrere in costi computazionali massicci.

Metodologia: L'Algoritmo di "Blending"

Gli autori propongono un nuovo metodo stocastico chiamato algoritmo di "blending" per approssimare efficientemente e senza bias i propagatori fermionici all-to-all.

• Costruzione dello Spazio di Blending: Il metodo costruisce uno "spazio di blending" per un tempo fisso combinando due sottospazi:
  1. Modi a Bassa Frequenza (LFM): Il sottospazio $L_1$ spaziato dai $N_e$ autovettori più bassi del Laplaciano spaziale invariante per gauge (simile alla distillazione). Questo cattura la fisica a lunga distanza e permette la costruzione di stati adronici spazialmente estesi.
  2. Modi ad Alta Frequenza (SHFM): Il complemento ortogonale $L_2$ è approssimato stocasticamente utilizzando $N_{st}$ vettori casuali ($\eta_j$) campionati uniformemente da $L_2$.
• Stima Unbiased (Non Distorta): L'operatore identità sullo spazio completo è approssimato da una somma pesata di proiezioni su questi vettori di base.
  • I proiettori dei modi bassi mantengono il loro peso esatto (1).
  • I proiettori stocastici dei modi alti sono assegnati un peso di $\dim(L_2)/N_{st}$ per garantire una stima non distorta dell'identità sul complemento.
  • Per operatori che coinvolgono prodotti di campi (ad esempio, nelle funzioni $N$-punto), viene applicato un tensore di peso $\Omega^{(2)}_{ij}$ per correggere il doppio conteggio e garantire che lo stimatore rimanga non distorto, particolarmente per i termini stocastici diagonali.
• Compressione: Il propagatore di quark all-to-all $S$ è compresso proiettandolo su questo spazio di blending. Ciò riduce i requisiti di memoria da $O((4N_c N_L^3 N_T)^2)$ a $O((4N_T(N_e + N_{st}))^2)$, consentendo il calcolo flessibile di funzioni di correlazione con arbitrarie combinazioni di operatori interpolanti e inserzioni di correnti.
• Controllo degli Stati Eccitati: Il metodo facilita l'uso dell'analisi del Problema del Valore Proprio Generalizzato (GEVP) con una base arricchita di interpolatori (inclusi operatori che coinvolgono la corrente come gli stati di scattering $N f_1$) senza costi di inversione aggiuntivi, poiché i propagatori compressi possono essere riutilizzati per molteplici combinazioni di operatori.

Contributi Chiave e Risultati

Gli autori dimostrano l'efficienza e l'accuratezza del metodo di blending attraverso diverse applicazioni su due ensemble di CLQCD: C24P29 ($m_\pi \approx 292$ MeV, $a \approx 0.105$ fm) e F64P13 ($m_\pi \approx 134$ MeV, $a \approx 0.077$ fm).

1. Validazione della Non Distorsione (Unbiasedness):

   • Il metodo è stato testato utilizzando l'elemento di matrice della corrente vettoriale conservata del pione ($\pi^+$). I risultati hanno mostrato che l'uso dei soli modi bassi ($N_{st}=0$) produceva un risultato distorto ($\approx 0.64$), mentre l'inclusione dei modi stocastici ($N_{st} \ge 50$) ha ripristinato il risultato al valore teorico atteso di 1.000(29).
   • Controlli di coerenza sono stati eseguiti per il fattore di forma elettrico del pione ($f_{\pi\pi}(Q^2)$) e il raggio di carica utilizzando sia funzioni di correlazione a 3 che a 4 punti. I risultati sono coerenti entro le incertezze, validando il calcolo della funzione a 4 punti.

2. Accoppiamento Assiale del Nucleone ($g_A$):

   • Iso-vettore ($g_A^{u-d}$): Utilizzando un fit congiunto a 3 stati attraverso vari conteggi di modi bassi ($\bar{N}_e$) e separazioni sorgente-sink, gli autori hanno ottenuto $g_A^{u-d} = 1.2337(84)$ alla massa fisica del pione. Questo risultato elimina l'ambiguità associata alla selezione di un $\bar{N}_e$ specifico e mostra un eccellente accordo con altri studi ad alta precisione.
   • Decomposizione Singoletto e di Flavor: Il metodo ha permesso il calcolo dei loop di quark disconnessi per la carica assiale singoletto ($g_A^{u+d+s}$) e i contributi di singolo flavor.
     • Utilizzando un'analisi GEVP con una base che include interpolatori $N$, $N f_1$ e $N f'_1$, hanno soppresso efficacemente la contaminazione da stati eccitati.
     • I risultati finali al punto fisico ($a=0.077$ fm) sono:
       • $g_A^u = 0.8408(86)$
       • $g_A^d = -0.3929(86)$
       • $g_A^s = -0.0381(57)$
       • $g_A^{u+d+s} = 0.410(20)$

3. Efficienza Computazionale:

   • Il metodo di blending riduce significativamente il numero di inversioni dell'operatore di Dirac richieste rispetto ai metodi tradizionali di sorgente sequenziale o stocastica.
   • Per l'isovettore $g_A$, il costo è stato ridotto di circa un ordine di grandezza rispetto a studi recenti ad alta precisione (es. ETM, RQCD, PNDME) che mirano a una precisione simile.
   • Per il singoletto $g_A$ (che coinvolge loop disconnessi), il metodo ha raggiunto una precisione del 4% con 0.70 milioni di inversioni, rispetto ai 1.37–34.1 milioni di inversioni in lavori precedenti comparabili. Questa efficienza deriva dalla capacità di riutilizzare il propagatore compresso per molteplici interpolatori e inserzioni di correnti senza nuove inversioni.

Significato e Rivendicazioni

Il paper sostiene che il metodo di blending fornisce un quadro robusto per i calcoli di LQCD ad alta precisione della struttura adronica, affrontando specificamente le doppie sfide del costo computazionale e della contaminazione da stati eccitati.

• Controllo Sistematico: Permettendo la costruzione di funzioni di correlazione con molteplici operatori interpolanti che condividono gli stessi numeri quantici, il metodo facilita l'estrazione robusta degli elementi di matrice dello stato fondamentale tramite analisi variazionale (GEVP), sopprimendo efficacementamente l'ESC.
• Scalabilità: Alla massa fisica del pione, il metodo rende le previsioni sub-percentuali per gli elementi di matrice del nucleone ottenibili con risorse computazionali significativamente ridotte.
• Versatilità: L'approccio è applicabile a varie funzioni $N$-punto, incluse le diagrammi disconnessi, le correzioni QED e gli effetti di ISB, che sono cruciali per testare il Modello Standard.
• Limitazioni: Gli autori notano che il costo computazionale rimane proibitivo per i fermioni Overlap, a causa del loro maggiore costo di inversione, suggerendo la necessità di ulteriori ottimizzazioni per una più ampia applicabilità.

In sintesi, il metodo di blending rappresenta un passo avanti significativo nel rendere fattibili i calcoli ad alta precisione di complessi osservabili adronici da primi principi, consentendo studi futuri di scattering multi-corpo e fisica dei partoni con incertezze sistematiche e statistiche ridotte.