1 引言
引力波的探测开启了研究黑洞(Black Holes, BHs)等致密天体物理现象的新窗口。黑洞并合最后阶段的引力场“铃宕”(ringdown)信号被认为包含了黑洞周围时空结构的“指纹”——准正则模(Quasinormal Modes, QNMs)，其谱线可用于推测黑洞的质量和角动量。理论上，这些QNMs被理解为在黑洞光子球(Photon Sphere, PS)区域附近被束缚的时空阻尼振荡。然而，光子在经典光子球轨道上是不稳定的，这引发了一个根本性问题：一个不稳定的区域如何能够局域化引力波模式？这需要对QNMs，特别是其空间分布，进行更精确的表征。
模拟引力(Analogue gravity)实验为在实验室可控环境中研究弯曲时空物理提供了一条有力途径。本文提出并实现了一种新颖的光学模拟方案：利用费马度规，将一个真实的4维(3+1D) Schwarzschild黑洞时空度规映射到一个2维(2+1D)曲面上，并严格保留了原时空中的零测地线（光线轨迹）。通过求解这个2D曲面上的波动方程（等价于弯曲时空中的克莱因-戈登方程），理论上可以模拟出黑洞的QNMs。基于此，研究团队设计并利用3D打印技术制备了掺染料的Schwarzschild曲面光学微腔，在实验上成功演示了在光子球处以及传统回音壁模式(Whispering Gallery Modes, WGMs)的激光激射。
2 理论预言
2.1 光子球及其稳定性
对于一个无电荷、无角动量的球对称Schwarzschild黑洞，其时空在赤道面上的费马度规（其空间部分）等同于一个嵌入三维欧氏空间的旋转曲面（称为Schwarzschild曲面）的诱导度规。这个曲面上存在一个唯一的圆形周期轨道，正好位于其腰部（即曲率半径最小的“细腰”处），对应于Schwarzschild黑洞的光子球半径rPS= 3GM/c2。通过雅可比方程和曲面的高斯曲率分析，可以证明在整个所考虑的外部区域（r> rh），高斯曲率均为负值，这意味着这个光子球轨道是经典不稳定的。这种不稳定性在图1c所示的庞加莱截面(Poincaré surface of section)中也得到了直观反映，光子球轨道（蓝色星号）位于一个鞍点上。
2.2 Schwarzschild曲面上的波动方程
在黑洞的铃宕阶段，其背景下的标量扰动演化由波动方程（无质量克莱因-戈登方程）描述。有趣的是，这个方程同样描述了限制在一个2D弯曲曲面（即一个(2+1)D时空）上的电磁波传播。这个等价性是建立引力波与弯曲曲面上光波之间类比的基础，使得在实验室中通过薄层弯曲波导来模拟黑洞QNMs成为可能。
2.3 光子球模的存在性
利用旋转对称性，可以将波动方程分离变量，得到一个关于径向坐标的类薛定谔方程。方程中的有效势Veff(r)在Schwarzschild曲面上呈现出一个“山丘”状，其峰值正好位于光子球rPS处，这与量子力学中的势垒类似。为了计算开放腔的本征模，研究人员将系统建模为一个“三明治”结构：中间的Schwarzschild曲面（区域II）连接着两个外侧的截头圆锥面（区域I和III，代表空气），并匹配界面处的边界条件。
通过运用WKB近似和Airy函数方法求解，研究人员发现，当频率足够高时（即有效势的峰值低于“能量”水平），有效势会在腔内形成两个“势阱”，分别位于两个边界附近，从而支持出类似于传统WGMs的模，其能量被全内反射约束在边界附近。而当频率低于势垒峰值时，有效势仅在光子球附近提供一个振荡区域，而两侧的边界附近为衰减区。在这种情况下，研究人员从行列式的零点中解析地求解出了一族新模式，其包络函数在光子球rPS处达到峰值，其空间强度分布也集中在光子球环附近，这被识别为基本的光子球模(Photon Sphere Modes, PSMs)。理论计算得到的PSMs和WGMs的径向强度分布与后续的数值模拟结果高度吻合。
2.4 准正则模的数值模拟
为了验证理论预言，研究团队进行了有限差分时域(Finite Difference Time Domain, FDTD)模拟。通过将Schwarzschild曲面保角变换到梯度折射率平面上进行仿真，研究人员获得了系统的共振谱。谱中的峰值对应于不同的准正则模，从它们的强度分布中可以清晰地识别出两个家族：强度峰值位于腰部（光子球）的PSM家族，和强度被约束在系统边界附近的WGM家族。提取的径向强度分布与理论预测高度一致，强有力地验证了理论。此外，模拟还给出了不同模式的质量因子(Quality Factor, Q)。由于WGMs主要通过边界的全内反射约束，其Q因子通常高于PSMs。然而，PSMs的Q因子与高阶WGMs相当甚至更高，这突显了曲率诱导的光约束在传播层中的有效性。
3 实验演示
3.1 光子球模的实验证据
研究团队利用双光子聚合3D打印技术，在掺有Pyrromethene 597染料的IP-G光刻胶中制备了Schwarzschild曲面微腔激光器。由于WGMs的阈值较低，在均匀泵浦下容易占主导地位。为了增强观测PSM激光的几率，研究人员采用空间光调制器将532 nm的泵浦光整形为一个窄条，并选择性地泵浦在微腔的腰部（光子球位置）。
在足够高的泵浦能量下，观察到了多模激光发射。通过对发射光谱进行Lomb–Scargle周期图分析，可以识别出与不同圆形轨道长度（即光学腔长）对应的谱峰。最强的峰值对应的轨道直径（44.3 μm）与微腔腰部的直径（44.5 μm）非常接近，这提供了PSM家族存在的第一个实验证据。相应的激光谱峰在光谱中形成了一个与方位角量子数m增加对应的等间距频率梳，与理论预言一致。在显微镜下，可以直接在腰部位置观察到两个激光光斑，直观地证实了光子球激光。当将泵浦条纹移动到边界附近以选择性地泵浦WGMs时，Lomb–Scargle周期图中只剩下对应边界轨道长度的谱峰，而腰部的激光发射则消失，这进一步区分了PSMs和WGMs。
3.2 光子球模的空间分布测量
为了精确表征一个特定PS激光模（波长为580.3 nm）的空间分布，研究人员将泵浦窄条沿微腔对称轴（z轴）扫描，并记录该波长处的峰值强度。在应用Lucy–Richardson反卷积算法后，获得了该PSM的空间强度分布，其被强烈地限制在腰部位置，有力地证明了其约束机制是由曲率诱导的势阱驱动，而非样品边界的反射。同样，对代表整个PSM家族的Lomb–Scargle谱峰幅度进行扫描测量，得到的平均空间轮廓与根据实验样品尺寸和平均方位角量子数m≈377计算的理论PSM基本模的预测表现出极好的一致性，在位置和空间延展上都吻合良好。
4 结论
总而言之，本文提出了一个维度约化的4D Schwarzschild黑洞光学模拟，它忠实地保留了其零测地线并共享相同的控制方程。通过理论分析、数值模拟和实验验证，首次在桌面条件下观测并证实了光学模拟Schwarzschild曲面上光子球模的存在和激射。这种由空间曲率（而非腔体边界）提供的全新模式约束机制，为研究引力物理中的黑洞现象（如铃宕谱、模式稳定性、黑洞光谱学等）提供了一个互补的光子学平台，并启发了基于非欧几里得曲面和势阱捕获原理的新型激光微腔的设计思路。
