Pith. sign in

REVIEW 4 major objections 2 references

Non-detections of receding relativistic jets cut parameter uncertainty by over 40% and expose strong observational bias against fast off-axis ejecta.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-12 14:07 UTC pith:JCVYI2BB

load-bearing objection Solid methodological paper: non-detection likelihood penalty tightens posteriors ~40% on MAXI J1535 and rest-frame maps show clear high-Γ bias; symmetry assumption is load-bearing but standard and partially checked on the two-sided source. the 4 major comments →

arxiv 2606.13806 v1 pith:JCVYI2BB submitted 2026-06-11 astro-ph.HE astro-ph.IM

Observational Biases and Improved Modelling of Off-axis Relativistic Jets

A. J. Cooper , A. P. Scott , L. Rhodes , F. Carotenuto , A. K. Hughes , J. H. Matthews , K. Savard , F. J. Cowie
show 3 more authors
E. L. Elley C. Lilje R. Fender
This is my paper
classification astro-ph.HE astro-ph.IM
keywords relativistic jetsDoppler boostingX-ray binariesoff-axis jetsobservational biaskinematic modellingnon-detectionsLorentz factor
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Relativistic Doppler boosting hides or dims large parts of the jet population that observers would otherwise detect. Using two Galactic X-ray binary case studies, the authors show that the mere absence of a receding jet can be turned into a quantitative constraint that shrinks the posterior uncertainties on jet parameters by more than 40 percent. They further recover a common rest-frame power-law light curve for the approaching and receding components of MAXI J1820+070 and use that template to map the full initial-Lorentz-factor versus viewing-angle plane, revealing that present radio strategies systematically miss ejecta launched at Lorentz factors greater than about 5 and most receding components. The practical message is that early-time, high-resolution observations together with sensitive late-time monitoring, plus deliberate use of non-detections in the likelihood, are required to remove the bias and to prepare for the coming flood of off-axis gravitational-wave and optical transients.

Core claim

Incorporating the non-detection of a receding jet into kinematic modelling of the one-sided ejecta of MAXI J1535-571 reduces average 1-sigma posterior uncertainties by more than 40 percent, while the rest-frame emission recovered from the double-sided jets of MAXI J1820+070 follows a common power-law decay that, when mapped across parameter space, demonstrates that current observing cadences strongly suppress detections of high-Gamma0 and receding components.

What carries the argument

Likelihood penalisation that converts predicted receding-jet flux into a soft constraint whenever that flux would have exceeded a realistic signal-to-noise and angular-separation threshold (Eq. 7), together with Doppler de-boosting that recovers a shared rest-frame power-law light curve used as a template across the Gamma0-theta_v plane.

Load-bearing premise

The two jets are assumed to be intrinsically identical so that only the Doppler factor differs; if the approaching and receding sides have different energies or ambient densities the conversion and the likelihood penalty both fail.

What would settle it

A second double-sided X-ray binary jet whose approaching and receding rest-frame light curves, after Doppler correction, do not lie on a common power-law, or a kinematic fit of a one-sided source in which adding the non-detection penalty does not shrink the posteriors.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 0 minor

Summary. The paper studies how relativistic Doppler boosting biases observations and modelling of discrete X-ray binary jet ejecta, using two case studies. For the one-sided ejecta of MAXI J1535−571, nested-sampling kinematic fits that penalise parameter sets predicting a detectable receding component (Eqs. 3–7; Figs. 2–3) reduce average 1σ posterior uncertainties by ~40% relative to kinematics alone. For the bipolar ejecta of MAXI J1820+070, deboosting and rest-frame time correction of both components (using the Carotenuto et al. 2024 deceleration profile) yield a common power-law rest-frame light curve (Eq. 12; Fig. 5). That template is then mapped across the Γ0–θv plane to produce detectability maps (Figs. 6–7, A2), arguing that current strategies miss high-Γ0 (≳5) approaching ejecta and receding components over much of parameter space, and that early- and late-time high-resolution observations plus non-detection constraints improve modelling of off-axis jetted transients.

Significance. If the methodology holds, the work supplies a practical, immediately usable way to fold receding-jet non-detections into kinematic inference and a concrete map of observational bias in the Γ0–θv plane for XRB-like ejecta. The side-by-side nested-sampling posteriors (Figs. 2–3), the collapse of both J1820 components onto a single rest-frame power law (Fig. 5), and the explicit detectability contours (Figs. 6–7, A2) are reproducible and falsifiable contributions. The framing toward GW-triggered and orphan off-axis afterglows is timely. The main numerical claims (∼40% uncertainty reduction; strong bias against Γ0 ≳ 5) rest on the intrinsic symmetry of approaching and receding sides and on a fixed Doppler exponent; quantifying robustness to those choices would make the result more durable for the community.

major comments (4)
  1. Sections 2.1 and 3.2 (Eqs. 3–5, 7, 11–12): both the J1535 likelihood penalty and the J1820 rest-frame template / Γ0–θv maps assume intrinsically identical approaching and receding jets (same E0, n0, opening angle, and rest-frame emission history), so that Doppler factor and light-travel time alone convert one light curve into the other. Fig. 5’s alignment is obtained under that same assumption and therefore does not independently validate it for the one-sided case or for the bias maps. A load-bearing sensitivity test is needed: e.g. allow O(1) side-to-side differences in energy or ambient density, and/or vary the Doppler exponent in the documented range ∼2–3 (and the GRS 1915 estimates cited), and show how the 40% uncertainty reduction and the detectability contours in Figs. 6–7 and A2 change.
  2. Section 2.1, Eq. 7 and surrounding text: the non-detection penalty for MAXI J1535 is applied to synthesised receding fluxes at the observer times of the approaching detections, with fixed SNR>8 (Fnoise=10 μJy) and θlim from a constant 5″ beam (Eq. 6), explicitly without requiring that real observations existed at those receding arrival times. That choice makes the reported 41.5% (49.7% excluding Γ0) uncertainty reduction specific to an idealised observing cadence rather than to the actual ATCA/MeerKAT campaign. Either re-run the penalty only at epochs when the source was observed (or when a receding component would have fallen in the FoV/beam), or clearly reframe the 40% figure as an upper bound under continuous monitoring and show the degradation when real epochs are used.
  3. Section 3.2–3.3 and Eq. 12: the rest-frame power-law template is fitted only for t_rest ≥ 25 d and then forced to zero flux at earlier rest-frame times when generating synthetic light curves. The paper notes that reverse-shock crossing depends on Γ0, E0 and n0, so a fixed 25 d cut-off is not self-consistent across the Γ0–θv grid. Because early rest-frame emission maps to late observer times for the receding jet (and to early times for low-θv approaching jets), this cut-off materially shapes the hatched non-detectable regions in Figs. 6–7 and A2. Provide at least one alternative (earlier cut-off, reverse-shock rise, or Γ0-dependent crossing time) and show which qualitative bias conclusions survive.
  4. Section 3.1 and Appendix Fig. A1: nested sampling leaves Γ0 essentially unconstrained for J1820 kinematics alone, while the same-separation flux argument yields only Γ(r_sep)≈1.28 with large systematic uncertainty from non-identical separations and frequencies. The subsequent bias maps nonetheless adopt the Carotenuto et al. (2024) Γ0≈2.6 and E_eff as the baseline. State explicitly how the detectability conclusions change if Γ0 is drawn from the broad nested-sampling posterior rather than fixed at the MCMC best-fit, or justify why the MCMC value is preferred for the template.

Circularity Check

0 steps flagged

No significant circularity: non-detection penalties and Γ0–θv bias maps are genuine model applications, not results forced by definition from their inputs.

full rationale

The paper’s two central numerical claims do not reduce to their inputs by construction. For MAXI J1535, the kinematic likelihood (Eq. 2) is first fit to approaching-jet separations alone; a separate penalty (Eq. 7) is then applied only when the Doppler-mapped receding prediction would have exceeded stated SNR and resolution thresholds. The reported ~40% shrinkage of posterior widths is an empirical comparison of two nested-sampling runs, not a tautology. For MAXI J1820, rest-frame fluxes are recovered by applying kinematic Γ(t) profiles (from position data) to independent flux measurements; the subsequent power-law fit (Eq. 12) is used only as a template for hypothetical jets at other (Γ0, θv). That is standard forward modelling, not a fitted quantity re-labelled as a prediction of the same data. The mild self-consistency that both de-boosted light curves land on a common power law is acknowledged by the authors as relying on the symmetry assumption already used in the joint kinematic fit; it is presented as a consistency check, not as an independent derivation. Self-citations (Cooper et al. 2025 for fixed θc, n0; Carotenuto et al. 2024 for the baseline kinematic solution) supply fixed parameters or prior results but are not load-bearing uniqueness theorems that force the present conclusions. The symmetry assumption itself is a modelling premise whose validity is a correctness/robustness question, not a circularity. Overall circularity is negligible.

Axiom & Free-Parameter Ledger

6 free parameters · 4 axioms · 0 invented entities

The central numerical claims rest on standard relativistic kinematics plus a small set of domain assumptions (jet symmetry, fixed Doppler exponent, simplified observability thresholds) and a handful of free parameters that are either fitted or fixed from earlier work. No new physical entities are introduced.

free parameters (6)
  • Γ0 (initial bulk Lorentz factor)
    Free parameter with uniform prior [1,100]; remains poorly constrained by kinematics alone and is only bounded from below by the non-detection penalty.
  • E_iso (isotropic-equivalent energy)
    Log-uniform prior; posterior becomes less bimodal once receding-jet constraints are added.
  • θv (viewing angle)
    Cosine prior; high-θv cut-off appears only after non-detection information is included.
  • D (source distance)
    Truncated normal prior centred on HI absorption measurement; low-distance cut-off enforced by receding-jet non-detection.
  • S0, γ (rest-frame power-law amplitude and index)
    Fitted to deboosted MAXI J1820 light curves (log10 S0 = 1.566 ± 0.07, γ = −3.579 +0.13/−0.14) and then used as the universal template for the bias maps.
  • SNR threshold and beam size
    Hand-chosen observability criteria (SNR > 8, 5-arcsec beam) that define when a predicted receding jet would have been detected.
axioms (4)
  • domain assumption Approaching and receding jets are intrinsically identical in energy, opening angle, ambient density and rest-frame emission history.
    Required to convert observed approaching fluxes into predicted receding fluxes via the Doppler factor alone (Sections 2.1, 3.2).
  • domain assumption Observed flux scales as δ^{3−α} with α ≈ −0.5 (or 2–3 range).
    Standard discrete-ejecta boosting formula (Eq. 3); the precise exponent is not marginalised.
  • domain assumption Blast-wave energy conservation and conical geometry (Huang et al. 1999 style) correctly describe the deceleration profile.
    Underpins both the jetsimpy kinematics and the Carotenuto et al. (2024) profiles used for deboosting.
  • ad hoc to paper Early-time (t_rest < 25 d) emission can be set to zero when extrapolating the power-law template.
    Motivated by reverse-shock crossing arguments but is a modelling choice that affects predicted early receding fluxes.

pith-pipeline@v1.1.0-grok45 · 22007 in / 3045 out tokens · 30470 ms · 2026-07-12T14:07:47.904093+00:00 · methodology

0 comments
read the original abstract

Relativistic Doppler boosting significantly affects the observed emission of astrophysical jets resulting in observational biases. In this work we investigate the observational biases and modelling opportunities which arise due to relativistic boosting using two X-ray binary case studies. Using the one-sided jet ejecta from MAXI J1535-571, we demonstrate that incorporating non-detections of the receding jet ejecta into kinematic modelling can significantly improve parameter estimation, reducing posterior uncertainties by over 40%. For the bipolar jets of MAXI J1820+070, we recover the intrinsic jet rest-frame emission of both approaching and receding jet components, demonstrating that they follow a common powerlaw evolution. Using this rest-frame emission profile as a base model, we show that current observational strategies strongly bias against detecting ejecta with high initial Lorentz factors >5 and receding ejecta components across a broad region of parameter space. These results highlight the importance of observational strategy selection, particularly early-time and late-time observations, and leveraging non-detections in the modelling of relativistic jets. More generally, quantifying observational biases and maximising modelling capabilities by incorporating the non-detection of receding jets can be employed to enhance interpretation of future gravitational-wave/optically-triggered observations of off-axis, extragalactic jetted transients.

Figures

Figures reproduced from arXiv: 2606.13806 by A. J. Cooper, A. K. Hughes, A. P. Scott, C. Lilje, E. L. Elley, F. Carotenuto, F. J. Cowie, J. H. Matthews, K. Savard, L. Rhodes, R. Fender.

Figure 1
Figure 1. Figure 1: Approaching (top panel) and predicted receding (bottom two pan￾els) jet flux from the initial kinematic fit of MAXI J1535, coloured by fre￾quency. As expected, the receding jet flux is lower and delayed as compared to approaching jet flux for both low-Γ0 and high-Γ0 solutions. The underlying data are published in Russell et al. (2019). given parameter set: Lpen (𝜃) = ln L (𝜃) − ∑︁ 𝑖∈I " max 0, 𝑓 rec 𝑖 (𝜃) … view at source ↗
Figure 4
Figure 4. Figure 4: Lorentz factor observer-frame temporal profile of the approaching (blue) and receding (red) components of MAXI J1820 jet ejecta (Bright et al. 2020) from kinematic modelling presented in Carotenuto et al. (2024). The kinematic data of both jets was fit by Carotenuto et al. (2024) utilising a trans-relativistic, conical blastwave model (e.g. constant half-opening angle, 𝜙jet), in which the ejecta shell dece… view at source ↗
Figure 2
Figure 2. Figure 2: Corner plot of posterior parameter space of kinematic-only fit to the separation data of MAXI J1535 without receding jet constraints. Orange lines show the best-fit parameters determined by the maximum likelihood, and blue contours highlight 1𝜎, 2𝜎, and 3𝜎 credible regions. Here, 𝐸iso = 16𝐸0/𝜃 2 c ≈ 104𝐸0 for our fixed value of 𝜃c = 2.25 deg [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Corner plot of posterior parameter space of kinematic-only fit to MAXI J1535 including receding jet constraints [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Panel 1: Observed fluxes of approaching jet (blue colours), receding jet (red colours), and unresolved AMI-LA (dark blue) from MAXI J1820. Panel 2: Same data scaled to a common frequency of 1.28 GHz. Panel 3: Same data corrected for (de)boosting of the flux due to bulk relativistic motion assuming the decelerating blast wave model fit. Panel 4: Same detections corrected to rest-frame emission time, with th… view at source ↗
Figure 7
Figure 7. Figure 7: Predicted MAXI J1820-like emission from the approaching (left panels) and receding (right panels) as a function of time and Γ0 for three view￾ing angles (𝜃v = [15, 45, 75] deg). Note the powerlaw colour scale changes for each row. Cyan dashed contours again correspond to 𝐹𝜈 = 150 𝜇Jy flux, with the white hatched region excluding regions of lower flux. The shape of the observable region shows how observatio… view at source ↗
Figure 6
Figure 6. Figure 6: Predicted MAXI J1820-like emission from the approaching (left panels) and receding (right panels) as a function of the initial Lorentz factor Γ0 and viewing angle 𝜃v for 50, 100, 200, and 400 days post-launch in the observer frame. Note that the colour scale changes for each row. Cyan dashed contours correspond to 𝐹𝜈 = 150 𝜇Jy flux, with the white hatched region excluding regions of lower flux. The green c… view at source ↗

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 · 1 linked inside Pith

  1. [1]

    D., et al., 2018, ApJ, 863, L18 Andreoni I., et al., 2019, PASP, 131, 068004 Atri P., et al., 2020, MNRAS, 493, L81 Bacon N

    Alexander K. D., et al., 2018, ApJ, 863, L18 Andreoni I., et al., 2019, PASP, 131, 068004 Atri P., et al., 2020, MNRAS, 493, L81 Bacon N. J., Cooper A. J., Kantzas D., Matthews J. H., Fender R., 2026, MNRAS, 546, stag080 Bahramian A., et al., 2023, ApJ, 948, L7 Blandford R. D., Payne D. G., 1982, MNRAS, 199, 883 Blandford R. D., Znajek R. L., 1977, MNRAS,...

    Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/0303339 2018

  2. [2]

    MNRAS000, 1–10 (2026)

    FigureA2.DetectabilityofMAXIJ1820-likeejectaat60,120,200,and280 days post launch at1.28GHz, as a function of the initial Lorentz factorΓ0 andthe viewingangle𝜃 v.Dark bluecorresponds toregions ofthe parameter spacewherebothjetsaresimultaneouslydetectable,lighterblueonlytheap- proachingjet,andgreenonlytherecedingjet.Theyellowregioncorresponds to the case wh...

    2026