---
title: "坡印廷—罗伯逊拖曳，怎样带走轨道能量"
description: "推导一阶拖曳、始终非正的轨道功率和圆轨道内迁定律，再与需要热滞后的反冲分开。"
canonical_url: "https://sandboxphysics.com/zh/stories/comet-dust/slow-loss-of-orbital-energy/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/comet-dust/slow-loss-of-orbital-energy/"
last_modified: "2026-10-28"
---

[从阳光推力读到真实彗尾](https://sandboxphysics.com/zh/stories/comet-dust/)

L3 · 彗尘动力学

# 坡印廷—罗伯逊拖曳，怎样带走轨道能量

Sandbox Physics · 2026-10-28

径向辐射压改变尘粒走哪一条圆锥轨道。还有一个小得多、取决于速度的力，能缓慢改变圆锥轨道本身。它含有轨道速度与光速之比，因此时钟和照片里可见的尘尾很不一样。

## 开始前需要什么

前两章的减弱引力、单位质量轨道能量和角动量。时间导数表示瞬时变化率。下文圆轨道参考要求粒子已按减弱引力下的圆轨道速度运动，并且受束缚。

## 保留速度带来的一阶修正

以日心速度 $\boldsymbol v$ 和径向分量 $v_r=\boldsymbol v\cdot\hat{\boldsymbol r}$ 描述尘粒。保留速度除以光速的一阶项，减弱引力方程外还要加一个与速度有关的修正。其中既有径向项，也有反向于完整速度的项；只随手画一根向后的箭头，会漏掉径向贡献。

此式采用固定压力效率，并按常见近似处理粒子静止系中的各向同性再辐射，不要求下午半球更热。吸收的太阳动量和再辐射都要在参考系之间一致变换；不能因为静止系里各向同性，就断言太阳系中没有相应动量损失。

$$
\ddot{\boldsymbol r}=-\mu(1-\beta)\frac{\boldsymbol r}{r^3}-\frac{\beta\mu}{cr^2}\left(v_r\hat{\boldsymbol r}+\boldsymbol v\right)
$$

第一项是短时尘尾采用的径向力，第二项才是 PR 修正，两项都是加速度。一天文单位圆轨道的速度与光速之比约为万分之一，因此修正远小于径向辐射压。

对应阅读：Burns · JPL Solar System Dynamics astronomical parameters

## 用能量和角动量检查方向

用减弱引力参数 $\mu'=\mu(1-\beta)$ 定义单位质量轨道能量，对动能和势能一起求导。中心力部分互相抵消，剩下速度与 PR 加速度的点积。最后，两个速度平方都乘着负系数。

力矩也给出一个直接检验：径向力对太阳不产生力矩，反速度项则削小轨道角动量。它们能检查积分器的物理方向，不必仅凭一幅螺旋图看起来是否可信。

$$
\varepsilon=\frac{v^2}{2}-\frac{\mu'}r,\qquad \frac{d\varepsilon}{dt}=-\frac{\beta\mu}{cr^2}(v_r^2+v^2)\le0
$$

这里讨论力比固定时的粒子轨道能量；完整能量计算还包含辐射场。轨道功率为负，不表示每颗初始逃逸粒子都会在离开前重新被束缚。

$$
\boldsymbol h=\boldsymbol r\times\boldsymbol v,\qquad \frac{d\boldsymbol h}{dt}=-\frac{\beta\mu}{cr^2}\boldsymbol h
$$

此近似下，角动量方向不变、大小减小。若模型只有中心辐射压，角动量反而必须守恒。

对应阅读：Burns · Finson & Probstein

## 推导受束缚圆轨道的时钟

轨道近圆且演化缓慢时，用半长轴代替距离，用减弱引力下的圆速度计算动能。单位质量能量为负的减弱引力参数除以两倍半长轴。对半长轴求导，再与刚才的轨道功率相等，减弱引力因子会消去，剩下一个简单的内迁速度。

将漂移方程乘以两倍半长轴再积分，得到轨道尺寸的平方随时间线性减小。它是受束缚轨道的平均近似；刚从母天体释放、还带着较快母体速度的粒子，未必满足这个初态。

$$
\frac{da}{dt}=-\frac{2\beta\mu}{ca},\qquad a(t)^2=a_0^2-\frac{4\beta\mu}{c}t
$$

圆轨道近似要求轨道尺寸变化远慢于一周公转。形式上到零半径的时间适合比较量级，不能当作穿过升华区和太阳表面的实际预测。

$$
t_{\rm PR,formal}=\frac{ca_0^2}{4\beta\mu}=\frac{400.494}{\beta}\left(\frac{a_0}{\mathrm{AU}}\right)^2\mathrm{yr}
$$

系数采用每年三百六十五点二五日与冻结太阳参数。初始轨道从一天文单位缩到半天文单位，参考时间就变成四分之一。

对应阅读：Burns · JPL Solar System Dynamics astronomical parameters

## 算算一百二十日能改变多少

十微米教学粒子的力比约零点零五七四，一天文单位圆轨道的形式时间约六千九百七十四年。在一百二十日内，轨道尺寸相对变化只有约百万分之二十四。同期形成的大尘尾主要来自释放后新轨道与释放速度，PR 内迁在这里缓慢得多。

下表与第一章使用相同密度、效率和圆轨道定律，可以独立计算。Lab 没有放大拖曳系数来挤进几个月的场景；独立图直接把横轴改为真实年。另用实际 PR 力积分十年，检查轨道平均公式。

| 半径（微米） | 力比 | 一天文单位圆轨道形式时间（年） |
| --- | --- | --- |
| 1 | 0.574 | 697 |
| 3 | 0.191 | 2090 |
| 10 | 0.0574 | 6970 |
| 100 | 0.00574 | 69700 |

对应阅读：Burns · JPL Solar System Dynamics astronomical parameters

## 热滞后需要再回答一个问题

旋转固体的下午侧可能比上午侧辐射更多热量。相应反冲取决于温度分布、自转方向、倾角、热导率与形状，这才进入昼夜雅可夫斯基效应；季节温度变化还会增加另一项。单凭力比无法决定这些作用。

交互场景因此只计算中心辐射压，并另列定义清楚的 PR 参考。原文导读再从热变量读到天体测量证据。对彗星，经验非引力加速度还可能描述不对称出气；拟合出的横向分量还需物理解释，分清出气与热反冲。

对应阅读：Vokrouhlický et al. · JPL Horizons API and observer/vector conventions

## 亲手观察

### [分开读两只时钟](https://sandboxphysics.com/zh/labs/comet-dust-tails/)

1. 对照运动的尘埃扇面与下方圆轨道参考，先看清 PR 横轴的年份再比较距离。
2. 把半径增大十倍。模型力比降十倍，圆轨道 PR 参考时间增十倍。

## 练习与答案

程序的 PR 修正在尘粒向外运动时增加了轨道能量，应该检查哪一条代数关系？

力比固定时，功率应为 $-\beta\mu(v_r^2+v^2)/(cr^2)$，径向速度取正取负都不应变成正值。检查径向符号、速度是否取自太阳参考系，以及势能是否采用同一个减弱引力参数。

十微米圆轨道参考，从一天文单位缩到半天文单位要多久？

尺寸平方定律需要形式时间的四分之三，约五千二百三十年。前提是粒子性质不变并一直存活；它不描述任意一颗从偏心彗星上逃逸的尘粒。

小行星存在横向轨道加速度，还需要什么信息才能支持昼夜热反冲解释？

自转方向、倾角、尺寸和热物性应能预测推力符号及量级，并由热模型解释温度滞后；还要比较可能的出气或其他辐射模型。系数描述加速度，其他观测才进一步检验成因。

[上一篇：继承的速度，怎样把光压变成尘尾](https://sandboxphysics.com/zh/stories/comet-dust/inherited-orbits-and-release-times/) [下一篇：把固定 JPL 尘埃几何投影到真实 HST 的 ISON 图像](https://sandboxphysics.com/zh/stories/comet-dust/ison-image-projection-tutorial/)

## 来源与计算范围

1. [Finson & Probstein, A Theory of Dust Comets I (1968)](https://adsabs.harvard.edu/pdf/1968ApJ...154..327F)

   ApJ 154，327—352，DOI 10.1086/149761。第 I—III 节，特别是第 329、334—335 页与图 4：脱离气体拖曳后的尘粒、减弱的太阳引力、等时线与等粒径线。教学场景另加明确写出的释放速度假设。

2. [Burns, Lamy & Soter, Radiation Forces on Small Particles (1979)](https://ntrs.nasa.gov/citations/19790070469)

   Icarus 40，1—48，DOI 10.1016/0019-1035(79)90050-2。第 III—IV 节推导一阶辐射力；第 VI 节式 (19) 与图 7 讨论压力效率；第 VII—VIII 节、式 (47)—(50) 讨论减弱引力的轨道与轨道平均 PR 演化。核读作者托管的 48 页全文。

3. [JPL Horizons API and observer/vector conventions](https://ssd.jpl.nasa.gov/horizons/manual.html)

   API 响应版本 1.2，API 文档版本 1.3。固定彗星记录 90004196、解 JPL#57（2025-Nov-21\_13:42:48）、DE441。对应 TIME\_TYPE、TLIST、ICRF/FRAME 向量约定和观测量 1、19、20、21、24、27、30；保存的请求与响应保留时标和几何选择。

4. [JPL Solar System Dynamics astronomical parameters](https://ssd.jpl.nasa.gov/astro_par.html)

   天文单位、日、光速以及相对力检验所用行星系统引力参数。太阳引力参数 132712440041.93938 立方千米每平方秒、光度三点八二八乘十的二十六次方瓦固定自 Horizons 太阳元数据。教学密度与压力效率是设定，不是 ISON 尘粒实测值。

5. [Vokrouhlický et al., The Yarkovsky and YORP Effects, arXiv:1502.01249v1](https://arxiv.org/abs/1502.01249v1)

   固定 2015 年 v1，共 24 页。第 2.1 节式 (1)—(3) 给出昼夜与季节性半长轴漂移；第 2.2 节讨论热物性和形状；第 3.1 节及图 2—3 讨论 Golevka、1992 BF 等天体测量探测。本页导读原文，不添加新热模拟，也不声称汇总最新探测。