[論文レビュー] Markov Chain Monte Carlo methods applied to measuring the fine structure constant from quasar spectroscopy

Julian A. King, D. Mortlock+2 more|ArXiv.org|Oct 14, 2009

Galaxies: Formation, Evolution, Phenomena被引用数 3

ひとこと要約

この論文は、クェーサー吸収線からの微細構造定数 α の宇宙論的変動を測定するためのスペクトルフィッティングソフトウェア VPFIT の信頼性を検証するために、マルコフ連鎖モンテカルロ法（MCMC）を適用している。結果として、VPFIT は Δα/α 及びその不確実性の一貫性があり正確な推定値をもたらすことが確認され、α の変動に関する先行研究の妥当性が裏付けられる。

ABSTRACT

Recent attempts to constrain cosmological variation in the fine structure constant, alpha, using quasar absorption lines have yielded two statistical samples which initially appear to be inconsistent. One of these samples was subsequently demonstrated to not pass consistency tests; it appears that the optimisation algorithm used to fit the model to the spectra failed. Nevertheless, the results of the other hinge on the robustness of the spectral fitting program VPFIT, which has been tested through simulation but not through direct exploration of the likelihood function. We present the application of Markov Chain Monte Carlo (MCMC) methods to this problem, and demonstrate that VPFIT produces similar values and uncertainties for (Delta alpha)/(alpha), the fractional change in the fine structure constant, as our MCMC algorithm, and thus that VPFIT is reliable.

研究の動機と目的

• クェーサー吸収線からの Δα/α の推定に広く用いられるスペクトルフィッティングツール VPFIT の信頼性を評価すること。
• 一部の研究で懸念されたように、VPFIT の最適化アルゴリズムが真の最尤解に収束しない可能性を検証すること。
• MCMC 法を用いて、VPFIT が生成する統計的不確実性およびパラメータ推定値の独立的検証を行うこと。
• Keck と VLT の Δα/α 測定値の乖離が、アルゴリズム的欠陥によるものか、物理的変動によるものかを特定すること。
• 他のパラメータ（例：吸収断面積）が正規分布でない場合でも、VPFIT が報告する Δα/α の不確実性が統計的に妥当であることを確認すること。

提案手法

• 吸収線系のモデルパラメータ（Δα/α を含む）の事後分布をサンプリングするために、メトロポリス＝ハスティングス MCMC アルゴリズムを用いた。
• 事後分布を定義するために、尤度関数 L(𝐱) ≡ exp(−χ²(𝐱)/2) を使用した。ここで χ² は標準的な自由度補正済みカイ二乗統計量である。
• Voigtプロファイルを用い、1つの遷移に対して最大3成分まで含める形で、クェーサー3例（LBQS 2206−1958、LBQS 0013−0029、Q 0551−366）の複数の吸収線（例：Si ii、Al iii、Fe ii、Mg i）をフィットした。
• 収束診断を実施し、Δα/α およびその他のパラメータの周辺事後分布を検討して、正規性および信頼性を評価した。
• MCMC で得られた Δα/α 及びその不確実性を、VPFIT の結果と直接比較して一貫性を検証した。
• MCMC 結果の重み付き平均を用いて、共通の α 値を仮定して、複数の吸収系にわたる統合された Δα/α を推定した。

実験結果

リサーチクエスチョン

• RQ1 VPFIT の最適化アルゴリズムは、クェーサー吸収線分光データにおける Δα/α の真の最尤解に信頼性を持って収束するか？
• RQ2 他のモデルパラメータ（例：吸収断面積）が非正規分布を示す場合でも、VPFIT が報告する Δα/α の統計的不確実性は正確か？
• RQ3 MCMC 法は、複雑なスペクトルフィッティング状況下で VPFIT のパラメータ推定値の信頼性を独立的に確認できるか？
• RQ4 Keck と VLT の Δα/α 測定値の乖離は、VPFIT のアルゴリズム的欠陥によるものか、それとも実際の宇宙論的変動によるものか？
• RQ5 非正規分布を示すパラメータ分布が Δα/α の推定値の信頼性に与える影響はどの程度か。また、Δα/α 自体は概ね正規分布に近いか？

主な発見

• 3つのクェーサー吸収線系（LBQS 2206−1958、LBQS 0013−0029、Q 0551−366）における MCMC の Δα/α 結果は、VPFIT の結果と一貫しており、VPFIT の信頼性が裏付けられた。
• 吸収断面積やドーピラー幅パラメータの事後分布に顕著な非正規性が見られたが、Δα/α の周辺事後分布は概ね正規分布に近い。
• MCMC 法は、スペクトルフィッティングにおける潜在的なアルゴリズム的欠陥を効果的に同定・是正でき、VPFIT の性能が妥当であることを裏付けた。
• 3つの吸収系からの統合的 Δα/α 評価は (−0.74 ± 0.59) × 10⁻⁵ であり、α の変化がないという仮説と統計的に整合的であった。
• 他のパラメータが非正規的または poorly constrained であっても、VPFIT が報告する Δα/α の不確実性が信頼できることが確認された。
• Murphy ら（2004）による初期の Δα/α ≈ (−0.57 ± 0.11) × 10⁻⁵ の検出結果が、VPFIT の最適化アルゴリズムの失敗によるものではないことが、本研究で裏付けられた。