RotCurves: A PYTHON package for efficient modelling and fitting of galactic rotation curves at high-z

Resumen Técnico: RotCurves – Un paquete de Python para el modelado y ajuste eficiente de curvas de rotación galáctica en alto desplazamiento al rojo ($z$)

Planteamiento del Problema
Las curvas de rotación son fundamentales para comprender la estructura, la dinámica y los halos de materia oscura de las galaxias. Sin embargo, analizar estas curvas en altos desplazamientos al rojo ($z$) presenta desafíos significativos debido al oscurecimiento por brillo superficial y a la resolución angular limitada. Estas limitaciones provocan el "ensanchamiento del haz" (beam smearing), un suavizado artificial de los gradientes de velocidad causado por la convolución de la luz de áreas físicas extensas sobre la fuente. Corregir este efecto es crítico pero computacionalmente costoso. Las herramientas existentes de modelado directo paramétrico, como dysmalpy, construyen cubos de datos 3D o 4D completos (incluyendo dimensiones espaciales y espectrales) para simular observaciones. Aunque son precisos, estos métodos son computacionalmente intensivos, requiriendo a menudo segundos o minutos por cada realización, lo que dificulta la exploración rápida de parámetros, la prueba de múltiples supuestos de modelos de masa y el procesamiento de grandes muestras de los próximos sondeos de Unidades de Campo Integral (IFU). Existe la necesidad de una herramienta que equilibre la precisión de la corrección del ensanchamiento del haz con la eficiencia computacional.

Metodología
Los autores presentan RotCurves, una herramienta de modelado directo paramétrico diseñada para ajustar curvas de rotación galáctica de manera eficiente mediante la aproximación del proceso de ensanchamiento del haz en 2D en lugar de 3D.

• Construcción del Modelo: El código construye un modelo de galaxia idealizado y axisimétrico en equilibrio centrífugo. Calcula la velocidad circular intrínseca en el plano medio galáctico sumando las contribuciones de varios componentes de masa:
  • Bariones: Modelados como discos (delgados o gruesos) y bulbos utilizando perfiles de Sérsic (incluyendo perfiles exponenciales y de de Vaucouleurs).
  • Materia Oscura: Modelada como halos esféricamente simétricos utilizando perfiles como NFW, $\alpha$-NFW, Burkert, Einasto y Dekel-Zhao. El código incluye una opción para la contracción adiabática del halo debido a la acumulación bariónica.
  • Soporte de Presión: El modelo corrige la velocidad de rotación observada por los gradientes de presión en el medio interestelar (ISM) turbulento, utilizando supuestos de equilibrio hidrostático, convirtiendo la velocidad circular intrínseca en velocidad de rotación observable basada en la dispersión de velocidad ($\sigma_0$).
• Aproximación del Ensanchamiento del Haz: A diferencia de las herramientas que rotan un cubo galáctico 3D para que coincida con el cielo y lo conllevan con una Función de Propagación de Punto (PSF), RotCurves proyecta el haz sobre el plano medio galáctico. Para un ángulo de inclinación $i$ dado, una PSF gaussiana circular se transforma en una gaussiana bivariante con una relación de ejes de $1/\cos(i)$. La velocidad observada y la dispersión de velocidad en un radio dado se calculan como promedios ponderados sobre este "haz efectivo", teniendo en cuenta el perfil de emisión (ponderación por luz) y la proyección en la línea de visión.
• Procedimiento de Ajuste: El paquete incluye un motor de inferencia bayesiana integrado utilizando el muestreador MCMC emcee. Ajusta el modelo a las curvas de rotación observadas minimizando una función de log-probabilidad que comprende una verosimilitud de mínimos cuadrados y priors definidos por el usuario (uniformes o gaussianos) sobre parámetros tales como la masa bariónica, el radio efectivo, la relación bulbo-a-total, la dispersión de velocidad intrínseca y la masa del halo de materia oscura.

Contribuciones Clave

1. Eficiencia Computacional: RotCurves logra un tiempo de ejecución de $\sim 10-20$ ms por realización para una rejilla típica, lo que representa una velocidad aproximadamente 250 veces mayor que dysmalpy (que tarda $\sim 3-4$ segundos para la misma tarea). Esta aceleración se logra evitando la construcción de hipercubos 4D completos.
2. Proyección de Haz en 2D: El artículo introduce y valida un método de proyección 2D para el ensanchamiento del haz que tiene en cuenta la inclinación sin requerir la integración de volumen 3D completa.
3. Evaluación Comparativa (Benchmarking): Los autores comparan rigurosamente RotCurves contra dysmalpy utilizando observaciones sintéticas de un modelo de galaxia fiduciaria en $z=2$.

Resultados

• Precisión vs. dysmalpy: Para sistemas bien resueltos (FWHM de la PSF $< R_{\text{eff}}$), RotCurves reproduce las curvas de rotación y dispersión observadas simuladas por dysmalpy con un error inferior al 5% hasta $3R_{\text{eff}}$. En sistemas marginalmente resueltos (FWHM de la PSF $\gtrsim 1.5 R_{\text{eff}}$), las discrepancias aumentan hasta un 15%.
  • Diferencias Sistemáticas: RotCurves tiende a subestimar ligeramente la dispersión de velocidad en las regiones internas (por $\sim 10\%$) porque solo considera los movimientos del plano medio, mientras que dysmalpy integra partículas fuera del plano medio que ensanchan la línea.
• Recuperación de Parámetros: Utilizando el ajuste MCMC en datos sintéticos, RotCurves recupera con éxito los parámetros intrínsecos del modelo (masa, radio, dispersión, fracción de materia oscura) a través de una amplia gama de propiedades galácticas.
  • Sesgos: Surgen sesgos sistemáticos para el radio efectivo ($R_{\text{eff}}$), que suele sobreestimarse por un factor de $\sim 1.5$, y para masas de disco bajas ($M_{\text{disk}} \lesssim 3 \times 10^9 M_\odot$), donde la masa del disco está mal restringida.
  • Relación Señal-Ruido (S/N): La recuperación de parámetros es excelente en S/N altos ($\gtrsim 25$). En S/N más bajos, las desviaciones aumentan, particularmente para la masa del disco y las fracciones de materia oscura. Sin embargo, la masa del bulbo y la dispersión de velocidad intrínseca se mantienen relativamente robustas incluso con S/N más bajos.
• Límites de Inclinación: La herramienta es más adecuada para inclinaciones entre $10^\circ$ y $80^\circ$.

Significancia y Reivindicaciones
El artículo posiciona a RotCurves principalmente como una herramienta exploratoria para el testeo rápido de supuestos de modelos de masa y estudios de parámetros, más que como un reemplazo del modelado 3D completo en todos los contextos. Su significancia radica en su capacidad para:

• Procesar eficientemente grandes muestras de datos observacionales de los actuales y futuros sondeos de IFU de gran escala (por ejemplo, VLT/MAVIS, JWST/NIRSpec, ELT/MOSAIC).
• Facilitar la prueba de diferentes priors, formas de perfil y efectos sistemáticos sin la carga computacional de las simulaciones 3D completas.
• Proporcionar estimaciones de parámetros rápidas y razonables para galaxias de alto desplazamiento al rojo, sirviendo como punto de partida para un análisis más detallado.

Los autores enfatizan que, si bien RotCurves es altamente eficiente y precisa para sistemas bien resueltos, los usuarios deben tener precaución al interpretar los resultados de sistemas mal resueltos o con bajo S/N, donde los efectos de ensanchamiento del haz son más severos y los sesgos sistemáticos (particularmente en $R_{\text{eff}}$) son más pronunciados. El código está disponible públicamente en GitHub.