Pith. sign in

REVIEW 4 major objections 6 minor 45 references

The magnetic field along SS 433's inner jet declines as H^-0.50, slowly enough that magnetic rigidity grows downstream and protons can reach ~2-3 PeV within a few hundred AU.

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 →

T0 review · deepseek-v4-flash

2026-08-04 09:27 UTC pith:JRKQJARS

load-bearing objection A genuinely new resolved B(H) measurement for SS 433's inner jet, but the 'rigidity grows with H' headline is an artifact of assuming a cone the data contradict. the 4 major comments →

arxiv 2608.02314 v1 pith:JRKQJARS submitted 2026-08-03 astro-ph.HE

Magnetic rigidity reveals the PeVatron acceleration region in SS433

classification astro-ph.HE
keywords SS 433PeVatronHillas limitmagnetic rigiditymicroquasar jetsequipartition magnetic fieldultra-high-energy gamma raysjet magnetic field profile
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.

This paper tries to locate the hidden PeVatron in SS 433, a microquasar whose >100 TeV gamma-ray emission requires protons accelerated to petaelectronvolt energies. Using archival multi-frequency radio interferometric images taken on a single night, the authors resolve the eastern jet into ten transverse slices spanning about 24-75 AU from the central object. From the local synchrotron luminosity, volume, and spectral index they derive an equipartition magnetic field of roughly 0.20-0.35 G that declines as B(H) ∝ H^-0.50±0.12. Because this is slower than 1/H, the speed-weighted Hillas energy B(H)H grows downstream, reaching ~2.2 PeV at 100 AU and crossing the ~2.6 PeV proton cutoff inferred from the >100 TeV gamma-ray spectrum at ~140 AU. The upshot is that the AU-scale baryonic ejecta, not the parsec-scale TeV-emitting lobes, are the region where PeV protons are accelerated and injected, with the radiation occurring later in the surrounding gas-rich W50 environment.

Core claim

The paper's central claim is that the compact radio jet of SS 433 itself clears the PeV bar: the equipartition magnetic field along the inner (eastern) jet falls as B(H) ∝ H^-0.50±0.12 over 24-75 AU, with normalization B100 ≈ 0.19 G at 100 AU. For a conical jet with effective acceleration radius R_acc = α_j H (α_j ≈ 0.1) and jet speed β_j ≈ 0.26c, the speed-weighted Hillas limit E_max ≈ β_j Z e B R_acc then rises downstream as H^0.50, reaching ≈2.2 PeV at 100 AU and crossing the 2.63 PeV cutoff (from a lepto-hadronic fit to the >30 TeV spectrum) at ≈140 AU. The authors therefore interpret the tens-to-hundreds of AU scale baryonic ejecta as the hidden PeV-hadron acceleration and injection reg

What carries the argument

The load-bearing identity is the speed-weighted Hillas condition E_max ≃ β_j Z e B R_acc, combined with the geometric closure R_acc ≃ α_j H for a conical jet, so E_max ∝ B(H)H. The magnetic field profile B(H) = B100 (H/100 AU)^(-q) with q = 0.50±0.12 is obtained from the classical equipartition/minimum-energy formula, which turns a measured synchrotron luminosity, emitting volume (from beam-deconvolved jet width), and spectral index into a field strength; the q < 1 condition is what makes B(H)H grow downstream and turns the inner jet into a PeV-capable accelerator. The measured width growth W ∝ H^0.44±0.12 (close to a toroidal field scaling B ∝ W^-1) is a supporting element.

Load-bearing premise

The load-bearing premise is that the equipartition (minimum-energy) magnetic field estimated from radio synchrotron luminosity equals the actual field that confines protons, and that the fitted power-law decline measured over 24-75 AU continues unchanged out to ~100-140 AU.

What would settle it

A future radio observation resolving SS 433's jet at deprojected distances of ~100-140 AU that finds B(H) falling more steeply than H^-1, or an independent measurement (e.g. via Faraday rotation or core shift) giving B100 well below 0.19 G, would falsify the claim that the Hillas energy crosses the ~2.6 PeV cutoff by 140 AU.

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

If this is right

  • The acceleration site for the >100 TeV protons in SS 433 is the inner baryonic jet on tens-to-hundreds of AU scales, not the parsec-scale TeV lobes.
  • The >100 TeV gamma rays are produced later, after the PeV protons escape or are advected into W50's neutral gas, naturally explaining why the UHE emission is centrally concentrated and overlaps the H I cloud.
  • The magnetic field configuration in the inner jet is predominantly toroidal (scaling as width^-1), which future jet-launching and collimation models must respect.
  • Only a small fraction (roughly 1-4% for k=20-60) of the one-sided jet power needs to be carried locally in fields and non-thermal particles, so the compact PeVatron is energetically comfortable.
  • PeV electrons are ruled out as the source of the >100 TeV photons from the compact ejecta, since they would cool in ~110 seconds; a hadronic injection picture is required.

Where Pith is reading between the lines

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

  • If the shallow B ∝ H^-0.5 decline is generic in microquasar jets, the Hillas ceiling should grow downstream in any conical, mildly relativistic jet; this predicts a correlation between the observed gamma-ray cutoff energy and the jet's opening angle.
  • The crossing at ~140 AU rests on extrapolating a power law fitted to 24-75 AU; a steeper decline beyond 75 AU (e.g. due to recollimation or field reconnection) would move the crossing farther out or erase it, so next-generation radio imaging of the 100 AU scale is a sharp test.
  • An independent field measurement at the same radii-available through core-shift analysis or Faraday rotation of the polarized jet-would either confirm the 0.19 G normalization or expose the equipartition assumption's limitations.
  • If this compact injection zone feeds the W50 hadronic reservoir, the same protons should also produce neutrinos and a characteristic pion-decay bump below the >100 TeV photons; deep gamma-ray/neutrino observations of the central region could test that multi-zone picture.

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 / 6 minor

Summary. The paper analyzes archival VLBA+Y1 observations of SS 433 at 5, 8.4, and 15 GHz. From ten transverse slices of the eastern jet spanning deprojected distances of about 24–75 AU, it measures flux densities, beam-deconvolved widths, and spectral indices, and under minimum-energy assumptions derives B_eq(H) ∝ H^{-0.50±0.12}. Using a conical prescription R_acc = α_j H with α_j = 0.1, it converts this profile into a speed-weighted Hillas energy E_max ∝ B H that grows downstream, reaching ~2.2 PeV at 100 AU and crossing the LHAASO-inferred proton cutoff E_p,cut ≈ 2.63 PeV at ~140 AU. The paper then fits the TeV–UHE spectrum with lepto-hadronic and pure hadronic models, argues that the compact ejecta are not the gamma-ray radiation zone, and proposes a multi-zone picture in which PeV protons are accelerated in the AU-scale baryonic ejecta and later produce the extended >100 TeV emission via pp interactions in the W50 environment. A projected diffusion model is used to argue that quasi-isotropic to mildly anisotropic transport of such protons into the HI reservoir matches the LHAASO morphology.

Significance. If the central claim held, the paper would identify the long-sought acceleration site of PeV protons in SS 433 and establish AU-scale microquasar jets as PeVatron injection regions. The observational work is valuable: it uses public multi-frequency VLBA data, propagates uncertainties with Monte Carlo, explicitly explores the k-dependence of the equipartition normalization, and includes energy-budget and acceleration-timescale checks. The authors are also unusually candid about some model dependencies. However, the advertised "rigidity grows with H" result is not supported by the paper's own measured jet widths: W ∝ H^0.44 implies B_eq W is roughly flat, so the PeV crossing at 140 AU disappears under the measured transverse scale. The paper's sensitivity analysis already contains the numbers that contradict its headline. The work can be reframed as a credible, more limited claim about PeV-scale confinement under motivated baryon loading, but as written the central conclusion is overstated and needs major revision.

major comments (4)
  1. [§3.3, Eq. (10); Table 1] The central claim that magnetic rigidity grows downstream rests on the conical closure R_acc = α_j H with α_j = 0.1. But the paper's own resolved widths give W ∝ H^{0.44±0.12} (§3.1, Table 1), so the jet is collimating, not conical. Using the measured transverse scale R_acc = W/2, the tabulated values give B_eq × W ≈ 400±50 mG mas with no trend, i.e. B_eq W ∝ H^{-0.04±0.17}. Equation (14) then yields E_max ∝ H^{-0.06}, not H^{0.50}. Numerically, with the paper's own W_{100}=13.1 AU extrapolation, B_{100}=0.19 G and β_j=0.26, the speed-weighted Hillas energy is ~1.45 PeV at 100 AU, not 2.2 PeV, and the 2.63 PeV crossing at ~140 AU disappears. The manuscript acknowledges this near Eq. (17), but the abstract, title, and Section 2 present the conical case as the result. Since this is the load-bearing step for the claimed identification of the PeVatron region, the measured-width case must be
  2. [§3.1, Eq. (2); §3.3, Eq. (13)] The equipartition field B_eq is a minimum-energy reference, not a direct measurement of the confining field. The normalization B_{100}=0.19 G assumes k=1, φ=1, and a specific cylindrical volume, and the quoted ±0.12 slope error does not include this normalization systematic. Because E_max and the crossing distance scale linearly with B (Eqs. 14–16), the quantitative PeV claim inherits this systematic. The text itself calls B_eq a 'reference scale rather than a direct field measurement', yet Eqs. (13)–(17) treat it as the actual field. The motivated baryon-loading range k=20–60 (B=0.37–0.50 G) is the appropriate benchmark for the PeV claim and should be presented as such rather than as a secondary sensitivity test.
  3. [§3.2, lepto-hadronic fit] The adopted proton cutoff E_p,cut = 2.63 PeV comes from a lepto-hadronic fit with BIC = 13.65 for the >30 TeV range, while the pure hadronic fit has BIC = 9.76 and E_p,max = 2.34 PeV. The paper notes this but adopts the higher-BIC model as fiducial. The claimed consistency with the LHAASO-inferred cutoff is sensitive to this choice. The conclusion should be framed for the full range of E_p,max across the fitted models (central values roughly 1.6–3.1 PeV) or the model preference should be justified from physical priors. This does not destroy the PeV-scale conclusion, but it affects the specific equality/crossing language used in the abstract and Section 2.
  4. [§3.3, Eq. (13)] The fit B(H) ∝ H^{-0.50±0.12} is anchored only at 24–75 AU, but the PeV crossing and the headline numbers rely on extrapolation to 100–140 AU. No data constrain B or W in this extrapolated range. Given that the measured width is still expanding as W ∝ H^{0.44±0.12} over the fitted range, the conical extrapolation to 140 AU is especially fragile. The 100–140 AU energies should be explicitly labeled as extrapolations, and the authors should test whether the conclusion survives at the extremes of the width and slope uncertainties.
minor comments (6)
  1. [Abstract and title] The phrase 'rigidity grows with H' is only true under the fiducial conical closure; as written it overstates what the data show. Please qualify it as a geometry-dependent scenario, or reframe around the measured width result.
  2. [§3.1; Table 1] The text says 'ten consecutive slices spanning H ≈ 3.5–12 mas', but Table 1 lists 3.56–11.51 mas. Please align the quoted range with the table.
  3. [Eq. (3)] The analytic expression for c12 has a removable singularity at α = 0.5. The linear interpolation between α = 0.44 and 0.56 is mentioned, but the implementation should be documented in more detail, including the interpolation grid.
  4. [Figure 2] The k = 0, 1, 100 points are nearly indistinguishable in the figure. Consider showing a single curve with the shaded k = 0–100 range, which would make the k-independence of the slope clearer.
  5. [References] The astrometric alignment and the collimation slope rely on Yan et al. [21], which is in press and shares authors with this paper. The dependence should be stated more explicitly, and the authors should confirm the final published calibration is used.
  6. [§3.7] The transport calculation is described as a 'conditional morphology diagnostic', which is appropriate. The abstract-level claim of quasi-isotropic morphology should not be stated more strongly than this conditional framing.

Circularity Check

0 steps flagged

No significant circularity: radio B(H) profile and gamma-ray cutoff are independent; the conical-geometry assumption is disclosed and tested, not a fitted input.

full rationale

The derivation chain is not circular. The magnetic-field profile B(H) is derived from VLBA radio observables through the equipartition formula (Eq. 2), using independently measured flux densities, widths, and spectral indices; no gamma-ray data enter this step. The fitted slope B(H) ∝ H^(-0.50±0.12) (Eqs. 6 and 13) is a least-squares fit to ten radio slices. The LHAASO-inferred proton cutoff Ep,cut = 2.63 PeV is obtained from a separate lepto-hadronic fit to the gamma-ray spectrum (Section 3.2), and the Hillas comparison (Eqs. 14–17) only compares B(H)·R_acc·β_j with that cutoff; the cutoff is not fed back into the radio fit or the B-normalization, so there is no fitted-input-called-prediction loop. The conical scale R_acc = α_j H is explicitly introduced as a 'geometric closure' (Eq. 10), and the paper tests the measured-width alternative, noting that the nominal crossing distance is 'geometry dependent' (after Eq. 17). This is a disclosed model assumption and extrapolation, not a circular definition. The citation to Yan et al. [21] is a self-citation in the sense of shared authors, but it supplies an external observational measurement (core-shift alignment, viewing angle, collimation slope) and is not used to assert uniqueness or to smuggle in the conical ansatz; an inaccurate alignment or geometry would be a calibration/correctness risk, not circularity. The morphology transport calculation is explicitly flux-calibrated (Section 3.7) and is not presented as an independent prediction. No load-bearing step reduces to its own input, so the appropriate finding is no significant circularity.

Axiom & Free-Parameter Ledger

4 free parameters · 8 axioms · 0 invented entities

No new particles, fields, or dimensions are introduced. The free parameters are the energy-partition ratio k, the geometric opening factor α_j, the filling factor, and the gamma-fit cutoff used as a comparison threshold. The axioms are dominated by the equipartition mapping and the jet geometry closure.

free parameters (4)
  • k = U_nonrad/U_rad (baryon/electron energy ratio) = 1 (reference), 0-100 bracket
    Multiplicative in Eq. 2; changes Beq normalization by up to a factor of 4 but not the slope; used for energy-budget checks.
  • α_j (semi-opening factor of magnetized jet envelope) = 0.1 (fiducial)
    Geometric closure R_acc=α_jH in Eq. 10; sets the Hillas normalization and the H=140 AU crossing; the paper tests W/2 and a Poynting-bound alternative.
  • Volume filling factor φ = 1
    Adopted in Eq. 2 following the standard equipartition convention; if φ<1, Beq increases as φ^(-2/7), leaving the slope unchanged.
  • Proton cutoff energy Ep,cut (gamma-ray fit) = 2.63 PeV (lepto-hadronic, >30 TeV; 1.64-3.08 PeV across models)
    Fitted to the LHAASO spectrum; used as the target for the Hillas comparison, not as an input to B(H). Broad uncertainties weaken the 'factor 1.2 agreement' claim.
axioms (8)
  • domain assumption Equipartition/minimum-energy conversion (Eq. 2) maps synchrotron luminosity and volume to magnetic field.
    Central for B(H); assumes radiating particles and field share the volume and are near energy minimum; normalization depends on k and filling factor.
  • domain assumption Retained slices are optically thin synchrotron emitters with 0.3 < α < 1.0.
    Justifies Eq. 2; supported by spectral indices in Fig. 6, but requires excluding innermost and outermost slices.
  • ad hoc to paper The effective transverse scale for proton confinement is R_acc = α_jH with α_j = 0.1.
    Explicit 'geometric closure' in Methods 3.3, not directly measured; the paper tests W/2 and Poynting-flux alternatives.
  • domain assumption SS 433's bulk jet speed β_j=0.26, viewing angle 57°, and distance 5.5 kpc from literature.
    Used for deprojection and timescales; sourced from [12,13,21,22].
  • standard math Near-Bohm first-order Fermi acceleration (η_I=1) bounds the acceleration time from below.
    If the diffusion coefficient is larger than Bohm, t_acc grows and PeV acceleration within the flow time may fail.
  • domain assumption The 1998 June 16 observation is representative of the jet's quasi-steady B(H) profile.
    Single epoch; SS 433's jet is precessing and shows time variability.
  • domain assumption The accelerated particles are protons (Z=1), not heavier nuclei.
    The pp gamma-ray interpretation fixes the charge; heavier nuclei would change the rigidity by Z.
  • standard math Hillas criterion r_L < R_acc is the correct confinement condition in a magnetized jet.
    Standard Larmor-radius argument; real confinement also depends on field geometry and diffusion, which the paper treats separately.

pith-pipeline@v1.3.0-daily-deepseek · 25834 in / 19025 out tokens · 161281 ms · 2026-08-04T09:27:43.229924+00:00 · methodology

0 comments
read the original abstract

PeVatrons are cosmic accelerators capable of driving particles to petaelectronvolt (PeV) energies. Recently, microquasar jets have emerged as compelling Galactic PeVatron candidates. This is especially the case for SS 433 as its $>100$ TeV gamma-ray emission is spatially coincident with an atomic cloud. However, the exact region where PeV protons are accelerated and injected within these jets remains unresolved. Here we report, using archival, multi-frequency VLBA observations, the magnetic field profile $B(H)$ along the SS 433 inner jet on tens of AU scale, where $H$ is the distance from the central compact object. We find that the field declines as $B(H) \propto H^{-0.50\pm0.12}$, demonstrating that the magnetic rigidity $B(H)R_{\rm acc}$ grows with $H$ for a conical jet. This implies the Hillas limit ($E_{\rm max} \propto BH$) to lie well beyond a PeV at a few hundred-AU scale, which becomes a highly potential site for accelerating protons to energies $E_{\rm cut} \simeq 2.6$ PeV inferred from the LHAASO gamma-ray spectrum. These results reveal a hidden PeVatron within the baryonic ejecta of microquasar SS 433, well upstream of the extended TeV-emitting lobes.

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

45 extracted references · 3 canonical work pages · 1 internal anchor

  1. [1]

    Astronomy & Astrophysics Review 21, 70 (2013) https://doi.org/10.1007/s00159-013-0070-7 arXiv:1311.7346 [astro- ph.HE]

    Blasi, P.: The origin of galactic cosmic rays. Astronomy & Astrophysics Review 21, 70 (2013) https://doi.org/10.1007/s00159-013-0070-7 arXiv:1311.7346 [astro- ph.HE]

  2. [2]

    Cambridge University Press, Cambridge (2016)

    Gaisser, T.K., Engel, R., Resconi, E.: Cosmic Rays and Particle Physics, 2nd edn. Cambridge University Press, Cambridge (2016). https://doi.org/10.1017/ CBO9781139192194

  3. [3]

    Bell, A.R.: The acceleration of cosmic rays in shock fronts. i. Monthly Notices of the Royal Astronomical Society 182, 147–156 (1978) https://doi.org/10.1093/ mnras/182.2.147

  4. [4]

    The Astrophysical Journal Letters 221, 29–32 (1978) https://doi.org/10.1086/182658

    Blandford, R.D., Ostriker, J.P.: Particle acceleration by astrophysical shocks. The Astrophysical Journal Letters 221, 29–32 (1978) https://doi.org/10.1086/182658

  5. [5]

    Reports on Progress in Physics 46, 973– 1027 (1983) https://doi.org/10.1088/0034-4885/46/8/002

    Drury, L.O.: An introduction to the theory of diffusive shock acceleration of energetic particles in tenuous plasmas. Reports on Progress in Physics 46, 973– 1027 (1983) https://doi.org/10.1088/0034-4885/46/8/002

  6. [6]

    Physics Reports 154, 1–75 (1987) https://doi.org/10.1016/ 0370-1573(87)90134-7

    Blandford, R., Eichler, D.: Particle acceleration at astrophysical shocks: A theory of cosmic ray origin. Physics Reports 154, 1–75 (1987) https://doi.org/10.1016/ 0370-1573(87)90134-7

  7. [7]

    Collaboration: Acceleration of petaelectronvolt protons in the galactic centre

    H.E.S.S. Collaboration: Acceleration of petaelectronvolt protons in the galactic centre. Nature 531, 476–479 (2016) https://doi.org/10.1038/nature17147

  8. [8]

    Nature 594, 33–36 (2021) https://doi.org/ 10.1038/s41586-021-03498-z

    LHAASO Collaboration: Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 gamma-ray galactic sources. Nature 594, 33–36 (2021) https://doi.org/ 10.1038/s41586-021-03498-z

  9. [9]

    A&A 698, 188 (2025) https://doi.org/10.1051/0004-6361/ 202452987 arXiv:2411.08762 [astro-ph.HE]

    Peretti, E., Petropoulou, M., Vasilopoulos, G., Gabici, S.: Particle acceleration and multi-messenger radiation from ultra-luminous X-ray sources: A new class of Galactic PeVatrons. A&A 698, 188 (2025) https://doi.org/10.1051/0004-6361/ 202452987 arXiv:2411.08762 [astro-ph.HE]

  10. [10]

    ApJL 989(2), 25 (2025) https://doi.org/10

    Wang, J., Reville, B., Aharonian, F.A.: Galactic Superaccreting X-Ray Bina- ries as Super-PeVatron Accelerators. ApJL 989(2), 25 (2025) https://doi.org/10. 3847/2041-8213/adf3a4 arXiv:2507.21048 [astro-ph.HE]

  11. [11]

    ApJ 997(1), 128 (2026) https: //doi.org/10.3847/1538-4357/ae246c arXiv:2507.02763 [astro-ph.HE]

    Wan, S.-Y., Wang, J.-s., Liu, R.-Y.: A Leptonic Interpretation of the Ultra-high- energy Gamma-Ray Emission from V4641 Sgr. ApJ 997(1), 128 (2026) https: //doi.org/10.3847/1538-4357/ae246c arXiv:2507.02763 [astro-ph.HE]

  12. [12]

    Astrophysics and Space Physics Reviews 12, 1–152 (2004) https://doi.org/10.48550/arXiv

    Fabrika, S.: The jets and supercritical accretion disk in ss 433. Astrophysics and Space Physics Reviews 12, 1–152 (2004) https://doi.org/10.48550/arXiv. 33 astro-ph/0603390

  13. [13]

    Annual Review of Astronomy and Astrophysics 22, 507–536 (1984) https://doi.org/10.1146/annurev.aa.22.090184

    Margon, B.: Observations of ss 433. Annual Review of Astronomy and Astrophysics 22, 507–536 (1984) https://doi.org/10.1146/annurev.aa.22.090184. 002451

  14. [14]

    Astronomy & Astrophysics 348, 910–916 (1999) arXiv:astro-ph/9907169

    Paragi, Z., Vermeulen, R.C., Fejes, I., Schilizzi, R.T., Spencer, R.E., Stirling, A.M.: The inner radio jet region and the complex environment of ss 433. Astronomy & Astrophysics 348, 910–916 (1999) arXiv:astro-ph/9907169

  15. [15]

    : Very-high-energy particle accel- eration powered by the jets of the microquasar ss 433

    Abeysekara, A.U., Albert, A., Alfaro, R., et al. : Very-high-energy particle accel- eration powered by the jets of the microquasar ss 433. Nature 562, 82–85 (2018) https://doi.org/10.1038/s41586-018-0565-5 arXiv:1810.01892 [astro-ph.HE]

  16. [16]

    Collaboration: Acceleration and transport of relativistic electrons in the jets of the microquasar ss 433

    H.E.S.S. Collaboration: Acceleration and transport of relativistic electrons in the jets of the microquasar ss 433. Science 383, 402–406 (2024) https://doi.org/10. 1126/science.adi2048 arXiv:2401.16019

  17. [17]

    National Science Review 12(12), 496 (2025) https://doi

    LHAASO Collaboration: Ultrahigh-energy gamma-ray emission associated with black hole–jet systems. National Science Review 12(12), 496 (2025) https://doi. org/10.1093/nsr/nwaf496 arXiv:2410.08988

  18. [18]

    Annual Review of Astronomy and Astrophysics 22, 425–444 (1984) https://doi.org/10.1146/ annurev.aa.22.090184.002233

    Hillas, A.M.: The origin of ultra-high-energy cosmic rays. Annual Review of Astronomy and Astrophysics 22, 425–444 (1984) https://doi.org/10.1146/ annurev.aa.22.090184.002233

  19. [19]

    Astrophysics and Space Science 309, 119–125 (2007) https://doi.org/10

    Rieger, F.M., Bosch-Ramon, V., Duffy, P.: Fermi acceleration in astrophysical jets. Astrophysics and Space Science 309, 119–125 (2007) https://doi.org/10. 1007/s10509-007-9466-z arXiv:astro-ph/0610141

  20. [20]

    Annual Review of Astronomy and Astrophysics 49, 119–153 (2011) https://doi

    Kotera, K., Olinto, A.V.: The astrophysics of ultrahigh-energy cosmic rays. Annual Review of Astronomy and Astrophysics 49, 119–153 (2011) https://doi. org/10.1146/annurev-astro-081710-102620 arXiv:1101.4256 [astro-ph.HE]

  21. [21]

    The First Measurement of Jet Collimation Profiles in an X-ray Binary: the Case of SS 433

    Yan, X., Cui, L., Paragi, Z., Frey, S.: The First Measurement of Jet Collimation Profiles in an X-ray Binary: the Case of SS 433. arXiv e-prints, 2607–28122 (2026) https://doi.org/10.48550/arXiv.2607.28122 arXiv:2607.28122 [astro-ph.HE]

  22. [22]

    The Astrophysical Journal Letters 616, 159–162 (2004) https://doi.org/10.1086/426542

    Blundell, K.M., Bowler, M.G.: Symmetry in the changing jets of ss 433: Implica- tions for distance determination. The Astrophysical Journal Letters 616, 159–162 (2004) https://doi.org/10.1086/426542

  23. [23]

    Pacholczyk, A.G.: Radio Astrophysics: Nonthermal Processes in Galactic and Extragalactic Sources. W. H. Freeman, San Francisco (1971)

  24. [24]

    Wiley-VCH, New York (1979) 34

    Rybicki, G.B., Lightman, A.P.: Radiative Processes in Astrophysics. Wiley-VCH, New York (1979) 34

  25. [25]

    Reviews of Modern Physics 42, 237–270 (1970) https://doi.org/10.1103/RevModPhys.42

    Blumenthal, G.R., Gould, R.J.: Bremsstrahlung, synchrotron radiation, and compton scattering of high-energy electrons traversing dilute gases. Reviews of Modern Physics 42, 237–270 (1970) https://doi.org/10.1103/RevModPhys.42. 237

  26. [26]

    Physical Review D 74, 034018 (2006) https://doi.org/10.1103/ PhysRevD.74.034018 arXiv:astro-ph/0606058

    Kelner, S.R., Aharonian, F.A., Bugayov, V.V.: Energy spectra of gamma rays, electrons, and neutrinos produced at proton–proton interactions in the very high energy regime. Physical Review D 74, 034018 (2006) https://doi.org/10.1103/ PhysRevD.74.034018 arXiv:astro-ph/0606058

  27. [27]

    Physical Review D 90, 123014 (2014) https://doi.org/10.1103/PhysRevD.90.123014 arXiv:1406.7369 [astro-ph.HE]

    Kafexhiu, E., Aharonian, F., Taylor, A.M., Vila, G.S.: Parametrization of gamma- ray production cross-sections for pp interactions in a broad proton energy range from the kinematic threshold to pev energies. Physical Review D 90, 123014 (2014) https://doi.org/10.1103/PhysRevD.90.123014 arXiv:1406.7369 [astro-ph.HE]

  28. [28]

    Astronomy & Astrophysics 594, 116 (2016) https:// doi.org/10.1051/0004-6361/201629178 arXiv:1610.06175 [astro-ph.GA]

    HI4PI Collaboration, Ben Bekhti, N., Flöer, L., et al.: Hi4pi: A full-sky h i survey based on ebhis and gass. Astronomy & Astrophysics 594, 116 (2016) https:// doi.org/10.1051/0004-6361/201629178 arXiv:1610.06175 [astro-ph.GA]

  29. [29]

    In: Heck, A

    Greisen, E.W.: AIPS, the VLA, and the VLBA. In: Heck, A. (ed.) Informa- tion Handling in Astronomy - Historical Vistas. Astrophysics and Space Science Library, vol. 285, p. 109 (2003). https://doi.org/10.1007/0-306-48080-8_7

  30. [30]

    In: Hunt, G., Payne, H

    Shepherd, M.C.: Difmap: an Interactive Program for Synthesis Imaging. In: Hunt, G., Payne, H. (eds.) Astronomical Data Analysis Software and Systems VI. Astronomical Society of the Pacific Conference Series, vol. 125, p. 77 (1997)

  31. [31]

    Publications of the Astronomical Society of the Pacific 109(732), 166–172 (1997)

    Condon, J.: Errors in elliptical gaussian fits. Publications of the Astronomical Society of the Pacific 109(732), 166–172 (1997)

  32. [32]

    Astronomische Nachrichten 326, 414–427 (2005) https://doi.org/10.1002/asna

    Beck, R., Krause, M.: Revised equipartition and minimum energy formula for magnetic field strength estimates from radio synchrotron observations. Astronomische Nachrichten 326, 414–427 (2005) https://doi.org/10.1002/asna. 200510366 arXiv:astro-ph/0507367

  33. [33]

    Reviews of Modern Physics 56(2), 255–351 (1984) https://doi.org/10

    Begelman, M.C., Blandford, R.D., Rees, M.J.: Theory of extragalactic radio sources. Reviews of Modern Physics 56(2), 255–351 (1984) https://doi.org/10. 1103/RevModPhys.56.255

  34. [34]

    The Astro- physical Journal 906(2), 105 (2021) https://doi.org/10.3847/1538-4357/abc42d arXiv:2010.14470 [astro-ph.HE]

    Chen, L., Zhang, B.: Analytical solution of magnetically dominated astrophys- ical jets and winds: Jet launching, acceleration, and collimation. The Astro- physical Journal 906(2), 105 (2021) https://doi.org/10.3847/1538-4357/abc42d arXiv:2010.14470 [astro-ph.HE]

  35. [35]

    Blundell, K.M., Laing, R., Lee, S., Richards, A.: SS433’s Jet Trace from ALMA 35 Imaging and Global Jet Watch Spectroscopy: Evidence for Post-launch Par- ticle Acceleration 867(2), 25 (2018) https://doi.org/10.3847/2041-8213/aae890 arXiv:1811.00760 [astro-ph.HE]

  36. [36]

    Nature Astronomy (2026) https://doi.org/10.1038/s41550-026-02922-6

    López-Miralles, J., Perucho, M., Vallés-Pérez, D., Martí, J.-M., Bosch-Ramon, V., Miller-Jones, J.C.A., Motta, S.E., Migliari, S., Marshall, H.L.: Magnetic field topology and colliding discrete ejecta in the precessing jets of ss 433. Nature Astronomy (2026) https://doi.org/10.1038/s41550-026-02922-6

  37. [37]

    PoS(ICRC2015), 922 (2015) arXiv:1509.03319

    Zabalza, V.: naima: A python package for inference of relativistic particle energy distributions from observed nonthermal spectra. PoS(ICRC2015), 922 (2015) arXiv:1509.03319

  38. [38]

    The Annals of Statistics 6(2), 461–464 (1978) https://doi.org/10.1214/aos/1176344136

    Schwarz, G.: Estimating the dimension of a model. The Annals of Statistics 6(2), 461–464 (1978) https://doi.org/10.1214/aos/1176344136

  39. [39]

    Publications of the Astronomical Society of Japan 73(3), 530–545 (2021) https://doi.org/10.1093/pasj/psab018 arXiv:2103.05578 [astro-ph.HE]

    Sakemi, H., Omae, R., Ohmura, T., Machida, M.: Energy estimation of high energy particles associated with the ss433/w50 system through radio observation at 1.4 ghz. Publications of the Astronomical Society of Japan 73(3), 530–545 (2021) https://doi.org/10.1093/pasj/psab018 arXiv:2103.05578 [astro-ph.HE]

  40. [40]

    Physical Review D 78(3), 034013 (2008) https://doi.org/10.1103/PhysRevD.78

    Kelner, S.R., Aharonian, F.A.: Energy spectra of gamma rays, electrons, and neu- trinos produced at interactions of relativistic protons with low energy radiation. Physical Review D 78(3), 034013 (2008) https://doi.org/10.1103/PhysRevD.78. 034013 arXiv:0803.0688 [astro-ph]

  41. [41]

    Astronomy & Astrophysics 669, 149 (2023) https://doi.org/10.1051/0004-6361/ 202243951 arXiv:2210.09390 [astro-ph.HE]

    Fogantini, F.A., García, F., Combi, J.A., Chaty, S., Martí, J., Luque Escamilla, P.: A nustar view of ss 433: Precessional evolution of the jet–disk system. Astronomy & Astrophysics 669, 149 (2023) https://doi.org/10.1051/0004-6361/ 202243951 arXiv:2210.09390 [astro-ph.HE]

  42. [42]

    The Astrophysical Journal 564(2), 941–952 (2002) https://doi.org/10.1086/324398 arXiv:astro-ph/0108206

    Marshall, H.L., Canizares, C.R., Schulz, N.S.: The high-resolution x-ray spectrum of ss 433 using the chandra hetgs. The Astrophysical Journal 564(2), 941–952 (2002) https://doi.org/10.1086/324398 arXiv:astro-ph/0108206

  43. [43]

    Annual Review of Astronomy and Astrophysics 47, 27–61 (2009) https://doi.org/10

    Kalberla, P.M.W., Kerp, J.: The h i distribution of the milky way. Annual Review of Astronomy and Astrophysics 47, 27–61 (2009) https://doi.org/10. 1146/annurev-astro-082708-101823

  44. [44]

    Monthly Notices of the Royal Astronomical Society 396, 1629–1639 (2009) https://doi.org/10.1111/j.1365-2966.2009.14832.x arXiv:0901.4549 [astro- ph.HE]

    Gabici, S., Aharonian, F.A., Casanova, S.: Broad-band non-thermal emission from molecular clouds illuminated by cosmic rays from nearby supernova rem- nants. Monthly Notices of the Royal Astronomical Society 396, 1629–1639 (2009) https://doi.org/10.1111/j.1365-2966.2009.14832.x arXiv:0901.4549 [astro- ph.HE]

  45. [45]

    Physical Review D 88, 023010 (2013) 36 https://doi.org/10.1103/PhysRevD.88.023010 arXiv:1306.3209 [astro-ph.HE] 37

    Giacinti, G., Kachelriess, M., Semikoz, D.V.: Anisotropic cosmic ray diffusion and its implications for gamma-ray astronomy. Physical Review D 88, 023010 (2013) 36 https://doi.org/10.1103/PhysRevD.88.023010 arXiv:1306.3209 [astro-ph.HE] 37