Bound-state decay of a charged scalar particle in an external Coulomb field

Resumo Técnico: Decaimento de Estado Ligado de uma Partícula Escalar Carregada em um Campo de Coulomb Externo

Enunciado do Problema
O artigo aborda o cálculo da taxa de decaimento $\Gamma$ para uma partícula ligada no campo de Coulomb de um núcleo, examinando especificamente a interação entre a dinâmica de decaimento, efeitos de ligação e interações de estado final. Embora estudos anteriores tenham analisado o decaimento de múons ligados, formulações numéricas que dependem de expansões de ondas parciais de funções de onda contínuas enfrentam desafios significativos para núcleos leves ($Z < 6$). Neste regime, o parâmetro de ligação $\alpha Z$ é pequeno, mas a série de momento angular converge lentamente. Cálculos diretos de ondas parciais são computacionalmente dispendiosos porque cada termo envolve integrais radiais altamente oscilatórias, e reconstruir uma razão de taxa de decaimento próxima de umidade com precisão suficiente para resolver pequenas correções de Coulomb é numericamente exigente.

Metodologia
Para investigar esses problemas de convergência, os autores empregam um modelo escalar simplificado: uma partícula escalar negativamente carregada ($K^-$) ligada a um núcleo estático, infinitamente pesado e pontual, decai em partículas escalares carregadas ($\pi^-$) e neutras ($\pi^0$) via uma interação escalar tripla local. Todas as partículas são tratadas como pontuais, e os píons finais são de massa nula. O sistema é descrito pela equação de Klein-Gordon em um potencial de Coulomb externo.

A inovação metodológica central é um método de subtração projetado para superar a lenta convergência da expansão de ondas parciais:

1. Separação de Contribuições: A taxa de decaimento total é decomposta em uma contribuição de "píon livre" (negligenciando as interações de Coulomb no estado final) e uma correção de "distorção de Coulomb" ($\Delta C$).
2. Redução Analítica: A contribuição do píon livre é avaliada separadamente usando uma representação integral unidimensional derivada da transformada de Fourier da função de onda do estado ligado. Isso evita a necessidade de uma soma de ondas parciais para o termo dominante.
3. Aplicação de Ondas Parciais: A expansão de ondas parciais é aplicada apenas à diferença entre a função de onda do píon interagente e a função de onda do píon livre (a distorção de Coulomb). Como esta correção é significativamente menor que a taxa total, a soma de ondas parciais converge muito mais rápido e requer um corte (cutoff) inferior ($l_{max}$).
4. Verificação Analítica: Os resultados numéricos são verificados cruzadamente usando uma expansão diagramática de pequeno $\alpha Z$ através da ordem $(\alpha Z)^4$. Isso envolve o cálculo da parte imaginária da autoenergia do kaon no campo externo, considerando inserções de fótons de Coulomb de zero, uma e duas linhas de píons, incluindo diagramas de seagull e de escada (ladder).

Principais Resultados

• Eficiência Numérica: Para núcleos leves ($1 \le Z \le 10$), o método de subtração reduz o corte de onda parcial necessário em aproximadamente 25–30% em comparação com cálculos diretos. Para $Z=1$, o corte é reduzido de $l_{max}=1136$ (direto) para $802$ (subtraído). Os resultados para $1 \le Z \le 5$ de ambos os métodos concordam até quatro dígitos significativos.
• Razão da Taxa de Decaimento: A razão entre as taxas de decaimento ligada e livre é determinada numericamente. O artigo apresenta a supressão da taxa $1 - \Gamma/\Gamma_0$ para $Z$ até 10.
• Expansão Analítica: A expansão diagramática fornece a seguinte série para a razão da taxa de decaimento:
  $$\frac{\Gamma}{\Gamma_0} = 1 - \frac{1}{2}(\alpha Z)^2 + \frac{17}{24}(\alpha Z)^4 + O((\alpha Z)^5)$$
  O termo líder $-(\alpha Z)^2/2$ provém inteiramente da contribuição do píon livre. A distorção de Coulomb da partícula de saída não contribui na ordem $(\alpha Z)^2$, mas fornece a correção líder na ordem $(\alpha Z)^4$.
• Coeficiente de Distorção de Coulomb: O coeficiente do termo $(\alpha Z)^4$ proveniente da distorção de Coulomb é encontrado como $4/3$. Isso é confirmado numericamente, conforme a quantidade $\Delta C / (\alpha Z)^4$ aproxima-se de $4/3$ conforme $Z$ diminui.
• Mecanismo de Cancelamento: Na ordem $(\alpha Z)^4$, as contribuições de dois fótons (gráfico seagull e duas inserções lineares) se cancelam exatamente ($+2(\alpha Z)^4 - 2(\alpha Z)^4 = 0$), deixando a contribuição de um fóton como a única fonte da correção de distorção nesta ordem.

Significância e Alegações
O artigo alega que o método de subtração resolve efetivamente a instabilidade numérica associada aos cálculos de decaimento de estado ligado de baixo $Z$, um problema que anteriormente levava pesquisadores a omitir núcleos leves ($Z < 6$) das análises de decaimento de múon ligado. Ao isolar a pequena correção de distorção de Coulomb, o método evita a dificuldade numérica de reconstruir uma razão próxima de unidade a partir de termos grandes e quase canceláveis.

Os autores enfatizam que, embora o coeficiente líder da supressão ($-1/2$) coincida com o resultado para o decaimento de múon ligado (o resultado de Überall), a origem física difere: neste modelo escalar, a distorção de Coulomb da partícula de saída só entra na ordem $(\alpha Z)^4$, enquanto no caso do múon de spin-1/2, ela contribui na ordem $(\alpha Z)^2$. O trabalho fornece um arcabouço validado e ferramentas (a estratégia de subtração e a expansão diagramática) para futuros cálculos mais realistas de decaimentos de estado ligado, incluindo o problema completo de spin-1/2 do decaimento de múon ligado.