变化磁场与引力：从麦克斯韦方程组到爱因斯坦场方程的推导

把磁场和引力放在一起推公式第一眼看上去像是什么“民科”选题。但这篇文章不推翻任何现有物理也不打算抛出“旋转磁场就能反重力”这种结论。我想做的事其实很朴素从麦克斯韦方程组和爱因斯坦场方程出发用第一性原理的方式走一遍“变化磁场如何影响时空”的推导全程看看能推出什么、能验证什么以及这件事在统一场论的大背景里到底意味着什么。核心会落在三个地方——变化磁场产生的感生电场方程、电磁场作为引力源时的度规扰动表达式以及这种“电磁-引力耦合”实际效应到底有多强、能不能被实验看到。适合有一定电动力学或广义相对论基础、但还没系统算过这类问题的读者零基础也没关系公式我会用白话拆开讲。1. 问题拆解变化磁场到底能“产生”什么1.1 先厘清一个常见误解第一次看到“变化磁场产生引力场”这种说法绝大多数人的第一反应是那是不是给电磁铁通上交流电周围物体就会变轻这个画面很吸引人但它和物理现实差得很远。先看标题里的另外一半——“变化磁场产生电场方程”这个没有任何争议就是法拉第电磁感应定律。变化的磁场会在周围空间激发涡旋电场这是发电机、变压器、无线充电所有技术的基础方程形式极其干净∇×E -∂B/∂t。磁场随时间一变电场立刻出现这是电磁学内部的事。但“变化磁场产生引力场”就完全不是一个层面的事了。引力在广义相对论里的来源不是“磁场”这个物理量本身而是能量动量张量。任何有能量的东西——静质量、动能、电磁场能量、压力、应力——都会弯曲时空。所以“变化磁场产生引力效应”的正确问题是变化磁场携带了能量动量这些能量动量会如何作为引力源去弯曲周围的时空这个区别非常重要。变化磁场“产生”电场是电磁场内部的动力学过程变化磁场“产生”引力效应是电磁场作为物质源去影响几何的过程。一个是方程左边和右边同属麦克斯韦系统另一个是物质场进入爱因斯坦方程。机制完全不同不能拿“磁生电”的直觉去套“磁生引力”。1.2 第一性原理推导的技术路线所谓第一性原理推导在我的理解里就是用最少的假设、从最基本的方程出发把结论一步一步算出来而不是靠类比或直觉。针对这个题目技术路线可以拆成四步。第一步写出变化磁场的电磁场构型。最常见也最容易处理的是均匀磁场随时间简谐变化比如B(t) B₀cos(ωt)或者更一般的缓变B(t)。这一步同时要把感生电场算出来因为一旦磁场变了空间里就会有涡旋电场电磁场的能量密度是电场和磁场一起贡献的。第二步计算这个电磁场的能量动量张量T^{μν}。这是连接电磁学和引力学的桥。T^{00}是能量密度T^{0i}是动量密度T^{ij}是应力张量。对纯磁场的特殊情况T^{00} B²/(2μ₀)就是磁场能量密度。第三步把T^{μν}代入线性化的爱因斯坦场方程或者说弱场近似下的引力场方程求出度规扰动h_{μν}。这一步得到的就是“磁场产生的引力场”的定量描述。电磁场的能量密度在引力方程里扮演的角色等价于一个“等效质量密度”ρ_eff u/c²。第四步做数量级估算。这一步很多理论推导会跳过但恰恰是最有价值的。算出1特斯拉的磁场、10特斯拉的MRI磁体、乃至10¹¹特斯拉的磁星它们产生的额外引力到底有多强和现有的测量精度差多少。只有把这件事落到数字上才能判断这个效应在什么尺度下有意义。整个路线不引入任何新物理假设用的全是教科书里成熟的内容。这样做的好处是结论可靠不会滑向玄学坏处是听起来不够“颠覆”但这本来就是物理学该有的样子。2. 基础工具麦克斯韦方程组、能动张量与弱场近似2.1 变化磁场与感生电场法拉第定律的微分形式从麦克斯韦方程组说起。在SI单位制下描写变化磁场和电场关系的方程是∇×E -∂B/∂t这是法拉第电磁感应定律的微分形式也是标题里“变化磁场产生电场方程”的直接来源。它说的是磁场如果随时间变化它周围会卷起一个涡旋电场电场线的环量正比于磁通量的变化率。对均匀轴向磁场B(t) B₀cos(ωt)感生电场是围绕对称轴的涡旋场。取一个半径为r的圆环用法拉第积分形式∮E·dl -d/dt ∫B·dS可以得到E_φ(r,t) -(r/2)dB/dt (ω r B₀/2)sin(ωt)这个式子很有意思。感生电场的幅值正比于频率ω、正比于离轴距离r而且和磁场存在90°相位差。当磁场达到峰值时电场刚好为零磁场过零的时候电场最大。后面算能量密度的时间平均时这个相位关系会直接决定交叉项怎么抵消、能量密度怎么分配。2.2 电磁场凭什么影响引力场能量动量张量引力在广义相对论里不和“磁场”直接作用它和作用在时空上的能量动量分布直接耦合。描述这种分布的数学对象就是能量动量张量T^{μν}它是一个4×4的对称矩阵每个分量的物理含义不同。电磁场的能量动量张量可以从拉格朗日密度L -¼F_{μν}F^{μν}出发用变分法推出T^{μν} F^{μα}F^{ν}{}{α} - ¼η^{μν}F{αβ}F^{αβ}这个式子看着抽象但拆开看非常直观。F_{μν}是电磁场张量它的分量里面装的就是电场E和磁场B。T^{00}是能量密度T^{0i}是沿i方向的动量密度坡印廷动量T^{ij}是电磁应力张量。在只有磁场、没有电场的简单情况下能量密度为T^{00} B²/(2μ₀)磁场B1特斯拉时能量密度大约为4×10⁵ J/m³。换算一下1立方米的1特斯拉磁场储能相当于40克汽油完全燃烧放热。这其实不小了但放到引力效应里就非常小因为引力耦合常数G/c⁴只有8.26×10⁻⁴⁵ s²/(kg·m)是一个小到反直觉的数字。值得强调的是T^{ij}的各向异性应力同样会弯曲时空不只是能量密度。磁场沿着某一方向会产生横向的磁压力和纵向的磁张力这些项在精确解里也会进入度规扰动。不过对数量级估算来说能量密度项通常起主导作用后面我在特例里会统一处理。2.3 弱场近似下的爱因斯坦方程一个“源”与“响应”的线性系统完整爱因斯坦场方程是非线性的直接拿来解这个题目会很吃力。好在日常生活和绝大多数实验室场景里引力都弱得可怜可以把时空度规写成平直度规加一个小扰动g_{μν} η_{μν} h_{μν}|h_{μν}| 1把爱因斯坦场方程线性化取谐和规范后可以得到一个形式上非常优美的方程□h̄_{μν} -(16πG/c⁴)T_{μν}这里h̄_{μν}是迹反演后的度规扰动□是达朗贝尔算符。看一眼这个式子的结构左边是度规扰动的动力学右边是物质的能量动量完全是一个线性响应系统。物质告诉时空怎么弯曲时空告诉物质怎么运动——在线性化下前半句被这个方程精确表达。再做一个牛顿近似场缓变忽略时间二阶导数物质运动速度远小于光速00分量会退化成经典的泊松方程∇²Φ 4πGρ_eff其中ρ_eff T_{00}/c²。这一步其实就说明了所有问题的核心在弱场低速条件下任何形式的能量密度都可以当成“等效质量密度”去产生牛顿引力势。磁场能量密度也不例外。所以要算“变化磁场产生的引力场”本质上就是在解一个泊松方程源项是磁场能量密度除以c²。这个路径极其清晰也是整篇推导里最容易被忽略的简化策略——很多人一想到引力就头大但一旦回到牛顿极限一切就都变成高中物理里的万有引力图景了。3. 核心推导从磁场到度规扰动的完整过程3.1 电磁场能量动量张量的显式展开先明确符号约定。采用平直时空度规η_{μν} diag(1, -1, -1, -1)电场E和磁场B的分量关系为F_{0i} -E_i/cF_{ij} ε_{ijk}B_k。在这个约定下电磁场能量动量张量的各分量展开为T^{00} ½(ε₀E² B²/μ₀)T^{0i} (1/μ₀c)(E × B)_iT^{ij} -ε₀E_iE_j - (1/μ₀)B_iB_j ½δ_{ij}(ε₀E² B²/μ₀)对纯磁场情况T^{00} B²/(2μ₀)。如果磁场还随时间变化法拉第定律保证空间里同时会有感生电场所以严格讲E不能忽略。但感生电场的能量密度相比磁场本身小多少取决于变化频率和系统尺度的关系这个我在3.2节专门分析。这里必须提一下T^{ij}的表达式中的“磁压”和“磁张力”。沿磁场方向看应力是“拉”的负压强垂直磁场方向看是“推”的正压强。这种各向异性意味着磁场不是一个简单的理想流体它的引力效应会比同能量的尘埃更复杂。在中子星或磁星的建模里这种各向异性应力会直接影响星体的形状和时空几何。3.2 简谐变化磁场的特例能量密度、等效质量密度与引力势考虑一个半径为R的柱形区域里面是均匀轴向磁场B(t) B₀cos(ωt)。这个模型虽然不是最复杂的但足够说明所有关键性质。根据3.1节电磁场能量密度为u(t) ½(ε₀E² B²/μ₀)其中感生电场E_φ (ω r B₀/2)sin(ωt)。代入后u(t) B₀²/(2μ₀)cos²(ωt) ε₀ω²r²B₀²/8·sin²(ωt)时间平均后⟨u⟩ B₀²/(4μ₀) ε₀ω²r²B₀²/16第一项是磁场能量的时间平均第二项是感生电场能量的时间平均。两个的比值是⟨u_E⟩/⟨u_B⟩ (ε₀ω²r²B₀²/16) / (B₀²/(4μ₀)) μ₀ε₀ω²r²/4 (ωr/c)²/4这个比值解释了什么时候能忽略电场能量。如果系统尺度r远小于电磁波的约化波长c/ω也就是变化足够缓慢、系统足够小那么电场能量远远小于磁场能量可以安全忽略。对实验室里一个1米尺度、频率在千赫兹到兆赫兹的磁场这个比值大约在10⁻¹⁰到10⁻⁶量级完全可忽略。但如果是光频段的电磁场电场和磁场能量就均分了。忽略电场能量后时间平均能量密度为ρ_eff ⟨u⟩/c² B₀²/(4μ₀c²)代入数值B₀ 1T时ρ_eff ≈ 2.2×10⁻¹² kg/m³。注意这是时间平均值如果是瞬时值系数会从1/4变回1/2。这个等效质量密度非常低一立方米的1特斯拉磁场“重量”只有约2×10⁻¹²千克比一粒灰尘轻好几个数量级。现在解泊松方程。为了估算方便把磁场区域近似看成一个半径为R的均匀密度球。球内外的引力势为Φ(r) (2π/3)Gρ_eff(r² - 3R²) (r ≤ R)Φ(r) -(4π/3)Gρ_effR³/r (r R)对应的额外引力加速度球内g -(4π/3)Gρ_eff·r方向指向球心球外g -GM_eff/r²M_eff (4π/3)R³ρ_eff。这个结果意味着磁场区域确实会像一团普通物质一样产生一个微弱的吸引性引力场。3.3 附加引力的数量级估算为什么实验这么难做把数字带进去感受一下现实。B₀ 1T的螺线管半径R 1m。ρ_eff ≈ 2.2×10⁻¹² kg/m³总有效质量M_eff ≈ (4π/3)×2.2×10⁻¹² ≈ 9.2×10⁻¹² kg。它在表面产生的附加引力加速度g_extra GM_eff/R² 6.67×10⁻¹¹×9.2×10⁻¹² ≈ 6.1×10⁻²² m/s²地球表面重力加速度是9.8 m/s²二者相差约10²²倍。这是什么概念即便把磁场拉到10T这是医用MRI磁体的两倍还多g_extra也只提升100倍到10⁻²⁰量级。目前最强的稳态磁场约45T实验室脉冲磁场能做到百T量级但即便100Tg_extra也只有10⁻¹⁸到10⁻¹⁷ m/s²。要观测这种信号需要加速度计灵敏度达到10⁻¹⁸g量级。现代超导重力仪的灵敏度大约在10⁻¹²g量级冷原子重力仪在低频段能做到10⁻¹¹到10⁻¹³g量级。也就是说差了至少5到7个数量级。想用传统实验来验证这个效应得先颠覆现有测量技术。那有没有办法把有效质量密度做大磁场能量密度正比于B²所以唯一的路是超级磁场。迄今宇宙中已知最强的宏观磁场来自磁星一种中子星的特例表面磁场可达B ≈ 10¹¹T。此时能量密度u ≈ (10¹¹)²/(2×4π×10⁻⁷) ≈ 4×10²⁷ J/m³等效质量密度ρ_eff ≈ 4.4×10¹⁰ kg/m³是水的密度的一千多万倍。磁星半径约10km磁场储能的等效质量约为M_eff ≈ 4.4×10¹⁰ × (4π/3)(10⁴)³ ≈ 1.8×10²³ kg这大约是地球质量的3%。所以磁星这种极端天体磁场对自身引力的贡献已经不可忽略了它会额外加深引力势、增大引力红移、影响星体平衡结构。这也是天体物理里“磁化中子星”模型需要用广义相对论而非牛顿理论严格处理的原因之一。数字告诉我们题目里的效应不是不存在而是太微小只有在宇宙极端尺度下才可能留下可观测痕迹。3.4 时变部分的物理边界振荡牛顿势不等于可探测引力波从3.2节可以看到由于B(t) ~ cos(ωt)能量密度里会出现cos²(ωt) ½(1cos(2ωt))的项。所以磁场产生的引力势并不是静态的它包含一个直流分量和一个以2ω振荡的分量。这个振荡分量在近场区域就是一个随时间变化的牛顿势会以2ω的频率周期性地调制周围粒子的运动。从引力本质讲这意味着磁场区域的等效质量在周期性地“脉动”。但这里有一个极其常见、又极其隐蔽的误区不能把“振荡的牛顿势”直接当成“引力波”。引力波的产生需要时变的四极矩或更高阶多极矩质量单极不辐射、偶极也不辐射动量守恒禁止。对于均匀轴对称的磁场能量密度的振荡只会改变局部的质量密度不会改变质量分布的四极矩所以不会产生引力波辐射。换句话说如果拿一个均匀螺线管通上交流电周围会感受到一个同频(确切说2ω)的很弱的引潮力调制但不会有引力波以这种频率辐射出去。想通过磁场变化产生引力波需要让磁场本身具有时变的高阶多极矩比如让两个强磁体高速旋转并保持质量分布不对称这种源产生的引力波辐射强度依然弱到现有探测器无法触及。区分“近场振荡”和“波辐射”这两个概念是这个题目里最值得记下的物理边界之一。4. 验证途径与现实边界4.1 实验室验证精密力测量与噪声困境门槛已经算清楚了实验室场景里这个效应被淹没在噪声之下的概率几乎是100%。但依然有人试过。历史上围绕“强磁场是否产生异常引力效应”做过不少实验比如旋转超导体的引力屏蔽实验、强磁场下的扭秤实验。这些实验的报告里偶尔出现“异常信号”但最终都在更精密的重复测量中被排除或者被归因于热效应、电磁干扰等系统误差。我见过很多初学者对这类实验抱有极大热情但我的态度一直是先做估算再谈实验。用上一节的数量级分析实验室尺度的强磁场产生的附加引力加速度在10⁻²⁰ m/s²量级而最好的重力仪灵敏度在10⁻¹³ m/s²左右中间隔着七个数量级。除非出现某一项原理性的测量技术突破否则这个通道没有现实可测性。这不是说实验室路线全无意义。如果把目标从“直接测引力加速度”改成“测度规扰动的相位效应”用原子干涉仪或超导谐振腔测量磁场区域内外的光速变化、时间延迟差异理论上更接近可测范围。但即便如此现有的干涉测量精度仍然离需求差几个数量级。这个方向可以作为长远的仪器科学挑战不宜抱短期的验证期待。4.2 天体物理场景磁星与宇宙磁场实验室测不到不代表自然界没有体现。上一节算过磁星磁场能量等效质量可以达到地球质量的3%。这意味着磁星的引力场和普通相同质量中子星的引力场并不完全相同。这种差异会带来两个可直接检验的效应。第一个是引力质量与转动惯量的关系变化。磁星因为储藏着巨大磁场能量其有效引力质量比同质量的非磁化星略大同时磁场应力各向异性会让星体的形变产生“磁椭率”。这种椭率会影响星体自转的岁差运动也能通过脉冲星计时观测加以限制。当前对磁星的观测显示磁场的引力贡献虽然没有被精确分离但广义相对论框架下的磁化中子星模型与观测整体相符。第二个是磁场能量密度的引力透镜效应。虽然在磁星周围这个效应依然很弱但相比实验室差了好几十个数量级至少是理论上有希望的方向。更普遍地说任何具有宇宙磁场的结构——星系团、活动星系核——磁场能量密度虽然远低于物质密度局部区域却可能达到与动力学平衡不可忽略的水平。在星系团尺度微高斯磁场在十万秒差距尺度上储存的能量对应的等效质量能与星系团内热气体的部分能量相比虽然通常不到百分之几但在冷团核心区域已经进入模型必须考虑的范围。这类天体物理场景本质上是“电磁-引力耦合”最现实的检验平台。4.3 统一性意义从爱因斯坦-麦克斯韦系统到卡鲁扎-克莱因回到标题里最大词——“统一性意义”。这个题目推导出来的电磁场引力耦合在物理学史上到底处于什么位置答案要从广义相对论诞生之前说起。在爱因斯坦方程的右边放上电磁场能量动量张量得到的就是爱因斯坦-麦克斯韦系统。这个系统在经典层面是完全自洽的电磁场弯曲时空弯曲时空反过来影响电磁场传播。也就是说“电磁和引力统一在经典框架下”这件事在1915年广义相对论完成以后就已经做到了——它不叫统一叫“所有物质场都是时空几何的源”。但物理学家并不满足于这种“被动统一”他们想要的是更根本的统一让电磁力本身变成时空几何的一部分。这个方向最著名的尝试是卡鲁扎-克莱因理论。1919年卡鲁扎提出在五维时空里写纯引力方程五维度规的某些分量在四维视角下恰好表现为一个矢量场而这个矢量场满足的方程正是麦克斯韦方程组。换句话说电磁场被“解释”成了第五维的几何效应。这个理论非常优雅但代价是第五维为什么卷曲得那么小、为什么看不见——这些问题直到现代弦论才给出新的回答。今天我们把电磁场理解为U(1)规范场把引力理解为自旋2的无质量场二者的统一需要在量子引力的框架内完成。标题推导的这个“经典耦合”只是这幅大图景里的第一块积木而不是终点。4.4 主流物理的边界经典耦合成功不等于量子统一完成一定要说清楚边界这篇文章从头到尾处理的都是经典场论电和磁、引力和时空都是经典对象。在这个层面上理论自洽、可以推导、可以估算。但物理学家真正头痛的统一是量子层面。电磁场是量子化的光子是U(1)规范玻色子引力场理论上也应该量子化对应自旋2的引力子。但把引力量子化和物质场放在同一个框架里会碰到紫外发散等问题。弦论、圈量子引力、渐近安全引力等方案都试图解决这个问题但目前没有一个得到决定性实验验证。电磁和引力的“经典动力学耦合”已经够用了但“量子统一”仍然悬而未决。理解这两者的区别比记住任何公式都重要。5. 推导中的常见坑与实操心得5.1 符号约定和单位制第一次算最容易被G/c⁴坑我刚开始做这类计算时最喜欢在G/c⁴这个因子上翻车。G ≈ 6.67×10⁻¹¹ N·m²/kg²c² ≈ 9×10¹⁶ m²/s²所以G/c⁴ ≈ 8.26×10⁻⁴⁵ s²/(kg·m)。这个数字小到不真实以至于算完总觉得是不是把π放错了位置。建议的做法是固定一套单位制再动笔。做这类半经典推导时我习惯用几何单位制(G c 1)把星号写在草稿最上方所有源项直接以能量形式进入方程最后需要代入实验数值时再统一还原成SI。这么做能避免大量不必要的因子混乱。另一个容易翻车的地方是度规符号约定。平直度规用(,-,-,-)还是(-,,,)会导致F_{0i}、T^{00}的符号定义不同E和B的对应关系差个负号。不要在两套约定之间反复横跳选定一套写到论文结束。我自己常用(,-,-,-)因为电磁学里的常规公式在那个约定下最自然。5.2 “产生”一词的使用边界物理学里用词必须准确。标题里“变化磁场产生引力场”这个说法和“变化磁场产生电场”严格讲不是同一种“产生”。法拉第定律里的“产生”是同层次场之间的动力学激发电场和磁场都是电磁场的一部分麦克斯韦方程组把它们耦合在一起。而引力场方程中的“产生”其实是“贡献”——变化磁场通过它的能量动量张量成为爱因斯坦方程右边的一个源项。它并没有像法拉第定律那样在磁场之外“拽”出一个新的基本场。磁场产生的引力效应在原理上和一根烧红的铁棍产生的引力效应没有本质区别都是“能量弯曲时空”。想明白这一层就不会被标题骗去期待某种新力。也就能理解为什么统一场论的真正难点不在经典层面——经典层面早就统一进能量动量张量里了难的是量子层面的重构。5.3 碰到这类跨领域推导我的工作流程这类从电磁学跨到引力学的推导我建议按固定流程走既省时间又不容易出错。第一步写清楚出发点和假设。这篇的假设是弱场、缓变、非相对论运动。任何一个假设不成立后续公式都要换。第二步求源项。把T^{μν}完整写出来不要只写能量密度。第三步做量纲分析。用G/c⁴的量纲检查方程两边是否平衡这一步能解决一大半低错。第四步代特例、算数量级。不要代完公式就结束一定要估算这个效应在现实条件下有多强这才是一个物理推导真正接地气的部分。最后一步回到假设检查结论适用范围。如果发现结论依赖某些极端条件就说清楚条件是什么。这套流程救过我很多次。特别是“代完公式必须算数量级”这条能挡住很多貌似正确、实则无用的推导。写到最后说一点我的体会。这类题目真正让人着迷的地方不是它给出了某个炫酷的结论而是它把两套看起来关系不大的物理理论干干净净地编织在一起。磁场和引力一个是电磁学的日常一个是时空几何的日常它们之间的耦合强度被G/c⁴压得低到离谱但在磁星表面、在宇宙尺度、在引力波和强磁场并存的极端环境里这种耦合终究会露出身影。做推导的时候别急着追求“大新闻”老老实实把每一步因子算对、把每一个数量级估准就已经比大多数网上流传的电磁引力文章高出一截了。如果你也想动手推一遍我建议从B(t) B₀cos(ωt)这个最简单的模型开始一步步走完整个流程你会对“场论里的源”这个概念有一个非常直观的新理解。