[論文レビュー] Rigorous "Rich Argument" in Microlensing Parallax

Andrew Gould|arXiv (Cornell University)|Feb 3, 2020

History and Developments in Astronomy参考文献 3被引用数 4

ひとこと要約

本論文は、離散的縮退する微小レンズ効果のパラメータ解の相対尤度について、厳密な解析的導出を提供する。尤度比は、はるかに大きな銀河的星数要因 $ H_i $ とヤコビアン要因 $ B_i = D_{L,i}/\pi_{{\rm E},i} $ に依存しており、これは『リッチの議論』を幾何的事前分布として形式化する。主な貢献は、大パラメータ解の抑制を表すヒューリスティックな $ \pi_{\rm E}^{-2} $ の一般化と正当化を可能にする閉形式の式を得たことである。この式は、微小レンズ効果パラメータ解のベイズ解析における堅牢な一貫性チェックを提供する。

ABSTRACT

I show that when the observables $(\vec π_{ m E},t_{ m E},θ_{ m E},π_s,\vec μ_s)$ are well measured up to a discrete degeneracy in the microlensing parallax vector $\vec π_{ m E}$, the relative likelihood of the different solutions can be written in closed form $P_i = K H_i B_i$, where $H_i$ is the number of stars (potential lenses) having the mass and kinematics of the inferred parameters of solution $i$ and $B_i$ is an additional factor that is formally derived from the Jacobian of the transformation from Galactic to microlensing parameters. The Jacobian term $B_i$ constitutes an explicit evaluation of the ``Rich Argument'', i.e., that there is an extra geometric factor disfavoring large-parallax solutions in addition to the reduced frequency of lenses given by $H_i$. Here $t_{ m E}$ is the Einstein timescale, $θ_{ m E}$ is the angular Einstein radius, and $(π_s,\vec μ_s)$ are the (parallax, proper motion) of the microlensed source. I also discuss how this analytic expression degrades in the presence of finite errors in the measured observables.

研究の動機と目的

• 離散的縮退する微小レンズ効果パラメータ解の相対尤度を、厳密に形式化し、『リッチの議論』が大パラメータ解を不利益にする理由を示す。
• 物理的星数密度と幾何的ヤコビアン要因を組み込んだ、離散的パラメータ解の相対尤度の閉形式解析的式を導出する。
• 文献に見られるヒューリスティックな $ \pi_{\rm E}^{-2} $ の大パラメータ解抑制が、物理的および幾何的事前分布の結果として正当化され、一般化されることを検証する。
• 特に $ \theta_{\rm E} $ が測定されない点レンズ事象において、数値的ベイズ解析の整合性チェックを提供する。

提案手法

• 相対尤度 $ P_i \propto K H_i B_i $ を導出する。ここで $ H_i $ は解 $ i $ の物理的パラメータと一致する銀河的星の数であり、$ B_i = D_{L,i}/\pi_{{\rm E},i} $ は微小レンズパラメータへの変換におけるヤコビアン要因である。
• レンズ-源相対運動の分布が解空間上で概ね線形であると仮定し、小さな誤差バーの下で有効性を保証する。
• 物理的星数の銀河的事前分布をモデル化し、パラメータ空間における幾何的歪みを補正するためのヤコビアン項を含める。
• 良好に測定された観測量にこの形式を適用し、$ B_i $ に現れる $ \pi_{\rm E}^{-1} $ の依存性が、『リッチの議論』の $ \pi_{\rm E}^{-2} $ ヒューリスティックを一般化していることを示す。
• $ \theta_{\rm E} $ が測定されない場合を検討し、$ D_L $ がベイズ的手法と銀河モデルを用いて推定される場合でも、式が依然有効であることを示す。
• 解析的結果が、完全な事後分布が計算された場合でさえも、数値的ベイズ統合の sanity check として機能することを示す。

実験結果

リサーチクエスチョン

• RQ1 離散的縮退する微小レンズ効果パラメータ解の相対尤度を、閉形式で厳密に導出する方法は何か？
• RQ2 大パラメータ解を抑制する『リッチの議論』の物理的および幾何的起源は何か？
• RQ3 文献に見られるヒューリスティックな $ \pi_{\rm E}^{-2} $ の抑制と、ヤコビアン要因 $ B_i = D_{L,i}/\pi_{{\rm E},i} $ との関係は何か？
• RQ4 この解析的式が破綻する状況は何か？ また、何時、数値的ベイズ手法に頼る必要があるか？
• RQ5 この形式は、$ \theta_{\rm E} $ が測定されない場合に、微小レンズ効果の数値的ベイズ解析の整合性チェックとして機能するか？

主な発見

• 縮退する微小レンズ効果パラメータ解の相対尤度は $ P_i = K H_i B_i $ で与えられ、ここで $ H_i $ は解 $ i $ の物理的パラメータと一致する銀河的星の数であり、$ B_i = D_{L,i}/\pi_{{\rm E},i} $ は幾何的ヤコビアン要因である。
• $ B_i $ に現れる $ \pi_{\rm E}^{-1} $ の依存性は、ヒューリスティックな $ \pi_{\rm E}^{-2} $ の大パラメータ解抑制を一般化・形式化し、『リッチの議論』を幾何的事前分布として正当化する。
• 解空間上で相対運動の分布が線形であるという仮定のもとで導出が成立し、小さな誤差バーの下ではこれが妥当である。
• 解析的結果は、$ \theta_{\rm E} $ が測定されず $ D_L $ が推定される点レンズ事象において特に、数値的ベイズ解析の堅牢な整合性チェックを提供する。
• $ \theta_{\rm E} $ が測定されない場合でも、$ D_L $ がベイズ的銀河モデルを用いて推定される限り、この式は依然有効であり、解析的公式が相互確認の手段として機能する。
• この手法は、パラメータ解の尤度分布の正規性からの逸脱、特に非対称または弧状の尤度等高線の状況に最も敏感であり、そのような場合には数値的評価が必要となる場合がある。