Radial velocity and atmospheric parameter calculations for the GaiaNIR spectrograph

Resumen Técnico: Cálculos de velocidad radial y parámetros atmosféricos para el espectrógrafo GaiaNIR

Planteamiento del problema
La misión Gaia ha revolucionado la astrometría y la espectroscopia en el régimen óptico, pero está severamente limitada por la extinción debida al polvo interestelar, lo que impide un mapeo exhaustivo del interior de la Galaxia, el bulbo y el lado lejano del disco de la Vía Láctea. La propuesta de la misión GaiaNIR tiene como objetivo superar esto operando en el infrarrojo cercano (NIR, 800–2300 nm), donde la extinción es significativamente menor. Aunque GaiaNIR planea realizar astrometría y fotometría, su retorno científico podría potenciarse exponencialmente añadiendo un espectrógrafo de resolución moderadamente alta ($R > 10,000$) para medir velocidades radiales ($v_{rad}$), parámetros atmosféricos ($T_{eff}$, $\log g$, $[M/H]$) y abundancias químicas para cientos de millones de estrellas. Sin embargo, las restricciones técnicas (tamaño del detector, telemetría y requisitos de enfriamiento) impiden que un único espectrógrafo cubra todo el rango de 800–2300 nm con alta resolución. En consecuencia, el desafío crítico consiste en identificar la ventana de longitud de onda óptima que equilibre una alta precisión en la velocidad radial con la capacidad de derivar parámetros atmosféricos y abundancias precisas, maximizando al mismo tiempo el beneficio de la reducción de la extinción interestelar.

Metodología
Los autores llevaron a cabo un estudio de simulación sistemático para identificar ventanas espectrales óptimas sin depender de una biblioteca empírica exhaustiva de espectros de alta resolución en el infrarrojo cercano, la cual actualmente no existe.

• Generación de datos sintéticos: El estudio utilizó la biblioteca de espectros sintéticos BOSZ, generando 10,000 espectros ficticios a través de un amplio rango de parámetros estelares ($T_{eff}$: 2800–10,000 K, $\log g$: 0–5.5 dex, $[M/H]$: -1.5 a 0.5 dex). Los modelos incorporaron aproximaciones LTE utilizando atmósferas MARCS y ATLAS9, opacidades continuas actualizadas y 23 listas de líneas moleculares.
• Parámetros de simulación: Los espectros se generaron a cuatro resoluciones ($R = 5,000, 10,000, 15,000, 20,000$) y se sometieron a ruido gaussiano aleatorio para simular relaciones señal-ruido (S/N) entre 5 y 100.
• Análisis de velocidad radial: Los espectros ficticios fueron desplazados Doppler y correlacionados cruzadamente con plantillas ideales utilizando la función crosscorrRV. El análisis escaneó ventanas de 100 nm en el rango de 850–2250 nm en incrementos de 5 nm. El rendimiento se evaluó utilizando la dispersión ($\sigma$), la desviación absoluta de la mediana (MAD) y el percentil 95 (P95) de los residuales de velocidad.
• Derivación de parámetros atmosféricos: Se empleó el código FERRE para ajustar $T_{eff}$, $\log g$ y $[M/H]$ a los espectros ficticios mediante la minimización de $\chi^2$ contra rejillas precomputadas.
• Estratificación de la muestra: El análisis se estratificó en cinco submuestras basadas en temperatura y gravedad superficial (Enanas Rojas, Gigantes Rojas, enanas/gigantes de tipo FGK y estrellas calientes de tipo A) para dar cuenta de la dominancia variable de las líneas atómicas frente a las moleculares.
• Validación: La metodología fue validada comparando las precisiones derivadas contra el rendimiento empírico de Gaia Data Release 3 (DR3) en la región óptica del triplete de Ca II, confirmando que el enfoque de simulación produce estimaciones de incertidumbre realistas para objetivos de alto S/N.

Resultos clave
El estudio identificó seis ventanas candidatas para un refinamiento posterior, estrechando finalmente el enfoque hacia dos regiones primarias:

1. La banda K (1926–2000 nm):

   • Rendimiento: Esta región ofrece la mayor precisión de velocidad radial para estrellas frías (tipos FGKM), con precisiones estimadas de 160–260 m s$^{-1}$ para $R \approx 16,000–20,000$ en los objetivos más brillantes.
   • Parámetros atmosféricos: Proporciona una precisión de $T_{eff}$ mejor que 100 K y una precisión de metalicidad de alrededor de 0.06–0.08 dex, comparable a las capacidades actuales de Gaia. Sin embargo, la precisión de $\log g$ es algo inferior a los resultados ópticos.
   • Extinción: La banda K ofrece una ventaja crítica con un enrojecimiento interestelar aproximadamente 7.4 veces menor que en el óptico, permitiendo el mapeo de regiones oscurecidas por el polvo.
   • Abundancias: Se pueden derivar diez especies (O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Ni), aunque las enanas M están limitadas principalmente a Ca debido al empaquetamiento de H$_2$O.
   • Limitaciones: El rendimiento se degrada significativamente para estrellas más calientes que 7000–8000 K debido a la falta de líneas metálicas, aunque la línea de hidrógeno Br-$\delta$ ofrece una utilidad limitada.

2. La transición J/H (1158–1202 nm):

   • Rendimiento: Esta región proporciona una precisión de velocidad radial ligeramente menor ($\approx 340$ m s$^{-1}$) comparada con la banda K, pero ofrece un rango de temperatura efectiva más amplio, extendiéndose hasta 12,000–16,000 K.
   • Parámetros atmosféricos: Las precisiones de $T_{eff}$ y metalicidad son competitivas, aunque la determinación de $\log g$ sigue siendo un desafío.
   • Extinción: La reducción de la extinción es marginal en comparación con el óptico, lo que limita su valor estratégico para sondear las regiones galácticas más oscurecidas.
   • Abundancias: Diez especies están disponibles (Mg, Si, K, Ca, Ti, V, Cr, Mn, Co, Ni). Una ventaja distintiva es la capacidad de medir abundancias para enanas M por debajo de 4000 K (Mg, Si, K, Ca, Ti, Cr) donde las líneas de H$_2$O son menos dominantes que en la banda K.

Conclusiones y significación
El artículo concluye que la ventana 1926–1968 nm (centrada en $\lambda_C \approx 1947$ nm) representa el compromiso estratégico óptimo para el espectrógrafo GaiaNIR. Mientras que la ventana más estrecha de 1984–2000 nm ofrece la mayor precisión de velocidad radial, sacrifica la exactitud de los parámetros atmosféricos. La región 1926–1968 nm, operando a resoluciones de hasta $R \approx 20,000$, entrega precisiones de velocidad radial de 160–260 m s$^{-1}$ y precisiones de parámetros atmosféricos comparables a las de Gaia, todo ello aprovechando la transmisión superior de la banda K a través del polvo.

Los autores enfatizan que esta selección permite a GaiaNIR cumplir con su mandato principal: mapear el interior de la Galaxia oscurecido por el polvo y el lado lejano del disco. Sin embargo, el artículo señala advertencias significativas:

• Aglomeración (Crowding): La banda K puede sufrir una severa aglomeración espectral en el plano galáctico, lo que potencialmente requiere resoluciones más bajas o tamaños de ventana menores, lo cual requiere más estudios.
• Enanas M: La incapacidad de derivar abundancias robustas para las enanas M más frías (debido al H$_2$O) es una desventaja señalada, ya que estas estrellas constituyen una gran fracción de la población objetivo.
• Trabajo futuro: Los autores afirman que la creación de una nueva lista de líneas de alta calidad para la banda K es esencial, ya que las observaciones actuales desde la Tierra están fuertemente contaminadas por líneas telúricas. Los estudios futuros deben abordar los compromisos entre resolución, tamaño de ventana y aglomeración estelar para finalizar el diseño del espectrógrafo.

El estudio no pretende haber resuelto todos los desafíos de diseño, pero proporciona un marco basado en datos para seleccionar la región de longitud de onda que maximice el retorno científico de la misión GaiaNIR dentro de sus restricciones técnicas.