跳到主要内容
Sandbox Physics
从阳光推力读到真实彗尾

L3 · 彗尘动力学

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

Sandbox Physics ·

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

保留速度带来的一阶修正

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

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

r¨=−μ(1−β)rr3−βμcr2(vrr^+v)\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

用能量和角动量检查方向

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

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

ε=v22−μ′r,dεdt=−βμcr2(vr2+v2)≤0\varepsilon=\frac{v^2}{2}-\frac{\mu'}r,\qquad \frac{d\varepsilon}{dt}=-\frac{\beta\mu}{cr^2}(v_r^2+v^2)\le0
这里讨论力比固定时的粒子轨道能量；完整能量计算还包含辐射场。轨道功率为负，不表示每颗初始逃逸粒子都会在离开前重新被束缚。
h=r×v,dhdt=−βμcr2h\boldsymbol h=\boldsymbol r\times\boldsymbol v,\qquad \frac{d\boldsymbol h}{dt}=-\frac{\beta\mu}{cr^2}\boldsymbol h
此近似下，角动量方向不变、大小减小。若模型只有中心辐射压，角动量反而必须守恒。

对应阅读：Burns · Finson & Probstein

推导受束缚圆轨道的时钟

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

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

dadt=−2βμca,a(t)2=a02−4βμct\frac{da}{dt}=-\frac{2\beta\mu}{ca},\qquad a(t)^2=a_0^2-\frac{4\beta\mu}{c}t
圆轨道近似要求轨道尺寸变化远慢于一周公转。形式上到零半径的时间适合比较量级，不能当作穿过升华区和太阳表面的实际预测。
tPR,formal=ca024βμ=400.494β(a0AU)2yrt_{\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 力积分十年，检查轨道平均公式。

半径（微米）力比一天文单位圆轨道形式时间（年）
10.574697
30.1912090
100.05746970
1000.0057469700

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

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

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

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

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

亲手观察

分开读两只时钟

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

练习与答案

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

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

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

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

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

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

来源与计算范围

  1. Finson & Probstein, A Theory of Dust Comets I (1968)

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

  2. Burns, Lamy & Soter, Radiation Forces on Small Particles (1979)

    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

    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

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

  5. Vokrouhlický et al., The Yarkovsky and YORP Effects, arXiv:1502.01249v1

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

反馈这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策