Search

Galactic cosmic rays (GCRs) constitute a significant part of the energy budget of our Galaxy, and the study of their accelerators is of high importance in modern astrophysics. Their main sources are likely supernova remnants (SNRs). These objects are capable to convert a part of their mechanical energy into accelerated charged particles. However, even though the mechanical energy reservoir of SNRs is promising, a conversion rate into particle energy of 10 to 20% is necessary to feed the population of GCRs. Such an efficiency is however not guaranteed. Complementary sources deserve thus to be investigated. This communication aims to address the question of the contribution to the acceleration of GCRs by pre-supernova massive stars in binary or higher multiplicity systems, Comment: 8 pages, 1 figure, Published in the Proceedings of the 3rd BINA Workshop on the Scientific Potential of the Indo-Belgian Cooperation

We model the currently available $\gamma$-ray data from the Fermi Large Area Telescope on Cyg X-3. Thanks to its very strong $\gamma$-ray activity during 2018--2021, the data quality has significantly improved. We study the strong orbital modulation of the $\gamma$-rays observed during at high $\gamma$-ray fluxes. The modulation, as found earlier, is well modeled by anisotropic Compton scattering of the donor blackbody emission by relativistic electrons in a jet strongly misaligned with respect to the orbital axis. We confirm that this model fits well both the average $\gamma$-ray modulation light curve and the spectrum. However, we find that if the jet is aligned with the spin axis of a rotating black hole, it would undergo geodetic precession with the period of $\sim$50 years. However, its presence is ruled out by both the $\gamma$-ray and radio data. Therefore, we consider an alternative model in which the average jet direction jet is aligned, but it is bent to outside the orbit owing to the thrust of the donor stellar wind, and thus precesses at the orbital period. The $\gamma$-ray modulation appears then owing to the variable Doppler boosting of synchrotron self-Compton jet emission. The model also fits well the data. However, the fitted bending angle is much larger than the theoretical one based on the binary and wind parameters as currently known. Thus, both models disagree with important aspects of our current theoretical understanding of the system. We discuss possible ways to find the correct model., Comment: ApJ, in press

We aim to interpret the radio emission of stellar bow shocks and assess under which conditions it could be either thermal or non-thermal, and how to use the observational data to infer their physical properties. We used an extended non-thermal emission model for stellar bow shocks for which we incorporated a consistent calculation of the thermal emission from the forward shock. We fitted this model to the available radio data, including largely unexplored data at low frequencies. In addition, we used a simplified one-zone model to estimate the gamma-ray emission from particles escaping the bow shocks. We can only explain the radio data from (BD+43{\deg}3654 and BD+60{\deg}2522) assuming a hard electron energy distribution below $\sim$1 GeV, a high efficiency of conversion of (shocked) wind kinetic power into relativistic electrons, and a relatively high magnetic-to-thermal pressure ratio. In the other systems, the interpretation of the observed flux density is ambiguous, although a non-thermal scenario is also favoured. We also show how complementary observations at other frequencies can allow us to place stronger constraints in the model. We also estimated the gamma-ray fluxes from the HII regions around the bow shocks of BD+43{\deg}3654 and BD+60{\deg}2522, and obtained luminosities at GeV energies of $\sim 10^{33}$~erg\,s$^{-1}$ and $10^{32}$~erg\,s$^{-1}$, respectively. Stellar bow shocks can potentially be very efficient particle accelerators. This work provides multi-wavelength predictions of their emission and demonstrates the key role of low-frequency radio observations in unveiling particle acceleration processes. The prospects of detections with next-generation observatories such as SKA and ngVLA are promising. Finally, BD+43{\deg}3654 may be detected in GeV in the near future, while bow shocks in general may turn out to be non-negligible sources of low-energy cosmic rays., Comment: 14 pages, 12 figures, 2 tables. Accepted for publication in Section 2 of Astronomy and Astrophysics: Astrophysical processes

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. We carried out a first dedicated exploration of the role of supernovae inside extragalactic jets in the production of ultra high-energy cosmic rays. 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. 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 $\gtrsim 10$~EeV. 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., Comment: 4 pages, Letter accepted for publication in Astronomy and Astrophysics (in press)

We report on the presence of very rapid hard X-ray variability in the $\gamma$-ray binary LS I +61 303. The results were obtained by analysing NuSTAR data, which show two achromatic strong flares on ks time-scales before apastron. The Swift-BAT orbital X-ray light curve is also presented, and the NuSTAR data are put in the context of the system orbit. The spectrum and estimated physical conditions of the emitting region indicate that the radiation is synchrotron emission from relativistic electrons, likely produced in a shocked pulsar wind. The achromaticity suggests that losses are dominated by escape or adiabatic cooling in a relativistic flow, and the overall behaviour in hard X-rays can be explained by abrupt changes in the size of the emitting region and/or its motion relative to the line of sight, with Doppler boosting potentially being a prominent effect. The rapid changes of the emitter could be the result of different situations such as quick changes in the intra-binary shock, variations in the re-accelerated shocked pulsar wind outside the binary, or strong fluctuations in the location and size of the Coriolis shock region. Although future multi-wavelength observations are needed to further constrain the physical properties of the high-energy emitter, this work already provides important insight into the complex dynamics and radiation processes in LS I +61-303., Comment: 8 pages, 7 figures. Accepted for publication in MNRAS

Gamma-ray binary systems, a subclass of high-mass X-ray binaries, show non-thermal emissions from radio to TeV. While efficient electron acceleration is considered to take place in them, the nature of the acceleration mechanism and the physical environments in these systems have been a long-standing question. In this work, we report on long-term recurrent patterns in the short-term variability of the soft X-ray emission of LS 5039, one of the brightest gamma-ray binary systems. The Neutron star Interior Composition Explorer (NICER) observed LS 5039 four times from 2018 to 2021. By comparing them with the previous Suzaku and NuSTAR long-exposure observations, we studied the long-term evolution of the orbital light curve in the soft X-ray band. Although the observations by NICER and Suzaku are separated by $\sim$14 years, i.e., more than 10^3 orbits, the orbital light curves show remarkable consistency after calculating their running averages with a window width 70 ks. Furthermore, all of the light curves show short-term variability with a time scale of $\sim$10 ks. Since the column density did not vary when the flux changed abruptly, such a short-term variability seems to be an intrinsic feature of the X-ray emission. We propose that the short-term variability is caused by clumps (or inhomogeneities) of the companion star wind impacting the X-ray production site. The observed time scale matches well with the lifetime of the clumps interacting with the pulsar wind and the dynamical time scale of the relativistic intrabinary shock in the pulsar wind scenario., Comment: 17 pages, 13 figures, accepted for publication in ApJ

Context: Massive colliding-wind binaries (CWBs) can be non-thermal sources. The emission produced in their wind-collision region (WCR) encodes information of both the shocks properties and the relativistic electrons accelerated in them. The recently discovered system Apep, a unique massive system hosting two Wolf-Rayet stars, is the most powerful synchrotron radio emitter among the known CWBs, being an exciting candidate to investigate the non-thermal processes associated with stellar wind shocks. Aims: We intend to break the degeneracy between the relativistic particle population and the magnetic field strength in the WCR of Apep by probing its hard X-ray spectrum, where inverse-Compton (IC) emission is expected to dominate. Methods: We observe Apep with NuSTAR for 60 ks and combine this with a re-analysis of a deep archival XMM-Newton observation to better constrain the X-ray spectrum. We use a non-thermal emission model to derive physical parameters from the results. Results: We detect hard X-ray emission consistent with a power-law component. This is compatible with IC emission produced in the WCR for a magnetic field of 100-160 mG and a fraction of ~1.5e-4 of the total wind kinetic power being converted into relativistic electron acceleration. Conclusions: This is the first time that the non-thermal emission from a CWB is detected both in radio and high energies. This allows us to derive the most robust constraints of the particle acceleration efficiency and magnetic field intensity in a CWB so far, reducing the typical uncertainty of a few orders of magnitude to just within a factor of two. This constitutes an important step forward in our characterisation of the physical properties of CWBs., Comment: 12 pages, 6 figures, accepted for publication in A&A

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. 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. 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. 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. 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., Comment: 10 pages, 9 figures Accepted in A&A on October 23, 2022

Should Primordial Black Holes (PBHs) exist in nature, they would inevitably accrete baryonic matter in their vicinity. In turn, the consequent emission of high-energy radiation could affect the thermal history of the universe to an extent that can be probed with a number of cosmological observables such as the Cosmic Microwave Background (CMB) anisotropies. However, our understanding of the accretion and radiation emission processes in the context of PBHs is still in its infancy, and very large theoretical uncertainties affect the resulting constraints on the PBH abundance. Building on state-of-the-art literature, in this work we take a step towards the development of a more realistic picture of PBH accretion by accounting for the contribution of outflows. Specifically, we derive CMB-driven constraints on the PBH abundance for various accretion geometries, ionization models and mass distributions in absence and in presence of mechanical feedback and non-thermal emissions due to the outflows. As a result, we show that the presence of such outflows introduces an additional layer of uncertainty that needs to be taken into account when quoting cosmological constraints on the PBH abundance, with important consequences in particular in the LIGO-Virgo-KAGRA observational window., Comment: 34 pages, 8 figures

MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA., Comment: 15 pages, 5 figures. Accepted for publication in MNRAS

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 modelling of the outflow interactions, spanning a broad range of spatial and temporal scales, is needed. Full 3-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 3-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 $\gtrsim 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., Comment: 18 pages, 11 figures

Context. There is a population of runaway stars that move at extremely high speeds with respect to their surroundings. The fast motion and the stellar wind of these stars, plus the wind-medium interaction, can lead to particle acceleration and non-thermal radiation. Aims. We characterise the interaction between the winds of fast runaway stars and their environment, in particular to establish their potential as cosmic-ray accelerators and non-thermal emitters. Methods. We model the hydrodynamics of the interaction between the stellar wind and the surrounding material. We self-consistently calculate the injection and transport of relativistic particles in the bow shock using a multi-zone code, and compute their broadband emission from radio to $\gamma$-rays. Results. Both the forward and reverse shocks are favourable sites for particle acceleration, although the radiative efficiency of particles is low and therefore the expected fluxes are in general rather faint. Conclusions. We show that high-sensitivity observations in the radio band can be used to detect the non-thermal radiation associated with bow shocks from hypervelocity and semi-relativistic stars. Hypervelocity stars are expected to be modest sources of sub-TeV cosmic rays, accounting perhaps for a $\sim 0.1$% of that of galactic cosmic rays., Comment: 11 pages, 8 figures, accepted for publication section 2. Astrophysical processes of Astronomy and Astrophysics

Gamma Cygni is a young supernova remnant located in the Cygnus region. MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes detected TeV emission (MAGIC J2019+408) to the north-west of this remnant, about 5 arcmin from its border. We want to identify the radio sources within the region encompassing gamma Cygni and MAGIC J2019+408 to shed light on their nature and investigate if these radio sources could be potential contributors to gamma-ray emission. We carried out a detailed study of the data we obtained with a survey of the Cygnus region at 325 and 610 MHz with the Giant Metrewave Radio Telescope (GMRT). We detected several radio sources in the region where the radio and the TeV emission overlap, as well as several areas of enhanced radio emission. In particular, two of these areas of diffuse enhanced emission may correspond to the supernova remnant interacting with a high density region, which seems to be the best candidate for the MAGIC source. Another two radio sources, which may or may not contribute to the gamma rays, are also spatially coincident with the emission peak of the MAGIC TeV source. One of them displays a rather peculiar extended morphology whose nature is completely unknown. We have identified the radio sources overlapping gamma Cygni and MAGIC J2019+408 and have shown that their potential gamma-ray contribution is likely not dominant. In addition, some of the studied sources show peculiar physical characteristics that deserve deeper multi-wavelength observations., Comment: 8 pages, 6 figures, Accepted for publication in Astronomy & Astrophysics

The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic fields during the first hundreds of seconds of jet propagation, focusing on the magnetized flow dynamics and the mechanisms of energy conversion. We have developed the code L\'ostrego v1.0, a new 3D RMHD code to simulate astrophysical plasmas in Cartesian coordinates. Using this tool, we performed the first 3D RMHD numerical simulations of relativistic magnetized jets propagating through the clumpy stellar wind in a HMXB. The overall morphology and dynamics of weakly magnetized jet models is similar to previous hydrodynamical simulations, where the jet head generates a strong shock in the ambient medium and the initial over-pressure with respect to the stellar wind drives one or more recollimation shocks. In the time scales of our simulations, these jets are ballistic and seem to be more stable against internal instabilities than jets with the same power in the absence of fields. However, moderate-to-strong toroidal magnetic fields favour the development of current-driven instabilities and the disruption of the jet within the binary. A detailed analysis of the energy distribution in the relativistic outflow and the ambient medium reveals that both magnetic and internal energies can contribute to the effective acceleration of the jet. We certify that the jet feedback into the ambient medium is highly dependent on the jet energy distribution at injection, where hotter, more dilute and/or more magnetized jets are more efficient, as anticipated by feedback studies in the case of jets in active galaxies., Comment: 28 pages, 23 figures. Abstract shortened for arXiv. Accepted for publication in A&A

Isolated stellar-mass BH accrete gas, often at supersonic speeds, and can form outflows that can influence the accreted gas. The latter process, known as mechanical feedback, can significantly affect the accretion rate. We use hydrodynamical simulations to assess the impact of mechanical feedback on the accretion rate when the BH moves supersonically through a uniform medium. We carry out 3D hydrodynamical simulations of outflows fueled by accretion that interact with a uniform medium, probing scales similar 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 are run for different outflow-medium motion angles and velocity ratios. The impact of different degrees of outflow collimation, accretor size, and resolution is also investigated. In general, the accretion rate is significantly affected by mechanical feedback. The effect is small for outflows perpendicular to the medium motion, and quickly grows for smaller angles. Moreover, the smaller the medium-to-outflow velocity ratio is, the more accretion decreases. On the other hand, the impact of outflow collimation seems moderate. The effect is enhanced when the accretor size is reduced. For a population of BH with random outflow orientations, the average accretion rate drops by (high-low 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. Our results strongly indicate that, on the considered scales, mechanical feedback can easily reduce the energy available from supersonic accretion by a factor of several. This should be taken into account when studying the mechanical, thermal and non-thermal output of isolated BH., Comment: 11 pages, 19 figures, astronomy and astrophysics, version after proofs

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\sim 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 $\sim 1$~AU., Comment: 16 pages, 13 figures, MNRAS accepted

The recently discovered massive binary system Apep is the most powerful synchrotron emitter among the known Galactic colliding-wind binaries. This makes this particular system of great interest to investigate stellar winds and the non-thermal processes associated with their shocks. This source was detected at various radio bands, and in addition the wind-collision region was resolved by means of very-long baseline interferometric observations. We use a non-thermal emission model for colliding-wind binaries to derive physical properties of this system. The observed morphology in the resolved maps allows us to estimate the system projection angle on the sky to be $\psi \approx 85^\circ$. The observed radio flux densities also allow us to characterise both the intrinsic synchrotron spectrum of the source and its modifications due to free--free absorption in the stellar winds at low frequencies; from this we derive mass-loss rates of the stars of $\dot{M}_\mathrm{WN} \approx 4\times10^{-5}$ $\mathrm{M}_\odot$ yr$^{-1}$ and $\dot{M}_\mathrm{WC} \approx 2.9\times10^{-5}$ $\mathrm{M}_\odot$ yr$^{-1}$. Finally, the broadband spectral energy distribution is calculated for different combinations of the remaining free parameters, namely the intensity of the magnetic field and the injected power in non-thermal particles. We show that the degeneracy of these two parameters can be solved with observations in the high-energy domain, most likely in the hard X-rays but also possibly in $\gamma$-rays under favourable conditions., Comment: 10 pages, 5 figures, 2 tables. This version has been accepted for publication in PASA after improvements and updates in the discussion of the results

Several supersonic runaway pulsar wind nebulae (sPWNe) with jet-like extended structures have been recently discovered in X-rays. If these structures are the product of electrons escaping the system and diffusing into the surrounding interstellar medium, they can produce a radio halo extending for several arcmin around the source. We model the expected radio emission in this scenario in the Lighthouse Nebula sPWN. We assume a constant particle injection rate during the source lifetime, and isotropic diffusion into the surrounding medium. Our predictions strongly depend on the low- and high-energy cutoffs given in the particle distribution. Our results indicate that extended radio emission can be detected from the Lighthouse Nebula without the need to invoke extreme values for the model parameters. We provide synthetic synchrotron maps that can be used to constrain these results with observations by current highly sensitive radio instruments., Comment: Accepted for publication in A&A

Binary systems of Population III can evolve to microquasars when one of the stars collapses into a black hole. When the compact object accretes matter at a rate greater than the Eddington rate, powerful jets and winds driven by strong radiation pressure should form. We investigate the structure of the jet-wind system for a model of Population III microquasar on scales beyond the jet-wind formation region. Using relativistic hydrodynamic simulations we find that the ratio of kinetic power between the jet and the disk wind determines the configuration of the system. When the power is dominated by the wind, the jet fills a narrow channel, collimated by the dense outflow. When the jet dominates the power of the system, part of its energy is diverted turning the wind into a quasi-equatorial flow, while the jet widens. From the results of our simulations, we implement semi-analytical calculations of the impact of the quasiequatorial wind on scales of the order of the size of the binary system. Our results indicate that Population III microquasars might inject gamma rays and relativistic particles into the early intergalactic medium, contributing to its reionization at large distances from the binary system., Comment: 12 pages, 9 figures, accepted for publication in Astrophysical and Space Science

LS~5039 is a powerful gamma-ray binary that probably hosts a non-accreting pulsar. Despite the wealth of data available, the power source of the non-thermal emitter is still unknown. We use a dynamical-radiative numerical model and multiwavelength data to constrain the properties of a pulsar wind that may power the non-thermal emitter in LS~5039. We ran simulations of an ultrarelativistic (low-$B$) cold $e^\pm$-wind that Compton scatters stellar photons and that dynamically interacts with the stellar wind. The effects of energy losses on the unshocked $e^\pm$-wind dynamics, and the geometry of the two-wind contact discontinuity, are computed for different wind models. The predicted unshocked $e^\pm$-wind radiation at periastron, when expected to be highest, is compared to LS~5039 data. The minimum possible radiation from an isotropic cold $e^\pm$-wind overpredicts the X-ray to gamma-ray fluxes at periastron by a factor of $\sim 3$. In the anisotropic wind case X-ray and $\gtrsim 100$ MeV data are not violated by wind radiation if the wind axis is at $\lesssim 20-40^\circ$ from the line of sight (probability of $\lesssim 6-24$\%), depending on the anisotropic wind model, or if the wind Lorentz factor $\in 10^2-10^3$, in which case the wind power can be higher, but it requires $e^\pm$-multiplicities of $\sim 10^6$ and $10^9$ for a $10^{-2}$~s and 10~s pulsar period, respectively. The studied model predicts that a low-$B$ cold pulsar $e^\pm$-wind in LS~5039 should be strongly anisotropic, with either a wind Lorentz factor $\in 10^2-10^3$ and very high multiplicities or with a fine-tuned wind orientation. A low-$B$, cold baryon-dominated wind would be possible, but then the multiplicities should be rather low, while the baryon-to-$e^\pm$ energy transfer should be very efficient at wind termination. A strongly magnetized cold wind seems to be the most favorable (least constrained) option., Comment: 6 pages, 2 figures, version after proofs, published in Astronomy and Astrophysics (abstract abridged), Eq. 10 corrected

Gamma-ray binaries are systems that emit non-thermal 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. We develop a semi-analytical model of the non-thermal emission of the colliding-wind structure including the dynamical effects of orbital motion. We apply the model to a general case and to LS 5039. 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, accounting for Doppler boosting and absorption processes associated to the massive star. 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. 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. The model is quite successful in reproducing the overall non-thermal behavior of LS 5039. Some improvements are suggested to better explain the phenomenology observed in this source, like accounting for particle reacceleration beyond the Coriolis turnover, unshocked pulsar wind emission, and the three-dimensional extension of the emitter., Comment: 12 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics

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 modelling of accretion onto these objects is required to obtain robust constraints on their abundance. We investigate the effect of mechanical feedback, i.e. the impact of outflows (winds and/or jets) on the medium, on primordial black hole accretion, and thereby on the associated radiation. Using analytical and numerical calculations, we study for the first time whether outflows can reduce the accretion rate of primordial black holes with masses similar to those detected by the LIGO-Virgo collaboration. Despite the complexity of the accretion rate evolution, mechanical feedback is able to significantly reduce the primordial black hole accretion rate, at least by an 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 still non-negligible. 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., Comment: 10 pages, 6 figures, to be published in Astronomy and Astrophysics (version after proofs; very minor changes)

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

Cosmic reionization put an end to the dark ages that came after the recombination era. Observations seem to favor the scenario where massive stars generating photons in low-mass galaxies were responsible for the bulk of reionization. Even though a possible contribution from accretion disks of active galactic nuclei (AGN) has been widely considered, they are currently thought to have had a minor role in reionization. Our aim is to study the possibility that AGN contributed to reionization not only through their accretion disks, but also through ionizing photons coming from the AGN jets interacting with the IGM. We adopt an empirically derived AGN luminosity function at $z\simeq6$, use X-ray observations to correct it for the presence of obscured sources, and estimate the density of jetted AGN. We then use analytical calculations to derive the fraction of jet energy that goes into ionizing photons. Finally, we compute the contribution of AGN jets to the H II volume filling factor at redshifts $z\simeq15-5$. We show that the contribution of the AGN jet lobes to the reionization of the Universe at $z\sim6$ might have been as high as $\gtrsim 10$\% of that of star-forming galaxies, under the most favorable conditions of jetted and obscuration fraction. The contribution of AGN to the reionization, while most likely not dominant, could have been higher than previously assumed, thanks to the radiation originated in the jet lobes., Comment: 7 pages, 3 figures, accepted for publication in A&A

The persistent gamma-ray source 4FGL J0647.7-4418 is tentatively associated in the latest Fermi catalogue with the sub-dwarf O-type X-ray binary HD 49798. However, an AGN candidate is also mentioned as an alternative identification in updated versions of the catalogue accompanying paper. If the first association were correct, this would add HD 49798 to the handful of currently known gamma-ray binaries and, therefore, represent a significant breakthrough not only because of a new member addition, but also because of the apparent white dwarf companion in this system. Despite these perspectives, here we show that the stellar association is likely wrong and that the proposed AGN object, well inside the Fermi 95% confidence ellipse, is a more conceivable counterpart candidate to the Fermi source due to its strong blazar similarities., Comment: 7 pages, 9 Figures, accepted for publication in MNRAS

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. 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. 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. Asymmetric light curves are obtained owing to the helical trajectory of the jets. 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., Comment: 11 pages, 12 figures, published in Astronomy & Astrophysics. v2: Typo corrected in Eq. 2

The innermost parts of powerful jets in active galactic nuclei are surrounded by dense, high-velocity clouds from the broad-line region, which may penetrate into the jet and lead to the formation of a strong shock. Such jet-cloud interactions are expected to have measurable effects on the $\gamma$-ray emission from blazars. We characterise the dynamics of a typical cloud-jet interaction scenario, and the evolution of its radiative output in the 0.1-30~GeV energy range, to assess to what extent these interactions can contribute to the $\gamma$-ray emission in blazars. We use semi-analytical descriptions of the jet-cloud dynamics, taking into account the expansion of the cloud inside the jet and its acceleration. Assuming that electrons are accelerated in the interaction and making use of the hydrodynamical information, we then compute the high-energy radiation from the cloud, including the absorption of $\gamma$-rays in the ambient photon field through pair creation. Jet-cloud interactions can lead to significant $\gamma$-ray fluxes in blazars with a broad-line region, in particular when the cloud expansion and acceleration inside the jet are taken into account. This is caused by 1) the increased shocked area in the jet, which leads to an increase in the energy budget for the non-thermal emission; 2) a more efficient inverse Compton cooling with the boosted photon field of the broad-line region; and 3) an increased observer luminosity due to Doppler boosting effects. For typical broad-line region parameters, either (i) jet-cloud interactions contribute significantly to the persistent $\gamma$-ray emission from blazars or (ii) the broad-line region is far from spherical or the fraction of energy deposited in non-thermal electrons is small., Comment: 9 pages, 9 figures. Accepted for publication in A&A

Blazars in elliptical hosts have a population of red giants surrounding their jet. These stars can carry large wind-blown bubbles into the jets, leading to gamma-ray emission through bubble-jet interactions. We study the interaction dynamics and the gamma-ray emission produced when the bubbles formed by red giant winds penetrate the jet of a blazar in an elliptical galaxy. First, we characterized the masses and penetration rates of the red giant wind bubbles that enter the jet. Then, the dynamical evolution of these bubbles under the jet impact was analysed analytically and numerically, and the radiation losses of the particles accelerated in the interaction were characterised. Finally, the synchrotron and the inverse Compton contributions above 100 MeV were estimated under different jet magnetic fields, powers, and Lorentz factors. We find that an analytical dynamical model is a reasonable approximation for the red giant wind bubble-jet interaction. The radiation produced by these wind bubbles interacting with a jet can have a duty cycle of up to 1. For realistic magnetic fields, gamma rays could be detectable from sources within the local universe, preferentially from those with high Lorentz factors (of about 10), and this could be a relatively common phenomenon for these sources. For magnetic fields in equipartition with the jet power, and high acceleration rates, synchrotron gamma rays may be detectable even for modest Lorentz factors (of about 3), but with a much lower duty cycle., Comment: 9 pages, 3 figures, 1 table, accepted for publication in A&A

Core-collapse supernovae are found in galaxies with ongoing star-formation. In a starburst galaxy hosting an active galactic nucleus with a relativistic jet, supernovae can take place inside the jet. The collision of the supernova ejecta with the jet flow is expected to lead to the formation of an interaction region, in which particles can be accelerated and produce high-energy emission. We study the non-thermal radiation produced by electrons accelerated as a result of a supernova explosion inside the jet of an active galactic nucleus within a star-forming galaxy. We first analyzed the dynamical evolution of the supernova ejecta impacted by the jet. Then, we explored the parameter space using simple prescriptions for the observed gamma-ray lightcurve. Finally, the synchrotron and the inverse Compton spectral energy distributions for two types of sources, a radio galaxy and a powerful blazar, are computed. For a radio galaxy, the interaction between a supernova and a jet of power $\sim 10^{43}-10^{44}$~erg~s$^{-1}$ can produce apparent gamma-ray luminosities of $\sim 10^{42}-10^{43}$~erg~s$^{-1}$, with an event duty cycle of supernova remnant (SNR) interacting with the jet close to one for one galaxy. For a blazar with a powerful jet of $\sim 10^{46}$~erg~s$^{-1}$, the jet-supernova ejecta interaction could produce apparent gamma-ray luminosities of $\sim 10^{43}-10^{44}$~erg~s$^{-1}$, but with a much lower duty cycle. The interaction of supernovae with misaligned jets of moderate power can be relatively frequent, and can result in steady gamma-ray emission potentially detectable for sources in the local universe. For powerful blazars much farther away, the emission would be steady as well, and it might be detectable under very efficient acceleration, but the events would be rather infrequent., Comment: 13 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Microquasars are binary systems consisting of a compact object and a star that present relativistic jets. When the companion is a high-mass star, significant interaction between the stellar wind and the jets is expected. In this work, an overview of the different effects that a strong stellar wind may have in microquasar jets is given. Both analytical estimates and numerical simulations performed in the last years are reviewed. The results of a model for the non-thermal emission of such jets at large scales are also shown and discussed. Observational studies of two high-mass microquasars, Cygnus X-1 and Cygnus X-3, are compared to the model predictions., Comment: 7 pages, 3 figures, proceedings of the Frontier Research in Astrophysics III workshop (FRAPWS 2018), Mondello, Palermo, Italy

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models., Comment: 9 pages, 3 figures. Accepted for publication in MNRAS the 2nd September 2018

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 suggest 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 lightcurve unless the jet is slow, $\beta < 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 on 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 $5\rm\,GeV$., Comment: 15 pages, 12 figures, MNRAS in press

The reionization of the Universe ends the dark ages that started after the recombination era. In the case of H, reionization finishes around $z\sim 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\sim 6-7$. We assumed that, as it has been proposed, active galactic nuclei at $z\sim 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\sim 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\sim 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\gtrsim 6$., Comment: Letter to the Editor accepted for publication in Astronomy & Astrophysics, very minor changes introduced

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. We calculate the broadband emission from a jet affected by the impact of the stellar wind and orbital motion in a high-mass microquasar. 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. 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. 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., Comment: 11 pages, 11 figures, accepted for publication in Astronomy & Astrophysics. [v2,v3]: typos corrected. [v4]: figures 2 and 3 corrected

Context: Runaway stars produce bowshocks that are usually observed at infrared (IR) wavelengths. Non-thermal radio emission has been detected so far only from the bowshock of BD+43{\deg}3654, whereas the detection of non-thermal radiation from these bowshocks at high energies remains elusive. Aims: We aim at characterising in detail the radio, X-ray, and gamma-ray emission from stellar bowshocks accounting for the structure of the region of interaction between the stellar wind and its environment. Methods: We develop a broadband-radiative, multi-zone model for stellar bowshocks that takes into account the spatial structure of the emitting region and the observational constraints. The model predicts the evolution and the emission of the relativistic particles accelerated and streaming together with the shocked flow. Results: We present broadband non-thermal spectral energy distributions for different scenarios, synthetic radio-cm synchrotron maps that reproduce the morphology of BD+43{\deg}3654, and updated predictions in X-ray and gamma-ray energy ranges. We also compare the results of the multi-zone model applied in this work with those of a refined one-zone model. Conclusions: A multi-zone model provides better constraints than a one-zone model on the relevant parameters, namely the magnetic field intensity and the amount of energy deposited in non-thermal particles. However, one-zone models can be improved by carefully characterising the intensity of the IR dust photon field and the escape rate of the plasma from the shocked region. Finally, comparing observed radio maps with those obtained from a multi-zone model enables constraints to be obtained on the direction of stellar motion with respect to the observer., Comment: 12 pages, 10 figures. Accepted for publication in A&A

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 quasi-plateau 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., Comment: 7 pages, 5 figures, submitted to MNRAS

The family links between radio galaxies and microquasars have been strongly strengthened thanks to a new common phenomenon: the presence of extended winged features. The first detection of such structures in a Galactic microquasar, recently reported in Nature Communications (http://rdcu.be/zgX8), widens the already known analogy between both kinds of outflow sources (Mart\'i et al. 2017). This observational result also has potential implications affecting the black hole merger scenarios that contribute to the gravitational wave background., Comment: 3 pages, 1 figure. Based on an invited talk presented to the workshop on "Interstellar Medium and Radio Astronomy in Catalonia: the Robert Estalella Legacy", held at the University of Barcelona on 21-22 December 2017

Gamma-ray observations of microquasars at high and very-high energies can provide valuable information of the acceleration processes inside the jets, the jet-environment interaction and the disk-jet coupling. Two high-mass microquasars have been deeply studied to shed light on these aspects: Cygnus X-1 and Cygnus X-3. Both systems display the canonical hard and soft X-ray spectral states of black hole transients, where the radiation is dominated by non-thermal emission from the corona and jets and by thermal emission from the disk, respectively. Here, we report on the detection of Cygnus X-1 above 60 MeV using 7.5 yr of Pass8 Fermi-LAT data, correlated with the hard X-ray state. A hint of orbital flux modulation was also found, as the source is only detected in phases around the compact object superior conjunction. We conclude that the high-energy gamma-ray emission from Cygnus X-1 is most likely associated with jets and its detection allow us to constrain the production site. Moreover, we include in the discussion the final results of a MAGIC long-term campaign on Cygnus X-1 that reaches almost 100 hr of observations at different X-ray states. On the other hand, during summer 2016, Cygnus X-3 underwent a flaring activity period in radio and high-energy gamma rays, similar to the one that led to its detection in the high-energy regime in 2009. MAGIC performed comprehensive follow-up observations for a total of about 70 hr. We discuss our results in a multi-wavelength context., Comment: Proceedings of the 35th International Cosmic Ray Conference (ICRC 2017), Bexco, Busan, Korea (arXiv:1708.05153)

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $\Gamma=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility., Comment: 15 pages, 5 figures. Accepted for publication in MNRAS

HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the only gamma-ray emitting high-mass binary known to host a non-accreting pulsar. In this Letter, we present hydrodynamical simulations of HESS J0632+057 in the context of a pulsar and a stellar wind interacting in an eccentric binary, and propose a scenario for the non-thermal phenomenology of the source. In this scenario, the non-thermal activity before and around apastron is linked to the accumulation of non-thermal particles in the vicinity of the binary, and the sudden drop of the emission before apastron is produced by the disruption of the two-wind interaction structure, allowing these particles to efficiently escape. In addition to providing a framework to explain the non-thermal phenomenology of the source, this scenario predicts extended, moving X-ray emitting structures similar to those observed in PSR B1259-63., Comment: 6 pages, 5 figures, accepted by MNRAS Letters

Dense stellar winds may mass-load the jets of active galactic nuclei, although it is unclear what are the time and spatial scales in which the mixing takes place. We study the first steps of the interaction between jets and stellar winds, and also the scales at which the stellar wind may mix with the jet and mass-load it. We present a detailed two-dimensional simulation, including thermal cooling, of a bubble formed by the wind of a star. We also study the first interaction of the wind bubble with the jet using a three-dimensional 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. The two-dimensional simulations point at quick wind bubble expansion and fragmentation after about one bubble shock crossing time. Three-dimensional simulations and stability analysis point at local mixing in the case of strong perturbations and relatively small 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 suggest highly unstable interaction structures and efficient wind-jet flow mixing on the scale of the jet interaction height, possibly producing strong inhomogeneities within the jet. In addition, the initial wind bubble shocked by the jet leads to a transient, large interaction surface. The interaction structure can be a source of significant non-thermal emission., Comment: Accepted for publication in Astronomy & Astrophysics

Fast radio bursts, or FRBs, are transient sources of unknown origin. Recent radio and optical observations have provided strong evidence for an extragalactic origin of the phenomenon and the precise localization of the repeating FRB 121102. Observations using the Karl G. Jansky Very Large Array (VLA) and very-long-baseline interferometry (VLBI) have revealed the existence of a continuum non-thermal radio source consistent with the location of the bursts in a dwarf galaxy. All these new data rule out several models that were previously proposed, and impose stringent constraints to new models. We aim to model FRB 121102 in light of the new observational results in the active galactic nucleus (AGN) scenario. We propose a model for repeating FRBs in which a non-steady relativistic $e^\pm$-beam, accelerated by an impulsive magnetohydrodynamic (MHD)-driven mechanism, interacts with a cloud at the centre of a star-forming dwarf galaxy. The interaction generates regions of high electrostatic field called cavitons in the plasma cloud. Turbulence is also produced in the beam. These processes, plus particle isotropization, the interaction scale, and light retardation effects, provide the necessary ingredients for short-lived, bright coherent radiation bursts. The mechanism studied in this work explains the general properties of FRB 121102, and may also be applied to other repetitive FRBs. Coherent emission from electrons and positrons accelerated in cavitons provides a plausible explanation of FRBs., Comment: 8 pages, 1 figure, accepted for publication in A&A

The central regions of galaxies are complex environments, rich in evolved and/or massive stars. For galaxies hosting an active galactic nucleus (AGN) with jets, the interaction of the jets with the winds of the stars within can lead to particle acceleration, and to extended high-energy emitting regions. We compute the non-thermal emission produced by the jet flow shocked by stellar winds on the jet scale, far from the jet-star direct interaction region. First, prescriptions for the winds of the relevant stellar populations in different types of galaxies are obtained. The scenarios adopted include galaxies with their central regions dominated by old or young stellar populations, and with jets of different power. Then, we estimate the available energy to accelerate particles in the jet shock, and compute the transport and energy evolution of the accelerated electrons, plus their synchrotron and inverse Compton emission, in the shocked flow along the jet. A significant fraction of the jet energy, $\sim 0.1-10$\%, can potentially be available for the particles accelerated in jet-wind shocks in the studied cases. The non-thermal particles can produce most of the high-energy radiation on jet scales, far from the jet shock region. This high-energy emission will be strongly enhanced in jets aligned with the line of sight due to Doppler boosting effects. The interaction of relativistic jets with stellar winds may contribute significantly to the persistent high-energy emission in some AGNs with jets. However, in the particular case of M87, this component seems too low to explain the observed gamma-ray fluxes., Comment: 15 pages, 11 figures. Accepted to be published in A&A

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., Comment: 9 pages, 7 figures, accepted for publication in A&A