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

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

Planteamiento del Problema
El paradigma cosmológico estándar postula que las perturbaciones de curvatura primordiales son casi gaussianas e invariantes de escala, en consistencia con las observaciones del Fondo Cósmico de Microondas (CMB). Sin embargo, en escalas pequeñas ($k \gg \text{Mpc}^{-1}$), estas perturbaciones permanecen mayormente sin restringir. Un desafío teórico significativo surge al intentar explicar la formación de Agujeros Negros Primordiales (PBH), particularmente como semillas para Agujeros Negros Supermasivos (SMBH) y los "Little Red Dots" observados por el JWST. Generar una abundancia suficiente de PBH requiere típicamente perturbaciones de curvatura grandes ($P_\zeta \sim 10^{-2}$), lo cual entra en tensión con las restricciones de distorsión espectral del CMB ($\mu$-distorsiones) que limitan la amplitud del espectro de potencia a $P_\zeta \lesssim \mathcal{O}(10^{-5})$ en las escalas relevantes.

Se ha propuesto que la no-gaussianidad (NG) fuerte es un mecanismo para evadir estas restricciones al mejorar la cola de alta amplitud de la función de densidad de probabilidad (PDF) sin aumentar el pico del espectro de potencia, pero los estudios previos a menudo han dependido de expansiones perturbativas (truncando la serie $\zeta = \zeta_G + \frac{3}{5}f_{NL}\zeta_G^2 + \dots$) o de ansätze ad hoc. Estos enfoques fallan en el régimen fuertemente no-gaussiano donde la serie perturbativa diverge o colapsa. Además, los modelos de curvatón existentes suelen asumir un potencial puramente cuadrático, omitiendo las auto-interacciones que pueden alterar significativamente la dinámica del campo en su transición de un estado congelado a oscilaciones.

Metodología
Los autores emplean el formalismo $\delta N$ para derivar el mapa local no-perturbativo de no-gaussianidad $\zeta(x) = F(\zeta_G(x))$ para una clase de potenciales de curvatón con auto-interacción. La metodología involucra varios pasos clave:

1. Dinámica Generalizada del Curvatón: El estudio va más allá del estándar potencial cuadrático para incluir potenciales monomiales ($V \propto |\phi|^n$), cuárticos ($V \propto \phi^4$) y de tipo axión (coseno). Los autores dan cuenta del hecho de que el inicio de las oscilaciones es dependiente del campo en estos modelos, introduciendo una nueva fuente de no-linealidad.
2. Acción Abreviada y Descongelamiento Instantáneo: Para manejar la transición no-adiabática de la fase de espectador congelado a la fase oscilatoria, los autores utilizan la acción abreviada $W(\rho_\phi)$. Esto les permite definir una ecuación de estado efectiva y emplear una aproximación de "descongelamiento instantáneo", acoplando directamente el régimen congelado con la descripción del fluido adiabático. Esto evita la necesidad de una integración numérica costosa de las ecuaciones completas de movimiento durante la fase no-adiabática.
3. Mapeo No-Perturbativo: En lugar de truncar la expansión de $\zeta$ en términos del campo gaussiano $\zeta_G$, los autores derivan la función de forma cerrada $F(\zeta_G)$ directamente de la dinámica del curvatón. Analizan las propiedades estadísticas de este mapeo, incluyendo la PDF de $\zeta$ y el espectro de potencia no-gaussiano resultante $P_\zeta(k)$.
4. Espectro de Potencia y Estadísticas de PBH: Los autores computan el espectro de potencia físico de forma no-perturbativa utilizando la integral bivariante del mapa no-lineal. Posteriormente, aplican estos resultados a la formación de PBH mediante estadísticas de umbral para la función de compactación $C$, calculando la fracción de masa $\beta$ y la abundancia de PBH resultante $f_{PBH}$.
5. Análisis de Restricciones: Los escenarios derivados se prueban contra los límites de $\mu$-distorsión de COBE/FIRAS para determinar si es posible lograr una abundancia viable de PBH sin violar las restricciones de distorsión espectral.

Contribuciones Clave y Resultados

• Derivación del Mapa de No-Gaussianidad: El artículo proporciona derivaciones explícitas del mapeo $F(\zeta_G)$ para potenciales cuadráticos, cuárticos, de tipo axión y lineales. Demuestra que las auto-interacciones pueden potenciar o suprimir la no-gaussianidad dependiendo de la forma del potencial y las condiciones iniciales.
• Ansatz Analítico: Para potenciales cuadráticos y de tipo axión, los autores proponen un 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 reproduce con precisión los resultados numéricos exactos, capturando el comportamiento logarítmico característico y la presencia de un mínimo (o cúspide) en el mapeo, lo cual es crucial para modelar las colas de la distribución.
• Espectros de Potencia No-Perturbativos: El estudio muestra que las expansiones perturbativas truncadas fallan al reproducir las colas ultravioletas (UV) correctas del espectro de potencia, incluso en órdenes altos ($N \sim 50$). El cálculo completo no-perturbativo revela que las NG fuertes pueden desplazar y potenciar el pico del espectro de potencia respecto al caso gaussiano, una característica omitida por los tratamientos perturbativos.
• Impacto de las Auto-interacciones: En el caso cuártico, se demuestra que las fuertes auto-interacciones suprimen las perturbaciones de curvatura para valores grandes de $\zeta_G$. Por el contrario, en el escenario de tipo axión, elecciones específicas del campo de fondo cerca de la cima del potencial pueden generar grandes NG positivas incluso para fracciones de decaimiento ($r_{dec}$) grandes.
• Formación de PBH y Restricciones de $\mu$-distorsión: Los autores demuestran un escenario "bottom-up" donde las NG fuertemente positivas permiten la formación de semillas de PBH supermasivos ($M_{PBH} \sim 10^5 M_\odot$) con amplitudes de espectro de potencia pico tan bajas como $A_{pk} \sim 10^{-3}$. Esta amplitud es compatible con los límites de $\mu$-distorsión de COBE/FIRAS, a diferencia del caso gaussiano que requeriría un $A_{pk} \sim 10^{-2}$ (violando los límites) para producir la misma abundancia.
• Viabilidad del Curvatón de Tipo Axión: El curvatón de tipo axión se identifica como un entorno particularmente natural para este mecanismo. Su potencial permite los parámetros de $\beta$ grande y negativo (en el ansatz) y $\epsilon$ pequeño necesarios para generar la no-gaussianidad requerida sin el regulador $\epsilon$ grande requerido en los modelos cuadráticos estándar, lo que de otro modo potenciaría demasiado el pico del espectro de potencia.

Significancia
El artículo sostiene que una comprensión profunda del mapa de no-gaussianidad no-perturbativo es central para estimaciones fiables de la abundancia de PBH en escenarios de curvatón. Los autores argumentan que las técnicas perturbativas estándar son insuficientes para los regímenes fuertemente no-gaussianos, llevando a predicciones incorrectas tanto para la forma del espectro de potencia como para la abundancia de PBH resultante.

Al incorporar la dinámica completa no-lineal de los curvatones con auto-interacción, el estudio abre un nuevo escenario primordial para el origen de los SMBH y los "Little Red Dots" observados por el JWST. Específicamente, sugiere que los PBH pesados pueden formarse directamente mediante el colapso gravitacional de perturbaciones de curvatura que son consistentes con las restricciones de distorsión espectral del CMB, evitando la necesidad de historias de acreción extrema o de fusiones jerárquicas. Los autores señalan que, aunque los PBH fueron la aplicación principal, el marco es general y aplicable a otros observables sensibles a la cola de alto-$\zeta$, tales como las ondas gravitacionales inducidas por escalares (SIGWs).