Precision Study of Semileptonic and Non-Leptonic $B_c$ Decays to $\eta_c$ and P Wave Charmonia

Resumo Técnico: Estudo de Precisão de Decaimentos Semilepônicos e Não Leptônicos de $B_c$ para $\eta_c$ e Charmônios de Onda P

Definição do Problema
O méson $B_c$, composto por dois quarks pesados ($b\bar{c}$), oferece um laboratório único para testar a dinâmica de quarks pesados e buscar física além do Modelo Padrão (SM), particularmente através da violação da Universalidade de Sabor Leptônico (LFU). Enquanto decaimentos semilepônicos para charmônios de onda S (como o $J/\psi$) são bem estudados, os decaimentos para charmônios de onda P ($\chi_{c0}, \chi_{c1}, h_c$) permanecem pouco compreendidos. Esses modos de onda P são críticos não apenas como potenciais fundos para análises de onda S, mas também como sondas independentes de efeitos não perturbativos de QCD e simetria de spin. Um grande obstáculo em prever esses decaimentos é a falta de determinações confiáveis de QCD em rede (lattice QCD) ou experimentais diretas dos fatores de forma de transição relevantes em todo o intervalo cinemático de $q^2$. As previsões existentes dependem fortemente de abordagens dependentes de modelos (por exemplo, QCD perturbativa, modelos de quarks, somas de QCD), levando a incertezas teóricas significativas.

Metodologia
Os autores empregam uma abordagem orientada por dados dentro do arcabouço de fatoração de QCD não-relativística (NRQCD) para restringir os fatores de forma de transição para decaimentos de $B_c \to$ charmônios de onda P. A metodologia procede em várias etapas fundamentais:

1. Extração de Parâmetros da Função de Onda: Em vez de depender de modelos teóricos para os inputs não-perturbativos, os autores extraem as derivadas das funções de onda radiais na origem ($\psi'_{M}(0)$) para $\chi_{c0}, \chi_{c1}$ e $h_c$, bem como as funções de onda para $J/\psi$ e $\eta_c$. Isso é alcançado realizando um ajuste global de $\chi^2$ aos dados experimentais de decaimentos radiativos (por exemplo, $\chi_{c1} \to V\gamma, \chi_{c0} \to \gamma\gamma, h_c \to \eta_c\gamma$) e decaimentos leptônicos, combinados com inputs de QCD em rede para constantes de decaimento. Parâmetros de incerteza (nuisance parameters), incluindo correções relativísticas e operadores de ordem superior, são tratados com priors gaussianos.
2. Cálculo do Fator de Forma no Recuo Máximo: Usando os parâmetros de função de onda extraídos e as expressões analíticas de NRQCD (incluindo Ordem Líder e Correções Relativísticas até $O(v^2)$) da Ref. [45], os autores calculam os fatores de forma de transição em $q^2=0$ (recuo máximo). Para $B_c \to \chi_{c1}$ e $h_c$, onde correções de Próxima Ordem (NLO) não estão disponíveis, uma incerteza conservadora de 20% é atribuída para contabilizar efeitos perturbativos de ordem superior ausentes.
3. Determinação da Forma e Extrapolação: Para descrever os fatores de forma em toda a região física de $q^2$, os autores utilizam dois métodos complementares:
   • Expansão de Polo: Os parâmetros são restringidos usando argumentos de simetria de spin, assumindo que os fatores de forma para transições para estados com o mesmo spin total $J$ compartilham os mesmos parâmetros de forma ($\alpha_i, \beta$). Esses parâmetros são primeiro determinados ajustando dados sintéticos gerados para $B_c \to J/\psi$ e $B_c \to \eta_c$ (baseados em estudos prévios de lattice e HQSS).
   • Expansão $z$ (Parametrização BCL): Os resultados da expansão de polo são usados para gerar pontos de pseudo-dados, que são então ajustados usando a parametrização Boyd-Caprini-Lellouch (BCL), que é independente de modelo. Esta abordagem garante que as restrições de analiticidade e unitariedade sejam satisfeitas, fornecendo uma descrição robusta da dependência de $q^2$.
4. Previsões Fenomenológicas: Os fatores de forma restringidos são usados para computar frações de branching semilepônicas, razões de LFU ($R(H)$), taxas de decaimento não-leptônicos exclusivos (usando fatorização de QCD com correções não-fatorizáveis) e seções de choque de produção de charmônios de onda P em decaimentos de bósons $Z$ e aniquilação $e^+e^-$.

Principais Contribuições e Resultados

• Fatores de Forma Orientados por Dados: O artigo fornece as primeiras estimativas de fatores de forma para $B_c \to \chi_{c0,1}, h_c$ restringidas por dados. Ao extrair as derivadas das funções de onda diretamente das taxas de decaimento radiativo experimental, os autores minimizam a dependência de modelos puramente teóricos.
• Razões de LFU: Os autores apresentam previsões precisas do Modelo Padrão para as razões de LFU:
  • $R(\chi_{c0}) = 0,185(3)$
  • $R(\chi_{c1}) = 0,147(26)$
  • $R(h_c) = 0,068(2)$
    Essas razões são teoricamente limpas, pois as incertezas hadrônicas largamente se cancelam entre o numerador e o denominador. A precisão é limitada primariamente pelas incertezas dos fatores de forma, que são estimadas em torno de 5% para $R(\chi_{c0})$ e $R(h_c)$.
• Frações de Branching Semilepônicas: Previsões atualizadas para frações de branching semilepônicas são fornecidas, por exemplo, $\mathcal{B}(B_c \to \eta_c \tau^- \bar{\nu}_\tau) = (1,90 \pm 0,49) \times 10^{-3}$ e $\mathcal{B}(B_c \to \eta_c \mu^- \bar{\nu}_\mu) = (6,51 \pm 1,97) \times 10^{-3}$. As incertezas permanecem relativamente grandes (25–30%) devido à atual falta de inputs de lattice precisos para os fatores de forma de onda P.
• Decaimentos Não-Leptônicos: O estudo estende-se a decaimentos não-leptônicos exclusivos de dois corpos ($B_c \to \chi_{cJ}, h_c + \pi, K, \rho, K^*$). Previsões são fornecidas para frações de branching, incluindo correções fatorizáveis de Próxima Ordem (NLO) e correções não-fatorizáveis de Ordem Líder (LO). Por exemplo, $\mathcal{B}(B_c \to \chi_{c0}\pi)$ é previsto como $5,3(60) \times 10^{-3}$ (com correções).
• Modos de Produção: O artigo atualiza as previsões para a produção de charmônio de onda P em decaimentos de bósons $Z$ (por exemplo, $Z \to \chi_{c0}\gamma$) e produção de charmônio duplo em aniquilação $e^+e^-$ (por exemplo, $e^+e^- \to J/\psi + \chi_{c0}$), comparando-as com estimativas da literatura anterior.

Significância e Alegações
Os autores alegam que este trabalho representa um dos primeiros esforços para restringir os fatores de forma de charmônio de onda P para $B_c$ usando dados experimentais em vez de depender exclusivamente de modelos teóricos. Ao combinar o formalismo de NRQCD com um ajuste global aos dados de decaimento radiativo, eles fornecem uma base "controlada e sistemática" para avaliar as larguras de decaimento.

O artigo enfatiza que, embora as previsões atuais tenham incertezas consideráveis (impulsionadas pelos fatores de forma hadrônicos mal restringidos), as razões de LFU derivadas são sondas teoricamente limpas. Os autores afirmam que suas previsões estão prontas para serem testadas por experimentos futuros no LHC e Belle II. Eles explicitamente observam que a análise atual é um passo intermediário, antecipando que futuros cálculos de QCD em rede e de Somas de Cone de Luz (LCSR) permitirão um tratamento independente de modelo para esses fatores de forma, possibilitando estimativas do SM mais precisas e buscas mais sensíveis por Nova Física. O trabalho não pretende resolver as tensões existentes em $R(J/\psi)$, mas destaca a necessidade de compreender os fundos de onda P e fornecer testes de precisão complementares no setor pesado-pesado.