[論文レビュー] Revisiting the quasinormal modes of the Schwarzschild black hole: Numerical analysis

Luis A. H. Mamani, Angel D. D. Masa+2 more|arXiv (Cornell University)|Jun 7, 2022

Astrophysical Phenomena and Observations参考文献 57被引用数 28

ひとこと要約

本稿では、スピン-0、1/2、1、3/2、2、5/2の場に対して、擬スペクトル法と漸近的反復法（AIM）を用いてシュバルツシルトブラックホールの準正規モード（QNMs）を再検討する。スピン-5/2場の完全に虚数のQNM周波数の初回の数値計算が報告され、スピン-1/2および3/2場についても解析的結果と整合することを確認した。また、高次のオーバートーンモードにおいても、両手法が非常に精度が高く、相補的な性質を有することを検証した。

ABSTRACT

We revisit the problem of calculating the quasinormal modes of spin $0$, $1/2$, $1$, $3/2$, $2$, and spin $5/2$ fields in the asymptotically flat Schwarzschild black hole spacetime. Our aim is to investigate the problem from the numerical point of view, by comparing some numerical methods available in the literature and still not applied for solving the eigenvalue problems arising from the perturbation equations in the Schwarzschild black hole spacetime. We focus on the pseudo-spectral and the asymptotic iteration methods. These numerical methods are tested against the available results in the literature, and confronting the precision between each other. Besides testing the different numerical methods, we calculate higher overtones quasinormal frequencies for all the investigated perturbation fields in comparison with the known results. In particular, we obtain purely imaginary frequencies for spin $1/2$ and $3/2$ fields that are in agreement with analytic results reported previously in the literature. The purely imaginary frequencies for the spin $1/2$ field are exactly the same as the frequencies obtained for the spin $3/2$ field. In turn, the quasinormal frequencies for the spin $5/2$ perturbation field are calculated for the very first time, and purely imaginary frequencies are found also in this case. We conclude that both methods provide accurate results and they complement each other.

研究の動機と目的

• シュバルツシルト時空におけるすべての整数および半整数スピン場（s = 0, 1/2, 1, 3/2, 2, 5/2）の準正規モード（QNMs）を数値的に計算すること。
• 摂動方程式の固有値問題を解くために、擬スペクトル法と漸近的反復法（AIM）の精度および性能を比較すること。
• 既存の結果を拡張し、特にこれまで数値的に研究されていなかったスピン-5/2場の高次のオーバートーンQNM周波数を計算すること。
• 特にスピン-1/2および3/2場における完全に虚数の周波数について、既知の解析的および文献上の結果と照合し、数値的妥当性を検証すること。
• 一般の2階線形微分方程式（QNM問題に関連）を解くために、AIMを実装したオープンソースのソフトウェアパッケージの開発および公開すること。

提案手法

• 擬スペクトル法は、直交多項式（例：チェビシェフ）の基底における波動関数の展開を用い、微分固有値問題を行列固有値問題に変換する。
• 漸近的反復法（AIM）は、QNMのシュレーディンガー型微分方程式を、反復関係を繰り返し適用することで固有値を特定する手法を用いる。
• 径数座標をトロイツォー座標 $ r_* $ に変換し、各スピン $ s $ に対してポテンシャル $ V_s(r) $ を導出する。その形は $ V_s(r) = f(r) \left[ \frac{\ell(\ell+1)}{r^2} + \frac{(1-s^2)2M}{r^3} \right] $ である。
• 境界条件として、ホライズン（$ r \to 2M $）では入射波、空間無限遠（$ r \to \infty $）では放射波を要求し、複素離散固有周波数 $ \omega $ が得られる。
• 数値解は、一貫性の確認のため、既知の文献結果（例：Leaverの連分数法）と照合される。
• 一般の2階線形微分方程式に適用可能なAIMのオープンソース実装が開発され、一般利用可能に公開された。

実験結果

リサーチクエスチョン

• RQ1 擬スペクトル法と漸近的反復法は、シュバルツシルト時空におけるすべてのスピン場に対して一貫性があり、高精度な準正規モード周波数をもたらすか？
• RQ2 スピン-1/2および3/2場における完全に虚数の準正規周波数は、解析的に予測された通りに、数値的に確認できるか？
• RQ3 漸近的反復法は、スピン-5/2を含むすべてのスピン場について、高次のオーバートーンモードを信頼性高く計算できるか？
• RQ4 高角運動量および高オーバートーン数モードに適用した際、両手法の精度および収束性の挙動はいかがなものか？
• RQ5 QNM固有値問題を解く一般用途のオープンソースソフトウェアパッケージをAIMに基づいて開発することは可能か？

主な発見

• 擬スペクトル法と漸近的反復法の両方が、非常に高い精度の準正規モード周波数をもたらし、互いに良好な一致が得られ、既知の文献結果とも一致した。
• スピン-1/2および3/2場では、完全に虚数の準正規周波数が数値的に確認され、解析的予測と正確に一致した：$ \ell=0 $ の場合、$ M\omega_n = i(n+1)/4 $（$ n=0,1,2,\dots $）。
• シュバルツシルト時空におけるスピン-5/2場の準正規モードの初回の数値計算が達成され、解析的期待と整合する完全に虚数の周波数が得られた。
• 高次のオーバートーンモード（$ n \geq 1 $）に対しては、漸近的反復法が選択肢となる。同法は、同じパrameter範囲で擬スペクトル法が捉えられないモードを捉えることができる。
• AIMを実装したオープンソースソフトウェアパッケージは、成功裏に開発および検証され、一般に重力物理学における2階線形微分方程式に応用可能となった。
• $ \ell=1 $、$ n=1 $ の場合、唯一AIMが解を提供した。これは、擬スペクトル法が特定の高オーバートーンモードを捉えられない限界を示しており、一方で両手法が $ \ell=1, n=0 $ モードを再現した。