SS 433: Prototypical Galactic Microquasar

• SS 433 is a Galactic microquasar defined by a high-mass X-ray binary with a black hole and persistent baryonic jets.
• Its supercritical accretion powers a thick, edge-on disk with precession and complex wind dynamics that illuminate jet-disk symbiosis.
• Jet interactions with the W50 supernova remnant accelerate particles to multi-TeV energies, offering insights into cosmic-ray origins.

SS 433 is a prototypical Galactic microquasar comprising a high-mass X-ray binary with a compact object, persistent supercritical accretion, and relativistic, precessing baryonic jets. The system serves as the canonical laboratory for super-Eddington accretion, jet-disk symbiosis, and jet-powered particle acceleration at parsec scales, with definitive implications for binary evolution, accretion physics, and cosmic-ray origin.

1. System Architecture and Dynamical Properties

SS 433 consists of a compact object (dynamically confirmed black hole with $M_x\gtrsim 7\,M_\odot$ based on mass-function constraints and radial-velocity measurements) and a Roche-lobe–filling evolved donor (A-type supergiant, $M_d\sim 11$–$15\,M_\odot$) in a 13.082-day eclipsing orbit with inclination $i\simeq 78\degree$–$79\degree$ (Picchi et al., 2020, Bowler, 2020, Cherepashchuk et al., 2019). The system is embedded within the W50 supernova remnant at a distance $d\simeq 4.5$–$5.5$ kpc. The donor overfills its Roche lobe, and mass transfer proceeds at super-Eddington rates, driving a thick, vertically extended accretion disk, with persistent baryonic jets (Paragi et al., 2012, Picchi et al., 2020, Cherepashchuk et al., 2019).

The binary separation $a \approx 0.27$ AU corresponds to $58\,R_\odot$. The donor's Roche radius is $R_L\sim 27\,R_\odot$; spectroscopic broadening and light curves indicate that the A-supergiant fills or slightly exceeds this radius (Picchi et al., 2020). The characteristic mass ratio $M_d\sim 11$0–$M_d\sim 11$1 is empirically constrained by both radial-velocity amplitudes and system stability (Picchi et al., 2020, Cherepashchuk et al., 2019).

2. Accretion Disk, Supercritical Outflow, and Wind Phenomenology

The accretion disk in SS 433 is geometrically and optically thick, viewed nearly edge-on and subject to both 162.375-day precession (half-opening angle $M_d\sim 11$2) and 6.3-day nutation (disc nodding) (Robinson et al., 2017, Paragi et al., 2012). The super-critical mass accretion ($M_d\sim 11$3) powers a persistent, ionization-stratified wind with complex P Cygni profiles, variable on both orbital and precessional cycles. Observed absorption edge velocities span $M_d\sim 11$4 to $M_d\sim 11$5 km\,s$M_d\sim 11$6 by species and ionization zone (Picchi et al., 2020).

Polarimetric and spectroscopic campaigns show that most disk and wind emission lines trace Keplerian and wind kinematics, with the stationary He II $M_d\sim 11$7 and C II $M_d\sim 11$8 lines providing reliable tracers of the compact object's motion, yielding semi-amplitudes $M_d\sim 11$9–$15\,M_\odot$0 km\,s$15\,M_\odot$1 and supporting a compact object mass $15\,M_\odot$2 for canonical donor masses (Bowler, 2020, Picchi et al., 2020). The wind exhibits clear ionization stratification, with absorption velocities (and centroid offsets) systematically tracking the disk orientation.

3. Jet Launching, Precession, and Multiphase Outflows

SS 433 launches persistent twin jets at $15\,M_\odot$3 (deprojected velocity), each precessing about the disk axis with a remarkably stable period $15\,M_\odot$4 d and negligible secular period change $15\,M_\odot$5 (Robinson et al., 2017). Jet opening angles are $15\,M_\odot$6, matching the disk precession cone. High-resolution VLBI and VLBA imaging resolves the inner "core–jet" at milliarcsecond scales, revealing self-absorbed core emission, proper motions $15\,M_\odot$7–$15\,M_\odot$8 mas/d, and equipartition brightness temperatures $15\,M_\odot$9–$i\simeq 78\degree$0 K (Marshall et al., 2013, Paragi et al., 2012).

X-ray and optical spectroscopy show that X-ray and optical jet emission regions are separated by $i\simeq 78\degree$1 cm (comparable to $i\simeq 78\degree$2 light-hours), with X-ray emission fading in thousands of seconds and optical components persisting for days. Jet Doppler-shift variations imply aperiodic direction perturbations, consistent with weak shocks due to jet–environment interaction (Marshall et al., 2013). The jet base density is $i\simeq 78\degree$3–$i\simeq 78\degree$4 cm$i\simeq 78\degree$5, with Ni abundances enhanced by $i\simeq 78\degree$6 over solar, indicative of an anomalous core-collapse SN nucleosynthetic origin (Marshall et al., 2013).

4. High-Energy Emission and Jet-Driven Particle Acceleration

Parsec-scale jets from SS 433 drive particle acceleration at the interface with the W50 SNR, resolved by ground-based TeV γ-ray instruments (HAWC, VERITAS, H.E.S.S.) (Collaboration et al., 2018, Rho et al., 2019, Aharonian et al., 2024, Kleiner, 25 Sep 2025, Collaboration et al., 23 Mar 2026). Both lobes exhibit extended TeV emission, spatially coincident with the jet–SNR interaction regions ("e1"/"w1"), with angular extents $i\simeq 78\degree$7°–$i\simeq 78\degree$8° corresponding to $i\simeq 78\degree$9–$79\degree$0 pc. Multi-TeV γ-ray spectra are power laws ($79\degree$1–$79\degree$2), with steady flux at $79\degree$3 cm$79\degree$4 s$79\degree$5 per lobe (Collaboration et al., 23 Mar 2026, Kleiner, 25 Sep 2025).

Energy-dependent centroid shifts, with higher-energy emission closer to the binary ($79\degree$6–$79\degree$7 pc), require an acceleration site at the jet collimation shocks, not at the jet terminal e3/w2 knots. Particle-transport modeling yields downstream advection speeds $79\degree$8, diffusion coefficients $79\degree$9 cm$d\simeq 4.5$0 s$d\simeq 4.5$1, and electron energy cutoffs $d\simeq 4.5$2 TeV, attaining near-Bohm-limit acceleration efficiency (Aharonian et al., 2024). Magnetic field estimates from SED fitting are $d\simeq 4.5$3–$d\simeq 4.5$4G (Collaboration et al., 2018).

Spectral energy distributions from radio through TeV are well-reproduced by leptonic models requiring a single electron population (index $d\simeq 4.5$5, cutoff $d\simeq 4.5$6 PeV), total $d\simeq 4.5$7 energy $d\simeq 4.5$8 erg, and high acceleration efficiency; hadronic scenarios require unrealistic proton energetics or ambient densities (Collaboration et al., 2018, Rho et al., 2019). UHE ($d\simeq 4.5$9 TeV) emission at larger distances and disconnected hot spots may arise from $5.5$0PeV neutron beams and in-flight β-decay secondary $5.5$1/γ-photons produced in explosive jet–disk–photon interactions, with neutron collimation preserved over $5.5$2 light-years (Fargion et al., 7 Dec 2025).

5. High-Energy Variability, Gamma-ray Periodicity, and Outflow Structure

Fermi-LAT detects GeV emission offset by $5.5$3 pc from the SS 433 core. Only one excess exhibits the jet precession period $5.5$4 d, a "gamma-ray heartbeat" providing evidence for long-range particle transport or magnetic flux tube-guided propagation from the central engine to ambient clouds (Li et al., 2020). No phase-resolved VHE variability is seen in VERITAS data: orbital/precessional modulations are $5.5$5\% Crab flux (Collaboration et al., 23 Mar 2026). Constraints on variability and extension disfavor transient or highly anisotropic core emission at VHE.

Optical/IR emission and polarimetric studies confirm a stratified wind structure, with stationary lines (O I, He I, C II) mapping to distinct disk/wind/irradiated-stream zones by Doppler-lag and phase. Outflow through the $5.5$6 point is linked to uneclipsed, highly Doppler-split H$5.5$7 flares, with GRAVITY Br-$5.5$8 mapping confirming the presence of rotating, expanding equatorial shells at $5.5$9 (binary separation), angular momenta consistent with ejection from $a \approx 0.27$0 rather than accretion disk rim (Bowler, 2021).

6. Binary Evolution, Population Synthesis, and Progenitor Constraints

Probabilistic reconstruction of the zero-age main sequence (ZAMS) progenitor space using simulation-calibrated Bayesian inference (simulation-based calibration with COSMIC, dynamic nested sampling in 10D—masses, period, eccentricity, mass transfer, natal kick) finds: $a \approx 0.27$1, $a \approx 0.27$2, $a \approx 0.27$3 d, $a \approx 0.27$4, $a \approx 0.27$5, $a \approx 0.27$6, $a \approx 0.27$7 km s$a \approx 0.27$8 (Steinle et al., 23 Oct 2025). These ranges allow for lower-mass black-hole progenitors than previously inferred and accommodate high-eccentricity, common-envelope evolution. Population synthesis predicts that a Hertzsprung-gap donor at RLOF onset, $a \approx 0.27$9 yr after BH formation and rapid thermal-timescale mass transfer ($58\,R_\odot$0–$58\,R_\odot$1), matches the age constraints set by the W50 SNR (Han et al., 2020, Steinle et al., 23 Oct 2025).

7. Broader Context: Jet-Disk Symbiosis and Future Directions

SS 433 exemplifies jet-disk symbiosis under super-Eddington accretion. The X-ray/jet emission budget is dominated by reprocessing; observed X-rays ($58\,R_\odot$2 erg s$58\,R_\odot$3) are four orders of magnitude lower than the intrinsic disk energy release ($58\,R_\odot$4 erg s$58\,R_\odot$5) because a thick Compton/wind funnel thermalizes original X-rays, re-emitting them in the UV/EUV (Fabrika et al., 2010). The analogy to ULXs is direct: the observed reflected/soft-excess spectral signatures, hard/soft X-ray variability, and the flat, luminous incident spectrum support the notion that face-on ULXs are SS 433 analogs with observable funnel emission (Fabrika et al., 2010).

SKA, CTA, e-VLBI, and future space-based X-ray and spectropolarimetric campaigns promise unprecedented spatial, temporal, and spectral resolution on accretion flows, wind structure, jet launching, and cosmic-ray feedback from microquasars (Paragi et al., 2012). SS 433 remains key to understanding not only the physics of supercritical disks and baryonic jets but also the evolutionary pathways to and from high-mass X-ray binaries, GW progenitors, and PeVatron cosmic-ray sources.