Modelling the expulsion of baryons from haloes from first principles: the role of feedback and of the cosmological constant

Aquí presento un resumen técnico detallado del artículo "Modelling the expulsion of baryons from haloes from first principles: the role of feedback and of the cosmological constant" (Modelado de la expulsión de bariones desde halos desde primeros principios: el papel de la retroalimentación y de la constante cosmológica), publicado en MNRAS (2026).

1. Planteamiento del Problema

La formación de galaxias y la estructura a gran escala dependen críticamente de cómo los procesos astrofísicos (como la retroalimentación estelar y de AGN) interactúan con el modelo cosmológico subyacente para expulsar bariones de los halos de materia oscura.

• El problema de los bariones perdidos: Se observa que los halos de masa inferior a $\sim 10^{12.8} M_\odot$ contienen una fracción de bariones significativamente menor que la fracción cósmica ($f_b = \Omega_b/\Omega_m$).
• El radio de cierre ($R_{closure}$): Se define como la distancia desde el centro del halo donde la fracción de bariones encerrada vuelve a igualar la fracción cósmica. Estudios previos (ej. Ayromlou et al. 2023) han encontrado una relación empírica entre este radio y la masa del halo, pero carecen de una explicación física derivada de primeros principios.
• Limitaciones actuales: Los modelos existentes son puramente empíricos, no explican la dependencia cosmológica (específicamente el papel de la energía oscura o $\Lambda$) y es difícil aislar los efectos de la retroalimentación de los parámetros cosmológicos en simulaciones numéricas estándar.

2. Metodología

Los autores desarrollan un modelo analítico basado en primeros principios y lo validan contra un conjunto de simulaciones hidrodinámicas cosmológicas.

A. Modelo Analítico

El modelo se basa en el equilibrio entre la gravedad, la presión hidrostática y la expansión cósmica (energía oscura).

1. Perfil de densidad de gas: Asumiendo simetría esférica, se deriva un perfil de densidad de gas fuera del radio virial ($R_{200}$) resolviendo la ecuación de equilibrio de fuerzas. En el límite de baja densidad ($\rho \ll \rho_{crit}$), la solución es una ley de potencia: $\rho(r) \propto r^{-\eta}$.
2. Derivación del radio de cierre: Se integra el perfil de densidad hasta que la masa bariónica encerrada iguala la fracción cósmica. Esto lleva a una expresión analítica para el radio de cierre normalizado ($R_{closure}/R_{200}$) que depende de:
   • La fracción de gas en el halo ($f_{gas}$).
   • El índice de la ley de potencia de la densidad ($\eta$).
   • Parámetros cosmológicos ($\Omega_m, \Omega_\Lambda, H_0$).
   • El desplazamiento al rojo ($z$).
   • La fracción de bariones del halo ($f_{b-halo}$).

La ecuación resultante (Eq. 21 en el texto) es:
$$\frac{R_{closure}}{R_{200}} = \left[ \xi(w_n, n) (f_{b-cosmic} - f_{b-halo}) \frac{M_{200}}{M_{200}} + R_{200}^{\frac{2}{n+3}} \right]^{\frac{n}{2+3n}}$$
(Nota: La forma simplificada y final se presenta en la ecuación 21 del texto, mostrando la dependencia explícita de $\Lambda$ y $z$).

B. Validación Numérica

• Simulaciones: Se utilizaron seis variantes de las simulaciones EAGLE (Evolution and Assembly of GaLaxies and their Environments) con diferentes valores de la constante cosmológica: $\Lambda = 0$ (EdS), $\Lambda_0$ (nuestro universo), $10\Lambda_0$, $30\Lambda_0$ y $100\Lambda_0$.
• Procedimiento: Se calcularon los perfiles de fracción bariónica encerrada para halos con $M_{200} > 10^{11} M_\odot$ en el rango de desplazamiento al rojo $0 < z < 3$.
• Ajuste: Se ajustaron los parámetros libres del modelo analítico ($\eta$ y $f_{gas}$) minimizando el chi-cuadrado ($\chi^2$) frente a los datos de las simulaciones.

3. Contribuciones Clave

1. Primera explicación física: Proporciona la primera derivación analítica de la relación entre el radio de cierre y la masa/fracción bariónica del halo, sin recurrir a ajustes empíricos puros.
2. Generalización cosmológica: El modelo es válido para un rango amplio de constantes cosmológicas ($0 \le \Lambda/\Lambda_0 \le 100$), demostrando robustez frente a variaciones en la energía oscura.
3. Conexión con modelos previos: Demuestra que la fórmula empírica de Ayromlou et al. (2023) es un caso límite de su modelo analítico completo bajo ciertas aproximaciones (linealización y límites de alto/bajo $z$).
4. Cuantificación del impacto de $\Lambda$: Aísla y cuantifica el papel de la energía oscura en la expulsión de bariones, diferenciándolo de la retroalimentación astrofísica.

4. Resultados Principales

• Acuerdo con simulaciones: El modelo analítico reproduce con excelente precisión los resultados de las simulaciones EAGLE para halos masivos ($M_{200} > 10^{11} M_\odot$) en el rango $0 < z < 3$.
• Evolución de parámetros:
  • El índice de densidad $\eta$ y la fracción de gas $f_{gas}$ muestran una evolución con el desplazamiento al rojo.
  • En el modelo linealizado (aproximación de primer orden), los parámetros ajustados $\eta_{approx} \approx 1.8$ y $f_{gas-approx} \approx 0.002$ resultan ser constantes a través de diferentes desplazamientos al rojo y valores de $\Lambda$. Esto sugiere que las diferencias en el radio de cierre entre cosmologías se deben principalmente a los parámetros cosmológicos ($\Omega_m, \Omega_\Lambda$) y no a cambios en las propiedades termodinámicas intrínsecas del gas.
• Impacto de la Energía Oscura:
  • Un aumento en $\Lambda$ reduce el radio crítico ($R_{crit}$) y hace que el perfil de densidad de gas sea menos concentrado (menor $\eta$), facilitando que los bariones sean expulsados más allá del radio virial.
  • Cuantificación: En el modelo cosmológico fiducial ($\Lambda = \Lambda_0$), el radio de cierre a $z < 2$ es aproximadamente un 30% mayor que en un universo Einstein-de Sitter ($\Lambda=0$). Esto implica que la energía oscura contribuye con un ~30% a la expulsión de bariones, siendo el resto dominado por la retroalimentación astrofísica.
  • En cosmologías con $\Lambda \ge 10\Lambda_0$, la energía oscura se convierte en el factor dominante en la evacuación de bariones.
• Relación con la fórmula A23: La fórmula empírica de Ayromlou et al. (2023) funciona bien en nuestro universo y a $z$ fijos, pero falla o requiere re-calibración en universos con $\Lambda$ alto. El modelo analítico generaliza esta relación.

5. Significado e Implicaciones

• Resolución del problema de los bariones perdidos: El modelo ofrece una explicación física coherente de por qué los bariones "faltantes" se encuentran a grandes distancias de los halos, integrando tanto la retroalimentación como la expansión acelerada del universo.
• Argumentos Antropicos: Los resultados tienen implicaciones profundas para el principio antrópico. En universos con una constante cosmológica significativamente mayor que la nuestra, la energía oscura expulsaría tan eficientemente los bariones que la formación de estrellas y galaxias se vería severamente suprimida. Esto limitaría la formación de elementos pesados y, por ende, la posibilidad de vida compleja.
• Herramienta predictiva: El modelo proporciona una herramienta robusta para predecir la distribución de bariones en diferentes entornos cosmológicos y épocas, útil para interpretar futuras observaciones de instrumentos como JWST y ELT, que podrán mapear la distribución de gas en el medio circumgaláctico e intergaláctico.

En conclusión, el trabajo establece un marco teórico sólido que unifica la física de la retroalimentación galáctica y la cosmología para explicar la distribución de la materia bariónica, demostrando que la energía oscura juega un papel no despreciable, y en algunos casos dominante, en la evolución de los halos de materia oscura.