Nonthermal emission from high-mass microquasar jets affected by orbital motion
Key Words.:X-rays: binaries - radiation mechanisms: nonthermal - stars: winds, outflows - stars: massive
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.
High-mass microquasars (HMMQ) are X-ray binaries that host a massive star and a compact object (CO) from which jets are produced. The stellar wind can strongly influence the jet propagation, both because the jet has to propagate surrounded by wind material, and because the wind lateral impact may significantly bend the jet away from the star. Several authors have used numerical and analytical methods to study, on the scales of the binary, the interaction of HMMQ jets with stellar winds and their radiative consequences (e.g., Romero et al. 2003; Romero & Orellana 2005; Perucho & Bosch-Ramon 2008; Owocki et al. 2009; Araudo et al. 2009; Perucho et al. 2010; Perucho & Bosch-Ramon 2012; Yoon & Heinz 2015; Bosch-Ramon & Barkov 2016; Yoon et al. 2016).
On scales larger than the binary system, orbital motion should also affect the dynamics of the jet, making it follow a helical trajectory (e.g., Bosch-Ramon 2013; Bosch-Ramon & Barkov 2016). There are a few well-established microquasars in which evidence of a helical jet has been observed; for example, SS 433, 1E 1740.72942 and Cygnus X-3 (Abell & Margon 1979; Mioduszewski et al. 2001; Stirling et al. 2002; Miller-Jones et al. 2004; Luque-Escamilla et al. 2015; see also Sell et al. 2010 for the possible case of Circinus X-1). In the case of SS 433, the jet helical geometry likely originates in the accretion disk (e.g., Begelman et al. 2006); in the case of 1E 1740.72942, the system likely hosts a low-mass star and the role of its wind could be minor; in the case of Cygnus X-3, the apparent jet helical shape has an unclear origin. Bosch-Ramon & Barkov (2016) proposed that Cygnus X-1 and Cygnus X-3 could be affected by orbital motion and the stellar wind, although the uncertainties on the jet and wind properties make quantitative predictions difficult. Unlike the case of Cygnus X-3, no clear evidence of a helical or bent jet has yet been found for Cygnus X-1 (Stirling et al. 2001).
In this work we study the implications of the impact of the stellar wind and orbital motion for the nonthermal emission of a HMMQ, from radio to gamma rays. Using an analytical prescription for the jet dynamics based on the results of Bosch-Ramon & Barkov (2016), the nonthermal radiation from a helical jet (spectra, light curves, and morphology) is computed numerically, considering a leptonic model with synchrotron and inverse Compton (IC) emission.
The paper is organized as follows: In Sect. 2, the interaction between the stellar wind and the jet is described. The technical aspects of the model are explained in Sect. 3. The results of the calculations are presented in Sect. 4 and a discussion is provided in Sect. 5. Unless stated otherwise, the convention is used throughout the paper, with in cgs units.
2 Jet-wind interaction
The physical scenario considered in this work consists of a massive star with a strong stellar wind and an accreting CO from which two jets are launched in opposite directions perpendicular to the orbital plane. The binary has a relatively compact orbit, taken circular for simplicity, with a period of days and a separation of cm (). The luminosity of the star is taken to be erg s, and its temperature K (). The stellar wind is assumed to be spherically symmetric, with a mass-loss rate and a velocity of M yr and cm s, respectively (typical for O-type stars; e.g., Muijres et al. 2012). Such a system would have similar parameters to those of Cygnus X-1, although it would be significantly wider and with a weaker wind than Cygnus X-3 (but the adopted wind and jet momentum rate relation will be similar; see below).
The orbit is assumed to lie on the -plane, with the star at the center of coordinates. The CO would be on the -axis, orbiting the star counter clockwise, and the jet (counter jet) would initially point along (). The axes are defined in a frame corotating with the CO around the companion star. The jets are assumed to have a conical shape, with an initial half-opening angle rad . In this work, we do not need to adopt a specific value for , although to make our predictions qualitatively valid, it should be smaller than the jet deflection angle imposed by the wind impact (see below). We take an initial jet Lorentz factor of , although the final results are not sensitive to this parameter as long as the jet is strongly supersonic and mildly relativistic. The jet power is taken as erg s.
Figure 1 sketches the scenario under study. The jet goes through three stages, the first of which is not present in the figure: Initially, the jet moves perpendicularly to the orbital plane. Then, the wind impact produces an asymmetric recollimation shock on the jet, which gets inclined away from the star (stage 2). Further from the binary, the combined effect of the stellar wind and orbital motion leads to a (orbit-related) force that bends the already inclined jet, now against orbital motion (clockwise in Fig. 1). In this way, the jet becomes helical, with asymmetric shocks being expected because of the orbit-related force (stage 3). Here it is assumed that electrons are efficiently accelerated in those shocks, ascribed here to one specific point: the onset of the helical jet, which neglects the fact that acceleration may occur along the helical jet region. Relativistic protons may also be generated in those shocks, but their radiation efficiency would be much lower in general, and therefore they are not considered in our calculations (see e.g., Bosch-Ramon & Khangulyan 2009, for a discussion of the different electron and proton cooling timescales in HMMQ).
For the radiation calculations, a fraction of is injected in the form of nonthermal electrons right where the jet starts its helical path. This fraction is only constrained by , although the predicted luminosities are proportional to its value and easy to translate if a different fraction is adopted. Nevertheless, adopting does not require particularly efficient acceleration processes. The inclination of the system with respect to the line of sight of the observer, assumed to be at a reference distance of kpc, is characterized by the angle . The values for the different fixed parameters of our model are listed in Table 1.
|Star luminosity||erg s|
|Wind speed||cm s|
|Wind mass-loss rate||M yr|
|Jet luminosity||erg s|
|Initial jet Lorentz factor||2|
|Half-opening angle||0.1 rad|
|Nonthermal energy fraction|
|Effective flow velocity||cm s|
|Radial flow velocity||cm s|
|Distance to the observer||3 kpc|
|Radio interferometer resolution||FWHM||1 mas|
2.1 Wind impact: jet bending
Interaction between the stellar wind and the jet can lead to bending of the latter in the -axis direction (stage 2 above). This bending can be characterized by an angle with respect to the jet initial direction (i.e., the -axis; see Fig. 1). Following Bosch-Ramon & Barkov (2016) we introduce a nondimensional parameter, , which corresponds to the ratio between the wind momentum rate intercepted by the jet, and the jet momentum rate:
where is in units of M yr in the rightmost expression, is the stellar wind momentum rate, and is the jet momentum rate, which can be expressed in terms of , , and the initial jet velocity in units :
An approximate computation of the bending angle , with a discrepancy of less than 10% with respect to the numerical case, yields the following expression (for more details, see the Appendix in Bosch-Ramon & Barkov 2016):
For bending angles , the lateral impact of the stellar wind may produce a weak recollimation shock or a sound wave in the jet, and the large-scale jet evolution should not be strongly affected (although the role of instability growth may still be important; see, e.g., Perucho & Bosch-Ramon 2008; Perucho et al. 2010; Perucho & Bosch-Ramon 2012, under uniform and clumpy stellar winds). However, when is significantly larger than , a strong asymmetric recollimation shock is expected, with the jet substantially deviating from its initial direction. In such a case, the jet is likely to turn into a helical-shaped structure at larger scales (stage 3, above). Consequently, bending due to stellar wind impact becomes important for the jet evolution when .
2.2 Orbital effects: helical jet
The orbital motion of the system combined with a jet inclined by the stellar wind (stage 2) leads, through the orbit-related force, to a jet helical pattern. This effect becomes clearer for significant bending angles. Otherwise, for , the conical jet expansion will dominate the geometry of the resulting structure. We note that a helical pattern would be expected from an inclined jet even under ballistic jet propagation, but the presence of the stellar wind makes jet propagation nonballistic. Instability growth should also play a role, but here we consider its effects only phenomenologically through a slower jet flow, caused by the wind mass entrainment and hence deceleration associated to instabilities.
One can estimate the distance that a nonballistic jet travels in the -axis direction before acquiring a negative velocity -component, i.e., bending against the orbital motion. With the additional assumptions that the half-opening angle of the jet remains constant after being bent (see below), and that is much larger than the orbital separation, we get (see Bosch-Ramon & Barkov 2016, for the derivation of this expression):
where is the orbital angular velocity, with the normalization corresponding to a period of several days (e.g., rad s), and is a constant of , introduced to account for the specific details of the wind-jet interaction. Here we fix .
Knowing and , we can compute the height at which the jet starts forming the helical structure: . For the parameters given above, and following Eqs. (3) and (4), the location vector of the beginning of the helical structure is cm, with . We note that cm, as required. The distance from to the star is cm and to the CO is cm. The jet curvature due to orbital motion for distances smaller than is neglected, that is, .
The point (the -helical- jet base hereafter) is the position where nonthermal electrons are injected. The evolution of these electrons is followed from this point onwards along the helical jet assuming that no more significant electron acceleration happens further downstream in the jet111We note that as the jet propagation is nonballistic, (weak) acceleration regions may also exist all along the jet due to jet kinetic energy dissipation.. Electron acceleration may also happen closer to the binary system (see Sect. 5), but here we focus specifically on the nonthermal emission produced in the helical jet region.
From onwards along the jet, the stellar wind should have a very important dynamical role. As the jet becomes helical due to the presence of wind plus orbital motion, the former is also pushing the latter in the radial direction. Therefore, the radial velocity () of the jet flow may be as slow as the stellar wind itself, or higher if the shocked stellar wind is effectively accelerated radially by the jet. In addition, effects of wind entrainment are likely to affect the jet flow itself, slowing the latter down to an effective velocity . We adopt here and as phenomenological parameters that fulfill . They provide a simplified prescription for the helical jet geometry, defining the pitch of the helical jet trajectory. It is also assumed that this trajectory is confined to the conical surface, characterized, in the nonrotating frame, by the direction determined by the vector along the orbit in that frame (i.e., a ring traced by the points , plus the radial direction from the star), which should be a relatively good approximation as long as .
The half-opening angle of the helical jet flow () need not be constant, as the nonballistic nature of the jet trajectory, with weak-shock heating, rarefaction waves, wind (re)confinement, and mass-loading, may lead to (re)collimation or (re)widening of the helical jet, depending on the details of the jet and wind interaction. For simplicity, at this stage we assume . In addition, we also take cm s and cm s; this is somewhat arbitrary, but except for the radio morphology (see Sect. 5) their actual values do not qualitatively change the results. Hydrodynamical simulations are required to properly characterize , , , the role of instability growth, mass-loading, and so on. Nevertheless, we consider our dynamical model realistic on scales of a few , as large-scale perturbations should grow up to a size to disrupt the helical structure.
3 Model description
We use a semi-analytical code to compute the energy and spatial distribution of the accelerated electrons. As the orbital period is significantly longer than the crossing time of the helical jet length (), the nonthermal electrons are assumed to be in the steady state at each orbital phase.
3.1 Energy losses
Analytical expressions for the energy losses are used in this work. For synchrotron cooling, assuming, for simplicity, an isotropic magnetic field in the flow frame, one obtains, in cgs units (Longair 1981):
where is the Thomson cross-section, , and is the Lorentz factor of an electron with energy . The magnetic field at the jet base is parametrized through the ratio of magnetic pressure to stellar photon energy density, , as
This way of normalizing allows an easy comparison between the expected luminosity of synchrotron and IC emission; for IC in the Thomson regime, the two processes contribute similarly for electrons of the same energy if . However, for electrons of energy (with being the Boltzmann constant), which scatter photons in the IC Klein-Nishina regime, already leads to global strong synchrotron dominance, which is corroborated by our calculations. One can relate with the ratio of magnetic pressure to jet energy density () for the same via: ( for the values adopted here).
Losses due to IC are computed as in Khangulyan et al. (2014) for an isotropic electron population embedded in the stellar photon field taken as a black body of temperature . This expression is valid in both the Thomson and the Klein-Nishina regimes and has an accuracy of :
with being the radius of the star, the distance from the star to the emitter, , and
For a conical jet, adiabatic losses are
where is the jet expansion velocity and is the helical jet radius.
The total energy-loss rate is
3.2 Electron distribution
The nonthermal electrons are injected in the jet base with an energy distribution that follows a power law of index -2, typical for efficient nonthermal sources (and for acceleration in nonrelativistic strong shocks; Drury 1983), plus a cutoff at high energies set by energy and escape losses:
We note that a much steeper would mean inefficient electron acceleration, whereas a much harder one would imply exceptionally efficient acceleration in which most of energy is in the highest-energy electrons.
The normalization of the electron injection is taken as a fraction of the jet power available for the acceleration of nonthermal electrons, :
The value of is fixed here to 1 MeV for simplicity, although a very high , greater than GeV for example, would affect the more compact radio emission. The value of is obtained from electron acceleration constraints.
Electrons are accelerated gaining energy at a rate , where is the acceleration efficiency, and the electron charge. We take , which would correspond, for example, to acceleration by a high-velocity (nonrelativistic) shock under Bohm diffusion (Drury 1983). Low acceleration efficiencies would reduce the maximum energy of synchrotron photons, for which the photon energy is , much more than that of IC photons, produced in the Klein-Nishina regime at the highest spectral end, for which (for IC in Thomson: ).
Equating yields the maximum energy that electrons can reach in the accelerator (at ): . However, electrons can diffuse away from the accelerator before they reach . Therefore, the diffusion timescale, , is also to be compared with the acceleration timescale, , where is the characteristic jet base radius, and is the diffusion coefficient assuming it proceeds in the Bohm regime. This gives a maximum energy value ; for a larger diffusion coefficient (still ), . The value is obtained as the smallest among and , whereas can be taken as several times .
To compute the electron energy distribution along the jet, the latter is divided in cylindrical segments with increasing radius. The segment lengths are determined (with some exception; see below) by multiplied by one fifth of the local shortest cooling time, which is derived from . In this way, the energy and spatial evolution of the fastest evolving electrons is reasonably well sampled. The magnetic field is assumed to be mostly perpendicular to the flow motion, typical for jets far from their origin, and therefore under frozen conditions and constant . Each individual jet division is treated as a homogeneous emitter, which in the worse case is correct within a % error (spatial scales grow segment by segment by a factor , where is the number of segments).
For , synchrotron cooling times become very short (, with being the advection timescale), and the large makes calculations very long. In those extreme cases, we have adopted an approximation which consists in limiting and then introducing a correction factor, , to the injected population at the first segment. This simplified approach can overestimate the high energy part of the electron distribution by a factor around two, but speeds up calculations by several orders of magnitude for extreme cases. Nevertheless, we indicate that is too high, as in the region of interest the jet power is likely largely dominated by kinetic energy, and is considered here for illustrative purposes only.
To find the energy evolution up to a given segment, one computes the initial energy that electrons of energy had when they were injected at the jet base. This is done iteratively by making small time steps back in time (and in space) to sample the energy evolution of electrons backwards from a given segment:
Here, is the energy that an electron had at a time , and is the electron energy at time . We note that due to the negative sign of . The linear approximation used in Eq. (13) is only valid as long as .
Once is known, the electron energy distribution at the jet base, , of a segment k with length , is computed as
with (if small enough; otherwise ). The electron energy distribution at a segment k is computed from the distribution when at the segment k-1:
where and are energies related by Eq. (13) in segments k and k-1, respectively222We note that Eq. (15) should be corrected for segment length and velocity if they are not constant.. By applying this procedure down to the jet base, one gets the general expression for the electron energy distribution at each jet segment:
keeping in mind that and .
3.3 Emission of radiation
The power per photon energy unit radiated by an electron of energy via synchrotron emission is computed following (Pacholczyk 1970)
where is the pitch angle between the electron velocity and vectors, and , with being the characteristic synchrotron photon energy defined as
The geometry of is not well known, so an isotropic distribution in the fluid frame is assumed, taking for simplicity . The synchrotron spectral energy distribution (SED) for an isotropic population of electrons with distribution is:
The IC emission is computed considering that relativistic electrons interact with a beam of photons with energy . The IC kernel for an angle with respect to the initial beam direction is (Aharonian & Atoyan 1981):
where , , and can vary in the limits . The (apparent) IC SED at a given direction is obtained by convolving Eq. (21) with the electron energy distribution and the energy distribution density of stellar photons :
Target photons other than those of stellar origin, like those coming from an accretion disk, are unimportant for the jet regions studied in this work. Synchrotron self-Compton is also negligible even when one assumes that most of the nonthermal energy is released as synchrotron radiation.
3.4 Absorption of radiation
The main absorption mechanism considered in this work is electron-positron pair production from gamma rays interacting with stellar photons. The cross section for this process is (Gould & Schréder 1967)
and are the energies of the high-energy and stellar photons, respectively, and is the angle between their propagation directions.
The optical depth of gamma-ray absorption is (Gould & Schréder 1967)
where is the distance to the observer covered by the gamma ray, and is the minimum energy of the gamma rays absorbed by the target stellar photons.
For the cases studied here, synchrotron self-absorption is estimated in radio by computing the optically thick-to-optically thin transition frequency (e.g., Bosch-Ramon 2009), and is found to be negligible in general; it has therefore not been included in the calculations. Free-free radio absorption in the stellar wind is also negligible at the scales of the emitting region. Nonthermal ultraviolet (UV) photons are also strongly absorbed on their way to the observer, but we have not included this effect as UV photons from the star would be largely dominant in any case. In a very compact system with a very powerful stellar wind, such as, for example, Cygnus X-3, radio absorption could still be an issue at a distance from the binary, although we defer specific source studies for future work.
The results are obtained using the parameter values given in Table 1, together with the sets of values presented in Table 2 for , the orbital phase , and , whose variations have strong effects on the nonthermal emission. Considering the free parameters, the electron energy distribution is affected only by , whereas SEDs, light curves, and the emitter geometry are also affected by and . Most SEDs are calculated for as a typical example, which corresponds to one of the orbit nodes. Varying , we also compute some SEDs and IC light curves to illustrate orbital variations of this component. A comparison of SEDs and light curves is done between the extended and the point-like (one-zone) emitter, to check which features arise when an extended (helical) jet is considered. Finally, maps at 5 GHz are computed to illustrate the expected radio morphologies from a helical jet.
4.1 Nonthermal electrons
The electron energy distribution for the jet333The jet and counter-jet electron energy distributions are equal in our model., to emphasize the energy content at different energy scales, is shown in Fig. 2 for and 1, and up to three different jet lengths starting from : , and (whole jet). As seen in the figure, a significant fraction of electrons below GeV survive beyond , and for the lowest , TeV electrons can also reach a larger distance from the binary (see also, e.g., Khangulyan et al. 2008). As expected, the energy distribution is steep for electron energies dominated by synchrotron and Thomson IC losses. When IC dominates, the SED becomes harder at GeV because of the Thomson-to-Klein-Nishina transition of the electron distribution, getting steeper again only at the highest energies due to synchrotron losses (and the effect of ).
4.2 Spectral energy distribution
Figure 3 shows the SEDs of the synchrotron and the IC emission for the same cases studied in Fig. 2, showing now within each panel the jet and the counter-jet contributions separately. An inclination of is adopted, and . As expected from the severe drop in radiation efficiencies, the contribution to the emission of the regions beyond is minor for all the cases studied, although radio emission is still non-negligible beyond that distance, as seen in Sect. 4.5. It is worth noting that most of the IC radiation comes from the counter-jet because of the larger IC -values. Doppler boosting is negligible for the adopted -value and is not included.
The effect of the magnetic field on the synchrotron and IC SEDs is shown in Fig. 4 for and , and for the jet and the counter-jet separately. As expected, the synchrotron to IC luminosity ratio depends strongly on . Figure 4 also shows the effect of gamma-ray absorption (stronger for than for ), which is quite modest given the relatively large distances from the star and the -value considered in this plot.
Figure 5 shows the effect of the orbital phase on the IC spectrum. With and , the SEDs for and are plotted. For , which corresponds to the CO superior conjunction (i.e., the CO is behind the star), there is a drop of up to five orders of magnitude in the counter-jet emission above 10 GeV due to gamma-ray absorption, which makes the jet dominate the radiation output in this range, effectively smoothing the overall impact of gamma-ray absorption for this orbital phase.
4.3 Light curve
Figure 6 shows the IC light curves for . The peak around and GeV for is explained by the larger IC -angles, and the small dip for manifests some level of gamma-ray absorption, which starts at GeV in a relatively hard IC SED (see Fig. 5). For energies GeV, absorption has a very strong impact on the fluxes by a factor of several, although the jet component smoothens somewhat the light curve around (superior conjunction of the CO) because of reduced gamma-ray absorption due to a smaller -value (see Fig. 5). On the other hand, the counter-jet still leads to significant IC fluxes around (inferior conjunction; see Fig. 5), unlike an emitter located very close to the CO, which would lead to stronger differences in flux.
4.4 Extended versus one-zone emitter
Figure 7 shows a comparison between the SEDs of an extended (helical) jet and a one-zone emitter, for the jet and the counter-jet, for , , and . The one-zone emitter is assumed to be just one segment located at and of length , which sets the electron escape time. Because of fast radiation losses, the one-zone approximation works well in the high-energy part of the synchrotron and IC SEDs, as expected. Qualitatively, the extended and the one-zone emitters look quite similar, although differences of up to a factor of 2 are seen, mainly but not only in the jet synchrotron and IC components at low energies, as low-energy electrons escaped from the one-zone region are still radiatively relevant. We note that inaccurate source knowledge and model simplifications imply systematic uncertainties likely larger than a factor of 2, which means that the SED alone cannot help to discriminate a helical jet model from that of a one-zone emitter.
Figure 8 compares the light curves of an extended and a one-zone emitter in the cases presented in Fig. 6. The fluxes are larger (lower) for the extended emitter for GeV ( GeV), and are symmetric around . The different behavior between the two energy bands is caused by the fact that electrons emitting GeV photons via IC reach farther from the star than those emitting GeV photons (see Sect. 4.1), both for the counter-jet and the jet, decreasing the IC flux as they meet a more diluted target photon field. We note that for higher jet powers, the helical jet would start further away (see Eqs. 1 and 4), in which case the light curve would become asymmetric for the extended jet.
4.5 Radio emission
Radio sky maps at 5 GHz for and different orbital phases are shown in Figs. 9 and 10 for = 30 and = 60, respectively. These maps are obtained following three steps: first, the helical jet structure is projected into the plane of the sky. Subsequently, the emission from each segment is convolved with a Gaussian with standard deviation to approximately distribute the emission over the actual sky-projected area of each segment. Finally, this emission is convolved again with a Gaussian with FWHM = 1 mas in order to mimic the response of a radio, very long baseline interferometer (see, e.g., Walker 1995, for a VLBA description). With the sensitivity of current instrumentation, typically of a few tens of Jy/beam, the extended emission could be resolved, but only from the parts of the jet that are relatively close to its base. The bending of the jet appears more dramatic the smaller the inclination due to projection effects. The total received flux for each map is 7.4 mJy. As the flux increases proportionally to , one could detect the jet up to larger distances for higher . To illustrate this, Fig. 11 shows a comparison between two sky maps with different nonthermal energy fractions, and . We note that the radio fluxes also depend on as .
5 Summary and discussion
As argued by Bosch-Ramon & Barkov (2016), a nonballistic helical jet region is the likely outcome of a HMMQ jet interacting with the stellar wind under the effect of orbital motion. In the present work, we have shown that such a region is predicted to produce significant fluxes from radio to gamma rays if erg s, for a source at a few kiloparsecs. These significant fluxes are explained by the high synchrotron and IC efficiencies at the helical jet onset location, , for typical HMMQ -values. Since realistic -values are expected to be well below equipartition with the radiation field, the IC component is likely to dominate the nonthermal emission, peaking around 10 GeV. Specific gamma-ray light curve features are also predicted: a non-negligible impact of gamma-ray absorption, combined with angular effects for this process and IC. Peculiar changing radio morphologies, which can trace the jet helical structure, are expected as well.
We note that hydrodynamical instabilities and wind-jet mixing make our results mostly valid for the inner helical jet region, say up to a few . Beyond that point, the jet helical geometry is likely to become blurred, eventually turning into a bipolar, relatively wide, and collimated supersonic outflow; a mixture of jet and wind material. Although the flow may still be mildly relativistic for very energetic ejections, in general the resulting ”jet” is expected to be much faster than the stellar wind, but nonrelativistic.
Our results are based on a number of strong simplifications, in particular regarding the stability of the helical jet. Therefore, detailed numerical simulations of the jet-wind interaction on middle to large scales are necessary for more accurate predictions. Nevertheless, despite the simplifications adopted, the main features of the scenario are expected to be rather robust, at least at a semi-quantitative level. These are: non-negligible fluxes; specific light curves characterized by the value of , the interplay of jet and counter-jet emission, and angular IC and gamma-ray absorption effects; and curved jet radio morphologies with characteristic orbital evolution. This is so because the first- or even zeroth-order level of the dynamical, radiative, and geometrical effects involved are taken into account in our calculations.
Two important system parameters, eccentricity and , have been assumed constant in our calculations. Despite being difficult to ascertain the impact of eccentricity in a high-mass microquasar without numerical calculations, simulations of pulsar-star wind colliding in eccentric binaries (Barkov & Bosch-Ramon 2016; Bosch-Ramon et al. 2017) suggest that the effect of eccentricity should be important only for rather eccentric binaries. In that case, the nature of the interaction structure may not be helical at all, being not only strongly asymmetric, and inclined towards the apastron side of the orbit, but also likely much more sensitive to disruption. The value of , and thus , may also change. This would leave the helical jet geometry unaffected: both the helical vertical step and are . On the other hand, the luminosity would change as except for emission produced under the fast radiative cooling regime of electrons (in our setup the highest-energy end of the synchrotron and IC emission), which is independent of . The helical jet radiation in a very wide system would be hardly detectable unless the flow is slow, which allows electrons to radiate more energy. On the other hand, in a very compact system, radio emission would be hard to resolve for a slow helical jet flow, as the source would be too small.
The emission on the scales of the binary, produced for instance by jet internal shocks or by the first wind-induced recollimation shock in the jet (stage 2), could be important. In which case: (i) assuming that the same nonthermal luminosity is injected in both regions, the radio emission from the binary scales would be strongly absorbed via synchrotron self-absorption and wind free-free absorption; and (ii) the X- and gamma-ray emission would be at a similar level if produced in the fast radiation cooling regime, and at a higher level otherwise because of the higher -value and IC-target density. Gamma-ray absorption could however strongly attenuate the gamma-ray luminosity GeV for most of the orbit. Therefore, the helical jet could contribute significantly to the overall nonthermal radiation of HMMQ, and particularly in radio and GeV, even if electrons are also accelerated closer to the jet base. This was the main motivation for this work: to perform a first exploration of specific radiation features of the helical jet region, so that it could be disentangled from the other emitting sites. For the HMMQ Cygnus X-1 and Cygnus X-3, their uncertain wind and jet parameters make it difficult to make concrete predictions; for example, a non-ballistic jet region may not form at all (see Yoon et al. 2016; Bosch-Ramon & Barkov 2016, and references therein), and an individual detailed source analysis is out of the scope of this work. In spite of this, our findings are not in contradiction with the radio morphology and the gamma-ray light curves of these sources (e.g., Stirling et al. 2001; Mioduszewski et al. 2001; Miller-Jones et al. 2004; Zanin et al. 2016; Zdziarski et al. 2018). Future detailed modeling together with deep, high-resolution radio observations and the high-energy light curves could provide hints of helical jet emission.
The wind-jet interaction in HMMQ resembles, to a significant extent, the wind-wind interaction in high-mass binaries hosting young pulsars. The latter sources are probably extended mostly on the orbital plane and along the orbit semi-major axis, whereas the bipolar helical-jet structure is focused mostly in a direction perpendicular to the orbit plane (Bosch-Ramon et al. 2015; Bosch-Ramon & Barkov 2016; Bosch-Ramon et al. 2017). Nevertheless, the interactions of the stellar wind with a relativistic pulsar wind or jet, on binary, middle, and large scales, share many qualitative properties, and this may mask a fundamentally different engine (accretion vs. a pulsar wind) when observing sources of yet unknown CO. These similarities are: extended radio structures with orbital evolution found at mas scales; strong radio absorption on the binary scales; similar hydrodynamics, and therefore impact of instabilities, adiabatic losses, and Doppler boosting effects444The much higher energy per particle of shocked pulsar winds can lead to stronger Doppler boosting effects than in a HMMQ jet, although mass-loading smoothens this difference., for the shocked flows outside the binary region; and gamma-ray light curves strongly affected by IC and gamma-ray absorption angular effects, in a likely extended/multi-zone emitter. Radio observations, high-quality multiwavelength data, numerical simulations, and realistic radiation calculations (e.g., Aharonian et al. 2006; Perucho & Bosch-Ramon 2008; Albert et al. 2009; Fermi LAT Collaboration et al. 2009; Perucho et al. 2010; Moldón et al. 2011b, a, 2012; Perucho & Bosch-Ramon 2012; Bosch-Ramon et al. 2015; Dubus et al. 2015; Zanin et al. 2016; Yoon & Heinz 2015; Yoon et al. 2016; de la Cita et al. 2017; Bosch-Ramon et al. 2017; Wu et al. 2018; Barkov & Bosch-Ramon 2018) are therefore needed to characterize radiation scenarios in high-mass binaries with unknown CO, such that empirical observables can help to disentangle the CO nature.
Acknowledgements.We would like to thank the referee for a constructive and useful report that helped to improve the manuscript. We acknowledge support by the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grants AYA2013-47447-C3-1-P and AYA2016-76012-C3-1-P, with partial support by the European Regional Development Fund (ERDF/FEDER), MDM-2014-0369 of ICCUB (Unidad de Excelencia ‘María de Maeztu’), and the Catalan DEC grant 2017 SGR 643. EM acknowledges support from MINECO through grant BES-2016-076342.
- Abell & Margon (1979) Abell, G. O. & Margon, B. 1979, Nature, 279, 701
- Aharonian et al. (2006) Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 460, 743
- Aharonian (2004) Aharonian, F. A. 2004, Very high energy cosmic gamma radiation: a crucial window on the extreme Universe (World Scientific Publishing Co)
- Aharonian & Atoyan (1981) Aharonian, F. A. & Atoyan, A. M. 1981, Ap&SS, 79, 321
- Albert et al. (2009) Albert, J., Aliu, E., Anderhub, H., et al. 2009, ApJ, 693, 303
- Araudo et al. (2009) Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2009, A&A, 503, 673
- Barkov & Bosch-Ramon (2016) Barkov, M. V. & Bosch-Ramon, V. 2016, MNRAS, 456, L64
- Barkov & Bosch-Ramon (2018) Barkov, M. V. & Bosch-Ramon, V. 2018, MNRAS, 479, 1320
- Begelman et al. (2006) Begelman, M. C., King, A. R., & Pringle, J. E. 2006, MNRAS, 370, 399
- Bosch-Ramon (2009) Bosch-Ramon, V. 2009, A&A, 493, 829
- Bosch-Ramon (2013) Bosch-Ramon, V. 2013, in European Physical Journal Web of Conferences, Vol. 61, European Physical Journal Web of Conferences, 03001
- Bosch-Ramon & Barkov (2016) Bosch-Ramon, V. & Barkov, M. V. 2016, A&A, 590, A119
- Bosch-Ramon et al. (2017) Bosch-Ramon, V., Barkov, M. V., Mignone, A., & Bordas, P. 2017, MNRAS, 471, L150
- Bosch-Ramon et al. (2015) Bosch-Ramon, V., Barkov, M. V., & Perucho, M. 2015, A&A, 577, A89
- Bosch-Ramon & Khangulyan (2009) Bosch-Ramon, V. & Khangulyan, D. 2009, International Journal of Modern Physics D, 18, 347
- de la Cita et al. (2017) de la Cita, V. M., Bosch-Ramon, V., Paredes-Fortuny, X., Khangulyan, D., & Perucho, M. 2017, A&A, 598, A13
- Drury (1983) Drury, L. O. 1983, Reports on Progress in Physics, 46, 973
- Dubus et al. (2015) Dubus, G., Lamberts, A., & Fromang, S. 2015, A&A, 581, A27
- Fermi LAT Collaboration et al. (2009) Fermi LAT Collaboration, Abdo, A. A., Ackermann, M., et al. 2009, Science, 326, 1512
- Gould & Schréder (1967) Gould, R. J. & Schréder, G. P. 1967, Physical Review, 155, 1408
- Khangulyan et al. (2008) Khangulyan, D., Aharonian, F., & Bosch-Ramon, V. 2008, MNRAS, 383, 467
- Khangulyan et al. (2014) Khangulyan, D., Aharonian, F. A., & Kelner, S. R. 2014, ApJ, 783, 100
- Longair (1981) Longair, M. S. 1981, High energy astrophysics (Cambridge University Press)
- Luque-Escamilla et al. (2015) Luque-Escamilla, P. L., Martí, J., & Martínez-Aroza, J. 2015, A&A, 584, A122
- Melrose (1980) Melrose, D. B. 1980, Plasma astrohysics. Nonthermal processes in diffuse magnetized plasmas (Gordon & Breach)
- Miller-Jones et al. (2004) Miller-Jones, J. C. A., Blundell, K. M., Rupen, M. P., et al. 2004, ApJ, 600, 368
- Mioduszewski et al. (2001) Mioduszewski, A. J., Rupen, M. P., Hjellming, R. M., Pooley, G. G., & Waltman, E. B. 2001, ApJ, 553, 766
- Moldón et al. (2011a) Moldón, J., Johnston, S., Ribó, M., Paredes, J. M., & Deller, A. T. 2011a, ApJ, 732, L10
- Moldón et al. (2011b) Moldón, J., Ribó, M., & Paredes, J. M. 2011b, A&A, 533, L7
- Moldón et al. (2012) Moldón, J., Ribó, M., & Paredes, J. M. 2012, A&A, 548, A103
- Muijres et al. (2012) Muijres, L. E., Vink, J. S., de Koter, A., Müller, P. E., & Langer, N. 2012, A&A, 537, A37
- Owocki et al. (2009) Owocki, S. P., Romero, G. E., Townsend, R. H. D., & Araudo, A. T. 2009, ApJ, 696, 690
- Pacholczyk (1970) Pacholczyk, A. G. 1970, Radio astrophysics. Nonthermal processes in galactic and extragalactic sources (W. H. Freeman & Company)
- Perucho & Bosch-Ramon (2008) Perucho, M. & Bosch-Ramon, V. 2008, A&A, 482, 917
- Perucho & Bosch-Ramon (2012) Perucho, M. & Bosch-Ramon, V. 2012, A&A, 539, A57
- Perucho et al. (2010) Perucho, M., Bosch-Ramon, V., & Khangulyan, D. 2010, A&A, 512, L4
- Romero & Orellana (2005) Romero, G. E. & Orellana, M. 2005, A&A, 439, 237
- Romero et al. (2003) Romero, G. E., Torres, D. F., Kaufman Bernadó, M. M., & Mirabel, I. F. 2003, A&A, 410, L1
- Sell et al. (2010) Sell, P. H., Heinz, S., Calvelo, D. E., et al. 2010, ApJ, 719, L194
- Stirling et al. (2002) Stirling, A. M., Jowett, F. H., Spencer, R. E., et al. 2002, MNRAS, 337, 657
- Stirling et al. (2001) Stirling, A. M., Spencer, R. E., de la Force, C. J., et al. 2001, MNRAS, 327, 1273
- Walker (1995) Walker, R. C. 1995, in Astronomical Society of the Pacific Conference Series, Vol. 82, Very Long Baseline Interferometry and the VLBA, ed. J. A. Zensus, P. J. Diamond, & P. J. Napier, 133
- Wu et al. (2018) Wu, Y. W., Torricelli-Ciamponi, G., Massi, M., et al. 2018, MNRAS, 474, 4245
- Yoon & Heinz (2015) Yoon, D. & Heinz, S. 2015, ApJ, 801, 55
- Yoon et al. (2016) Yoon, D., Zdziarski, A. A., & Heinz, S. 2016, MNRAS, 456, 3638
- Zanin et al. (2016) Zanin, R., Fernández-Barral, A., de Oña Wilhelmi, E., et al. 2016, A&A, 596, A55
- Zdziarski et al. (2018) Zdziarski, A. A., Malyshev, D., Dubus, G., et al. 2018, MNRAS, 479, 4399