Recalibration of the H$\alpha$ surface brightness-radius relation for planetary nebulae using Gaia DR3: new distances and the Milky Way oxygen radial gradient

Resumen Técnico: Recalibración de la Relación Brillo Superficial Hα–Radio para Nebulosas Planetarias utilizando Gaia DR3

Enunciado del Problema
La distribución espacial de los elementos químicos en el disco galáctico sirve como una restricción crítica para los modelos de evolución de galaxias. Sin embargo, los estudios que utilizan nebulosas planetarias (NPs) como trazadores han estado históricamente limitados por grandes incertidumbres en la determinación de distancias. Si bien las escalas de distancia estadísticas, como la relación brillo superficial Hα–radio establecida por Frew et al. (2016, aquíafter F16), han mejorado las estimaciones de distancia, siguen sujetas a incertidumbres que a menudo superan factores de dos. La llegada de la misión Gaia, particularmente el tercer lanzamiento de datos (DR3), ofrece paralajes astrométricos precisos para cientos de estrellas centrales de nebulosas planetarias (ECNPs), proporcionando una oportunidad para medir distancias directamente y recalibrar las escalas estadísticas para reducir estas incertidumbres.

Metodología
Los autores presentan una recalibración exhaustiva de la relación brillo superficial Hα–radio utilizando paralajes de Gaia DR3. La metodología implica varios pasos clave:

1. Identificación y Selección de Fuentes: Se compiló un catálogo de 1.200 NPs con coordenadas fiables y diámetros angulares. Los autores consultaron el archivo de Gaia para fuentes dentro de un radio de búsqueda definido por la geometría de la nebulosa. Para garantizar la fiabilidad astrométrica, las fuentes se filtraron basándose en el Error de Peso Unitario Renormalizado (RUWE < 1.4) y el número de periodos de visibilidad (> 8).
2. Emparejamiento de ECNP: Se empleó un método color-distancia, mejorado respecto a trabajos anteriores, para identificar las ECNPs más probables. Esto implicó aplicar restricciones de temperatura (Teff > 13.000 K) y relaciones color-temperatura. Las identificaciones resultantes se clasificaron en grupos de fiabilidad (A, B, C) y se cruzaron con el catálogo de Chornay & Walton (2021). Se excluyeron fuentes con valores de extinción visual que se desviaban significativamente de los mapas de polvo galáctico, resultando en una muestra final de 436 NPs con asociaciones de alta probabilidad.
3. Derivación de Distancias:
   • Recalibración: Utilizando 135 NPs con incertidumbres fraccionales de paralaje ($f = \sigma_\varpi/\varpi$) $\le 0.20$, los autores realizaron una regresión de Mínimos Cuadrados Ordinarios (OLS) con simulaciones de Monte Carlo (MC) para recalibrar la relación $\log R_{pc}$ vs. $\log S_{H\alpha}$.
   • Distancias Bayesianas: Para una submuestra de 415 NPs con paralajes de Gaia, se derivaron distancias bayesianas. Este enfoque utiliza la distancia estadística recalibrada como priori y el paralaje de Gaia como verosimilitud. Para las 711 NPs restantes sin contrapartes de Gaia, se adoptó la distribución posterior como distribución priori.
4. Análisis de Abundancia Química: Las abundancias de oxígeno (O/H) se compilaron de Stanghellini & Haywood (2018) y Bucciarelli & Stanghellini (2023, aquíafter BS23), recalcularon utilizando un factor de corrección de ionización (ICF) consistente. La muestra se separó en poblaciones de disco delgado, disco grueso, bulbo y halo basándose en velocidades peculiares ($V_{pec}$) derivadas de los movimientos propios y velocidades radiales de Gaia.
5. Ajuste de Gradientes: El gradiente radial O/H para 231 NPs de disco se analizó utilizando regresión lineal por tramos. Se probaron tres modelos: un ajuste lineal único, un ajuste segmentado con una ruptura y un ajuste segmentado con dos rupturas. La selección del modelo se guió mediante el Criterio de Información de Akaike (AIC). Además, se generó un mapa de abundancia 2D utilizando interpolación de kriging universal.

Contribuciones y Resultados Clave

• Escala de Distancia Recalibrada: La nueva calibración produce una relación más plana en comparación con F16, resultando en distancias estadísticas más cortas. La relación adoptada es $\log R_{pc} = (-0.257 \pm 0.008) \times \log S_{H\alpha} - (1.51 \pm 0.03)$. Esta escala demuestra métricas de bondad de ajuste mejoradas ($\langle k \rangle = 1.020$, $\langle \sigma \rangle = 0.142$) en comparación con la calibración original de F16 aplicada a los mismos datos.
• Comparaciones de Distancia: Las distancias bayesianas derivadas son, en promedio, un 16% más pequeñas que las de Chornay & Walton (2021) y un 28% más pequeñas que las distancias de F16. El método bayesiano reduce significativamente la dispersión y el sesgo para NPs con incertidumbres fraccionales de paralaje $f < 0.5$.
• Gradiente Radial de Oxígeno: El análisis de 231 NPs de disco revela que un gradiente lineal único no es el mejor ajuste. El modelo con una ruptura es estadísticamente preferido (AIC más bajo), indicando un cambio en la pendiente cerca del radio solar ($R \approx 7.25 \pm 1.00$ kpc).
  • Disco Interior ($R \lesssim 8$ kpc): El gradiente es plano o ligeramente positivo ($0.026 \pm 0.028$ dex/kpc).
  • Disco Exterior ($R \gtrsim 8$ kpc): El gradiente se vuelve más pronunciado y negativo ($-0.034 \pm 0.009$ dex/kpc).
  • También se probó un modelo de dos rupturas, pero arrojó resultados inconclusos sobre la segunda ruptura, aunque sugirió una estructura tipo escalón cerca del Sol.
• Disco Delgado vs. Disco Grueso: Al separar la muestra por cinemática, las NPs del disco delgado exhiben un gradiente más pronunciado ($-0.035 \pm 0.012$ dex/kpc) en comparación con el disco grueso ($-0.005 \pm 0.017$ dex/kpc). Los autores sugieren que la ruptura aparente en la muestra combinada puede resultar de la superposición de estas poblaciones estelares distintas.
• Asimetría Azimutal: El mapa de kriging 2D revela una asimetría azimutal modesta, con abundancias aumentadas cerca de la barra galáctica en longitudes positivas. El círculo solar parece actuar como una región de transición, creando una estructura de abundancia bimodal entre las regiones interiores y exteriores.

Significado y Afirmaciones
El artículo afirma que la recalibración de la relación brillo superficial Hα–radio utilizando paralajes de Gaia DR3 proporciona un indicador de distancia más robusto para las NPs galácticas, abordando incertidumbres de larga data que han obstaculizado el estudio de la evolución química.

El significado principal radica en la evidencia de un gradiente radial de oxígeno no lineal en el disco de la Vía Láctea. La detección de una ruptura cerca del radio solar desafía la suposición de un gradiente lineal constante. Los autores proponen varios mecanismos físicos para esta característica:

1. Efectos Dinámicos: La ruptura puede reflejar cambios en la eficiencia de formación estelar impulsados por la barra galáctica o la resonancia de corrotación de los brazos espirales.
2. Poblaciones Estelares: La ruptura observada puede ser un artefacto de la superposición de los discos delgado y grueso, que poseen diferentes gradientes de metalicidad radial.
3. Evolución Química: Los resultados proporcionan nuevas restricciones para los modelos de evolución química galáctica, particularmente con respecto al impacto de las estructuras no axisimétricas (barras y espirales) en los flujos de gas y las distribuciones de metalicidad.

Los autores concluyen que, aunque los datos apoyan un cambio en la pendiente cerca del radio solar, se requiere una investigación adicional con muestras más grandes y homogéneas y con incertidumbres reducidas para confirmar definitivamente el origen de estas rupturas y refinar los modelos de la evolución química de la Vía Láctea.