扩展引力现象

Learning objectives / 学习目标

完成本课后，你应当能够：

1. 从引力场的一阶 Taylor 展开得到潮汐张量，并解释径向拉伸、横向压缩与自由落体中的可测相对加速度。
2. 区分简单差分引力估计、刚性自引力模型和同步流体 Roche 模型的系数及假设。
3. 从逐层组装推导球体引力自能，说明密度剖面对系数的影响。
4. 写出 Newtonian \(N\) 体方程、摄动加速度与圆形限制性三体问题的基本结构。
5. 解释 Lagrange 点、轨道摄动和混沌敏感性，而不把它们误写成静态“引力相消”或能量不守恒。
6. 用 \(GM/(rc^2)\) 和 \(v/c\) 判断 Newtonian 近似，并识别需要广义相对论的现象。

Prerequisites / 先修知识

• Newtonian Gravitation：球对称场、势与自引力能量语言。
• Orbits and Kepler's Laws 与 Central-Force Methods：二体轨道、有效势和摄动前的基准解。
• 会使用一阶 Taylor 展开、矩阵本征方向、常微分方程和基本数值误差概念；广义相对论只作尺度边界，不推导场方程。

Notation and modeling contract / 符号与建模契约

0. What changes beyond the point-mass two-body model? / 超出点质量二体模型后改变什么

Knowledge 1 — Tidal forces and Roche limits / 潮汐力与洛希极限

Knowledge ID: P01-U07-L04-K01

1.1 Tides are differential gravity

自由落体质心处仍可有很强的共同加速度；有限物体真正“感觉”到的是不同位置的场差：

\[\delta\mathbf g =\mathbf g(\mathbf R+\boldsymbol\xi)-\mathbf g(\mathbf R) \approx(\boldsymbol\xi\cdot\boldsymbol\nabla)\mathbf g.\]

对点质量外场

\[\mathbf g(\mathbf R)=-\frac{GM}{R^3}\mathbf R.\]

定义潮汐张量

\[T_{ij}=\frac{\partial g_i}{\partial x_j} =\frac{GM}{R^3}\left(3n_in_j-\delta_{ij}\right), \qquad \delta g_i\approx T_{ij}\xi_j,\]

其中 \(\mathbf n=\mathbf R/R\) 从源指向质心。沿径向的本征值为 \(+2GM/R^3\)，两个横向本征值各为 \(-GM/R^3\)：径向拉伸、横向压缩。在源外真空中迹为零；在有质量密度处，Poisson 方程改变这一结果。

1.2 Exact difference and linear range

若两点沿径向相隔小距离 \(\Delta r\)，

\[\Delta g =GM\left[\frac1{R^2}-\frac1{(R+\Delta r)^2}\right] \approx\frac{2GM}{R^3}\Delta r.\]

误差量级为 \(O[(\Delta r/R)^2]\)。当物体大小不再远小于 \(R\)，必须使用精确场差，甚至求解物体自身形变后的引力与应力。

1.3 Worked Example 1 — lunar tidal acceleration across Earth

已知。 月球标准引力参数 \(GM_{\mathrm M}=4.9049\times10^{12}\ \mathrm{m^3/s^2}\)、地月中心距 \(R=3.844\times10^8\ \mathrm m\)、地球半径 \(R_\oplus=6.371\times10^6\ \mathrm m\)。

研究对象、系统与参考系。 比较地球近月表面点与地心的月球引力；先在月心惯性系求两场之差，等价地在随地心共同自由落体的局部系中求相对加速度。

坐标与正方向。 径向单位矢量 \(\mathbf n\) 从月心指向地心；近月点的内部坐标为 \(\boldsymbol\xi=-R_\oplus\mathbf n\)。下式报告相对加速度的大小。

模型及理由。 月球视为点质量，且 \(R_\oplus/R=0.0166\ll1\)，可把场差保留到 \(\boldsymbol\xi/R\) 的一阶潮汐项。

方程与计算。 径向潮汐本征值为 \(2GM_{\mathrm M}/R^3\)，所以

\[\delta a_r\approx\frac{2GM_{\mathrm M}R_\oplus}{R^3} =\frac{ 2\left(4.9049\times10^{12}\ \mathrm{m^3\,s^{-2}}\right) \left(6.371\times10^6\ \mathrm m\right)} {\left(3.844\times10^8\ \mathrm m\right)^3} =1.10\times10^{-6}\ \mathrm{m/s^2}.\]

SI、有效数字与检查。 \([GM R_\oplus/R^3]=\mathrm{m/s^2}\)；近月点的相对加速度方向为 \(-\mathbf n\)，即相对地心更朝月球。以点质量场直接作差可得

\[\delta a_{\rm exact} =GM_{\mathrm M}\left[ \frac{1}{(R-R_\oplus)^2}-\frac{1}{R^2} \right] =1.13\times10^{-6}\ \mathrm{m/s^2}.\]

令 \(\epsilon=R_\oplus/R=0.0166\)，精确展开为背景场乘 \(2\epsilon+3\epsilon^2+\cdots\)；所以下一项相对领先潮汐项约为 \(3\epsilon/2=2.49\%\)，一阶值相对精确场差低约 \(2.48\%\)。因此 \(1.10\times10^{-6}\ \mathrm{m/s^2}\) 是受控的百分之几精度估计，而不是三位数都精确的测量值。

物理解释。 该量远小于月球对地心的共同加速度，却可持续驱动潮汐；潮汐是减去质心共同自由落体后留下的场梯度，不是“月球只拉近侧”。

1.4 Roche limits are model-dependent

对主天体半径 \(R_M\)、密度 \(\rho_M\) 和卫星密度 \(\rho_s\)，把卫星半径处的线性潮汐加速度与自身表面引力相等，得到简单估计

\[d_{\mathrm{simple}} =2^{1/3}R_M \left(\frac{\rho_M}{\rho_s}\right)^{1/3} \approx1.26R_M \left(\frac{\rho_M}{\rho_s}\right)^{1/3}.\]

经典同步、不可压、无材料强度的流体卫星模型给出更大的系数：

\[d_{\mathrm{fluid}} \approx2.44R_M \left(\frac{\rho_M}{\rho_s}\right)^{1/3}.\]

\(1.26\) 常被称为保持球形的刚性自引力估计，\(2.44\) 是可形变同步流体极限。真实破坏还依赖抗拉强度、孔隙率、黏度、自转、轨道偏心率、内部结构和通过极限的时间；小而坚固的航天器可在“流体 Roche 极限”内完好存在。

1.5 Worked Example 2 — two Roche estimates for the Moon

已知。 主天体取地球：\(R_\oplus=6371\ \mathrm{km}=6.371\times10^6\ \mathrm m\)、\(\rho_\oplus=5514\ \mathrm{kg/m^3}\)；卫星取月球平均密度 \(\rho_{\mathrm M}=3344\ \mathrm{kg/m^3}\)。

研究对象、系统与参考系。 研究一个月球密度卫星绕地球的潮汐破坏中心距；在随圆轨道转动的框架中比较两个不同内部响应模型。

坐标与正方向。 \(d\) 从地心量到卫星质心，径向向外为正；两结果都是中心距，不是离地高度。

模型及理由。 分别采用保持球形的刚性自引力估计和同步、不可压、无材料强度的流体模型；相同密度比但不同系数正是要比较的物理假设。

方程与计算。 密度因子为

\[\left(\frac{\rho_\oplus}{\rho_{\mathrm M}}\right)^{1/3} =\left( \frac{5514\ \mathrm{kg\,m^{-3}}}{3344\ \mathrm{kg\,m^{-3}}} \right)^{1/3} =1.181.\]

于是

\[d_{\mathrm{rigid}} \approx1.26\left(6.371\times10^6\ \mathrm m\right)(1.181) =9.48\times10^6\ \mathrm m,\]

\[d_{\mathrm{fluid}} \approx2.44\left(6.371\times10^6\ \mathrm m\right)(1.181) =1.84\times10^7\ \mathrm m.\]

SI、有效数字与检查。 密度比及其立方根无量纲，故 \(d\) 的单位继承 \(R_\oplus\)；结果取三位有效数字。两值均大于 \(R_\oplus\) 且流体极限更远，符合可形变体更易被潮汐破坏的预期；也都远小于当前地月距离。

物理解释。 数值差异不是“哪个公式错了”，而是内部强度、形变和同步自转假设不同；真实天体还需材料与结构模型。

Knowledge 2 — Gravitational self-energy / 引力自能

Knowledge ID: P01-U07-L04-K02

2.1 Assemble a sphere shell by shell

引力自能是把连续质量从无穷远逐步组装成给定分布时，外力所做的准静态功。已组装的包围质量为 \(M(r)\)，再放入薄壳 \(dM\)：

\[dW=-\frac{GM(r)}{r}\,dM.\]

所以一般球对称分布

\[W=-\int_0^M\frac{GM(r)}{r}\,dM.\]

等价的连续形式为

\[W=\frac12\int\rho(\mathbf r)\Phi(\mathbf r)\,d^3r,\]

其中 \(1/2\) 防止质量对被重复计数。\(W<0\) 表示束缚；完全拆散到无穷远至少需要输入 \(-W\)，若过程中结构和内能不变。

2.2 Uniform-density sphere

对常密度球

\[M(r)=M\frac{r^3}{R^3}, \qquad dM=3M\frac{r^2}{R^3}\,dr.\]

因此

\[W=-\int_0^R \frac{GMr^3/R^3}{r} \left(3M\frac{r^2}{R^3}\,dr\right) =-\frac{3GM^2}{5R}.\]

系数 \(3/5\) 只属于均匀密度。中心凝聚越强，通常束缚能的绝对值越大；不能对任意恒星、行星都机械套用 \(3/5\)。

2.3 Worked Example 3 — uniform-Earth estimate

已知。 \(G=6.67430\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}}\)、\(M_\oplus=5.9722\times10^{24}\ \mathrm{kg}\)、\(R_\oplus=6.371\times10^6\ \mathrm m\)，\(c=2.9979\times10^8\ \mathrm{m/s}\)。

研究对象、系统与参考系。 研究地球全部质量元组成的自引力系统；在地心球坐标中从无穷远准静态组装各薄壳，整体平动参考系不影响内部自能。

坐标与正方向。 壳半径 \(r\) 从地心向外增加；组装使系统势能降低，故自能应为负。

模型及理由。 把地球近似成静止、均匀密度球，才能使用系数 \(3/5\)；该模型用于尺度估计而非精密地球结构计算。

方程与计算。 均匀球逐壳积分给出

\[W=-\frac{3G M_\oplus^2}{5R_\oplus} =-\frac{ 3\left(6.67430\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}}\right) \left(5.9722\times10^{24}\ \mathrm{kg}\right)^2} {5\left(6.371\times10^6\ \mathrm m\right)} =-2.2419\times10^{32}\ \mathrm J \approx-2.24\times10^{32}\ \mathrm J.\]

相对静质量能的比例为

\[\frac{\lvert W\rvert}{M_\oplus c^2} =\frac{2.2419\times10^{32}\ \mathrm J} {\left(5.9722\times10^{24}\ \mathrm{kg}\right) \left(2.9979\times10^8\ \mathrm{m\,s^{-1}}\right)^2} =4.18\times10^{-10}.\]

SI、有效数字与检查。 \(GM^2/R\) 的单位为 \(\mathrm{kg\,m^2/s^2}=\mathrm J\)，结果按输入取三位有效数字。负号表示束缚；若保持质量而增大 \(R\)，\(\lvert W\rvert\propto1/R\) 应减小。极小的无量纲比值确认整体弱场。

物理解释。 拆散理想均匀地球至少需输入约 \(2.24\times10^{32}\ \mathrm J\)；真实地球中心更致密，均匀球只给第一尺度估计。

Knowledge 3 — Many-body orbital dynamics and perturbations / 多体轨道与摄动

Knowledge ID: P01-U07-L04-K03

3.1 Newtonian \(N\)-body equations

孤立点质量系统满足

\[\ddot{\mathbf r}_i =G\sum_{j\ne i}m_j \frac{\mathbf r_j-\mathbf r_i} {\lvert\mathbf r_j-\mathbf r_i\rvert^3}.\]

总能量、总线动量和总角动量仍守恒，但 \(N\ge3\) 一般没有类似 Kepler 二体问题的统一闭式解。守恒量降低自由度，却不等于每条轨道可用初等函数写出。

3.2 Perturbing acceleration and osculating orbits

在日心坐标中，行星 \(p\) 对测试体相对太阳运动的摄动加速度为

\[\mathbf a_{\mathrm{pert}} =Gm_p\left[ \frac{\mathbf r_p-\mathbf r}{\lvert\mathbf r_p-\mathbf r\rvert^3} -\frac{\mathbf r_p}{r_p^3} \right].\]

第一项是行星对测试体的直接引力，第二项是行星使坐标原点太阳加速的 indirect term。漏掉第二项就没有写出相对运动。摄动理论把真实轨道视为一系列瞬时相切的 osculating Kepler orbits，其轨道根数随时间缓慢变化。

3.3 Restricted three-body structure and Lagrange points

在 circular restricted three-body problem 中，两主质量作圆运动，第三体质量可忽略。在随两主质量转动的坐标系中，定义

\[\Omega(x,y) =\frac{Gm_1}{r_1}+\frac{Gm_2}{r_2} +\frac12n^2(x^2+y^2).\]

方程为

\[\ddot x-2n\dot y=\frac{\partial\Omega}{\partial x}, \qquad \ddot y+2n\dot x=\frac{\partial\Omega}{\partial y}.\]

\(\boldsymbol\nabla\Omega=0\) 给五个 Lagrange 平衡点。\(L_1,L_2,L_3\) 线性不稳定；\(L_4,L_5\) 在主质量比足够小时可线性稳定。定义 \(\eta=m_2/(m_1+m_2)\le1/2\)，三角点的线性稳定条件为

\[27\eta(1-\eta)<1, \qquad \eta<\eta_R=\frac12\left(1-\sqrt{\frac{23}{27}}\right)\approx0.03852.\]

它们是旋转系中的动力学平衡，包含离心结构；不能解释为两引力简单相消。实际航天器常在 Lagrange 点附近的 halo 或 Lissajous 轨道运行，并需 station keeping。

3.4 Chaos and numerical prediction

某些三体和多体轨道对初值敏感：早期相邻轨迹差可近似

\[\lVert\delta\mathbf X(t)\rVert \sim\lVert\delta\mathbf X(0)\rVert e^{\lambda t}, \qquad \lambda>0.\]

这叫确定性混沌，不等于随机力，也不取消守恒定律。长期预测应报告初值不确定度、模型误差和数值误差；保守轨道长期积分常选 symplectic integrator，以减少能量系统漂移。

3.5 Worked Example 4 — Jupiter as a perturbation on Earth

已知。 Earth–Sun 距离 \(a=1.00\ \mathrm{AU}\)，Sun–Jupiter 距离取 \(D=5.2\ \mathrm{AU}\)；最近共线时 Earth–Jupiter 距离为 \(D-a=4.2\ \mathrm{AU}\)。质量比 \(M_J/M_\odot=9.54\times10^{-4}\)。

研究对象、系统与参考系。 研究 Earth 相对 Sun 的运动；在日心相对坐标中，必须从 Jupiter 对 Earth 的直接加速度减去 Jupiter 对 Sun 的加速度。

坐标与正方向。 取瞬时 Sun–Jupiter 连线为径向轴，并以 Sun 指向 Jupiter 为正；最近共线时 Earth 位于两者之间，径向潮汐本征值带系数 \(2\)。

模型及理由。 三天体视为点质量；先把 Jupiter 场围绕 Sun 在尺度 \(a\) 上作一阶 Taylor 展开，因此展开分母必须是 Sun–Jupiter 距离 \(D\)，径向摄动为 \(2GM_Ja/D^3\)。由于 \(a/D=0.192\) 并非极小，随后再用精确的 direct-minus-indirect 差值检查一阶估计。

方程与计算。 相对于太阳主加速度 \(a_\odot=GM_\odot/a^2\)，线性径向潮汐估计为

\[\frac{a_{\mathrm{pert}}}{a_\odot} \approx 2\frac{M_J}{M_\odot}\left(\frac{a}{D}\right)^3 =2\left(9.54\times10^{-4}\right) \left(\frac{1.00\ \mathrm{AU}}{5.2\ \mathrm{AU}}\right)^3 =1.36\times10^{-5}.\]

最近共线时，精确的直接项减间接项给出

\[\left(\frac{a_{\mathrm{pert}}}{a_\odot}\right)_{\rm exact} =\frac{M_J}{M_\odot}a^2 \left[\frac{1}{(D-a)^2}-\frac{1}{D^2}\right] =\left(9.54\times10^{-4}\right)(1.00\ \mathrm{AU})^2 \left[ \frac{1}{(4.2\ \mathrm{AU})^2} -\frac{1}{(5.2\ \mathrm{AU})^2} \right] =1.88\times10^{-5}.\]

SI、有效数字与检查。 每个距离都以同一单位 AU 代入，质量比与整个加速度比均无量纲。线性值比最近共线的精确差值低约 \(28\%\)，与 \(a/D=0.192\) 只提供量级级别的一阶展开相符。若 \(M_J\to0\) 或 \(D\to\infty\)，比值趋零；若漏掉 indirect term，远处行星造成的共同加速度不会相消，便得不到 \(D^{-3}\) 的潮汐极限。

物理解释。 瞬时比值虽小，周期性摄动仍可长期累积为近日点进动、轨道根数振荡或共振。一阶潮汐式揭示 \(D^{-3}\) 标度，精确差值则提醒：当 \(a/D\) 不很小时，不能把量级展开冒充精密数值。

Knowledge 4 — Boundary of validity of Newtonian gravity / 牛顿引力的适用边界

Knowledge ID: P01-U07-L04-K04

4.1 Two small parameters

Newtonian 引力通常要求

\[\epsilon_g=\frac{GM}{rc^2}\ll1, \qquad \beta=\frac vc\ll1.\]

\(\epsilon_g\) 衡量引力势相对 \(c^2\) 的大小，\(\beta\) 衡量运动相对论性。弱场、低速、无需传播延迟或极高精度时，Newtonian 模型极其成功。以下现象必须使用广义相对论或其后 Newtonian 展开：强场致密天体、黑洞视界、光线弯曲、引力时间膨胀、引力波、精密近心进动和高速双星。

4.2 Schwarzschild scale and a warning

定义

\[r_s=\frac{2GM}{c^2}.\]

把 Newtonian 逃逸速度形式设成 \(c\) 会碰巧得到同一代数尺度，但不能由此推导黑洞：Newtonian 理论没有时空曲率、光锥、事件视界或正确的光传播。\(r_s\) 的物理意义来自 Schwarzschild 几何。

4.3 Relativistic perihelion advance

对围绕球对称非旋转质量的弱场椭圆轨道，GR 的一阶额外近心进动为每圈

\[\Delta\varpi =\frac{6\pi GM}{a(1-e^2)c^2}.\]

这不是 Newtonian 多体摄动；实际比较观测前必须先扣除其他行星、太阳扁率等经典效应。

4.4 Worked Example 5 — Mercury's relativistic scale

已知。 Mercury 的 \(a=5.791\times10^{10}\ \mathrm m\)、\(e=0.2056\)，太阳 \(GM_\odot=1.3271\times10^{20}\ \mathrm{m^3/s^2}\)，\(c=2.9979\times10^8\ \mathrm{m/s}\)；每世纪约 \(415\) 圈。

研究对象、系统与参考系。 研究 Sun–Mercury 的相对轨道，在日心惯性系中求球对称、非旋转太阳给出的一阶 GR 近心进动；其他行星与太阳扁率暂从模型中剔除。

坐标与正方向。 轨道平面内顺公转取正角方向；\(\Delta\varpi>0\) 表示近心点沿运动方向进动。

模型及理由。 \(GM_\odot/(ac^2)\ll1\) 且 \(v/c\ll1\)，故使用 Schwarzschild 弱场的一阶后 Newtonian 修正，而非完整数值相对论。

方程与计算。 每圈进动为

\[\Delta\varpi =\frac{6\pi GM_\odot}{a(1-e^2)c^2} =\frac{ 6\pi\left(1.3271\times10^{20}\ \mathrm{m^3\,s^{-2}}\right)} {\left(5.791\times10^{10}\ \mathrm m\right) \left[1-(0.2056)^2\right] \left(2.9979\times10^8\ \mathrm{m\,s^{-1}}\right)^2} =5.02\times10^{-7}\ \mathrm{rad/orbit} =0.1035\ \mathrm{arcsec/orbit}.\]

Mercury 每世纪约运行 \(415\) 圈，因此额外进动约

\[\Delta\varpi_{\mathrm{century}} =\left(0.1035\ \mathrm{arcsec/orbit}\right) \left(415\ \mathrm{orbit/century}\right) \approx43.0\ \mathrm{arcsec/century}.\]

同时

\[\epsilon_g =\frac{GM_\odot}{ac^2} =\frac{1.3271\times10^{20}\ \mathrm{m^3\,s^{-2}}} {\left(5.791\times10^{10}\ \mathrm m\right) \left(2.9979\times10^8\ \mathrm{m\,s^{-1}}\right)^2} =2.55\times10^{-8}\ll1.\]

SI、有效数字与检查。 进动式分子、分母都有 \(\mathrm{m^3/s^2}\)，故弧度无量纲；\(0.1035\times415=43.0\)，世纪值与每圈值一致。修正为正，且在 \(c\to\infty\) 或 \(a\to\infty\) 时趋零，极限正确。

物理解释。 这是弱场中微小却可测的一阶修正，不表示 Newtonian 引力在太阳系整体失效；与观测比较前仍须扣除经典多体和太阳扁率进动。

5. Lesson synthesis / 全课整合

5.1 Model-selection ladder

1. 尺寸远小于场变化尺度且只有两体：点质量二体模型。
2. 物体有限但 \(\lvert\boldsymbol\xi\rvert/R\ll1\)：潮汐张量线性化。
3. 接近 Roche 破坏：加入自引力、旋转、形变与材料本构；声明采用哪一个系数模型。
4. 研究整体束缚：按密度剖面积分自能，不能只背均匀球系数。
5. 三体或多体：写完整耦合方程；小作用可做摄动，长期或共振问题通常数值积分。
6. 强场、高速、光传播或精密相对论效应：转向 GR 或受控 post-Newtonian 模型。

5.2 Formula sheet

5.3 Model boundaries / 模型边界

• 潮汐线性式只保留 \(\lvert\boldsymbol\xi\rvert/R\) 的一阶项；接近源或跨越大尺度必须计算精确场差。
• Roche 系数取决于刚性、流体性、同步自转、密度结构与材料强度，不能脱离模型引用单一数值。
• \(N\) 体摄动和 Lagrange 点在指定参考系、初值与质量比下解释；混沌不破坏确定性或守恒律，却限制长期逐点预报。
• \(GM/(Rc^2)\ll1\)、\(v/c\ll1\) 只是 Newtonian 近似的必要尺度检查；所需测量精度可能更早迫使采用相对论模型。

5.4 What remains invariant across models

清晰的系统边界、参考系、尺度分析和守恒量仍是核心。扩展模型不是“再加一条神奇公式”，而是辨认先前被忽略的尺度：有限尺寸对应梯度，内部结构对应自能和材料响应，多体对应耦合与共振，高精度强场对应时空几何。

6. Misconception clinic / 易错概念诊所

7. Practice and retrieval / 练习与检索

A. Conceptual questions

1. 为什么自由落体实验室能消去共同重力，却不能消去潮汐？
2. 点质量外场潮汐张量的三个本征值为何和为零？
3. \(1.26\) 与 \(2.44\) 两个 Roche 系数为什么不能混用？
4. 引力自能公式中的负号和连续形式中的 \(1/2\) 各表示什么？
5. 为什么 \(N\) 体系统有守恒量仍不保证可积？
6. 混沌系统的短期高精度预测与长期统计描述是否矛盾？

B. Quantitative and derivation questions

7. 推导点质量场在径向和横向的潮汐本征值。
8. 若主天体密度不变、卫星密度变为原来的 8 倍，Roche 距离如何变化？
9. 均匀球半径在质量不变时收缩为一半，自能如何变化？释放多少能量，以初始 \(\lvert W_i\rvert\) 表示？
10. 写出行星摄动加速度中的 direct 与 indirect 两项，并解释第二项来源。
11. 某轨道有 \(GM/(rc^2)=10^{-9}\)、\(v/c=10^{-4}\)。它是否处于典型 Newtonian 区？
12. 太阳 Schwarzschild 半径约为多少？取 \(GM_\odot=1.3271\times10^{20}\ \mathrm{m^3/s^2}\)。

Exit ticket

为下列四种任务各选一个最低充分模型并写出判据：地球实验室两米尺度自由落体、接近 Roche 破坏的流体卫星、木星扰动下的小行星、Mercury 近日点的高精度进动。

8. References and source notes / 参考资料与使用说明

8.1 Source-to-content map

8.2 Bibliographic and modeling notes

OpenStax §§13.6–13.7 provide the undergraduate bridge from tides to relativity. NASA sources supply direct spacecraft, Roche, Lagrange and chaos context; MIT 8.901 independently supplies the self-energy derivation. Roche coefficients are explicitly separated by model rather than presented as one universal boundary. All prose, examples, equations and figures are original syntheses.