Probing the Physical Origin of the Balmer Decrement in the Broad-line Region of Nearby Active Galactic Nuclei via Spectral Variability

Suyeon Son, Minjin Kim, Luis C. Ho, and Ruancun Li

Suyeon Son et al 2025 ApJ 995 37(新しいタブで開きます)

※ 月曜雑誌会での論文紹介のために作成したページです。

※ 出典のない図やグラフは紹介する論文から引用しています。ここにある記述には誤りも含まれているかもしれませんので、ご注意ください。

Abstract

AGNBLRからのBalmer Decrementの原因を探るために、SDSSで複数回分光観測されたlow-zのAGNのデータを用いて、fluxの時間的な変化を調べた。mean spectraを用いた解析で、Balmer decrementとAGN luminosityには相関がないことが明らかになり、Eddington比とは負の相関関係を持つことが分かった。一方で、個々のAGNのfluxの時間変動を調べるとluminosityとの間に更に強い負の相関関係があることが判明し、AGNのBalmer decrementを決定する上でluminosityの変化が重要であることが示唆された。Balmer decrementの時間進化と連続光のカラーを比べることで、AGN自身によるreddening(赤化)がBalmer decrementを決定する主要因ではなく、輻射輸送と励起のメカニズムが極めて重要であることが明らかになった。この結果はchanging-look AGNや、"little red dots"(LRDs)などのhigh-z AGNsの物理にも有益な洞察を与えることになる。(LRDsはBalmer decrementが極端に大きく、ダストの証拠として誤って解釈される。)

1. Introduction

AGNにおいて、水素の再結合線(Hα、Hβ、Hγ、...)はBroad Line Region(BLR)やNarrow Line Region(NLR)の大きな特徴の一つ。(降着円盤からの電離光子が水素ガスを電離して輝線が出ている。)
BL: Δv ~ 5000 km/s、Δλ ~ 100 Å、   NL: Δv ~ 1000 km/s、 Δλ ~ 20 Å

Emmaalexander unified agn
AGNのイメージ図(https://commons.wikimedia.org/wiki/File:Emmaalexander_unified_agn.pngより引用。リサイズしています。)
Type1,Type2 AGN の composite スペクトル
Type1,Type2 AGN それぞれの composite スペクトル(Shu Wang et al 2024 ApJ 966 128 Figure 1.より抜粋)
Hα、HβのBroad Line(Type1のみ),Narrow Line(Type1,2両方)が見えている。(Hγ、HδのBLも見える。)

Balmer Decrementはバルマー系列のfluxの強度比のことで、ダスト減光の指標として用いられていた。
特に、Balmer decrementはgasの温度や密度の条件に対する感度が比較的鈍いのでよく利用され、特にHα/Hβはそれぞれの輝線強度が強く、観測しやすいので使われることが多い。

NLR,planetary nebulae,HII regionのような密度の低い環境からの放射でBalmer Decrementが大きくなるのは、ダスト吸収によるreddening(赤化)が原因と考えられる。

一方で、BLRはガス密度が高い(~108-1011-3)ので、Balmer decrementがガス密度や温度に影響を受けやすい。
BLRのBalmer decrementはType1 AGNでは、Hα/Hβ ~ 2.7-3.5 と報告されていて、Photoionizing modelsで説明可能。
BLRのBalmer decrementの原因がdust reddeningなのか、他の要因によるものなのか議論が続いている。

continuum slopeとX線のphoton indexを調べて、Balmer decrementをdust extinctionの効果で説明できると主張する論文もあれば、dust吸収の兆候が小さいAGNでもBalmer decrementが大きくなると主張する論文もあり、dust吸収の効果と降着円盤の活動状態の効果を切り分けることが難しい。

そこでBalmer decrementの変動を調べることが新たな解決の手段となる。
しかし、Balmer decrementはAGNの活動状態にも依存することが分かっているので、その点も考慮することが必要。

その点で活動状態が劇的に変化するChanging-Look AGNは研究対象として適しているが、32個のCLAGNを用いた調査(S. Panda & M. Sneigowska(2024))ではBalmer decrementとAGN の活動状態の間に相関関係がないという結果が出てしまっているので、さらなる系統的な調査が必要
→ 本研究では多くのType 1 AGNサンプルを用いた系統的な調査を行う。

2. Data

SDSS DR16 quasar catalogからType 1 AGNsを選ぶ。

Balmer decrementの時間変動を見るために、複数epochの分光観測データがある天体を選ぶ。(ただし、BOSS分光器で撮られたスペクトルのみ使用)

z≤0.45の天体を選出(Hα線(6563Å)がスペクトルの中に入るように)

目で見てHβやHα付近のスペクトルが不自然なものは排除

Type 1 AGN 1403天体、1天体ごとに作ったスペクトルのペアが計69275ペア
(~86%くらいの天体は2epochから作った1ペアだけ、SDSSのreverbration mapping projectの対象になっていたいくつかの天体は1000ペア以上)

3. Spectral Fitting

PyQSFit(SDSS Quasarのスペクトルフィッティング用のPythonライブラリ)を利用。

Galactic extinctionを補正(E.L. Fitzpatrick(1999)のextinction lawとE.F. Schlafly & D.P. Finkbeiner (2011)のE(B-V)を利用)

スペクトルを静止波長系に変換(redshiftはB.W. Lyke et al.(2020)の値を利用)

continuum fitting

を取り除き、輝線成分を抽出

特に、母銀河成分の分離に最適化された手法を利用した。(W. Ren et al. 2024)
∵降着円盤から放射されるcontinuumと、バルマーラインのフラックスを正確に計測するために極めて重要だから。

輝線fit

最後に、目で見て(visually inspect)うまくfittingできていないものは除外

fittingプロセスの1例が下図

スペクトルフィッティングの図
Figure 1. ある天体のスペクトルfitの図

最終的なサンプル数は、1116 targets、スペクトルペアは44679組
~86%程度の天体はスペクトル2つだけ、いくらかの天体は2000以上のスペクトルの組をもつ

4. Results

4.1. Balmer Decrement versus AGN Properties

Balmer decrement(f/f)の平均値 : 3.66±0.95   >   3.1 (AGNのNLRでCase B recombinationを考えた場合(?)の値)
3%程度の天体は測定誤差かBLとNLの分離の失敗でBalmer decrementが2.7(普通のCase Bの場合)を下回った。

Balmer decrementと以下のAGNの物理量との相関関係を調べた。

最終的なサンプルの物理量は以下のように分布

解析には平均の物理量を利用。また、extinctionの原因が分かっていないので、一切の補正は行っていない。

Balmer decrementと各種物理量の相関関係は以下のとおり

Balmer decrementとAGNの物理量の相関関係
Figure 2. Balmer decrementとbolometric luminosity(左)、Edington比(右)の相関関係
density contourはlogスケール

いずれも先行研究と一致 (Balmer decrementとEdington比の相関関係は、low-luminosity AGNに限るともっと強い負の相関が見られる。)

4.2. Temporal Variation of Balmer Decrement

Balmer decrementの時間変化 Δf/f とAGNの性質との関係を調べる。(時間変化を調べるので、[OIII]を使って追加でflux calibrationを実施。)

Δf/f と Δf5100 の間には負の相関がみられた。(相関係数r = -0.24)

Balmer decrementの変化とbrightnessの変化の相関関係
Figure 3. Balmer decrementの変化とAGN brightnessの変化(Δ log f5100)の間の相関関係
density contourは対数スケール

直交距離回帰で分布を回帰すると、 Δlog(f/f) = -0.40 Δlog f5100
原点を通る直線になっているので、Balmer decrementの変化とbrightnessの変化の相関関係はAGNの明るさによらない。

5. Implications

5.1. Physical Origin of Balmer Decrement

4.2.でBalmer decrementの変化とAGN luminosityの変化に負の相関が見られることが分かったが、その原因が視線上のdustによる減光か否か調べるために、continuumのカラーの変動とBalmer decrementの変動の相関関係を調査↓

Balmer decrementの変動とカラーの変動の相関関係
Figure 4. Balmer decrementとcontinuum colorの変動の関係図
青色の3つの線はAGNのintrinsic extinctionから予想される関係
赤線と緑線はそれぞれstar-forming galaxyとAGNの減光曲線から予想される関係
黒色の破線が観測データを直交回帰した結果

赤化するほど(f5100/f3800が大きくなるほど)Balmer decrementは大きくなっている。

←extinctionがBalmer decrementの物理的な要因の可能性 (ただし、降着円盤とBLRが同程度のextinctionを受けているという前提)
しかし、天の川銀河star-forming galaxyAGNのextinction curveから予想されるBalmer decrementより遥かに大きい傾きが見て取れる。
したがって、dust extinctionのみがBalmer decrementの原因ではないとも考えられるが、AGNの"bluer when brighter"という性質にも注意が必要。
∵ AGNのfluxが下がる → Balmer decrementが大きくなる(4.2.の結果) & 赤化(log f5100 - log f3800が大きくなる) → Figure 4.のようになる。

←降着円盤とBLRが受けるextinctionの量が違うと考えて解釈できるかも
ただし、これはAGN luminosityの減少に対してtorusのcovering factorが系統的に増加して、かつdiskよりBLRの方がよりextinctionの影響を受けている場合の話。

← よりあり得るシナリオは、BLR中のHαとHβのcontinuumの変化に対する感度の違い。(Hαの方がHβよりも鈍い)
BLRのような高密度のガス中では、lineの光学的厚みの影響が大きくなってくる。(Hα:n=3→2、Hβ:n=4→2だから、Hαの方がガス中で吸収されやすい。)
Hαは衝突励起もあるので、AGN luminosityが下がってもHαはfluxの減少が鈍いという効果も。

実際、個別のType 1 AGNやChanging-Look AGNに対するBalmer decrementの変動を調査した先行研究でも以下のような傾向が見つかっている。

他にもBalmer decrementの変動とHα fluxの変動に正の相関があることも見つかった。

Balmer decrementの変化とHα fluxの変化の関係
Figure 5. Balmer decrementの変動とHα fluxの変動の関係

MBHを求める際に使われるLとL5100の関係式によるとLが増加するとL5100も増加するということになっているが、これはFigure 3.とFigure 5.に矛盾する。
← BLR中の輻射輸送効果と励起メカニズムの関係で生じているのだろう。
← Balmer decrementも考慮した関係式が以下のようになる。(luminosityの単位は erg s-1)

\[\log L_{5100} = 0.99 \log L_{H\alpha} - 0.48 \log(f_{H\alpha}/f_{H\beta}) + 1.98\]

さらに、この効果によってHαの等価幅の変動がBalmer decrementの変動と強い相関関係を示すことがわかる。(Figure 6.)

Hα輝線の等価幅とHα luminosity、Balmer decrementの変動の関係
Figure 6. 左図はHαのluminosityと等価幅の関係。赤い点はLRDs。右図はHαの等価幅の変動とBalmer decrementの変動の関係

物理量の変化ではないが、Balmer decrementとEdington比の間に弱い負の相関があることも分かった。(先行研究ともconsistentではある。)

5.2. Application to High-z AGNs

近年、James Webb宇宙望遠鏡による観測で"Little Red Dot(LRD)"と呼ばれるAGNの新たな種族が見つかっている。(D. D. Kocevski et al. 2023; R. Maiolino et al. 2024)
特に、V字型のSEDをもち、high-zに多くあることが特徴(Y. Harikane et al. 2023; J. E. Greene et al. 2024; V. Kokorev et al. 2024)

LRDは比較的大きいHαの等価幅とsteepなBalmer decrementを持つという特徴がある。(J. Matthee et al. 2024; M. Brooks et al. 2025; R. Maiolino et al. 2025)

Balmer decrementの変動とHαの等価幅の関係は本研究で扱ったサンプルにも見て取れる。(Figure 6.)
この関係がBLR中の輻射輸送効果と励起メカニズムによるものだった(Figure 5.)ことは、LRDsのようなhigh-z AGNの物理への示唆となる。(LRDsのsteepなBalmer decrementが赤化によるものではないかも)

6. Summary

SDSSから1000天体以上のlow-z AGNsを選び、multiepochスペクトルデータを用いて、1型AGNのBLRのBalmer decrementについて、その物理的な起源を調査した。

  1. mean spectraを用いた解析から、Balmer decrementはAGNのluminosityとの相関がなく、エディントン比とは弱い負の相関を持つことが分かった。(同様の結果が先行研究でも得られている。)
  2. 個々の天体のスペクトルを解析することで、Balmer decrementの時間変動がAGNのluminosityの時間変動と負の相関を持つことが分かった。→brightnessの変動がBalmer decrementを決めるうえで極めて重要
  3. continuum colorの変化とBalmer decrementの変動を比べると、先行研究の候補となるようなdust extinction curveでは極端なケースを考えても観測結果と合わなかった。→AGN中のintrinsicなreddeningはBalmer decrementの変動においては支配的な要因でないことを示唆
  4. BLRのBalmer decrementはBLRの高密度に起因するlineに対するoptical depthの効果や衝突励起が影響していそう。特にionizing continuumに対するHαとHβの感度の違いがBalmer decrementの変動に大きく影響しているだろう。

この研究でBLRのBalmer decrementは視線上のdustによる赤化の効果というよりは、主にBLRの物理状態によって決定されていることが明らかになった。
この解釈はchanging-look AGNやlittle red dotsのようなsteepなBalmer decrementを示すhigh-z AGNの物理についても有益な示唆を与えることになる。