Theoretical framework for lattice QCD computations of $B\to K \ell^+ \ell^-$ and $\bar{B}_s\to \ell^+\ell^- \gamma$ decays rates, including contributions from charming penguin diagrams

Il paper sviluppa una strategia basata su metodi di densità spettrale per calcolare in QCD su reticolo i contributi complessi alle ampiezze di decadimento $B\to K\ell^+\ell^-$ e $\bar{B}_s\to\gamma\ell^+\ell^-$, inclusi i "pappagalli charm" e le divergenze ultraviolette, fornendo un quadro teorico per ridurre le incertezze fenomenologiche legate agli stati intermedi reali.

L'essenziale

Il Grande Puzzle delle Particelle: Come i Fisici "Vedono" l'Invisibile

Immagina di voler capire come funziona un orologio antico e complesso, ma non puoi aprirlo. Puoi solo osservare le lancette muoversi dall'esterno. Questo è un po' quello che fanno i fisici quando studiano il mondo subatomico: le particelle sono così piccole e veloci che non possiamo "vederle" direttamente mentre accadono.

In questo articolo, un gruppo di ricercatori italiani e britannici ha sviluppato un nuovo metodo per guardare dentro l'orologio, in particolare per studiare due tipi di "orologi" molto rari e delicati: il decadimento di una particella chiamata B in una K (con due leptoni, come elettroni o muoni) e il decadimento di una particella Bs in un fotone e due leptoni.

Ecco i concetti chiave, spiegati con analogie quotidiane.

1. Il Problema: Il "Fantasma" nel Mezzo

Quando una particella pesante (come la B) decade, a volte succede qualcosa di strano. Non va direttamente dal punto A al punto B. Prima, crea un "fantasma" temporaneo: una coppia di particelle (in questo caso, un quark charm e il suo antiparticella) che appare e scompare molto velocemente.

• L'analogia: Immagina di lanciare una palla da tennis contro un muro. Di solito, rimbalza e torna indietro. Ma a volte, la palla colpisce il muro, crea un'eco che rimbalza dentro la stanza per un attimo prima di uscire. Quell'eco è il "fantasma".
• Il problema: In fisica, questi "fantasmi" (chiamati pinguini affascinanti o charming penguins) sono pericolosi per i calcoli. Creano un'ambiguità: non sappiamo se il risultato è "reale" o "immaginario" (nel senso matematico). I metodi tradizionali di calcolo (come quelli usati finora) sono come tentare di misurare l'eco in una stanza buia: non riescono a distinguere il suono reale dal rumore di fondo. Questo crea un'incertezza enorme: i fisici non sanno se ciò che vedono è una nuova legge della natura (Nuova Fisica) o solo un'illusione causata da questi fantasmi.

2. La Soluzione: La "Fotografia Sfumata" (Metodo SFR)

Gli autori del paper hanno detto: "Non possiamo misurare l'eco istantaneamente, ma possiamo ricostruirla". Hanno usato una tecnica chiamata Metodo di Ricostruzione della Densità Spettrale (SFR).

• L'analogia: Immagina di voler capire la forma di un oggetto che si muove troppo velocemente per essere fotografato. Invece di cercare di congelare il tempo, scatti molte foto con un obiettivo sfocato (un "filtro" che rende tutto morbido). Ogni foto sfocata ti dà un'idea diversa. Se poi combini tutte queste foto sfocate con un algoritmo intelligente, riesci a ricostruire la forma precisa dell'oggetto, anche se non l'hai mai visto nitido.
• Come funziona: I ricercatori usano un supercomputer (la "Griglia" o Lattice) per calcolare come si comportano le particelle in uno spazio "finto" (lo spazio euclideo, che è come il nostro tempo ma girato di 90 gradi). Poi, usano un trucco matematico per "ripristinare" il tempo reale e vedere cosa succede con i fantasmi. È come se prendessero le foto sfocate e usassero un software per mettere a fuoco l'immagine finale, separando il segnale vero dal rumore.

3. Il "Rumore" e la Pulizia (Renormalizzazione)

C'è un altro ostacolo. Quando si fanno questi calcoli al computer, appaiono dei "difetti" matematici, come se il calcolo producesse numeri infiniti o senza senso quando due particelle si toccano troppo da vicino.

• L'analogia: È come se stessimo cercando di ascoltare una conversazione delicata in una stanza piena di gente che urla. Se due persone urlano troppo vicine, il microfono va in saturazione (diventa infinito).
• La soluzione: Gli autori spiegano come "pulire" il microfono. Hanno sviluppato un metodo per sottrarre matematicamente quel "rumore" infinito (chiamato divergenze di potenza) prima di ascoltare la conversazione vera. In pratica, dicono: "Sappiamo che il microfono fa questo rumore quando due persone sono vicine, quindi lo togliamo dal calcolo per vedere cosa dicono davvero".

4. Il Primo Test: Un Esperimento di Prova

Per vedere se il loro nuovo metodo funziona davvero, hanno fatto un "esperimento pilota".

• Hanno usato un computer per simulare il decadimento di una particella B in una K.
• Hanno usato una massa per la particella B leggermente più leggera di quella reale (per rendere il calcolo più veloce e meno costoso), ma hanno mantenuto le altre particelle con la loro massa reale.
• Il risultato: È stato un successo! Hanno potuto vedere i "fantasmi" (i pinguini affascinanti) e misurare quanto contribuiscono al decadimento. Hanno confrontato i loro risultati con un modello teorico semplice (come una mappa approssimativa) e hanno visto che i loro dati "veri" corrispondevano bene alla mappa, confermando che il metodo funziona.

Perché è Importante?

Perché tutto questo?
Perché oggi c'è un mistero nella fisica: alcuni esperimenti mostrano che le particelle si comportano in modo leggermente diverso da come prevede la teoria standard (il Modello Standard). Alcuni pensano che ci sia una "Nuova Fisica" (particelle o forze sconosciute). Altri pensano che sia solo un errore di calcolo dovuto a questi "fantasmi" che non siamo riusciti a gestire bene.

Questo articolo ci dice: "Non dobbiamo più indovinare!".
Ora abbiamo la mappa e gli strumenti per calcolare esattamente quanto pesano questi "fantasmi". Se dopo averli calcolati con precisione i dati sperimentali non combaciano ancora con la teoria, allora sì, avremo scoperto una Nuova Fisica. Se invece combaciano, allora il mistero era solo un calcolo impreciso.

In Sintesi

Gli autori hanno inventato un nuovo modo di "fotografare" eventi quantistici complessi che prima erano impossibili da calcolare con precisione. Hanno dimostrato che, usando i supercomputer e un po' di magia matematica (i filtri sfocati e la pulizia del rumore), possiamo finalmente vedere chiaramente cosa succede quando le particelle decadono. Questo ci avvicina molto di più a scoprire se l'universo nasconde segreti oltre le nostre conoscenze attuali.

Sommario tecnico

1. Il Problema: Calcolo delle Ampiezze di Decadimento FCNC

Il lavoro affronta la sfida teorica di calcolare con precisione le ampiezze di decadimento per i processi a Corrente Neutra che Cambia Sapore (FCNC) nel Modello Standard, in particolare:

• $B \to K \ell^+ \ell^-$
• $\bar{B}_s \to \gamma \ell^+ \ell^-$

Questi decadimenti sono soppresi nel Modello Standard e rappresentano finestre sensibili per la ricerca di Nuova Fisica (NP). Tuttavia, il confronto tra dati sperimentali (es. LHCb) e previsioni teoriche è limitato dall'incertezza nella valutazione degli elementi di matrice adronici.
In particolare, le contribuzioni complesse (con parte immaginaria) derivanti da stati intermedi "on-shell" (che possono propagarsi tra l'operatore debole e la corrente elettromagnetica) non possono essere calcolate con le tecniche standard della QCD su reticolo. Queste includono:

• I "Charming Penguins": Diagrammi contenenti loop di quark charm ($c\bar{c}$) che possono propagarsi a lunghe distanze, specialmente vicino alle risonanze del charmonio ($J/\psi$, $\psi(2S)$).
• L'operatore chromomagnetico ($O_8$).

Attualmente, queste contribuzioni sono stimate tramite modelli fenomenologici (es. Dominio dei Vettori di Vettore), introducendo incertezze significative che ostacolano test precisi del Modello Standard.

2. Metodologia: Il Framework SFR/HLT

Gli autori sviluppano una strategia completa per calcolare queste ampiezze utilizzando la QCD su reticolo combinata con metodi di densità spettrale.

A. Rappresentazione Spettrale e Ordinamento Temporale

Il cuore del metodo risiede nel trattamento delle funzioni di correlazione euclidee per recuperare le ampiezze di Minkowski complesse.

• L'ampiezza è espressa come integrale di una densità spettrale $\rho(E)$ su uno stato intermedio.
• A causa della presenza di stati intermedi on-shell con energie inferiori allo stato esterno, l'integrale presenta poli reali sull'asse reale, rendendo impossibile il calcolo diretto in spazio euclideo senza tecniche speciali.
• Il metodo SFR (Spectral Function Reconstruction), proposto in lavori precedenti [24], permette di ricostruire sia la parte reale che quella immaginaria dell'ampiezza.
• Viene utilizzata la tecnica HLT (Hansen-Lupo-Tantalo) [23] per approssimare il kernel di smearing (che regolarizza il polo) tramite una somma di esponenziali, permettendo il calcolo tramite funzioni di correlazione euclidee.

B. Generalizzazione a Operatori Multi-locali

Il paper generalizza il metodo SFR per gestire:

1. Decadimenti $B \to K \ell^+ \ell^-$: Richiedono la matrice di due operatori locali (corrente elettromagnetica e operatore debole $O_{1,2}^{(c)}$). Vengono analizzati i due ordinamenti temporali ($t<0$ e $t>0$); solo $t>0$ contribuisce alla parte immaginaria e richiede SFR.
2. Decadimenti $\bar{B}_s \to \gamma \ell^+ \ell^-$: Richiedono la matrice di tre operatori locali (due correnti elettromagnetiche e l'operatore debole). Il framework viene esteso per gestire sei possibili ordinamenti temporali, tre dei quali contengono poli e richiedono l'uso di densità spettrali generalizzate (integrale su due variabili di energia).

C. Rinormalizzazione e Sottrazione delle Divergenze

Un aspetto cruciale trattato è la gestione delle divergenze ultraviolette (UV):

• Termini di Contatto: Quando la corrente elettromagnetica e l'operatore debole si avvicinano, sorgono divergenze UV. Sebbene la conservazione della corrente elettromagnetica riduca la divergenza da quadratica a logaritmica nel risultato totale, ogni singolo ordinamento temporale è quadrativamente divergente.
• Proposta di Sottrazione: Gli autori propongono di separare le divergenze introducendo identità algebriche nel kernel di integrazione (Eq. 75) per isolare i termini divergenti (che possono essere calcolati con metodi standard) dalla parte fisica a lunga distanza (che richiede SFR).
• Divergenze di Potenza (Power Divergences): Gli operatori $O_{1,2}^{(c)}$ si mescolano con operatori di dimensione inferiore (densità scalari e pseudoscalari) con coefficienti che divergono come potenze inverse del passo reticolare ($1/a^2$, $1/a$).
  • Viene mostrata una procedura di sottrazione non perturbativa basata sulle simmetrie della QCD con fermioni twisted-mass (simmetrie CS, $S_3$, $R_5^{sp}$).
  • Si utilizza un modello ausiliario con quark valenza replicati per eliminare le divergenze di potenza imponendo condizioni di rinormalizzazione su elementi di matrice specifici (es. $\langle 0 | \tilde{O} | B \rangle = 0$).

3. Risultati Chiave e Studio Esplorativo

Gli autori hanno eseguito un calcolo numerico esplorativo per testare la fattibilità del metodo:

• Setup: Un singolo ensemble di gauge ETMC ($N_f = 2+1+1$, fermioni twisted-mass Wilson-Clover), passo reticolare $a \approx 0.08$ fm.
• Configurazione: Massa del quark $b$ più leggera del fisico ($m_b = 2m_c$) per ridurre gli artefatti di reticolo, masse degli altri quark fisiche.
• Processo Studiato: Il diagramma "charming penguin" in cui il fotone virtuale è emesso dal loop di charm (Fig. 1a).
• Analisi:
  • Calcolo delle funzioni di correlazione euclidee per $B \to K \ell^+ \ell^-$ e un processo ausiliario $B_s \to \eta_{ss'} \ell^+ \ell^-$ (dove $\eta_{ss'}$ è un mesone fittizio composto da due quark strange) per ridurre il rumore statistico.
  • Applicazione del metodo HLT per ricostruire l'ampiezza "smearata" $\tilde{H}(\epsilon)$ per diversi valori del parametro di smearing $\epsilon$.
  • Confronto con Modelli: I risultati sono stati confrontati con un modello basato sull'approssimazione di saturazione del vuoto (VSA) e distribuzioni di Breit-Wigner per le risonanze.
• Risultati Numerici:
  • Le componenti reale e immaginaria dell'ampiezza sottratta sono state determinate con successo.
  • Si osserva una struttura risonante chiara (picchi a $J/\psi$ e $\psi(2S)$) nella parte immaginaria, coerente con le aspettative fisiche.
  • Per l'operatore $O_2$, si osserva una soppressione significativa rispetto a $O_1$ (maggiore di quanto previsto dal VSA), dovuta alla cancellazione parziale tra componenti vettoriali e assiali.
  • Sfida dell'Estrapolazione: La precisione attuale non permette ancora di raggiungere pienamente il regime asintotico ($\epsilon \to 0$) necessario per un'estrapolazione polinomiale semplice vicino alle risonanze strette. Tuttavia, i risultati sono incoraggianti e suggeriscono che con statistiche maggiori e tecniche di riduzione del rumore, l'estrapolazione è fattibile.

4. Contributi e Significatività

• Framework Teorico Completo: Il paper fornisce la prima descrizione dettagliata di come calcolare rigorosamente le contribuzioni dei "charming penguins" e dell'operatore $O_8$ in QCD su reticolo, risolvendo il problema della continuazione analitica per stati intermedi on-shell.
• Risoluzione delle Divergenze: Offre una soluzione sistematica e non perturbativa per la sottrazione delle divergenze di potenza e dei termini di contatto, un prerequisito essenziale per qualsiasi calcolo di precisione.
• Generalizzabilità: Il metodo è applicabile a una vasta classe di decadimenti semileptonici e radiativi con stati finali stabili o instabili, aprendo la strada a calcoli di precisione per testare la Nuova Fisica.
• Impatto Fenomenologico: Permetterà di ridurre drasticamente le incertezze teoriche nelle previsioni per $B \to K \ell^+ \ell^-$ e $\bar{B}_s \to \gamma \ell^+ \ell^-$, permettendo di distinguere tra effetti di QCD non perturbativa (come i penguin charm) e segnali genuini di Nuova Fisica nelle anomalie osservate finora.

In sintesi, questo lavoro getta le basi per calcoli "ab-initio" di precisione delle ampiezze di decadimento FCNC, superando uno dei principali colli di bottiglia teorici nella fisica dei sapori moderna.