Memoirs of the curvaton: non-perturbative non-Gaussianity and supermassive primordial black holes

Resumo Técnico: "Memoirs of the curvaton: non-perturbative non-Gaussianity and supermassive primordial black holes"

Enunciado do Problema
O paradigma cosmológico padrão postula que as perturbações de curvatura primordial são quase gaussianas e invariantes de escala, consistentes com as observações da Radiação Cósmica de Fundo (CMB). No entanto, em escalas menores ($k \gg \text{Mpc}^{-1}$), essas perturbações permanecem amplamente não restringidas. Um desafio teórico significativo surge ao tentar explicar a formação de Buracos Negros Primordiais (PBHs), particularmente como sementes para Buracos Negros Supermassivos (SMBHs) e os "Little Red Dots" observados pelo JWST. Gerar uma abundância suficiente de PBHs tipicamente requer grandes perturbações de curvatura ($P_\zeta \sim 10^{-2}$), o que está em tensão com as restrições de distorção espectral do CMB ($\mu$-distortions) que limitam a amplitude do espectro de potência a $P_\zeta \lesssim \mathcal{O}(10^{-5})$ em escalas relevantes.

A não-gaussianidade (NG) forte tem sido proposta como um mecanismo para contornar essas restrições, aumentando a cauda de alta amplitude da função de densidade de probabilidade (PDF) sem aumentar o pico do espectro de potência, mas estudos anteriores frequentemente basearam-se em expansões perturbativas (truncando a série $\zeta = \zeta_G + \frac{3}{5}f_{NL}\zeta_G^2 + \dots$) ou ansatzes ad hoc. Essas abordagens falham no regime fortemente não-gaussiano, onde a série perturbativa diverge ou colapsa. Modelos de curvaton existentes muitas vezes assumem um potencial puramente quadrático, negligenciando autointerações que podem alterar significativamente a dinâmica da transição do campo de um estado congelado para oscilações.

Metodologia
Os autores empregam o formalismo $\delta N$ para derivar o mapa de não-gaussianidade local exato e não-perturbativo $\zeta(x) = F(\zeta_G(x))$ para uma classe de potenciais de curvaton com autointeração. A metodologia envolve várias etapas principais:

1. Dinâmica Generalizada do Curvaton: O estudo vai além do padrão de potencial quadrático para incluir potenciais monomiais ($V \propto |\phi|^n$), quárticos ($V \propto \phi^4$) e do tipo axion (cosseno). Os autores levam em conta o fato de que o início das oscilações é dependente do campo nesses modelos, introduzindo uma nova fonte de não-linearidade.
2. Ação Abreviada e Descongelamento Instantâneo: Para lidar com a transição não-adiabática da fase de espectador congelada para a fase oscilatória, os autores utilizam a ação abreviada $W(\rho_\phi)$. Isso permite definir uma equação de estado efetiva e empregar uma aproximação de "descongelamento instantâneo", combinando diretamente o regime congelado com a descrição de fluido adiabático. Isso evita a necessidade de integração numérica computacionalmente cara das equações completas de movimento durante a fase não-adiabática.
3. Mapeamento Não-Perturbativo: Em vez de truncar a expansão de $\zeta$ em termos do campo gaussiano $\zeta_G$, os autores derivam a função de forma fechada $F(\zeta_G)$ diretamente da dinâmica do curvaton. Eles analisam as propriedades estatísticas desse mapeamento, incluindo a PDF de $\zeta$ e o espectro de potência não-gaussiano resultante $P_\zeta(k)$.
4. Espectro de Potência e Estatística de PBH: Os autores computam o espectro de potência físico de forma não-perturbativa usando a integral bivariada do mapa não-linear. Eles aplicam então esses resultados à formação de PBH usando estatísticas de limiar da função de compacidade $C$, calculando a fração de massa $\beta$ e a abundância de PBH resultante $f_{PBH}$.
5. Análise de Restrições: Os cenários derivados são testados contra as restrições de $\mu$-distorção do COBE/FIRAS para determinar se uma abundância viável de PBH pode ser alcançada sem violar os limites de distorção espectral.

Principais Contribuições e Resultados

• Derivação do Mapa Não-Gaussiano: O artigo fornece derivações explícitas do mapeamento $F(\zeta_G)$ para potenciais quadráticos, quárticos, de cosseno e lineares. Demonstra que as autointerações podem aumentar ou suprimir a não-gaussianidade dependendo da forma do potencial e das condições iniciais.
• Ansatz Analítico: Para potenciais quadráticos e do tipo axion, os autores propõem um ansatz logarítmico regularizado:
  $$\zeta \approx \frac{1+\epsilon^2}{2\beta} \ln\left[(1+\beta\zeta_G)^2 + \epsilon^2\right] + F_0$$
  Este ansatz reproduz com precisão os resultados numéricos exatos, capturando o comportamento logarítmico característico e a presença de um mínimo (ou cúspide) no mapeamento, o que é crucial para modelar as caudas da distribuição.
• Espectros de Potência Não-Perturbativos: O estudo mostra que as expansões perturbativas truncadas falham em reproduzir as caudas ultravioletas (UV) corretas do espectro de potência, mesmo em ordens elevadas ($N \sim 50$). O cálculo não-perturbativo completo revela que NGs fortes podem deslocar e aumentar o pico do espectro de potência em relação ao caso gaussiano, uma característica perdida pelos tratamentos perturbativos.
• Impacto das Autointerações: No caso quártico, mostra-se que fortes autointerações suprimem as perturbações de curvatura em grandes valores de $\zeta_G$. Por outro lado, no cenário do tipo axion, escolhas específicas do campo de fundo próximo ao topo do potencial podem gerar grandes NGs positivos mesmo para grandes frações de decaimento ($r_{dec}$).
• Viabilidade do Curvaton do Tipo Axion: O curvaton do tipo axion é identificado como um cenário particularmente natural para este mecanismo. Seu potencial permite os parâmetros necessários de grande $\beta$ negativo (no ansatz) e pequeno $\epsilon$ para gerar a não-gaussianidade necessária sem o grande regulador $\epsilon$ exigido nos modelos quadráticos padrão, que de outra forma aumentaria demais o pico do espectro de potência.
• Formação de PBH e Restrições de $\mu$-distorção: Os autores demonstram um cenário "bottom-up" onde NGs fortemente positivas permitem a formação de sementes de PBH supermassivas ($M_{PBH} \sim 10^5 M_\odot$) com amplitudes de espectro de potência de pico tão baixas quanto $A_{pk} \sim 10^{-3}$. Esta amplitude é compatível com as restrições de $\mu$-distorção do COBE/FIRAS, ao contrário do caso Gaussiano, que exigiria $A_{pk} \sim 10^{-2}$ (violando os limites) para produzir a mesma abundância.

Significância
O artigo afirma que a compreensão profunda do mapa não-gaussiano não-perturbativo é central para estimativas confiáveis da abundância de PBH em cenários de curvaton. Os autores argumentam que as técnicas perturbativas padrão são insuficientes para regimes fortemente não-gaussianos, levando a previsões incorretas tanto para a forma do espectro de potência quanto para a abundância de PBH resultante.

Ao incorporar a dinâmica completa não-linear de curvatons com autointeração, o estudo abre um novo cenário primordial para a origem de SMBHs e os "Little Red Dots" observados pelo JWST. Especificamente, sugere que PBHs pesados podem se formar diretamente via colapso gravitacional de perturbações de curvatura que são consistentes com as restrições de distorção espectral do CMB, evitando a necessidade de processos extremos de acreção ou histórias de fusão hierárquica. Os autores observam que, embora os PBHs tenham sido a aplicação primária, o framework é geral e aplicável a outras observáveis sensíveis à cauda de alto $\zeta$, como ondas gravitacionais induzidas por escalar (SIGWs).