Context.
LS 5039 is a powerful high-mass gamma-ray binary that probably hosts a young non-accreting pulsar. However, despite the wealth of data available, the means by which the non-thermal emitter ...is powered are still unknown.
Aims.
We use a dynamical-radiative numerical model, and multiwavelength data, to constrain the properties of a hypothetical pulsar wind that would power the non-thermal emitter in LS 5039.
Methods.
We ran simulations of an ultrarelativistic (weakly magnetized) cold
e
±
-wind that Compton scatters stellar photons and that dynamically interacts with the stellar wind. The effects of energy losses on the unshocked
e
±
-wind dynamics, and the geometry of the two-wind contact discontinuity, are computed for different wind models. The predicted unshocked
e
±
-wind radiation at periastron, when expected to be the highest, is compared to LS 5039 data.
Results.
The minimum possible radiation from an isotropic cold
e
±
-wind overpredicts the X-ray to gamma-ray fluxes at periastron by a factor of ∼3. In the anisotropic (axisymmetric) wind case X-ray and ≳100 MeV data are not violated by wind radiation if the wind axis is at ≲20−40° from the line of sight (chance probability of ≲6−24%), depending on the anisotropic wind model, or if the wind Lorentz factor ∈10
2
− 10
3
, in which case the wind power can be higher, but it requires
e
±
-multiplicities of ∼10
6
and 10
9
for a 10
−2
s and 10 s pulsar period, respectively.
Conclusions.
The studied model predicts that a weakly magnetized cold pulsar
e
±
-wind in LS 5039 should be strongly anisotropic, with either a wind Lorentz factor ∈10
2
− 10
3
and very high multiplicities or with a fine-tuned wind orientation. A weakly magnetized, cold baryon-dominated wind would be a possible alternative, but then the multiplicities should be rather low, while the baryon-to-
e
±
energy transfer should be very efficient at wind termination. A strongly magnetized cold wind seems to be the most favorable option as it is consistent with recent research on pulsar winds and does not require fine-tuning of the pulsar wind orientation, and the wind multiplicity and Lorentz factor are less constrained.
Context.
Dark matter may consist, at least partially, of primordial black holes formed during the radiation-dominated era. The radiation produced by accretion onto primordial black holes leaves ...characteristic signatures on the properties of the medium at high redshift, before and after hydrogen recombination. Therefore, reliable modeling of accretion onto these objects is required to obtain robust constraints on their abundance.
Aims.
We investigate the effect of mechanical feedback, that is, the impact of outflows (winds and– or –jets) on the medium, on primordial black hole accretion, and thereby on the associated radiation.
Methods.
Using analytical and numerical calculations, we studied for the first time the possibility that outflows can reduce the accretion rate of primordial black holes with masses similar to those detected by the LIGO-Virgo collaboration.
Results.
Despite the complexity of the accretion rate evolution, mechanical feedback is useful in to significantly reducing the primordial black hole accretion rate, at least by one order of magnitude, when outflows are aligned with the motion of the compact object. If the outflow is perpendicular to the direction of motion, the effect is less important, but it is still non-negligible.
Conclusions.
Outflows from primordial black holes, even rather weak ones, can significantly decrease the accretion rate, effectively weakening abundance constraints on these objects. Our results motivate further numerical simulations with a more realistic setup, which would yield more precise quantitative predictions.
Context.
High-mass gamma-ray binaries are powerful nonthermal galactic sources, some of them hosting a pulsar whose relativistic wind interacts with a likely inhomogeneous stellar wind. So far, ...modeling these sources including stellar wind inhomogeneities has been done using either simple analytical approaches or heavy numerical simulations, none of which allow for an exploration of the parameter space that is both reasonably realistic and general.
Aims.
Applying different semi-analytical tools together, we study the dynamics and high-energy radiation of a pulsar wind colliding with a stellar wind with different degrees of inhomogeneity to assess the related observable effects.
Methods.
We computed the arrival of clumps to the pulsar wind-stellar wind interaction structure using a Monte Carlo method and a phenomenological clumpy-wind model. The dynamics of the clumps that reach deep into the pulsar wind zone was computed using a semi-analytical approach. This approach allows for the characterization of the evolution of the shocked pulsar wind region in times much shorter than the orbital period. With this three-dimensional information about the emitter, we applied analytical adiabatic and radiative models to compute the variable high-energy emission produced on binary scales.
Results.
An inhomogeneous stellar wind induces stochastic hour-timescale variations in the geometry of the two-wind interaction structure on binary scales. Depending on the degree of stellar wind inhomogeneity, 10–100% level hour-scale variability in the X-rays and gamma rays is predicted, with the largest variations occurring roughly once per orbit.
Conclusions.
Our results, based on a comprehensive approach, show that present X-ray and future very-high-energy instrumentation can allow us to trace the impact of a clumpy stellar wind on the shocked pulsar wind emission in a gamma-ray binary.
Context.
Gamma-ray binaries are systems that emit nonthermal radiation peaking at energies above 1 MeV. One proposed scenario to explain their emission consists of a pulsar orbiting a massive star, ...with particle acceleration taking place in shocks produced by the interaction of the stellar and pulsar winds.
Aims.
We develop a semi-analytical model of the nonthermal emission of the colliding-wind structure, which includes the dynamical effects of orbital motion. We apply the model to a general case and to LS 5039.
Methods.
The model consists of a one-dimensional emitter, the geometry of which is affected by Coriolis forces owing to orbital motion. Two particle accelerators are considered: one at the two-wind standoff location and the other one at the turnover produced by the Coriolis force. Synchrotron and inverse Compton emission is studied taking into account Doppler boosting and absorption processes associated to the massive star.
Results.
If both accelerators are provided with the same energy budget, most of the radiation comes from the region of the Coriolis turnover and beyond, up to a few orbital separations from the binary system. Significant orbital changes of the nonthermal emission are predicted in all energy bands. The model allows us to reproduce some of the LS 5039 emission features, but not all of them. In particular, the MeV radiation is probably too high to be explained by our model alone, the GeV flux is recovered but not its modulation, and the radio emission beyond the Coriolis turnover is too low. The predicted system inclination is consistent with the presence of a pulsar in the binary.
Conclusions.
The model is quite successful in reproducing the overall nonthermal behavior of LS 5039. Some improvements are suggested to better explain the phenomenology observed in this source, such as accounting for particle reacceleration beyond the Coriolis turnover, unshocked pulsar wind emission, and the three-dimensional extension of the emitter.
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.
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.
Context. Galactic and extragalactic relativistic jets are surrounded by rich environments that are full of moving objects, such as stars and dense medium inhomogeneities. These objects can enter into ...the jets and generate shocks and non-thermal emission. Aims. We characterize the emitting properties of the downstream region of a standing shock formed due to the interaction of a relativistic jet with an obstacle. We focus on the case of red giants interacting with an extragalactic jet. Methods. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet meeting an obstacle of very large inertia. The results are interpreted in the framework of a red giant whose dense and slow wind interacts with the jet of an active galactic nucleus. Assuming that particles are accelerated in the standing shock generated in the jet as it impacts the red giant wind, we compute the non-thermal particle distribution, the Doppler boosting enhancement, and the non-thermal luminosity in gamma rays. Results. The available non-thermal energy from jet-obstacle interactions is potentially enhanced by a factor of ~100 when accounting for the whole surface of the shock induced by the obstacle, instead of just the obstacle section. The observer gamma-ray luminosity, including the effective obstacle size, the flow velocity and Doppler boosting effects, can be ~300 (γj/10)2 times higher than when the emitting flow is assumed at rest and only the obstacle section is considered, where γj is the jet Lorentz factor. For a whole population of red giants inside the jet of an active galactic nucleus, the predicted persistent gamma-ray luminosities may be potentially detectable for a jet pointing approximately to the observer. Conclusions. Obstacles interacting with relativistic outflows, for instance clouds and populations of stars for extragalactic jets, or stellar wind inhomogeneities in microquasar jets and in winds of pulsars in binaries, should be taken into account when investigating the origin of the non-thermal emission from these sources.
Context. The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary system suffers periodic and random variations of orbital and non-linear dynamical origins. ...The characterization of the evolution of the wind interaction region is necessary for understanding the rich phenomenology of these sources. Aims. For the first time, we simulate in 3 dimensions the interaction of isotropic stellar and relativistic pulsar winds along one full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolution and size, and of different state equations: a γ̂-constant ideal gas, and an ideal gas with γ̂ dependent on temperature. Methods. We used the code PLUTO to carry out relativistic hydrodynamical simulations in 2 and 3 dimensions of the interaction between a slow dense wind and a mildly relativistic wind with Lorentz factor 2, along one full orbit in a region up to ~100 times the binary size. The different 2-dimensional simulations were carried out with equal and larger grid resolution and size, and one was done with a more realistic equation of state than in 3 dimensions. Results. The simulations in 3 dimensions confirm previous results in 2 dimensions, showing: a strong shock induced by Coriolis forces that terminates the pulsar wind also in the opposite direction to the star; strong bending of the shocked-wind structure against the pulsar motion; and the generation of turbulence. The shocked flows are also subject to a faster development of instabilities in 3 dimensions, which enhances shocks, two-wind mixing, and large-scale disruption of the shocked structure. In 2 dimensions, higher resolution simulations confirm lower resolution results, simulations with larger grid sizes strengthen the case for the loss of the general coherence of the shocked structure, and simulations with two different equations of state yield very similar results. In addition to the Kelvin-Helmholtz instability, discussed in the past, we find that the Richtmyer-Meshkov and the Rayleigh-Taylor instabilities are very likely acting together in the shocked flow evolution. Conclusions. Simulations in 3 dimensions confirm that the interaction of stellar and pulsar winds yields structures that evolve non-linearly and become strongly entangled. The evolution is accompanied by strong kinetic energy dissipation, rapid changes in flow orientation and speed, and turbulent motion. The results of this work strengthen the case for the loss of the coherence of the whole shocked structure on large scales, although simulations of more realistic pulsar wind speeds are needed.
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 emission.