从光的动量，推到尘粒感受到的减弱引力

先约定单位与正方向

把 $\boldsymbol r$ 定义为从太阳指向尘粒的位置向量，长度为 r，向外的径向单位向量为 $\hat{\boldsymbol r}$ 。太阳光度是每秒发出的能量，把它铺在半径为 r 的球面上，就得到尘粒处的能流密度。

本章先用米、千克、秒。太阳光度固定为三点八二八乘十的二十六次方瓦，太阳引力参数固定为一点三二七一二四四零零四一九三九三八乘十的二十次方立方米每平方秒。模拟改用天文单位和日，只改变数值单位。

S(r)=\frac{L_\odot}{4\pi r^2},\qquad \boldsymbol F_g=-\frac{GM_\odot m}{r^2}\hat{\boldsymbol r}

能流密度的单位是每平方米多少瓦。引力前的负号表示向内。光照和引力都按距离平方反比减弱，因此稍后求比值时，日心距离会消去。

光携带能量 E，也携带大小为 $E/c$ 的动量。半径为 s 的尘粒用横截面 $\pi s^2$ 截获阳光。把截获功率乘上无量纲压力效率 $Q_{\rm pr}$ ，再除以光速，就得到向外的力。这个效率同时概括吸收与散射造成的动量转移，并不等于吸光比例。

球形尘粒的质量等于密度乘体积。把质量代入，再让辐射力除以引力，横截面、日心距离和部分圆周率便消去了；分母最后留下一次方的粒子半径。

F_{\rm rad}=\frac{L_\odot Q_{\rm pr}s^2}{4cr^2},\qquad m=\frac{4\pi}{3}\rho s^3,\qquad \beta=\frac{3L_\odot Q_{\rm pr}}{16\pi cGM_\odot\rho s}

力比没有单位。它与日心距离无关，前提是粒子性质和按太阳光谱平均的压力效率保持不变。

取 $\rho=1000\,\mathrm{kg\,m^{-3}}$ 、 $Q_{\rm pr}=1$ 。三微米半径是百万分之三米，代入后得到力比约零点一九一。向外的辐射力抵消了约百分之十九的太阳引力，剩下约百分之八十一的向内加速度。

下表统一使用同一套太阳常数、球形密度与效率，直接代公式即可重算。半径加倍，力比减半；密度加倍也一样。真实压力效率会随尺寸和波长变化，因此不能把这张反比表一路延伸到分子尺度。

\beta=0.574236761\,Q_{\rm pr}\left(\frac{1000\,\mathrm{kg\,m^{-3}}}{\rho}\right)\left(\frac{1\,\mu\mathrm m}{s}\right)

系数来自前述太阳光度、引力参数与精确光速。这是在给定密度下的模型换算，不是对 ISON 尘粒密度的测量。

把两个径向力相加，就得到引力参数减小的开普勒问题。这个替换改变加速度，却不会重置释放时已有的速度。最清楚的第一步，是先把气体带来的释放速度设为零：粒子与彗核从同一点、带着同一速度出发，之后才因加速度不同而分开。

力比小于一，只保证简化模型里的中心加速度向内，不能保证轨道受束缚。轨道能量还取决于继承的速度。下一章会求出逃逸门槛，并说明常见的二分之一门槛为什么只适用于圆形母轨道。

\ddot{\boldsymbol r}=-GM_\odot(1-\beta)\frac{\boldsymbol r}{r^3}\equiv-\mu'\frac{\boldsymbol r}{r^3}

短时模型假定力比固定、尘粒之间无碰撞、气体加速已经结束，省略相对量级为速度除以光速的 PR 项。交互粒径都让力比小于一，因而可用有吸引力的椭圆、近抛物与双曲圆锥轨道传播。