[论文解读] 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

一句话总结

本文使用数值方法——伪谱法和渐近迭代法（AIM）——重新研究了史瓦西黑洞的准正规模（QNMs），针对自旋为0、1/2、1、3/2、2和5/2的场。该研究首次实现了对自旋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）在求解微分方程特征值问题时的精度与性能
• 通过计算更高倍频的QNM频率，扩展现有结果，特别是对尚未被数值研究过的自旋5/2场
• 将数值结果与已知的解析结果及文献结果进行验证，尤其针对自旋1/2和3/2场的纯虚数频率
• 开发并发布一个开源软件包，实现用于一般二阶常微分方程（与QNM问题相关）的AIM

提出的方法

• 通过将波函数展开为正交多项式基（如切比雪夫多项式），应用伪谱法，将微分特征值问题转化为矩阵特征值问题
• 使用渐近迭代法（AIM）求解QNM的薛定谔型微分方程，通过迭代应用递推关系以检测特征值
• 将径向坐标变换为tortoise坐标 $ 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连分数法）进行比较，以进行一致性检验
• 开发并发布一个开源的AIM实现，用于一般性求解具有复边界条件的二阶常微分方程

实验结果

研究问题

• RQ1 伪谱法与渐近迭代法在史瓦西时空中对所有自旋场的准正规模频率是否能产生一致且准确的结果？
• RQ2 自旋1/2和3/2场的纯虚数准正规模频率是否在数值上得到确认，如解析预测所示？
• RQ3 渐近迭代法是否可可靠地扩展至计算所有自旋场（包括自旋5/2）的高倍频模式？
• RQ4 当应用于高角动量和高倍频数模式时，两种数值方法的精度与收敛性行为如何？
• RQ5 是否可行开发一个通用的、开源的软件包，基于AIM求解黑洞微扰理论中的QNM特征值问题？

主要发现

• 伪谱法与渐近迭代法均能产生高度准确的准正规模频率，两者之间以及与已知文献结果之间均表现出极佳的一致性
• 对于自旋1/2和3/2场，纯虚数准正规模频率在数值上得到确认，与解析预测完全一致：$ M\omega_n = i(n+1)/4 $（当 $ \ell=0 $，$ n=0,1,2,\dots $）
• 首次实现了对史瓦西时空中自旋5/2场的准正规模的数值计算，得到的频率为纯虚数，与解析预期一致
• 对于高倍频模式（$ n \geq 1 $），渐近迭代法为首选方法，其能捕捉到伪谱法在相同参数范围内未能发现的模式
• 基于AIM的开源软件包已成功开发并验证，可广泛应用于引力物理中二阶常微分方程的求解
• 对于 $ \ell=1 $，$ n=1 $，仅AIM能获得解，表明伪谱法在捕捉某些高倍频模式时存在局限性，而两种方法均能复现 $ \ell=1, n=0 $ 模式