Examining Active Galactic Nucleus UV/Optical Variability beyond the Simple Damped Random Walk. II. Insights from 22 yr Observations of SDSS, PS1, and ZTF

Weixiang Yu (于伟翔), Gordon T. Richards, John J. Ruan, Michael S. Vogeley, Franz E. Bauer, and Matthew J. Graham

Weixiang Yu et al 2025 ApJ 992 130(新しいタブで開きます)

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

※ 出典のない図は全て紹介する論文より引用したものです。

2.の論文の概要

22年にわたるAGNのUV/Opticalライトカーブに対するモデリングから以下の2点が分かった。

  1. モデルパラメータとクエーサーの物理的な性質の相関
    →ライトカーブモデリングのみで(分光観測をしなくても)物理的な性質を調べられる。
    (1.の論文で調べきれなかったことに再挑戦したみたい。)
  2. あるタイプのライトカーブを持つクエーサーの光度変動は2つの要因で説明できる。

AGN（Active Galactic Nuclei） について

AGN:母銀河全体を凌駕するほど強力な電磁波を放射している銀河中心部領域

AGN structure

中心に超大質量ブラックホール(Black Hole)、そのまわりに降着円盤(Accretion Disk)がある。 ガスやダストなどの物質がSMBHに降着することで開放される莫大な重力エネルギーが強力な電磁波の放射源と考えられている。

AGNの特徴
UV/Opticalのフラックスが数日~数十年で10%程度変動する。
←物理的なメカニズムは未解明
←ライトカーブのモデリングで、物理的な性質とモデルパラメータの相関関係を調査

https://heasarc.gsfc.nasa.gov/docs/objects/agn/agn_model.htmlより引用

UV/Optical ライトカーブのモデリング

Damped Random Walk(DRW) モデル

DRWモデルの式:\[\frac{dx}{dt}+\alpha x=\beta \epsilon(t)\] \(x\)は等級(mag)、\(\epsilon(t)\)はガウシアンホワイトノイズ

DRWモデルでは、ライトカーブの時間変動の振幅(\(\sigma\))、 タイムスケール(\(\tau\))それぞれと エディントン比 \(L/L_{\text{Edd}}\)、 ブラックホール質量 \(M_{\text{BH}}\) の間に相関があることが分かる。

DRW correlation

左の2つのグラフから、重たくて明るい星ほど、変動が小さい。

DRW correlation

重たい星ほど変動のタイムスケールが長い。(左のグラフ)
明るい星ほど変動のタイムスケールが短い。(右のグラフ)

B. C. Kelly et al. 2009より引用

DRWモデル：\(M_{\text{BH}}=10^8 M_{\odot} \rightarrow \tau\sim100 \text{days}\)
→ 時間変動は降着円盤の thermal fluctuation によって駆動されていることを示唆(Kelly et al. 2009; Sun et al. 2020; Burke et al. 2021)
※「thermal fluctuation time」は、diskの局所的な非平衡が熱平衡に戻るまでの時間

しかし、Kepler telescope（high cadence & S/N）での観測で見つかったより短時間での時間変動成分を考慮すると
DRWモデル × DHOモデル 〇 （B. C. Kelly et al. 2014; V. P. Kasliwal et al. 2017）

Damped Harmonic Oscillator(DHO)モデル

DHOモデルの式 \[\frac{d^2x}{dt^2}+2\xi\omega_0\frac{dx}{dt}+\omega_0^2 x=\sigma_{\epsilon}\epsilon(t)+\tau_{\text{perturb}}\epsilon\frac{d\epsilon}{dt}(t)\] \(x\)は等級(mag)、\(\epsilon(t)\)はガウシアンホワイトノイズ

上述のパラメータからAGNライトカーブの変動を定量的に表すパラメータがいくつか求められるが、今回は以下の4つを対象として扱う。

特に一つ目の\(\tau_{\text{decay}}\)は今回の論文で22年にわたる長期のライトカーブを扱ったことで調べられるようになった。

DHO_images

左側のグラフはデルタ関数的な外力を与えた際の減衰振動のグラフ。3組のパラメータに対するグラフが描画されている。
右側のグラフは、そのパラメータ値で上記のDHOモデル式からシミュレーションされたライトカーブ

W. Yu et al (2022)より引用

使用した測光データ

3つとも外れ値は除外、同一夜のデータは平均をとる。合計22年分の光度曲線が得られる。

recalibrationとphotometric offset の調整

異なる望遠鏡で得られたデータを一つにまとめるには、それぞれの望遠鏡の特性に起因するズレを補正してやる必要がある。 (ある望遠鏡では18等の星が、別の望遠鏡では18.1等になるといったふうに。)
SDSS-I以外の各サーベイに対して、標準星のライトカーブでガウス過程回帰を行い、pipeline photometric errorで説明しきれていないエラーを推定してその分を補正。
左図は等級に対するreduced-\(\chi^2\)の図。縦軸の値が1より大きいとエラーが過小評価されており、1より小さいとエラーが過大評価されている。
望遠鏡によってエラーバーの幅が大きく出やすかったり小さく出やすかったりするために、その程度を揃えているということかも？
ライトカーブのフィッティングの際に、エラーが大きい方に結果が引っ張られてしまうのを防げる。
SED fittingを行って、PS1、ZTFのライトカーブをSDSSのライトカーブに統合した。
左図は例としてランダムに選んだ3天体の統合したライトカーブ。（縦軸はgバンド等級）

対象クエーサーの絞り込み

Initial Sample

Clean Sample

Core Sample

Results

Structure Function

Structure Function(SF)：ある時間差に対応する2時点間の等級の差

クエーサーライトカーブのSF; \[SF(\Delta t) = A(\frac{\Delta t}{\Delta t_0})^{\gamma}\] \(\gamma\)がSFの形を特徴付ける。

DHO モデルの SF（青実線）がDRW モデルのSF（赤破線）よりも観測データのSF（黒点）をうまく説明できる。

Power Spectral Density Function

パワースペクトル密度関数：各振動数成分がどれくらいの強さで含まれているかを表す。

下図は 振動数×PSD と(補完された)ライトカーブの例

2成分ローレンチアンでfit (X線でのPSDの先行研究に倣う)
A：narrow + broad Lorentzian
B：うまくfitできず (for future work)
C：broad Lorentzian×2

変動の要因(Class C)

image of two component constituting light curve variability

Class CのPSDの2つのローレンチアンを以下のように物理的に解釈;

左上側の図ではSMBH付近のコロナから放射されたX線が降着円盤でリプロセスされて放射されている光と降着円盤そのものから放射されている光が描かれており、 それぞれが左下側のPSDの図では、低周波数側にピークを持つローレンチアンと高周波数側にピークを持つローレンチアンで説明されている。
accretion rate fluctuation は降着円盤を外側から内側に向かって伝播する。
これはX-ray binaryと少数のAGNのX線ライトカーブでの解析から得た発想 (Nowak 2000、McHardy et al.(2007)など)
※accretion rate fluctuationはmagnetorotational instabilityというものによって降着円盤の各小領域で局所的に生み出され、 それが内側に運ばれるにつれてどんどん増幅されるらしい。 （Ingram & van der Klis(2013)）

Discussion

tmp

accretion rate fluctuationとreprocessed X-ray variabilityを仮定すると、DHOモデルのパラメータとクエーサーの物理量の間の相関がどのように説明されるか。 Hagen et al. 2024