Search

[Context] γ 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, ~5 from its border. [Aims] We want to identify the radio sources within the region encompassing γ 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. [Methods] 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. [Results] 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. [Conclusions] We have identified the radio sources overlapping γ 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., P.B. acknowledges support from ANPCyT PICT 0773-2017. C.H.I.C. acknowledges the support of the Department of Atomic Energy, Government of India, under the project 12-R&D-TFR-5.02-0700. J.M.P. and V.B.R. acknowledge financial support from the Spanish Agencia Estatal de Investigación through the grant PID2019-105510GB-C31/ AEI / 10.13039/501100011033 and ‘Unit of Excellence María de Maeztu 20202023’ award to the Institute of Cosmos Sciences (CEX2019-000918-M).

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

Context LS 5039 is a powerful high-mass gamma-ray binary that probably hosts a young non-accreting pulsar. However, despite the wealth of data available, the means by which the non-thermal emitter is powered are still unknown. Aims. We use a dynamical-radiative numerical model, and multiwavelength data, to constrain the properties of a hypothetical pulsar wind that would power the non-thermal emitter in LS 5039. Methods. We ran simulations of an ultrarelativistic (weakly magnetized) cold e±-wind that Compton scatters stellar photons and that dynamically interacts with the stellar wind. The effects of energy losses on the unshocked e±-wind dynamics, and the geometry of the two-wind contact discontinuity, are computed for different wind models. The predicted unshocked e±-wind radiation at periastron, when expected to be the highest, is compared to LS 5039 data. Results. The minimum possible radiation from an isotropic cold e±-wind overpredicts the X-ray to gamma-ray fluxes at periastron by a factor of ∼3. In the anisotropic (axisymmetric) wind case X-ray and ≳ 100 MeV data are not violated by wind radiation if the wind axis is at ≲ 20-40° from the line of sight (chance probability of ≲ 6-24%), depending on the anisotropic wind model, or if the wind Lorentz factor ∈ 102 - 103, in which case the wind power can be higher, but it requires e±-multiplicities of ∼106 and 109 for a 10-2 s and 10 s pulsar period, respectively. Conclusions. The studied model predicts that a weakly magnetized cold pulsar e±-wind in LS 5039 should be strongly anisotropic, with either a wind Lorentz factor ∈ 102 - 103 and very high multiplicities or with a fine-tuned wind orientation. A weakly magnetized, cold baryon-dominated wind would be a possible alternative, but then the multiplicities should be rather low, while the baryon-to-e± energy transfer should be very efficient at wind termination. A strongly magnetized cold wind seems to be the most favorable option as it is consistent with recent research on pulsar winds and does not require fine-tuning of the pulsar wind orientation, and the wind multiplicity and Lorentz factor are less constrained., Acknowledgements. We thank the anonymous referee for constructive and useful comments that helped to improve the manuscript. We are grateful to Dmitry Khangulyan for insightful comments on this work. V.B-R. acknowledges financial support by the Spanish Ministerio de Economía, Industria y Competitivi-dad (MINEICO/FEDER, UE) under grant AYA2016-76012-C3-1-P, from the State Agency for Research of the Spanish Ministry of Science and Innovation under grant PID2019-105510GB-C31 and through the ”Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M), and by the Catalan DEC grant 2017 SGR 643. V.B-R. is Correspondent Researcher of CONICET, Argentina, at the IAR.

[Context]: Gamma-ray binaries are systems that emit nonthermal radiation peaking at energies above 1 MeV. One proposed scenario to explain their emission consists of a pulsar orbiting a massive star, with particle acceleration taking place in shocks produced by the interaction of the stellar and pulsar winds. [Aims]: We develop a semi-Analytical model of the nonthermal emission of the colliding-wind structure, which includes the dynamical effects of orbital motion. We apply the model to a general case and to LS 5039. [Methods]: The model consists of a one-dimensional emitter, the geometry of which is affected by Coriolis forces owing to orbital motion. Two particle accelerators are considered: one at the two-wind standoff location and the other one at the turnover produced by the Coriolis force. Synchrotron and inverse Compton emission is studied taking into account Doppler boosting and absorption processes associated to the massive star. [Results]: If both accelerators are provided with the same energy budget, most of the radiation comes from the region of the Coriolis turnover and beyond, up to a few orbital separations from the binary system. Significant orbital changes of the nonthermal emission are predicted in all energy bands. The model allows us to reproduce some of the LS 5039 emission features, but not all of them. In particular, the MeV radiation is probably too high to be explained by our model alone, the GeV flux is recovered but not its modulation, and the radio emission beyond the Coriolis turnover is too low. The predicted system inclination is consistent with the presence of a pulsar in the binary. [Conclusions]: The model is quite successful in reproducing the overall nonthermal behavior of LS 5039. Some improvements are suggested to better explain the phenomenology observed in this source, such as accounting for particle reacceleration beyond the Coriolis turnover, unshocked pulsar wind emission, and the three-dimensional extension of the emitter., We acknowledge support by the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grant AYA2016-76012-C3-1-P, with partial support by the European Regional Development Fund (ERDF/FEDER), and from the Catalan DEC grant 2017 SGR 643. EM acknowledges support from MINECO through grant BES-2016-076342.

The persistent gamma-ray source 4FGL J0647.7−4418 is tentatively associated in the latest Fermi catalogue with the subdwarf O-type X-ray binary HD 49798. However, an active galactic nucleus (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 per cent confidence ellipse, is a more conceivable counterpart candidate to the Fermi source due to its strong blazar similarities., This work was supported by the Agencia Estatal de Investigación grants AYA2016-76012-C3-1-P and AYA2016-76012-C3-3-P from the Spanish Ministerio de Economía y Competitividad (MINECO), by Consejería de Economía, Innovación, Ciencia y Empleo, Junta de Andalucía under research group FQM-322, by grant MDM-2014-0369 of the Institute of Cosmos Sciences (ICCUB, Unidad de Excelencia ‘María de Maeztu’), and by the Catalan DEC grant 2017 SGR 643, as well as Fondo Europeo de Desarrollo Regional (FEDER) funds. RC acknowledges support from National Aeronautics and Space Administration (NASA) Fermi grant NNX15AU83G.

[Context] Dark matter may consist, at least partially, of primordial black holes formed during the radiation-dominated era. The radiation produced by accretion onto primordial black holes leaves characteristic signatures on the properties of the medium at high redshift, before and after hydrogen recombination. Therefore, reliable modeling of accretion onto these objects is required to obtain robust constraints on their abundance., [Aims] We investigate the effect of mechanical feedback, that is, the impact of outflows (winds and– or –jets) on the medium, on primordial black hole accretion, and thereby on the associated radiation., [Methods] Using analytical and numerical calculations, we studied for the first time the possibility that outflows can reduce the accretion rate of primordial black holes with masses similar to those detected by the LIGO-Virgo collaboration., [Results] Despite the complexity of the accretion rate evolution, mechanical feedback is useful in to significantly reducing the primordial black hole accretion rate, at least by one order of magnitude, when outflows are aligned with the motion of the compact object. If the outflow is perpendicular to the direction of motion, the effect is less important, but it is still non-negligible., [Conclusions] Outflows from primordial black holes, even rather weak ones, can significantly decrease the accretion rate, effectively weakening abundance constraints on these objects. Our results motivate further numerical simulations with a more realistic setup, which would yield more precise quantitative predictions., V.B-R. acknowledges support by the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grant AYA2016-76012-C3-1-P, with partial support by the European Regional Development Fund (ERDF/FEDER), MDM-2014-0369 of ICCUB (Unidad de Excelencia “María de Maeztu”), and the Catalan DEC grant 2017 SGR 643. NB is supported by the Spanish MINECO under grant BES-2015-073372.

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

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., 13 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

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

The stellar wind in high-mass microquasars should interact with the jet. This interaction, coupled with orbital motion, is expected to make the jet follow a helical, nonballistic trajectory. The jet energy dissipated by this interaction, through shocks for example, could lead to nonthermal activity on scales significantly larger than the system size. We calculate the broadband emission from a jet affected by the impact of the stellar wind and orbital motion in a high-mass microquasar. We employ a prescription for the helical trajectory of a jet in a system with a circular orbit. Subsequently, assuming electron acceleration at the onset of the helical jet region, we compute the spatial and energy distribution of these electrons, and their synchrotron and inverse Compton emission including gamma-ray absorption effects. For typical source parameters, significant radio, X- and gamma-ray luminosities are predicted. The scales on which the emission is produced may reduce, but not erase, orbital variability of the inverse Compton emission. The wind and orbital effects on the radio emission morphology could be studied using very long baseline interferometric techniques. We predict significant broadband emission, modulated by orbital motion, from a helical jet in a high-mass microquasar. This emission may be hard to disentangle from radiation of the binary itself, although the light curve features, extended radio emission, and a moderate opacity to very high-energy gamma rays, could help to identify the contribution from an extended (helical) jet region., Comment: 11 pages, 11 figures, accepted for publication in Astronomy & Astrophysics. [v2,v3]: typos corrected. [v4]: figures 2 and 3 corrected

[Context]: Cosmic reionization put an end to the dark ages that came after the recombination era. Observations seem to favor the scenario of massive-star photons generated in low-mass galaxies being responsible for the bulk of reionization and, whereas a possible contribution from AGN accretion disks has been widely considered, they are currently thought to have had a minor role in reionization. [Aims]: We aim to study the possibility of AGN having contributed to reionization not only through their accretion disks, but also through ionizing photons coming from the AGN jets interacting with the intergalactic medium. Methods. We adopt an empirically derived AGN luminosity function at z ≃ 6, 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 ≃ 15-5. [Results]: We show that the contribution of the AGN jet lobes to the reionization of the Universe at z ∼ 6 might have been as high as ? 10% of that of star-forming galaxies, under the most favorable conditions of jetted and obscuration fraction. [Conclusions]: 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., We acknowledge support by the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grant AYA2016-76012-C3-1-P, with partial support by the European Regional Development Fund (ERDF/FEDER), MDM-2014-0369 of ICCUB (Unidad de Excelencia “María de Maeztu”), and the Catalan DEC Grant 2017 SGR 643. N. T.-A. acknowledges support from MINECO through FPU14/04887 Grant.

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

Aims. We probe the high-energy (>60 MeV) emission from the black hole X-ray binary system, Cygnus X-1, and investigate its origin.Methods. We analyzed 7.5 yr of data by Fermi -LAT with the latest Pass 8 software version.Results. We report the detection of a signal at ~8σ statistical significance that is spatially coincident with Cygnus X-1 and has a luminosity of 5.5 × 1033 erg s-1 , above 60 MeV. 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 ~20 GeV without any sign of spectral break, is well fit by a power-law function with a photon index of 2.3 ± 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 indicates the anisotropic inverse-Compton on stellar photons as the mechanism at work, thus constraining the emission region to a distance 1011 –1013 cm from the black hole.

The active galactic nucleus PG 1553+113 was observed by the MAGIC telescope in July 2006 during a multiwavelength campaign, in which telescopes in the optical, X-ray, and very high energies participated. Although the MAGIC data were affected by strong atmospheric absorption (calima), they were analyzed after applying a correction. In 8.5 hours, a signal was detected with a significance of 5.0 sigma. The integral flux above 150 GeV was (2.6 +/- 0.9)*10^{-7} ph/s/m^2. By fitting the differential energy spectrum with a power law, a spectral index of -4.1 +/- 0.3 was obtained., 4 pages, 3 figures, accepted for publication in A&A

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 γ-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 γ-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., Instituto Argentino de Radioastronomía, Facultad de Ciencias Astronómicas y Geofísicas

LS 5039 is an X-ray binary that presents non-thermal radio emission. The radiation at $\sim 5$ GHz is quite steady and optically thin, consisting on a dominant core plus an extended jet-like structure. There is a spectral turnover around 1 GHz, and evidence of variability at timescales of 1 yr at 234 MHz. We investigate the radio emitter properties using the available broadband radio data, and assuming two possible scenarios to explain the turnover: free-free absorption in the stellar wind, or synchrotron self-absorption. We use the relationships between the turnover frequency, the stellar wind density, the emitter location, size and magnetic field, and the Lorentz factor of the emitting electrons, as well as a reasonable assumption on the energy budget, to infer the properties of the low-frequency radio emitter. Also, we put this information in context with the broadband radio data. The location and size of the low-frequency radio emitter can be restricted to $\ga$ few AU from the primary star, its magnetic field to $\sim 3\times 10^{-3}-1$ G, and the electron Lorentz factors to $\sim 10-100$. The observed variability of the extended structures seen with VLBA would point to electron bulk velocities $\ga 3\times 10^8$ cm s$^{-1}$, whereas much less variable radiation at 5 GHz would indicate velocities for the VLBA core $\la 10^8$ cm s$^{-1}$. The emission at 234 MHz in the high state would mostly come from a region larger than the dominant broadband radio emitter. We suggest a scenario in which secondary pairs, created via gamma-ray absorption and moving in the stellar wind, are behind the steady broadband radio core, whereas the resolved jet-like radio emission would come from a collimated, faster, outflow., accepted for publication in A&A, 5 pages, 2 figures, 1 table

Strong interactions between jets and stellar winds at binary system spatial scales could occur in high-mass microquasars. We study here, mainly from a dynamical but also a radiative point of view, the collision between a dense stellar wind and a mildly relativistic hydrodynamical jet of supersonic nature.}{We have performed numerical 2-dimensional simulations of jets, with cylindrical and planar (slab) symmetry, crossing the stellar wind material. From the results of the simulations, we derive estimates of the particle acceleration efficiency, using first order Fermi acceleration theory, and give some insight on the possible radiative outcomes. We find that, during jet launching, the jet head generates a strong shock in the wind. During and after this process, strong recollimation shocks can occur due to the initial overpressure of the jet with its environment. The conditions in all these shocks are convenient to accelerate particles up to $\sim$ TeV energies, which can lead to leptonic (synchrotron and inverse Compton) and hadronic (proton-proton) radiation. In principle, the cylindrical jet simulations show that the jet is stable, and can escape from the system even for relatively low power. However, when accounting for the wind ram pressure, the jet can be bent and disrupted for power $\la 10^{36}$ erg s$^{-1}$., Accepted for publication in Astronomy & Astrophysics