具有拓扑缺陷的普拉默暗物质黑洞：影子、灰体因子、准正则模与热力学

《Annals of Physics》：Plummer dark matter black hole with topological defects: Shadow, greybody factors, quasinormal modes, and thermodynamics

　　

摘要

我们研究了一种嵌入于核球型 Plummer 暗物质晕并环绕弦云的静态球对称黑洞时空。该晕的密度轮廓通过相对论切向速度（旋转曲线）关系进入度规的构造，而非作为直接源项出现，我们对由此产生的几何结构求解爱因斯坦场方程。支撑物质是一种两组分各向异性流体，其能量–动量张量可加性地分解为晕部分与 Letelier 弦云部分。我们明确表明，以此方式得到的晕能量密度虽然以 Plummer 密度为种子，但与该种子 Plummer 密度并不相同：两者共享相同的 $r^{-4}$ 尾部，在 $r = r_0$ 处重合，而在 $\rho_0$ 的领头阶以因子 $(r^2 + r_0^2)/(2r^2)$ 相差。我们对物质内容进行验证，借助 Maple 进行符号运算，并在相关参数范围内进行数值计算，确认在外部区域 $r > r_h$ 内处处满足全部四个标准的逐点能量条件（零能条件、弱能条件、强能条件与优势能条件）。我们通过计算光子球与阴影半径以及利用 Gauss–Bonnet 定理得到的弱偏折角来分析其光学与动力学性质。我们也确定最内稳定圆轨道。然后将阴影预测与事件视界望远镜（EHT）对 M87* 与 Sgr A* 的测量结果进行对比，从而给出弦云张力与 Plummer 晕密度的参数上限。标量微扰通过有效势进行分析，灰体因子通过 Boonserm-Visser 方法给出上界，同时研究霍金辐射谱以及由切比雪夫谱方法求解并经三阶 WKB 交叉验证得到的准正则模频率。我们还进行了热力学分析，包括霍金温度、熵、热容、类 Smarr 关系与 Gibbs 自由能。在所有参数范围内热容始终为负，表明热力学不稳定，且不存在任何 Davies 型相变。总体而言，一个清晰的层级关系浮现出来：弦云参数对所有可观测量提供主导修正，而 Plummer 晕密度贡献次要效应。

引言

广义相对论（GR）迄今为止通过了所有高精度检验，从太阳系实验到 LIGO/Virgo 对引力波的直接探测 [1]，再到事件视界望远镜（EHT）获得的首张黑洞阴影图像 [2]。尽管取得了这些成就，仍有两个重大悬而未决的问题：星系尺度上暗物质的本质，以及拓扑缺陷——宇宙弦（CS）、全局单极子（GM）以及弦云（CoS）——在早期宇宙对称性破缺相变中产生的作用 [3]。

多条独立的观测证据支持暗物质晕与拓扑缺陷在天体物理黑洞附近同时存在。EHT 对 M87* [2] 与 Sgr A* [4] 的测量对事件视界附近的几何结构施加了严格约束，任何偏离纯 Schwarzschild 轮廓的部分都可归因于环境物质。星系旋转曲线、弱引力透镜巡天以及星系团的 X 射线观测提供了强有力的证据，表明星系中心存在暗物质晕，而这些区域正是超大质量黑洞的所在地 [5]、[6]。虽然对宇宙弦的直接天文探测仍在进行，但最新的 CMB 温度与偏振测量将无量纲弦张力限制在 $G\mu \lesssim 10^{-7}$ [3]，这一水平与本工作中探索的 $\alpha \in [0, 0.3]$ 参数范围完全一致。脉冲星计时阵列对宇宙弦网络产生的随机引力波背景的搜寻已经对弦云张力参数施加了约束，而近期的基于 EHT 的现象学研究已开始利用黑洞阴影观测来联合约束弦云与暗物质参数 [7]。综合而言，这些观测窗口表明，Plummer-CoS 黑洞模型不仅是一个数学练习，而是一个具有现象学意义的研究对象。

在研究规则黑洞或源于非线性电动力学、各向异性流体的黑洞时，一个核心问题是：这些物质构型能否在现实的天体物理环境中实现。实际上，黑洞从来都不是完全孤立的；它们嵌入在大尺度结构中，而这些结构往往由暗物质主导。来自星系与星系团的观测证据强烈支持暗物质晕的存在 [5]、[6]，尽管暗物质的微观本质仍不确定。为了对暗物质晕建模，星系动力学与宇宙学模拟中通常采用各种现象学密度轮廓，例如 Hernquist [8]、Einasto [9]、Dehnen [10] 以及 Dekel-Zhao 模型 [11]、[12] 提出的方案。值得注意的是，这些轮廓指定了密度分布，而将压强留作未指定项，为与黑洞解兼容的有效流体描述提供了灵活性。在本工作中，我们考虑另一种已知为 Plummer 轮廓 [13]、[14] 的暗物质晕轮廓，并研究一个被弦云环绕的黑洞解。Plummer 轮廓的密度函数 [13]、[14] 为：

$$\rho(r) = \frac{\rho_0}{[1 + (r/r_0)^2]^2}$$

该轮廓结合了有限的中心密度、陡峭的外部衰减 $\rho \sim r^{-4}$ 以及可解析处理的包络质量函数。Senjaya 等人 [15] 最近构造了精确的 Schwarzschild-Plummer 黑洞解，并分析了其光子动力学、eikonal 极限下的准正则模以及热力学性质 [15]。他们的度规函数为：

$$h_P(r) = \exp\left[-\frac{4\pi \rho_0 r_0^3}{r} \arctan\left(\frac{r}{r_0}\right)\right] - \frac{r_s}{r}$$

其中 $r_s = 2M$ 为 Schwarzschild 半径。该解捕捉了核球型暗物质晕对中心黑洞几何的引力效应，但未考虑在现实星系环境中可能与晕共存的拓扑缺陷。

拓扑缺陷为时空修正提供了独立的来源。Letelier [16] 首先提出的弦云代表了穿过时空的一维物体分布。其能量–动量张量由单个张力参数 $\alpha \in [0, 1)$ 参数化，该参数作为加性常数 $-\alpha$ 进入度规函数 [16]、[17]。另一方面，全局单极子引入一个乘性实体亏损因子 $(1 - 2\eta^2)$ [18]，而由 Nambu-Goto 作用量产生的宇宙弦产生类似的锥形亏损修正 [19]。这些缺陷在统一大理论相变中自然出现，并可能穿过星系中心超大质量黑洞的时空，而这些区域同时存在暗物质晕。

近年来，暗物质晕与拓扑缺陷对黑洞可观测量的联合效应引起了越来越多的关注 [7]、[20]、[21]、[22]、[23]、[24]。然而，大多数现有研究将两者效应分别处理：要么单独研究暗物质晕，要么将拓扑缺陷叠加在真空黑洞之上。将核球型暗物质轮廓与弦云缺陷同时纳入单一黑洞几何，并研究其完整现象学后果的工作，对于 Plummer 晕而言尚未开展。

黑洞热力学研究引起了相当大的兴趣，因为它深化了我们对黑洞及其与支配物理系统基本原理之间联系的理解。此外，它对量子引力理解具有重要意义。黑洞热力学的基础概念在开创性著作 [25]、[26] 中引入，其中建立了黑洞力学四定律，与经典热力学定律类比。后续研究表明黑洞发射热辐射 [27]、[28]，并拥有与其事件视界面积成正比的熵 [29]、[30]。在黑洞热力学中，已研究了多种类型的相变。例如，Davies 型相变源于热容发散 [31]，而 Hawking–Page 转变描述热 Anti-de Sitter（AdS）时空与黑洞之间的相变 [32]。其他例子包括极端相变 [33]、[34]，以及扩展相空间中的 Van der Waals 类行为，其中宇宙常数被视为热力学压强，黑洞质量被解释为焓 [35]、[36]。

黑洞通过阻尼振荡返回平衡态，其特征由其准正则模（QNMs）刻画，使这些复频率成为强场引力有价值的探针。准正则模的实部对应微扰的振荡频率，而虚部决定衰减速率 [37]、[38]。在数学上，准正则模表现为具有耗散边界条件——在事件视界处入射、在空间无穷远处出射——的类 Schr?dinger 波动方程的本征值 [39]、[40]。由于对于任意势函数精确解通常不可行，准正则模通常使用半解析或数值方法计算，例如 WKB 近似、Padé 求和、连分数或时域积分 [41]、[42]、[43]。WKB 方法最初由 Refs. [44]、[45]、[45] 应用于 Schwarzschild 黑洞，提供了一种相对简单但有效的半解析方法来研究黑洞微扰。后来，该技术被扩展到更一般的时空，包括旋转（Kerr）和带电（Reissner–Nordstr?m）黑洞 [46]、[47]。此后，准正则模在广泛的黑洞几何和周围物质场中被广泛研究，包括近期对规则黑洞和量子修正黑洞的研究 [48]、[49]。

在本文中，我们通过构造 Plummer-CoS 黑洞解来填补这一空白——该解通过将 Letelier 弦云参数 $\alpha$ 添加到 Ref. [15] 中的 Plummer-Schwarzschild 度规而获得——并研究其在七个相互关联领域内的性质：物质内容与能量条件、事件视界结构、光子测地线与阴影、弱引力透镜、类时测地线与 ISCO、标量微扰（灰体因子、准正则模与霍金辐射）以及热力学。我们的分析表明，弦云张力 $\alpha$ 与 Plummer 晕密度 $\rho_0$ 以相关但层级有序的方式影响所有可观测量：$\alpha$ 主导事件视界半径、光子球、阴影大小、ISCO、准正则模频率和热力学量的偏移，而 $\rho_0$ 提供次要的加性修正。本工作通过显式推导爱因斯坦场方程确立的一个关键结构特征是，该度规假设由各向异性流体支撑，其能量–动量张量清晰地分解为晕部分与 Letelier 弦云部分，后者精确再现 Letelier 的原始张量结构 [16] $T_{tt} = T_{rr} = -\alpha/(8\pi r^2)$，$T_{\theta\theta} = T_{\phi\phi} = 0$，采用第 2.2 节中固定的符号与场方程约定。我们在开篇强调，该源中的晕部分并非原始 Plummer 密度：Plummer 轮廓仅作为切向速度关系的输入，而爱因斯坦方程返回的密度与它相差一个可计算的因子，如第 2.3 节所量化。我们表明，组合物质源在整个外部区域 $r > r_h$ 内满足全部四个标准的逐点能量条件（NEC、WEC、SEC、DEC），该性质已通过 Maple 2024 进行符号与数值双重验证。该解的另一个结构特征是，对所有 $\alpha < 1$，它仅存在单一非简并事件视界，不存在极端或多视界构型，这是径向方程中不存在排斥势垒的结果；我们通过基于 $A(r)$ 严格单调性的视界唯一性闭式证明来强化这一观察。我们还表明，热容在所有参数取值下严格为负，确认 Plummer-CoS 黑洞在正则系综中热力学不稳定，且不存在 Davies 型相变。此外，我们所得阴影半径与 EHT 对 M87* 和 Sgr A* 的约束进行了对比，提供了经观测校准的 $(\rho_0, \alpha)$ 上限。

论文组织如下。第 2 节构造 Plummer-CoS 度规，从爱因斯坦场方程推导物质源，检验能量条件，分析渐近结构，并对全参数空间进行分类事件视界构型，给出闭式唯一性证明。第 3 节研究零测地线、光子球、黑洞阴影半径、光子有效径向力以及通过 Gauss–Bonnet 定理得到的弱偏折角，并将阴影预测与 EHT 对 M87* 和 Sgr A* 的数据进行对比。第 4 节处理类时测地线并确定 ISCO 半径。第 5 节推导标量微扰势，使用 Boonserm-Visser 方法计算灰体因子上界，给出霍金辐射谱，并通过 WKB 近似获得准正则模频率。第 6 节涵盖热力学分析：霍金温度、Bekenstein-Hawking 熵、热容、类 Smarr 关系与 Gibbs 自由能。第 7 节给出结论。全文使用自然单位 $G = c = \hbar = k_B = 1$，并采用几何化单位，其中所有长度均归一化到黑洞质量 $M$；因此，图中出现的参数为无量纲比值，如 $\rho_0 M^2$、$r_0/M$、$r/M$、$\omega M$、$T_H M$ 等。

章节片段

度规构造与视界结构

位于星系中心的天体物理黑洞被暗物质晕环绕，其引力印记改变了真空 Schwarzschild 几何 [24]、[50]、[51]。为了解析地对此环境建模，将一个中心黑洞嵌入到物质分布中，其密度轮廓 $\rho(r)$ 通过包络质量函数进入度规。

在文献中研究的密度轮廓中——尖峰型Navarro-Frenk-White（NFW）形式ρ∝r?1(1+r/rs)?2 [5]，[6]，以及带核的伪等温球——
零测地线：阴影与弱透镜效应
时空的因果结构可通过研究试探粒子与光子的测地线运动来探测[39]，[60]。这一分析可获取光子球（PS）、黑洞阴影边界以及弱场偏折角——这三个量编码了Plummer暗物质晕与弦宇宙学（CoS）对光子传播的联合引力印记。将测地线方程(ds/dλ)2=?（有质量粒子?=?1，光子?=0）应用于度规(10)，得到?A(r)t?2+r?2/A(r)+r2θ?2+r2sin2θφ?2=?，
类时测地线与ISCO分析
围绕黑洞的有质量试探粒子运动决定了吸积盘结构、能量提取效率以及X射线双星中观测到的准周期振荡频率[71]，[72]。在此背景下，一个中心量是ISCO（最内稳定圆轨道），它标记了存在稳定圆轨道的最内半径。在ISCO以内，不存在束缚圆运动，下落物质将直接坠入事件视界（EH）。对于Schwarzschild黑洞，ISCO位于r_ISCO=6M；
标量扰动、灰体因子与准正规范畴
黑洞时空的线性稳定性及其霍金辐射的谱性质，均由扰动场所经历的有效势控制。在本节中，我们推导了无质量标量场在Plummer-CoS度规(10)中传播的Schr?dinger型径向方程，考察了所得势垒的形状，使用Boonserm-Visser方法[73]，[74]计算了灰体因子的下界，并用切比雪夫谱方法确定了准正规模谱。
热力学
黑洞的热力学行为编码在一组小的状态函数中——质量M、霍金温度T_H、熵S、热容C_H和Gibbs自由能G——均从视界处的度规函数导出[27]，[29]。对于Plummer-CoS黑洞，附加参数ρ0、r0和α通过其对视界方程(31)的影响进入所有这些量。本节中我们求出M、T_H、S、C_H和G作为视界半径r_h的函数的显式形式，以及
结论
在本文中，我们研究了一个嵌入带核Plummer暗物质晕中并由弦宇宙学（CoS）包围的静态球对称黑洞，提供了一个统一框架，使星系尺度物质分布与拓扑缺陷的引力特征可以在同等条件下进行比较。度规函数A(r)对所有允许的参数值α∈[0,1)都存在唯一且明确定义的事件视界，视界半径随ρ0和α单调增长，而α≥1的情形对应于
CRediT作者贡献声明
Ahmad Al-Badawi：撰写–审阅与编辑，撰写–原稿，可视化，监督，资源，方法论，研究，形式分析，数据管理，构思。Faizuddin Ahmed：撰写–审阅与编辑，撰写–原稿，可视化，验证，监督，软件，资源，方法论，研究，资金获取，形式分析，数据管理，构思。?zzet Sakall?：撰写–审阅与编辑，撰写–原稿，
利益冲突声明
作者声明，他们没有已知的竞争性经济利益或个人关系可能影响本文所报告的工作。
致谢
作者感谢主编、处理编辑和匿名审稿人对稿件的仔细阅读和富有建设性的意见，这些意见改进了本工作的呈现。F.A.感谢印度浦那大学间天体与天体物理学中心（IUCAA）在他访问期间给予的亲切接待，部分工作在该访问期间完成。?.S.感谢TüB?TAK、ANKOS和SCOAP3的支持，以及COST行动CA22113、CA21106、CA23130的支持。
Ahmad Al-Badawi|Faizuddin Ahmed|?zzet Sakall?