--从 MOC 统一曲率极值原理推导万有引力定律

作者：张苏杭  洛阳

一、MOC 统一曲率极值原理（公设）

设时空流形配备统一曲率标量 \mathcal{R}_{\text{总}} ，它由引力曲率标量 R_{\text{grav}} 与规范场曲率贡献 R_{\text{gauge}} 共同构成：

\mathcal{R}_{\text{总}} = R_{\text{grav}} + R_{\text{gauge}}

MOC 极值原理表述为：物理上允许的场构型使总曲率标量的时空积分取极值：

\boxed{\delta \int \mathcal{R}_{\text{总}} \, \sqrt{-g} \, d^4x = 0}

该变分对全部自由度（度规、规范场、物质场）进行，无需额外假设作用量形式。

---

二、投影到纯引力子流形（规范场为零）

在只考虑引力相互作用的区域（如星体外部真空，电磁/弱/强力可忽略），规范曲率贡献为零： R_{\text{gauge}} = 0 ，则

\mathcal{R}_{\text{总}} = R_{\text{grav}} \equiv R

 R 为时空的黎曼曲率标量。MOC 极值原理退化为：

\delta \int R \, \sqrt{-g} \, d^4x = 0

---

三、变分推导爱因斯坦场方程

对上述泛函关于度规 g_{\mu\nu} 变分，标准计算给出：

\int \left( R_{\mu\nu} - \frac12 g_{\mu\nu} R \right) \delta g^{\mu\nu} \sqrt{-g} \, d^4x = 0

考虑到物质场作用量时，右端增加 8\pi G \, \delta S_m / \delta g^{\mu\nu} ，最终得到：

\boxed{R_{\mu\nu} - \frac12 g_{\mu\nu} R = 8\pi G \, T_{\mu\nu}}

即爱因斯坦场方程。其中 G 为牛顿引力常数，由与物质场的耦合尺度自然出现。

---

四、真空球对称静态解（史瓦西度规）

在星体外部真空 T_{\mu\nu}=0 ，且要求球对称、静态，史瓦西度规为方程的解：

ds^2 = -\left(1 - \frac{2GM}{c^2 r}\right) c^2 dt^2 + \left(1 - \frac{2GM}{c^2 r}\right)^{-1} dr^2 + r^2 d\Omega^2

---

五、弱场低速牛顿极限

当 \frac{2GM}{c^2 r} \ll 1 （弱场），且粒子速度远小于光速，度规时间分量近似为：

g_{00} \approx -\left(1 + \frac{2\phi}{c^2}\right), \quad \phi(r) = -\frac{GM}{r}

其中 \phi(r) 为牛顿引力势。

---

六、导出万有引力公式

在牛顿极限下，试验粒子在引力场中的加速度为：

\boldsymbol{a} = -\nabla \phi = -\frac{d\phi}{dr} \hat{\boldsymbol{r}} = -\frac{GM}{r^2} \hat{\boldsymbol{r}}

试探质量 m 所受引力为：

\boldsymbol{F} = m \boldsymbol{a} = -\,\frac{GMm}{r^2} \hat{\boldsymbol{r}}

取大小即得万有引力定律：

\boxed{F = G \frac{Mm}{r^2}}

---

总结：完整逻辑链

\begin{aligned}

&\text{MOC 极值原理: } \delta \int \mathcal{R}_{\text{总}} \sqrt{-g}\,d^4x = 0 \\

&\quad \xrightarrow{\text{纯引力子流形投影}} \delta \int R \sqrt{-g}\,d^4x = 0 \\

&\quad \xrightarrow{\text{变分}} \text{爱因斯坦场方程} \\

&\quad \xrightarrow{\text{真空球对称静态解}} \text{史瓦西度规} \\

&\quad \xrightarrow{\text{弱场低速极限}} \text{牛顿势 } \phi = -GM/r \\

&\quad \xrightarrow{\text{加速度 } a = -\nabla\phi} \text{万有引力定律 } F = G\frac{Mm}{r^2}

\end{aligned}