[論文レビュー] Lane-Emden Equation: perturbation method

Z. F. Seidov|arXiv (Cornell University)|Feb 5, 2004

Mathematics and Applications参考文献 3被引用数 4

ひとこと要約

本稿では、既知の解が存在する $ n = 0, 1, 5 $ の場合を越えて、任意の多相指数 $ n $ に対してラン・エムデン方程式を数値的摂動法で解く手法を提示する。$ \delta = n - n_0 $ のべき級数に解を展開することで、次元なし半径 $ X $、質量 $ \mu $、中心密度対平均密度比 $ \rho_c/\rho_m $ を高い精度で計算し、既知の数値データと一致することを確認した。また、$ 0 < n < 5 $ の範囲で有効であることを示した。

ABSTRACT

The perturbation method is applied to numerical solution of the Lane-Emden Equation of arbitrary index n, and the global parameters of polytropes are found as function of polytropic index n.

研究の動機と目的

• 任意の多相指数 $ n $ に対して、既知の解析的解が得られる $ n = 0, 1, 5 $ の3つの場合を越えて、ラン・エムデン方程式（LEE）を解くための数値的摂動法を開発すること。
• 次元なし半径 $ X $、質量 $ \mu $、中心密度対平均密度比 $ \rho_c/\rho_m $ といった、全構造的パラメータを $ n $ の関数として高精度に計算すること。
• ジャバール（1984年）、ホレッド（1987年、1990年）など、既存の文献値と比較することで、摂動法の数値的妥当性を検証すること。
• 特に $ n \approx 5 $ の近傍で摂動項が元の解と同程度の大きさをとる領域において、この手法の有効性と限界を評価すること。

提案手法

• 解を $ \delta = n - n_0 $ のべき級数に展開し、$ n_0 \in \{0, 1, 5\} $ として2次まで $ y = y_0 + \delta y_1 + \delta^2 y_2 $ とする。
• 非線形項 $ y^n $ の摂動展開を用いて、$ y_0 $、$ y_1 $、$ y_2 $ に対する連立2階常微分方程式を導出し、式（6）～（8）を得る。
• 補助変数 $ z_i = -x^2 y_i' $ を導入することで、6つの1階常微分方程式系に変換し、Mathematica の NDSolve を用いて安定な数値積分を実行する。
• ラン・エムデン解の既知のテイラー級数を用いて、$ x = 0 $ 近傍での初期条件を高次（$ x^{24} $ まで）の級数展開で得る。
• 中間点 $ x = X_0 $ で解をマッチングし、$ \mu = -X^2 y'(X) $ を用いて $ X $、$ \mu $、$ \rho_c/\rho_m $ を決定する。
• ジャバール（1984年）およびホレッド（1987年、1990年）の発表値と比較することで、計算された $ X(n) $、$ \mu(n) $、$ \rho_c/\rho_m $ の妥当性を検証し、6桁程度の一致を確認した。

実験結果

リサーチクエスチョン

• RQ1 任意の $ n \in (0,5) $ に対して、数値的摂動法が多相体の全構造パラメータを正確に計算できるか、特に $ n \approx 5 $ の近傍でも有効か？
• RQ2 $ n $ が 0、1、5 に近い場合に、2次摂動補正項 $ y_1 $、$ y_2 $ が $ X(n) $ および $ \mu(n) $ をどれほど正確に予測できるか？
• RQ3 摂動項 $ y_1 $、$ y_2 $ が元の解 $ y_0 $ と同程度の大きさに達する領域では、摂動法がどの程度有効に保たれるか？
• RQ4 計算された $ X(n) $、$ \mu(n) $、$ \rho_c/\rho_m $ の値は、ジャバール（1984年）およびホレッド（1987年、1990年）の既存の数値結果とどの程度一致するか？
• RQ5 特に $ n \approx 5 $ の場合に、大域的な $ x $-区間を渡る積分において、摂動スキームの数値的安定性と収束性はいかがなものか？

主な発見

• $ n = 3.1 $ の場合、計算された半径 $ X(n=3.1) = 7.3073764235 $ はジャバールの値と6桁一致し、ホレッドの値とも5桁一致した。
• $ n = 2.9 $ の場合、$ X(n=2.9) = 6.5273727715 $ はジャバールの結果と6桁一致し、ホレッドの値とも5桁一致した。
• 質量パラメータ $ \mu(n=3.1) = 1.9897152342 $ はホレッドの推定値 $ 1.9897028553 $ と一致し、差は $ 10^{-5} $ 未満であった。
• $ n = 3.6 $ の場合、$ X(n=3.6) = 10.2979713253 $ はジャバールの値と6桁一致し、ホレッドの値とも5桁一致した。
• $ n = 4.6 $ の場合、$ X(n=4.6) = 40.0619871651 $ はジャバールの値 $ 40.4132296343 $ と比べて $ 0.86\% $ の差で、$ \mu(n=4.6) = 1.7306833002 $ はホレッドの値 $ 1.7307484954 $ と $ 3.8 \times 10^{-5} $ の差で一致した。
• 本手法は $ n \in (0,5) $ 全体にわたり $ \mu(n) $ と $ X(n) $ を正確に計算でき、テストされたほとんどの値で $ \mu $ の相対誤差が $ 10^{-4} $ 未満であった。これにより、天体物理学的モデル化における実用的有用性が裏付けられた。