# Dynamical Friction in a Gaseous Medium

## 1 Introduction

The process of dynamical friction, defined as momentum loss by a massive moving object due to its gravitational interaction with its own gravitationally-induced wake, arises in many astronomical systems. Examples of systems in which such effects are well known to be important include stars in clusters or galaxies, galaxies in galaxy clusters, and binary star cores in the common envelope phase of evolution. In the first two examples, the surrounding background medium in general consists of a combination of collisionless matter (stars, galaxies, dark matter) and gas, while in the third example the surrounding medium is entirely gaseous. As a corollary to the dynamical friction process, in all these cases the background medium is heated at an equal and opposite rate to the energy lost by the perturber.

The analytic theory for the gravitational drag in collisionless systems was developed by Chandrasekhar (1943), and over the decades since has enjoyed widespread theoretical application, extensive verification by numerical experiments, and well-documented embodiment in observed astronomical systems. The variety of important consequences of gravitational drag in collisionless astronomical systems includes mass segregation in star clusters, sinking satellites in dark matter galaxy halos, orbital decay of binary supermassive black holes after galaxy mergers, etc.; see e.g. Binney & Tremaine (1987).

Less well-developed is the corresponding theory of dynamical friction in a gaseous (i.e., collisional) medium.<sup>1</sup><sup>1</sup>1Instead, most studies of the gravitational interaction between a moving massive body and the surrounding gaseous medium have focused on the problem of accretion, following on the analysis of Hoyle and Lyttleton (1939) and Bondi and Hoyle (1944), and the early numerical work of Hunt (1971). For supersonic motion, analytic linear-theory estimates of the gravitational drag under assumption of a steady state were obtained by Dokuchaev (1964), Ruderman & Spiegel (1971), and Rephaeli & Salpeter (1980). The resulting estimates for the drag force in the steady supersonic case take the form

$$F_{SS}=\frac{4\pi (GM_p)^2\rho _0}{V^2}\mathrm{ln}\left(\frac{r_{max}}{r_{min}}\right)$$
(1)   
where $`r_{max}`$ and $`r_{min}`$ correspond respectively to the effective linear sizes of the surrounding medium and the perturbing object, similar to the drag formula obtained for a collisionless medium. Although there is some ambiguity in the definition of $`r_{max}`$ and $`r_{min}`$, the estimate (1) appears consistent with calculations of the gravitational drag that are obtained as a by-product of numerical hydrodynamic investigations focused on the Bondi-Hoyle-Lyttleton accretion problem – see, e.g. Shima et al (1985), Shankar et al (1993), and Ruffert and collaborators (see Ruffert (1996) and references therein).

For steady-state, subsonic motion of the perturber, the front-back symmetry of the perturbed density distribution about the perturber led Rephaeli & Salpeter (1980) to argue that gravitational drag is absent in the subsonic, inviscid case. Considering that the drag force for the supersonic case increases proportional to $`V^2`$ with decreasing perturber speed, $`V`$, it seems counterintuitive for the dynamical drag to become exactly zero when $`V`$ becomes infinitesimally smaller than the sound speed. In this paper, we reconsider the linear-theory drag as a time-dependent rather than steady-state problem, and arrive instead at a nonzero value for the dynamical friction for the subsonic case, while still verifying that the drag force is maximized for perturbers with $`Vc_s`$. In §2, we derive results for the perturbed density distributions created by, and dynamical drag on, a massive perturber on a constant-velocity trajectory through a uniform, infinite medium. In §3 we relate our results to the classical (collisionless) dynamical friction formula, and briefly consider applications to a few astronomical systems.

## 2 Analysis

### 2.1 Wave Equation and Formal Solution using Green Functions

We begin with the linearized equations for the perturbed density $`\rho \rho _0[1+\alpha (𝐱,t)]`$ and velocity $`𝐯c_s\beta (𝐱,t)`$ of an adiabatic gaseous medium that is subject to an external gravitational potential $`\mathrm{\Phi }_{ext}(𝐱,t)`$:

$$\frac{1}{c_s}\frac{\alpha }{t}+\text{}\beta =0$$
(2)   
and

$$\frac{1}{c_s}\frac{\beta }{t}+\text{}\alpha =\frac{1}{c_s^2}\text{}\mathrm{\Phi }_{ext},$$
(3)   
where $`\rho _0`$ is the unperturbed density, $`c_s(P_0/\rho _0)^{1/2}`$ is the sound speed, and the perturbation amplitudes $`\alpha ,|\beta |<<1`$. By substituting equation (2) in the divergence of equation (3), we have

$$^2\alpha \frac{1}{c_s^2}\frac{^2\alpha }{t^2}=\frac{1}{c_s^2}^2\mathrm{\Phi }_{ext}4\pi f(𝐱,t)$$
(4)   
where $`\rho _{ext}(𝐱,t)c_s^2f(𝐱,t)`$ is the mass density of the perturber.

In the absence of any disturbance prior to the action of $`\mathrm{\Phi }_{ext}`$, the solution to equation (4) is found using the retarded Green function $`G^+`$ for the 3D wave equation (Jackson (1975) §6.6) as

$$\alpha (𝐱,t)=d^3x^{}𝑑t^{}\frac{\delta [t^{}(t|𝐱𝐱^{}|/c_s)]f(𝐱^{},t^{})}{|𝐱𝐱^{}|}$$
(5)   
### 2.2 Perturbed Density Distributions for Constant-Velocity Perturbers

We now specialize to the case of a point mass $`M_p`$ on a straight-line trajectory with velocity $`V\widehat{z}`$, passing at time $`t=0`$ through $`𝐱=0`$. If $`(t)`$ describes the time over which the perturber is active, then $`f(𝐱,t)=(GM_p/c_s^2)\delta (zVt)\delta (x)\delta (y)(t)`$. We define $`szVt`$ as the distance along the line of motion relative to the perturber, $`wz^{}z`$, and $`V/c_s`$ as the Mach number; $`w>0(<0)`$ corresponds to the perturbation from a backward (forward) - propagating wave. Then

$$\alpha (𝐱,t)=\frac{GM_p}{c_s^2}_{\mathrm{}}^{\mathrm{}}𝑑w\frac{\delta [w+s+(R^2+w^2)^{1/2}]((w+z)/V)}{(R^2+w^2)^{1/2}}.$$
(6)   
Here $`R=(x^2+y^2)^{1/2}`$ is the cylindrical radius.

To evaluate the integral (6), one expands the argument of the $`\delta `$ function about its possible roots

$$w_\pm =\frac{s\pm [s^2+R^2(1^2)]^{1/2}}{^21};$$
(7)   
for $`<1`$, only $`w_+`$ is a valid root, whereas for $`>1`$ both roots are valid as long as $`s<0`$ and $`|s|/R>(^21)^{1/2}`$, and neither otherwise. Using $`\delta ((ww_\pm )A)=\delta (ww_\pm )/|A|`$ for $`A=1+w_\pm /(R^2+w_\pm ^2)^{1/2}`$ and substituting the solution (7), the result is

$$\alpha =\frac{GM_p/c_s^2}{[s^2+R^2(1^2)]^{1/2}}\underset{\mathrm{roots}\mathrm{w}_0}{}((z+w_0)/V).$$
(8)   
For a purely steady solution over all time, $``$ is unity for all arguments, with the result that

$$\alpha _S=\frac{GM_p/c_s^2}{[s^2+R^2(1^2)]^{1/2}}\times \{\begin{array}{cc}1\hfill & \text{if }<1\hfill \\ 2\hfill & \text{if }>1\text{ and }s/R<(^21)^{1/2}\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}$$
(9)   
This confirms the analytic, linear-theory results of previous authors for completely steady flow created by a point mass on a straight-line, constant-speed trajectory: (a) A subsonic perturber generates a density distribution centered on the perturber at $`s=0`$, with contours of constant density corresponding to similar ellipses in the $`sR`$ plane with eccentricity $`e=`$, and the short axis along the line of motion of the perturber (i.e. the 3D density distribution consists of concentric similar oblate spheroids). The density contrast $`\rho _1/\rho _0\alpha `$ is unity for an elliptical section with semiminor axis $`GM_p/c_s^2`$; the linearization is sensible only outside of this elliptical section. This density distribution is a generalization of the far-field limit of the hydrostatic envelope $`\rho /\rho _0=exp(GM_p/(c_s^2r))`$ that surrounds a stationary perturber.<sup>2</sup><sup>2</sup>2We note that for subsonic flow, the far-field density enhancement from equation (9) is within 10% of Hunt’s (1971) $`=0.6`$ numerical solution outside of $`r=GM_p/c_s^2`$. Also, the rms anisotropy of 7% from the far-field steady solution (9) with $`=0.6`$ is comparable to the 5% mean anisotropy cited for Hunt’s near-field numerical solution. (b) A supersonic perturber generates a density wake only within the rear Mach cone of half-opening angle $`\mathrm{sin}\theta =1/`$ defined by $`s/R<(^21)^{1/2}`$; the surfaces of constant density within the wake correspond to hyperbolae in the $`sR`$ plane, with eccentricity $`e=`$.

Now consider the case where the perturber is “turned on” at $`t=0`$, so that $``$ is a Heaviside function. For a subsonic perturber $`<1`$ the only root is $`w_+`$ (eq. 7). Algebraic calculation shows that $`z+w_+>0`$ when $`R^2+z^2<(c_st)^2`$, for $`<1`$. Thus, the region of perturbation is the sphere centered on the original position of the perturber, within which a sound wave has traveled in time $`t`$. Within this region of perturbation, the density distribution has reached the value given by the steady solution (eq. 9); outside the causal region for sound waves, the density remains unperturbed.

For a supersonic perturber, any density disturbance must be confined within the rear Mach cone $`s/R<(^21)^{1/2}`$. Algebraic calculation shows that within the sphere $`R^2+z^2<(c_st)^2`$, $`z+w_+>0`$ and $`z+w_{}<0`$; hence only $`w_+`$ contributes in equation (8). Within the Mach cone and to the right ($`R<|zc_st|/(^21)^{1/2}`$, $`z>c_st/`$) of this sphere ($`R^2+z^2>(c_st)^2`$), both $`z+w_+`$ and $`z+w_{}`$ are real and positive, hence contribute in eq. (8).

The results for the perturbed density for this finite-time perturbation are summarized as follows: <sup>3</sup><sup>3</sup>3These results appear to have been obtained previously by Just & Kegel (1990) via an alternative mathematical formalism; the two calculations serve as independent checks of the formulae.

$$\alpha (t)=\frac{GM_p/c_s^2}{[s^2+R^2(1^2)]^{1/2}}\times \{\begin{array}{cc}1\hfill & \text{if }R^2+z^2<(c_st)^2\hfill \\ 2\hfill & \text{if }>1\text{}R^2+z^2>(c_st)^2\text{}s/R<(^21)^{1/2}\text{, and }z>c_st/\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}$$
(10)   
The region of perturbed density has the shape of a loaded ice-cream-cone dragged by its point by the perturber $`M_p`$; only $`w_+`$ contributes in the ice-cream region (“region 1”), while both $`w_\pm `$ contribute in the cone region (“region 2”). The cone shrinks in size as $``$ decreases for $`>1`$, and is nonexistant for $`<1`$. Figures 1 and 2 show examples of the perturbed density distributions for subsonic and supersonic perturbers, respectively. Because of the linear-theory assumptions made at the outset, we note that equation (10) is properly valid only for $`\alpha <<1`$, i.e. $`s^2+R^2(1^2)>>(GM_p/c_s^2)^2`$.

### 2.3 Gravitational Drag Formulae

To compute the gravitational drag on the perturber, we need to evaluate the gravitational force between the perturber and its wake,

$$F_{df}=2\pi GM_p\rho _0𝑑s𝑑RR\frac{\alpha (t)s}{(s^2+R^2)^{3/2}}.$$
(11)   
To perform the volume integral, it is convenient to transform to a spherical polar coordinate system ($`r,\theta `$) centered on the massive perturber, with $`Rr\mathrm{sin}\theta `$ and $`sr\mathrm{cos}\theta `$. Defining $`\mu =\mathrm{cos}\theta `$, $`xr/c_st`$, we have

$$F_{df}=I,\frac{4\pi (GM_p)^2\rho _0}{V^2},$$
(12)   
where

$$I=\frac{1}{2}\frac{dx}{x}𝑑\mu \frac{\mu ^2𝒮_{}}{(1^2+\mu ^2^2)^{1/2}}.$$
(13)   
Here $`𝒮_{}`$ represents the sum in equation (8).

For the purely steady-state density perturbation given in equation (9), in the subsonic case $`𝒮_{}=1`$ everywhere in space and front-back antisymmetry in the angle integral (due to the symmetric spheroidal density distribution) argues that there is zero net force on the perturber. For the steady-state, supersonic case, $`𝒮_{}=2`$ for all $`\mu `$ between $`1`$ and $`\mu _M\sqrt{^21}/`$ (the boundary of the Mach cone), so that $`I=𝑑x/x\mathrm{ln}(r_{max}/r_{min})`$. This is consistent with previous results, and yields an identical formula to that representing the dynamical friction force in a collisionless medium when $`V`$ is much larger than the background particle velocity dispersion (e.g. Binney & Tremaine (1987), eq. 7-18).

For the finite-time case, based on the perturbed density distribution in equation (10), there is nonzero contribution to the integral $`I`$ only from a finite region in space: region 1, in which $`𝒮_{}=1`$, and region 2, in which $`𝒮_{}=2`$. Region 2 is nonexistant for subsonic perturbers. For the subsonic case, the perturber is surrounded by a concentric distribution of similar ellipsoids, and is displaced by $`Vt`$ forward from the center of the sonic sphere (radius $`c_st`$) surrounding its initial position. The nearby, complete ellipsoids exert no net force on the perturber, but the larger, cut-off ones with semiminor axes between $`(c_sV)t`$ and $`(c_s+V)t`$ lag behind the perturber (see figure 1) and exert a gravitational drag. Thus, the radial integral in equation (13) has upper/lower limits $`x=1\pm `$, and the angular integral has limits $`\mu =1,\mu _C`$ with $`\mu _C(1^2x^2)/(2x)`$. The integrals are straightforward; the result is

$$I_{subsonic}=\frac{1}{2}\mathrm{ln}\left(\frac{1+}{1}\right).$$
(14)   
The implicit assumptions in deriving this equation are that $`(c_sV)t`$ exceeds the effective size of the perturber ($`r_{min}`$), and that $`(c_s+V)t`$ is smaller than the effective size of the surrounding gaseous medium ($`r_{max}`$). Under these conditions, the dynamical drag is time-independent, and non-zero. The steady-state result that zero net force results from the front-back symmetry of the density distribution is misleading; because of the long-range nature of the Coulomb potential, the total drag force at any time depends on the unchanging ratio $`(c_st+vt)/(c_stvt)=(1+)/(1)`$ between the semiminor axes of the furthest and closest perturbing partial spheroids. The gravitational drag is always dominated by the far-field, and at any time the perturber is located ahead of center of the sonic sphere. Physically, we can associate the energy loss arising from the drag force with the rate at which the expanding sound wave does work on the background medium. In the limit of a very slow perturber $`<<1`$, $`I_{subsonic}^3/3`$, so that the drag force is proportional to the perturber’s speed $`V`$.

For the supersonic case, the whole of the perturbed density distribution lags the perturber. The angular integration limits are $`\mu =1,\mu _{}`$ for $`x=r_{min}/(c_st)`$ to $`\sqrt{^21}`$, and $`\mu =1,\mu _C`$ for $`x=\sqrt{^21}`$ to $`+1`$; $`𝒮_{}`$ takes on values 2 and 1 in regions 2, 1. The result of the integration is

$$I_{supersonic}=\frac{1}{2}\mathrm{ln}\left(\frac{+1}{1}\right)+\mathrm{ln}\left(\frac{1}{r_{min}/c_st}\right)=\frac{1}{2}\mathrm{ln}\left(1\frac{1}{^2}\right)+\mathrm{ln}\left(\frac{Vt}{r_{min}}\right).$$
(15)   
We have assumed that $`Vtc_st>r_{min}`$, and that the effective size of the background medium exceeds $`Vt+c_st`$ . In the limit $`>>1`$, we have $`I_{supersonic}\mathrm{ln}(Vt/r_{min})`$; with $`Vtr_{max}`$, this recovers the steady-state result.<sup>4</sup><sup>4</sup>4As pointed out by the referee Scott Tremaine, the notion that the maximum impact parameter increases as $`Vt`$ was earlier introduced by J. Ostriker & Davidsen (1968) in a time-dependent analysis of stellar relaxation.

In Figure 3, we plot the dynamical friction force as a function of the Mach number for several values of $`c_st/r_{min}`$.

## 3 Discussion

The main formal result of this paper is the evaluation, in linear perturbation theory, of the gravitational drag force $`F_{df}`$ on a massive perturber $`M_p`$ moving on a straight-line trajectory through an infinite, homogeneous, gaseous medium of density $`\rho _0`$ and sound speed $`c_s`$. Together with the definition in equation (12), the expressions (14) and (15) give the dynamical friction (df) drag forces on subsonic and supersonic perturbers, respectively. Figure 3 presents the same results graphically, showing how the drag force varies with the Mach number $``$ of the perturber, and the time over which the perturber has been moving with fixed speed $`V=c_s`$.

For comparison, we have also included in Figure 3 the result for the gravitational drag on a particle of mass $`M_p`$ moving through a collisionless medium with the same density $`\rho _0`$ as the gaseous medium we have considered, and with a Maxwellian distribution of particle velocities with $`\sigma =c_s`$. From equation (7-18) of Binney & Tremaine (1987), the collisionless df drag force is given by expression (12) with

$$I_{collisionless}=\mathrm{ln}\left(\frac{r_{max}}{r_{min}}\right)[\mathrm{erf}(X)\frac{2X}{\sqrt{\pi }}e^{X^2}],$$
(16)   
where $`XV/(\sigma \sqrt{2})`$. From Figure 3, it is clear that (a) for $`>>1`$, the collisionless and gaseous df forces are identical (as has been previously noted); (b) for $`<1`$, the drag force is generally larger in a collisionless medium than in a gaseous medium, because in the latter case pressure forces create a symmetric distribution in the background medium in the vicinity of the perturber; (c) the functional form of the gaseous df drag is much more sharply peaked near $`=1`$ than it is for the collisionless df drag; perturbers moving at speeds near Mach 1 resonantly interact with the pressure waves that they launch in the background medium; (d) for a given value of $`\mathrm{ln}(\mathrm{\Lambda })\mathrm{ln}(r_{max}/r_{min})`$, the peak value of the gaseous df force is much larger than the corresponding peak value of the collisionless df force; at $`1`$, there is a factor of four difference in the force between the two cases. Below, we explore some potential consequences of these results in a variety of astronomical systems.

As a consequence of the stronger gaseous df force than collisionless df drag force for supersonic motion, massive objects may make their way more rapidly to the center of a star cluster or galaxy if they arrive at the outer edge before, rather than after, the gas is turned into stars. For a particular example, we consider the decay of a massive perturber’s near-circular orbit in a spherical density distribution with a singular isothermal sphere profile, $`\rho (r)=c_s^2/(2\pi Gr^2)`$ so $`M(r)=2c_s^2r/G`$ (here $`c_s`$ respectively denotes the sound speed or velocity dispersion for a gaseous or stellar distribution). For this density profile, the circular speed is constant, $`V=c_s\sqrt{2}`$ . By equating the rate of decrease of angular momentum $`d(M_pVr)/dt`$ to the torque $`\tau _{df}=rF_{df}`$, one finds that the time for the perturber’s orbit to decay from $`r_{init}`$ to $`r<<r_{init}`$ is given approximately by

$$\frac{t_{df}(r_{init})}{t_{orb}(r_{init})}\frac{M(r_{init})}{4\pi M_p\mathrm{ln}\left(\frac{r}{r_{min}}\right)}$$
(17)   
for gaseous df. This time is a factor of $`0.428`$ shorter than the corresponding decay time under stellar df (cf. Binney & Tremaine (1987), eq. 7-25). For galactic applications, the implication is that condensed objects that form early (e.g. globular clusters) could spiral into the galactic center from a factor 1.5 further in a galaxy than would be predicted by stellar df theory (Tremaine, Ostriker, & Spitzer (1975)), within the epoch over which the baryonic distribution remains gaseous. For star cluster applications, the shorter df time for gaseous distributions may help explain why young, embedded stellar clusters like the Orion Nebular Cluster (Hillenbrand (1997), Hillenbrand & Hartmann (1998)) appear significantly more relaxed than expected from stellar df alone; e.g., for the ONC, n-body simulations show that stellar df would require a time a factor 3-4 longer than the best estimate of the cluster age to establish the observed mass segregation (Marshall, Ostriker, & Teuben (1998)).

As mentioned above (see Fig. 3), the gaseous df force is strongly depressed for subsonic perturbers. For modeling the global evolution of combined star-gas systems in which the particle velocity dispersion and gas sound speed are comparable, the strong cutoff of $`F_{df}`$ for subsonic perturbers implies that setting $`F_{df}=0`$ for $`V<c_s`$ (e.g. Saiyadpour, Deiss, & Kegel (1997)) should yield realistic results. However, in other circumstances it is interesting to enquire how the small, but nonzero, df drag on subsonic perturbers (computed in this paper) can affect their orbit evolution.

As a model problem, we consider the decay of a subsonically-moving mass $`M_p`$ on a near-circular orbit embedded within a constant-density gaseous sphere of radius $`r_0`$. For constant background density, the angular orbit frequency $`\mathrm{\Omega }=\sqrt{4\pi G\rho _0/3}`$ and orbital period $`t_{orb}=2\pi /\mathrm{\Omega }`$ are independent of the distance $`r`$ from the center. If this constant-density region represents the core of a nonsingular isothermal sphere with core radius $`r_0`$ and sound speed $`c_s`$, then $`r_0\mathrm{\Omega }/c_s=\sqrt{3}`$; thus, the Mach number for a circular orbit at $`r`$ is $`=\sqrt{3}r/r_0`$. Just as above for decay in a power-law density profile, we can compute the time for the perturber’s orbit to decay from $`r_{init}=r_0/\sqrt{3}`$ to $`r_f`$ as

$$\frac{t(r_f)}{t_{orb}}=\frac{M_{core}}{\pi 3^{5/2}M_p}_{M_f=\frac{\sqrt{3}r_f}{r_0}}^1\frac{d}{I_{subsonic}()};$$
(18)   
the numerical coefficient equals $`0.0204`$.

Figure 4a shows how the decay of the orbital radius depends on time and on the mass of the perturber relative to the whole core. In Figure 4b, we verify that the assumption of a near-circular orbit is valid provided $`M_p/M_{core}<<1`$, since

$$\frac{v_r}{v_\phi }=\frac{3^{5/2}M_p}{2M_{core}}\frac{I_{subsonic}()}{^3},$$
(19)   
which has the limit $`2.598M_p/M_{core}`$ for $`<<1`$. In arriving at the results shown in Fig. 4, we have used equation (14). Its validity depends, however, on the size of the uniform-density core exceeding that of the perturbed-density region. This assumption must fail, and the df drag force consequently decrease, when the forward wave defining the disturbed-density region (cf. Fig. 1) reaches a distance $`r_0`$ ahead of the perturber.<sup>5</sup><sup>5</sup>5Because the perturber follows a circular rather than straight-line orbit, the df drag at late times should still be nonzero: Since the direction of $`\widehat{V}`$ changes by $`\pi /2`$ four times per orbit, the forward-wave propagation effectively “restarts” as well. When $`r<<r_0`$, over each quarter orbit, each new expansion wave propagates to just $`r_0`$; over much of the quarter-orbit, an unbalanced trailing density enhancement will remain within the core. Thus, we expect that the df drag will still be some fraction of the value found using eq. (14), although a more refined calculation is needed to predict the fraction quantitatively. The position $`r(t_{fw})`$ of the perturber at this time is plotted as a function of $`M_p/M_{core}`$ in Figure 5. Based on this figure, only perturbers of mass $`M_p>0.2M_{core}`$ are predicted to reach within one-tenth of a core radius before the df drag decreases. For lower-mass perturbers, this implies that the decay of orbits to very small radii may stall, if the df drag is sharply reduced after the time $`t_{fw}`$. This result may have relevance for models of QSO evolution (Silk & Rees (1998)) in which primordial black holes are formed away from the centers of galaxies, later to be driven there by df during mergers. If, as assumed in this scenario, these events occur before the advent of star formation, than the relevant df drag is the gaseous df examined in this paper. If the df drag becomes inefficient when the orbit reaches the protogalaxy core and becomes subsonic, then these massive black holes may have more time to grow by accretion before they finally sink to the centers of their host galaxies.

Finally, we comment on the applicability of our results to the interaction of protoplanets with a gaseous circumstellar nebula in which they may grow. Protoplanets orbit faster than the surrounding gaseous nebula, with the Mach number of the drift speed $`v_{rel}/c_sc_s/(\mathrm{\Omega }_Kr)<<1`$ (where $`\mathrm{\Omega }_K`$ is the local Kepler angular rotation frequency and $`r`$ is the protoplanet’s semimajor axis). Drag forces occur as a result of this relative motion, leading to inward radial migration of protoplanets as they lose angular momentum to the surrounding nebula (e.g. Ward (1997) and references therein). Using parameters based on solid bodies in the protosolar nebula, it can be shown that the nominal df drag from equations (12) and (14), $`F_{df}(4\pi /3)(GM_p)^2\rho _0/(c_s\mathrm{\Omega }_Kr)(2\pi /3)(GM_p)^2\mathrm{\Sigma }/(c_s^2r)`$, exceeds the large-Reynolds number aerodynamic drag $`F_{aer}(1/4)v_{rel}^2\pi R_p^2\rho _0`$ (Landau & Lifshitz (1987), §45) for protoplanets of radius $`R_p`$ greater than a few $`100\mathrm{k}\mathrm{m}`$.

It is sometime argued that gravitationally-enhanced drag on protoplanets (cf. Takeda et al (1985), Ohtsuki, Nakagawa, & Nakazawa (1988), but note the different scaling from our results) may enhance migration rates over those predicted to arise from differential resonant torques (Goldreich & Tremaine (1980)) between the protoplanet and the surrounding nebula. In fact, provided that the difference of gas orbits from Kepler orbits is included so that the perturber does not corotate with the nebular gas at the same radial distance, then the df drag is automatically incorporated in the net resonant torque. Indeed, calculations show that the net resonant torque from a thin disk (e.g. Korycansky & Pollack (1993), Ward (1997)), taking this velocity difference into account, is of the same magnitude (and scaling) as the df drag estimated above. Thus, while the gaseous df drag may have implications for the late stages of planet formation, its effects are naturally incorporated within a full resonant formalism. Although such a calculation has not yet been performed in three dimensions, it is reassuring that our simple estimate of the drag – which neglects gradients of velocity, density, and temperature in the disk – and the estimates from 2D resonant torque differences – which neglect the disk thickness – nevertheless yield comparable answers.

The author is grateful to the referee Scott Tremaine for helpful comments on the manuscript.