Cosmology with one galaxy: An analytic formula relating $Ω_{\rm m}$ with galaxy properties

Resumen Técnico: Cosmología con una Galaxia

Enunciado del Problema
Los análisis cosmológicos estándar suelen desacoplar la formación de galaxias de la inferencia de parámetros cosmológicos, recurriendo a estadísticas a nivel de población como el agrupamiento, el lente gravitacional o las abundancias de halos para restringir el parámetro de densidad de materia, $\Omega_m$. Aunque estudios recientes de aprendizaje automático sobre la suite de simulaciones CAMELS han demostrado que las propiedades de galaxias individuales pueden codificar información cosmológica, el origen físico de esta señal ha permanecido incierto, y los enfoques anteriores dependían de redes neuronales de "caja negra" carentes de interpretabilidad física. Este trabajo aborda el desafío de derivar una relación analítica físicamente transparente que vincule $\Omega_m$ directamente con observables intrínsecos de galaxias individuales, evitando la necesidad de estadísticas de agrupamiento a gran escala.

Metodología
Los autores utilizan el proyecto CAMELS (Cosmología y Astrofísica con Simulaciones de Aprendizaje Automático), específicamente los conjuntos de datos de Hipercubo Latino (LH) y el conjunto de datos extendido TNG-SB28. Estas suites contienen más de 16.000 simulaciones hidrodinámicas que varían parámetros cosmológicos ($\Omega_m$, $\sigma_8$, $\Omega_b$, $H_0$, etc.) y parámetros astrofísicos subgrid (eficiencias de retroalimentación de supernovas y AGN).

La metodología central emplea Regresión Simbólica (utilizando el algoritmo PySR) para descubrir formas funcionales analíticas compactas que relacionan las propiedades de las galaxias con $\Omega_m$.

1. Preparación de Datos: El conjunto de entrenamiento consiste en galaxias muestreadas de simulaciones IllustrisTNG en $z=0$. Se construyen variables clave adimensionales, incluyendo la metalicidad estelar ($Z_\star$), la fracción de masa estelar ($f_\star$), la velocidad dinámica efectiva ($V_{\text{eff}}$) y un parámetro de compacidad ($R_{\text{compact}}$).
2. Búsqueda Simbólica: El algoritmo busca una expresión que minimice el error cuadrático medio entre el $\Omega_m$ predicho y el real, sujeto a penalizaciones por parsimonia y restricciones sobre la complejidad de los operadores (por ejemplo, limitando la profundidad de anidamiento y los rangos de exponentes).
3. Validación: La fórmula analítica resultante, entrenada exclusivamente en IllustrisTNG, se prueba contra suites de simulación independientes (ASTRID, SIMBA, Swift-EAGLE) y a través de diferentes corrimientos al rojo ($z=0$ a $z=3$).
4. Espacio de Parámetros Extendido: El marco se prueba adicionalmente en el conjunto de datos TNG-SB28, que varía 28 parámetros simultáneamente, incluidos $\Omega_b$ y $H_0$, para evaluar la robustez bajo variaciones cosmológicas y de retroalimentación más amplias.

Contribuciones Clave

1. Primer Rastreador Analítico de Galaxia Única: El artículo presenta la primera fórmula analítica físicamente interpretable que relaciona $\Omega_m$ directamente con propiedades de galaxias individuales. La expresión derivada toma la forma:
   $$\Omega_m \approx \ln \left[ \sigma \left( \frac{Z_\star}{M_\star} \frac{\Omega_b V_{\text{eff}}^3}{G H_0} \frac{\Omega_b}{k} \right) + c_0 \right] - \frac{a_0}{R_{\text{compact}}}$$
   donde $\sigma(x)$ es una función sigmoide.
2. Robustez Transversal a Suites: La forma funcional derivada de IllustrisTNG recupera con éxito $\Omega_m$ en las simulaciones ASTRID, SIMBA y Swift-EAGLE con solo una recalibración menor de coeficientes numéricos ($k, c_0, a_0$). Esto sugiere que el acoplamiento físico subyacente es universal entre diferentes implementaciones numéricas de la física de formación de galaxias.
3. Interpretación Física: Los autores proporcionan una interpretación física transparente de la variable latente. El término dentro de la sigmoide codifica un acoplamiento barion-cosmología mediado por la Relación Masa-Metalicidad (MZR). Representa la interacción entre la profundidad del potencial gravitatorio (determinada por $V_{\text{eff}}$ y la cosmología) y la eficiencia de retención de bariones y enriquecimiento químico regulada por la retroalimentación.
4. Evolución con el Corrimiento al Rojo: Los autores derivan una extensión dependiente del corrimiento al rojo de la fórmula utilizando escalas de ley de potencia de $(1+z)$, lo que mantiene la precisión predictiva hasta $z=3$ para IllustrisTNG y ASTRID.
5. Generalización a Espacios de Parámetros de Alta Dimensión: En la suite TNG-SB28 (donde varían $\Omega_b$ y $H_0$), la forma analítica se adapta para incluir una dependencia explícita de la relación del parámetro de Hubble ($H_0/H_{0,\text{fid}}$) y un término refinado de eficiencia de retención $E(Z_\star, f_\star)$, consistente con un marco de regulador de gas.

Resultados

• Rendimiento Predictivo: En $z=0$, el modelo alcanza coeficientes de determinación ($R^2$) de 0.77 para IllustrisTNG y 0.73 para ASTRID. SIMBA y Swift-EAGLE muestran un rendimiento ligeramente degradado ($R^2 \approx 0.6-0.7$) pero aún exhiben una fuerte correlación monótona, lo que indica que la señal no es un artefacto de sobreajuste a una simulación específica.
• Comportamiento de Saturación: En el conjunto de datos de alta dimensión TNG-SB28, el predictor exhibe un régimen de saturación a altos $\Omega_m$ ($\Omega_m \gtrsim 0.43$). Los autores atribuyen esto a fuertes degeneraciones entre la física de retroalimentación y la cosmología en el espacio de parámetros de alta dimensión, donde los cambios en la cosmología son cada vez más compensados por procesos bariónicos, creando un techo de información.
• Mecanismo Físico: El análisis confirma que la "rodilla" de la relación $\Omega_m$–metalicidad corresponde a la transición entre enriquecimiento dominado por cosmología (baja metalicidad) y enriquecimiento saturado por retroalimentación (alta metalicidad). La pendiente de esta transición está gobernada por la prescripción de retroalimentación, lo que explica por qué diferentes simulaciones requieren diferentes calibraciones de coeficientes.

Significado y Afirmaciones
El artículo afirma establecer un puente directo e interpretable entre la física de galaxias a pequeña escala y la cosmología a gran escala. Al demostrar que $\Omega_m$ está codificado localmente en el ciclo bariónico de galaxias individuales, el trabajo abre una vía complementaria para la inferencia cosmológica que no depende de las estadísticas tradicionales de estructura a gran escala.

Los autores enfatizan que este enfoque no pretende reemplazar las sondas estándar, sino proporcionar un marco novedoso para comprender cómo la información cosmológica sobrevive al procesamiento bariónico. El descubrimiento de una relación analítica compacta sugiere que la huella de la densidad de materia en la formación de galaxias es robusta y está gobernada por principios fundamentales de retención de bariones y profundidad del potencial gravitatorio. El trabajo destaca que, aunque el aprendizaje automático puede identificar correlaciones, la regresión simbólica es esencial para descubrir las leyes físicas subyacentes que gobiernan la conexión galaxia-cosmología.