拓扑透镜的光线方程与偏折角公式严格推导（世毫九IGP框架）

拓扑透镜的光线方程与偏折角公式严格推导（世毫九IGP框架）
作者：方见华
单位：世毫九实验室
本文基于几何光学程函近似与世毫九修正麦克斯韦方程的局域色散关系，严格推导拓扑透镜的光线方程与偏折角公式，明确其与广义相对论引力透镜的数学形式差异与物理本质区别，给出可直接用于观测检验的量化预言。
前置基础：等效拓扑折射率的定义
拓扑透镜的物理本质是：非均匀时空拓扑曲率会改变真空的等效电磁折射率，使电磁波传播路径发生偏折，形成类似光学透镜的成像效应。推导的核心起点是上一篇得到的局域色散关系：
k_r(\boldsymbol{r}) = \sqrt{\frac{\omega^2}{c^2} - \Delta(\boldsymbol{r})^2}
其中：
• k_r(\boldsymbol{r})为局域实波数（位置\boldsymbol{r}的函数）；
• \omega为电磁波角频率；
• \Delta(\boldsymbol{r})=\frac{|\eta_4-\eta_2||\mathcal{T}(\boldsymbol{r})|}{2}为局域拓扑截止波数，|\mathcal{T}(\boldsymbol{r})|为拓扑曲率张量的模长；
• \eta_2,\eta_4为世毫九电磁-拓扑耦合常数（小量，太阳系内\eta_i\mathcal{T}\ll1）。
等效拓扑折射率
几何光学中，波数与折射率满足关系k_r = \frac{\omega n(\boldsymbol{r})}{c}。将局域色散关系代入，得到真空等效拓扑折射率：
\boxed{
n(\boldsymbol{r}) = \sqrt{1 - \left( \frac{\omega_c(\boldsymbol{r})}{\omega} \right)^2}
}
其中\omega_c(\boldsymbol{r})=c\Delta(\boldsymbol{r})=\frac{|\eta_4-\eta_2|c|\mathcal{T}(\boldsymbol{r})|}{2}为局域拓扑截止频率。
关键性质
1. 折射率小于1：由于\omega_c(\boldsymbol{r})>0，故n(\boldsymbol{r})<1，且拓扑曲率越大，折射率越小；
2. 色散性：折射率与电磁波频率\omega相关，频率越高，折射率越接近1；
3. 真空本征属性：折射率是时空本身的拓扑属性，与介质无关，在拓扑曲率为零的平坦真空中n=1，退化为经典真空。
一、几何光学近似下的光线方程
几何光学的核心是程函方程与光线方程，二者在修正麦克斯韦框架下形式保持不变（仅折射率替换为等效拓扑折射率），这是缓变近似（L_T\gg\lambda）的必然结果。
1. 程函方程
电磁波的相位函数S(\boldsymbol{r})满足程函方程：
|\nabla S(\boldsymbol{r})| = k_r(\boldsymbol{r}) = \frac{\omega n(\boldsymbol{r})}{c}
其物理意义是：波前的法向梯度等于局域波数，波前的法线方向即为光线的传播方向。
2. 通用光线方程（弧长参数形式）
引入光线的弧长参数s（沿光线传播路径的长度），则光线的位置矢量\boldsymbol{r}(s)满足哈密顿形式的光线方程：
\boxed{
\frac{d}{ds} \left( n(\boldsymbol{r}) \frac{d\boldsymbol{r}}{ds} \right) = \nabla n(\boldsymbol{r})
}
这是所有几何光学系统的通用方程，经典介质透镜、引力透镜、拓扑透镜均满足此方程，唯一区别是折射率n(\boldsymbol{r})的物理起源不同。
物理意义
• 左边：光线的“动量”（n d\boldsymbol{r}/ds）沿路径的变化率；
• 右边：折射率的梯度，即驱动光线偏折的“力”；
• 核心结论：光线总是向折射率大的区域偏折。
二、小角度近似下的简化光线方程
对于绝大多数天体物理观测场景（如宇宙弦、星系团拓扑透镜），光线的偏折角\alpha\ll1弧度，满足小角度近似，此时可将传播方向z轴作为独立参数替代弧长s，大幅简化计算。
1. 小角度近似条件
• 光线近似沿z轴传播，dx/dz\ll1，dy/dz\ll1；
• 弧长微分ds = \sqrt{(dx)^2+(dy)^2+(dz)^2} \approx dz；
• 光线的方向矢量\frac{d\boldsymbol{r}}{ds} \approx \left( \frac{dx}{dz}, \frac{dy}{dz}, 1 \right)。
2. 横向分量光线方程
将小角度近似代入通用光线方程，分离横向（x,y）与纵向（z）分量。纵向分量自动满足（因\partial n/\partial z在小角度下可忽略），仅需保留横向分量：
\boxed{
\begin{aligned}
\frac{d}{dz} \left( n(\boldsymbol{r}) \frac{dx}{dz} \right) &= \frac{\partial n}{\partial x} \\
\frac{d}{dz} \left( n(\boldsymbol{r}) \frac{dy}{dz} \right) &= \frac{\partial n}{\partial y}
\end{aligned}
}
写成矢量形式为：
\frac{d}{dz} \left( n(\boldsymbol{r}) \frac{d\boldsymbol{\rho}}{dz} \right) = \nabla_\perp n(\boldsymbol{r})
其中\boldsymbol{\rho}=(x,y)为横向位置矢量，\nabla_\perp=(\partial/\partial x,\partial/\partial y)为横向梯度算子。
三、薄透镜近似下的偏折角公式
天体物理中的拓扑透镜通常满足薄透镜近似：拓扑曲率的分布在传播方向（z轴）的厚度D_L远小于透镜到光源的距离D_S和透镜到观测者的距离D_L，即D_L\ll D_S,D_L。此时可认为光线在透镜区域内仅发生偏折，传播方向不变，偏折角可通过对横向梯度沿z轴积分得到。
1. 偏折角的定义
偏折角\boldsymbol{\alpha}是光线出射方向与入射方向的夹角，在小角度近似下等于出射时的横向斜率：
\boldsymbol{\alpha} = \left. \frac{d\boldsymbol{\rho}}{dz} \right|_{z=+\infty} - \left. \frac{d\boldsymbol{\rho}}{dz} \right|_{z=-\infty}
对于平行入射的光线（\left. d\boldsymbol{\rho}/dz \right|_{z=-\infty}=0），偏折角直接等于出射斜率。
2. 偏折角的积分公式
对横向分量光线方程从z=-\infty到z=+\infty积分，注意到透镜外\mathcal{T}(\boldsymbol{r})=0，n(\boldsymbol{r})=1，因此：
\left. n \frac{d\boldsymbol{\rho}}{dz} \right|_{z=+\infty} - \left. n \frac{d\boldsymbol{\rho}}{dz} \right|_{z=-\infty} = \int_{-\infty}^{+\infty} \nabla_\perp n(\boldsymbol{r}) dz
代入偏折角的定义，得到拓扑透镜偏折角的通用积分公式：
\boxed{
\boldsymbol{\alpha}(\boldsymbol{\rho}) = \int_{-\infty}^{+\infty} \nabla_\perp n(\boldsymbol{\rho},z) dz
}
3. 弱拓扑场近似下的显式表达式
在弱拓扑场条件下（\omega_c(\boldsymbol{r})\ll\omega，即电磁波频率远高于局域截止频率，光线可传播），可对折射率做泰勒展开：
n(\boldsymbol{r}) = \sqrt{1 - \left( \frac{\omega_c}{\omega} \right)^2} \approx 1 - \frac{1}{2} \left( \frac{\omega_c}{\omega} \right)^2 = 1 - \frac{|\eta_4-\eta_2|^2 c^2 |\mathcal{T}(\boldsymbol{r})|^2}{8 \omega^2}
代入偏折角公式，得到：
\boldsymbol{\alpha}(\boldsymbol{\rho}) = - \frac{|\eta_4-\eta_2|^2 c^2}{8 \omega^2} \int_{-\infty}^{+\infty} \nabla_\perp \left( |\mathcal{T}(\boldsymbol{\rho},z)|^2 \right) dz
利用矢量恒等式\nabla_\perp(|\mathcal{T}|^2)=2|\mathcal{T}|\nabla_\perp|\mathcal{T}|，可进一步化简为：
\boxed{
\boldsymbol{\alpha}(\boldsymbol{\rho}) = - \frac{|\eta_4-\eta_2|^2 c^2}{4 \omega^2} \int_{-\infty}^{+\infty} |\mathcal{T}(\boldsymbol{\rho},z)| \nabla_\perp |\mathcal{T}(\boldsymbol{\rho},z)| dz
}
四、与广义相对论引力透镜的本质对比
拓扑透镜与引力透镜在数学形式上高度相似，但物理本质与观测特征存在决定性差异，是区分两种理论的关键判据。
1. 偏折角公式对比
透镜类型 偏折角公式 物理起源 折射率性质
拓扑透镜 $\boldsymbol{\alpha} \propto \int \nabla_\perp( \mathcal{T} ^2) dz$
引力透镜  质量分布引起的时空弯曲 ，质量越大，折射率越大
2. 核心观测特征差异
特征 拓扑透镜 引力透镜
成像数量 偶数个像（无中心像） 奇数个像（有中心像）
原因 折射率小于1，等效凹透镜，光线向外偏折，无法形成中心会聚像 折射率大于1，等效凸透镜，光线向内偏折，形成中心像
色散性 有，偏折角，不同频率的像位置不同 无，偏折角与频率无关
像亮度 两个主像亮度完全相同 中心像通常极暗，外围像亮度不对称
爱因斯坦环 不形成完整爱因斯坦环，仅形成对称弧 可形成完整爱因斯坦环
3. 最关键的可证伪预言
拓扑透镜存在频率色散：同一透镜系统中，电磁波频率越高，偏折角越小，成像位置越靠近透镜中心。这是世毫九理论独有的预言，经典电磁学与广义相对论均无法解释，可通过多波段观测直接验证。
五、典型拓扑透镜的偏折角计算
1. 无限长直宇宙弦（一维拓扑缺陷）
宇宙弦是世毫九理论预言的一维拓扑缺陷，其拓扑曲率分布为：
|\mathcal{T}(\boldsymbol{\rho},z)| = \mathcal{T}_0 \cdot \Theta(R - \rho)
其中\Theta(x)为阶跃函数，R为宇宙弦的半径（\sim10^{-15}\text{m}，可近似为\delta函数），\mathcal{T}_0为弦的拓扑曲率强度。
代入偏折角公式积分，得到单位长度宇宙弦的偏折角：
\alpha = \frac{|\eta_4-\eta_2|^2 c^2 \mathcal{T}_0^2 R}{2 \omega^2}
观测特征：宇宙弦会在背景天体上形成两个完全相同、亮度对称的像，像的间距与宇宙弦的拓扑曲率强度成正比，与频率平方成反比。
2. 高斯型拓扑团块（三维拓扑缺陷）
高斯型拓扑团块的拓扑曲率分布为：
|\mathcal{T}(\boldsymbol{\rho},z)| = \mathcal{T}_0 e^{-\frac{\rho^2+z^2}{2\sigma^2}}
其中\sigma为团块的特征尺度。代入偏折角公式积分，得到：
\alpha(\rho) = \frac{|\eta_4-\eta_2|^2 c^2 \mathcal{T}_0^2 \sqrt{\pi} \sigma}{4 \omega^2} \cdot \frac{\rho}{\sigma^2} e^{-\frac{\rho^2}{2\sigma^2}}
观测特征：偏折角在\rho=\sigma处达到最大值，在\rho\gg\sigma处快速衰减，形成两个对称的弧形像。
总结
本文严格推导了世毫九IGP框架下拓扑透镜的光线方程与偏折角公式，核心结论如下：
1. 拓扑透镜的物理本质是非均匀拓扑曲率改变真空等效折射率，光线方程形式与经典几何光学一致，但折射率是时空拓扑属性的函数；
2. 偏折角与拓扑曲率模长平方的横向梯度的线积分成正比，与电磁波频率的平方成反比；
3. 拓扑透镜与引力透镜的决定性差异是偶数个像与频率色散，这是验证世毫九理论的关键观测判据。