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

Resumen Técnico: Decaimientos inclusivos $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ desde QCD en el retículo

Planteamiento del Problema
La determinación del elemento de la matriz de Cabibbo–Kobayashi–Maskano (CKM) $|V_{cb}|$ es un desafío central en la física de sabor, actualmente caracterizado por una tensión significativa entre las determinaciones experimentales exclusivas e inclusivas. Mientras que los decaimientos semileptónicos inclusivos de mesones $B$ ofrecen una vía teóricamente robusta mediante la Expansión de Productos Operadores (OPE), el cálculo de primeros principios mediante QCD en el retículo (lattice QCD) de la tasa de decaimiento inclusivo $\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}]$ ha sido históricamente considerado inaccesible debido al "problema de las dos escalas" (satisfacer simultáneamente $a m_b \ll 1$ y $\Lambda_{QCD} L \gg 1$) y a la rápida degradación de la relación señal-ruido en los correladores de cuatro puntos a medida que aumenta la masa del quark pesado.

Metodología y Estrategia Computacional
Los autores presentan una estrategia novedosa para calcular esta tasa de decaimiento inclusivo desde primeros principios, utilizando conjuntos de medida (ensembles) de Masa Torcida Extendida (ETMC) con $N_f = 2+1+1$ con masas de quarks ligeros, extraños y de charm físicas. La metodología se basa en tres pilares principales:

1. Interpolación de Masas mediante OPE: En lugar de simular directamente en la masa física del quark $b$ (lo cual es computacionalmente prohibitivo debido a los efectos de corte), los autores simulan en masas de quark pesado más ligeras que las físicas ($M_I < M_{\bar{B}_s}$). Mantienen una razón fija entre las masas del mesón inicial y final ($r_{IF} = M_F/M_I = M_{D_s}/M_{\bar{B}_s}$) para mantener constante el volumen del espacio de fase. El resultado físico se obtiene interpolando estos puntos de datos no perturbativos en el retículo hacia el punto físico utilizando las predicciones de la OPE, las cuales se vuelven exactas en el límite estático ($m_b \to \infty$). Esto transforma el problema de extrapolación en un problema de interpolación.

2. Nueva Estrategia Numérica para Correladores de Cuatro Puntos: Una innovación algorítmica crítica es la introducción de una nueva estrategia para computar los correladores euclidianos de cuatro puntos requeridos. Los autores observan que la "estrategia antigua" (fijar las separaciones temporales entre los operadores de fuente y sumidero) conduce a una pérdida catastrófica de la relación señal-ruido para mesones pesados. La nueva estrategia fija las separaciones temporales entre las inserciones de la corriente débil y los operadores de fuente/sumidero ($t_{snk} - t = 0 - t_{src}$) mientras varía el tiempo de inserción de la corriente $t$. Este enfoque, aunque es computacionalmente más costoso (requiere más inversiones), supera con éxito el problema de la señal-ruido, permitiendo la resolución de correladores en masas de mesones pesados grandes donde el método antiguo falló.

3. Reconstrucción Espectral y Análisis: La tasa de decaimiento se extrae de los correladores euclidianos amputados utilizando el algoritmo de Hansen–Lupo–Tantalo (HLT). Los autores realizan un análisis riguroso de las incertidumbres sistemáticas variando el orden de los límites (continuo $a \to 0$, suavizado $\sigma \to 0$, e integración numérica) y empleando dos enfoques de interpolación distintos hacia el punto físico:

   • Enfoque de Función de Ajuste: Utilizando una función racional para la contribución no perturbativa combinada con términos de la OPE perturbativa.
   • Enfoque de Espacio Funcional: Utilizando Procesos Gaussianos con una distribución previa anclada por el resultado de la OPE perturbativa.

Resultados Clave
El cálculo se realizó en tres conjuntos de ETMC (B64, C80, D96) con espaciados de red $a \approx 0.08, 0.07, 0.06$ fm y volúmenes espaciales $L \approx 5.5$ fm. Los autores obtuvieron una determinación preliminar de la tasa de decaimiento inclusivo:
$$\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}] = 2.47(13)(10)[17] \times 10^{-11} \text{ GeV}$$
donde los errores son estadísticos, sistemáticos (provenientes de los pasos de análisis) y totales, respectivamente. El error relativo total es aproximadamente del 7%.

Utilizando este resultado y despreciando la contribución de $X_{\bar{s}u}$ (justificada por la supresión $|V_{ub}|^2/|V_{cb}|^2$), los autores derivan un valor para $|V_{cb}|$:
$$|V_{cb}|_{\text{latt-incl}} = 40.9(1.2)_{\text{latt}}(1.7)_{\text{exp}}[2.1]_{\text{tot}} \times 10^{-3}$$
Este valor es consistente tanto con la determinación exclusiva ($39.62(47) \times 10^{-3}$) como con la determinación inclusiva basada en OPE ($41.64(47) \times 10^{-3}$), aunque la precisión estadística actual es insuficiente para resolver la tensión existente entre estos métodos.

Significancia y Reivindicaciones
El artículo reclama varias contribuciones significativas al campo:

• Avance Algorítmico: La nueva estrategia numérica para correladores de cuatro puntos resuelve el prolongado problema de la relación señal-ruido para mesones pesados, permitiendo el estudio de decaimientos inclusivos de $B$ en el retículo.
• Validación Metodológica: Este trabajo demuestra que las tasas de decaimiento semileptónico inclusivo pueden computarse desde primeros principios mediante QCD en el retículo a un nivel de precisión fenomenológicamente relevante, combinando datos no perturbativos con restricciones de la OPE.
• Verificación Independiente: Al emplear una discretización de retículo diferente (Masa Torcida) y una técnica de análisis distinta (algoritmo HLT) en comparación con esfuerzos concurrentes (por ejemplo, Ref. [15]), este trabajo proporciona un resultado de primeros principios totalmente independiente para $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$.
• Perspectiva Futura: Los autores enfatizan que el error actual del 7% está dominado por el número limitado de conjuntos de medida y la falta de un estudio de volumen finito basado en datos. Afirman que con la adición de espaciados de red más finos ($a \approx 0.05$ fm) y volúmenes más grandes, el error total puede reducirse significativamente, permitiendo potencialmente una resolución definitiva del enigma de $|V_{cb}|$.

Los autores mantienen un tono modesto respecto al resultado actual, señalando que es un "primer resultado físico", y que las incertidumbres sistemáticas asociadas con los efectos de volumen finito y el límite continuo requieren un mayor refinamiento en publicaciones próximas.