Showing 1–50 of 173 results for author: Bosch-Ramon, V

Investigating the Role of Pre-supernova Massive Stars in the Acceleration of Galactic Cosmic Rays

Authors: Michael De Becker, Santiago del Palacio, Paula Benaglia, Anandmayee Tej, Benito Marcote, Gustavo Esteban Romero, Valenti Bosch-Ramon, C. H. Ishwara-Chandra

Abstract: 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… ▽ More 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 △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 8 pages, 1 figure, Published in the Proceedings of the 3rd BINA Workshop on the Scientific Potential of the Indo-Belgian Cooperation

Journal ref: Bulletin de la Société Royale des Sciences de Liège, 2024, Volume 93, No 2, 536-543

Two models for the orbital modulation of $γ$-rays in Cyg X-3

Authors: Anton Dmytriiev, Andrzej A. Zdziarski, Denys Malyshev, Valenti Bosch-Ramon, Maria Chernyakova

Abstract: We model the currently available $γ$-ray data from the Fermi Large Area Telescope on Cyg X-3. Thanks to its very strong $γ$-ray activity during 2018--2021, the data quality has significantly improved. We study the strong orbital modulation of the $γ$-rays observed during at high $γ$-ray fluxes. The modulation, as found earlier, is well modeled by anisotropic Compton scattering of the donor blackbo… ▽ More We model the currently available $γ$-ray data from the Fermi Large Area Telescope on Cyg X-3. Thanks to its very strong $γ$-ray activity during 2018--2021, the data quality has significantly improved. We study the strong orbital modulation of the $γ$-rays observed during at high $γ$-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 $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: ApJ, submitted, comments are welcome

Probing the non-thermal physics of stellar bow shocks using radio observations

Authors: Javier Rodrigo Martinez, Santiago del Palacio, Valentí Bosch-Ramon

Abstract: 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 mode… ▽ More 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°3654 and BD+60°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°3654 and BD+60°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°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. △ Less

Submitted 28 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures, 2 tables. Accepted for publication in Section 2 of Astronomy and Astrophysics: Astrophysical processes

Journal ref: A&A 680, A99 (2023)

The role of supernovae inside AGN jets in UHECR acceleration

Authors: Valenti Bosch-Ramon

Abstract: 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 c… ▽ More 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. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 4 pages, Letter accepted for publication in Astronomy and Astrophysics (in press)

Journal ref: A&A 677, L14 (2023)

Achromatic rapid flares in hard X-rays in the $γ$-ray binary LS I +61-303

Authors: Enzo A. Saavedra, Gustavo E. Romero, Valenti Bosch-Ramon, Elina Kefala

Abstract: We report on the presence of very rapid hard X-ray variability in the $γ$-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 t… ▽ More We report on the presence of very rapid hard X-ray variability in the $γ$-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. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 8 pages, 7 figures. Accepted for publication in MNRAS

Unveiling properties of the non-thermal X-ray production in the gamma-ray binary LS 5039 using the long-term pattern of its fast X-ray variability

Authors: Hiroki Yoneda, Valenti Bosch-Ramon, Teruaki Enoto, Dmitry Khangulyan, Paul S. Ray, Tod Strohmayer, Toru Tamagawa, Zorawar Wadiasingh

Abstract: 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… ▽ More 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. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 17 pages, 13 figures, accepted for publication in ApJ

Evidence for non-thermal X-ray emission from the double WR colliding-wind binary Apep

Authors: S. del Palacio, F. García, M. De Becker, D. Altamirano, V. Bosch-Ramon, P. Benaglia, B. Marcote, G. E. Romero

Abstract: 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, bei… ▽ More 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. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 12 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 672, A109 (2023)

Modeling the effects of clumpy winds in the high-energy light curves of γ-ray binaries

Authors: E. Kefala, V. Bosch-Ramon

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 2 November, 2022; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 10 pages, 9 figures Accepted in A&A on October 23, 2022

Journal ref: A&A 669, A21 (2023)

The effect of outflows on CMB bounds from Primordial Black Hole accretion

Authors: Lorenzo Piga, Matteo Lucca, Nicola Bellomo, Valentì Bosch-Ramon, Sabino Matarrese, Alvise Raccanelli, Licia Verde

Abstract: 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 accret… ▽ More 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. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 34 pages, 8 figures

Report number: ULB-TH/22-14, UTWI-11-2022

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: 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 expe… ▽ More 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. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

The major role of eccentricity in the evolution of colliding pulsar-stellar winds

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: 18 pages, 11 figures

Journal ref: July 2021, Universe, vol. 7, issue 8, p. 277

Non-thermal emission in hyper-velocity and semi-relativistic stars

Authors: J. R. Martinez, S. del Palacio, V. Bosch-Ramon, G. E. Romero

Abstract: 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 the… ▽ More 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 $γ$-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. △ Less

Submitted 23 March, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 11 pages, 8 figures, accepted for publication section 2. Astrophysical processes of Astronomy and Astrophysics

Journal ref: A&A 661, A102 (2022)

Exploring the region encompassing gamma Cygni SNR and MAGIC J2019+408 with the GMRT at 325 and 610 MHz

Authors: Josep M. Paredes, Paula Benaglia, C. H. Ishwara-Chandra, Valenti Bosch-Ramon, Marcel Strzys

Abstract: 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… ▽ More 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. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: 8 pages, 6 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A73 (2022)

3D RMHD simulations of jet-wind interactions in High Mass X-ray Binaries

Authors: J. López-Miralles, M. Perucho, J. M. Martí, S. Migliari, V. Bosch-Ramon

Abstract: 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 field… ▽ More 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óstrego 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. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 28 pages, 23 figures. Abstract shortened for arXiv. Accepted for publication in A&A

Journal ref: A&A 661, A117 (2022)

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

Authors: Valenti Bosch-Ramon

Abstract: 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 hydro… ▽ More 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. △ Less

Submitted 28 March, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 11 pages, 19 figures, astronomy and astrophysics, version after proofs

Journal ref: A&A 660, A5 (2022)

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

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

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 stud… ▽ More 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. △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: 16 pages, 13 figures, MNRAS accepted

The non-thermal emission from the colliding-wind binary Apep

Authors: S. del Palacio, P. Benaglia, M. De Becker, V. Bosch-Ramon, G. E. Romero

Abstract: 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 v… ▽ More 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 $ψ\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 $γ$-rays under favourable conditions. △ Less

Submitted 17 December, 2021; v1 submitted 14 November, 2021; originally announced November 2021.

Comments: 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

Journal ref: Publ. Astron. Soc. Aust. 39 (2022) e004

Signatures of extended radio emission from escaping electrons in the Lighthouse Nebula

Authors: Pol Bordas, Xiying Zhang, Valentí Bosch-Ramon, Josep Maria Paredes

Abstract: 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 Lighth… ▽ More 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. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 654, A4 (2021)

Equatorial outflows driven by jets in Population III microquasars

Authors: Pablo Sotomayor Checa, Gustavo E. Romero, Valenti Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 12 pages, 9 figures, accepted for publication in Astrophysical and Space Science

Properties of a hypothetical cold pulsar wind in LS~5039

Authors: V. Bosch-Ramon

Abstract: 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-… ▽ More 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. △ Less

Submitted 16 April, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 6 pages, 2 figures, version after proofs, published in Astronomy and Astrophysics (abstract abridged), Eq. 10 corrected

Journal ref: A&A 645, A86 (2021)

A dynamical and radiation semi-analytical model of pulsar-star colliding winds along the orbit: Application to LS 5039

Authors: Edgar Molina, Valentí Bosch-Ramon

Abstract: 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 incl… ▽ More 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. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 12 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A84 (2020)

Mechanical feedback effects on primordial black hole accretion

Authors: Valenti Bosch-Ramon, Nicola Bellomo

Abstract: 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 constra… ▽ More 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. △ Less

Submitted 9 June, 2020; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: 10 pages, 6 figures, to be published in Astronomy and Astrophysics (version after proofs; very minor changes)

Journal ref: A&A 638, A132 (2020)

The high-energy emission from HD~93129A near periastron

Authors: S. del Palacio, F. García, D. Altamirano, R. H. Barbá, V. Bosch-Ramon, M. Corcoran, M. De Becker, K. Hamaguchi, J. Maíz Apellániz, P. Munar Adrover, J. M. Paredes, G. E. Romero, H. Sana, M. Tavani, A. ud-Doula

Abstract: 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 $γ$-ray satellite \textit{AGILE}… ▽ More 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 $γ$-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 $γ$-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. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

AGN jets versus accretion as reionization sources

Authors: N. Torres-Albà, V. Bosch-Ramon, K. Iwasawa

Abstract: 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 ro… ▽ More 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. △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: 7 pages, 3 figures, accepted for publication in A&A

Journal ref: A&A 635, A57 (2020)

A blazar as the likely counterpart to 4FGL J0647.7-4418 instead of a gamma-ray binary

Authors: Josep Martí, Estrella Sánchez-Ayaso, Pedro L. Luque-Escamilla, Josep M. Paredes, Valentí Bosch-Ramon, Robin H. D. Corbet

Abstract: 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 bina… ▽ More 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. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Comments: 7 pages, 9 Figures, accepted for publication in MNRAS

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

Authors: Edgar Molina, Santiago del Palacio, Valentí Bosch-Ramon

Abstract: 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 Dopp… ▽ More 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. △ Less

Submitted 5 May, 2021; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 11 pages, 12 figures, published in Astronomy & Astrophysics. v2: Typo corrected in Eq. 2

Journal ref: A&A 629, A129 (2019)

Gamma rays from jets interacting with BLR clouds in blazars

Authors: S. del Palacio, V. Bosch-Ramon, G. E. Romero

Abstract: 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 $γ$-ray emission from blazars. We characterise the dynamics of a typical cloud-jet interaction scenario, and t… ▽ More 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 $γ$-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 $γ$-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 $γ$-rays in the ambient photon field through pair creation. Jet-cloud interactions can lead to significant $γ$-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 $γ$-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. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 9 pages, 9 figures. Accepted for publication in A&A

Gamma rays from red giant wind bubbles entering the jets of elliptical host blazars

Authors: N. Torres-Albà, V. Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: 9 pages, 3 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 623, A91 (2019)

Non-thermal emission resulting from a supernova explosion inside an extragalactic jet

Authors: F. L. Vieyro, V. Bosch-Ramon, N. Torres-Albà

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: 13 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 622, A175 (2019)

The influence of the optical star on the jets of high-mass microquasars

Authors: Edgar Molina, Valentí Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 10 September, 2018; originally announced September 2018.

Comments: 7 pages, 3 figures, proceedings of the Frontier Research in Astrophysics III workshop (FRAPWS 2018), Mondello, Palermo, Italy

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee , et al. (133 additional authors not shown)

Abstract: 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 o… ▽ More 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-$σ$ 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 $σ$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 9 pages, 3 figures. Accepted for publication in MNRAS the 2nd September 2018

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 481, Issue 2, 1 December 2018, Pages 2479-2486

Inverse Compton Emission from Relativistic Jets in Binary Systems

Authors: D. Khangulyan, V. Bosch-Ramon, Y. Uchiyama

Abstract: 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… ▽ More 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, $β< 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$. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 15 pages, 12 figures, MNRAS in press

The role of AGN jets in the reionization epoch

Authors: V. Bosch-Ramon

Abstract: 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 con… ▽ More 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$. △ Less

Submitted 7 September, 2018; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: Letter to the Editor accepted for publication in Astronomy & Astrophysics, very minor changes introduced

Journal ref: A&A 617, L3 (2018)

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

Authors: Edgar Molina, Valentí Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 12 September, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: 11 pages, 11 figures, accepted for publication in Astronomy & Astrophysics. [v2,v3]: typos corrected. [v4]: figures 2 and 3 corrected

Journal ref: A&A 618, A146 (2018)

Multi-zone non-thermal radiative model for stellar bowshocks

Authors: S. del Palacio, V. Bosch-Ramon, A. L. Müller, G. E. Romero

Abstract: 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°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 accoun… ▽ More 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°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°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. △ Less

Submitted 1 July, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 12 pages, 10 figures. Accepted for publication in A&A

Journal ref: A&A 617, A13 (2018)

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

Authors: Maxim Barkov, Valenti Bosch-Ramon

Abstract: 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… ▽ More 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. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: 7 pages, 5 figures, submitted to MNRAS

GRS 1758-258: the first winged microquasar

Authors: Josep Martí, Pedro L. Luque-Escamilla, Valentí Bosch-Ramon, Josep M. Paredes

Abstract: 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í et al. 2017). This obser… ▽ More 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í et al. 2017). This observational result also has potential implications affecting the black hole merger scenarios that contribute to the gravitational wave background. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: 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 rays from microquasars Cygnus X-1 and Cygnus X-3

Authors: A. Fernández-Barral, O. Blanch, E. de Oña Wilhelmi, D. Galindo, J. Herrera, M. Ribó, J. Rico, A. Stamerra, F. Aharonian, V. Bosch-Ramon, R. Zanin

Abstract: 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… ▽ More 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. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: Proceedings of the 35th International Cosmic Ray Conference (ICRC 2017), Bexco, Busan, Korea (arXiv:1708.05153)

Report number: MAGIC-ICRC/2017/11

Journal ref: PoS(ICRC2017)734

Search for very-high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (122 additional authors not shown)

Abstract: 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 t… ▽ More 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 $Γ=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. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 472, Issue 3, 11 December 2017, Pages 3474-3485

HESS J0632+057: hydrodynamics and nonthermal emission

Authors: Valentí Bosch-Ramon, Maxim V. Barkov, Andrea Mignone, Pol Bordas

Abstract: 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… ▽ More 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. △ Less

Submitted 31 July, 2017; originally announced August 2017.

Comments: 6 pages, 5 figures, accepted by MNRAS Letters

The impact of red giant/AGB winds on AGN jet propagation

Authors: Manel Perucho, Valentí Bosch-Ramon, Maxim V. Barkov

Abstract: 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,… ▽ More 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. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

A model for the repeating FRB 121102 in the AGN scenario

Authors: Florencia L. Vieyro, Gustavo E. Romero, Valentí Bosch-Ramon, Benito Marcote, María V. del Valle

Abstract: 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-therma… ▽ More 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. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: 8 pages, 1 figure, accepted for publication in A&A

Journal ref: A&A 602, A64 (2017)

Collective non-thermal emission from an extragalactic jet interacting with stars

Authors: Florencia L. Vieyro, Núria Torres-Albà, Valentí Bosch-Ramon

Abstract: 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 je… ▽ More 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. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 15 pages, 11 figures. Accepted to be published in A&A

Journal ref: A&A 604, A57 (2017)

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

Authors: V. M. de la Cita, S. del Palacio, V. Bosch-Ramon, X. Paredes-Fortuny, G. E. Romero, D. Khangulyan

Abstract: 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.… ▽ More 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. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: 9 pages, 7 figures, accepted for publication in A&A

Non-thermal radiation from a pulsar wind interacting with an inhomogeneous stellar wind

Authors: Víctor M. de la Cita, Valentí Bosch-Ramon, Xavier Paredes-Fortuny, Dmitry Khangulyan, Manel Perucho

Abstract: Binaries hosting a massive star and a non-accreting pulsar are powerful non-thermal emitters due to the interaction of the pulsar and the stellar wind. The winds of massive stars are thought to be inhomogeneous, which could have an impact on the non-thermal emission. We study numerically the impact of the presence of inhomogeneities or clumps in the stellar wind on the high-energy non-thermal radi… ▽ More Binaries hosting a massive star and a non-accreting pulsar are powerful non-thermal emitters due to the interaction of the pulsar and the stellar wind. The winds of massive stars are thought to be inhomogeneous, which could have an impact on the non-thermal emission. We study numerically the impact of the presence of inhomogeneities or clumps in the stellar wind on the high-energy non-thermal radiation of high-mass binaries hosting a non-accreting pulsar. We compute the trajectories and physical properties of the streamlines in the shocked pulsar wind without clumps, with a small clump, and with a large one. This information is used to compute the synchrotron and inverse Compton emission from the non-thermal populations, accounting also for the effect of gamma-ray absorption through pair creation. A specific study is done for PSR B1259-63/LS2883. When stellar wind clumps perturb the two-wind interaction region, the associated non-thermal radiation in the X-ray band,of synchrotron origin, and in the GeV-TeV band, of inverse Compton origin, is affected by several effects: (i) strong changes in the the plasma velocity direction that result in Doppler boosting factor variations, (ii) strengthening of the magnetic field that mainly enhances the synchrotron radiation, (iii) strengthening of the pulsar wind kinetic energy dissipation at the shock, potentially available for particle acceleration, and (iv) changes in the rate of adiabatic losses that affect the lower energy part of the non-thermal particle population. The radiation above 100 GeV detected, presumably, during the post-periastron crossing of the Be star disc in PSR B1259-63/LS2883, can be roughly reproduced assuming that the crossing of the disc is modeled as the encounter with a large inhomogeneity. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 13 pages, 12 figures, submitted to A&A

Journal ref: A&A 598, A13 (2017)

Detection of gamma rays of likely jet origin in Cygnus X-1

Authors: Roberta Zanin, Alba Fernández-Barral, Emma de Oña-Wilhelmi, Felix Aharonian, Oscar Blanch, Valentí Bosch-Ramon, Daniel Galindo

Abstract: Aims: Probe the high-energy ($>$60 MeV) emission from the black hole X-ray binary system, Cygnus X-1, and investigate its origin. Methods: We analysed 7.5 yr of data by Fermi/LAT with the latest PASS8 software version. Results: We report the detection of a signal at $\sim$8 $σ$ statistical significance spatially coincident with Cygnus X-1 and a luminosity above 60 MeV of 5.5$\times$10$^{33}$ erg s… ▽ More Aims: Probe the high-energy ($>$60 MeV) emission from the black hole X-ray binary system, Cygnus X-1, and investigate its origin. Methods: We analysed 7.5 yr of data by Fermi/LAT with the latest PASS8 software version. Results: We report the detection of a signal at $\sim$8 $σ$ statistical significance spatially coincident with Cygnus X-1 and a luminosity above 60 MeV of 5.5$\times$10$^{33}$ erg s$^{-1}$. The signal is correlated with the hard X-ray flux: the source is observed at high energies only during the hard X-ray spectral state, when the source is known to display persistent, relativistic radio emitting jets. The energy spectrum, extending up to $\sim$20 GeV without any sign of spectral break, is well fitted by a power-law function with a photon index of 2.3$\pm$0.2. There is a hint of orbital flux variability, with high-energy emission mostly coming around the superior conjunction. Conclusions: We detected GeV emission from Cygnus X-1 and probed that the emission is most likely associated with the relativistic jets. The evidence of flux orbital variability points to the anisotropic inverse Compton on stellar photons as the mechanism at work, thus constraining the emission region to a distance $10^{11}-10^{13}$ cm from the black hole. △ Less

Submitted 19 September, 2016; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: accepted A\&A (9 pages and 7 figures)

Journal ref: A&A 596, A55 (2016)

A model for the non-thermal emission of the very massive colliding-wind binary HD 93129A

Authors: Santiago del Palacio, Valentí Bosch-Ramon, Gustavo E. Romero, Paula Benaglia

Abstract: The binary stellar system HD 93129A is one of the most massive known binaries in our Galaxy. This system presents non-thermal emission in the radio band, which can be used to infer its physical conditions and predict its emission in the high-energy band. We intend to constrain some of the unknown parameters of HD 93129A through modelling the non-thermal emitter, and also to analyse the detectabili… ▽ More The binary stellar system HD 93129A is one of the most massive known binaries in our Galaxy. This system presents non-thermal emission in the radio band, which can be used to infer its physical conditions and predict its emission in the high-energy band. We intend to constrain some of the unknown parameters of HD 93129A through modelling the non-thermal emitter, and also to analyse the detectability of this source in hard X-rays and $γ$-rays. We develop a broadband radiative model for the wind-collision region taking into account the evolution of the accelerated particles streaming along the shocked region, the emission by different radiative processes, and the attenuation of the emission propagating through the local matter and radiation fields. From the analysis of the radio emission, we find that the binary HD~93129A is more likely to have a low inclination and a high eccentricity. The minimum energy of the non-thermal electrons seems to be between $\sim 20 - 100$MeV, depending on the intensity of the magnetic field in the wind-collision region. The latter can be in the range $\sim 20 - 1500$ mG. Our model is able to reproduce the observed radio emission, and predicts that the non-thermal radiation from HD~93129A will increase in the near future. With instruments such as \textit{NuSTAR}, \textit{Fermi}, and CTA, it will be possible to constrain the relativistic particle content of the source, and other parameters such as the magnetic field strength in the wind collision zone, which in turn can be used to obtain upper-limits of the magnetic field on the surface of the very massive stars, thereby inferring whether magnetic field amplification is taking place in the particle acceleration region. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 11 pages, 9 figures. Accepted for publication in A&A

Journal ref: A&A 591, A139 (2016)

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

Authors: V. Bosch-Ramon, M. V. Barkov

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

Submitted 21 April, 2016; originally announced April 2016.

Comments: 7 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 590, A119 (2016)

Coupling hydrodynamics and radiation calculations for star-jet interactions in AGN

Authors: Víctor M. de la Cita, Valentí Bosch-Ramon, Xavier Paredes-Fortuny, Dmitry Khangulyan, Manel Perucho

Abstract: Stars and their winds can contribute to the non-thermal (NT) emission in extragalactic jets. Given the complexity of jet-star interactions, the properties of the resulting emission are strongly linked to those of the emitting flows. We simulate the interaction between a stellar wind and a relativistic extragalactic jet and use the hydrodynamic results to compute the NT emission under different con… ▽ More Stars and their winds can contribute to the non-thermal (NT) emission in extragalactic jets. Given the complexity of jet-star interactions, the properties of the resulting emission are strongly linked to those of the emitting flows. We simulate the interaction between a stellar wind and a relativistic extragalactic jet and use the hydrodynamic results to compute the NT emission under different conditions. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet interacting with a supersonic, non-relativistic stellar wind. We compute the corresponding streamlines out of the simulation results, and calculate the injection, evolution, and emission of NT particles accelerated in the jet shock, focusing on electrons or $e^\pm$-pairs. Several cases are explored, considering different jet-star interaction locations, magnetic fields and observer lines of sight. The jet luminosity and star properties are fixed, but the results are easily scalable under changes of these parameters. Individual jet-star interactions produce synchrotron and inverse Compton (IC) emission that peaks from X-rays to MeV energies (depending on the magnetic field), and at $\sim 100-1000$ GeV (depending on the stellar type), respectively. The radiation spectrum is hard in the scenarios explored here due to non-radiative cooling dominance, as low-energy electrons are efficiently advected even under relatively high magnetic fields. Interactions of jets with cold stars lead to an even harder IC spectrum because of the Klein-Nishina effect in the cross-section. Doppler boosting has a strong impact on the observer luminosity. The emission levels for individual interactions found here are in the line of previous, more approximate, estimates, strengthening the hypothesis that collective jet-star interactions could significantly contribute at high energies under efficient particle acceleration. △ Less

Submitted 15 November, 2016; v1 submitted 7 April, 2016; originally announced April 2016.

Comments: 13 pages, 2 tables, 14 figures. Language corrections mostly. Matches the published version in A&A 591, A15 (2016)

Journal ref: A&A 591, A15 (2016)

The origin of the X-ray-emitting object moving away from PSR B1259-63

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: A mysterious X-ray-emitting object has been detected moving away from the high-mass gamma-ray binary PSR B1259-63, which contains a non-accreting pulsar and a Be star whose winds collide forming a complex interaction structure. Given the strong eccentricity of this binary, the interaction structure should be strongly anisotropic, which together with the complex evolution of the shocked winds, coul… ▽ More A mysterious X-ray-emitting object has been detected moving away from the high-mass gamma-ray binary PSR B1259-63, which contains a non-accreting pulsar and a Be star whose winds collide forming a complex interaction structure. Given the strong eccentricity of this binary, the interaction structure should be strongly anisotropic, which together with the complex evolution of the shocked winds, could explain the origin of the observed moving X-ray feature. We propose here that a fast outflow made of a pulsar-stellar wind mixture is always present moving away from the binary in the apastron direction, with the injection of stellar wind occurring at orbital phases close to periastron passage. This outflow periodically loaded with stellar wind would move with a high speed, and likely host non-thermal activity due to shocks, on scales similar to those of the observed moving X-ray object. Such an outflow is thus a very good candidate to explain this X-ray feature. This, if confirmed, would imply pulsar-to-stellar wind thrust ratios of $\sim 0.1$, and the presence of a jet-like structure on the larger scales, up to its termination in the ISM. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 5 pages, 5 figures