Context. High-resolution Very Long Baseline Interferometry (VLBI) observations of active galactic nuclei revealed traveling and stationary or quasi-stationary radio components in several blazar jets. ...The traveling radio components are, in general, interpreted as shock waves generated by pressure perturbations injected at the jet nozzle. The stationary features can be interpreted as recollimation shocks in nonpressure matched jets if they show a quasi-symmetric bump in the spectral index distribution. In some jets there may be interactions between the two kinds of shocks. These shock-shock interactions are observable with VLBI techniques and their signature should also be imprinted on the single-dish light curves. Aims. In this paper, we investigate the spectral evolution produced by the interaction between a recollimation shock with traveling shock waves to address the question of whether these interactions contribute to the observed flares and what their signature in both single-dish and VLBI observations looks like. Methods. We performed relativistic hydrodynamic simulations of overpressured and pressure-matched jets. To simulate the shock interaction we injected a perturbation at the jet nozzle once a steady state was reached. We computed the nonthermal emission, including adiabatic and synchotron losses, resulting from the simulation. Results. We show that the injection of perturbations in a jet can produce a bump in emission at GHz frequencies previous to the main flare, which is produced when the perturbation fills the jet in the observer’s frame. The detailed analysis of our simulations and the nonthermal emission calculations show that interaction between a recollimation shock and traveling shock produce a typical and clear signature in both the single-dish light curves and in the VLBI observations: the flaring peaks are higher and delayed with respect to the evolution of a perturbation through a conical jet. This fact can allow us to detect such interactions for stationary components lying outside of the region in which the losses are dominated by inverse Compton scattering.
Context. The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary system suffers periodic and random variations of orbital and non-linear dynamical origins. ...The characterization of the evolution of the wind interaction region is necessary for understanding the rich phenomenology of these sources. Aims. For the first time, we simulate in 3 dimensions the interaction of isotropic stellar and relativistic pulsar winds along one full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolution and size, and of different state equations: a γ̂-constant ideal gas, and an ideal gas with γ̂ dependent on temperature. Methods. We used the code PLUTO to carry out relativistic hydrodynamical simulations in 2 and 3 dimensions of the interaction between a slow dense wind and a mildly relativistic wind with Lorentz factor 2, along one full orbit in a region up to ~100 times the binary size. The different 2-dimensional simulations were carried out with equal and larger grid resolution and size, and one was done with a more realistic equation of state than in 3 dimensions. Results. The simulations in 3 dimensions confirm previous results in 2 dimensions, showing: a strong shock induced by Coriolis forces that terminates the pulsar wind also in the opposite direction to the star; strong bending of the shocked-wind structure against the pulsar motion; and the generation of turbulence. The shocked flows are also subject to a faster development of instabilities in 3 dimensions, which enhances shocks, two-wind mixing, and large-scale disruption of the shocked structure. In 2 dimensions, higher resolution simulations confirm lower resolution results, simulations with larger grid sizes strengthen the case for the loss of the general coherence of the shocked structure, and simulations with two different equations of state yield very similar results. In addition to the Kelvin-Helmholtz instability, discussed in the past, we find that the Richtmyer-Meshkov and the Rayleigh-Taylor instabilities are very likely acting together in the shocked flow evolution. Conclusions. Simulations in 3 dimensions confirm that the interaction of stellar and pulsar winds yields structures that evolve non-linearly and become strongly entangled. The evolution is accompanied by strong kinetic energy dissipation, rapid changes in flow orientation and speed, and turbulent motion. The results of this work strengthen the case for the loss of the coherence of the whole shocked structure on large scales, although simulations of more realistic pulsar wind speeds are needed.
ABSTRACT This work presents the first characterization of the internal structure of overpressured, steady superfast-magnetosonic relativistic jets in connection with their dominant type of energy. To ...this aim, relativistic magnetohydrodynamic simulations of different jet models threaded by a helical magnetic field have been analyzed covering a wide region in the magnetosonic Mach number-specific internal energy plane. The merit of this plane is that models dominated by different types of energy (internal energy: hot jets; rest-mass energy: kinetically dominated jets; magnetic energy: Poynting-flux-dominated jets) occupy well-separated regions. The analyzed models also cover a wide range of magnetizations. Models dominated by the internal energy (i.e., hot models, or Poynting-flux-dominated jets with magnetizations larger than but close to one) have a rich internal structure characterized by a series of recollimation shocks and present the largest variations in the flow Lorentz factor (and internal energy density). Conversely, in kinetically dominated models, there is not much internal or magnetic energy to be converted into kinetic, and the jets are featureless with small variations in the flow Lorentz factor. The presence of a significant toroidal magnetic field threading the jet produces large gradients in the transversal profile of the internal energy density. Poynting-flux-dominated models with high magnetization ( 10 or larger) are prone to be unstable against magnetic pinch modes, which sets limits on the expected magnetization in parsec-scale active galactic nucleus jets or constrains their magnetic field configuration.
Context.
IceCube has reported a very-high-energy neutrino (IceCube-170922A) in a region containing the blazar
TXS 0506+056
. Correlated gamma-ray activity has led to the first high-probability ...association of a high-energy neutrino with an extragalactic source. This blazar has been found to be in a radio outburst during the neutrino event.
Aims.
Our goal is to probe the sub-milliarcsecond properties of the radio jet right after the neutrino detection and during the further evolution of the radio outburst.
Methods.
We performed target of opportunity observations at 43 GHz frequency using very long baseline interferometry imaging, corresponding to 7 mm in wavelength, with the Very Long Baseline Array two and eight months after the neutrino event.
Results.
We produced two images of the radio jet of TXS 0506+056 at 43 GHz with angular resolutions of (0.2 × 1.1) mas and (0.2 × 0.5) mas, respectively. The source shows a compact, high brightness temperature core, albeit not approaching the equipartition limit and a bright and originally very collimated inner jet. Beyond approximately 0.5 mas from the millimeter-VLBI core, the jet loses this tight collimation and expands rapidly. During the months after the neutrino event associated with this source, the overall flux density is rising. This flux density increase happens solely within the core. Notably, the core expands in size with apparent superluminal velocity during these six months so that the brightness temperature drops by a factor of three despite the strong flux density increase.
Conclusions.
The radio jet of TXS 0506+056 shows strong signs of deceleration and/or a spine-sheath structure within the inner 1 mas, corresponding to about 70–140 pc in deprojected distance, from the millimeter-VLBI core. This structure is consistent with theoretical models that attribute the neutrino and gamma-ray production in TXS 0506+056 to interactions of electrons and protons in the highly relativistic jet spine with external photons originating from a slower moving jet region. Proton loading due to jet-star interactions in the inner host galaxy is suggested as the possible cause of deceleration.
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.
Aims.
The collimation of relativistic jets in galaxies is a poorly understood process. Detailed radio studies of the jet collimation region have been performed so far in a few individual objects, ...providing important constraints for jet formation models. However, the extent of the collimation zone as well as the nature of the external medium possibly confining the jet are still debated.
Methods.
In this article, we present a multifrequency and multiscale analysis of the radio galaxy NGC 315, including the use of mm-VLBI data up to 86 GHz, aimed at revealing the evolution of the jet collimation profile. We then consider results from the literature to compare the jet expansion profile in a sample of 27 low-redshift sources, mainly comprising radio galaxies and BL Lacs, which were classified based on the accretion properties as low-excitation (LEG) and high-excitation (HEG) galaxies.
Results.
We propose that the jet collimation in NGC 315 is completed on sub-parsec scales. A transition from a parabolic to conical jet shape is detected at
z
t
= 0.58 ± 0.28 parsecs or ∼5 × 10
3
Schwarzschild radii (
R
S
) from the central engine, a distance which is much smaller than the Bondi radius,
r
B
∼ 92 pc, estimated based on X-ray data. The jet in this and in a few other LEG in our sample may be initially confined by a thick disk extending out to ∼10
3
− 10
4
R
S
. A comparison between the mass-scaled jet expansion profiles of all sources indicates that jets in HEG are surrounded by thicker disk-launched sheaths and collimate on larger scales with respect to jets in LEG. These results suggest that disk winds play an important role in the jet collimation mechanism, particularly in high-luminosity sources. The impact of winds on the origin of the FRI and FRII dichotomy in radio galaxies is also discussed.
Supermassive black holes (SMBH) are essential for the production of jets in radio-loud active galactic nuclei (AGN). Theoretical models based on (Blandford & Znajek 1977, MNRAS, 179, 433) extract the ...rotational energy from a Kerr black hole, which could be the case for NGC 1052, to launch these jets. This requires magnetic fields on the order of 103G to 104G. We imaged the vicinity of the SMBH of the AGN NGC 1052 with the Global Millimetre VLBI Array and found a bright and compact central feature that is smaller than 1.9 light days (100 Schwarzschild radii) in radius. Interpreting this as a blend of the unresolved jet bases, we derive the magnetic field at 1 Schwarzschild radius to lie between 200 G and ~ 8.3 × 104 G consistent with Blandford & Znajek models.
Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or subrelativistic speeds on kiloparsec scales. Several mechanisms have been proposed ...to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code ratpenat. We study the evolution of Fanaroff–Riley Class I (FR I) jets, with kinetic powers L
j ∼ 1041–1044 erg s−1, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FR I jets within the first kiloparsec. However, powerful FR I jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kiloparsec scales is initiated by mass loading in the weaker FR I jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.
Context.
With the increase in the sensitivity and resolution of radio interferometry within recent years, the study of the collimation and acceleration region of extragalactic jets in active galactic ...nuclei (AGN) has come into focus. Whereas a large fraction of AGN jets show a change from parabolic to conical collimation profile around the Bondi radius, there is a small number of sources that display a deviation from this standard picture, including the radio galaxy NGC 1052.
Aims.
We study the jet width profile, which provides valuable information about the interplay between the central engine and accretion disk system and the formation, acceleration, and collimation of the jets.
Methods.
We observed the double-sided, low-radio-power active galaxy NGC 1052 at six frequencies with the VLBA in 2017 and at 22 GHz with RadioAstron in 2016. These data are combined with archival 15, 22, and 43 GHz multi-epoch VLBA observations. From ridge-line fitting we obtained width measurements along the jet and counter-jet which were fitted with single and broken power laws.
Results.
We find a clear break point in the jet collimation profile at ∼10
4
R
S
(Schwarzschild radii). Downstream of the break, the collimation is conical with a power-law index of 1.0 − 1.2 (cylindrical 0; parabolic 0.5; conical 1) for both jets. On the other hand, the upstream power-law index of 0.36 for the approaching jet is neither cylindrical nor parabolic, and the value of 0.16 for the receding jet suggests this latter is close-to cylindrical. For both jets we find a large opening angle of ∼30° at a distance of ∼10
3
R
S
and well-collimated structures with an opening angle of < 10° downstream of the break.
Conclusions.
There are significant differences between the upstream collimation profiles of the approaching (eastern) and receding (western) jets. Absorption or scattering in the surrounding torus as well as an accretion wind may mimic a cylindrical profile. We need to increase the observing frequencies, which do not suffer from absorption, in order to find the true jet collimation profile upstream of 10
4
R
S
.
Context. Extragalactic jets are formed close to supermassive black-holes in the center of galaxies. Large amounts of gas, dust, and stars cluster in the galaxy nucleus, and interactions between this ...ambient material and the jet base should be frequent, having dynamical as well as radiative consequences. Aims. This work studies the dynamical interaction of an obstacle, a clump of matter or the atmosphere of an evolved star, with the innermost region of an extragalactic jet. Jet mass-loading and the high-energy outcome of this interaction are briefly discussed. Methods. Relativistic hydrodynamical simulations with axial symmetry have been carried out for homogeneous and inhomogeneous obstacles inside a relativistic jet. These obstacles may represent a medium inhomogeneity or the disrupted atmosphere of a red giant star. Results. Once inside the jet, an homogeneous obstacle expands and gets disrupted after few dynamical timescales, whereas in the inhomogeneous case, a solid core can smoothen the process, with the obstacle mass-loss dominated by a dense and narrow tail pointing in the direction of the jet. In either case, matter is expected to accelerate and eventually get incorporated to the jet. Particles can be accelerated in the interaction region, and produce variable gamma-rays in the ambient matter, magnetic and photon fields. Conclusions. The presence of matter clumps or red giants into the base of an extragalactic jet likely implies significant jet mass-loading and slowing down. Fast flare-like gamma-ray events, and some level of persistent emission, are expected due to these interactions.