Reheating matters: Starobinsky inflation in light of joint CMB+BAO results and gravitational-wave forecasts

Resumen Técnico: El recalentamiento importa: La inflación de Starobinsky a la luz de los resultados conjuntos de CMB+BAO y los pronósticos de ondas gravitacionales

Planteamiento del problema
Análisis conjuntos recientes de datos del Fondo Cósmico de Microondas (CMB) del Telescopio Cosmológico de Atacama (ACT) y de Oscilaciones Acústicas de Bariones (BAO) del Instrumento de Espectroscopía de Energía Oscura (DESI) han arrojado un índice espectral escalar de $n_s = 0.9752 \pm 0.0030$. Este valor representa una tensión de $\sim 2\sigma$ con las predicciones estándar de la inflación de Starobinsky, que típicamente arrojan un $n_s \approx 0.960\text{--}0.967$ bajo el supuesto de un episodio de recalentamiento estándar dominado por un fluido de tipo materia ($N_* \approx 50\text{--}60$).

La literatura sugiere que un episodio de recalentamiento post-inflacionario dominado por un fluido "rígido" (ecuación de estado $w_{\text{reh}} > 1/3$) podría aumentar el número de e-folds ($N_*$) y desplazar la predicción de $n_s$ hacia arriba para alinearse con estos nuevos conjuntos de datos. Sin embargo, un episodio rígido induce un sesgo azul (blue-tilt) en el espectro de ondas gravitacionales (GW) primordiales. El problema central abordado es si tales escenarios dominados por un fluido rígido pueden satisfacer simultáneamente las restricciones de CMB+BAO sobre $n_text{s}$ sin violar los límites de la Nucleosíntesis del Big Bang (BBN) sobre la densidad de radiación ($\Delta N_{\text{eff}}$), y si el resultante espectro de GW con sesgo azul es detectable por futuros interferómetros.

Metodología
Los autores realizan un análisis estadístico bayesiano utilizando el marco de trabajo Cobaya y una versión modificada del código de Boltzmann CLASS. Utilizan el conjunto de datos exhaustivo P-ACT-LB2-BK18 (Planck 2018 + ACT DR6 + DESI DR2 + BICEP/Keck 2018).

El análisis procede en tres etapas:

1. Exploración del espacio de parámetros: Los autores muestrean el espacio de parámetros del modelo de Starobinsky, específicamente la escala del potencial inflacionario ($\ln V_0$) y el parámetro de recalentamiento reescalado ($\ln R_{\text{reh}}$), para derivar restricciones sobre la ecuación de estado efectiva ($w_{\text{reh}}$) y la temperatura de recalentamiento ($T_{\text{reh}}$).
2. Modelado de escenarios: Se investigan tres historias de recalentamiento distintas:
   • Rígido puro: El recalentamiento está dominado enteramente por un fluido con una $w_{\text{reh}}$ constante $> 1/3$.
   • Materia + Rígido: El recalentamiento comienza con una fase dominada por materia (imitando las oscilaciones del inflatón alrededor de un potencial cuadrático) seguida de una fase dominada por un fluido rígido.
   • Radiación + Rígido: El recalentamiento comienza con una fase dominada por radiación seguida de una fase dominada por un fluido rígido.
     En los escenarios de ley de potencia rota (Materia/Rígido y Radiación/Rígido), el espectro depende de un parámetro adicional, $w_{\text{sd}}$, que representa la EoS durante el episodio rígido, donde $w_{\text{reh}} \leq w_{\text{sd}} \leq 1$.
3. Restricciones y pronósticos: Los espectros de GW derivados se prueban contra:
   • Límites de BBN: La densidad de energía integrada de las GW está restringida por el límite $\Delta N_{\text{eff}} \lesssim 0.4$.
   • Detectabilidad de GW: Se calcula la Relación Señal-Ruido (SNR) para futuros interferómetros: LISA, Einstein Telescope (ET), DECIGO y BBO, adoptando un criterio de detección de SNR > 10.

Contribuciones clave y resultados

• Exclusión de la región 1$\sigma$: El análisis demuestra que toda la región $1\sigma$ del espacio de parámetros $(w_{\text{reh}}, T_{\text{reh}})$ que concuerda con los datos de ACT+DESI está excluida por las restricciones de BBN. La $w_{\text{reh}}$ rígida requerida para elevar $n_s$ al valor observado produce inevitablemente un espectro de GW con una amplitud que viola el límite de $\Delta N_{\text{eff}}$. Este resultado es válido para los tres escenarios de recalentamiento considerados.
• Impacto de las leyes de potencia rotas: Contrario a la expectativa ingenua de que un episodio temprano de dominancia de materia o radiación aliviaría las restricciones de BBN mediante un sesgo rojo o un aplanamiento del espectro, los autores encuentran lo opuesto. Debido a que la EoS promedio ($w_{\text{reh}}$) debe permanecer alta para satisfacer los datos de CMB, el episodio rígido subsiguiente debe tener una EoS aún más alta ($w_{\text{sd}}$) para compensar. Esto conduce a un sesgo azul más pronunciado en altas frecuencias, expandiendo la región excluida en el espacio de parámetros en comparación con el caso de rígido puro.
• Región 2$\sigma$ y detectabilidad: Aunque la región $1\sigma$ está descartada, una parte de la región $2\sigma$ permanece permitida por BBN.
  • Detectores de tercera generación (LISA, ET): En el escenario de rígido puro, estos instrumentos solo pueden sondear regiones ya excluidas por BBN. Sin embargo, en los escenarios de dos pasos (Materia+Rígido o Radiación+Rígido), ganan sensibilidad hacia partes de la región $2\sigma$ permitida, sondeando temperaturas de recalentamiento hasta $T_{\text{reh}} \sim 10^3$ GeV.
  • Detectores futuros (DECIGO, BBO): Estos experimentos pueden sondear una fracción significativa de la región $2\sigma$ restante, incluso en el escenario de rígido puro, accediendo a temperaturas de recalentamiento de hasta $T_{\text{reh}} \sim 10^5\text{--}10^6$ GeV y $w_{\text{reh}} \gtrsim 0.35$.
• Reevaluación de la tensión: Cuando se aplica el filtro de exclusión de BBN a las distribuciones posteriores, la tensión entre la inflación de Starobinsky y los datos conjuntos de CMB+BAO persiste al nivel de $>2\sigma$ ($n_s \leq 0.9688$). La inclusión del recalentamiento rígido reduce la tensión de $\sim 2\sigma$ a $\sim 1.4\sigma$, pero no logra el pleno acuerdo con los datos sin violar los límites de BBN.

Significancia
El artículo concluye que, si bien un episodio de recalentamiento dominado por un fluido rígido puede teóricamente desplazar la predicción de Starobinsky para $n_s$ hacia las mediciones recientes de ACT+DESI, lo hace a costa de violar las restricciones de BBN en la región de concordancia de $1\sigma$. Por consiguiente, el modelo no puede reconciliarse con los datos al nivel de $1\sigma$ independientemente de los detalles del recalentamiento.

La significancia del trabajo radica en su demostración de que los futuros observatorios de ondas gravitacionales sirven como sondas críticas para estos escenarios. Específicamente, la no detección del espectro con sesgo azul predicho por DECIGO o BBO restringiría aún más el espacio de parámetros de recalentamiento, potencialmente descartando las regiones restantes permitidas por $2\sigma$. Los autores destacan que sigue existiendo una ventana específica ($T_{\text{reh}} \in (10^7, 10^{10})$ GeV) donde la tensión es $<3\sigma$, los límites de BBN se cumplen y la señal es actualmente indetectable; notablemente, este rango se solapa con la temperatura de recalentamiento predicha por el recalentamiento puramente gravitacional en los modelos de Starobinsky. Así, los experimentos de GW están posicionados para probar definitivamente la viabilidad de estas historias específicas de recalentamiento inflacionario.