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

Resumo Técnico: A Forma da Distribuição de Ação Vertical Localiza os Espalhadores que Aquecem o Disco Galáctico

Enunciado do Problema
A dispersão de velocidade vertical das estrelas no disco da Via Láctea aumenta com a idade, um fenômeno conhecido como aquecimento do disco. Embora a existência desse aquecimento seja bem estabeleada, a identidade e a distribuição espacial dos "espalhadores" responsáveis permanecem sem resolução. Existem duas classes primárias de modelos:

1. Espalhadores de preenchimento de volume: Perturbações como a estrutura espiral, a Barra Galáctica ou subhalos de matéria escura que atuam sobre uma estrela durante toda a sua órbita.
2. Espalhadores confinados ao plano médio: Perturbações, como Nuvens Moleculares Gigantes (GMCs), que são confinadas a uma camada fina ao redor do plano médio galáctico e interagem com as estrelas apenas durante suas passagens verticais.

Tentativas anteriores de distinguir entre esses mecanismos usando a relação dispersão de velocidade–idade ou a forma do elipsoide de velocidade foram inconclusivas devido a degenerescências entre a abundância, a massa do espalhador e a história do aquecimento. O artigo argumenta que a ação vertical ($J_z$), uma invariante adiabática, oferece um diagnóstico superior porque rotula a órbita em vez da velocidade instantânea e é insensível ao crescimento do potencial do disco ao longo do tempo.

Metodologia
Os autores derivam a distribuição teórica da ação vertical para uma população coeval de estrelas aquecida por espalhadores com uma distribuição de altura vertical arbitrária.

• Derivação Teórica:

  • O processo de aquecimento é modelado como uma difusão no espaço da ação vertical usando a equação de Fokker–Planck.
  • Os autores derivam a difusividade $D(J_z)$ como uma função do perfil de altura $q(z)$ do espalhador. Eles mostram que a difusividade escala como uma lei de potência $D \propto J_z^b$, onde o expoente $b$ é determinado unicamente pela geometria dos espalhadores.
  • Caso de preenchimento de volume: Se os espalhadores preencherem o volume, a difusividade cresce linearmente com a ação ($b=1$), resultando em uma distribuição exponencial $p(J_z) \propto \exp(-J_z/J_0)$.
  • Caso de camada no plano médio: Se os espalhadores estiverem confinados a uma camada fina, o tempo de interação diminui à medida que a velocidade de travessia da estrela aumenta. Isso reduz o crescimento da difusividade para $D \propto J_z^{1/2}$ ($b=1/2$), resultando em uma distribuição de "exponencial estendida" $p(J_z) \propto \exp[-(J_z/J_0)^{3/2}]$.
  • A solução geral para uma população coeval é $p(J_z) \propto \exp[-(J_z/J_0)^{2-b}]$, onde a forma (expoente da cauda) depende apenas da geometria ($b$), enquanto a escala ($J_0$) codifica a amplitude e a história do aquecimento.

• Dados e Ajuste:

  • O modelo é ajustado a uma amostra de 7.589 estrelas de clump vermelho de baixo $[\alpha/\text{Fe}]$ do levantamento APOGEE DR14, cruzados com astrometria do Gaia DR2.
  • A amostra é restrita a estrelas com raios efetivos $5 \le R_{\text{eff}} \le 10$ kpc e idades $2 \le \tau \le 8$ Gyr.
  • Um componente crítico da metodologia é o tratamento rigoroso da função de seleção, que perde preferencialmente estrelas de alta ação próximas ao limite de magnitude do levantamento. Os autores aplicam uma função de seleção por estrela derivada de Ting & Rix (2019) diretamente à verossimilhança do modelo.
  • O procedimento de ajuste marginaliza os erros de medição em idade, ação e raio usando quadratura gaussiana. Crucialmente, o parâmetro de escala $J_0$ é deixado livre em cada idade (parametrizado por quatro nós e um declive radial). Isso permite a medição do expoente de forma $b$ sem assumir uma história de aquecimento específica ou dispersão de nascimento, efetivamente desacoplando a geometria da taxa de aquecimento.

Principais Resultados

1. Expoente Medido: O ajuste produz um declive de difusividade de $b = 0,51^{+0,06}_{-0,07}$.

   • Este valor é consistente com a previsão para uma camada confinada ao plano médio ($b=0,5$), mas fortemente inconsistente com a previsão de preenchimento de volume ($b=1$).
   • Testes sintéticos confirmam que o método pode recuperar o expoente correto sem viés e separa os dois geometrias com uma significância de $5,6\sigma$.
   • O resultado mantém-se em diferentes bins radiais e cortes de idade, e é robusto contra variações nos erros de ação assumidos.

2. Massa Efetiva do Espalhador:

   • Usando o declive de difusividade medido e os valores de escala livre derivados ($b \approx 0,5$), os autores calculam a amplitude do aquecimento.
   • Para o bin de idade de 2–4 Gyr, comparando a amplitude de aquecimento com a densidade de superfície molecular atual, obtém-se uma massa efetiva do espalhador de $2,7 \times 10^6 M_\odot$.
   • Este valor cai dentro da faixa ($1,4\text{--}3,4 \times 10^6 M_\odot$) derivada de catálogos de nuvens independentes e funções de massa (por exemplo, Miville-Deschênes et al. 2017).
   • Para estrelas mais velhas, a massa efetiva aparente aumenta se comparada à densidade de gás atual. No entanto, quando corrigida pela declinante história de formação estelar (e, portanto, pela declinante densidade de gás) ao longo dos últimos vários giganos, a massa efetiva para todos os bins de idade converge para a faixa de $1,9\text{--}2,8 \times 10^6 M_\odot$.

3. História de Aquecimento:

   • As escalas livres crescem com a idade como $J_0 \propto \tau^{1,12}$. Isso é mais rápido do que o esperado ($\tau^{0,67}$) para uma taxa de aquecimento constante com um nascimento frio, indicando que o disco aqueceu mais fortemente no passado.

Significância e Alegações
O artigo afirma resolver o debate de longa data sobre a localização dos espalhadores do disco Galáctico ao utilizar a forma da distribuição da ação vertical em vez de apenas sua média ou dispersão de velocidade.

• Localização dos Espalhadores: O expoente medido $b \approx 0,5$ fornece evidência direta de que o principal mecanismo de aquecimento vertical é localizado no plano médio Galáctico. Isso descarta mecanismos de preenchimento de volume (como subhalos de matéria escura, ondas de flexão ou perturbações espirais diretas) como a fonte dominante de aquecimento difusivo. Os autores estimam que qualquer componente de preenchimento de volume contribui com menos de 7% para o aquecimento.
• Identificação dos Espalhadores: A combinação da geometria do plano médio e da massa efetiva derivada ($\sim 10^6 M_\odot$) identifica as Nuvens Moleculares Gigantes (GMCs) como os espalhadores. A massa é consistente com o limite superior da função de massa das GMCs, e a geometria coincide com a camada fina de gás molecular.
• Resolução de Discrepâncias Anteriores: O artigo explica por que estudos anteriores (ex: Binney & Lacey 1988) inferiram $b=1$ (preenchimento de volume). Esses estudos analisaram populações de idades mistas ou dispersões de velocidade sem resolver a forma em idades fixas. Os autores demonstram que misturar populações com diferentes escalas de aquecimento artificialmente achata a cauda da distribuição, mimetizando uma geometria de preenchimento de volume ($b \approx 1$), mesmo que a verdadeira geometria subjacente seja uma camada fina ($b=0,5$).
• Mecanismo vs. Fonte de Energia: Os autores esclarecem que, embora o espalhamento ocorra no plano médio (via nuvens), a energia pode originar-se de movimentos no plano (ex: impulsionados por estrutura espiral) que são redirecionados verticalmente durante as passagens pelas nuvens. O expoente mede a localização do redirecionamento, não a fonte última da energia.

Em conclusão, o artigo estabelece que o disco da Via Láctea é aquecido próximo ao plano por uma população evolutiva de objetos com a massa de Nuvens Moleculares Gigantes, e que a forma da distribuição da ação vertical é um diagnóstico robusto e independente de modelo para a distribuição espacial dos agentes de espalhamento.