[論文レビュー] How Cumbersome is a Tenth Order Polynomial?: The Case of Gravitational Triple Lens Equation

Sun Hong Rhie|arXiv (Cornell University)|Feb 15, 2002

Relativity and Gravitational Theory参考文献 1被引用数 13

ひとこと要約

本稿では、重力的三重レンズ方程式が1つの複素変数における10次多項式に簡略化可能であることを示しており、Numerical Recipes や Mathematica などの標準的な計算ツールを用いて効率的な数値解法が可能である。著者らは多項式の明示的な係数を提供しており、従来の複雑さの認識とは対照的に、この系は取り扱えることが証明されている。

ABSTRACT

Three point mass gravitational lens equation is a two-dimensional vector equation that can be embedded in a tenth order analytic polynomial equation of one complex variable, and we can solve the one variable equation on the source trajectories using recipies for Fortran or $C$ (portable for $C$++ or $C_{jj}$) in Numerical Recipes, or using packages such as Mathemetica, Matlab, etc. This ready solvability renders fitting microlensing light curves including triple lenses a normal process, and such was done in a circumbinary planet fit for MACHO-97-BLG-41. Subsequently, there was a claim that converting the triple lens equation into the analytic equation was rather cumbersome, and the impressionable judgement has caused an effect of mysterious impedance around the perfectly tractable lens equation. There are judgements. Then, there is nature. We looked up for one of the quantities of highest precision measurements: electron $g$-factor correction $a_e \equiv g/2-1$. The current best experimental values of $a_e$ agree to eight significant digits with the theoretical value, and the theoretical calculation involves more than one thousand Feynman diagrams -- many orders of magnitude messier than the triple lens equation coefficients. We seem to have only choice to be compliant to nature and its appetite for elegant mess and precision numerics. In fact, the triple lens equation coefficients take up less than a page to write out and are presented here for users' convenience.

研究の動機と目的

• 重力的三重レンズ方程式を解くことが計算的に取り扱えないのだとする誤解を解きほぐすこと。
• 三重レンズ方程式に相当する10次多項式の明示的かつ解析的に導出された係数を提供すること。
• この系が、扱いやすい代数的構造を持つため、標準的な数値的手法に適していることを示すこと。
• 三重レンズ方程式の複雑さは、電子の g 因子補正のような高精度量子場理論計算ほどではないことを強調すること。
• 特に宇宙望遠鏡観測による系外惑星探査において、正確な解析解の使用を提唱すること。

提案手法

• 代数的代入を用いて共役変数を消去することで、三重レンズ方程式を1変数複素多項式に再定式化する。
• レンズ方程式を、複素変数 z および源位置 ω に関する多項式 G と H の累乗を含む形に変換する。
• 得られた方程式を、∑cff(k)·z^k = 0 の形で表し、cff(k) はレンズパラメータの対称多項式から導かれる係数である。
• 係数 H_ijk は組み合わせ的恒等式を用いて計算され、次数10まで明示的にリストアップされており、コードへの直接実装が可能である。
• Laguerre の方法などの既知の数値根発見アルゴリズム（Numerical Recipes より）を活用し、多項式の安定な解法を実現する。
• 対称的ケースのテストにより、既知の画像挙動と一致し、先行するミクロレンズ研究の結果と比較することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 一部の研究者が主張するように、重力的三重レンズ方程式は本当に解くことのできない複雑さを有するのか、それとも解ける多項式形に簡略化可能なのか？
• RQ2 三重レンズ方程式に相当する10次多項式の明示的代数的係数は何か？
• RQ3 三重レンズ方程式を解く複雑さは、電子の g 因子補正のような他の高精度理論的計算と比べてどの程度か？
• RQ4 標準的な数値ライブラリおよび根発見アルゴリズムは、ミクロレンズ光曲線モデリングのための完全な三重レンズ系を効果的に解けるか？
• RQ5 三重レンズ方程式の解析解を用いることで、ミクロレンズ惑星探査の精度と信頼性はどの程度向上するのか？

主な発見

• 重力的三重レンズ方程式は、複素変数 z における単一の10次多項式に簡略化可能であり、解析的に解ける。
• 10次多項式の係数は明示的に導出されており、1ページ未満に収まり、境界項として H_39 = 1 および H_30 = c^3 が与えられている。
• 三重レンズ多項式の複雑さは、1000以上のフェニマン図を含む電子 g 因子の計算よりもはるかに単純である。
• 本手法により、ミクロレンズ光曲線の強固なフィッティングが可能であることが、MACHO-97-BLG-41 の circumbinary 行星解析で実証された。
• Laguerre の方法とニュートン移動法による画像曲線追跡は、カウスティック領域外では実行可能であるが、カウスティック付近ではヤコビアン行列式が大きいため失敗する。
• 有限光源効果は、画像平面の画素化と、各画素ごとの輝度プロファイルの重み付けにより組み込むことができ、カウスティック通過時の正確な光曲線モデリングが可能になる。