$K \pi$ scattering as a step towards $B \to K^* \ell^+ \ell^-$ from Lattice QCD

Questo studio presenta i primi risultati esplorativi di un calcolo reticolare che combina la determinazione variazionale degli stati $K\pi$ in volume finito con il formalismo $1+J\to2$ e una strategia duale per i quark pesanti, al fine di ottenere input adronici precisi per il decadimento $B \to K^* \ell^+ \ell^-$, concentrandosi inizialmente sull'analisi cinematica ad alto $q^2$.

L'essenziale

Immagina di essere un detective che sta cercando di risolvere il caso più importante della fisica moderna: c'è qualcosa di nuovo e misterioso oltre le regole conosciute dell'universo?

Per trovare la risposta, i fisici osservano un "reato" molto raro che avviene nelle particelle subatomiche: un quark pesante (chiamato quark b) che si trasforma in un quark più leggero (quark s) emettendo due particelle di luce (elettroni o muoni). Questo processo è come un'auto che cambia colore mentre passa attraverso un tunnel, rilasciando due fari luminosi.

Il problema è che questo "tunnel" non è vuoto. È pieno di traffico caotico fatto di altre particelle (adroni) che si scontrano e si fondono. Per capire se il "cambio di colore" dell'auto è dovuto a una nuova legge fisica o semplicemente al caos del traffico, dobbiamo essere in grado di prevedere esattamente come si comportano queste particelle nel tunnel.

Ecco di cosa parla questo documento: è un aggiornamento sui lavori in corso di un gruppo di scienziati (tra cui ricercatori del CERN e dell'Università di Edimburgo) che stanno costruendo un "simulatore di traffico" ultra-preciso al computer per risolvere questo mistero.

Ecco i punti chiave spiegati con metafore semplici:

1. Il Problema: Il "Tunnel" Instabile

Nella fisica delle particelle, spesso trattiamo le particelle finali come se fossero solide e stabili, come dei mattoni. Ma in questo caso specifico, il prodotto finale è una particella chiamata $K^*$ (K-asterisco).
Immagina la $K^*$ non come un mattone, ma come una bolla di sapone. È instabile, vive per un istante brevissimo e poi scoppia in due pezzi (un kaone e un pione).
Fino a poco tempo fa, i computer facevano fatica a simulare bolle di sapone che esplodono mentre si muovono velocemente. I vecchi metodi trattavano la bolla come se fosse un sasso solido, il che va bene per oggetti stabili, ma è un errore grave per le bolle. Questo documento dice: "Abbiamo imparato a simulare le bolle di sapone che esplodono!".

2. La Sfida: Il Peso dell'Auto e la Dimensione del Tunnel

Per simulare questo fenomeno, devono usare un "tunnel" (chiamato volume finito nel linguaggio tecnico) abbastanza grande da non schiacciare le particelle, ma abbastanza piccolo da farci stare il computer.
C'è un altro problema: il quark pesante ($b$) è come un elefante in una stanza piena di topolini (i quark leggeri).

• Se il tunnel è troppo piccolo, l'elefante sbatte contro i muri (errori di calcolo).
• Se il tunnel è troppo grande, il computer impiega un'eternità a calcolare come si muove l'elefante.
• Inoltre, per simulare l'elefante in modo realistico, serve una griglia di calcolo molto fitta (come una rete da pesca con maglie piccolissime), altrimenti l'elefante sembra un'ombra sfocata.

3. La Soluzione: Due Strategie in Una

Gli scienziati hanno adottato un approccio geniale, come se avessero due diversi tipi di occhiali per guardare lo stesso oggetto:

1. Occhiali per Elefanti: Usano un metodo specifico calibrato direttamente sulla massa reale dell'elefante (il quark $b$).
2. Occhiali per Topi Giganti: Usano un altro metodo simulando prima un "topo gigante" (un quark pesante ma non quanto l'elefante, simile al quark charm) e poi usano la matematica per "stirare" i risultati fino a raggiungere la massa dell'elefante.
   Usando entrambi i metodi sullo stesso set di dati, possono incrociare le informazioni e ottenere un risultato molto più sicuro, come due testimoni che confermano la stessa storia.

4. La Tecnica Segreta: La "Distillazione"

Per vedere cosa succede dentro questo caos, usano una tecnica chiamata Distillazione.
Immagina di avere una stanza piena di fumo (i dati grezzi del computer). La distillazione è come un filtro speciale che lascia passare solo i "profumi" più importanti (le particelle che ci interessano) e blocca il resto. Questo permette di costruire un elenco di "ingredienti" (operatori) che possono creare sia la particella stabile che la bolla di sapone instabile, tutto in un unico esperimento. È come avere un unico set di ingredienti per cucinare sia una torta semplice che un soufflé complesso.

5. Dove Siamo Ora? (Lo Stato dell'Arte)

Il documento è un "aggiornamento di lavoro" (non il risultato finale).

• Cosa hanno fatto: Hanno costruito il primo "motore" del simulatore. Hanno verificato che le loro "bolle di sapone" ($K^*$) si comportano come previsto e che i dati che stanno raccogliendo sono puliti.
• Cosa manca: Hanno solo una piccola quantità di dati (come se avessero guardato il traffico per 5 minuti). Hanno bisogno di raccogliere dati per giorni e giorni (migliaia di configurazioni) per avere una statistica solida.
• L'obiettivo: Una volta raccolti tutti i dati, potranno calcolare con precisione come l'auto cambia colore (la particella decade) in diverse condizioni.

Perché è importante?

Se i loro calcoli precisi non coincidono con ciò che vedono gli esperimenti reali (come quelli al CERN o in Giappone), allora abbiamo scoperto una nuova fisica. Potrebbe esserci una particella o una forza sconosciuta che sta influenzando il traffico nel tunnel.

In sintesi: questo gruppo sta costruendo il simulatore di guida più preciso mai creato per capire se le leggi della natura sono quelle che pensiamo, o se c'è un "fantasma" che guida l'auto. Al momento, stanno solo testando se il motore del simulatore funziona e se le "bolle di sapone" esplodono correttamente.

Sommario tecnico

Titolo: Scattering $K\pi$ come passo verso $B \to K^*\ell^+\ell^-$ dalla QCD su reticolo

1. Il Problema e la Motivazione

Le decadimenti rari $b \to s \ell^+\ell^-$ rappresentano uno dei test più sensibili del Modello Standard (SM) e una finestra privilegiata per la ricerca di nuova fisica. In particolare, il canale fenomenologicamente cruciale $B \to K^* (\to K\pi) \ell^+\ell^-$ richiede un trattamento rigoroso degli effetti adronici.

• Sfida principale: Lo stato finale $K^*$ è una risonanza (non una particella stabile) che decade in $K\pi$. I calcoli precedenti sulla QCD su reticolo per stati finali vettoriali (come $K^*$ e $\phi$) sono spesso datati o si basano sull'approssimazione di larghezza stretta, trattando gli stati finali come stabili. Questa approssimazione è accettabile per il $\phi$ (molto stretto), ma inadeguata per il $K^*$, che è molto più largo.
• Complessità dinamica: Il calcolo deve gestire simultaneamente:
  1. La dinamica di un sistema a due adroni risonante ($K\pi$) in un volume finito.
  2. La dinamica pesante-leggera che coinvolge il quark $b$, richiedendo un controllo preciso degli errori di discretizzazione.
• Obiettivo: Fornire un input adronico non perturbativo, preciso e sistematicamente migliorabile per le previsioni del Modello Standard, superando le limitazioni delle approssimazioni attuali.

2. Metodologia

Il lavoro presenta un calcolo esplorativo che combina tecniche avanzate di QCD su reticolo:

• Formalismo di Volume Finito:
  • Utilizzo del formalismo Lüscher per lo scattering elastico $K\pi$ per determinare le energie e le fasi di scattering.
  • Applicazione del formalismo generale $1+J \to 2$ (analogia del fattore Lellouch-Lüscher) per collegare gli elementi di matrice in volume finito agli ampiezze fisiche in volume infinito per le transizioni elettrodeboli.
• Strategia dei Quark Pesanti (Dual Heavy-Quark):
  Per gestire il quark $b$ su un reticolo con spaziatura $a^{-1} \approx 2.7$ GeV, viene adottata una strategia duale sulla stessa ensemble:
  1. Azione RHQ (Relativistic Heavy Quark): Tarata direttamente sulla massa fisica del quark $b$.
  2. Quark Domain-Wall pesanti: Tre masse aggiuntive nel range $m_c \le m_h \lesssim 0.5 m_b$, che forniscono un "braccio di leva" per un'interpolazione esplorativa tra il charm e il bottom.
• Tecnica di Distillazione (Stochastic Distillation):
  • Tutti i correlatori (a due e tre punti) sono calcolati utilizzando la distillazione, che proietta i campi di quark su un sottospazio a rango ridotto (autovettori dell'operatore Laplaciano).
  • Setup Ibrido: Utilizzo di distillazione esatta su alcuni slice temporali e distillazione stocastica (con diluizione) per le linee quark sorgente-pozzo. Questo permette di gestire efficientemente operatori estesi e contrazioni multi-adroniche.
  • Inserimento di Corrente: Per le correnti elettrodeboli, viene mantenuta un'inserzione locale non "smearata" (non spalmata), necessaria per la precisione, a costo di una maggiore complessità I/O.
• Ensemble di Reticolo:
  • Utilizzo dell'ensemble F1M di RBC/UKQCD (fermioni a muro di dominio).
  • Spaziatura inversa: $a^{-1} \approx 2.7$ GeV.
  • Dimensioni: $48^3 \times 96$.
  • Masse dei quark leggeri: $M_\pi = 232$ MeV, $M_K = 510$ MeV (vicino al punto fisico, ma non esattamente al punto fisico).

3. Risultati Preliminari e Stato Attuale

Il documento presenta i primi risultati di uno studio esplorativo:

• Analisi dei Correlatori a Due Punti:
  • Costruzione di una matrice di correlazione $5 \times 5$ nel sistema di riferimento a riposo, utilizzando una base variazionale che include operatori quark-antiquark ($V$) e operatori a due adroni ($K_i\pi_i$).
  • I dati attuali provengono da 8 configurazioni di campo di gauge (con 24 sorgenti ciascuna, raggruppate in 3 bin per l'analisi preliminare).
  • Le incertezze statistiche sono attualmente sottostimate a causa del basso numero di configurazioni; l'obiettivo è raggiungere 30-60 configurazioni per un'analisi GEVP (Generalized Eigenvalue Problem) stabile.
  • I risultati mostrano una buona qualità del segnale e la capacità di risolvere i livelli energetici nella regione della risonanza $K^*$ (circa 960 MeV), ben separata dalla soglia $K\pi$ (circa 742 MeV).
• Accesso Cinematico:
  • Lo studio si concentra attualmente sulla regione ad alto $q^2$ (basso rinculo), dove il controllo reticolare è più forte.
  • Raggiungere $q^2 \approx 0$ (alto rinculo) richiederebbe momenti molto elevati non fattibili su questo reticolo specifico, oltre a entrare in regimi con soglie inelastiche multiple.

4. Contributi Chiave

1. Integrazione di Tecniche Avanzate: È uno dei primi tentativi di combinare il formalismo variazionale per stati risonanti $K\pi$ con il formalismo $1+J \to 2$ per decadimenti di mesoni pesanti.
2. Strategia Duale per il Quark $b$: L'uso simultaneo di un'azione RHQ tarata al fisico e di masse domain-wall intermedie offre un approccio robusto per controllare gli errori di discretizzazione e l'interpolazione di massa.
3. Versatilità del Dataset: L'uso della distillazione stocastica permette di calcolare una vasta gamma di transizioni elettrodeboli ($b \to s, d, u$ e $c \to s, d, u$) in stati finali risonanti all'interno di un unico dataset, massimizzando l'efficienza computazionale.
4. Gestione I/O e Scalabilità: Nonostante le sfide tecniche legate all'inserimento di correnti locali (che richiedono buffer di memoria temporanei dell'ordine di 500 TB), il setup dimostra la fattibilità di calcoli complessi su larga scala.

5. Significato e Prospettive Future

• Impatto sulla Fisica delle Particelle: Questo lavoro getta le basi per calcoli di precisione delle forme di fattore per $B \to K^*\ell^+\ell^-$, fondamentali per interpretare le attuali anomalie sperimentali e cercare nuova fisica.
• Passi Successivi:
  • Aumentare le statistiche fino a 40-60 configurazioni.
  • Stabilizzare l'estrazione degli stati in volume finito tramite GEVP.
  • Incorporare i risultati dello scattering $K\pi$ per mappare gli elementi di matrice in volume finito alle ampiezze fisiche.
  • Estendere il calcolo a regioni di $q^2$ più basse (richiedendo ensemble con spaziatura più grossolana o strategie combinate) e trattare gli effetti delle risonanze di charmonio nella regione intermedia di $q^2$.

In sintesi, questo documento segna l'avvio di un programma ambizioso per calcolare da primi principi i decadimenti rari dei mesoni B in stati finali risonanti, superando le limitazioni delle approssimazioni di larghezza stretta e fornendo un input teorico cruciale per la fisica delle alte energie.