Thermal axion production in QCD from the lattice

Resumo Técnico: Produção Térmica de Áxions em QCD a partir da Rede (Lattice)

Definição do Problema
A taxa de produção de áxions térmicos no Universo primordial é governada pela função espectral da densidade de carga topológica da QCD, $\sigma(k, k)$, avaliada na luz-cone. Enquanto métodos perturbativos funcionam em temperaturas muito altas e a teoria de perturbação quiral se aplica abaixo da região de transição, a região de crossover ($T \sim 150\text{--}300$ MeV) é fortemente acoplada, tornando essas técnicas não confiáveis. As abordagens padrão de QCD em rede enfrentam um obstáculo fundamental: elas computam funções de correlação euclidianas, enquanto a taxa de produção requer a função espectral de Minkowski. Tentativas anteriores de preencher essa lacuna basearam-se em continuação analítica, um processo propenso a grandes erros sistemáticos devido à natureza mal posta da inversão. Este artigo aborda o desafio de computar integrais ponderadas da taxa de produção de áxions diretamente a partir de dados de rede euclidianos sem a necessidade de continuação analítica.

Metodologia
Os autores propõem um método inspirado em trabalhos recentes sobre taxas de produção de fótons (Harvey Meyer et al.), que permite o cálculo de integrais ponderadas de funções espectrais usando corretores euclidianos.

1. Estrutura Teórica:

   • A taxa de produção de áxions está relacionada à função espectral $\sigma(\omega, \omega)$.
   • Os autores definem um correlador euclidiano $H_E(\omega_n)$ envolvendo momento espacial imaginário, que se relaciona com a função espectral via uma integral de kernel.
   • A avaliação direta de $H_E(\omega_n)$ é problemática porque o kernel cresce exponencialmente com a distância, e a integral é divergente no ultravioleta (UV) e potencialmente em altas frequências.
   • Para resolver os problemas de divergência, os autores constroem combinações lineares de modos de Matsubara, especificamente $H_E^{(0,1)}$ e $H_E^{(0,1,2)}$. Essas combinações introduzem kernels de subtração que cancelam as principais divergências UV e garantem a convergência em altas frequências.

2. Implementação em Rede (Lattice):

   • Configuração: O estudo é realizado em teoria de Yang-Mills $SU(3)$ pura (aproximação quenched) a uma temperatura $T = 1.5 T_c$ em uma rede com $N_\tau = 16$ (aspect ratio 4).
   • Fluxo de Gradiente (Gradient Flow): Para lidar com a renormalização multiplicativa e o ruído de curta distância do operador de densidade de carga topológica $q(x)$, os autores empregam o fluxo de gradiente melhorado (Zeuthen). Isso suaviza o operador sobre um raio $r \sim \sqrt{8t_f}$, estabilizando a carga topológica.
   • Ajuste de Cauda (Tail Fitting): O correlator euclidiano $G(\omega_n, r)$ decai exponencialmente a grandes distâncias. No entanto, a função de peso $\sinh(\omega_n r)/(\omega_n r)$ cresce exponencialmente. Para computar a integral, os autores ajustam a "cauda" do correlator (grande $r$) a uma soma de exponenciais ($e^{-mr}/r$) derivada da representação espectral da teoria em uma geometria compactificada. A integral é então computada combinando dados precisos de rede em pequenos $r$ com a função ajustada em grandes $r$.
   • Extrapolação: Os resultados são obtidos em tempos de fluxo finitos $t_f$. Os autores argumentam que para $t_f \gtrsim 0.4a^2$, a função de suavização é efetivamente invariante sob $O(4)$, o que implica que a dependência do tempo de fluxo da quantidade integrada $H_E$ deve ser desprezível (escalando como $t_f^2$ ou superior). Eles realizam um ajuste constante nos dados através de diferentes tempos de fluxo para extrapolar para $t_f \to 0$.

Principais Contribuições

• Método Livre de Continuação Analítica: O artigo demonstra um caminho viável para extrair integrais ponderadas da taxa de produção de áxions diretamente de simulações de QCD em rede euclidiana, evitando a necessidade de continuação analítica instável.
• Tratamento de Ruído Exponencial: Os autores desenvolvem uma estratégia numérica robusta combinando fluxo de gradiente com técnicas de ajuste de cauda para superar a degradação da relação sinal-ruído causada pela função de peso exponencialmente crescente no integral.
• Esquemas de Subtração: A formulação de combinações específicas de modos de Matsubara ($H_E^{(0,1)}$ e $H_E^{(0,1,2)}$) que tornam os integrais finitos e bem definidos na rede.

Resultos

• Valores Numéricos: Em QCD de glúons puros a $T = 1.5 T_c$, os autores extraem com sucesso:
  • $H_E^{(0,1)} = -0.0357(58) T^4$
  • $H_E^{(0,1,2)} = -0.0152(91) T^4$
    O sinal negativo é consistente com a natureza anti-hermitiana do operador de densidade de carga topológica.
• Comparação com Teoria de Perturbação: Os resultados são comparados com cálculos de perturbação de Próxima Ordem (NLO). Os resultados da rede são consistentes com as previsões de NLO (usando $\alpha_s \approx 0.3$), enquanto os resultados de Ordem Líder (LO) são significativamente menores.
• Viabilidade: O estudo confirma que a extração de pelo menos um, e possivelmente dois, integrais ponderados é viável com um esforço numérico modesto (aprox. 60.000 horas-núcleo de CPU) em um único espaçamento de rede.

Significância e Alegações
O artigo afirma abrir um novo caminho para determinar as taxas de produção térmica de áxions no regime fortemente acoplado da QCD, onde outros métodos falham. Ao computar com sucesso esses integrais em teoria de glúons puros, os autores estabelecem a viabilidade da abordagem. Eles postulam que este método pode ser estendido para a QCD completa em temperaturas ao redor de e ligeiramente acima do crossover ($T \in [200, 1000]$ MeV), um intervalo crítico para restrições cosmológicas sobre a matéria escura de áxions.

Os autores permanecem modestos quanto ao impacto fenomenológico imediato, observando que seu método fornece momentos integrais em vez da taxa de produção diferencial $d\Gamma/d^3k$ como função do momento. Eles reconhecem que traduzir essas restrições integrais em limites precisos sobre a massa e o acoplamento do áxon requer mais desenvolvimento. Além disso, destacam que o sucesso do método depende de o estado $0^{-+}$ sem sabor mais leve ser suficientemente pesado em comparação com a escala de temperatura para garantir que o correlator decaia mais rápido do que a função de peso crescente; esta condição pode tornar-se desafiadora em temperaturas muito altas, onde o acoplamento se torna fraco.