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

Technische Samenvatting: Inclusieve $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ vervalprocessen vanuit Lattice QCD

Probleemstelling
De bepaling van het Cabibbo–Kobayashi–Maskawa (CKM) matrixelement $|V_{cb}|$ is een centrale uitdaging in de flavor-fysica, die momenteel wordt gekenmerkt door een significante spanning tussen exclusieve en inclusieve experimentele bepalingen. Hoewel inclusieve semileptische vervalprocessen van $B$-mesonen een theoretisch robuust pad bieden via de Operator Product Expansion (OPE), werd een eerste-principes lattice QCD-berekening van de inclusieve verval-snelheid $\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}]$ historisch gezien als ontoegankelijk beschouwd vanwege het "twee-schaalprobleem" (het simultaan voldoen aan $a m_b \ll 1$ en $\Lambda_{QCD} L \gg 1$) en de snelle verslechtering van de signaal-ruisverhouding van vier-punts correlatoren naarmate de zware quarkmassa toeneemt.

Methodologie en Computationele Strategie
De auteurs presenteren een nieuwe strategie om deze inclusieve vervalsnelheid vanuit eerste principes te berekenen, gebruikmakend van $N_f = 2+1+1$ Extended Twisted Mass (ETMC) gauge-ensembles met fysieke lichte, strange en charm quarkmassa's. De methodologie rust op drie hoofdpijlers:

1. Massa-interpolatie via OPE: In plaats van direct te simuleren bij de fysieke $b$-quarkmassa (wat computationeel onhaalbaar is vanwege cutoff-effecten), simuleren de auteurs bij lichtere-dan-fysische zware quarkmassa's ($M_I < M_{\bar{B}_s}$). Zij behouden een vaste ratio tussen de initiële en finale mesonmassa's ($r_{IF} = M_F/M_I = M_{D_s}/M_{\bar{B}_s}$) om het faseruimte-volume constant te houden. Het fysieke resultaat wordt verkregen door deze niet-perturbatieve lattice-datapunten te interpoleren naar het fysieke punt met behulp van OPE-voorspellingen, die exact worden in het statische limiet ($m_b \to \infty$). Dit transformeert het extrapolatieprobleem in een interpolatieprobleem.

2. Nieuwe Numerieke Strategie voor Vier-punts Correlatoren: Een cruciale algoritmische innovatie is de introductie van een nieuwe strategie voor het berekenen van de vereiste vier-punts Euclidische correlatoren. De auteurs observeren dat de "oude strategie" (het fixeren van de tijdsafstanden tussen bron- en sink-operatoren) leidt tot een catastrofaal verlies van de signaal-ruisverhouding voor zware mesonen. De nieuwe strategie fixeert de tijdsafstanden tussen de zwakke stroom-inserties en de bron/sink-operatoren ($t_{snk} - t = 0 - t_{src}$) terwijl de tijd van de stroom-insertie $t$ wordt gevarieerd. Deze aanpak, hoewel computationeel duurder (vereist meer inversies), overwint succesvol het signaal-ruisprobleem, waardoor correlatoren op grote zware mesonmassa's kunnen worden opgelost waar de oude methode faalde.

3. Spectrale Reconstructie en Analyse: De vervalsnelheid wordt geëxtraheerd uit de geamputeerde Euclidische correlatoren met behulp van het Hansen–Lupo–Tantalo (HLT) algoritme. De auteurs voeren een rigoureuze analyse uit van systematische onzekerheden door de volgorde van limieten te variëren (continuüm $a \to 0$, smearing $\sigma \to 0$, en numerieke integratie) en twee verschillende interpolatiebenaderingen naar het fysieke punt toe te passen:

   • Fitting-functie Benadering: Gebruikmakend van een rationale functie voor de niet-perturbatieve bijdrage gecombineerd met perturbatieve OPE-termen.
   • Functionele Ruimte Benadering: Gebruikmakend van Gaussian Processes met een prior-distributie die verankerd is door het perturbatieve OPE-resultaat.

Belangrijkste Resultaten
De berekening werd uitgevoerd op drie ETMC-ensembles (B64, C80, D96) met lattice-afstanden $a \approx 0.08, 0.07, 0.06$ fm en ruimtelijke volumes $L \approx 5.5$ fm. De auteurs verkregen een voorlopige bepaling van de inclusieve vervalsnelheid:
$$\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}] = 2.47(13)(10)[17] \times 10^{-11} \text{ GeV}$$
waarbij de fouten respectievelijk statistisch, systematisch (uit analyse-stappen) en totaal zijn. De totale relatieve fout is ongeveer 7%.

Gebruikmakend van dit resultaat en door de $X_{\bar{s}u}$ bijdrage te verwaarlozen (gerechtvaardigd door de $|V_{ub}|^2/|V_{cb}|^2$ onderdrukking), leiden de auteurs een waarde voor $|V_{cb}|$ af:
$$|V_{cb}|_{\text{latt-incl}} = 40.9(1.2)_{\text{latt}}(1.7)_{\text{exp}}[2.1]_{\text{tot}} \times 10^{-3}$$
Deze waarde is consistent met zowel de exclusieve bepaling ($39.62(47) \times 10^{-3}$) als de OPE-gebaseerde inclusieve bepaling ($41.64(47) \times 10^{-3}$), hoewel de huidige statistische precisie onvoldoende is om de bestaande spanning tussen deze methoden op te lossen.

Betekenis en Claims
Het artikel claimt enkele significante bijdragen aan het vakgebied:

• Algoritmische Doorbraak: De nieuwe numerieke strategie voor vier-punts correlatoren lost het langdurige signaal-ruisprobleem voor zware mesonen op, wat de studie van inclusieve $B$-verval op de lattice mogelijk maakt.
• Methodologische Validatie: Dit werk demonstreert dat inclusieve semileptische vervalsnelheden vanuit eerste principes met lattice QCD kunnen worden berekend op een fenomenologisch relevant nauwkeurigheidsniveau door niet-perturbatieve data te combineren met OPE-restricties.
• Onafhankelijke Verificatie: Door een andere lattice-discretisatie (Twisted Mass) en analyse-techniek (HLT algoritme) te gebruiken vergeleken met gelijktijdige inspanningen (bijv. Ref. [15]), biedt dit werk een volledig onafhankelijk eerste-principes resultaat voor $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$.
• Toekomstperspectief: De auteurs benadrukken dat de huidige 7% fout wordt gedomineerd door het beperkte aantal lattice-ensembles en het gebrek aan een data-gedreven finite-volume studie. Zij stellen dat met de toevoeging van fijnere lattice-afstanden ($a \approx 0.05$ fm) en grotere volumes, de totale fout aanzienlijk kan worden verminderd, wat potentieel een definitieve oplossing kan bieden voor de $|V_{cb}|$ puzzel.

De auteurs hanteren een bescheiden toon met betrekking tot het huidige resultaat, waarbij zij opmerken dat het een "eerste fysiek resultaat" is en dat de systematische onzekerheden geassocieerd met finite-volume effecten en de continuümlimiet verdere verfijning vereisen in toekomstige publicaties.