Finding the boundary: Using galaxy membership to inform galaxy cluster extent through machine learning

Resumen Técnico: Determinación de la Extensión de los Cúmulos de Galaxias mediante Aprendizaje Automático

1. El Problema

La definición de los límites físicos de los cúmulos de galaxias es un desafío fundamental en la cosmología y la astrofísica. Tradicionalmente, se han utilizado umbrales rígidos basados en la densidad de materia oscura (como $R_{200,mean}$, $R_{200,crit}$) o dinámicos (como el radio de splashback, $R_{sp}$) para separar las galaxias de campo de las galaxias miembros del cúmulo. Sin embargo, estos límites no capturan adecuadamente la complejidad de la transición entre ambos entornos.

El entorno del cúmulo influye en la evolución de las galaxias a través de procesos como el ram-pressure stripping (despojo por presión de arrastre), la estrangulación y las interacciones de marea. Estos procesos no ocurren instantáneamente en una frontera definida, sino que generan una región de transición difusa donde las propiedades intrínsecas de las galaxias (estelares, gaseosas y dinámicas) cambian gradualmente. La falta de un método cuantitativo para mapear esta región de transición probabilística limita la comprensión de la "pre-procesamiento" de las galaxias antes de entrar en el núcleo del cúmulo.

2. Metodología

Los autores emplean un marco de trabajo impulsado por datos utilizando Redes Neuronales Profundas (DNN) supervisadas para clasificar galaxias y mapear la extensión espacial de la influencia del entorno del cúmulo.

• Datos: Se utilizaron simulaciones cosmológicas magnetohidrodinámicas de IllustrisTNG, específicamente el volumen más grande y de mayor resolución, TNG300-1.

  • Muestras: Se seleccionaron cúmulos con masas $M_{200,mean} > 10^{13} M_{\odot}$ a $z=0$.
  • Galaxias: Se analizaron subhalos con masa estelar $M_{\star} > 10^9 M_{\odot}$ y masa de materia oscura $M_{dm} > 10^{9.5} M_{\odot}$.
  • Distancia: Se calculó la distancia de cada galaxia al cúmulo más cercano, normalizada por el radio $R_{200,mean}$ del host.

• Etiquetado (Labels):

  • Se definieron dos esquemas de etiquetado binario (0 = cúmulo, 1 = campo):
    1. Basado en distancia: Galaxias dentro de un umbral inferior ($r/R_x < 0.5$) como miembros del cúmulo y fuera de un umbral superior ($r/R_x > 5.0$) como de campo. Se probaron tres definiciones de radio: $R_{sp}$, $R_{200,mean}$ y $R_{200,crit}$.
    2. Basado en árbol de fusión (SUBLINK): Galaxias que alguna vez estuvieron asociadas a un grupo huésped masivo a lo largo de su rama progenitora principal.

• Modelo de Aprendizaje Automático:

  • Se entrenó una Red Neuronal de Perceptrón Multicapa (MLP) totalmente conectada.
  • Entradas: Se utilizaron 15 propiedades intrínsecas de las galaxias, agrupadas en tres categorías físicas: Estelares (masa, color, metalicidad), Gaseosas (fracción de gas, metalicidad, SFR) y Dinámicas (velocidad, dispersión, espín).
  • Preprocesamiento: Se aplicaron transformaciones arcsinh y logarítmicas para simetrizar las distribuciones de las características.
  • Entrenamiento: Se dividió el conjunto de datos por cúmulos (70% entrenamiento, 15% validación, 15% prueba) para evitar la contaminación ambiental entre conjuntos. Se utilizó regularización (Dropout, Batch Normalization, Weight Decay) y optimización con AdamW.

• Análisis de la Transición:

  • El modelo predice la probabilidad de pertenencia al cúmulo.
  • Se identificó la región de transición como el radio donde el logit de salida del modelo cruza cero (probabilidad = 0.5).
  • Se analizaron perfiles apilados y medianas de cero-cruce por bin de masa para estudiar la dependencia de la masa y la dispersión intrínseca.

3. Contribuciones Clave

1. Primer Mapa Probabilístico Cuantitativo: Presentan el primer mapa que define la región de transición entre galaxias de campo y de cúmulo no como una frontera física aguda, sino como una región amplia y dispersa intrínsecamente, basada en la probabilidad de membresía derivada de propiedades físicas.
2. Desacoplamiento de Procesos Físicos: Al entrenar modelos específicos para diferentes categorías de propiedades (gas, estelar, dinámica), logran aislar cómo y dónde ocurren las transiciones asociadas a diferentes mecanismos físicos.
3. Validación de la Naturaleza Difusa: Demuestran que las definiciones de límites tradicionales ($R_{200}$, $R_{sp}$) no separan limpiamente a las poblaciones de galaxias según sus propiedades intrínsecas, revelando una población mixta significativa en los bordes del cúmulo.

4. Resultados Principales

• Naturaleza de la Transición: La transición no es un límite duro. Existe una población "ambigua" (logits cercanos a 0) que se distribuye ampliamente, indicando que las galaxias experimentan un cambio gradual en sus propiedades a medida que se acercan al cúmulo.
• Dependencia de la Masa:
  • El radio de transición ($r_0$) aumenta con la masa del cúmulo ($r_0 \propto M_{200}^{0.1}$).
  • Esta tendencia es inversa a la relación masa-concentración de los halos de materia oscura. Los cúmulos de mayor masa tienen perfiles de densidad más planos en sus regiones exteriores, permitiendo que procesos como el ram-pressure stripping actúen a mayores fracciones del radio virial.
• Diferencias por Tipo de Propiedad:
  • Propiedades Dinámicas: Muestran la transición más profunda (más cerca del centro), ya que la dinámica orbital cambia rápidamente al entrar en el potencial del cúmulo.
  • Propiedades Estelares y Gaseosas:
    • En cúmulos de baja masa, las transiciones de gas y estrellas ocurren en ubicaciones similares.
    • En cúmulos de alta masa, la transición de las propiedades gaseosas ocurre más lejos (en las regiones más externas) que la estelar. Esto se debe a que el gas se despoja o calienta antes de que la formación estelar se apague (retraso de quenching) y la morfología cambie.
• Ineficacia de los Límites Tradicionales: Los modelos entrenados con etiquetas basadas en $R_{200}$ o $R_{sp}$ no logran una separación perfecta de las galaxias, confirmando que estas definiciones geométricas no coinciden con los cambios físicos reales en las galaxias.

5. Significado e Implicaciones

• Revisión de Definiciones de Cúmulos: El estudio sugiere que la definición de los límites de los cúmulos debe pasar de ser un radio fijo a una frontera probabilística que depende de las propiedades físicas de las galaxias y la masa del cúmulo.
• Pre-procesamiento Ambiental: Se confirma que el entorno del cúmulo comienza a afectar a las galaxias mucho antes de que estas alcancen el radio virial o el radio de splashback, especialmente en cúmulos masivos donde el gas se ve afectado en las regiones más externas.
• Aplicabilidad Observacional: Este marco metodológico puede adaptarse a datos observacionales (usando, por ejemplo, redMaPPer para identificar cúmulos y propiedades fotométricas/espectroscópicas) para refinar la selección de miembros de cúmulos y reducir la contaminación por galaxias de campo en estudios de evolución galáctica.
• Futuro: Los autores proponen validar estos hallazgos comparando sus límites basados en propiedades con mediciones independientes de ondas de choque de acreción (efecto SZ) y radios de splashback (lente gravitacional), lo que podría revelar si la transición está dominada por procesos gaseosos o dinámicos.

En conclusión, el artículo establece que la "frontera" de un cúmulo de galaxias es una zona de mezcla física compleja, y el aprendizaje automático ofrece una herramienta robusta para cuantificar esta transición de manera más precisa que los métodos geométricos tradicionales.