Pith. Sign in

REVIEW

Causal self-consistency of the Blandford--McKee self-similar solution

High Energy Astrophysical Phenomena · Submitted 31 Aug 2026 · 296 references

Pith's verdict: A genuine and mostly sound resolution of the BM self-consistency gap, conditional on the simplified equation of state used to build the sonic barrier.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 9, 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 →

· 2 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.00116 (version 1)
Pith IDDELIM24N, Pith's permanent record of this paper and its review
CitationDoron Kushnir, Eli Waxman, Gilad Sadeh, Tamar Faran. "Causal self-consistency of the Blandford--McKee self-similar solution." arXiv:2609.00116 (2026). https://pith.science/paper/DELIM24N
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
The paper

The paper claims

The Blandford-McKee solution near the shock is causally protected from the deep interior, closing a gap in its derivation.

The reading

The Blandford-McKee (BM) solution describes the ultra-relativistic flow behind a strong shock, but only holds near the shock; deeper downstream its assumptions fail. This raises a consistency problem: if the interior could send acoustic signals to the shock, the BM near-shock profiles would depend on unknown interior conditions. The paper shows that the flow itself prevents such communication. For shallow external density profiles (k < k_g ≈ 2.062), the BM characteristic structure alone blocks signals. For steeper profiles (k_g < k < 4), the flow develops a cooling layer that ends at a sonic point, beyond which no forward acoustic signals can reach the shock. Thus the BM near-shock solution is causally self-consistent.

What carries the argument

The analysis uses two main tools: (i) characteristic analysis of the hot BM solution in the similarity coordinate χ, giving the drift of forward characteristics relative to the breakdown coordinate; (ii) construction of a cooling self-similar solution following the composite method of Pan & Sari (2009), using the equation of state p = (ε−ρ)/3 to capture the transition from hot to cold fluid. The cooling solution is reduced to an autonomous dynamical system in the variables (z, η) (where z = qδ ξ g/[2(4−k)] and η = f√g/h), and the sonic point is identified where the forward-characteristic drift D+ vanishes. Smooth crossing of the sonic line imposes a regularity condition that selects the physical branch.

Extended reading notes

Core claim

The paper establishes that the BM solution near the shock is causally protected from the deep interior. For density index k < k_g ≈ 2.062, forward acoustic characteristics launched from the first breakdown point of the BM solution reach the shock only after the blast wave has become non-relativistic. For k_g < k < 4, the paper constructs a new self-similar cooling solution, matched asymptotically to the hot BM shell. This cooling solution terminates at a sonic point (either at finite sound speed or in the zero-pressure limit), which acts as an acoustic barrier: signals from beyond cannot propagate forward to the shock. The result is that the BM near-shock profiles are determined solely by the shock jump conditions, with no dependence on the unknown interior.

Load-bearing premise

The cooling solution uses the equation of state p = (ε − ρ)/3, which fixes the adiabatic index at 4/3 and adds rest-mass inertia but does not capture the full trans-relativistic transition to a 5/3 adiabatic index as the gas cools, potentially affecting the location and nature of the sonic-point barrier.

Editorial extensions

If this is right

  • The BM solution is causally self-consistent, closing a long-standing gap in its derivation for all k < 4.
  • The near-shock profiles are determined by shock conditions alone, independent of the structure of the deep interior.
  • A fully analytic example at k = 7/2 (g = ξ⁻¹, f = ξ⁻², h = 1) makes the causal shielding mechanism explicit and provides a benchmark for numerical work.
  • For gamma-ray burst afterglow modeling, the BM solution can be used confidently for the early ultra-relativistic phase, as the interior cannot affect the shock or the hot shell.

Reading between the lines

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

  • The causal barrier mechanism may apply to other self-similar shock solutions where downstream assumptions break down, potentially simplifying a class of otherwise ambiguous problems.
  • The cooling-layer profiles derived here could be used to compute synchrotron emission from the trans-relativistic layer, which may affect afterglow light curves at intermediate times.
  • The discontinuity in the cooling solution at k = 3 and k_c suggests that the afterglow from steep external media (e.g., a stellar wind) may have qualitatively different properties as the blast wave decelerates.
Original abstract

The Blandford--McKee (BM) solution describes the ultra-relativistic self-similar flow behind a strong spherical blast wave propagating into an external density profile $\propto r^{-k}$, where $k<4$ and $r$ is the distance from the center, but it applies only to the hot shell adjacent to the shock, sufficiently deep behind it, the fluid leaves the BM regime. Since the similarity profiles are determined solely by the shock conditions, with no conditions imposed at the inner end, the validity of the solution near the shock depends on whether the flow beyond the BM regime can imprint on this hot shell. For shallow density profiles, $k<k_g\simeq 2.062$, the characteristic structure of the BM solution already prevents forward-going acoustic information from the non-BM interior from reaching the shock during the ultra-relativistic stage. For steeper profiles, $k_g<k<4$, the hot equation of state fails while the flow is still relativistic and shock-connected. There, we derive a new self-similar solution for the cooling relativistic flow, which has its own similarity scale, equations, and characteristic structure. It overlaps with the hot BM solution toward the shock, and reaches a sonic point beyond this overlap. This critical point separates the shock-connected BM-plus-cooling composite from the deeper downstream region, so acoustic signals generated beyond it cannot propagate toward the shock. The mechanism is explicit at $k=7/2$, for which the cooling solution is obtained analytically. This causal structure is therefore what enables the BM solution, which remains self-consistent near the shock even though it does not globally describe the full downstream flow.

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.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,132 input tokens.

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

The abstract, sentence by sentence

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

The Blandford--McKee (BM) solution describes the ultra-relativistic self-similar flow behind a strong spherical blast wave propagating into an external density profile $\propto r^{-k}$, where $k<4$ and $r$ is the distance from the center, but it applies only to the hot shell adjacent to the shock, sufficiently deep behind it, the fluid leaves the BM regime. Since the similarity profiles are determined solely by the shock conditions, with no conditions imposed at the inner end, the validity of the solution near the shock depends on whether the flow beyond the BM regime can imprint on this hot shell. For shallow density profiles, $k<k_g\simeq 2.062$, the characteristic structure of the BM solution already prevents forward-going acoustic information from the non-BM interior from reaching the shock during the ultra-relativistic stage. For steeper profiles, $k_g<k<4$, the hot equation of state fails while the flow is still relativistic and shock-connected. There, we derive a new self-similar solution for the cooling relativistic flow, which has its own similarity scale, equations, and characteristic structure. It overlaps with the hot BM solution toward the shock, and reaches a sonic point beyond this overlap. This critical point separates the shock-connected BM-plus-cooling composite from the deeper downstream region, so acoustic signals generated beyond it cannot propagate toward the shock. The mechanism is explicit at $k=7/2$, for which the cooling solution is obtained analytically. This causal structure is therefore what enables the BM solution, which remains self-consistent near the shock even though it does not globally describe the full downstream flow.

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 The Blandford--McKee (BM) solution describes the ultra-relativistic self-similar flow behind a strong spherical blast wave propagating into an external density profile $\propto r^{-k}$, where $k<4$ and $r$ is the distance from the center, but it applies only to the hot shell adjacent to the shock, sufficiently deep behind it, the fluid leaves the BM regime. background 0.99 0.10 0.08 asserted 1.00
s1 Since the similarity profiles are determined solely by the shock conditions, with no conditions imposed at the inner end, the validity of the solution near the shock depends on whether the flow beyond the BM regime can imprint on this hot shell. background 0.59 0.13 0.06 asserted 0.90
s2 For shallow density profiles, $k<k_g\simeq 2.062$, the characteristic structure of the BM solution already prevents forward-going acoustic information from the non-BM interior from reaching the shock during the ultra-relativistic stage. result 0.99 0.46 0.18 asserted 1.00
s3 For steeper profiles, $k_g<k<4$, the hot equation of state fails while the flow is still relativistic and shock-connected. result 0.85 0.54 0.10 asserted 1.00
s4 There, we derive a new self-similar solution for the cooling relativistic flow, which has its own similarity scale, equations, and characteristic structure. contribution 0.97 0.93 0.07 asserted 1.00
s5 It overlaps with the hot BM solution toward the shock, and reaches a sonic point beyond this overlap. result 0.86 0.83 0.07 asserted 1.00
s6 This critical point separates the shock-connected BM-plus-cooling composite from the deeper downstream region, so acoustic signals generated beyond it cannot propagate toward the shock. result 0.51 0.81 0.23 asserted 1.00
s7 The mechanism is explicit at $k=7/2$, for which the cooling solution is obtained analytically. method 0.36 0.85 0.07 asserted 1.00
s8 This causal structure is therefore what enables the BM solution, which remains self-consistent near the shock even though it does not globally describe the full downstream flow. interpretation 0.96 0.79 0.84 asserted 1.00

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. FIG. 1. Shock Lorentz factor,
Figure 2. FIG. 2. Physical roots of the finite-sound-speed regularity condition,
Figure 5. FIG. 5. Location of the inner sonic point,

Showing 3 of the paper's 5 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.00116 by the authors.

References

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.