pith. sign in

arxiv: 2604.22665 · v1 · submitted 2026-04-24 · 🌌 astro-ph.GA · astro-ph.CO· astro-ph.EP· astro-ph.IM· astro-ph.SR

The Sound of the Universe: A Resonant Gravitational Instability Driven by Baryon-Dark Matter Relative Drift

Mohamad Shalaby , Avery Broderick This is my paper

Pith reviewed 2026-05-08 10:53 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.COastro-ph.EPastro-ph.IMastro-ph.SR
keywords resonant gravitational instabilitybaryon-dark matter relative driftJeans scalesound wavescollisionless dragastrophysical instabilitiesdark matter probes
0
0 comments X

The pith

Relative drift between baryons and dark matter after decoupling triggers a resonant gravitational instability that amplifies baryonic sound waves faster than standard cold dark matter growth on sub-Jeans scales.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that baryons and dark matter develop a relative velocity after early-universe decoupling. When this drift is slower than the sound speed in the baryons, their stable oscillatory modes resonate with the Doppler-shifted dark matter perturbations, driving exponentially growing density waves. The growth exceeds the usual rate from dark matter gravity alone and peaks below the baryon Jeans scale, with a window of full stability in between in baryon-rich settings. Timescales range from years to tens of millions of years across planets, disks, stars, clouds, galaxies, and clusters, often shorter than the systems' ages. The coupling also produces a collisionless drag that transfers momentum between the components and, in an expanding universe, makes growth depend on the orientation of the mode relative to the drift direction.

Core claim

Dark matter and baryons acquire a relative velocity after decoupling in the early Universe. Their relative drift triggers a resonant gravitational instability that drives sound waves in baryons. When the projected DM drift is subsonic, the stable oscillatory branch of baryons resonates with the Doppler-shifted DM mode, producing exponentially growing perturbations whose growth rates exceed the intrinsic CDM growth rate. The instability peaks below the baryon Jeans scale and, in baryon-dominated environments, opens a window of complete stability between the Jeans scale and the resonance.

What carries the argument

Resonant coupling between the stable baryon oscillatory sound-wave branch and the Doppler-shifted cold dark matter density mode driven by their relative drift velocity.

If this is right

  • The instability enhances baryon density perturbations for modes oriented appropriately to the drift in an expanding universe while suppressing those aligned with the dark matter stream.
  • Momentum transfer creates a non-viscous collisionless drag between the species.
  • Growth timescales from years to tens of millions of years across planets, protoplanetary disks, stars, molecular clouds, galaxies, and clusters are typically much shorter than system ages.
  • Supersonic relative drifts suppress the growth, and the mechanism may explain the persistence of spiral arms and heating of the intracluster medium.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Cosmological simulations that include relative velocities would show increased small-scale baryon power for certain orientations relative to the dark matter flow.
  • The resonance could be tested by searching for velocity-aligned overdensities in high-resolution maps of gas and dark matter in nearby galaxies or clusters.
  • In protoplanetary disks the short growth times suggest the instability may drive rapid density variations that affect early stages of planet formation.

Load-bearing premise

Perturbations remain small enough for the linear calculation to hold and the relative speed between baryons and dark matter stays below the sound speed long enough for the waves to grow in the environments examined.

What would settle it

A simulation or observation showing no faster-than-CDM growth of baryon perturbations when the relative drift is subsonic, or no directional dependence in perturbation growth aligned with or against the dark matter flow in an expanding universe.

Figures

Figures reproduced from arXiv: 2604.22665 by Avery Broderick, Mohamad Shalaby.

Figure 1
Figure 1. Figure 1: Graphical representation for the physical origin of the driven resonant instability. This figure shows the solutions of Equa￾tion (1) in various limiting cases. The oscillation frequency for the baryon-only case (ρDM = 0) is shown in blue, while the cold DM￾only case is shown in black. For the non-drifting cold DM case (rel￾ative to the baryon rest frame), the normal modes are non-oscillating (ω = 0), repr… view at source ↗
Figure 2
Figure 2. Figure 2: The impact of dark matter (DM) drift on the gravitational field perturbation growth rates for two mass density ratios: ΩDM = 10−2Ωb (baryon-dominated, left panel, with logarithmic y-axes) and ΩDM = Ωb (comparable densities, right panel, with linear y-axes). For each density ratio, growth rates are shown as a function of wavenumber for different relative drift velocities: subsonic drift vr = 0.9 (blue), sup… view at source ↗
Figure 3
Figure 3. Figure 3: The typical time scale for the resonantly driven sound waves as given in Equation (5), where we use vr = 0.5. We nor￾malize the both cold DM baryon mass densities with the Interstellar medium average density of 1 proton (with mass mp) in cm3 . The the growth time can be much shorter if the value of vr is chosen differently; see discussion at the end of section 4.1. White contours show the lines of equal gr… view at source ↗
Figure 4
Figure 4. Figure 4: Regime of unstable sound waves (light blue shaded regime) due to resonant instabilities in case of supersonic DM rela￾tive drift, i.e., vDM > cs along the blue vertical direction. In case of subsonic drift, the cone of stable sound waves closes, and sound waves are resonantly unstable in all direction. Modes are stable on scales below the Jeans scale of Baryons. When vr ≪ 1, the resonance appears on scales… view at source ↗
Figure 5
Figure 5. Figure 5: Cosmological evolution of baryon density perturbation δb with k = 104 Mpc−1 ). Colors indicate: no drift in dashed red, modes along DM streaming direction: solid red (cos θ = 1, super￾sonic drift, density growth suppression), early subsonic drift at re￾combination (black, cos θ = 1/5, early density amplification), late subsonic drift: z ≳ 200, blue (cos θ = 0.47, extra amplification at late times). subsoni… view at source ↗
Figure 6
Figure 6. Figure 6: The top panel shows the oscillatory behavior, while the bottom panel shows the growth, of gravitational perturbations in the presence of cold DM (cDM = 0) with various densities. we empha￾size here that the solutions in top and bottom panels are different branches, that is for real solution shown in the top penal, ℑ{ω} = 0, and for the growing solutions shown in the bottom panel ℜ{ω} = 0. As we show in App… view at source ↗
Figure 7
Figure 7. Figure 7: Same as in view at source ↗
Figure 9
Figure 9. Figure 9: Validating the approximate solution for the resonantly driven growth rate given in Equation (C38): The dependence of the growth rate on R at various values of ur = vr/ √ 1 − v 2 r . This, as expected, does not agree very well with the imagi￾nary part of the full solution for values of R close to 1. How￾ever, inspired by this solution, we can construct a numerical formula that agrees much better with the nu… view at source ↗
Figure 8
Figure 8. Figure 8: Solutions for the dispersion relation (Equation (C30)). Vertical gray lines in the last 2 columns show the expected reso￾nance scale given by Equation (C31), and the right-most column is a zoom of the results in the middle column. This show that the streaming-driven resonance has growth rate which is maximized at the resonance wavelength given by kcs/Ωb = 1/ √ 1 − v 2 r . The solution of Equation (C32) is … view at source ↗
Figure 10
Figure 10. Figure 10: Validating the approximate solution for the resonant growth rate given in Equation (C38): The dependence of the growth rate on ur at various values of R. C.2.1. Resonance Location Assuming R ≪ 1, the location of the fastest growth rates can be found as follows (x − yvr) = ±ycr → x = y(vr ± cr) = p y 2 − 1 → y 2 = 1 1 − (vr ± cr) 2 (C41) From solving the dispersion relation in this case (restricted to the … view at source ↗
Figure 11
Figure 11. Figure 11: Top panel show resonant instability growth rates (vr − cr = 0.9 for all cases) at various sound speed ratios (cr = cDM/cs, and DM to baryon density ratios, R = ΩDM/Ωb. Bottom panel shows the case of supersonic relative drift (vr − cr = 2 for all cases). For comparison, we show for all cases, the growth rates for the cases of cold non-drifting DM are shown in green in all panels. (Equation 11 of Tseliakhov… view at source ↗
Figure 12
Figure 12. Figure 12: The cosmological evolution of the DM–baryon relative drift, normalized to the evolving baryon sound speed, is shown in red. The black curve indicates the minimum angle θmin between the wave mode and the DM streaming direction for which the projected speed equals the sound speed, i.e., vr = vDM cos θ/cs = 1. All wave modes with angles θ ∈ [θmin, π/2) can be resonantly driven. That is while the range for re… view at source ↗
read the original abstract

Dark matter and baryons acquire a relative velocity after decoupling in the early Universe. Baryons are gravitationally unstable only above their Jeans scale, while cold dark matter (CDM) is unstable on all scales. We show for the first time that their relative drift triggers a resonant gravitational instability that drives sound waves in baryons. When the projected DM drift is subsonic, the stable oscillatory branch of baryons resonates with the Doppler-shifted DM mode, producing exponentially growing perturbations whose growth rates exceed the intrinsic CDM growth rate. The instability peaks below the baryon Jeans scale and, in baryon-dominated environments, opens a window of complete stability between the Jeans scale and the resonance. Supersonic drift suppresses growth, as previously noted. The resonant coupling also transfers momentum between the species, creating a non-viscous, collisionless drag. We derive an accurate analytical approximation for the growth rate at resonance and show that the associated timescales range from years to tens of millions of years across diverse environments -- planets, protoplanetary disks, stars, molecular clouds, galaxies, and galaxy clusters -- typically much shorter than their ages. In an expanding FLRW universe, the instability enhances baryon density perturbations at different redshifts for appropriately oriented modes while suppressing the growth of those aligned with the DM stream. The universe thus sings across all scales, and this resonant mechanism provides the means to listen: it offers a novel probe of dark matter through its seismic imprint on astrophysical objects and may explain long-standing puzzles such as the persistence of spiral arms and the heating of the intracluster medium in galaxy clusters.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript claims that the relative drift velocity between baryons and cold dark matter acquired after decoupling triggers a resonant gravitational instability. When the projected DM drift is subsonic, the stable baryon acoustic branch resonates with the Doppler-shifted DM mode, producing exponentially growing perturbations whose growth rates exceed the intrinsic CDM rate. The instability peaks below the baryon Jeans scale; an accurate analytical approximation for the resonant growth rate is derived, yielding timescales from years to tens of Myr across planets, disks, stars, clouds, galaxies and clusters. Supersonic drifts suppress growth. Momentum transfer induces collisionless drag. In an expanding FLRW universe the mechanism enhances baryon perturbations for appropriately oriented modes while suppressing those aligned with the DM stream.

Significance. If the resonance derivation and linear-regime persistence are confirmed, the result supplies a new, scale-independent probe of dark-matter properties through its seismic imprint on baryonic structures and offers a potential explanation for the persistence of spiral arms and intracluster-medium heating. The short growth timescales relative to system ages and the provision of an analytical growth-rate formula are notable strengths that would make the mechanism observationally testable across diverse environments.

major comments (2)
  1. Abstract: The central claim of sustained exponential growth (exceeding CDM rates) over timescales from years to ~10 Myr rests on the assumption that the linear perturbation analysis remains valid and that the mean relative drift velocity stays subsonic. No estimate is supplied for the time required to reach δ ≈ 1 or for the evolution of the bulk relative velocity under the collisionless drag that the resonance itself generates; if either condition fails before significant amplification, the resonance window closes and the mechanism cannot operate as stated.
  2. Abstract: The manuscript states that an 'accurate analytical approximation for the growth rate at resonance' is derived, yet supplies neither the dispersion relation nor the steps leading to that approximation. Without these, it is impossible to verify whether the resonance condition is obtained from first principles or incorporates post-hoc assumptions that affect the claimed growth-rate excess over CDM.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and insightful review of our manuscript. We appreciate the opportunity to clarify the points raised regarding the validity of the linear regime and the presentation of our analytical results. Below, we provide point-by-point responses to the major comments and outline the revisions we will implement.

read point-by-point responses

  1. Referee: Abstract: The central claim of sustained exponential growth (exceeding CDM rates) over timescales from years to ~10 Myr rests on the assumption that the linear perturbation analysis remains valid and that the mean relative drift velocity stays subsonic. No estimate is supplied for the time required to reach δ ≈ 1 or for the evolution of the bulk relative velocity under the collisionless drag that the resonance itself generates; if either condition fails before significant amplification, the resonance window closes and the mechanism cannot operate as stated.

    Authors: We agree that explicit estimates for the time to reach δ ≈ 1 and the evolution of the mean relative drift velocity under collisionless drag are necessary to substantiate sustained growth. In the revised manuscript we will add a new subsection that integrates the resonant growth rate to estimate the saturation time and solves the coupled mean-velocity equations to track the deceleration of the bulk drift. These calculations show that, for subsonic drifts, δ reaches order unity within the quoted growth timescales while the mean velocity remains subsonic throughout the linear phase, thereby keeping the resonance window open. revision: yes

  2. Referee: Abstract: The manuscript states that an 'accurate analytical approximation for the growth rate at resonance' is derived, yet supplies neither the dispersion relation nor the steps leading to that approximation. Without these, it is impossible to verify whether the resonance condition is obtained from first principles or incorporates post-hoc assumptions that affect the claimed growth-rate excess over CDM.

    Authors: The two-fluid dispersion relation with relative drift is stated in the main text, and the resonant growth-rate approximation is obtained from it by imposing the resonance condition and performing a perturbative expansion. To improve transparency and allow direct verification, the revised manuscript will include the explicit quartic dispersion relation, the resonance matching condition, and the full sequence of algebraic steps used to isolate the imaginary frequency. This will confirm that the excess growth over standard CDM follows directly from the first-principles equations without additional assumptions. revision: yes

Circularity Check

0 steps flagged

No significant circularity: growth rate derived from standard two-fluid dispersion relation

full rationale

The paper applies linear perturbation theory to the coupled baryon-DM fluid equations with an imposed relative drift velocity. The resonant instability and its analytical growth-rate approximation follow directly from solving the resulting dispersion relation under the subsonic-drift condition; no equation defines the output growth rate in terms of itself, no fitted parameter is relabeled as a prediction, and no load-bearing premise rests on a self-citation chain. The claim that growth exceeds the CDM rate is a direct consequence of the resonance term in the dispersion relation rather than a tautology. External references (e.g., to supersonic suppression) are not required to establish the new subsonic resonant branch.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on linear two-fluid perturbation theory in an expanding background plus the existence of a post-decoupling relative velocity; no new particles or forces are introduced, but the resonance condition is derived rather than postulated.

axioms (2)
  • domain assumption Linear perturbation theory remains valid for the coupled baryon-DM system on the scales of interest
    Invoked implicitly when stating exponentially growing modes and resonance between oscillatory and drift branches
  • domain assumption A relative velocity between baryons and dark matter persists after decoupling and can be treated as approximately constant over the instability growth time
    Stated in the opening sentence and used to define subsonic vs supersonic regimes

pith-pipeline@v0.9.0 · 5612 in / 1564 out tokens · 25470 ms · 2026-05-08T10:53:31.020864+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    For the force terms, we linearize equation (A5)

    Then, we add the first order perturbation, i.e., we take ρb =ρ b,0 +ρ b,1 &⃗ vb =⃗ vb,0 +⃗ vb,1 pb =p b,0 +p b,1 &f DM =f DM,0 +f DM,1 ⇒⃗ a=⃗ a0 +⃗ a1 orϕ=ϕ 0 +ϕ 1 (A6) For the baryons, we linearize equations (A1) and (A3), and for the DM, we linearize equation (A4). For the force terms, we linearize equation (A5). We assume that only the linear terms dep...

    work page 2008

  2. [2]

    × 10-4 0.001 0.005 0.010 Figure 10.Validating the approximate solution for the resonant growth rate given in Equation (C38): The dependence of the growth rate onu r at various values ofR. C.2.1.Resonance Location AssumingR≪1, the location of the fastest growth rates can be found as follows (x−yv r) =±yc r →x=y(v r ±c r) = p y2 −1 →y 2 = 1 1−(v r ±c r)2 (C...

    work page 2010