Pith. Sign in

REVIEW

Dynamical Friction in a Gaseous Medium

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Eve C. Ostriker

Astrophysics · Submitted 20 Oct 1998 · 24 references

Pith's verdict: A clean, checkable linear-theory calculation that settles the subsonic gaseous drag question; the soft spots are the turn-on bookkeeping and the order-of-magnitude applications, both of which the paper mostly owns.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 1 major objection · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/9810324 (version 1)
Pith IDPTBTMIIF, Pith's permanent record of this paper and its review
CitationEve C. Ostriker. "Dynamical Friction in a Gaseous Medium." arXiv:astro-ph/9810324 (1998). https://pith.science/paper/PTBTMIIF
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

A massive object moving through gas feels gravitational drag even below the sound speed, and the drag is strongest near Mach 1.

The reading

The paper asks whether a massive object moving through a uniform gas loses momentum by gravity even when it travels slower than the sound speed. Using time-dependent linear perturbation theory, the paper computes the wake and finds a nonzero subsonic drag, $I_{subsonic}=1/2 ln[(1+M)/(1-M)]-M$, which falls to $M^3/3$ as $M$ goes to zero; the supersonic result recovers the familiar $ln(Vt/r_{min})$ law. The point matters because the drag sets how fast gas-rich protogalaxies, proto-star clusters, and protoplanetary disks can pull massive bodies inward, and because it contradicts the older steady-state argument that subsonic inviscid gas produces no drag. The paper also compares the gaseous force with the collisionless formula and finds the gas is more efficient above Mach 1 but less efficient below it.

What carries the argument

The machinery is the retarded Green function for the three-dimensional wave equation sourced by the perturber's gravitational potential. It turns the linearized continuity and momentum equations into an integral for the fractional density perturbation $alpha(x,t)$, whose roots and Heaviside cutoffs define the causal wake: a sphere of radius $c_s t$ for $M<1$, and a Mach cone plus a sonic sphere for $M>1$. The drag integral then reduces to an angular and radial integral over that wake, with the finite-time cutoff supplying the asymmetry that a purely steady calculation misses. The same Green-function bookkeeping yields the closed forms in equations (14) and (15) and makes the comparison with the collisionless formula direct.

Extended reading notes

Core claim

The central claim is that gravitational dynamical friction in a gaseous medium is a finite-time causal effect, not a steady-state symmetry. For a point mass $M_p$ moving at constant speed $V$ through gas of density $rho_0$ and sound speed $c_s$, the drag is $F_{df}=-I 4pi(GM_p)^2rho_0/V^2$. In the subsonic case the wake is a set of oblate spheroids inside the sonic sphere; the near, complete spheroids pull symmetrically, but the larger ones cut off behind the perturber leave a net backward pull, giving $I_{subsonic}=1/2 ln[(1+M)/(1-M)]-M$, which is nonzero for all $0<M<1$. In the supersonic case the wake is confined to the trailing Mach cone and $I_{supersonic}=1/2 ln(1-1/M^2)+ln(Vt/r_{min})$, recovering $ln(Vt/r_{min})$ for $M$ much greater than 1. The result overturns the steady-state conclusion that subsonic inviscid drag vanishes, and it shows the drag peaks near $M=1$.

Load-bearing premise

The subsonic result treats the perturber as switched on at $t=0$ and then assumes the drag becomes time-independent once $(c_s-V)t$ exceeds the body size and $(c_s+V)t$ stays inside the medium, so the force rests on the ratio of those cutoffs rather than on a demonstrated history-independent steady state.

Editorial extensions

If this is right

  • Subsonic perturbers do sink, but slowly: setting $F_{df}=0$ for $V<c_s$ is a reasonable approximation only when the drag is negligible compared with other orbital forces.
  • For supersonic motion in an isothermal $rho$ proportional to $r^{-2}$ halo, the gaseous decay time is about 0.428 times the stellar-dynamical-friction time, so gas-rich systems can sink objects from farther out in the same epoch.
  • The drag is sharply peaked near $M=1$, so perturbers with speeds close to the sound speed lose angular momentum fastest.
  • In a constant-density core, the forward sound wave eventually leaves the core and the drag drops; the paper estimates that only perturbers with $M_p$ greater than about $0.2M_{core}$ reach within one-tenth of a core radius before that happens.
  • For protoplanets, the gravitational drag exceeds ordinary aerodynamic drag for bodies larger than a few hundred kilometers, although in a full disk calculation it is part of the resonant torque.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The subsonic formula is tied to a switch-on at $t=0$; a perturber with a different formation history, or one moving indefinitely through a finite medium, could in principle show a different force, which makes history-dependent drag a testable prediction rather than a nuisance.
  • Because the effect comes from unequal inner and outer causal cutoffs, any process that truncates the wake, such as finite core size, stratification, or orbital motion, should modify the drag in a calculable way, suggesting extensions to disks and cores.
  • The Mach-1 peak and the $M^3$ subsonic tail are clean targets for moving-mesh or grid hydrodynamics of a gravitating perturber in a uniform gas.
  • If three-dimensional resonant-torque calculations reproduce the same $Sigma/(c_s^2 r)$ scaling for protoplanets, that would tie the local gravitational-wake picture to the global disk-torque picture.
Original abstract

Using time-dependent linear perturbation theory, we evaluate the dynamical friction force on a massive perturber M_p traveling at velocity V through a uniform gaseous medium of density rho_0 and sound speed c_s. This drag force acts in the direction -\hat V, and arises from the gravitational attraction between the perturber and its wake in the ambient medium. For supersonic motion (M=V/c_s>1), the enhanced-density wake is confined to the Mach cone trailing the perturber; for subsonic motion (M<1), the wake is confined to a sphere of radius c_s t centered a distance V t behind the perturber. Inside the wake, surfaces of constant density are hyperboloids or oblate spheroids for supersonic or subsonic perturbers, respectively, with the density maximal nearest the perturber. The dynamical drag force has the form F_df= - I 4\pi (G M_p)^2\rho_0/V^2. We evaluate I analytically; its limits are I\to M^3/3 for M<<1 and I\to ln(Vt/r_min) for M>>1. We compare our results to the Chandrasekhar formula for dynamical friction in a collisionless medium, noting that the gaseous drag is generally more efficient when M>1 but less efficient when M<1. To allow simple estimates of orbit evolution in a gaseous protogalaxy or proto-star cluster, we use our formulae to evaluate the decay times of a (supersonic) perturber on a near-circular orbit in an isothermal \rho\propto r^{-2} halo, and of a (subsonic) perturber on a near-circular orbit in a constant-density core. We also mention the relevance of our calculations to protoplanet migration in a circumstellar nebula.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,314 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.02 · states limitation 0.27 · human participants 0.01 · controlled comparison 0.24 · randomized assignment 0.01 · states quantitative result 0.89 · compares with prior results 0.98 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Using time-dependent linear perturbation theory, we evaluate the dynamical friction force on a massive perturber M_p traveling at velocity V through a uniform gaseous medium of density rho_0 and sound speed c_s. This drag force acts in the direction -\hat V, and arises from the gravitational attraction between the perturber and its wake in the ambient medium. For supersonic motion (M=V/c_s>1), the enhanced-density wake is confined to the Mach cone trailing the perturber; for subsonic motion (M<1), the wake is confined to a sphere of radius c_s t centered a distance V t behind the perturber. Inside the wake, surfaces of constant density are hyperboloids or oblate spheroids for supersonic or subsonic perturbers, respectively, with the density maximal nearest the perturber. The dynamical drag force has the form F_df= - I 4\pi (G M_p)^2\rho_0/V^2. We evaluate I analytically; its limits are I\to M^3/3 for M<<1 and I\to ln(Vt/r_min) for M>>1. We compare our results to the Chandrasekhar formula for dynamical friction in a collisionless medium, noting that the gaseous drag is generally more efficient when M>1 but less efficient when M<1. To allow simple estimates of orbit evolution in a gaseous protogalaxy or proto-star cluster, we use our formulae to evaluate the decay times of a (supersonic) perturber on a near-circular orbit in an isothermal \rho\propto r^{-2} halo, and of a (subsonic) perturber on a near-circular orbit in a constant-density core. We also mention the relevance of our calculations to protoplanet migration in a circumstellar nebula.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Using time-dependent linear perturbation theory, we evaluate the dynamical friction force on a massive perturber M_p traveling at velocity V through a uniform gaseous medium of density rho_0 and sound speed c_s. method 0.63 0.95 0.05 asserted 0.99
s1 This drag force acts in the direction -\hat V, and arises from the gravitational attraction between the perturber and its wake in the ambient medium. interpretation 0.41 0.52 0.06 asserted 1.00
s2 For supersonic motion (M=V/c_s>1), the enhanced-density wake is confined to the Mach cone trailing the perturber; for subsonic motion (M<1), the wake is confined to a sphere of radius c_s t centered a distance V t behind the perturber. result 1.00 0.78 0.05 asserted 1.00
s3 Inside the wake, surfaces of constant density are hyperboloids or oblate spheroids for supersonic or subsonic perturbers, respectively, with the density maximal nearest the perturber. result 1.00 0.83 0.04 asserted 1.00
s4 The dynamical drag force has the form F_df= - I 4\pi (G M_p)^2\rho_0/V^2. result 0.99 0.91 0.26 asserted 1.00
s5 We evaluate I analytically; its limits are I\to M^3/3 for M<<1 and I\to ln(Vt/r_min) for M>>1. result 0.98 0.96 0.12 asserted 1.00
s6 We compare our results to the Chandrasekhar formula for dynamical friction in a collisionless medium, noting that the gaseous drag is generally more efficient when M>1 but less efficient when M<1. interpretation 0.61 0.80 0.46 asserted 0.90
s7 To allow simple estimates of orbit evolution in a gaseous protogalaxy or proto-star cluster, we use our formulae to evaluate the decay times of a (supersonic) perturber on a near-circular orbit in an isothermal \rho\propto r^{-2} halo, and of a (subsonic) perturber on a near-circular orbit in a constant-density core. contribution 0.58 0.95 0.07 asserted 0.99
s8 We also mention the relevance of our calculations to protoplanet migration in a circumstellar nebula. interpretation 0.87 0.50 0.03 unclear 0.39

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — Density perturbation profiles for subsonic perturb
Figure 2. — Density perturbation profiles for supersonic pertu
Figure 3. — Solid curves: Dynamical friction force in a gaseous medium as a function of Mach number M = V/cs. Curves correspond to ln(cst/rmin) = 4, 6, 8, ..., 16. Dashed curves: Corresponding dynamical friction force in a collisionless medium with particle velocity dispersion σ = cs and rmax ≡ V t = Mcst
Figure 4. — (a) Decay in time of radial distance r of massive perturber Mp on near-circular orbit about the center of uniform-density core of mass Mcore, radius r0, and sound speed cs = (GMcore/(3r0))1/2 ; (b) Radial-to-azimuthal velocity ratio for same situation as in (a)
Figure 5. — For same situation as in Fig. 4, perturber’s radial d

Figures reproduced from arXiv:astro-ph/9810324 by the authors.

References
How later papers used it

223 later papers cite it

Evidence: substantial. 222 of the citing papers share no author with it, and 82 of those carry a confident label from Pith's citation classifier.

Pith read 31 of the 421 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 6 of 55 papers with this label
Compares with it 1 of 22 papers with this label
How the literature uses it

82 arXiv papers cite this in 171 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 101 sentences without a confident label

Contests it 1 paper
  • Since the original $I(ℳ)$ of [ref] is singular at $ℳ = 1$, we adopt the regularized form of [ref] , [equation] where the Coulomb logarithm is set to $lnΛ= 1.4 ∼ln(R_⋆/R_HL)$, the ratio of the largest to the smallest radial scales contributing to the drag.

    Formation of black hole stars via star--black hole collisions · Model setups and initial condition · contest 0.77

Builds on it 1 paper
Uses its method or tool 13 papers
Compares with it 8 papers
Cites it as support 3 papers
Cites it as background 21 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.