Lattice Calculation of Short-Range Contributions to Neutrinoless Double-Beta Decay $\pi^-\to\pi^+ ee$ at Physical Pion Mass

Questo lavoro presenta un calcolo in QCD reticolare degli elementi di matrice $\pi^- \to \pi^+ ee$ alla massa fisica del pione utilizzando fermioni a muro di dominio, impiegando un metodo innovativo per sottrarre gli effetti "around-the-world" e una rinormalizzazione non perturbativa per ridurre le incertezze e riconciliare le discrepanze negli studi precedenti sui contributi a corto raggio del doppio decadimento beta senza neutrini.

L'essenziale

Il Quadro Generale: Risolvere un Mistero Cosmico

Immagina l'universo come un gigantesco puzzle e uno dei pezzi mancanti sia la natura dei neutrini. Queste sono particelle minuscole e fantasmatiche che attraversano tutto. Gli scienziati vogliono sapere: sono fermioni "Dirac" (come gli elettroni normali, dove una particella e la sua antiparticella sono distinte) o fermioni "Majorana" (dove una particella è la propria antiparticella)?

L'unico modo per risolvere questo mistero è osservare un evento molto raro chiamato decadimento doppio beta senza neutrini. È come guardare due atomi trasformarsi spontaneamente in atomi diversi ed espellere due elettroni, ma senza rilasciare alcun neutrino. Se vediamo che questo accade, dimostra che i neutrini sono le proprie antiparticelle.

Il Problema: Un Segnale Rumoroso

Per prevedere se questo evento avverrà, i fisici devono fare calcoli complessi. Scompongono il calcolo in due parti:

1. La Parte a Lunga Distanza: Come un sussurro che attraversa una stanza.
2. La Parte a Breve Distanza: Come una grida che avviene proprio accanto al tuo orecchio.

Questo documento si concentra sulla Parte a Breve Distanza. Nello specifico, calcolano come due pioni (particelle composte da quark) interagiscono per produrre due elettroni. Pensa a questo come misurare il "volume" di quella grida.

Il Conflitto: Due diversi gruppi di scienziati avevano precedentemente tentato di misurare questo "volume" utilizzando supercomputer (chiamati QCD su reticolo). Tuttavia, i loro risultati differivano di un fattore due. Era come se un team dicesse che la grida era a 60 decibel e l'altro dicesse che era a 120 decibel. Questo enorme disaccordo rendeva difficile fidarsi delle previsioni per il mistero dei neutrini.

La Soluzione: Un Nuovo Modo per Pulire i Dati

Gli autori di questo documento hanno deciso di condurre il proprio esperimento per mettere fine alla questione. Hanno utilizzato un massiccio supercomputer per simulare il mondo subatomico. Ma si sono trovati di fronte a un problema tecnico specifico: "Effetti Intorno al Mondo".

L'Analogia: Immagina di registrare una conversazione in una piccola stanza rimbombante con una parete circolare. Se batti le mani, il suono viaggia in avanti, colpisce il muro, si avvolge attorno alla stanza e torna da te da dietro. Nella simulazione al computer, la "stanza" è la griglia dello spaziotempo. Poiché la griglia è finita, le particelle possono viaggiare per tutto il percorso del loop e interferire con la misurazione, creando un confuso "eco" che rovina i dati.

L'Innovazione: I metodi precedenti cercavano di indovinare come cancellare questi echi. Questo team ha inventato un nuovo metodo di sottrazione.

• Invece di indovinare, hanno isolato il segnale dell'"eco" direttamente dai dati.
• Hanno calcolato esattamente quanto forte fosse l'eco e l'hanno sottratto dal segnale principale.
• Il Risultato: Il "rumore" è scomparso, lasciando un segnale pulito e stabile (un "plateau") di cui potevano fidarsi.

La Verifica: Controllare il Righello

Per assicurarsi che il loro nuovo metodo non fosse difettoso, hanno verificato il loro lavoro confrontandolo con uno standard noto. Hanno calcolato un valore specifico (chiamato "parametro di sacca") che era stato misurato da altri team in precedenza.

• Il loro risultato corrispondeva perfettamente allo standard affidabile.
• Quando hanno confrontato il loro risultato con il team che aveva la differenza di "fattore due", hanno scoperto che i loro numeri erano esattamente il doppio di quelli dell'altro team.
• La Conclusione: Si è scoperto che l'altro team aveva probabilmente utilizzato un "righello" leggermente diverso (convenzione di normalizzazione) per le loro misurazioni. Una volta tenuto conto di quella differenza, i punti dati si allineano effettivamente molto bene. Il metodo degli autori conferma che il loro calcolo è corretto e risolve la confusione.

Il Risultato Finale

Il team ha calcolato con successo il contributo a "breve raggio" al processo di decadimento doppio beta senza neutrini con una precisione molto superiore rispetto al passato.

• Hanno rimosso gli "echi" (effetti intorno al mondo) che stavano rovinando i dati precedenti.
• Hanno utilizzato due diverse "lenti" matematiche (schemi di rinormalizzazione) per garantire che la loro matematica fosse solida.
• Hanno fornito un numero definitivo e ad alta precisione che aiuta i fisici a prevedere se vedremo infine questo raro decadimento negli esperimenti reali.

In breve: Hanno costruito un microscopio migliore, hanno pulito il rumore statico e hanno confermato che il precedente disaccordo era solo una questione di utilizzo di diversi metri di misura. Ora, la comunità scientifica ha un numero affidabile per aiutare a risolvere il mistero della massa dei neutrini.

Sommario tecnico

1. Enunciazione del Problema

Il decadimento doppio beta senza neutrini ($0\nu\beta\beta$) è un processo critico per determinare se i neutrini siano fermioni di Majorana e per sondare la violazione del numero leptonico (LNV). Il tasso di decadimento dipende dagli elementi di matrice nucleare (NME). Nella Teoria di Campo Effettiva del Modello Standard (SMEFT), i contributi a corto raggio alla $0\nu\beta\beta$ derivano da operatori di dimensione 9. Questi contributi sono mediati dal processo $\pi^- \to \pi^+ ee$ tramite scambio di pioni tra nucleoni.

Il Problema Centrale: I precedenti calcoli reticolari QCD degli elementi di matrice a corto raggio per $\pi^- \to \pi^+ ee$ (in particolare i Ref. [30] e [31]) hanno riportato risultati significativamente discordanti. Per il parametro di bag $B_\pi^{MS}$ a $\mu=3$ GeV, il Ref. [30] ha riportato $\approx 0.42$, mentre il Ref. [31] ha riportato $\approx 0.20$—una discrepanza di circa un fattore due. Questa incoerenza ostacola previsioni teoriche affidabili per i tassi di decadimento $0\nu\beta\beta$. Inoltre, gli studi precedenti hanno sofferto di:

• Incertezze Sistemiche: Significativi effetti "intorno al mondo" (ATW) dovuti alla propagazione inversa dei pioni nei reticoli temporali periodici, che compromettono l'estrazione degli elementi di matrice dello stato fondamentale.
• Estrapolazione Chirale: Alcuni studi si sono basati su masse di pioni non fisiche, richiedendo un'estrapolazione al punto fisico, introducendo un'ulteriore dipendenza dal modello.

2. Metodologia

Gli autori hanno eseguito un calcolo QCD reticolare ad alta precisione utilizzando gli ensemble di Fermioni a Muro di Dominio (DWF) generati dalla collaborazione RBC/UKQCD.

A. Configurazione Reticolare

• Ensemble: Sono stati utilizzati ensemble $N_f = 2+1$, inclusi due con masse di pioni fisiche (48I e 64I) con diversi spazi reticolari ($a^{-1} \approx 1.73$ e $2.36$ GeV) per consentire un'estrapolazione diretta al continuo. Sono stati utilizzati anche ensemble con masse non fisiche per controlli incrociati.
• Operatori: Sono stati calcolati gli elementi di matrice per cinque operatori a quattro quark di dimensione 9 ($O_1, O_2, O_3, O'_1, O'_2$) rilevanti per la transizione $\pi^- \to \pi^+ ee$.

B. Mitigazione degli Effetti "Intorno al Mondo" (ATW)

Una delle principali innovazioni di questo lavoro è un nuovo metodo per gestire gli effetti ATW, che si manifestano quando i pioni propagano attraverso il confine temporale periodico, contaminando le funzioni di correlazione.

• Analisi: Gli autori hanno analizzato la struttura temporale delle funzioni di correlazione a tre punti, identificando i contributi dal diagramma desiderato (A) e dai diagrammi ATW (B, C, D).
• Nuovo Metodo di Sottrazione: Invece di affidarsi a metodi di rapporto (R1 o R2) che sopprimono solo parzialmente gli effetti ATW, gli autori hanno proposto una tecnica di sottrazione diretta.
  • Hanno identificato le regioni temporali in cui il contributo ATW da un singolo pione che attraversa il confine (Diagramma B) domina.
  • Hanno ricostruito il contributo ATW direttamente dai dati reticolari utilizzando separazioni temporali specifiche ($t_{\pi\pi}, t_\pi$).
  • Formula: L'elemento di matrice corretto è ottenuto sottraendo il termine ATW ricostruito:
    $$O^{(sub)}_i(t) = N_\pi^{-2} e^{2m_\pi t} \left[ C^i_3(t,t) - 2 C^i_3(t_{\pi\pi}, T - t_{\pi\pi} - t_\pi) e^{-m_\pi(T-t_\pi-2t_{\pi\pi})} \right]$$
  • Questo metodo ha ripristinato con successo plateau stabili nei grafici della massa efficace, riducendo significativamente gli errori sistematici rispetto ai precedenti metodi di rapporto.

C. Rinormalizzazione

• Schema: La rinormalizzazione non perturbativa è stata eseguita utilizzando il metodo RI/SMOM (Roma-Southampton) con cinematiche non eccezionali.
• Schemi: Sono stati impiegati due schemi per stimare gli errori di troncamento perturbativo: $(\gamma_\mu, \gamma_\mu)$ e $(\not{q}, \not{q})$.
• Accoppiamento: I risultati sono stati convertiti nello schema $\overline{MS}$ a $\mu = 3$ GeV utilizzando coefficienti di accoppiamento a un loop. La differenza tra i due schemi RI/SMOM è stata utilizzata per stimare l'incertezza sistematica dovuta al troncamento perturbativo.

3. Risultati Chiave

A. Elementi di Matrice Nudi

Lo studio ha estratto elementi di matrice nudi $\langle \pi^+ | O_i | \pi^- \rangle$ alla massa fisica del pione.

• Controllo Incrociato con Parametri di Bag: Gli autori hanno verificato la loro normalizzazione confrontando l'elemento di matrice di $O_3$ (relativo al parametro di bag del mesone neutro $B_K$) con i risultati precedenti revisionati dal FLAG. I risultati concordavano perfettamente con il Ref. [38] (che utilizzava la stessa convenzione di normalizzazione) ma erano approssimativamente due volte più grandi rispetto ai risultati nel Ref. [31].
• Discrepanza di Normalizzazione: Il rapporto tra i risultati degli autori e il Ref. [31] è costantemente $\approx 2.0$ per tutti e cinque gli operatori. Ciò suggerisce un errore globale di fattore di normalizzazione nel Ref. [31], piuttosto che una discrepanza fisica.

B. Elementi di Matrice Rinormalizzati e Parametro di Bag

I risultati finali nello schema $\overline{MS}$ ($\mu = 3$ GeV) sono:

• Parametro di Bag: $B_\pi^{MS}(3 \text{ GeV}) = 0.3769(24)_{stat}(76)_{PT}(1)_L$.
  • Questo valore è coerente con il Ref. [30] ($0.421$) entro le incertezze ma significativamente diverso dal Ref. [31] ($0.197$).
  • Il rapporto con il Ref. [31] è $1.91(18)$, confermando il problema di normalizzazione di un fattore 2.
• Estrapolazione al Continuo: Utilizzando gli ensemble a massa fisica (48I e 64I), gli autori hanno eseguito un'estrapolazione diretta al continuo, eliminando le incertezze associate all'estrapolazione chirale.
• Riduzione delle Incertezze: Le incertezze totali in questo lavoro sono significativamente ridotte rispetto agli studi precedenti grazie a:
  1. Calcolo diretto alla massa fisica del pione.
  2. Migliore precisione statistica.
  3. Rimozione efficace degli errori sistematici ATW.
  4. Stima robusta degli errori di troncamento perturbativo.

4. Significato e Conclusione

• Risoluzione della Discrepanza: Questo lavoro fornisce un controllo incrociato indipendente che supporta fortemente la normalizzazione utilizzata nel Ref. [30] e nella revisione FLAG, identificando al contempo un probabile errore di normalizzazione nel Ref. [31]. La discrepanza di un fattore 2 è risolta come una questione di convenzione di normalizzazione piuttosto che come un disaccordo fisico fondamentale.
• Avanzamento Metodologico: Il metodo di sottrazione proposto per gli effetti ATW è una tecnica robusta che può essere applicata ad altri calcoli reticolari che coinvolgono mesoni leggeri in volumi periodici, migliorando l'affidabilità dei futuri calcoli degli elementi di matrice $0\nu\beta\beta$.
• Impatto sulla Fisica: Fornendo elementi di matrice precisi, estrapolati al continuo e a masse di pioni fisiche, questo lavoro riduce le incertezze QCD non perturbative nella previsione teorica dei tassi di decadimento $0\nu\beta\beta$. Ciò è cruciale per interpretare i risultati sperimentali dai prossimi esperimenti di nuova generazione (ad es. LEGEND, nEXO, CUPID) per determinare la natura della massa del neutrino e la scala della nuova fisica.

In sintesi, questo lavoro stabilisce un nuovo punto di riferimento per i contributi a corto raggio al decadimento doppio beta senza neutrini, risolvendo le precedenti incoerenze attraverso una metodologia rigorosa, simulazioni a massa fisica e un trattamento innovativo degli artefatti di volume finito.