Inclusive $\bar B_s\mapsto X_{\bar sc} \ell \bar \nu$ decays from lattice QCD: computational strategy and a first physical result

Sintesi Tecnica: Decimenti inclusivi $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ da QCD su reticolo

Definizione del Problema
La determinazione dell'elemento della matrice di Cabibbo–Kobayashi–Maska (CKM) $|V_{cb}|$ rappresenta una sfida centrale nella fisica delle particelle di sapore, attualmente caratterizzata da una significativa tensione tra le determinazioni sperimentali esclusive e inclusive. Sebbene i decadimenti semileptonici inclusivi dei mesoni $B$ offrano un percorso teoricamente robusto tramite l'Espansione dei Prodotti Operatori (OPE), il calcolo da primi principi del tasso di decadimento inclusivo $\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}]$ è stato storicamente considerato inaccessibile a causa del "problema delle due scale" (soddisfare simultaneamente $a m_b \ll 1$ e $\Lambda_{QCD} L \gg 1$) e del rapido deterioramento del rapporto segnale-rumore nei correlatori a quattro punti all'aumentare della massa del quark pesante.

Metodologia e Strategia Computazionale
Gli autori presentano una nuova strategia per calcolare questo tasso di decadimento inclusivo da primi principi, utilizzando ensemble di gauge $N_f = 2+1+1$ Extended Twisted Mass (ETMC) con masse dei quark leggeri, strani e charm fisiche. La metodologia si basa su tre pilastri principali:

1. Interpolazione di Massa via OPE: Invece di simulare direttamente alla massa fisica del quark $b$ (computazionalmente proibitivo a causa degli effetti di cutoff), gli autori simulano a masse di quark pesanti inferiori alla massa fisica ($M_I < M_{\bar{B}_s}$). Mantengono un rapporto fisso tra le masse del mesone iniziale e finale ($r_{IF} = M_F/M_I = M_{D_s}/M_{\bar{B}_s}$) per mantenere costante il volume dello spazio delle fasi. Il risultato fisico viene ottenuto interpolando questi punti di dati non perturbativi sul reticolo al punto fisico utilizzando le predizioni OPE, che diventano esatte nel limite statico ($m_b \to \infty$). Ciò trasforma il problema dell'estrapolazione in un problema di interpolazione.

2. Nuova Strategia Numerica per i Correlatori a Quattro Punti: Un'innovazione algoritmica critica è l'introduzione di una nuova strategia per calcolare i richiesti correlatori euclidei a quattro punti. Gli autori osservano che la "vecchia strategia" (fissare le separazioni temporali tra gli operatori sorgente e pozzo) porta a una catastrofica perdita del rapporto segnale-rumore per i mesoni pesanti. La nuova strategia fissa le separazioni temporali tra le inserzioni della corrente debole e gli operatori sorgente/pozzo ($t_{snk} - t = 0 - t_{src}$) variando il tempo di inserzione della corrente $t$. Questo approccio, sebbene computazionalmente più costoso (richiede più inversioni), supera con successo il problema del segnale-rumore, permettendo la risoluzione dei correlatori a grandi masse di mesoni pesanti dove il vecchio metodo falliva.

3. Ricostruzione Spettrale e Analisi: Il tasso di decadimento viene estratto dai correlatori euclidei amputati utilizzando l'algoritmo di Hansen–Lupo–Tantalo (HLT). Gli autori eseguono un'analisi rigorosa delle incertezze sistematiche variando l'ordine dei limiti (continuum $a \to 0$, smearing $\sigma \to 0$ e integrazione numerica) e impiegando due approcci di interpolazione distinti al punto fisico:

   • Approccio della Funzione di Fitting: Utilizzando una funzione razionale per il contributo non perturbativo combinata con termini OPE perturbativi.
   • Approccio nello Spazio Funzionale: Utilizzando Processi Gaussiani con una distribuzione a priori ancorata al risultato OPE perturbativo.

Risultati Chiave
Il calcolo è stato eseguito su tre ensemble ETMC (B64, C80, D96) con spaziature di reticolo $a \approx 0.08, 0.07, 0.06$ fm e volumi spaziali $L \approx 5.5$ fm. Gli autori hanno ottenuto una determinazione preliminare del tasso di decadimento inclusivo:
$$\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}] = 2.47(13)(10)[17] \times 10^{-11} \text{ GeV}$$
dove gli errori sono statistici, sistematici (dalle fasi di analisi) e totali, rispettivamente. L'errore relativo totale è circa del 7%.

Utilizzando questo risultato e trascurando il contributo $X_{\bar{s}u}$ (giustificato dalla soppressione $|V_{ub}|^2/|V_{cb}|^2$), gli autori derivano un valore per $|V_{cb}|$:
$$|V_{cb}|_{\text{latt-incl}} = 40.9(1.2)_{\text{latt}}(1.7)_{\text{exp}}[2.1]_{\text{tot}} \times 10^{-3}$$
Questo valore è coerente sia con la determinazione esclusiva ($39.62(47) \times 10^{-3}$) che con la determinazione inclusiva basata su OPE ($41.64(47) \times 10^{-3}$), sebbene l'attuale precisione statistica sia insufficiente per risolvere la tensione esistente tra questi metodi.

Significatività e Rivendicazioni
Il lavoro rivendica diversi contributi significativi nel campo:

• Svolta Algoritmica: La nuova strategia numerica per i correlatori a quattro punti risolve il problema di segnale-rumore di lunga data per i mesoni pesanti, consentendo lo studio dei decadimenti inclusivi di $B$ sul reticolo.
• Validazione Metodologica: Il lavoro dimoste che i tassi di decadimento semileptonici inclusivi possono essere calcolati con la QCD su reticolo da primi principi a un livello di accuratezza fenomenologicamente rilevante, combinando dati non perturbativi con vincoli OPE.
• Verifica Indipendente: Utilizzando una diversa discretizzazione del reticolo (Twisted Mass) e una tecnica di analisi differente (algoritmo HLT) rispetto agli sforzi concorrenti (es. Rif. [15]), questo lavoro fornisce un risultato totalmente indipendente da primi principi per $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$.
• Prospettive Future: Gli autori sottolineano che l'attuale errore del 7% è dominato dal numero limitato di ensemble di reticolo e dalla mancanza di uno studio del volume finito guidato dai dati. Essi affermano che, con l'aggiunta di spaziature di reticolo più fini ($a \approx 0.05$ fm) e volumi maggiori, l'errore totale può essere ridotto significativamente, permettendo potenzialmente una risoluzione definitiva del puzzle di $|V_{cb}|$.

Gli autori mantengono un tono modesto riguardo al risultato attuale, notando che si tratta di un "primo risultato fisico" e che le incertezze sistematiche associate agli effetti di volume finito e al limite del continuo richiedono un ulteriore perfezionamento nelle pubblicazioni future.