3D hydrodynamical simulations of the impact of mechanical feedback on accretion in supersonic stellar-mass black holes

Context. Isolated stellar-mass black holes accrete gas from their surroundings, often at supersonic speeds, and can form outflows that may influence the accreted gas. The latter process, known as mechanical feedback, can significantly affect the accretion rate. Aims. We use hydrodynamical simulations to assess the impact of mechanical feedback on the accretion rate when the black hole moves supersonically through a uniform medium. Methods. We carried out three-dimensional (3D) hydrodynamical simulations of outflows fueled by accretion that interact with a uniform medium, probing scales equivalent to and larger than the accretor gravitational sphere of influence. In the simulations, the accretor is at rest and the medium moves at supersonic speeds. The outflow power is assumed to be proportional to the accretion rate. The simulations were run for different outflow-medium motion angles and velocity ratios. We also investigated the impact of different degrees of outflow collimation, accretor size, and resolution. Results. In general, the accretion rate is significantly affected by mechanical feedback. There is a minor reduction in accretion for outflows perpendicular to the medium motion, but the reduction quickly becomes more significant for smaller angles. Moreover, the decrease in accretion becomes greater for smaller medium-to-outflow velocity ratios. On the other hand, the impact of outflow collimation seems moderate. Mechanical feedback is enhanced when the accretor size is reduced. For a population of black holes with random outflow orientations, the average accretion rate drops by (low-high resolution) ∼0.2−0.4 and ∼0.1−0.2 for medium-to-outflow velocity ratios of 1/20 and 1/100, respectively, when compared to the corresponding cases without outflow. Conclusions. Our results strongly indicate that on the considered scales, mechanical feedback can easily reduce the energy available from supersonic accretion by at least a factor of a few. This aspect should be taken into account when studying the mechanical, thermal, and non-thermal output of isolated black holes

The role of AGN jets in the reionization epoch

The reionization of the Universe ends the dark ages that started after the recombination era. In the case of H, reionization finishes around z ~ 6. Faint star-forming galaxies are the best candidate sources of the H-ionizing radiation, although active galactic nuclei may have also contributed. We have explored whether the termination regions of the jets from active galactic nuclei may have contributed significantly to the ionization of H in the late reionization epoch, around z ~ 6−7. We assumed that, as it has been proposed, active galactic nuclei at z ~ 6 may have presented a high jet fraction, accretion rate, and duty cycle, and that non-thermal electrons contribute significantly to the pressure of jet termination regions. Empirical black-hole mass functions were adopted to characterize the population of active galactic nuclei. From all this, estimates were derived for the isotropic H-ionizing radiation produced in the jet termination regions, at z ~ 6, through inverse Compton scattering off CMB photons. We find that the termination regions of the jets of active galactic nuclei may have radiated most of their energy in the form of H-ionizing radiation at z ~ 6. For typical black-hole mass functions at that redshift, under the considered conditions (long-lasting, common, and very active galactic nuclei with jets), the contribution of these jets to maintain (and possibly enhance) the ionization of H may have been non-negligible. We conclude that the termination regions of jets from active galactic nuclei could have had a significant role in the reionization of the Universe at z ≳ 6

The role of supernovae inside AGN jets in UHECR acceleration

Context. Jets of active galactic nuclei are potential accelerators of ultra high-energy cosmic rays. Supernovae can occur inside these jets and contribute to cosmic ray acceleration, particularly of heavy nuclei, but that contribution has been hardly investigated so far. Aims. We carried out a first dedicated exploration of the role of supernovae inside extragalactic jets in the production of ultra high-energy cosmic rays. Methods. We characterized the energy budget of supernova-jet interactions, and the maximum possible energies of the particles accelerated in those events, likely dominated by heavy nuclei. This allowed us to assess whether these interactions can be potential acceleration sites of ultra high-energy cosmic rays, or at least of their seeds. For that, we estimated the cosmic ray luminosity for different galaxy types, and compared the injection rate of cosmic ray seeds into the jet with that due to galactic cosmic ray entrainment. Results. Since the supernova is fueled for a long time by the luminosity of the jet, the energy of a supernova-jet interaction can be several orders of magnitude greater than that of an isolated supernova. Thus, despite the low rate of supernovae expected to occur in the jet, they could still provide more seeds for accelerating ultra high-energy particles than cosmic ray entrainment from the host galaxy. Moreover, these interactions can create sufficiently efficient accelerators to be a source of cosmic rays with energies ≳10 EeV. Conclusions. Supernova-jet interactions can contribute significantly to the production of ultra high-energy cosmic rays, either directly by accelerating these particles themselves or indirectly by providing pre-accelerated seeds

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

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 Lj∼1037(M˙w/10−6M⊙yr−1) erg s−1, where 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

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

A hydrodynamics-informed, radiation model for HESS J0632 + 057 from radio to gamma-rays

Relativistic hydrodynamical simulations of the eccentric gamma-ray binary HESS J0632+057 show that the energy of a putative pulsar wind should accumulate in the binary surroundings between periastron and apastron, being released by fast advection close to apastron. To assess whether this could lead to a maximum of the non-thermal emission before apastron, we derive simple prescriptions for the non-thermal energy content, the radiation efficiency, and the impact of energy losses on non-thermal particles, in the simulated hydrodynamical flow. These prescriptions are used to estimate the non-thermal emission in radio, X-rays, GeV, and TeV, from the shocked pulsar wind in a binary system simulated using a simplified 3-dimensional scheme for several orbital cycles. Lightcurves at different wavelengths are derived, together with synthetic radio images for different orbital phases. The dominant peak in the computed lightcurves is broad and appears close to, but before, apastron. This peak is followed by a quasiplateau shape, and a minor peak only in gamma rays right after periastron. The radio maps show ejection of radio blobs before apastron in the periastron-apastron direction. The results show that a scenario with a highly eccentric high-mass binary hosting a young pulsar can explain the general phenomenology of HESS J0632+057: despite its simplicity, the adopted approach yields predictions that are robust at a semi-quantitative level and consistent with multiwavelength observations

The Major Role of Eccentricity in the Evolution of Colliding Pulsar-Stellar Winds

Binary systems that host a massive star and a non-accreting pulsar can be powerful non-thermal emitters. The relativistic pulsar wind and the non-relativistic stellar outflows interact along the orbit, producing ultrarelativistic particles that radiate from radio to gamma rays. To properly characterize the physics of these sources, and better understand their emission and impact on the environment, careful modeling of the outflow interactions, spanning a broad range of spatial and temporal scales, is needed. Full three-dimensional approaches are very computationally expensive, but simpler approximate approaches, while still realistic at the semi-quantitative level, are available. We present here the results of calculations done with a quasi three-dimensional scheme to compute the evolution of the interacting flows in a region spanning in size up to a thousand times the size of the binary. In particular, we analyze for the first time the role of different eccentricities in the large scale evolution of the shocked flows. We find that the higher the eccentricity, the closer the flows behave like a one-side outflow, which becomes rather collimated for eccentricity values ≳0.75. The simulations also unveil that the pulsar and the stellar winds become fully mixed within the grid for low eccentricity systems, presenting a more stochastic behavior at large scales than in the highly eccentric systems

Inverse Compton emission from relativistic jets in binary systems

The gamma-ray emission detected from several microquasars can be produced by relativistic electrons emitting through inverse Compton scattering. In particular, the GeV emission detected from Cygnus X-3, and its orbital phase dependence, strongly suggests that the emitting electrons are accelerated in a relativistic jet, and that the optical companion provides the dominant target. Here, we study the effects related to particle transport in the framework of the relativistic jet scenario. We find that even in the most compact binary systems, with parameters similar to Cygnus X-3, particle transport can have a substantial influence on the GeV light curve unless the jet is slow, β < 0.7. In more extended binary systems, strong impact of particle transport is nearly unavoidable. Thus, even for a very compact system such as Cygnus X-3, particle transport significantly affects the ability of one-zone models to infer the properties of the gamma-ray production site based on the shape of the GeV lightcurve. We conclude that a detailed study of the gamma-ray spectrum can further constrain the structure and other properties of the gamma-ray emitter in Cygnus X-3, although such a study should account for gamma-gamma attenuation, since it may strongly affect the spectrum above 5GeV⁠

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 lightcurves 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

Impact of red giant/AGB winds on active galactic nucleus jet propagation

Context. Dense stellar winds may mass-load the jets of active galactic nuclei, although it is unclear on what time and spatial scales the mixing takes place. Aims. Our aim is to study the first steps of the interaction between jets and stellar winds, and also the scales on which the stellar wind mixes with the jet and mass-loads it. Methods. We present a detailed 2D simulation - including thermal cooling - of a bubble formed by the wind of a star designed to study the initial stages of jet-star interaction. We also study the first interaction of the wind bubble with the jet using a 3D simulation in which the star enters the jet. Stability analysis is carried out for the shocked wind structure to evaluate the distances over which the jet-dragged wind, which forms a tail, can propagate without mixing with the jet flow. Results.The 2D simulations point to quick wind bubble expansion and fragmentation after about one bubble shock crossing time. Three-dimensional simulations and stability analysis point to local mixing in the case of strong perturbations and relatively low density ratios between the jet and the jet dragged-wind, and to a possibly more stable shocked wind structure at the phase of maximum tail mass flux. Analytical estimates also indicate that very early stages of the star jet-penetration time may be also relevant for mass-loading. The combination of these and previous results from the literature suggests highly unstable interaction structures and efficient wind-jet flow mixing on the scale of the jet interaction height. Conclusions. The winds of stars with strong mass loss can efficiently mix with jets from active galactic nuclei. In addition, the initial wind bubble shocked by the jet leads to a transient, large interaction surface. The interaction between jets and stars can produce strong inhomogeneities within the jet. As mixing is expected to be effective on large scales, even individual asymptotic giant branch stars can significantly contribute to the mass-load of the jet and thus affect its dynamics. Shear layer mass-entrainment could be important. The interaction structure can be a source of significant non-thermal emissio