High-energy emission from jet-clump interactions in microquasars
High-mass microquasars are binary systems consisting of a massive star and an accreting compact object from which relativistic jets are launched. There is considerable observational evidence that winds of massive stars are clumpy. Individual clumps may interact with the jets in high-mass microquasars to produce outbursts of high-energy emission. Gamma-ray flares have been detected in some high-mass X-ray binaries, such as Cygnus X-1, and probably in LS 5039 and LS I+61 303. We predict the high-energy emission produced by the interaction between a jet and a clump of the stellar wind in a high-mass microquasar. Assuming a hydrodynamic scenario for the jet-clump interaction, we calculate the spectral energy distributions produced by the dominant non-thermal processes: relativistic bremsstrahlung, synchrotron and inverse Compton radiation, for leptons, and for hadrons, proton-proton collisions. Significant levels of emission in X-rays (synchrotron), high-energy gamma rays (inverse Compton), and very high-energy gamma rays (from the decay of neutral pions) are predicted, with luminosities in the different domains in the range ~ 10^{32}-10^{35} erg/s. The spectral energy distributions vary strongly depending on the specific conditions. Jet-clump interactions may be detectable at high and very high energies, and provide an explanation for the fast TeV variability found in some high-mass X-ray binary systems. Our model can help to infer information about the properties of jets and clumpy winds by means of high-sensitivity gamma-ray astronomy.
💡 Research Summary
The paper investigates the high‑energy radiation that can be produced when a dense clump of a massive‑star wind collides with the relativistic jet of a high‑mass microquasar (MQ). The authors start by noting that massive‑star winds are known to be highly inhomogeneous, consisting of dense clumps embedded in a more tenuous outflow. When such a clump penetrates the jet, a strong shock forms at the interface, generating turbulence and a shear layer that are conducive to particle acceleration. The acceleration mechanism is assumed to be diffusive shock acceleration (DSA), which yields power‑law energy distributions (∝ E⁻²) for both electrons and protons up to a maximum energy limited by radiative cooling and the size of the interaction region.
The paper models the jet as a mildly relativistic flow (≈0.3 c) with a radius of order 10⁹ cm and an internal magnetic field of 10³–10⁴ G. The clumps are characterized by radii between 10⁹ and 10¹¹ cm and densities ranging from 10⁻¹³ to 10⁻¹¹ g cm⁻³, corresponding to mass fractions of 10⁻⁴–10⁻² of the total stellar wind. When the clump enters the jet, a bow shock propagates into the jet material while a reverse shock travels into the clump. The shocked region is assumed to be homogeneous for the purpose of calculating the non‑thermal emission.
Electrons accelerated at the shock lose energy through three dominant processes: (i) synchrotron radiation in the jet magnetic field, producing X‑ray photons (keV energies); (ii) external inverse‑Compton (EIC) scattering of the intense UV photon field from the massive companion star, generating high‑energy γ‑rays (GeV); and (iii) relativistic bremsstrahlung in the dense clump material, contributing to the MeV band. Protons, on the other hand, interact with the clump’s dense gas via proton‑proton (pp) collisions, leading to the production of neutral pions that decay into very‑high‑energy (VHE) γ‑rays (TeV).
The authors perform a series of numerical calculations to obtain spectral energy distributions (SEDs) for a range of parameter sets. In most cases the SED exhibits a double‑humped structure: a low‑energy hump from synchrotron emission and a high‑energy hump that is a combination of EIC and π⁰‑decay components. The relative strength of these components depends sensitively on clump size, density, and the magnetic field strength. Small, very dense clumps favor pp interactions, yielding TeV luminosities up to ∼10³⁵ erg s⁻¹, whereas larger, less dense clumps lead to dominant electron processes and stronger X‑ray synchrotron output. The predicted luminosities span 10³²–10³⁵ erg s⁻¹ across the X‑ray, GeV, and TeV bands, making them potentially detectable with current and upcoming facilities.
Temporal evolution is also addressed. The interaction time is set by the clump crossing time, typically seconds to minutes, which produces rapid γ‑ray flares. The model naturally explains the fast TeV variability observed in systems such as Cygnus X‑1, LS 5039, and LS I +61 303. By fitting observed flare light curves and spectra, one can infer jet parameters (e.g., magnetic field, particle acceleration efficiency) and clump properties (size distribution, density contrast).
The discussion highlights observational prospects. The predicted fluxes are within the sensitivity of instruments like Fermi‑LAT (GeV) and ground‑based Cherenkov arrays (H.E.S.S., MAGIC, and the future CTA) for TeV energies. The authors argue that systematic monitoring of microquasars could reveal a population of short‑duration γ‑ray flares that are signatures of jet‑clump collisions.
Finally, the paper acknowledges limitations. The hydrodynamic treatment neglects magnetic field amplification and anisotropic particle transport, which could modify acceleration efficiencies and cooling rates. The clump formation mechanisms and their statistical properties remain uncertain, and a full three‑dimensional magnetohydrodynamic (MHD) simulation would be required for a more realistic description. Nonetheless, the work provides a solid framework linking wind clumpiness to high‑energy variability in microquasars and offers a pathway to use γ‑ray observations as diagnostics of both jet physics and stellar wind structure.