Acoustic instability at shock-wave precursors

Antonio Capanema; Emanuele Sobacchi; Pasquale Blasi

arxiv: 2604.12514 · v1 · pith:42P2ZMO7new · submitted 2026-04-14 · 🌌 astro-ph.HE · hep-ph· physics.plasm-ph

Acoustic instability at shock-wave precursors

Antonio Capanema , Pasquale Blasi , Emanuele Sobacchi This is my paper

Pith reviewed 2026-05-10 15:27 UTC · model grok-4.3

classification 🌌 astro-ph.HE hep-phphysics.plasm-ph

keywords acoustic instabilitycosmic-ray modified shocksshock precursorsmagnetic field amplificationturbulence generationparticle accelerationsupernova remnantsMHD simulations

0 comments

The pith

With realistic Mach numbers and cosmic-ray efficiencies, the acoustic instability grows small density perturbations into large nonlinear structures as plasma crosses a shock precursor.

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

The paper investigates how small density variations in the upstream region of a cosmic-ray-modified shock interact with the cosmic-ray pressure gradient. It shows that under parameters closer to those expected in supernova remnants than in earlier models, this interaction drives the acoustic instability to nonlinear amplitudes before the fluid reaches the shock. The resulting vortical motions generate turbulence capable of amplifying magnetic fields. This process is studied through two-dimensional magnetohydrodynamic simulations that impose a fixed cosmic-ray pressure profile. The work also examines the spectrum of resulting magnetic fluctuations and notes possible coupling with streaming instabilities.

Core claim

By performing two-dimensional magnetohydrodynamic simulations with an imposed cosmic-ray pressure gradient, we demonstrate that the acoustic instability grows small upstream density perturbations into large-amplitude nonlinear structures as the plasma flows through the precursor of a cosmic-ray-modified shock, for Mach numbers and cosmic-ray acceleration efficiencies closer to realistic values than in earlier work.

What carries the argument

The acoustic instability driven by the coupling between density perturbations and the cosmic-ray pressure gradient in the shock precursor, which induces vorticity and turbulence.

If this is right

The instability produces a spectrum of turbulent magnetic fluctuations that can scatter energetic particles.
Vorticity generated by the nonlinear structures contributes to magnetic-field amplification upstream of the shock.
Constructive interference may occur between the acoustic instability and non-resonant streaming instabilities.
Current numerical resolutions limit access to the smallest scales where turbulence becomes fully nonlinear.
Future work can explore three-dimensional evolution and self-consistent cosmic-ray feedback.

Where Pith is reading between the lines

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

If the mechanism operates, it supplies an additional channel for magnetic-field amplification that could raise the maximum energy reachable by cosmic rays in supernova remnants.
The resulting turbulence spectrum might imprint observable features on radio synchrotron emission from young remnants.
Extending the simulations to include particle-in-cell or hybrid kinetic treatments could test whether the amplified fields alter the cosmic-ray acceleration efficiency itself.

Load-bearing premise

An externally prescribed cosmic-ray pressure gradient accurately represents the self-consistent precursor of a cosmic-ray-modified shock and two-dimensional simulations capture the essential nonlinear evolution.

What would settle it

High-resolution observations or simulations that show the precursor density perturbations remain linear in amplitude throughout the crossing time for realistic Mach numbers and acceleration efficiencies would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.12514 by Antonio Capanema, Emanuele Sobacchi, Pasquale Blasi.

**Figure 1.** Figure 1: Real and imaginary parts of ω. The CR pressure gradient is constant, |∇PCR| = ξCRρ0u 2 0 /L, and the perturbation wave vector k is aligned with ∇PCR. We assume ξCR = 0.1, Ms = 100, and Mms = 100/ √ 5. The dashed and dotted lines show the asymptotic behavior of ω. find −π/2 < ϕ < 0 if we take waves in the opposite direction (k · ∇PCR = −k|∇PCR| or equivalently ωR, ωI < 0), making pressure precede velocity … view at source ↗

**Figure 2.** Figure 2: Quasi-stationary x-profiles of the dimensionless density, pressure, velocity, and magnetic field perturbations from kx modes. We show the simulation results as solid lines and the expected growths for AI derived in Sect. 2 as dashed lines. In the bottom left panel, the initial magnetic field pointed along y, and the lighter gray line corresponds to the expected growth in the absence of magnetic fields. The… view at source ↗

**Figure 3.** Figure 3: Stationary profiles of the dimensionless density (left) and velocity (right) perturbations from ky modes. On the right, red/blue/black lines show the simulation results for δux/cs at a fixed y coordinate corresponding to the peak/trough/node of δρ/ρ0 oscillations. Dotted and dashed lines correspond to the analytical predictions of Eq. (20). Top: No magnetic fields. Bottom: Magnetic field perpendicular to t… view at source ↗

**Figure 4.** Figure 4: Linear-regime simulation results for initial density fluctuations following Eq. (23) with a flat power spectrum in a 2048 × 256 grid. Left: Perturbation profiles for (from top to bottom) density, velocity x- and y-components, and magnetic field x- and y-components. Note that δux/cs has been multiplied by 1/3 and δuy/cs by 2, so all quantities can be visualized using the same color bar scale. Right: dimensi… view at source ↗

**Figure 5.** Figure 5: 2D snapshots of quasi-stationary density, velocity, magnetic field, and magnetic energy density profiles, for different parameter choices, as indicated above each column, in a 4096 × 512 simulation grid. All quantities have been normalized so as to become dimensionless. For better contrast, colorbars use linear scaling for ρ/ρ0 and ux/cs , logarithmic scaling for (B/B0) 2 , and symmetric logarithmic scalin… view at source ↗

**Figure 7.** Figure 7: Left: Turbulent kinetic (thick lines) and magnetic (thin lines) power spectra in units of ρ0c 2 s for different grid resolutions (see right panel’s legend). Right: MFA as a function of x for different grid resolutions. To quantitatively analyze our results, we calculate both kinetic and magnetic omnidirectional turbulent energy density spectra, Eturb,kin(k) and Eturb,mag(k) respectively, as defined in Ap… view at source ↗

**Figure 8.** Figure 8: Left: Turbulent kinetic (thick lines) and magnetic (thin lines) power spectra in units of ρ0c 2 s for 2D (1024 × 128 grid; blue solid) and 3D (1024 × 128 × 128 grid; yellow dashed) simulations. Center: MFA as a function of x for 2D and 3D simulations. In 3D, we perform a yz-averaging rather than just a y-averaging. Right: Inverse eddy turnover timescales for the 2D (blue solid) and 3D (yellow dashed) runs,… view at source ↗

read the original abstract

Magnetic field amplification is an integral part of the process of particle acceleration at non-relativistic shocks. It is necessary to reach the maximum energies required by observations, especially in supernova remnants, thought to be sources of the bulk of Galactic cosmic rays. Such amplification can be caused by the acoustic instability that develops when small density perturbations interact with the cosmic-ray pressure gradient in the upstream of a cosmic-ray-modified shock. The vorticity induced by the nonlinear development of the instability may lead to turbulence, which amplifies the pre-existing magnetic fields. To study this phenomenon, we use the PLUTO code to carry out 2D (and some 3D) magnetohydrodynamical simulations of the evolution of small density perturbations in the presence of an assigned cosmic-ray pressure gradient. Adopting more realistic values of Mach number and cosmic-ray acceleration efficiency than previously assumed in the literature, we show that the acoustic instability can transform small density perturbations into large nonlinear structures while the fluid crosses the precursor region of a cosmic-ray-modified shock. We study the power spectrum of turbulent magnetic fluctuations that may be important to scatter particles. We comment on the possible constructive interference between acoustic and non-resonant streaming instabilities. We discuss limitations of previous and current numerical investigations in accessing spatial scales where turbulence is expected to turn nonlinear, and outline perspectives for future investigations.

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.

Desk Editor's Note private letter to a colleague

This paper runs the acoustic instability with more realistic Mach numbers and cosmic-ray efficiencies than earlier work, showing growth to nonlinear structures and a turbulent magnetic spectrum, but the fixed pressure gradient is a clear limitation once amplitudes get large.

read the letter

The new element is the choice of parameters closer to supernova remnant shocks. In 2D PLUTO MHD runs with an imposed cosmic-ray pressure gradient, small density perturbations grow into large nonlinear features over the precursor crossing time, and the authors extract the resulting magnetic fluctuation spectrum. They also flag possible overlap with the non-resonant streaming instability and are open about the scale and dimensionality limits of the current runs. That quantitative step with realistic numbers is the main advance over prior literature that used milder conditions where growth stayed weaker. The honest discussion of what the numerics cannot yet reach is also useful. The soft spot is the assigned gradient. At the Mach numbers and efficiencies they adopt, the precursor is thinner, so the time window for growth is shorter while the chance that growing vorticity and density fluctuations feed back on the cosmic rays is higher. The simulations omit that back-reaction by construction, and the stress-test concern lands here. The paper notes its own limitations on scales, which helps, but the self-consistency issue remains for the nonlinear regime. This is for modelers of cosmic-ray acceleration and magnetic amplification at shocks who need concrete spectra or parameter checks. A reader working on supernova remnant maximum energies or turbulence would find the results worth looking at, even if they treat the gradient as an approximation. It deserves peer review. The numerical setup is described clearly enough for referees to assess the growth claims and suggest how to test the back-reaction in follow-up work.

Referee Report

3 major / 2 minor

Summary. The manuscript presents 2D (with some 3D) MHD simulations using the PLUTO code of small density perturbations evolving in the presence of an imposed cosmic-ray pressure gradient representing the precursor of a cosmic-ray-modified shock. Adopting more realistic Mach numbers and cosmic-ray acceleration efficiencies than prior literature, the authors claim that the acoustic instability grows these perturbations into large nonlinear structures during precursor transit, inducing vorticity that can drive turbulence and amplify magnetic fields. They analyze the power spectrum of the resulting magnetic fluctuations and discuss possible constructive interference with the non-resonant streaming instability, while noting limitations in resolving fully nonlinear turbulence scales.

Significance. If the numerical results hold under the stated approximations, the work is significant for models of magnetic-field amplification and particle acceleration at non-relativistic shocks, such as those in supernova remnants. It advances prior studies by demonstrating the instability's efficacy with observationally motivated parameters, showing that small seeds can reach nonlinear amplitudes within the available transit time and providing quantitative spectra of turbulent fluctuations. A clear strength is the explicit adoption of realistic Mach numbers and efficiencies together with the exploration of 2D-to-3D differences. The imposed-gradient approach, however, limits self-consistency, so the significance hinges on whether back-reaction effects remain negligible.

major comments (3)

[§2] §2 (Numerical setup): The cosmic-ray pressure gradient is externally assigned rather than obtained from a self-consistent CR transport equation. This assumption is load-bearing for the central claim, because any back-reaction from the growing density and vorticity fluctuations on CR streaming or diffusion is omitted by construction. For the realistic (higher) Mach numbers and acceleration efficiencies adopted, the precursor is thinner, shrinking the growth window while increasing the potential impact of omitted back-reaction; explicit sensitivity tests (e.g., gradient amplitude variations or a minimal coupled CR module) are needed to support the result.
[Results sections] Results on nonlinear evolution (likely §3–4): The assertion that small perturbations become 'large nonlinear structures' while the fluid crosses the precursor rests on the simulations reaching nonlinear amplitudes within the transit time. No quantitative diagnostics—such as measured growth factors, time to nonlinearity relative to precursor crossing time, or direct comparison with linear acoustic-instability theory—are reported in the provided description, making it difficult to verify that the claimed transformation occurs under the stated conditions.
[Methods and results] Dimensionality and resolution: Primary results are 2D, with only limited 3D runs mentioned. In 2D MHD, inverse energy cascades can artificially enhance large-scale structures; the manuscript should demonstrate that the reported nonlinear structures and magnetic power spectra remain robust when the same parameters are evolved in 3D at comparable or higher resolution.

minor comments (2)

[Abstract] Abstract: The strong phrasing 'we show that the acoustic instability can transform...' should be tempered to reflect that the evidence is numerical and subject to the imposed-gradient approximation.
[Throughout] Figure captions and text: Ensure all initial perturbation amplitudes, grid resolutions, and exact functional form of the imposed CR gradient (including how it is held fixed) are stated explicitly so that the setup can be reproduced.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report, which highlights both the potential significance of the work and areas where the presentation and supporting evidence can be strengthened. We address each major comment below, indicating the revisions we intend to incorporate in the revised manuscript.

read point-by-point responses

Referee: [§2] §2 (Numerical setup): The cosmic-ray pressure gradient is externally assigned rather than obtained from a self-consistent CR transport equation. This assumption is load-bearing for the central claim, because any back-reaction from the growing density and vorticity fluctuations on CR streaming or diffusion is omitted by construction. For the realistic (higher) Mach numbers and acceleration efficiencies adopted, the precursor is thinner, shrinking the growth window while increasing the potential impact of omitted back-reaction; explicit sensitivity tests (e.g., gradient amplitude variations or a minimal coupled CR module) are needed to support the result.

Authors: We agree that the externally imposed cosmic-ray pressure gradient constitutes an approximation that omits possible back-reaction of the developing fluid perturbations on cosmic-ray transport. This choice follows the approach used in earlier studies of the acoustic instability and permits isolation of the instability mechanism under controlled conditions. Nevertheless, the referee correctly notes that the higher Mach numbers and efficiencies we adopt result in a thinner precursor, which could amplify the importance of omitted effects. In the revised manuscript we will add a dedicated sensitivity study in which the imposed gradient amplitude is varied over a range consistent with the adopted acceleration efficiency; we will also expand the discussion in §2 to quantify the expected regime of validity of the approximation and to outline why a fully coupled CR module lies beyond the scope of the present work. revision: partial
Referee: [Results sections] Results on nonlinear evolution (likely §3–4): The assertion that small perturbations become 'large nonlinear structures' while the fluid crosses the precursor rests on the simulations reaching nonlinear amplitudes within the transit time. No quantitative diagnostics—such as measured growth factors, time to nonlinearity relative to precursor crossing time, or direct comparison with linear acoustic-instability theory—are reported in the provided description, making it difficult to verify that the claimed transformation occurs under the stated conditions.

Authors: The referee is correct that the manuscript would be strengthened by explicit quantitative diagnostics supporting the claim that nonlinearity is reached within the precursor transit time. In the revised version we will report (i) the measured exponential growth factor of the density contrast as a function of time, (ii) the simulation time at which the density perturbation first reaches order-unity amplitude, normalized to the precursor crossing time, and (iii) a direct comparison of the observed growth rate with the analytic prediction from linear acoustic-instability theory for the adopted Mach number and cosmic-ray pressure gradient. These diagnostics will be added to the results sections and to the associated figures. revision: yes
Referee: [Methods and results] Dimensionality and resolution: Primary results are 2D, with only limited 3D runs mentioned. In 2D MHD, inverse energy cascades can artificially enhance large-scale structures; the manuscript should demonstrate that the reported nonlinear structures and magnetic power spectra remain robust when the same parameters are evolved in 3D at comparable or higher resolution.

Authors: We acknowledge that two-dimensional MHD simulations are susceptible to an inverse energy cascade that can artificially sustain large-scale structures. Although the manuscript already contains a limited set of three-dimensional runs, these are insufficient to fully address the referee’s concern. In the revision we will present additional three-dimensional simulations performed at resolutions comparable to the two-dimensional cases, together with side-by-side comparisons of the density and magnetic-field structures and of the magnetic power spectra. We will also add a brief discussion of the expected impact of dimensionality on the inverse cascade and on the robustness of the reported spectral slopes. revision: partial

Circularity Check

0 steps flagged

No significant circularity: results follow from direct numerical integration of standard equations with imposed gradient

full rationale

The paper conducts 2D/3D MHD simulations in PLUTO with an externally assigned cosmic-ray pressure gradient to evolve small density perturbations. The central result—that the acoustic instability grows into nonlinear structures for realistic Mach numbers and acceleration efficiencies—is obtained by integrating the governing equations forward in time rather than by algebraic reduction, fitting, or self-referential definition. No parameter is fitted to a subset of data and then relabeled a prediction; no uniqueness theorem or ansatz is imported via self-citation; and the imposed gradient is explicitly stated as an approximation whose limitations are discussed. The derivation chain is therefore self-contained against external benchmarks (standard MHD + prescribed source term) and receives a score of 0.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard MHD evolution under an externally imposed cosmic-ray pressure gradient; no new physical entities are introduced and the only adjustable inputs are the adopted Mach number and acceleration efficiency taken from the literature.

free parameters (2)

Mach number
Chosen as more realistic than prior studies; value not specified in abstract but treated as an input parameter.
cosmic-ray acceleration efficiency
Chosen as more realistic than prior studies; value not specified in abstract but treated as an input parameter.

axioms (2)

standard math Standard ideal MHD equations govern the plasma flow
Implemented via the PLUTO code for the simulations.
domain assumption An externally assigned cosmic-ray pressure gradient adequately represents the precursor of a cosmic-ray-modified shock
Used to drive the acoustic instability in the upstream region.

pith-pipeline@v0.9.0 · 5535 in / 1604 out tokens · 74831 ms · 2026-05-10T15:27:33.448103+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

[1]

For density fluctuations, the field isf=δρ/ρ 0, and we chooseC=1 such thatX n Eδρ/ρ0(¯kn)∆kn = X ab ∆2k P(kab)=ξ(0) = 1 A X i j ∆x∆y " δρ ρ0 (xi j) #2 = δρrms ρ0 !2 ,(B.12) whereδρ rms /ρ0 is the RMSδρ/ρ 0 over all scales in the box

work page
[2]

In other words, δBrms(¯kn)= q 8πEturb,mag(¯kn)∆kn (B.14) is the RMS value of magnetic field fluctuations with charac- teristic scaleℓ=2π/ ¯kn

For magnetic fields fluctuations,f=δB, we chooseC= 1/8πwhereby the turbulent magnetic energy density power spectrumE δB(k)≡E turb,mag(k) satisfies X n Eturb,mag(¯kn)∆kn = 1 A X i j ∆x∆y |δB(xi j)|2 8π = δB2 rms 8π (B.13) is the average magnetic energy density in perturbations at all scales, whileE turb,mag(¯kn)∆kn =δB 2 rms(¯kn)/8πis that coming from pert...

work page
[3]

For velocity fields,e.g.δu, we chooseC=ρ 0/2 such that X n Eδu(¯kn)∆kn = 1 A X i j ∆x∆y ρ0|δu(xi j)|2 2 = ρ0δu2 rms 2 . (B.15) Analogously to the magnetic field case, the RMS amplitude of velocity perturbations at a characteristic scaleℓ=2π/ ¯kn can be found through δurms(¯kn)= s 2Eδu(¯kn)∆kn ρ0 .(B.16)

work page
[4]

For studies of compressible turbulence, it is common to consider the density-weighted velocityw= √ρuto de- fine kinetic energy spectra (Kida & Orszag 1990; Grete et al. 2017). Letδwbe its fluctuations around a mean background profile. We can then chooseC=1/2 such that its corre- sponding power spectrumE δw(k)≡E turb,kin(k) is the turbu- lent kinetic energ...

work page 1990

[1] [1]

For density fluctuations, the field isf=δρ/ρ 0, and we chooseC=1 such thatX n Eδρ/ρ0(¯kn)∆kn = X ab ∆2k P(kab)=ξ(0) = 1 A X i j ∆x∆y " δρ ρ0 (xi j) #2 = δρrms ρ0 !2 ,(B.12) whereδρ rms /ρ0 is the RMSδρ/ρ 0 over all scales in the box

work page

[2] [2]

In other words, δBrms(¯kn)= q 8πEturb,mag(¯kn)∆kn (B.14) is the RMS value of magnetic field fluctuations with charac- teristic scaleℓ=2π/ ¯kn

For magnetic fields fluctuations,f=δB, we chooseC= 1/8πwhereby the turbulent magnetic energy density power spectrumE δB(k)≡E turb,mag(k) satisfies X n Eturb,mag(¯kn)∆kn = 1 A X i j ∆x∆y |δB(xi j)|2 8π = δB2 rms 8π (B.13) is the average magnetic energy density in perturbations at all scales, whileE turb,mag(¯kn)∆kn =δB 2 rms(¯kn)/8πis that coming from pert...

work page

[3] [3]

For velocity fields,e.g.δu, we chooseC=ρ 0/2 such that X n Eδu(¯kn)∆kn = 1 A X i j ∆x∆y ρ0|δu(xi j)|2 2 = ρ0δu2 rms 2 . (B.15) Analogously to the magnetic field case, the RMS amplitude of velocity perturbations at a characteristic scaleℓ=2π/ ¯kn can be found through δurms(¯kn)= s 2Eδu(¯kn)∆kn ρ0 .(B.16)

work page

[4] [4]

For studies of compressible turbulence, it is common to consider the density-weighted velocityw= √ρuto de- fine kinetic energy spectra (Kida & Orszag 1990; Grete et al. 2017). Letδwbe its fluctuations around a mean background profile. We can then chooseC=1/2 such that its corre- sponding power spectrumE δw(k)≡E turb,kin(k) is the turbu- lent kinetic energ...

work page 1990