The Shape of the Vertical Action Distribution Locates the Scatterers that Heat the Galactic Disc

Resumen Técnico: La forma de la distribución de la acción vertical localiza a los dispersores que calientan el disco galáctico

Planteamiento del problema
La dispersión de velocidades verticales de las estrellas en el disco de la Vía Láctea aumenta con la edad, un fenómeno conocido como calentamiento del disco. Si bien la existencia de este calentamiento está bien establecida, la identidad y la distribución espacial de los "dispersores" responsables sigue sin resolverse. Existen dos clases principales de modelos:

1. Dispersores de volumen completo (volume-filling): Perturbaciones tales como la estructura espiral, la barra Galáctica o subhalos de materia oscura que actúan sobre una estrella durante todo su órbita.
2. Dispersores confinados al plano medio: Perturbaciones, como las Nubes Moleculares Gigantes (GMCs, por sus siglas en inglés), que están confinadas a una capa delgada alrededor del plano medio Galáctico e interactúan con las estrellas solo durante sus cruces verticales.

Los intentos previos para distinguir entre estos mecanismos utilizando la relación edad–dispersión de velocidad o la forma del elipsoide de velocidades han sido inconclusos debido a las degeneraciones entre la abundancia, la masa de los dispersores y la historia del calentamiento. El artículo argumenta que la acción vertical ($J_z$), un invariante adiabático, ofrece un diagnóstico superior porque etiqueta la órbita en lugar de la velocidad instantánea y es insensible al crecimiento del potencial del disco a lo largo del tiempo.

Metodología
Los autores derivan la distribución teórica de la acción vertical para una población coeval de estrellas calentadas por dispersores con una distribución de altura vertical arbitraria.

• Derivación Teórica:

  • El proceso de calentamiento se modela como una difusión en el espacio de la acción vertical utilizando la ecuación de Fokker–Planck.
  • Los autores derivan la difusividad $D(J_z)$ como una función del perfil de altura del dispersor $q(z)$. Demuestran que la difusividad escala como una ley de potencia $D \propto J_z^b$, donde el exponente $b$ está determinado únicamente por la geometría de los dispersores.
  • Caso de volumen completo: Si los dispersores llenan el volumen, la difusividad crece linealmente con la acción ($b=1$), resultando en una distribución exponencial $p(J_z) \propto \exp(-J_z/J_0)$.
  • Caso de capa en el plano medio: Si los dispersores están confinados a una capa delgada, el tiempo de interacción disminuye a medida que la velocidad de cruce de la estrella aumenta. Esto reduce el crecimiento de la difusividad a $D \propto J_z^{1/2}$ ($b=1/2$), resultando en una distribución de "exponencial estirada" $p(J_z) \propto \exp[-(J_z/J_0)^{3/2}]$.
  • La solución general para una población coeval es $p(J_z) \propto \exp[-(J_z/J_0)^{2-b}]$, donde la forma (el exponente de la cola) depende solo de la geometría ($b$), mientras que la escala ($J_0$) codifica la amplitud y la historia del calentamiento.

• Datos y Ajuste:

  • El modelo se ajusta a una muestra de 7,589 estrellas de la rama de la agrupación roja (red clump) de bajo $\alpha$ del sondeo APOGEE DR14, cruzada con la astrometría de Gaia DR2.
  • La muestra se restringe a estrellas con radios efectivos $5 \le R_{\text{eff}} \le 10$ kpc y edades $2 \le \tau \le 8$ Gyr.
  • Un componente crítico de la metodología es el tratamiento riguroso de la función de selección, que omite preferencialmente las estrellas de alta acción cerca del límite de magnitud del sondeo. Los autores aplican una función de selección por estrella derivada en Ting & Rix (2019) directamente a la verosimilitud del modelo.
  • El procedimiento de ajuste marginaliza los errores de medición en edad, acción y radio mediante cuadratura gaussiana. Crucialmente, la escala del parámetro $J_0$ se deja libre en cada edad (parametrizada por cuatro nudos y una pendiente radial). Esto permite la medición del exponente de forma $b$ sin asumir una historia de calentamiento específica o una dispersión de nacimiento, desacoplando efectivamente la geometría de la tasa de calentamiento.

Resultados Clave

1. Exponente Medido: El ajuste arroja un exponente de difusividad de $b = 0.51^{+0.06}_{-0.7}$.

   • Este valor es consistente con la predicción para una capa confinada al plano medio ($b=0.5$) pero es fuertemente inconsistente con la predicción de volumen completo ($b=1$).
   • Pruebas sintéticas confirman que el método puede recuperar el exponente correcto sin sesgos y separa las dos geometrías con una significancia de $5.6\sigma$.
   • El resultado se mantiene a través de diferentes contenedores radiales y cortes de edad, y es robusto frente a variaciones en los errores de acción asumidos.

2. Masa Efectiva del Dispersor:

   • Utilizando el exponente $b \approx 0.5$ y las escalas libres medidas, los autores calculan la amplitud del calentamiento.
   • Para el contenedor de edad de 2–4 Gyr, al comparar la amplitud del calentamiento con la densidad superficial molecular actual, se obtiene una masa efectiva del dispersor de $2.7 \times 10^6 M_\odot$.
   • Este valor cae dentro del rango ($1.4\text{--}3.4 \times 10^6 M_\odot$) derivado de catálogos de nubes independientes y funciones de masa (por ejemplo, Miville-Deschênes et al. 2017).
   • Para las estrellas más viejas, la masa efectiva aparente aumenta si se compara con la densidad de gas actual. Sin embargo, cuando se corrige por la historia decreciente de formación estelar (y, por ende, la densidad de gas decreciente) durante los últimos varios gigaños, la masa efectiva para todos los contenedores de edad converge al rango de $1.9\text{--}2.8 \times 10^6 M_\odot$.

3. Historia del Calentamiento:

   • Las escalas libres crecen con la edad como $J_0 \propto \tau^{1.12}$. Esto es más rápido que el $\tau^{0.67}$ esperado para un calentamiento constante con un nacimiento frío, lo que indica que el disco se calentó más fuertemente en el pasado.

Significado y Reivindicaciones
El artículo afirma resolver el debate de larga data sobre la ubicación de los dispersores del disco Galáctico mediante el uso de la forma de la distribución de la acción vertical en lugar de solo su media o dispersión de velocidad.

• Localización de los Dispersores: El exponente medido $b \approx 0.5$ proporciona evidencia directa de que el mecanismo principal para el calentamiento vertical está localizado cerca del plano Galáctico. Esto descarta los mecanismos de volumen completo (como subhalos de materia oscura, ondas de flexión o perturbaciones espirales directas) como la fuente dominante de calentamiento difusivo. Los autores estiman que cualquier componente de volumen completo contribuye con menos del 7% al calentamiento.
• Identificación de los Dispersores: La combinación de la geometría del plano medio y la masa efectiva derivada ($\sim 10^6 M_\odot$) identifica a las Nubes Moleculares Gigantes (GMCs) como los dispersores. La masa es consistente con la parte superior de la función de masa de las GMC y la geometría coincide con la capa delgada de gas molecular.
• Resolución de Discrepancias Previas: El artículo explica por qué estudios previos (por ejemplo, Binney & Lacey 1988) infirieron $b=1$ (volumen completo). Esos estudios analizaron poblaciones de edades mixtas o dispersiones de velocidad sin resolver la forma a una edad fija. Los autores demuestran que mezclar poblaciones con diferentes escalas de calentamiento aplana artificialmente la cola de la distribución, simulando una geometría de volumen completo ($b \approx 1$) incluso si la verdadera geometría subyacente es una capa delgada ($b=0.5$).
• Mecanismo vs. Fuente de Energía: Los autores aclaran que, si bien el dispersión ocurre en el plano medio (vía nubes), la energía puede originarse en movimientos dentro del plano (por ejemplo, impulsados por la estructura espiral) que se redirigen verticalmente durante los cruces de las nubes. El exponente mide la localización de la redirección, no la fuente última de la energía.

En conclusión, el artículo establece que el disco de la Vía Láctea es calentado cerca del plano por una población evolutiva de objetos con la masa de las Nubes Moleculares Gigantes, y que la forma de la distribución de la acción vertical es un diagnóstico robusto e independiente del modelo para la distribución espacial de los agentes de dispersión.