Improved Big Bang Nucleosynthesis constraints on decaying massive relics

Resumen Técnico: Mejoras en las Restricciones de la Nucleosíntesis del Big Bang sobre Reliquias Masivas Decaentes

Enunciado del Problema
La Nucleosíntesis del Big Bang (BBN) sirve como una sonda rigurosa de la historia térmica del Universo temprano y de la física más allá del Modelo Estándar (BSM). Específicamente, las reliquias pesadas y de larga vida ($\phi$) producidas en el Universo temprano pueden decaer en partículas del Modelo Estándar (SM) durante o después de la época de la BBN (temperaturas $T \lesssim$ pocos MeV). Estos decaimientos inyectan energía en el plasma primordial, alterando las abundancias de elementos ligeros ($^1$H, D, $^3$He, $^4$He, $^7$Li) mediante mecanismos distintos: modificando la relación neutrón-protón ($n/p$) a través de interconversiones hadrónicas, desintegrando núcleos ligeros por hadrodisintegración y desintegrando núcleos mediante cascadas electromagnéticas. Análisis anteriores han proporcionado restricciones sobre tales reliquias, pero a menudo se basaron en mediciones de abundancia desactualizadas, tasas de reacción nuclear simplificadas o aproximaciones en el tratamiento de la radiación del estado final y la hadronización. Además, la interacción entre las interconversiones tempranas y los procesos de desintegración tardía requería un marco teórico más unificado.

Metodología
Los autores presentan una actualización exhaustiva de las restricciones de la BBN para reliquias pesadas ($m_\phi \gtrsim \text{pocos GeV}$) que decaen en varios canales SM de dos cuerpos (excluyendo neutrinos, que se reservan para trabajos futuros). El análisis se basa en una cadena de cálculo significativamente refinada:

1. Desarrollo de Código (AlterAlterBBN & ACROPOLIS): Los autores adaptaron el código de BBN AlterBBN a una versión extendida, AlterAlterBBN, capaz de manejar cosmologías de fondo arbitrarias e inyectar material hadrónico. Este código rastrea las interconversiones $p \leftrightarrow n$ y la hadrodisintegración durante la época activa de la BBN. Para la evolución tardía ($t \gtrsim 10^4$ s), donde la fotodesintegración se vuelve dominante, el código se interconecta con ACROPOLIS. Un logro técnico crítico es la coincidencia perfecta de estos dos regímenes a $T \approx 5$ keV, asegurando un tratamiento consistente de la fusión nuclear y la desintegración sin doble conteo.
2. Espectros de Inyección Hadrónica (PYTHIA 8): Para modelar con precisión la composición de los productos de decaimiento, los autores utilizan PYTHIA 8 para simular la radiación del estado final (FSR) y la hadronización. Esto proporciona entradas precisas para el número de mesones (piones, kaones) y nucleones inyectados, sus energías cinéticas promedio y la fracción de energía depositada en el sector electromagnético ($\zeta_{EM}$). Este enfoque reemplaza las antiguas aproximaciones de espectros "universales" menos precisos.
3. Aproximación de Equilibrio Dinámico: Los autores implementan un tratamiento refinado de las abundancias hadrónicas utilizando una aproximación de equilibrio dinámico. Dado que las tasas de interacción hadrónica (dispersión, aniquilación) son órdenes de magnitud más rápidas que la tasa de expansión de Hubble durante la BBN, las abundancias de mesones inyectados y antinucleones se resuelven algebraicamente como función de los términos fuente y las tasas de agotamiento. Esto permite el seguimiento eficiente de cadenas de interconversión complejas que involucran piones, kaones ($K^\pm, K_L^0$) y pares nucleón-antinucleón.
4. Entradas de Física Actualizadas:
   • Datos Observacionales: El análisis incorpora las últimas mediciones de abundancia primordial, notablemente la fracción de masa actualizada de $^4$He de la Ref. [31] y las relaciones revisadas de deuterio a hidrógeno.
   • Tasas Nucleares: Se implementan tasas de reacción nuclear actualizadas para la BBN estándar y nuevas secciones eficaces dependientes de la temperatura para interconversiones hadrónicas (específicamente $p \leftrightarrow n$ mediante kaones y piones) de la Ref. [57].
   • Tratamiento de Errores: Las incertidumbres teóricas en las tasas nucleares, las secciones eficaces de interconversión y las eficiencias de hadrodisintegración (parámetros $\xi$) se propagan utilizando un enfoque tipo Monte Carlo (ejecutando escenarios de tasa alta, media y baja) combinado con errores observacionales.
5. Cosmología de Fondo: El código cuenta con las modificaciones a la tasa de Hubble y a la temperatura de desacoplamiento de neutrinos causadas por la densidad de energía de la reliquia. También aborda la abundancia de "congelación" (freeze-in) irreducible de reliquias producidas mediante decaimientos inversos, estableciendo un límite físico inferior para la abundancia de la reliquia $Y_\phi$.

Contribuciones Clave

• Marco Unificado: El artículo unifica el tratamiento de las interconversiones hadrónicas (tiempos tempranos), la hadrodisintegración (tiempos intermedios) y la fotodesintegración (tiempos tardíos) en una sola cadena de cálculo coherente, corrigiendo limitaciones anteriores donde estos regímenes se trataban por separado o con límites inconsistentes.
• Tratamiento Hadrónico Refinado: La inclusión de interconversiones inducidas por kaones y el uso de PYTHIA 8 para FSR y hadronización proporcionan una descripción más precisa de los espectros de inyección, particularmente para canales de decaimiento electromagnético (por ejemplo, $\phi \to e^+e^-$) que anteriormente estaban subrestringidos en cuanto a efectos hadrónicos.
• Corrección de Errores Sistemáticos: Los autores corrigen un error en la curva de pérdida de energía electromagnética identificado en trabajos anteriores, que había aumentado artificialmente las tasas de hadrodisintegración por protones energéticos. Esta corrección suaviza las restricciones en regiones específicas de parámetros.
• Cobertura Exhaustiva de Canales: El estudio cubre una amplia gama de canales de decaimiento de dos cuerpos ($e^+e^-, u\bar{u}, \gamma\gamma, \mu^+\mu^-, \tau^+\tau^-, b\bar{b}, t\bar{t}, gg, WW, ZZ, hh$), proporcionando contornos de exclusión para cada uno.

Resultados
Los autores presentan contornos de exclusión actualizados en el plano $(m_\phi Y_\phi, \tau_\phi)$ para varios canales de decaimiento.

• Canales de Decaimiento:
  • Electromagnéticos ($e^+e^-, \gamma\gamma$): Las restricciones se ven impulsadas por la sobreproducción de $^4$He a vidas cortas (vía interconversiones) y la sobreproducción de D/$^3$He a vidas más largas (vía desintegración). La inclusión de hadrones inducidos por FSR endurece significativamente los límites para los decaimientos $e^+e^-$ en comparación con la literatura anterior.
  • Hadrónicos ($u\bar{u}, gg$): Estos canales producen las restricciones más fuertes debido a la inyección de energía hadrónica no suprimida. Para $m_\phi = 10$ GeV, las restricciones son tan severas que excluyen incluso la abundancia de congelación irreducible para vidas $0.1 \text{ s} \lesssim \tau_\phi \lesssim 10^5 \text{ s}$.
  • Leptones ($\mu^+\mu^-, \tau^+\tau^-$): El canal $\tau^+\tau^-$ muestra características únicas debido a la gran rama de decaimiento en piones, lo que conduce a fuertes límites de interconversión a bajas masas.
• Dependencia de la Masa: Para reliquias pesadas ($m_\phi \gtrsim 1$ TeV), las restricciones tienden a volverse independientes de la masa a vidas largas debido a la universalidad del espectro de cascada electromagnética. Sin embargo, en el régimen dominado por hadrones, las restricciones pueden debilitarse ligeramente para masas muy altas debido al crecimiento sublineal en el número de hadrones inyectados.
• Comparación con la Literatura: Los resultados coinciden cualitativamente en general con estudios anteriores (por ejemplo, Kawasaki et al.), pero muestran diferencias cuantitativas significativas. A vidas cortas, los nuevos límites son mucho más fuertes, principalmente debido a las tasas de interconversión actualizadas y la inclusión adecuada de hadrones inducidos por FSR. A vidas largas, las diferencias surgen del tratamiento refinado de la fotodesintegración y la abundancia actualizada de $^4$He.

Significado
El artículo afirma proporcionar un tratamiento de vanguardia de las restricciones de la BBN sobre reliquias pesadas decayeras. Al integrar datos observacionales actualizados, secciones eficaces nucleares precisas y una simulación sofisticada de lluvias de partículas, los autores establecen límites de exclusión robustos y mejorados. Estos resultados son particularmente significativos para restringir modelos BSM que involucran estados pesados y de larga vida, ya que las restricciones actualizadas pueden descartar espacios de parámetros que anteriormente se consideraban viables, incluidas regiones donde la abundancia de la reliquia es tan baja como la producción de congelación irreducible. El trabajo sirve como una referencia crítica para interpretar los datos de la BBN en el contexto de nueva física, mientras señala explícitamente que el tratamiento de los canales de neutrinos y las reliquias más ligeras (escala de MeV) sigue siendo un tema para futuras investigaciones.