Gas drag on a planetesimal binary

For linear gas drag with aerodynamic stopping times $t_1,t_2$, masses $m_1,m_2$, total mass $m=m_1+m_2$, and relative velocity $\mathbf{v}$, defineIf $\mathbf{w}$ is the gas velocity relative to the centre of mass, the relative drag acceleration is $-\kappa \mathbf{v}+\Delta \mathbf{w}$. The first term dissipates internal specific orbital energy; the second is a differential headwind that vanishes for equal aerodynamic stopping times.

For a weakly perturbed circular binary planetesimal, the constant differential headwind does zero work averaged over the unperturbed orbit. Orbit-averaged drag work givesWith constant aerodynamic stopping times, $a(t)=a_0{e}^{-t/\tau }$. This averaging assumes drag is slow compared with the orbital period and does not appreciably distort the circular orbit within one period.

For a nearly circular stellar orbit of radius $a_b$ in gas moving at $(1-\eta )$ times the local circular speed, the averaged tangential gas drag on the centre of mass is $-\eta {\kappa }_b{v}_K$. Orbit-averaged drag work givesIf $m_2/m_1=ϵ\ll 1$ and aerodynamic stopping time scales as $m^k$, $0<k<1$, then $\tau /{\tau }_b\simeq \eta {ϵ}^k\ll 1$. Internal contraction can outpace stellar migration. Tightly gas-coupled drift requires solving radial and azimuthal motion together instead of this weak-drag approximation.