High-mass Microquasars Research Articles

Effects of radiative losses on the relativistic jets of high-mass microquasars

Context.Relativistic jets are ubiquitous in astrophysics. High-mass microquasars (HMMQs) are useful laboratories for studying these jets because they are relatively close and evolve over observable timescales. The ambient medium into which the jet propagates, however, is far from homogeneous. Corresponding simulation studies to date consider various forms of a wind-shaped ambient medium, but typically neglect radiative cooling and relativistic effects.Aims.We investigate the dynamical and structural effects of radiative losses and system parameters on relativistic jets in HMMQs, from the jet launch to its propagation over several tens of orbital separations.Methods.We used 3D relativistic hydrodynamical simulations including parameterized radiative cooling derived from relativistic thermal plasma distribution to carry out parameter studies around two fiducial cases inspired by Cygnus X-1 and Cygnus X-3.Results.Radiative losses are found to be more relevant in Cygnus X-3 than Cygnus X-1. Varying jet power, jet temperature, or the wind of the donor star tends to have a larger impact at early times, when the jet forms and instabilities initially develop, than at later times when the jet has reached a turbulent state.Conclusions.Radiative losses may be dynamically and structurally relevant at least for Cygnus X-3 and thus should be examined in more detail.

Relativistic hydrodynamical simulations of the effects of the stellar wind and the orbit on high-mass microquasar jets

ABSTRACT High-mass microquasar jets, produced in an accreting compact object in orbit around a massive star, must cross a region filled with stellar wind. The combined effects of the wind and orbital motion can strongly affect the jet properties on binary scales and beyond. The study of such effects can shed light on how high-mass microquasar jets propagate and terminate in the interstellar medium. We study for the first time, using relativistic hydrodynamical simulations, the combined impact of the stellar wind and orbital motion on the properties of high-mass microquasar jets on binary scales and beyond. We have performed 3-dimensional relativistic hydrodynamic simulations, using the PLUTO code, of a microquasar scenario in which a strong weakly relativistic wind from a star interacts with a relativistic jet under the effect of the binary orbital motion. The parameters of the orbit are chosen such that the results can provide insight on the jet-wind interaction in compact systems like for instance Cyg X-1 or Cyg X-3. The wind and jet momentum rates are set to values that may be realistic for these sources and lead to moderate jet bending, which together with the close orbit and jet instabilities could trigger significant jet precession and disruption. For high-mass microquasars with orbit size a ∼ 0.1 AU, and (relativistic) jet power $L_j\sim 10^{37}(\dot{M}_w/10^{-6}\, {\rm M}_\odot \, {\rm yr}^{-1})$ erg s−1, where $\dot{M}_w$ is the stellar wind mass rate, the combined effects of the stellar wind and orbital motion can induce relativistic jet disruption on scales ∼1 AU.

HAWC Search for High-mass Microquasars

Abstract Microquasars with high-mass companion stars are promising very high energy (VHE; 0.1–100 TeV) gamma-ray emitters, but their behaviors above 10 TeV are poorly known. Using the High Altitude Water Cerenkov (HAWC) observatory, we search for excess gamma-ray emission coincident with the positions of known high-mass microquasars (HMMQs). No significant emission is observed for LS 5039, Cyg X-1, Cyg X-3, and SS 433 with 1523 days of HAWC data. We set the most stringent limit above 10 TeV obtained to date on each individual source. Under the assumption that HMMQs produce gamma rays via a common mechanism, we have performed source-stacking searches, considering two different scenarios: (I) gamma-ray luminosity is a fraction ϵ γ of the microquasar jet luminosity, and (II) VHE gamma rays are produced by relativistic electrons upscattering the radiation field of the companion star in a magnetic field B. We obtain ϵ γ < 5.4 × 10−6 for scenario I, which tightly constrains models that suggest observable high-energy neutrino emission by HMMQs. In the case of scenario II, the nondetection of VHE gamma rays yields a strong magnetic field, which challenges synchrotron radiation as the dominant mechanism of the microquasar emission between 10 keV and 10 MeV.

A model for high-mass microquasar jets under the influence of a strong stellar wind

Context. High-mass microquasars (HMMQs) are systems from which relativistic jets are launched. At the scales of several times the binary system size, the jets are expected to follow a helical path caused by the interaction with a strong stellar wind and orbital motion. Such a trajectory has its influence on the non-thermal emission of the jets, which also depends strongly on the observing angle due to Doppler boosting effects.Aims. We explore how the expected non-thermal emission of HMMQ jets at small scales is affected by the impact of the stellar wind and the orbital motion on the jet propagation.Methods. We studied the broadband non-thermal emission, from radio to gamma rays, produced in HMMQ jets up to a distance of several orbital separations, taking into account a realistic jet trajectory, different model parameters, and orbital modulation. The jet trajectory is computed by considering momentum transfer with the stellar wind. Electrons are injected at the position where a recollimation shock in the jets is expected due to the wind impact. Their distribution along the jet path is obtained assuming local acceleration at the recollimation shock, and cooling via adiabatic, synchrotron, and inverse Compton processes. The synchrotron and inverse Compton emission is calculated taking into account synchrotron self-absorption within the jet, free-free absorption with the stellar wind, and absorption by stellar photons via pair production.Results. The spectrum is totally dominated by the jet over the counter-jet due to Doppler boosting. Broadband emission from microwaves to gamma rays is predicted, with radio emission being totally absorbed. This emission is rather concentrated in the regions close to the binary system and features strong orbital modulation at high energies. Asymmetric light curves are obtained owing to the helical trajectory of the jets.Conclusions. The presence of helical shaped jets could be inferred from asymmetries in the light curves, which become noticeable only for large jet Lorentz factors and low magnetic fields. Model parameters could be constrained if accurate phase-resolved light curves from GeV to TeV energies were available. The predictions for the synchrotron and the inverse Compton radiation are quite sensitive of the parameters determining the wind-jet interaction structure.

Nonthermal emission from high-mass microquasar jets affected by orbital motion

Context. The stellar wind in high-mass microquasars should interact with the jet. This interaction, coupled with orbital motion, is expected to make the jet follow a helical, nonballistic trajectory. The jet energy dissipated by this interaction, through shocks for example, could lead to nonthermal activity on scales significantly larger than the system size.Aims. We calculate the broadband emission from a jet affected by the impact of the stellar wind and orbital motion in a high-mass microquasar.Methods. We employ a prescription for the helical trajectory of a jet in a system with a circular orbit. Subsequently, assuming electron acceleration at the onset of the helical jet region, we compute the spatial and energy distribution of these electrons, and their synchrotron and inverse Compton emission including gamma-ray absorption effects.Results. For typical source parameters, significant radio, X- and gamma-ray luminosities are predicted. The scales on which the emission is produced may reduce, but not erase, orbital variability of the inverse Compton emission. The wind and orbital effects on the radio emission morphology could be studied using very long baseline interferometric techniques.Conclusions. We predict significant broadband emission, modulated by orbital motion, from a helical jet in a high-mass microquasar. This emission may be hard to disentangle from radiation of the binary itself, although the light curve features, extended radio emission, and a moderate opacity to very high-energy gamma rays, could help to identify the contribution from an extended (helical) jet region.

Gamma rays from clumpy wind-jet interactions in high-mass microquasars

The stellar winds of the massive stars in high-mass microquasars are thought to be inhomogeneous. The interaction of these inhomogeneities, or clumps, with the jets of these objects may be a major factor in gamma-ray production. Our goal is to characterize a typical scenario of clump-jet interaction, and calculate the contribution of these interactions to the gamma-ray emission from these systems. We use axisymmetric, RHD simulations to model the emitting flow in a typical clump-jet interaction. Using the simulation results we perform a numerical calculation of the high-energy emission from one of these interactions. The radiative calculations are done for relativistic electrons locally accelerated at the jet shock, and the synchrotron and IC radiation spectra are computed for different stages of the shocked clump evolution. We also explore different parameter values, such as viewing angle and magnetic field strength. We generalize phenomenologically the results derived from one clump-jet interaction to multiple interactions under different wind models, estimating the clump-jet interaction rates, and the resulting luminosities in the GeV range. If particles are efficiently accelerated in clump-jet interactions, the apparent gamma-ray luminosity through IC scattering with the stellar photons can be significant even for rather strong magnetic fields and thus efficient synchrotron cooling. Moreover, despite the standing nature or slow motion of the jet shocks for most of the interaction stage, Doppler boosting in the postshock flow is relevant even for mildly relativistic jets. For clump-to-average wind density contrasts $\gtrsim 10$, clump-jet interactions could be bright enough to match the observed GeV luminosity in Cyg~X-1 and Cyg~X-3 when a jet is present in these sources, with required non-thermal-to-total available power fractions $\sim 0.01$ and $0.1$, respectively.

The effects of the stellar wind and orbital motion on the jets of high-mass microquasars

High-mass microquasar jets propagate under the effect of the wind from the companion star, and the orbital motion of the binary system. The stellar wind and the orbit may be dominant factors determining the jet properties beyond the binary scales.} An analytical study is performed to characterize the effects of the stellar wind and the orbital motion on the jet properties. Accounting for the wind thrust transferred to the jet, we derive analytical estimates to characterize the jet evolution under the impact of the stellar wind. We include the Coriolis force effect, induced by orbital motion and enhanced by the wind presence. Large-scale evolution of the jet is sketched accounting for wind-to-jet thrust transfer, total energy conservation, and wind-jet flow mixing. If the angle of the wind-induced jet bending is larger than its half-opening angle, the following is expected: (i) a strong recollimation shock; (ii) bending against orbital motion, caused by Coriolis forces and enhanced by the wind presence; and (iii) non-ballistic helical propagation further away. Even if disrupted, the jet can reaccelerate due to ambient pressure gradients, but wind entrainment can weaken this acceleration. On large scales, the opening angle of the helical structure is determined by the wind-jet thrust relation, and the wind-loaded jet flow can be rather slow. The impact of stellar winds on high-mass microquasar jets can yield non-ballistic helical jet trajectories, jet partial disruption and wind mixing, shocks, and possibly non-thermal emission. Among several observational diagnostics at different bands, the radio morphology on milliarcsecond scales can be particularly insightful regarding the wind-jet interaction.

ORIGIN OF MULTI-BAND EMISSION FROM THE MICROQUASAR CYGNUS X-1

We study the origin of non-thermal emissions from the Galactic black hole X-ray binary Cygnus X--1 which is a confirmed high mass microquasar. By analogy with those methods used in studies of active galactic nuclei, we propose a two-dimensional, time-dependent radiation model from the microquasar Cygnus X--1. In this model, the evolution equation for relativistic electrons in a conical jet are numerically solved by including escape, adiabatic and various radiative losses. The radiative processes involved are synchrotron emission, its self-Compton scattering, and inverse Compton scatterings of an accretion disk and its surrounding stellar companion. This model also includes an electromagnetic cascade process of an anisotropic $\gamma$-$\gamma$ interaction. We study the spectral properties of electron evolution and its emission spectral characteristic at different heights of the emission region located in the jet. We find that radio data from Cygnus X--1 are reproduced by the synchrotron emission, the \emph{Fermi} LAT measurements by the synchrotron emission and Comptonization of photons of the stellar companion, the TeV band emission fluxes by the Comptonization of the stellar photons. Our results show that: (1) Radio emission region extends from the binary system scales to the termination of the jet. (2) The GeV band emissions should originate from the distance close to the binary system scales. (3) The TeV band emissions could be inside the binary system, and these emissions could be probed by the upcoming CTA telescope. (4) The MeV tail emissions, which produces a strongly linearly polarized signal, are emitted inside the binary system. The location of the emissions is very close to the inner region of the jet.

3D simulations of microquasar jets in clumpy stellar winds

High-mass microquasars consist of a massive star and a compact object, the latter producing jets that will interact with the stellar wind. The evolution of the jets, and ultimately their radiative outcome, could depend strongly on the inhomogeneity of the wind, which calls for a detailed study. The hydrodynamics of the interaction between a jet and a clumpy wind is studied, focusing on the global wind and single clump-jet interplay. We have performed, using the code \textit{Ratpenat}, three-dimensional numerical simulations of a clumpy wind interacting with a mildly relativistic jet, and of individual clumps penetrating into a jet. For typical wind and jet velocities, filling factors of about > 0.1 are already enough for the wind to be considered as clumpy. An inhomogeneous wind makes the jet more unstable when crossing the system. Kinetic luminosities of the order 1.e37 erg/s allow the jet to reach the borders of a compact binary with an O star, as in the smooth wind case, although with a substantially higher degree of disruption. When able to enter into the jet, clumps are compressed and heated during a time of about their size divided by the sound speed in the shocked clump. Then, clumps quickly disrupt, mass-loading and slowing down the jet. We conclude that moderate wind clumpiness makes already a strong difference with the homogeneous wind case, enhancing jet disruption, mass-loading, bending, and likely energy dissipation in the form of emission. All this can have observational consequences at high-energies and also in the large scale radio jets.

TeV activity of Cygnus X-3 high mass binary system

Cygnus X-3 is a high mass X-ray binary and microquasar, with a compact object, which is either a neutron star or may be a black hole, and a companion object, which is a Wolf-Rayet star. The nature of the compact object is still uncertain. Cygnus X-3 galactic binary system has been regularly observed since a 1995 by SHALON Atmospheric Cherenkov telescope with the average gamma-ray flux (6.8 ± 0.7) × 10−13 cm−2 s−1. The observations of very high-energy gamma-radiation from the sources of this type would be important for understanding the nature of this astrophysical object.

Gamma-ray absorption and the origin of the gamma-ray flare in Cygnus X-1

The high-mass microquasar Cygnus X-1, the best-established candidate for a stellar-mass black hole in the Galaxy, has been detected in a flaring state at very high energies (VHE), E > 200 GeV, by the Atmospheric Cherenkov Telescope MAGIC. The flare occurred at orbital phase 0.91, where phase 1 is the configuration with the black hole behind the companion high-mass star, when the absorption of gamma-ray photons by photon-photon annihilation with the stellar field is expected to be highest. We aim to set up a model for the high-energy emission and absorption in Cyg X-1 that can explain the nature of the observed gamma-ray flare. We study the gamma-ray opacity due to pair creation along the whole orbit, and for different locations of the emitter. Then we consider a possible mechanism for the production of the VHE emission. We present detailed calculations of the gamma-ray opacity and infer from these calculations the distance from the black hole where the emitting region was located. We suggest that the flare was the result of a jet-clump interaction where the decay products of inelastic proton-proton collisions dominate the VHE outcome. We are able to reproduce the spectrum of Cyg X-1 during the observed flare under reasonable assumptions. The flare may be the first event of jet-cloud interaction ever detected at such high energies.

HIGH-ENERGY EMISSION AND ABSORPTION IN CYGNUS X-1

The high-mass microquasar Cygnus X-1 has been detected during a flaring state at very high energies, E > 200 GeV . The observation was performed by the Atmospheric Cherenkov Telescope MAGIC. It constitutes the first experimental evidence of very-high energy (VHE) emission produced by a Galactic stellar-mass black hole. The observed signal was detected in coincidence with an X–ray flare. The gamma-ray flare occurred when the compact object was located behind the companion star with respect to the line of sight (superior conjunction of the compact object). In this configuration the absorption of VHE photons by annihilation with the stellar photons is expected to be maximum. This suggests that the emission has been originated far above the compact object. The energy spectrum is well fitted by a relatively soft power-law. We present a model for the absorption and the emission of VHE gamma-rays in Cyg X-1. Detailed calculations of the gamma-ray opacity due to pair creation are provided and used to establish constraints on the emitting region. We propose that the high energy flare was the result of the interaction between the jet and a very dense clump from the stellar wind.

3D simulations of wind-jet interaction in massive X-ray binaries

High-mass microquasars may produce jets that will strongly interact with surrounding stellar winds on binary system spatial scales. We study the dynamics of the collision between a mildly relativistic hydrodynamical jet of supersonic nature and the wind of an OB star. We performed numerical 3D simulations of jets that cross the stellar wind with the code Ratpenat. The jet head generates a strong shock in the wind, and strong recollimation shocks occur due to the initial overpressure of the jet with its environment. These shocks can accelerate particles up to TeV energies and produce gamma-rays. The recollimation shock also strengthens jet asymmetric Kelvin-Helmholtz instabilities produced in the wind/jet contact discontinuity. This can lead to jet disruption even for jet powers of several times $10^{36}$ erg s$^{-1}$. High-mass microquasar jets likely suffer a strong recollimation shock that can be a site of particle acceleration up to very high energies, but also eventually lead to the disruption of the jet.

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.

STUDYING THE INTERACTION BETWEEN MICROQUASAR JETS AND THEIR ENVIRONMENTS

In high-mass microquasars (HMMQ), strong interactions between jets and stellar winds at binary system scales could occur. In order to explore this possibility, we have performed numerical two-dimensional hydrodynamical simulations of jets crossing the dense stellar material to study how the jet will be affected by these interactions. We find that the jet head generates strong shocks in the wind. These shocks reduce the jet advance speed, and compress and heat up the jet and wind material. In addition, strong recollimation shocks can occur where pressure balance between the jet side and the surrounding medium is reached. All this, together with jet bending, could lead to the destruction of jets with power < 1036 erg/s . The conditions around the outflow shocks would be convenient for accelerating particles up to ~ TeV energies. These accelerated particles could emit via synchrotron and inverse Compton (IC) scattering if they were leptons, and via hadronic processes if they were hadrons.

On the interaction of microquasar jets with stellar winds

Strong interactions between jets and stellar winds at binary system spatial scales could occur in high-mass microquasars. We study here, mainly from a dynamical but also a radiative point of view, the collision between a dense stellar wind and a mildly relativistic hydrodynamical jet of supersonic nature.}{We have performed numerical 2-dimensional simulations of jets, with cylindrical and planar (slab) symmetry, crossing the stellar wind material. From the results of the simulations, we derive estimates of the particle acceleration efficiency, using first order Fermi acceleration theory, and give some insight on the possible radiative outcomes. We find that, during jet launching, the jet head generates a strong shock in the wind. During and after this process, strong recollimation shocks can occur due to the initial overpressure of the jet with its environment. The conditions in all these shocks are convenient to accelerate particles up to $\sim$ TeV energies, which can lead to leptonic (synchrotron and inverse Compton) and hadronic (proton-proton) radiation. In principle, the cylindrical jet simulations show that the jet is stable, and can escape from the system even for relatively low power. However, when accounting for the wind ram pressure, the jet can be bent and disrupted for power $\la 10^{36}$ erg s$^{-1}$.

A microquasar model applied to unidentified gamma-ray sources

the date of receipt and acceptance should be inserted later Abstract. Among unidentified gamma-ray sources in the galactic plane, there are some that present significant variability and have been proposed to be high-mass microquasars. To deepen the study of the possible associa- tion between variable low galactic latitude gamma-ray sources and microquasars, we have applied a leptonic jet model based on the microquasar scenario that reproduces the gamma-ray spectrum of three unidentified gamma- ray sources, 3EG J1735 1500, 3EG J1828+0142 and GRO J1411 64, and is consistent with the observational constraints at lower energies. We conclude that if these sources were generated by microquasars, the particle acceleration processes could not be as efficient as in other objects of this type that present harder gamma-ray spectra. Moreover, the dominant mechanism of high-energy emission should be synchrotron self-Compton (SSC) scattering, and the radio jets may only be observed at low frequencies. For each particular case, further predictions of jet physical conditions and variability generation mechanisms have been made in the context of the model. Although there might be other candidates able to explain the emission coming from these sources, microquasars cannot be excluded as counterparts. Observations performed by the next generation of gamma-ray instruments, like GLAST, are required to test the proposed model.

Microquasar Models for 3EG J1828+0142 and 3EG J1735-1500

Microquasars are promising candidates to emit high-energy gamma-rays. Moreover, statistical studies show that variable EGRET sources at low galactic latitudes could be associated with the inner spiral arms. The variable nature and the location in the Galaxy of the high-mass microquasars, concentrated in the galactic plane and within 55 degrees from the galactic center, give to these objects the status of likely counterparts of the variable low-latitude EGRET sources. We consider in this work the two most variable EGRET sources at low-latitudes: 3EG J1828+0142 and 3EG J1735-1500, proposing a microquasar model to explain the EGRET data in consistency with the observations at lower energies (from radio frequencies to soft gamma-rays) within the EGRET error box.

On the Nature of the Variable Gamma-Ray Sources at Low Galactic Latitudes

Population studies of EGRET gamma-ray sources indicate that there is a distinctive population of bright sources at low galactic latitudes. The sources have a distribution consistent with that of young galactic objects, with a concentration toward the inner spiral arms. There is a subgroup that displays strong variability with timescales from days to months. Following an earlier suggestion by Kaufman Bernado et al. (2002), we explore the possibility that these sources could be high-mass microquasars. Detailed models for the gamma-ray emission that include inverse Compton interactions of electrons in the relativistic jets and photons from all local fields (stellar UV photons, synchrotron photons, soft X-ray photons from the accretion disk, and hard X-ray photons from a corona) are presented. We conclude that microquasars are excellent candidates for the parent population of the subgroup of variable low-latitude EGRET sources.

High-mass microquasars and low-latitude gamma-ray sources

Population studies of unidentified EGRET sources suggest that there exist at least three different populations of galactic gamma-ray sources. One of these populations is formed by young objects distributed along the galactic plane with a strong concentration toward the inner spiral arms of the Galaxy. Variability, spectral and correlation analysis indicate that this population is not homogeneous. In particular, there is a subgroup of sources that display clear variability in their gamma-ray fluxes on timescales from days to months. Following the proposal by Kaufman Bernadó et al. (2002), we suggest that this group of sources might be high-mass microquasars, i.e. accreting black holes or neutron stars with relativistic jets and early-type stellar companions. We present detailed inhomogeneous models for the gamma-ray emission of these systems that include both external and synchrotron self-Compton interactions. We have included effects of interactions between the jet and all external photon fields to which it is exposed: companion star, accretion disk, and hot corona. We make broadband calculations to predict the spectral energy distribution of the emission produced in the inner jet of these objects up to GeV energies. The results and predictions can be tested by present and future gamma-ray instruments like INTEGRAL, AGILE, and GLAST.