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

Resumo Técnico: Decaimentos Inclusivos $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ a partir de QCD em Rede

Problema
A determinação do elemento da matriz Cabibbo–Kobayashi–Maskawa (CKM) $|V_{cb}|$ é um desafio central na física de sabor, atualmente caracterizada por uma tensão significativa entre as determinações exclusivas e inclusivas experimentais. Embora os decaimentos semilepônicos inclusivos de mésons $B$ ofereçam um caminho teoricamente robusto via Expansão de Produtos Operadores (OPE), o cálculo de primeiros princípios via QCD em rede da taxa de decaimento inclusivo $\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}]$ tem sido historicamente considerado inacessível devido ao "problema das duas escalas" (satisfazer simultaneamente $a m_b \ll 1$ e $\Lambda_{QCD} L \gg 1$) e à rápida deterioração da razão sinal-ruído dos correladores de quatro pontos conforme a massa do quark pesado aumenta.

Metodologia e Estratégia Computacional
Os autores apresentam uma estratégia inovadora para computar esta taxa de decaimento inclusivo a partir de primeiros princípios, utilizando conjuntos de gauge $N_f = 2+1+1$ de Massa Twist Estendida (ETMC) com massas de quarks leves, estranhos e de charm fisicamente reais. A metodologia baseia-se em três pilares principais:

1. Interpolação de Massa via OPE: Em vez de simular diretamente na massa física do quark $b$ (o que é computacionalmente proibitivo devido aos efeitos de cutoff), os autores simulam em massas de quarks pesados mais leves do que as físicas ($M_I < M_{\bar{B}_s}$). Eles mantêm uma razão fixa entre as massas do mensageiro inicial e final ($r_{IF} = M_F/M_I = M_{D_s}/M_{\bar{B}_s}$) para manter o volume do espaço de fase constante. O resultado físico é obtido interpolando estes pontos de dados não perturbativos da rede para o ponto físico usando previsões de OPE, que se tornam exatas no limite estático ($m_b \to \infty$). Isso transforma o problema de extrapolação em um problema de interpolação.

2. Nova Estratégia Numérica para Correladores de Quatro Pontos: Uma inovação algorítmica crítica é a introdução de uma nova estratégia para computar os correladores euclidianos de quatro pontos necessários. Os autores observam que a "estratégia antiga" (fixar as separações temporais entre os operadores de fonte e sumidouro) leva a uma perda catastrófica da razão sinal-ruído para mésons pesados. A nova estratégia fixa as separações temporais entre as inserções de corrente fraca e os operadores de fonte/sumidouro ($t_{snk} - t = 0 - t_{src}$) enquanto varia o tempo de inserção da corrente $t$. Esta abordagem, embora computacionalmente mais cara (exigindo mais inversões), supera com sucesso o problema do sinal-ruído, permitindo a resolução de correladores em grandes massas de mésons pesados onde o método antigo falhou.

3. Reconstrução Espectral e Análise: A taxa de decaimento é extraída dos correladores euclidianos amputados usando o algoritmo Hansen–Lupo–Tantalo (HLT). Os autores realizam uma análise rigorosa das incertezas sistemáticas variando a ordem dos limites (contínuo $a \to 0$, smearing $\sigma \to 0$ e integração numérica) e empregando duas abordagens de interpolação distintas para o ponto físico:

   • Abordagem de Função de Ajuste: Utilizando uma função racional para a contribuição não perturbativa combinada com termos de OPE perturbativos.
   • Abordagem de Espaço Funcional: Utilizando Processos Gaussianos com uma distribuição a priori ancorada pelo resultado de OPE perturbativo.

Resultados Principais
O cálculo foi realizado em três conjuntos ETMC (B64, C80, D96) com espaçamentos de rede $a \approx 0.08, 0.07, 0.06$ fm e volumes espaciais $L \approx 5.5$ fm. Os autores obtiveram uma determinação preliminar da taxa de decaimento inclusivo:
$$\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}] = 2.47(13)(10)[17] \times 10^{-11} \text{ GeV}$$
onde os erros são estatísticos, sistemáticos (da análise de passos) e totais, respectivamente. O erro relativo total é de aproximadamente 7%.

Usando este resultado e negligenciando a contribuição $X_{\bar{s}u}$ (justificada pela supressão $|V_{ub}|^2/|V_{cb}|^2$), os autores derivam um 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 é consistente tanto com a determinação exclusiva ($39.62(47) \times 10^{-3}$) quanto com a determinação inclusiva baseada em OPE ($41.64(47) \times 10^{-3}$), embora a precisão estatística atual seja insuficiente para resolver a tensão existente entre estes métodos.

Significância e Alegações
O artigo reivindica várias contribuições significativas para o campo:

• Avanço Algorítmico: A nova estratégia numérica para correladores de quatro pontos resolve o problema de sinal-ruído de longa data para mésons pesados, permitindo o estudo de decaimentos inclusivos de $B$ em rede.
• Validação Metodológica: O trabalho demonstra que taxas de decaimento semilepônicos inclusivos podem ser computadas a partir de primeiros princípios em QCD em rede a um nível de precisão fenomenológico, combinando dados não perturbativos com restrições de OPE.
• Verificação Independente: Ao empregar uma discretização de rede diferente (Twisted Mass) e uma técnica de análise distinta (algoritmo HLT) em relação a esforços concomitantes (ex: Ref. [15]), este trabalho fornece um resultado de primeiros princípios totalmente independente para $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$.
• Perspectiva Futura: Os autores enfatizam que o erro atual de 7% é dominado pelo número limitado de conjuntos de rede e pela falta de um estudo de volume finito baseado em dados. Eles afirmam que, com a adição de espaçamentos de rede mais finos ($a \approx 0.05$ fm) e volumes maiores, o erro total pode ser significativamente reduzido, potencialmente permitindo uma resolução definitiva do enigma do $|V_{cb}|$.

Os autores mantêm um tom modesto em relação ao resultado atual, observando que é um "primeiro resultado físico", e que as incertezas sistemáticas associadas aos efeitos de volume finito e ao limite contínuo requerem mais refinamento em publicações futuras.