MIGHTEE-HI: The direct detection of neutral hydrogen in galaxies at $z>0.25$

1. El Problema

El hidrógeno neutro atómico (HI) es el reservorio de gas primordial del cual se forma el gas molecular y, posteriormente, las estrellas en las galaxias. Comprender su relación con la evolución galáctica a lo largo del tiempo cósmico es crucial, pero ha sido extremadamente difícil de observar directamente más allá del universo local.

• Limitaciones actuales: La detección de la línea de emisión de 21 cm del HI es difícil debido a su baja probabilidad de transición en comparación con otras líneas (como Lyman-α). Históricamente, las encuestas no dirigidas (blind) han sido sensibles solo a $z \lesssim 0.1$.
• Desafíos en alto desplazamiento al rojo: A $z > 0.25$, la línea de 21 cm se desplaza a frecuencias donde la detección es aún más desafiante. Los telescopios de plato único (como FAST) han logrado detecciones, pero su baja resolución angular puede llevar a sobrestimaciones de masa debido a la confusión con galaxias vecinas. Los interferómetros ofrecen mejor resolución, pero la sensibilidad necesaria para detectar galaxias individuales a estas distancias ha sido un obstáculo hasta la llegada de telescopios modernos como MeerKAT.

2. Metodología

El estudio presenta una búsqueda dirigida (targeted) de emisión de HI en galaxias individuales utilizando los datos del MIGHTEE (MeerKAT International Tiered Extragalactic Exploration) y el catálogo espectroscópico de DESI (Dark Energy Spectroscopic Instrument).

• Datos de Radio: Se utilizó el cubo espectral de la "Data Release 1" (DR1) de MIGHTEE sobre el campo COSMOS. Este cubo tiene un rango espectral de 960–1150 MHz, un ancho de canal de 104.5 kHz y una resolución angular de ~15 arcsec.
• Selección de Objetivos: Se cruzaron las posiciones y redshifts espectroscópicos ópticos de galaxias seleccionadas por líneas de emisión del catálogo DESI (rango $0.23 < z < 0.5$) con el cubo de radio. Se extrajeron "cubelets" (subcubos) de 80x80 arcsec² centrados en cada galaxia óptica.
• Proceso de Detección:
  1. Búsqueda Visual: Se colapsaron los cubos 3D a imágenes 2D y se examinaron espectros 1D para identificar candidatos a líneas de emisión. Se requirió una señal prominente en la imagen colapsada (dentro de 4 arcsec de la posición óptica) y emisión significativa en el dominio espectral (al menos 5 canales).
  2. Validación Estadística: Se midió la relación señal-ruido (SNR) integrando el flujo sobre la línea. El ruido se estimó midiendo ~500 regiones aleatorias cercanas a la fuente para capturar la variación espacial y espectral intrínseca. Se aplicó una prueba de Anderson-Darling para confirmar la naturaleza gaussiana del ruido.
  3. Criterios de Selección: Se retuvieron galaxias con un SNR integrado > 4.
  4. Controles de Sesgo: Se realizó una búsqueda de control en posiciones desplazadas (50 arcmin) para verificar falsos positivos. También se verificó la presencia de emisión de radio continuo para descartar artefactos de sustracción.
• Análisis de Propiedades: Para las 11 galaxias confirmadas, se realizaron ajustes de Distribución de Energía Espectral (SED) utilizando el código Bagpipes con datos multibanda (desde UV hasta infrarrojo lejano) para derivar masas estelares ($M_*$) y tasas de formación estelar (SFR).
• Relación Tully-Fisher: Se midieron los anchos de línea ($W_{50}$), se corrigieron por inclinación (usando la elipticidad en banda $g$) y se construyó la relación de Tully-Fisher Bariónica (bTFr).

3. Contribuciones Clave

• Primera Muestra Directa: Este trabajo presenta la primera muestra de 11 detecciones directas de HI en galaxias individuales a $z > 0.25$ utilizando un interferómetro (MeerKAT). La detección más lejana se encuentra en $z = 0.3841$.
• Superación de Sesgos: Al utilizar redshifts espectroscópicos ópticos conocidos, el estudio supera las dificultades de las búsquedas no dirigidas, permitiendo detectar galaxias que de otro modo serían demasiado débiles para ser encontradas en un escaneo ciego.
• Medición de bTFr a Alto Redshift: Se realiza la primera medición de la relación de Tully-Fisher bariónica (bTFr) utilizando HI a $z > 0.25$, permitiendo comparar directamente con la relación local sin depender de trazadores de gas diferentes (como CO o H$\alpha$).
• Caracterización de Galaxias Ricas en Gas: Se derivan propiedades físicas completas (masa estelar, SFR, masa de HI) para una muestra de galaxias en un volumen cósmico significativo, llenando un vacío en la comprensión de la evolución del gas neutro.

4. Resultados Principales

• Detecciones: Se identificaron 11 galaxias con HI en el rango $0.25 \lesssim z \lesssim 0.39$. No se detectaron galaxias en el rango $0.39 < z < 0.5$, lo cual se atribuye a la distribución de redshift del catálogo padre DESI (menos galaxias objetivo en ese rango) y a la profundidad de los datos.
• Masas de HI: Las galaxias detectadas tienen masas de HI significativamente mayores ($\log(M_{HI}/M_\odot) \approx 10.2 - 10.8$) en comparación con sus contrapartes seleccionadas por HI a bajo redshift ($z < 0.1$) para una masa estelar dada.
  • Interpretación: Esto no implica necesariamente una evolución en la función de masa de HI. Se explica por el mayor volumen cósmico probed a estos redshifts, que permite detectar la cola superior (más brillante) de la función de masa, que es rara en el universo local y requiere grandes áreas de sondeo para ser encontrada.
• Secuencia Principal de Formación Estelar: La mayoría de las galaxias se sitúan sobre o por encima de la secuencia principal de formación estelar (SFMS) a $z \sim 0.35$, lo cual es consistente con la selección basada en galaxias con líneas de emisión (alta formación estelar).
• Relación Tully-Fisher Bariónica (bTFr):
  • Las galaxias son consistentes con la relación bTFr local, pero muestran una tendencia a situarse por debajo de la relación local en el régimen de altas masas bariónicas y altas velocidades de rotación.
  • Esto sugiere un posible "aplanamiento" de la relación a altas masas/redshifts.
  • Explicaciones posibles: Podría indicar una evolución real en la relación masa bariónica-masa gravitacional (alineado con simulaciones hidrodinámicas como simba), o podría deberse a que la masa de gas molecular (no medida directamente) es significativa en estas galaxias masivas. Si se asume que la masa molecular es comparable a la de HI, los puntos se alinean mejor con la relación local.
• Tiempo de Agotamiento del Gas: Los tiempos de agotamiento ($t_{dep} = M_{HI}/SFR$) tienden a ser menores que en muestras de bajo redshift, aunque con una gran dispersión.

5. Significado e Impacto

• Validación de MeerKAT: El estudio demuestra que MeerKAT, combinado con datos ópticos profundos, es una herramienta poderosa para detectar HI en galaxias individuales más allá del universo local, superando las limitaciones de los telescopios de plato único en términos de resolución y precisión de masa.
• Preparación para el SKA: Este trabajo sienta las bases para futuros estudios del HI en el universo distante. Muestra que las búsquedas dirigidas son viables y que, con la sensibilidad del Square Kilometre Array (SKA), será posible realizar búsquedas ciegas y mapear la evolución de la función de masa de HI y la bTFr hasta $z \sim 0.5$ y más allá.
• Comprensión de la Evolución Galáctica: Proporciona evidencia observacional clave sobre cómo el reservorio de gas neutro evoluciona en galaxias masivas, sugiriendo que la relación entre el gas y la dinámica galáctica puede ser más compleja a altas masas, posiblemente influenciada por la fracción de gas molecular.

En resumen, este artículo marca un hito en la astronomía de radio al establecer la primera muestra robusta de galaxias ricas en HI a $z > 0.25$, validando metodologías híbridas (óptico-radio) y abriendo la puerta a una nueva era de estudios de evolución galáctica impulsados por el gas neutro.