查看“︁弯曲时空中的麦克斯韦方程组”︁的源代码

{{NoteTA|G1=物理学}}
[[Image:Gravitation space source.png|thumb|200px|引起时空弯曲]]
[[物理学]]中，'''弯曲时空中的麦克斯韦方程组'''（Maxwell's equations in curved spacetime）制约着弯曲[[时空]]（其间的[[度规]]可能不是[[闵可夫斯基时空|闵可夫斯基性]]的）中的[[电磁场]]的动力学。它们可以被认为是真空中的[[麦克斯韦方程组]]在[[广义相对论]]框架中的扩展，而真空中的麦克斯韦方程组只是一般化的麦克斯韦方程组在局部平直时空中的特殊形式。但由于在广义相对论中电磁场本身的存在也会引起时空的弯曲，因此真空中的麦克斯韦方程组应被理解为一种出于方便的近似形式。

然而，这种形式的麦克斯韦方程组仅仅对真空情形下的麦克斯韦方程组有用，这也被称作“微观”麦克斯韦方程组。对于宏观上与各向异性的物质相关的麦克斯韦方程组，物质的存在会建立一个参考系从而使方程组不再是[[协变]]的。

阅读本条目需要读者了解[[经典电磁理论的协变形式|平直时空中电磁理论的四维形式]]。

电磁场本身要求其几何描述与坐标选取无关，而麦克斯韦方程组在任何时空中的几何描述都是一样的，而不管这个时空是否是平直的。同时，当使用[[笛卡尔坐标|非笛卡尔]]的局部坐标时平直闵可夫斯基空间中的方程组会做同样的修改。例如本条目中方程组可以写成[[球坐标]]中的麦克斯韦方程组的形式。基于上述原因，更好的理解方法是将闵可夫斯基空间中的麦克斯韦方程组理解为一种特殊形式，而非将弯曲时空中的麦克斯韦方程组理解为一种相对论化的推广。

== 摘要 ==

[[广义相对论]]中真空中的电磁理论的方程为

:<math>F_{\alpha \beta} \, = \, \partial_{\alpha} A_{\beta} \, - \, \partial_{\beta} A_{\alpha} \,</math>
:<math>\mathcal{D}^{\mu\nu} \, = \, \frac{1}{\mu_{0}} \, g^{\mu\alpha} \, F_{\alpha\beta} \, g^{\beta\nu} \, \sqrt{-g} \,</math>
:<math>J^{\mu} \, = \, \partial_\nu \mathcal{D}^{\mu \nu} \,</math>
:<math>f_\mu \, = \, F_{\mu\nu} \, J^\nu \,</math>

其中<math>g^{\alpha\beta}\,</math>是度规张量<math>g_{\alpha\beta}\,</math>的倒数，而<math>g\,</math>是[[度规张量]]的行列式，<math>A_{\alpha}\,</math>是电磁场的[[四维势]]，<math>F^{\alpha\beta}\,</math>是[[电磁张量|电磁场的四维协变张量]]，<math>D^{\mu\nu}\,</math>是[[位移电流]]张量，<math>f_{\mu}\,</math>是[[洛伦兹力]]的密度，<math>J_{\mu}\,</math>是[[四维电流密度]]。尽管方程组中使用了[[偏导数]]，这些方程仍然在任意曲面坐标变换下是协变的。也就是说如果将偏导数换成[[协变导数]]，引入的附加项会自动消去从而保持形式不变。

== 电磁四维势 ==

电磁场的[[电磁四维势|四维势]]<math>A_{\alpha} \,,</math>是一个协变矢量，它的坐标变换规则为

:<math>\bar{A}_{\beta} = \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} A_{\gamma} \,.</math>

== 电磁场四维张量 ==

[[电磁张量|电磁场]]是一个协变的二阶反对称张量，它用电磁势可以定义为

:<math>F_{\alpha \beta} \, = \, \partial_{\alpha} A_{\beta} \, - \, \partial_{\beta} A_{\alpha} \,.</math>

为证明它的[[洛伦兹不变性]]，我们对其进行坐标变换

:<math>\bar{F}_{\alpha \beta} \, = \, \frac{\partial \bar{A}_{\beta}}{\partial \bar{x}^{\alpha}} \, - \, \frac{\partial \bar{A}_{\alpha}}{\partial \bar{x}^{\beta}} \,</math>
:<math>= \, \frac{\partial}{\partial \bar{x}^{\alpha}} \left( \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} A_{\gamma} \right) \, - \,  \frac{\partial}{\partial \bar{x}^{\beta}} \left( \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} A_{\delta} \right) \,</math>
:<math>= \, \frac{\partial^2 x^{\gamma}}{\partial \bar{x}^{\alpha} \, \partial \bar{x}^{\beta}} A_{\gamma} \, + \, \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} \frac{\partial A_{\gamma}}{\partial \bar{x}^{\alpha}} \, - \, \frac{\partial^2 x^{\delta}}{\partial \bar{x}^{\beta} \, \partial \bar{x}^{\alpha}} A_{\delta} \, - \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \frac{\partial A_{\delta}}{\partial \bar{x}^{\beta}} \,</math>
:<math>= \, \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \frac{\partial A_{\gamma}}{\partial x^{\delta}} \, - \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} \frac{\partial A_{\delta}}{\partial x^{\gamma}}  \,</math>
:<math>= \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} \, \left( \frac{\partial A_{\gamma}}{\partial x^{\delta}} \, - \, \frac{\partial A_{\delta}}{\partial x^{\gamma}} \right) \,</math>
:<math>= \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \frac{\partial x^{\gamma}}{\partial \bar{x}^{\beta}} \, F_{\delta \gamma} \,.</math>

这一定义暗示了电磁场张量满足关系

:<math>\partial_\lambda F_{\mu \nu} + \partial _\mu F_{\nu \lambda} + \partial_\nu F_{\lambda \mu} = 0 \,</math>

这一关系包含了[[法拉第电磁感应定律]]和[[磁场的高斯定理]]，因此也叫做法拉第-高斯方程。具体而言，

:<math>\partial_\lambda F_{\mu \nu} + \partial _\mu F_{\nu \lambda} + \partial_\nu F_{\lambda \mu} \,</math>
:<math>= \, \partial_\lambda \partial_\mu A_\nu - \partial_\lambda \partial_\nu A_\mu + 
\partial_\mu \partial_\nu A_\lambda - \partial_\mu \partial_\lambda A_\nu + 
\partial_\nu \partial_\lambda A_\mu - \partial_\nu \partial_\mu A_\lambda \, = 0 \,.</math>

虽然这个关系包含了64个分量方程，但只有四个是独立的。借助电磁场张量的反对称性，可知只有<math>\lambda, \mu, \nu\,</math>等于1,2,3或2,3,0或3,0,1或0,1,2时的方程是彼此独立的。

法拉第-高斯方程有时也写作下面的形式

:<math>F_{[\mu \nu ; \lambda]} \, = \, F_{[\mu \nu , \lambda]} \, = \, \frac{1}{6} \left( \partial_\lambda F_{\mu \nu} + \partial _\mu F_{\nu \lambda} + \partial_\nu F_{\lambda \mu} - \partial_\lambda F_{\nu \mu} - \partial _\mu F_{\lambda \nu} - \partial_\nu F_{\mu \lambda} \right) \,</math>
:<math>= \, \frac{1}{3} \left( \partial_\lambda F_{\mu \nu} + \partial _\mu F_{\nu \lambda} + \partial_\nu F_{\lambda \mu} \right) = 0 \,</math>

其中按照惯例用分号表示[[协变导数]]，逗号表示偏导数，方括号表示反对称形式。电磁场张量的协变导数为

:<math>F_{\alpha \beta ; \gamma} \, = \, F_{\alpha \beta , \gamma} - {\Gamma^{\mu}}_{\alpha \gamma} F_{\mu \beta} - {\Gamma^{\mu}}_{\beta \gamma} F_{\alpha \mu} \,</math>

其中<math>\Gamma^{\alpha}_{\beta\gamma}\,</math>是[[克里斯托费尔符号]]，两个下标是对称的。

== 电磁位移张量 ==

[[位移电场]]<math>\mathbf{D} \,,</math>和[[磁场|附屬磁场]]<math>\mathbf{H} \,,</math>构成一个反对称的逆变二阶张量。真空中这个张量为

:<math>\mathcal{D}^{\mu \nu} \, = \, \frac{1}{\mu_{0}} \, g^{\mu \alpha} \, F_{\alpha \beta} \, g^{\beta \nu} \, \sqrt{-g} \,.</math>

注意这个方程是电磁理论中唯一有度规（即引力）存在的方程。并且这一方程具有[[尺度不变性]]，即将度规乘以一个常数不会改变方程的形式。也就是说，引力只能通过改变全局坐标中的[[光速]]来影响相应的电磁场，由于光在引力场中会发生偏折，其效应等同于质量的引力场增加了周围时空的[[折射率]]，从而影响了对应的位移电场和附屬磁场。

更一般地，当物质中的磁化-极化张量不为零时，我们有

:<math>\mathcal{D}^{\mu \nu} \, = \, \frac{1}{\mu_{0}} \, g^{\mu \alpha} \, F_{\alpha \beta} \, g^{\beta \nu} \, \sqrt{-g} \, - \, \mathcal{M}^{\mu \nu} \,.</math>

电磁位移张量的坐标变换规则为

:<math>\bar{\mathcal{D}}^{\mu \nu} \, = \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \,</math>

其中使用了[[雅可比行列式]]。如果磁化-极化张量存在，它和电磁位移张量具有同样的变换规则。

== 电流 ==

电磁位移张量的[[散度]]被定义为[[电流]]，在真空中

:<math>J^{\mu} \, = \, \partial_{\nu} \mathcal{D}^{\mu \nu} \,.</math>

如果磁化-极化张量存在，则上式替换为
:<math>J^{\mu}_{\text{free}} \, = \, \partial_{\nu} \mathcal{D}^{\mu \nu} \,.</math>

这一方程包含了[[高斯定理]]和[[安培环路定理]]。

无论磁化-极化张量是否存在，电磁位移张量是反对称的这一事实暗示了电流是一个守恒量：

:<math>\partial_{\mu} J^{\mu} \, = \, \partial_{\mu} \partial_{\nu} \mathcal{D}^{\mu \nu} = 0 \,</math>

这个二阶偏导数为零是由于偏导是[[对易]]的。

电流的安培-高斯定义并不能决定它的大小，因为电磁四维势的大小并未确定。相反地，通常的步骤是使电流等于一个用其他场表示的表达式（主要是电荷），然后解得电磁位移、电磁场和电磁势。

电流是一个[[逆变]]的矢量，因而其变换规则是

:<math>\bar{J}^{\mu} \, = \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, J^{\alpha} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \,.</math>

变换规则的验证如下：

:<math>\bar{J}^{\mu} \, = \, \frac{\partial}{\partial \bar{x}^{\nu}} \left( \bar{\mathcal{D}}^{\mu \nu} \right) \, = \, \frac{\partial}{\partial \bar{x}^{\nu}} \left( \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \right) \,</math>
:<math>= \, \frac{\partial^2 \bar{x}^{\mu}}{\partial \bar{x}^{\nu} \partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \,</math>
::<math>\frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \frac{\partial \mathcal{D}^{\alpha \beta}}{\partial \bar{x}^{\nu}} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \frac{\partial}{\partial \bar{x}^{\nu}} \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \,</math>
:<math>= \, \frac{\partial^2 \bar{x}^{\mu}}{\partial x^{\beta} \partial x^{\alpha}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \,</math>
::<math>\frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \mathcal{D}^{\alpha \beta}}{\partial x^{\beta}} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \frac{\partial \bar{x}^{\rho}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial \bar{x}^{\nu} \partial \bar{x}^{\rho}} \,</math>
:<math>= \, 0 \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \,</math>
::<math>\frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, J^{\alpha} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \frac{\partial \bar{x}^{\rho}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\rho}} \,</math>
:<math>= \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, J^{\alpha} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \, + \, \frac{\partial \bar{x}^{\mu}}{\partial x^{\alpha}} \, \mathcal{D}^{\alpha \beta} \, \det \left[ \frac{\partial x^{\sigma}}{\partial \bar{x}^{\rho}} \right] \left( \frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, + \, \frac{\partial \bar{x}^{\rho}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\rho}} \right) \,.</math>

下面只需证明

:<math>\frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, + \, \frac{\partial \bar{x}^{\rho}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\rho}} \, = \, 0 \,</math>

这是[[微分学]]中一条已知定理的应用：

:<math>\frac{\partial^2 \bar{x}^{\nu}}{\partial \bar{x}^{\nu} \partial x^{\beta}} \, + \, \frac{\partial \bar{x}^{\rho}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\rho}} \, = \, \frac{\partial x^{\sigma}}{\partial \bar{x}^{\nu}} 
\frac{\partial^2 \bar{x}^{\nu}}{\partial x^{\sigma} \partial x^{\beta}} \, + \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\sigma}} \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\nu}} \,</math>
:<math>= \, \frac{\partial x^{\sigma}}{\partial \bar{x}^{\nu}} 
\frac{\partial^2 \bar{x}^{\nu}}{\partial x^{\beta} \partial x^{\sigma}} \, + \, \frac{\partial^2 x^{\sigma}}{\partial x^{\beta} \partial \bar{x}^{\nu}} \frac{\partial \bar{x}^{\nu}}{\partial x^{\sigma}} \, = \, \frac{\partial}{\partial x^{\beta}} \left( \frac{\partial x^{\sigma}}{\partial \bar{x}^{\nu}} \, \frac{\partial \bar{x}^{\nu}}{\partial x^{\sigma}} \right) \,</math>
:<math>= \, \frac{\partial}{\partial x^{\beta}} \left(  \, \frac{\partial \bar{x}^{\nu}}{\partial \bar{x}^{\nu}} \right) \, = \, \frac{\partial}{\partial x^{\beta}} \left( \mathbf{4} \right) \, = \, 0 \,.</math>

== 洛伦兹力 ==

[[洛伦兹力]]的密度是一个协变矢量：

:<math>f_{\mu} \, = \, F_{\mu \nu} \, J^{\nu} \,.</math>

如果在一个测试粒子上的作用力只有引力和电磁力，则有

:<math>\frac{d p_{\alpha}}{d t} \, = \, \Gamma^{\beta}_{\alpha \gamma} \, p_{\beta} \, \frac{d x^{\gamma}}{d t} \, +  \, q \, F_{\alpha \gamma} \, \frac{d x^{\gamma}}{d t} \,</math>

其中<math>p^{\alpha}\,</math>是粒子的[[四维动量]]，<math>t\,</math>是将粒子[[世界线]]参数化的任意时间坐标。式中第一项含有克里斯托费尔符号，表示引力项，第二项含有电荷，表示电磁力项。

方程具有洛伦兹不变性，对时间坐标变换的验证方法为将方程乘以<math>\frac{d t}{d \bar{t}} \,</math>并使用[[链式法则]]。

对空间坐标变换，通过对克里斯托费尔符号进行坐标变换

:<math>\bar{\Gamma}^{\beta}_{\alpha \gamma} \, = \, 
\frac{\partial \bar{x}^{\beta}}{\partial x^{\epsilon}} \,
\frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \,
\frac{\partial x^{\zeta}}{\partial \bar{x}^{\gamma}} \,
\Gamma^{\epsilon}_{\delta \zeta} \,
+ 
\frac{\partial \bar{x}^{\beta}}{\partial x^{\eta}}\, 
\frac{\partial^2 x^{\eta}}{\partial \bar{x}^{\alpha} \partial \bar{x}^{\gamma}} \,</math>

我们得到

:<math>\frac{d \bar{p}_{\alpha}}{d t} \, - \, \bar{\Gamma}^{\beta}_{\alpha \gamma} \, \bar{p}_{\beta} \, \frac{d \bar{x}^{\gamma}}{d t} \, -  \, q \, \bar{F}_{\alpha \gamma} \, \frac{d \bar{x}^{\gamma}}{d t} \,</math>
:<math>= \, \frac{d}{d t} \left( \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \, p_{\delta} \right) \, - \, 
\left( 
\frac{\partial \bar{x}^{\beta}}{\partial x^{\theta}} \,
\frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \,
\frac{\partial x^{\iota}}{\partial \bar{x}^{\gamma}} \,
\Gamma^{\theta}_{\delta \iota} \,
+ \, \frac{\partial \bar{x}^{\beta}}{\partial x^{\eta}}\, 
\frac{\partial^2 x^{\eta}}{\partial \bar{x}^{\alpha} \partial \bar{x}^{\gamma}}  
\right) \, \frac{\partial x^{\epsilon}}{\partial \bar{x}^{\beta}} \, p_{\epsilon} \, \frac{\partial \bar{x}^{\gamma}}{\partial x^{\zeta}} \, \frac{d x^{\zeta}}{d t} \, -  \, q \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \, F_{\delta \zeta} \, \frac{d x^{\zeta}}{d t} \,</math>
:<math>= \, \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \, \left(
\frac{d p_{\delta}}{d t} \, - \, \Gamma^{\epsilon}_{\delta \zeta} \, p_{\epsilon} \, \frac{d x^{\zeta}}{d t} \, -  \, q \, F_{\delta \zeta} \, \frac{d x^{\zeta}}{d t} \right) \, + \,</math>
::<math>\frac{d}{d t} \left( \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \right) \, p_{\delta} \, - \, 
\left( \frac{\partial \bar{x}^{\beta}}{\partial x^{\eta}}\, 
\frac{\partial^2 x^{\eta}}{\partial \bar{x}^{\alpha} \partial \bar{x}^{\gamma}}  
\right) \, \frac{\partial x^{\epsilon}}{\partial \bar{x}^{\beta}} \, p_{\epsilon} \, \frac{\partial \bar{x}^{\gamma}}{\partial x^{\zeta}} \, \frac{d x^{\zeta}}{d t} \,</math>
:<math>= \, 0 \, + \, \frac{d}{d t} \left( \frac{\partial x^{\delta}}{\partial \bar{x}^{\alpha}} \right) \, p_{\delta} \, - \, 
\frac{\partial^2 x^{\epsilon}}{\partial \bar{x}^{\alpha} \partial \bar{x}^{\gamma}} p_{\epsilon} \, \frac{d \bar{x}^{\gamma}}{d t} \, = \, 0 \,</math>

== 拉格朗日量 ==

在真空中经典电磁场的[[拉格朗日量]]（在密度意义下的单位为焦耳/米<sup>3</sup>）是一个标量：

:<math> \mathcal{L} \, = \, - \frac{1}{4 \mu_0} \, F_{\alpha \beta} \, F^{\alpha \beta} \, \sqrt{- g} \, + \, A_{\alpha} \, J^{\alpha} \,</math>

其中<math>F^{\alpha \beta} = g^{\alpha \gamma} F_{\gamma \delta} g^{\delta \beta} \,.</math>。这里的四维电流项应理解为各种对电流有贡献的场的总和。

如果我们将自由电流与束缚电流分离开来，拉格朗日量则变为

:<math> \mathcal{L} \, = \, - \frac{1}{4 \mu_0} \, F_{\alpha \beta} \, F^{\alpha \beta} \, \sqrt{- g} \, + \,A_{\alpha} \, J^{\alpha}_{\text{free}} \, + \, \frac12 \, F_{\alpha \beta} \, \mathcal{M}^{\alpha \beta} \,.</math>

== 电磁应力-能量张量 ==
{{main|电磁应力-能量张量}}

作为[[爱因斯坦引力场方程]]的源，电磁场的[[应力-能量张量]]是一个协变的对称张量

:<math>T_{\mu \nu} \, = \, - \frac{1}{\mu_0} ( F_{\mu \alpha} g^{\alpha \beta} F_{\beta \nu} \, - \, \frac{1}{4} g_{\mu \nu} \, F_{\sigma \alpha} g^{\alpha \beta} F_{\beta \rho} g^{\rho \sigma} ) \,</math>

并且它是无迹的

:<math>T_{\mu \nu} g^{\mu \nu} \, = \, 0 \,</math>

这是由于电磁场的传播速度在不同参考系下是不变的。

为了表达能量和动量的守恒律，电磁应力-能量张量的最佳表示方法是

:<math>\mathfrak{T}_{\mu}^{\nu} = T_{\mu \gamma} \, g^{\gamma \nu} \, \sqrt{-g}.</math>

对于上式可以证明

:<math>{\mathfrak{T}_{\mu}^{\nu}}_{; \nu} \, + \, f_{\mu} \, = \, 0 \,</math>

并可重写为

:<math>- {\mathfrak{T}_{\mu}^{\nu}}_{, \nu} \, = \, - \Gamma^{\sigma}_{\mu \nu} \mathfrak{T}_{\sigma}^{\nu} \, + \, f_{\mu} \,</math>

这个方程的意义是，电磁场能量的减少相当于电磁场对引力场以及通过洛伦兹力对物质所做的功，类似地，电磁场动量的减少率相当于作用在引力场上的电磁力以及作用在物质上的洛伦兹力。

守恒律的推导过程为

:<math>{\mathfrak{T}_{\mu}^{\nu}}_{; \nu} \, + \, f_{\mu} \, = \, - \frac{1}{\mu_0} ( F_{\mu \alpha ; \nu} g^{\alpha \beta} F_{\beta \gamma} g^{\gamma \nu} \, + \, F_{\mu \alpha} g^{\alpha \beta} F_{\beta \gamma ; \nu} g^{\gamma \nu} \, - \, \frac12 \delta_{\mu}^{\nu} \, F_{\sigma \alpha ; \nu} g^{\alpha \beta} F_{\beta \rho} g^{\rho \sigma} ) \sqrt{- g} \, + \,</math>
::<math>\frac{1}{\mu_{0}} \, F_{\mu \alpha} \, g^{\alpha \beta} \, F_{\beta \gamma ; \nu} \, g^{\gamma \nu} \, \sqrt{-g} \,</math>
:<math>= \, - \frac{1}{\mu_0} ( F_{\mu \alpha ; \nu} F^{\alpha \nu} \, - \, \frac12 F_{\sigma \alpha ; \mu} F^{\alpha \sigma} ) \sqrt{- g} \,</math>
:<math>= \, - \frac{1}{\mu_0} ( (- F_{\nu \mu ; \alpha} - F_{\alpha \nu ; \mu}) F^{\alpha \nu} \, - \, \frac12 F_{\sigma \alpha ; \mu} F^{\alpha \sigma} ) \sqrt{- g} \,</math>
:<math>= \, - \frac{1}{\mu_0} ( F_{\mu \nu ; \alpha} F^{\alpha \nu} - F_{\alpha \nu ; \mu} F^{\alpha \nu} \, + \, \frac12 F_{\sigma \alpha ; \mu} F^{\sigma \alpha} ) \sqrt{- g} \,</math>
:<math>= \, - \frac{1}{\mu_0} ( F_{\mu \alpha ; \nu} F^{\nu \alpha} - \frac12 F_{\alpha \nu ; \mu} F^{\alpha \nu} ) \sqrt{- g} \,</math>
:<math>= \, - \frac{1}{\mu_0} (- F_{\mu \alpha ; \nu} F^{\alpha \nu} \, + \, \frac12 F_{\sigma \alpha ; \mu} F^{\alpha \sigma} ) \sqrt{- g} \,</math>

所得的结果是为零的，因为它等于自身的负值（对照推导的第二步）。

== 电磁波方程 ==

从[[麦克斯韦方程组在狭义相对论中的形式|电磁理论的狭义相对论形式]]可以借助电磁场张量修改得到[[非齐次的电磁波方程]]

:<math> \Box F_{ab} \ \stackrel{\mathrm{def}}{=}\  F_{ab;}{}^d{}_d = \, -2 R_{acbd} F^{cd} + R_{ae}F^e{}_b - R_{be}F^e{}_a + J_{a;b} - J_{b;a}  </math>

其中

:<math>  R_{acbd} \,</math>

是[[黎曼张量]]的协变形式，而<math>  \Box \,</math>是[[达朗贝尔算符]]在协变导数情形下的推广。

波方程可以写成四维势的形式（参见ref 2, p. 569）

:<math>  \Box A^{a}  =  {{A^{a; }}^{b}}_{ b}  = - \mu_0 J^{ a }  + {R^{ a }}_{ b } A^{ b } </math>

其中 

:<math>  {R^{ a }}_{ b }  \ \stackrel{\mathrm{def}}{=}\  {R^{ s }}_{ a s b } </math>

是[[里奇张量]]。而这里假设了[[洛伦茨规范]]在弯曲时空中的推广具有下面的形式

:<math>  {A^{ a  }}_{ ; a } = 0 </math>

这个波方程与平直时空中的波方程具有十分类似的形式，除了导数被替换为协变导数，以及在表示源的项中多了一项描述时空曲率的项。这个方程和弯曲时空中的洛伦兹力也有相似之处，这里的四维势<math>  {A^{ a  }} \, </math>相当于洛伦兹力中的四维坐标。

== 麦克斯韦方程组在动态时空中的非线性 ==

在考虑麦克斯韦方程组本身与背景时空无关时，时空度规在广义相对论中被认为是受电磁场影响的动态变化的变量，这导致了电磁波方程和麦克斯韦方程组是非线性的。这一点可以从[[爱因斯坦场方程]]中的曲率张量依赖于[[应力-能量张量]]看出。

:<math>   G_{ab} =  \frac{8 \pi G} {c^4}   T_{ab}    </math>
 
其中
:<math>  {G}_{a  b} \ \stackrel{\mathrm{def}}{=}\        {R}_{a b} - {1 \over 2} {R} g_{a b}   </math> 

是[[爱因斯坦张量]]，<math>G\,</math>是[[万有引力常数]]，<math>R\,</math> （[[标量曲率]]）是里奇张量的迹。应力-能量张量来源于粒子的[[应力]]和能量，但同时也来源于电磁场，这造就了非线性。

== 参见 ==

* [[电磁波方程]]
* [[非齐次的电磁波方程]]
* [[麦克斯韦方程组在狭义相对论中的形式]]
* [[麦克斯韦方程组的历史]]
* [[广义相对论的理论动机]]
* [[广义相对论的数学方法]]
* [[电真空解]]

== 参考文献 ==

*{{lang|en|{{cite book | author=Einstein, A.  | title=Relativity: The Special and General Theory | url=https://archive.org/details/relativityspecia00eins_0  | location= New York | publisher=Crown| year=1961 | isbn=0-517-02961-8}}}}
*{{lang|en|{{cite book | author=Misner, Charles; Thorne, Kip S. & Wheeler, John Archibald | title=Gravitation | location=San Francisco | publisher=W. H. Freeman | year=1973 | isbn=0-7167-0344-0}}}}
*{{lang|en|{{cite book | author=Landau, L. D. and Lifshitz, E. M.| title=Classical Theory of Fields (Fourth Revised English Edition) | url=https://archive.org/details/classicaltheoryo0000land_k6k2| location=Oxford | publisher=Pergamon | year=1975 | isbn=0-08-018176-7}}}}
*{{lang|en|{{cite book | author=R. P. Feynman, F. B. Moringo, and W. G. Wagner | title=Feynman Lectures on Gravitation | url=https://archive.org/details/feynmanlectureso0000feyn_g4q1 | publisher=Addison-Wesley | year=1995 | isbn=0-201-62734-5}}}}

== 外部链接 ==
* [http://www.iop.org/EJ/abstract/0264-9381/22/2/011/ 弯曲时空中的麦克斯韦场] 

{{Physics-footer}}

[[Category:基本物理概念|W]]
[[Category:广义相对论|W]]
[[Category:偏微分方程|W]]