Time-Dependent Cosmic Ray Halos from Bursty Star Formation and Active Galactic Nuclei: Semi-Analytic Formalism and Galaxy Formation Implications

Resumen Técnico: Halos de Rayos Cósmicos Dependientes del Tiempo provenientes de la Formación Estelar Intermitente y de Núcleos Galácticos Activos

Planteamiento del Problema
La retroalimentación de los rayos cósmicos (CR, por sus siglas en inglés) es cada vez más reconocida como un mecanismo crítico en la evolución de las galaxias, particularmente para regular la formación estelar en galaxias masivas ($M_{\text{halo}} > 10^{13} M_{\odot}$) mediante Núcleos Galácticos Activos (AGN) y supernovas (SNe). Sin embargo, persisten incertidumbres significativas con respecto al transporte de los CR a través del medio interestelar y el medio circumgaláctico (ISM y CGM). Los esfuerzos de modelado actuales se enfrentan a una dicotomía: las simulaciones completas de magnetohidrodinámica de rayos cósmicos (CR-MHD) son computacionalmente prohibitivas para volúmenes cosmológicos grandes, mientras que los tratamientos analíticos y semianalíticos existentes suelen basarse en supuestos de estado estacionario. Estos supuestos de estado estacionario, aunque válidos para escenarios de inyección continua como la Vía Láctea, pueden fallar al capturar la física de las historias de formación estelar (SFH) intermitentes y la acreción episódica de AGN, las cuales son características de la formación de galaxias de alto desplazamiento al rojo y de galaxias masivas. Los autores argumentan que omitir la dependencia temporal y los efectos de tiempo de viaje finito puede conducir a predicciones inexactas de los perfiles de presión de los CR ($P_{\text{CR}}$) y sus gradientes, los cuales impulsan la retroalimentación dinámica.

Metodología
Los autores desarrollan un formalismo semianalítico y un marco numérico para modelar la evolución dependiente del tiempo de la presión de los CR en los halos galácticos, relajando el supuesto de inyección en estado estacionario.

1. Marco Analítico: El estudio resuelve la ecuación de difusión-advección 3D para la presión de los CR, tratando el transporte en el límite de fluido relativista. Los autores parametrizan la ignorancia de la dispersión microfísica mediante coeficientes de difusión efectiva ($\kappa_0$) y de streaming/advección ($v_{\text{eff}}$). Derivan soluciones de la función de Green para tres regímenes:

   • Caso 1: Difusión pura ($\kappa_0$ constante, $v_{\text{eff}} \to 0$).
   • Caso 2: Advección/streaming puro ($\kappa_0 \to 0$, $v_{\text{eff}}$ constante).
   • Caso 3: Combinación de difusión y advección/streaming. Para este caso, proponen una solución semianalítica aproximada (Ec. 7) que convoluciona una función fuente dependiente del tiempo con un núcleo de la función de Green espacial, normalizada para conservar la energía total.

2. Implementación Numérica: Los autores implementan soluciones numéricas exactas a la ecuación de transporte utilizando un enfoque de volumen finito en una rejilla radial logarítmica con un método de Runge-Kutta implícito (Radau IIA) para manejar problemas de difusión rígidos (stiff).

3. Historias de Inyección: Los modelos se prueban contra cuatro historias de inyección de CR distintas:

   • Una inyección de función $\delta$ única (pico en $z \sim 3$).
   • Una inyección constante de estado estacionario (promediada sobre un tiempo de Hubble).
   • Una historia de inyección de decaimiento lento motivada empíricamente (basada en el modelo Trinity).
   • Una historia intermitente y dependiente del tiempo derivada de una simulación de zoom de CR-MHD cosmológica (suite FIRE-3).

4. Validación: Los resultados semianalíticos y numéricos se comparan con una simulación de zoom CR-MHD totalmente cosmológica y autoconsistente de una galaxia masiva ($M_{\text{halo}} \approx 10^{13} M_{\odot}$) en $z \approx 1.3$.

Contribuciones Clave y Resultados

• La Dependencia del Tiempo Altera los Perfiles de Presión: El estudio demuestra que la inyección dependiente del tiempo altera significamente los perfiles de $P_{\text{CR}}$ en comparación con las predicciones de estado estacionario. Específicamente, la inyección intermitente o episódica eleva la presión en el halo interno ($r \lesssim 100$ kpc) debido a eventos de inyección tardíos, y aplana los perfiles en radios grandes ($r \gtrsim R_{\text{vir}}$). Este aplanamiento ocurre incluso en historias de inyección suavizadas empíricamente, lo que indica que la evolución lenta de las tasas de acreción promedio influye en las distribuciones de los halos de bajo desplazamiento al rojo.
• Efectos de Tiempo de Viaje Finito: Los autores muestran que los supuestos de estado estacionario, que a menudo asumen un "coeficiente de difusión efectivo" ($\kappa_{\text{eff}} \sim r$) para subsumir el comportamiento de streaming, fallan al capturar los efectos de tiempo de viaje finito. En escenarios de variabilidad temporal, estos supuestos pueden sobrepredecir $P_{\text{CR}}$ en radios grandes (al asumir un tiempo de propagación infinito) o subpredecirla (al perder las colas extendidas de la propagación difusiva), dependiendo del régimen de transporte específico y la historia de inyección.
• Ruptura de la Degeneración: El trabajo sugiere que la dependencia temporal puede romper la degeneración entre los comportamientos de transporte tipo "difusión" y tipo "streaming". Aunque estos regímenes pueden parecer degenerados en las soluciones de estado estacionario, sus predicciones divergen significativamente cuando se consideran la evolución de la fuente y los tiempos de viaje finitos.
• Validez Semianalítica: La aproximación semianalítica propuesta (Ec. 7) reproduce las soluciones numéricas "reales" y los resultados de la simulación CR-MHD con una precisión sorprendente, particularmente en radios grandes ($r > R_{\text{vir}}$). Aunque ocurren sobreestimaciones leves en radios intermedios y subestimaciones en radios internos, el enfoque captura las características de transporte principales que pasan desapercibidas para las formulaciones de estado estacionario.
• Impacto en los Gradientes: La evolución dependiente del tiempo conduce a gradientes de presión ($\nabla P_{\text{CR}}$) fuertemente variables que difieren de las expectativas de estado estacionario. Dado que estos gradientes impulsan los vientos y la retroalimentación inducida por CR, los autores advierten que los modelos de sub-malla que ignoran la dependencia temporal pueden predecir incorrectamente la magnitud y la ubicación de los efectos dinámicos.

Significancia y Reivindicaciones
El artículo afirma que la consideración cuidadosa de los efectos de "suavizado" dependientes del tiempo es esencial para las rutinas de sub-malla en simulaciones cosmológicas de gran volumen. Los autores argumentan que los modelos de sub-malla de estado estacionario actuales (por ejemplo, los de Hopkins et al.) pueden sobrepredecir las contribuciones de presión de los CR de fuentes de edad finita, mientras que simultáneamente omiten o desplazan los gradientes de presión fuertes. Esto podría conducir a errores no lineales en la predicción del crecimiento de las galaxias, las condiciones de salida (outflow) y la regulación de la formación estelar en halos masivos.

Los autores proponen su formalismo semianalítico como una vía computacionalmente eficiente para implementar la retroalimentación de CR dependiente del tiempo en volúmenes cosmológicos, permitiendo potencialmente la inclusión de contribuciones de AGN e inyecciones intermitentes de formación estelar sin el costo prohibitivo de una CR-MHD completa. Además, sugieren que las características distintivas derivadas de la inyección intermitente transportada mediante streaming/advección (como los aumentos transitorios de presión en radios grandes) podrían vincularse con fenómenos observables como los "Círculos de Radio Impares" (ORCs), ofreciendo un posible control observacional para constreñir los mecanismos de transporte de CR en el CGM exterior. En última instancia, este trabajo motiva un cambio hacia el modelado dependiente del tiempo para comprender mejor los efectos emergentes de la dispersión de los CR en la formación de galaxias a través de las épocas cósmicas.