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

技術要約：垂直作用量分布の形状が銀河円盤を加熱する散乱体の所在を特定する

問題提起
天の川銀河の円盤における星の垂直速度分散は、年齢とともに増加する。これは「円盤加熱（disc heating）」として知られる現象である。この加熱の存在自体は確立されているが、その原因となる「散乱体（scatterers）」の正体とその空間分布は未解決のままである。主に2つのモデルクラスが存在する：

1. 体積充填型散乱体（Volume-filling scatterers）： 螺旋構造、銀河棒、あるいはダークマター・サブハロのように、星の軌道全体にわたって作用する摂動。
2. 中間面限定型散乱体（Midplane-confined scatterers）： 巨大分子雲（GMC）のように、銀河中間面の薄い層に限定されており、星が垂直方向を横切る際のみ相互作用する摂動。

年齢–速度分散関係や速度楕円体の形状を用いてこれらのメカニズムを区別しようとするこれまでの試みは、散乱体の存在量、質量、および加熱履歴の間の縮退により、決定的な結論を得るには至らなかった。本論文は、**垂直作用量（vertical action, $J_z$）**が、瞬時的な速度ではなく軌道をラベル付けし、時間の経過に伴う円盤ポテンシャルの成長にも影響を受けないため、優れた診断指標であることを主張している。

手法
著者らは、任意の垂直高さ分布を持つ散乱体によって加熱された共時集団（coeval population）の垂直作用量の理論的分布を導出した。

• 理論的導出:

  • 加熱プロセスは、垂直作用量空間における拡散としてフォッカー・プランク方程式を用いてモデル化される。
  • 著者らは、散乱体の高さプロファイル $q(z)$ の関数として、拡散係数 $D(J_z)$ を導出した。拡散係数はべき乗則 $D \propto J_z^b$ に従い、指数 $b$ は散乱体の幾何学的形状のみによって決定されることを示している。
  • 体積充填型の場合： 散乱体が体積を満たしている場合、拡散係数は作用に対して線形に増大し（$b=1$）、その結果、指数分布 $p(J_z) \propto \exp(-J_z/J_0)$ が得られる。
  • 中間面層型の場合： 散乱体が薄い層に限定されている場合、星の横切り速度が増すにつれて相互作用時間が減少する。これにより、拡散の増大は $D \propto J_z^{1/2}$ （$b=1/2$）へと抑制され、「引き伸ばされた指数分布（stretched exponential）」$p(J_z) \propto \exp[-(J_z/J_0)^{3/2}]$ が得られる。
  • 共時集団に対する一般解は $p(J_z) \propto \exp[-(J_z/J_0)^{2-b}]$ であり、その形状（テールの指数）は幾何学（$b$）にのみ依存し、スケール（$J_0$）は加熱の振幅と履歴をコード化する。

• データとフィッティング:

  • モデルは、APOGEE DR14サーベイから得られた、Gaia DR2のアストロメトリとクロスマッチされた7,589個の低$\alpha$レッドクランプ星のサンプルに適合させている。
  • サンプルは、有効半径 $5 \le R_{\text{eff}} \le 10$ kpc、年齢 $2 \le \tau \le 8$ Gyrの星に限定されている。
  • 手法の重要な要素は、サーベイの等級限界付近で高作用量の星を優先的に取りこぼしてしまう**選択関数（selection function）**の厳密な取り扱いである。著者らは、Ting & Rix (2019) による個々の星に対する選択関数を、モデルの尤度（likelihood）に直接適用している。
  • フィッティング手順では、ガウス求積法を用いて、年齢、作用、および半径の測定誤差を周辺化している。極めて重要な点として、$J_0$ は各年齢において自由パラメータとして扱われている（4つのノットと半径方向の傾斜によってパラメータ化）。これにより、特定の加熱履歴や誕生時の分散を仮定することなく、形状指数 $b$ を測定でき、幾何学と加熱率を効果的に分離できる。

主な結果

1. 測定された指数: フィットの結果、拡散係数の傾きは $b = 0.51^{+0.06}_{-0.07}$ となった。

   • この値は、中間面限定層（$b=0.5$）の予測と一致しているが、体積充填型の予測（$b=1$）とは強く矛盾している。
   • 合成テストにより、本手法が偏りなく正しい指数を回収できること、および2つの幾何学的形状を $5.6\sigma$ の有意差で分離できることが確認された。
   • 結果は異なる半径ビンや年齢カットに対しても頑健であり、仮定された作用誤差の変化に対しても堅牢である。

2. 実効的な散乱体質量:

   • 測定された自由スケール（$J_0$）と導出された指数（$b \approx 0.5$）を用いて、加熱振幅を計算した。
   • 2–4 Gyrの年齢ビンにおいて、加熱振幅を現在の分子表面密度と比較すると、実効的な散乱体質量は $2.7 \times 10^6 M_\odot$ となる。
   • この値は、独立した雲カタログや質量関数（例：Miville-Deschênes et al. 2017）から導出された範囲（$1.4\text{--}3.4 \times 10^6 M_\odot$）に収まっている。
   • より古い星については、現在のガス密度と比較した場合、見かけ上の実効質量が増加する。しかし、過去数ギガ年にわたる星形成史（ひいてはガス密度の減少）による補正を行うと、すべての年齢ビンにおける実効質量は $1.9\text{--}2.8 \times 10^6 M_\odot$ の範囲に収束する。

3. 加熱履歴:

   • 自由スケールは $J_0 \propto \tau^{1.12}$ として年齢とともに増大する。これは、冷たい誕生分散を伴う一定の加熱率から予想される $\tau^{0.67}$ よりも速く、過去に円盤がより強く加熱されていたことを示している。

意義と主張
本論文は、垂直作用量の分布の「形状」を利用することで、垂直速度分散の「平均」や「分散」のみを分析する従来の手法を凌駕し、銀河円盤の散乱体の所在に関する長年の論争を解決したと主張している。

• 散乱体の局在化: 測定された指数 $b \approx 0.5$ は、垂直加熱の主要なメカニズムが銀河中間面に局在していることの直接的な証拠を提供している。これにより、体積充填型のメカニズム（ダークマター・サブハロ、曲げ波、あるいは直接的な螺旋摂動など）が、拡散的加熱の支配的な源であるという説を否定している。著者らは、体積充填成分が加熱に寄与するのは7%未満であると推定している。
• 散乱体の特定: 中間面の幾何学と導出された実効質量（$\sim 10^6 M_\odot$）の組み合わせにより、**巨大分子雲（GMC）**が散乱体であることを特定している。質量はGMC質量関数の上限域と一致しており、幾何学も分子ガスの薄い層と一致している。
• 従来の不一致の解決: なぜ以前の研究（例：Binney & Lacey 1988）が $b=1$（体積充填型）を推論したのかを説明している。それらの研究は、混合された年齢集団や速度分散を分析しており、固定された年齢における分布の形状を解明していなかった。著者らは、異なる加熱スケールを持つ集団を混合すると、分布の裾が人工的に平坦化され、たとえ真の幾何学が薄い層（$b=0.5$）であっても、体積充填型の幾何学（$b \approx 1$）を模倣してしまうことを示している。
• メカニズムとエネルギー源: 著者らは、散乱自体は（雲を介して）中間面で行われるが、エネルギーの起源は（螺旋構造などによって）面内の運動から供給されている可能性があることを明確にしている。指数は、エネルギーの究極の源ではなく、エネルギーが「再指向（redirection）」される場所を測定している。

結論として、本論文は、天の川銀河の円盤が、巨大分子雲の質量を持つ進化する集団によって、平面付近で加熱されていることを確立しており、垂直作用量の分布の形状が、散乱因子の空間分布を知るための頑健でモデルに依存しない診断指標であることを示している。