Since Fanaroff & Riley (1974) reported the morphological and brightness dichotomy of radiogalaxies, and it became clear that the symmetric emission from jets and counter-jets in the centre-brightened, less powerful, FRI sources could be caused by jet deceleration, many works have addressed different mechanisms that could cause this difference. Recent observational results seem to indicate that the deceleration must be caused by the development of small-scale instabilities that force mixing at the jet boundary. According to these results, the mixing layer expands and propagates down to the jet axis along several kiloparsecs, until it covers the whole jet cross-section. Several candidate mechanisms have been proposed as the initial trigger for the generation of such mixing layer. However, the instabilities proposed so far do not fully manage to explain the observations of FRI jets and/or require a triggering mechanism. Therefore, there is not still a satisfactory explanation for the original cause of jet deceleration. In this letter, I show that the penetration (and exit) of stars from jets could give the adequate explanation by means of creating a jet-interstellar medium mixing layer that expands across the jet.

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. We have performed target-of-opportunity VLBI imaging observations at 43 GHz frequency with the VLBA two and eight months, respectively, after the neutrino event. We produced two images of TXS 0506+056 with angular resolutions of (0.2x1.1) mas and (0.2x0.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 about 0.5 mas from the mm-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. The core expands in size with apparent superluminal velocity during these six months so that the brightness temperature drops by a factor of three in spite of the strong flux density increase. 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 pc to 140 pc in deprojected distance) from the mm-VLBI core. This structure is consistent with theoretical models that attribute the neutrino and {\gamma}-ray production 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

Detailed studies of relativistic jets in active galactic nuclei (AGN) require high-fidelity imaging at the highest possible resolution. This can be achieved using very long baseline interferometry (VLBI) at radio frequencies, combining worldwide (global) VLBI arrays of radio telescopes with a space-borne antenna on board a satellite. We present multiwavelength images made of the radio emission in the powerful quasar S5 0836+710, obtained using a global VLBI array and the antenna Spektr-R of the RadioAstron mission of the Russian Space Agency, with the goal of studying the internal structure and physics of the relativistic jet in this object. The RadioAstron observations at wavelengths of 18cm, 6cm, and 1.3cm are part of the Key Science Program for imaging radio emission in strong AGN. The internal structure of the jet is studied by analyzing transverse intensity profiles and modeling the structural patterns developing in the flow. The RadioAstron images reveal a wealth of structural detail in the jet of S5 0836+710 on angular scales ranging from 0.02mas to 200mas. Brightness temperatures in excess of $10^{13}$\,K are measured in the jet, requiring Doppler factors of $\ge 100$ for reconciling them with the inverse Compton limit. Several oscillatory patterns are identified in the ridge line of the jet and can be explained in terms of the Kelvin-Helmholtz (KH) instability. The oscillatory patterns are interpreted as the surface and body wavelengths of the helical mode of the KH instability. The interpretation provides estimates of the jet Mach number and of the ratio of the jet to the ambient density, which are found to be $M_\mathrm{j}\approx 12$ and $\eta\approx 0.33$. The ratio of the jet to the ambient density should be conservatively considered an upper limit because its estimate relies on approximations.

A simple look at the steady high-energy Universe reveals a clear correlation with outflows generated around compact objects (winds and jets). In the case of relativistic jets, they are thought to be produced as a consequence of the extraction of rotational energy from a Kerr black hole (Blandford-Znajek), or from the disc (Blandford-Payne). A fraction of the large energy budget provided by accretion and/or black hole rotational energy is invested into jet formation. After formation, the acceleration and collimation of these outflows allow them to propagate to large distances away from the compact object. The synchrotron cooling times demand that re-acceleration of particles takes place along the jets to explain high-energy and very-high-energy emission from kiloparsec scales. At these scales, jets in radio galaxies are divided in two main morphological/luminosity types, namely, Fanaroff-Riley type I and II (FRI, FRII), the latter being more luminous, collimated and edge-brightened than the former, which show clear hints of decollimation and deceleration. In this contribution, I summarise a set of mechanisms that may contribute to dissipate magnetic and kinetic energy: Magnetohydrodynamic instabilities or jet-obstacle interactions trigger shocks, shearing and mixing, which are plausible scenarios for particle acceleration. I also derive an expression for the expected distance in which the entrainment by stellar winds starts to be relevant, which is applicable to FRI jets. Finally, I discuss the differences in the evolutionary scenarios and the main dissipative mechanisms that take place in extragalactic and microquasar jets.

Our goal is to study the termination of an AGN jet in the young universe and to deduce physical parameters of the jet and the intergalactic medium. We use LOFAR to image the long-wavelength radio emission of the high-redshift blazar S5 0836+710 on arcsecond scales between 120 MHz and 160 MHz. The LOFAR image shows a compact unresolved core and a resolved emission region about 1.5 arcsec to the southwest of the radio core. This structure is in general agreement with previous higher-frequency radio observations with MERLIN and the VLA. The southern component shows a moderately steep spectrum with a spectral index of about $\gtrsim -1$ while the spectral index of the core is flat to slightly inverted. In addition, we detect for the first time a resolved steep-spectrum halo with a spectral index of about $-1$ surrounding the core. The arcsecond-scale radio structure of S5 0836+710 can be understood as an FR II-like radio galaxy observed at a small viewing angle. The southern component can be interpreted as the region of the approaching jet's terminal hotspot and the halo-like diffuse component near the core can be interpreted as the counter-hotspot region. From the differential Doppler boosting of both features, we can derive the hotspot advance speed to $(0.01-0.036)$ c. At a constant advance speed, the derived age of the source would exceed the total lifetime of such a powerful FR II-like radio galaxy substantially. Thus, the hotspot advance speed must have been higher in the past in agreement with a scenario in which the originally highly relativistic jet has lost collimation due to the growth of instabilities and has transformed into an only mildly relativistic flow. Our data suggest that the density of the intergalactic medium around this distant ($z=2.22$) AGN could be substantially higher than the values typically found in less distant FR II radio galaxies.

Particle acceleration in relativistic jets to very high energies occurs at the expense of the dissipation of magnetic or kinetic energy. Therefore, understanding the processes that can trigger this dissipation is key to the characterization of the energy budgets and particle acceleration mechanisms at action in active galaxies. Instabilities and entrainment are two obvious candidates to trigger dissipation. On the one hand, supersonic, relativistic flows threaded by helical fields, as expected from the standard formation models of jets in supermassive black-holes, are unstable to a series of magnetohydrodynamical instabilities, such as the Kelvin-Helmholtz, current-driven, or possibly the pressure-driven instabilities. Furthermore, in the case of expanding jets, the Rayleigh-Taylor and centrifugal instabilities may also develop. With all these destabilizing processes at action, a natural question is how can some jets keep their collimated structure along hundreds of kiloparsecs. On the other hand, the interaction of the jet with stars and clouds of gas that cross the flow in their orbits around the galactic centers provides another scenario in which kinetic energy can be efficiently converted into internal energy and particles can be accelerated to non-thermal energies. In this contribution, I review the conditions under which these processes occur and their role both in jet evolution and propagation and energy dissipation.

A number of extragalactic jets show periodic structures at different scales that can be associated with growing instabilities. The wavelengths of the developing instability modes and their ratios depend on the flow parameters, so the study of those structures can shed light on jet physics at the scales involved. In this work, we use the fits to the jet ridgeline obtained from different observations of S5 B0836$+$710 and apply stability analysis of relativistic, sheared flows to derive an estimate of the physical parameters of the jet. Based on the assumption that the observed structures are generated by growing Kelvin-Helmholtz (KH) instability modes, we have run numerical calculations of stability of a relativistic, sheared jet over a range of different jet parameters. We have spanned several orders of magnitude in jet-to-ambient medium density ratio, and jet internal energy, and checked different values of the Lorentz factor and shear layer width. This represents an independent method to obtain estimates of the physical parameters of a jet. By comparing the fastest growing wavelengths of each relevant mode given by the calculations with the observed wavelengths reported in the literature, we have derived independent estimates of the jet Lorentz factor, specific internal energy, jet-to-ambient medium density ratio and Mach number. We obtain a jet Lorentz factor $\gamma \simeq 12$, specific internal energy of $\varepsilon \simeq 10^{-2}\,c^2$, jet-to-ambient medium density ratio of $\eta\approx 10^{-3}$, and an internal (classical) jet Mach number of $M_\mathrm{j}\approx 12$. We also find that the wavelength ratios are better recovered by a transversal structure with a width of $\simeq 10\,\%$ of the jet radius. This method represents a powerful tool to derive the jet parameters in all jets showing helical patterns with different wavelengths.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

There is compelling evidence showing that extragalactic jets are a crucial ingredient in the evolution of host galaxies and their environments. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics at sub-parsec and parsec scales. They generate strong shocks in the ambient medium, associated with observed hotspots in FRII radio galaxies, and carve cavities that are filled with the shocked jet flow, dragging a large fraction of the interstellar gas along, in the form of slow, massive outflows within the host galaxies. In this paper, I discuss relevant processes associated to jet evolution in the frame of FRI-FRII dichotomy. In particular, I focus on the role of 1) the interaction between galactic atmospheres and the jet head on global FRII jet kinematics, and 2) mass load by stellar winds or small-scale instabilities on jet deceleration in FRI jets. The results presented are based on 3D relativistic hydrodynamical (RHD) and/or 2D axisymmetric, time-independent relativistic magnetohydrodynamical (RMHD) simulations.

Context: Few Active Galactic Nuclei (AGN) reveal a double sided jet system. However, these are crucial to understand basic physical properties of extragalactic jets. Aims: We address the questions whether jets in AGN are symmetric in nature, how well they are collimated on small scales, and how do they evolve with time. Methods: We monitored the sub-parsec scale morphology of NGC1052 with the VLBA at 43GHz from 2005 to 2009. Results: A detailed study of 29 epochs show an asymmetry between both jets. A kinematic analysis of the outflows reveals higher apparent velocities for the eastern jet as compared to the western one: beta_ej = 0.529+-0.038 and beta_wj = 0.343+-0.037, respectively. Contradictory to previous studies, we find higher flux densities for the western jet. The distribution of brightness temperature and jet width features well collimated jets up to 1mas distance to the dynamic center and a nearly conical outflow further outwards. Combining flux density ratios and velocities of the jet flow, we were not able to find a combination of intrinsic velocity and inclination angle of the jets that is consistent for all observations, which is contradictory to symmetrically evolving jets. Spectral index maps between quasi-simultaneous 22GHz and 43GHz observations support the existence of an optically thick absorber covering the innermost 1.6 mas around the 43GHz central feature and an optically thin jet emission with a spectral index of <= -1. Conclusions: Our results fit into a picture in which we expect larger internal energy and/or magnetic flux in the western jet and higher kinetic energy in the eastern jet. Previous observations at lower frequencies found slower velocities of the moving jet features as compared to this work. Considering the different velocities at different areas we suggest a spine-sheath structure with a faster inner layer and a slower outer layer.

We present a long-term numerical three-dimensional simulation of a relativistic outflow designed to be compared with previous results from axisymmetric, two-dimensional simulations, with existing analytical models and state-of-art observations. We follow the jet evolution from 1~kpc to 200~kpc, using a relativistic gas equation of state and a galactic profile for the ambient medium. We also show results from smaller scale simulations aimed to test convergence and different three-dimensional effects. We conclude that jet propagation can be faster than expected from axisymmetric simulations, covering tens of kiloparsecs in a few million years, until the dentist drill effect produced by the growth of helical instabilities slows down the propagation speed of the jet head. A comparison of key physical parameters of the jet structure as obtained from the simulations with values derived from observations of FRII sources reveals good agreement. Our simulations show that shock heating can play a significant role in the feedback from active galaxies, confirming previous 2D results. A proper description of galactic jets as a relativistic scenario, both dynamical and thermodynamical, reveals an extremely fast and efficient feedback process reheating the ICM, and therefore, with dramatic consequences on the galactic evolution. Our results point towards FRII jets as the source of the energetic electrons observed in radio relics.

We present relativistic magnetohydrodynamic (RMHD) simulations of stationary overpressured magnetized relativistic jets which are characterized by their dominant type of energy, namely internal, kinetic, or magnetic. Each model is threaded by a helical magnetic field with a pitch angle of $45^\circ$ and features a series of recollimation shocks produced by the initial pressure mismatch, whose strength and number varies as a function of the dominant type of energy. We perform a study of the polarization signatures from these models by integrating the radiative transfer equations for synchrotron radiation using as inputs the RMHD solutions. These simulations show a top-down emission asymmetry produced by the helical magnetic field and a progressive confinement of the emission into a jet spine as the magnetization increases and the internal energy of the non-thermal population is considered to be a constant fraction of the thermal one. Bright stationary components associated with the recollimation shocks appear presenting a relative intensity modulated by the Doppler boosting ratio between the pre-shock and post-shock states. Small viewing angles show a roughly bimodal distribution in the polarization angle due to the helical structure of the magnetic field, which is also responsible for the highly stratified degree of linear polarization across the jet width. In addition, small variations of the order of $26^\circ$ are observed in the polarization angle of the stationary components, which can be used to identify recollimation shocks in astrophysical jets.

Radio emission from protostellar jets is usually dominated by free-free emission from thermal electrons. However, in some cases, it has been proposed that non-thermal emission could also be present. This additional contribution from non-thermal emission has been inferred through negative spectral indices at centimeter wavelengths in some regions of the radio jets. In the case of HH 80-81, one of the most powerful protostellar jets known, linearly polarized emission has also been detected, revealing that the non-thermal emission is of synchrotron nature from a population of relativistic particles in the jet. This result implies that an acceleration mechanism should be taking place in some parts of the jet. Here, we present new high sensitivity and high angular resolution radio observations at several wavelengths (in the 3-20 cm range) of the HH80-81 radio jet. These new observations represent an improvement in sensitivity and angular resolution by a factor of $\sim$10 with respect to previous observations. This allows us to resolve the morphology of the radio jet, and to study the different emission mechanisms involved through spectral index maps. We conclude that synchrotron emission in this jet arises from an extended component detected at low frequencies and from the termination points of the jet, where strong shocks against the ambient medium can produce efficient particle acceleration.

We aim to better understand the dynamics within relativistic magneto-hydrodynamical flows in the extreme environment and close vicinity of supermassive black holes. To do so, we analyze the peculiar radio galaxy 3C 111, for which long-term polarimetric observations are available. We make use of the high spatial resolution of the VLBA network and the MOJAVE monitoring program, which provides high data quality also for single sources and allows us to study jet dynamics on parsec scales in full polarization with an evenly sampled time-domain. We additionally consider data from the IRAM 30-m Telescope as well as the SMA. Jet properties such as the electric vectors, the (polarized) flux density, feature size, and brightness temperature, describe a complex evolution of the polarized jet. The electric vector position angles (EVPAs) of features traveling down the jet perform a large and smooth rotation of $\gtrsim 180^{\circ}$ across a distance of about 20 pc. In contrast, the EVPAs are strongly variable within the first parsecs of the jet. We find a tendency towards transverse EVPAs across the jet with a local anomaly of aligned vectors in between. The transverse extent of the flow decreases coincident with a jump in brightness temperature around where we observe the EVPAs to turn into alignment with the jet flow. Also the gradients of the feature size and particle density with distance steepen in that region. We interpret the propagating polarized features with shocks and the observed local anomalies with the interaction of these shocks with a recollimation shock of the underlying flow. Together with a sheared magnetic field, this shock-shock interaction can explain the large rotation of the EVPA. The superimposed variability of the EVPAs close to the core is likely related to a clumpy Faraday screen, which also contributes significantly to the observed EVPA rotation in that region.

Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydrodynamics. As a first step, we have studied the structure of steady-state configurations of jets by using numerical simulations. We have used a helical field configuration and have changed different relevant parameters that control the way in which the energy flux is distributed in jets (namely, the proportion of the energy flux carried by internal, kinetic or magnetic energy). Our results show significant differences among the different kinds of jets. Finally, we also report on results based on synthetic maps of our simulated jets.

We address the linear stability of a discontinuous surface of a relativistic flow in the context of a jet that oscillates radially as it propagates. The restoring force of the oscillation is expected to drive a Rayleigh-Taylor instability (RTI) at the interface between the jet and its cocoon. We perform a linear analysis and numerical simulations of the growth of the RTI in the transverse plane to the jet flow with a uniform acceleration. In this system, an inertia force due to the uniform acceleration acts as the restoring force for the oscillation. We find that not only the difference in the inertia between the two fluids separated by the interface but also the pressure at the interface helps to drive the RTI because of a difference in the Lorenz factor across the discontinuous surface of the jet. The dispersion relation indicates that the linear growth rate of each mode becomes maximum when the Lorentz factor of the jet is much larger than that of the cocoon and the pressure at the jet interface is relativistic. By comparing the linear growth rates of the RTI in the analytical model and the numerical simulations, the validity of our analytically derived dispersion relation for the relativistic RTI is confirmed.

Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this Letter, we compare the pressure evolution in the hotspot and the cavity generated by relativistic and classical jets. Our results show that the classical approach underestimates the cavity pressure by a factor larger or equal to 2 for a given shocked volume during the whole active phase. The tension between both approaches can only be alleviated by unrealistic jet flow densities or gigantic jet areas in the classical case. As a consequence, the efficiency of a relativistic jet heating the ambient is typically around 20% larger compared with a classical jet, and the heated volume is 2 to 10 times larger during the time evolution. This conflict translates into two substantially disparate manners, both spatially and temporal, of heating the ambient medium. These differences are expected to have relevant implications on the star formation rates of the host galaxies and their evolution.

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.

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.

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.

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.

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 1) 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 nor magnetic energy to be converted into kinetic one 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) are prone to be unstable against magnetic pinch modes, which sets limits to the expected magnetization in parsec-scale AGN jets {and/or constrains their magnetic field configuration}.

Supermassive black holes (SMBH) are essential for the production of jets in radio-loud active galactic nuclei (AGN). Theoretical models based on Blandford & Znajek extract the rotational energy from a Kerr black hole, which could be the case for NGC1052, to launch these jets. This requires magnetic fields of the order of $10^3\,$G to $10^4\,$G. We imaged the vicinity of the SMBH of the AGN NGC1052 with the Global Millimetre VLBI Array and found a bright and compact central feature, 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 ~80000 G consistent with Blandford & Znajek models.

High resolution Very Long Baseline Interferometry (VLBI) observations of Active Galactic Nuclei (AGN) revealed traveling and stationary or quasi-stationary radio-components in several blazar jets. The traveling ones 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 non-pressure 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. We performed relativistic hydrodynamic (RHD) simulations of over-pressured 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 non-thermal emission (including adiabatic and synchotron losses) resulting from the simulation. 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 non-thermal 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 to detect such interactions for stationary components lying outside of the region in where the losses are dominated by inverse Compton scattering.

Baldwin (1982) wrote that \emph{the distribution of sources in the radio luminosity, $P$, overall physical size, $D$, diagram} could be considered as \emph{the radio astronomer's $H-R$ diagram}. However, unlike the case of stars, not only the intrinsic properties of the jets, but also those of the host galaxy and the intergalactic medium are relevant to explain the evolutionary tracks of radio radio sources. In this contribution I review the current status of our understanding of the evolution of radio sources from a theoretical and numerical perspective, using the $P-D$ diagram as a framework. An excess of compact (linear size $\leq 10$~kpc) sources could be explained by low-power jets being decelerated within the host galaxy, as shown by recent numerical simulations. Finally, I discuss the possible tracks that radio sources may follow within this diagram, and the physical processes that can explain the different tracks.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

The jet-counterjet system of the closest radio-loud active galaxy Centaurus A (Cen A) can be studied with Very Long Baseline Interferometry (VLBI) on unprecedented small linear scales of ~0.018pc. These high-resolution observations provide essential information on jet emission and propagation within the inner parsec of an AGN jet. We present the results of a kinematic study performed within the framework of the Southern-hemisphere AGN monitoring program TANAMI. Over 3.5years, the evolution of the central-parsec jet structure of Cen A was monitored with VLBI. These observations reveal complex jet dynamics which are well explained by a spine-sheath structure supported by the downstream acceleration occurring where the jet becomes optically thin. Both moving and stationary jet features are tracked. A persistent local minimum in surface brightness suggests the presence of an obstacle interrupting the jet flow, which can be explained by the interaction of the jet with a star at a distance of ~0.4pc from the central black hole. We briefly discuss possible implications of such an interaction regarding the expected neutrino and high-energy emission and the effect on a putative planet.

(abridged)The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary suffers periodic and random variations of orbital and non-linear dynamical origin. For the 1st time, we simulate in 3 D the interaction of isotropic stellar and relativistic pulsar winds along 1 full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolution and size, and of different EoOs: a gamma-constant ideal gas, and an ideal gas with gamma dependent on temperature. We carry out, with the code PLUTO, relativistic HD simulations in 2 and 3 D of the interaction of a slow wind and a relativistic wind with Gamma=2 along 1 full orbit up to ~100 x the binary size. The different 2-D simulations are carried out with equal and larger grid resolution and size, and 1 of them is done with a more realistic equation of state, than in 3 D. The simulations in 3 D confirm previous results in 2 D. The shocked flows are subject to a faster instabilities growth in 3 D, which enhances the presence of shocks, mixing, and large-scale disruption. In 2 D, higher resolution simulations confirm lower resolution results, simulations with larger grid sizes strengthen the case for the loss of global coherence of the shocked-wind structure, and simulations with 2 different EoOs yield very similar results. In addition to the KHI, we find that the Richtmyer-Meshkov and the RTI are likely acting together in the shocked flow evolution. Simulations in 3 D confirm that the interaction of stellar and pulsar winds yields structures that evolve non-linearly and get strongly entangled. The evolution is accompanied by strong kinetic energy dissipation, rapid changes in flow orientation and speed, and turbulent motion. The results strengthen the case for the loss of global coherence of the shocked structure on large scales, although higher pulsar wind speed simulations are needed.

Binary systems containing a massive star and a non-accreting pulsar present strong interaction between the stellar and the pulsar winds. The properties of this interaction, which largely determine the non-thermal radiation in these systems, strongly depend on the structure of the stellar wind, which can be clumpy or strongly anisotropic, as in Be stars. We study numerically the influence of inhomogeneities in the stellar wind on the structure of the two-wind interaction region. We carried out for the first time axisymmetric, relativistic hydrodynamical simulations, with Lorentz factors of ~6 and accounting for the impact of instabilities, to study the impact in the two-wind interaction structure of an over-dense region of the stellar wind. We also followed the evolution of this over-dense region or clump as it faces the impact of the pulsar wind. For typical system parameters, and adopting a stellar wind inhomogeneity with a density contrast >~10, clumps with radii of a few percent of the binary size can significantly perturb the two-wind interaction region, pushing the two-wind interface to <~40% of the initial distance to the pulsar. After it is shocked, the inhomogeneity quickly expands and is disrupted when it reaches the smallest distance to the pulsar. It eventually fragments, being advected away from the binary system. The whole interaction region is quite unstable, and the shocked pulsar wind can strongly change under small perturbations. We confirm the sensitive nature of the two-wind interaction structure to perturbations, in particular when the stellar wind is inhomogeneous. For realistic over-dense regions of the stellar wind, the interaction region may shrink by a factor of a few, with the shocked flow presenting a complex spatial and temporal pattern. This can lead to strong variations in the non-thermal radiation.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

Blazars are among the most powerful extragalactic objects, as a sub-class of active galactic nuclei. They launch relativistic jets and their emitted radiation shows strong variability across the entire electro-magnetic spectrum. The mechanisms producing the variability are still controversial and different models have been proposed to explain the observed variations in multi-frequency blazar light curves.We investigate the capabilities of the classical shock-in-jet model to explain and reconstruct the observed evolution of flares in the turnover frequency turnover flux density plane and their frequency-dependent light curve parameters. With a detailed parameter space study we provide the framework for future, detailed comparisons of observed flare signatures with the shock-in-jet scenario. Based on the shock model we compute synthetic single-dish light curves at different radio frequencies (2.6 to 345 GHz) and for different physical conditions in a conical jet (e.g. magnetic field geometry and Doppler factor). From those we extract the slopes of the different energy loss stages within the $\nu_\mathrm{m}$-$S_\mathrm{m}$ plane and deduce the frequency-dependence of different light curve parameters such as flare amplitude, time scale and cross-band delays. The evolution of the Doppler factor along the jet has the largest influence on the evolution of the flare and on the frequency-dependent light curve parameters. The synchrotron stage can be hidden in the Compton or in the adiabatic stage, depending mainly on the evolution of the Doppler factor, which makes it difficult to detect its signature in observations. In addition, we show that the time lags between different frequencies can be used as an efficient tool to better constrain the physical properties of these objects.

Extragalactic jets are a common feature of radio-loud active galaxies. The nature of the observed jets in relation to the bulk flow is still unclear. In particular it is not clear whether the observations of parsec-scale jets using the very long baseline interferometric technique (VLBI) reveal wave-like structures that develop and propagate along the jet, or trace the jet flow itself. In this contribution I review the evidence collected during the last years showing that the ridge-lines of helical radio-jets do not correspond to observational artifacts. This conclusion was reached by studying a number of VLBI observations of the radio jet in the quasar S5~0836+710 at different frequencies and epochs. The ridge-line of the emission in the jet coincides at all frequencies within the errors. Moreover, small differences between the ridge-lines as observed at different epochs reveal wave-like motion transversal to the jet propagation axis. I also discuss similar results, albeit with different interpretations, obtained by other authors. The current challenge is to measure the propagation velocities of these waves and to try to characterise them in terms of simple perturbations or Kelvin-Helmholtz instability, which would help understanding the physical conditions of the flow where the waves develop. This problem can only be tackled by high-resolution observations such as those that can be achieved by the space radio-antenna Radioastron.

Relativistic jets in active galactic nuclei represent one of the most powerful phenomena in the Universe. They form in the surroundings of the supermassive black holes as a by-product of accretion onto the central black hole in active galaxies. The flow in the jets propagates at velocities close to the speed of light. The distance between the first part of the jet that is visible in radio images (core) and the black hole is still a matter of debate. Only very-long-baseline interferometry observations resolve the innermost compact regions of the radio jet. Those can access the jet base, and combining data at different wavelenghts, address the physical parameters of the outflow from its emission. We have performed an accurate analysis of the frequency-dependent shift of the VLBI core location for a multi-wavelength set of images of the blazar CTA 102 including data from 6 cm down to 3 mm. The measure of the position of the central black hole, with mass $\sim 10^{8.93}\,M_\odot$, in the blazar CTA 102 reveals a distance of $\sim 8\times10^4$ gravitational radii to the 86 GHz core, in agreement with similar measures obtained for other blazars and distant radio galaxies, and in contrast with recent results for the case of nearby radio galaxies, which show distances between the black hole and the radio core that can be two orders of magnitude smaller.

Centaurus A is the closest radio-loud active galaxy. Very Long Baseline Interferometry (VLBI) enables us to study the jet-counterjet system on milliarcsecond (mas) scales, providing essential information for jet emission and propagation models. We study the evolution of the central parsec jet structure of Cen A over 3.5 years. The proper motion analysis of individual jet components allows us to constrain jet formation and propagation and to test the proposed correlation of increased high energy flux with jet ejection events. Cen A is an exceptional laboratory for such detailed study as its proximity translates to unrivaled linear resolution, where 1 mas corresponds to 0.018 pc. The first 7 epochs of high-resolution TANAMI VLBI observations at 8 GHz of Cen A are presented, resolving the jet on (sub-)mas scales. They show a differential motion of the sub-pc scale jet with significantly higher component speeds further downstream where the jet becomes optically thin. We determined apparent component speeds within a range of 0.1c to 0.3c, as well as identified long-term stable features. In combination with the jet-to-counterjet ratio we can constrain the angle to the line of sight to ~12{\deg} to 45{\deg}. The high resolution kinematics are best explained by a spine-sheath structure supported by the downstream acceleration occurring where the jet becomes optically thin. On top of the underlying, continuous flow, TANAMI observations clearly resolve individual jet features. The flow appears to be interrupted by an obstacle causing a local decrease in surface brightness and a circumfluent jet behavior. We propose a jet-star interaction scenario to explain this appearance. The comparison of jet ejection times with high X-ray flux phases yields a partial overlap of the onset of the X-ray emission and increasing jet activity, but the limited data do not support a robust correlation.

The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolution of jets with different powers in terms of the spatial and temporal scales of energy deposition. According to our results, the observed morphology in X-ray cavities requires that an important fraction of the jet's energetic budget is in the form of internal energy. Thus, light, lepton-dominated jets are favoured. In all cases, heating is mainly produced by shocks. Cavity overpressure is sustained by an important population of thermal particles. Our simulations reproduce the cool-core structure in projected, luminosity-weighted temperature. We have performed an additional simulation of a slow, massive jet and discuss the differences with its relativistic counterparts. Important qualitative and quantitative differences are found between the non-relativistic and the relativistic jets. Our conclusions point towards a dual-mode of AGN kinetic feedback, depending on the jet power.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic 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 (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, 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 FRI jets within the first kiloparsec. However, powerful FRI 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 kpc scales is initiated by mass loading in the weaker FRI 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.

Extragalactic jets in active galactic nuclei (AGN) are divided into two morphological types, namely Fanaroff-Riley I (FRI) and Fanaroff-Riley II (FRII). The former show decollimated structure at the kiloparsec scales and are thought to be decelerated by entrainment within the first kiloparsecs of evolution inside the host galaxy. The entrainment and deceleration can be, at least partly, due to the interaction of jets with stellar winds and gas clouds that enter in the jet as they orbit around the galactic centre. In this contribution, I review recent simulations to study the dynamic effect of entrainment from stellar winds in jets and the direct interaction of jets with gas clouds and stellar winds. I also briefly describe the importance of these interactions as a possible scenario of high-energy emission from extragalactic jets.

To obtain a better understanding of the location and mechanisms for the production of the gamma-ray emission in jets of AGN we present a detailed study of the HST-1 structure, 0.8 arcsec downstream the jet of M87, previously identified as a possible candidate for TeV emission. HST-1 shows a very peculiar structure with superluminal as well as possible stationary sub-components, and appears to be located in the transition from a parabolic to a conical jet shape, presumably leading to the formation of a recollimation shock. This scenario is supported by our new RHD simulations in which the interaction of a moving component with a recollimation shock leads to the appearance of a new superluminal component. To discern whether HST-1 is produced by a recollimation shock or some other MHD instability, we present new polarimetric 2.2 and 5 GHz VLBA, as well as 15, 22 and 43 GHz JVLA observations obtained between November 2012 and March 2013.

The dynamics and stability of extragalactic jets may be strongly influenced by small (and probable) differences in pressure between the jet and the ambient and within the jet itself. The former give rise to expansion and recollimation of the jet. This occurs in the form of conical shocks, or Mach disks, if the pressure difference is large enough. Pressure asymmetries within the jet may trigger the development of helical patterns via coupling to kink current-driven instability, or to helical Kelvin-Helmholtz instability, depending on the physical conditions in the jet. I summarize here the evidence collected during the last years on the presence of recollimation shocks and waves in jets. In the jet of CTA 102 evidence has been found for (traveling)shock-(standing)shock interaction in the core-region (0.1 mas from the core), using information from the light-curve of the source combined with VLBI data. The conclusions derived have been confirmed by numerical simulations combined with emission calculations that have allowed to study the spectral evolution of the perturbed jet. Helical structures can also be identified in radio-jets. The ridge-line of emission of the jet of S5~0836+710 has been identified as a physical structure corresponding to a wave developing in the jet flow. I review here the evidence that has allowed to reach this conclusion, along with an associated caveat. Current data do not allow to distinguish between magnetic or hydrodynamical instabilities. I finally discuss the importance of these linear and non-linear waves for jet evolution.

One of the open questions in extragalactic jet Astrophysics is related to the nature of the observed radio jet, namely whether it traces a pattern or the flow structure itself. In this paper I summarize the evidence collected for the presence of waves in extragalactic jets. The evidence points towards the peak of emission in helical jets corresponding to pressure-maxima of a wave that is generated within the core region and propagates downstream. Making use of a number of very long baseline interferometry (VLBI) observations of the radio jet in the quasar S5~0836+710 at different frequencies and epochs, Perucho et al. (2012) were able to observe wave-like behavior within the observed radio-jet. The ridge-line of the emission in the jet coincides within the errors at all frequencies. Moreover, small differences between epochs at 15 GHz reveal wave-like motion of the ridge-line transversal to the jet propagation axis. The authors conclude that the helicity is a real, physical structure. I report here on those results and discuss this result in the light of new results recently announced by other authors that confirm the presence of waves in the close-by object BL Lac (Cohen et al., in preparation), and the possible magnetic or hydrodynamic nature.

The temporal and spatial spectral evolution of the jets of AGN can be studied with multi-frequency, multi-epoch VLBI observations. The combination of both, morphological and spectral parameters can be used to derive source intrinsic physical properties such as the magnetic field and the non-thermal particle density. In the first two papers of this series, we analyzed the single-dish light curves and the VLBI kinematics of the blazar CTA 102 and suggested a shock-shock interaction between a traveling and a standing shock wave as a possible scenario to explain the observed evolution of the component associated to the 2006 flare. In this paper we investigate the core-shift and spectral evolution to test our hypothesis of a shock-shock interaction. We used 8 multi-frequency VLBA observations to analyze the temporal and spatial evolution of the spectral parameters during the flare. We observed CTA 102 between May 2005 and April 2007 using the VLBA at six different frequencies spanning from 2 - 86 GHz. After the calibrated VLBA images were corrected for opacity, we performed a detailed spectral analysis. From the derived values we estimated the magnetic field and the density of the relativistic particles. The detailed analysis of the opacity shift reveals that the position of the jet core is proportional to nu^-1 with some temporal variations. The value suggests possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions. From the variation of the physical parameters we deduced that the 2006 flare in CTA 102 is connected to the ejection of a new traveling feature (t=2005.9) and the interaction between this shock wave and a stationary structure around 0.1 mas from the core. The source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet.

The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis of the single--dish light curves revealed a possible shock-shock interaction as a scenario for the 2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz - 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent position of the core slightly deviates from the typical $\nu^{-1}$-dependence during the flaring state. This behaviour may reflect a mismatch between magnetic energy density and energy density of the relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic field and the particle density. The overall picture of the source could be best interpreted as an over-pressured jet perturbed by traveling shock waves.

We make use of VLBI observations of the radio jet in the quasar S5 0836+710 at different frequencies and epochs to study its properties. The jet shows helical structure at all frequencies. The ridge-line of the emission in the jet coincides at all frequencies and epochs, within the errors. We conclude that the helicity is a real, physical structure. Small differences between epochs reveal wave-like motion of the ridge-line transversal to the jet propagation axis. These transversal motions are measured to be superluminal. This unphysical result could correspond to a possible small amplitude oscillation of the ridge-line at the radio-core and to large errors in the determination of the positions. In addition, higher resolution images at 15 GHz show that the ridge-line does not coincide exactly with the centre of the radio jet. At arc-second scales, this powerful jet shows non-collimated, irregular structure and a lack of a hot-spot. Following this collection of evidence, we conclude that the ridge-line could be related to a pressure maximum within the jet cross-section, associated with the observed helical pattern that could lead to jet disruption at longer scales. We discuss the implications regarding the nature of the instability and the jet flow itself.

We performed multifrequency multiepoch Very Long Baseline Array (VLBA) observations of the blazar CTA 102 during its 2006 radio flare, the strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical properties of blazars, especially during these flaring events. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA\,102 (Paper I). In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 43 GHz. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible recollimation shocks around 0.1 mas (deprojected 18 pc at a viewing angle of 2.6 degrees) and 6.0 mas (deprojected 1 kpc) from the core. The 43 GHz observations reveal a feature ejected at epoch $t_\mathrm{ej}=2005.9\pm0.2$, which could be connected to the 2006 April radio flare. Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which gives rise to the observed double peak structure in the single-dish light curves reported in Paper I.

The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet structures and obtain information from these structures. For example, it was recently shown that the observed helical structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz instability. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before interaction with the intergalactic medium. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction has provided a complete image of the object at arc-second scales. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation between the VLBI region and the region of interaction with the ambient medium. This result, together with the previous identification of the helical structures in the jet with helical pressure waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not correspond to this classification. The implications of this fact are discussed.

The MOJAVE survey contains 101 quasars with a total of 354 observed radio components that are different from the radio cores, among which 95% move with apparent projected superluminal velocities with respect to the core, and 45% have projected velocities larger than 10c (with a maximum velocity 60c). Doppler boosting effects are analyzed to determine the statistics of the superluminal motions. We integrate over all possible values of the Lorentz factor the values of the kinetic energy corresponding to each component. The calculation of the mass in the ejection is carried out by assuming the minimum energy state. This kinetic energy is multiplied by the frequency at which the portions of the jet fluid identified as "blobs" are produced. Hence, we estimate the average total power released by the quasars in the form of kinetic energy in the long term on pc-scales. RESULTS. A selection effect in which both the core and the blobs of the quasar are affected by huge Doppler-boosting enhancement increases the probability of finding a jet ejected within 10 degrees of the line of sight >~40 times above what one would expect for a random distribution of ejection, which explains the ratios of the very high projected velocities given above. The average total kinetic power of each MOJAVE quasar should be very high to obtain this distribution: ~7E47 erg/s. This amount is much higher than previous estimates of kinetic power on kpc-scales based on the analysis of cavities in X-ray gas or radio lobes in samples of objects of much lower radio luminosity but similar black hole masses. The kinetic power is a significant portion of the Eddington luminosity, on the order of the bolometric luminosity, and proportional on average to square root of the radio luminosity, although this correlation might be induced by Malmquist-like bias.

The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. Potential sites of non-thermal activity are investigated. Relativistic hydrodynamical simulations in two dimensions, performed with the code PLUTO and using the adaptive mesh refinement technique, are used to model interacting stellar and pulsar winds on scales ~80 times the distance between the stars. The hydrodynamical results suggest the suitable locations of sites for particle acceleration and non-thermal emission. In addition to the shock formed towards the star, the shocked and unshocked components of the pulsar wind flowing away from the star terminate by means of additional strong shocks produced by the orbital motion. Strong instabilities lead to the development of turbulence and an effective two-wind mixing in both the leading and trailing sides of the interaction structure, which starts to merge with itself after one orbit. The adopted moderate pulsar-wind Lorentz factor already provides a good qualitative description of the phenomena involved in high-mass binaries with pulsars, and can capture important physical effects that would not appear in non-relativistic treatments. Simulations show that shocks, instabilities, and mass-loading yield efficient mass, momentum, and energy exchanges between the pulsar and the stellar winds. This renders a rapid increase in the entropy of the shocked structure, which will likely be disrupted on scales beyond the simulated ones. Several sites of particle acceleration and low- and high-energy emission can be identified. Doppler boosting will have significant and complex effects on radiation.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

Helical structures are common in extragalactic jets. They are usually attributed in the literature to periodical phenomena in the source (e.g., precession). In this work, we use VLBI data of the radio-jet in the quasar S5 0836+710 and hypothesize that the ridge-line of helical jets like this corresponds to a pressure maximum in the jet and assume that the helically twisted pressure maximum is the result of a helical wave pattern. For our study, we use observations of the jet in S5 0836+710 at different frequencies and epochs. The results show that the structures observed are physical and not generated artificially by the observing arrays. Our hypothesis that the observed intensity ridge-line can correspond to a helically twisted pressure maximum is confirmed by our observational tests. This interpretation allows us to explain jet misalignment between parsec and kiloparsec scales when the viewing angle is small, and also brings us to the conclusion that high-frequency observations may show only a small region of the jet flow concentrated around the maximum pressure ridge-line observed at low frequencies. Our work provides a potential explanation for the apparent transversal superluminal speeds observed in several extragalactic jets by means of transversal shift of an apparent core position with time.

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. 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. 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. 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. 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.

We present the results from an ultra-high-resolution 7mm Very Long Baseline Array (VLBA) study of the relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. Analysis of the image sequence reveals a sharp jet-position-angle swing by >100 deg. during [2004,2006], as viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our images also show for the first time a prominent erratic wobbling behavior of the innermost ~0.4mas of the jet with fluctuations in position angle of up to ~40 deg. over time scales ~2yr. This is accompanied by highly superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as possible origins of the phenomenon on the scales probed by our observations, although such processes may cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet.

The termination structures of the jets of Fanaroff & Riley galaxies are observed to produce extended non-thermal emission in a wide frequency range. The study of these structures can provide valuable insights on the conditions for particle acceleration and radiation at the shock fronts.We have studied the thermal and non-thermal emission that can be expected from the jet termination regions of Fanaroff & Riley type I sources. The broadband emission from these galaxies has been recently extended to include the high-energy gamma-ray domain, owing to the Fermi detection of Cen A lobes. Exploring the physics behind the jet/medium interactions in FRI can provide valuable insights on the conditions for particle acceleration and radiation in the jet termination shocks. Making use of the results of a fully relativistic numerical simulation code of the evolution of a FRI jet we model the expected radiative output and predict spectra and lightcurves of both thermal and non-thermal emission at different source ages.

High-mass microquasars consist of a massive star and a compact object, the latter producing jets that will interact with the stellar wind. The evolution of the jets, and ultimately their radiative outcome, could depend strongly on the inhomogeneity of the wind, which calls for a detailed study. The hydrodynamics of the interaction between a jet and a clumpy wind is studied, focusing on the global wind and single clump-jet interplay. We have performed, using the code \textit{Ratpenat}, three-dimensional numerical simulations of a clumpy wind interacting with a mildly relativistic jet, and of individual clumps penetrating into a jet. For typical wind and jet velocities, filling factors of about > 0.1 are already enough for the wind to be considered as clumpy. An inhomogeneous wind makes the jet more unstable when crossing the system. Kinetic luminosities of the order 1.e37 erg/s allow the jet to reach the borders of a compact binary with an O star, as in the smooth wind case, although with a substantially higher degree of disruption. When able to enter into the jet, clumps are compressed and heated during a time of about their size divided by the sound speed in the shocked clump. Then, clumps quickly disrupt, mass-loading and slowing down the jet. We conclude that moderate wind clumpiness makes already a strong difference with the homogeneous wind case, enhancing jet disruption, mass-loading, bending, and likely energy dissipation in the form of emission. All this can have observational consequences at high-energies and also in the large scale radio jets.

The analysis of the radio light curves of the blazar CTA102 during its 2006 flare revealed a possible interaction between a standing shock wave and a traveling one. In order to better understand this highly non-linear process, we used a relativistic hydrodynamic code to simulate the high energy interaction and its related emission. The calculated synchrotron emission from these simulations showed an increase in turnover flux density, $S_{m}$, and turnover frequency, $\nu_{m}$, during the interaction and decrease to its initial values after the passage of the traveling shock wave.

Relativistic jets carry energy and particles from compact to very large scales compared with their initial radius. This is possible due to their remarkable collimation despite their intrinsic unstable nature. In this contribution, I review the state-of-the-art of our knowledge on instabilities growing in those jets and several stabilising mechanisms that may give an answer to the question of the stability of jets. In particular, during the last years we have learned that the limit imposed by the speed of light sets a maximum amplitude to the instabilities, contrary to the case of classical jets. On top of this stabilising mechanism, the fast growth of unstable modes with small wavelengths prevents the total disruption and entrainment of jets. I also review several non-linear processes that can have an effect on the collimation of extragalactic and microquasar jets. Within those, I remark possible causes for the decollimation and decelleration of FRI jets, as opposed to the collimated FRII's. Finally, I give a summary of the main reasons why jets can propagate through such long distances.

Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our results show no buoyant bubbles, but large cocoon regions compatible with the observed X-ray cavities. The reheating is more efficient and faster than in previous scenarios, and it is produced by the shock wave driven by the jet, that survives for several hundreds of Myrs. Therefore, the X-ray cavities in clusters produced by powerful relativistic jets would remain confined by weak shocks for extremely long periods, whose detection could be an observational challenge.

The Compact Symmetric Objects (CSOs) are small (<1 kiloparsec) and powerful extragalactic radio sources showing emission on both sides of an active galactic nucleus and no signs of strong relativistic beaming. They may be young radio sources, progenitors of large FRII radio galaxies. We aim to study the statistical properties of CSOs by constructing and investigating a new large sample of CSO candidates on the basis of dual-frequency, parsec-scale morphology. For the candidate selection we utilized VLBI data for 4170 extragalactic objects obtained simultaneously at 2.3 and 8.6 GHz (S and X band) within the VLBA Calibrator Survey 1-6 and the Research and Development - VLBA projects. Properties of their broad-band radio spectra were characterized by using RATAN-600 observations. Numerical modeling was applied in an attempt to explain the observed effects. A sample of 64 candidate CSOs is identified. The median two-point S-X band spectral index of parsec-scale hotspots is found to be -0.52; with the median brightness temperature ~10^9 K at X band. Statistical analysis reveals a systematic difference between positions of brightest CSO components (associated with hotspots) measured in the S and X bands. The distance between these components is found to be on average 0.32+/-0.06 mas greater at 8.6 GHz than at 2.3 GHz. This difference in distances cannot be explained by different resolutions at the S and X bands. It is a manifestation of spectral index gradients across CSO components, which may potentially provide important physical information about them. Despite our detailed numerical modeling of a CSO hotspot, the model was not able to reproduce the magnitude of the observed positional difference. A more detailed modeling may shed light on the origin of the effect.

High energy emission can be produced in the interaction sites of both galactic and extragalactic jets with the surrounding medium. We have developed a radiative model that accounts for the continuous injection of relativistic electrons in the forward, reverse and recollimation shocks developed in the shell, cocoon and reconfinement interaction regions, respectively. We also performed hydrodynamical simulations to establish the physical properties in both galactic and extragalactic systems. The resulting non-thermal emission is studied assuming different values for the jet power, the external mass density and the source age for both FR I galaxies and galactic microquasars. The obtained fluxes are compared to current instrument sensitivities at radio, X-ray and gamma-ray bands.

The morphologies of detected jets in X-ray binaries are almost as diverse as their number. This is due to different jet properties and ambient media that these jets encounter. It is important to understand the physics of these objects and to obtain information about possible sites suitable for particle acceleration in order to explain the observations at very high energies. Here I present the results obtained from the first relativistic hydrodynamical simulations of jets in high-mass microquasars. Our results allow us to make estimates for the emission originated in different sites of the whole structure generated by the jets. These works represent a first step in trying to obtain a deeper understanding of the physics and emission processes related with jets in high-mass microquasars.

Context: The blazar CTA 102 (z=1.037) underwent a historical radio outburst in April 2006. This event offered a unique chance to study the physical properties of the jet. Aims: We used multifrequency radio and mm observations to analyze the evolution of the spectral parameters during the flare as a test of the shock-in-jet model under these extreme conditions. Methods: For the analysis of the flare we took into account that the flaring spectrum is superimposed on a quiescent spectrum. We reconstructed the latter from archival data and fitted a synchrotron self-absorbed distribution of emission. The uncertainties of the derived spectral parameters were calculated using Monte Carlo simulations. The spectral evolution is modeled by the shock-in-jet model, and the derived results are discussed in the context of a geometrical model (varying viewing angle) and shock-shock interaction. Results: The evolution of the flare in the turnover frequency-turnover flux density plane shows a double peak structure. The nature of this evolution is dicussed in the frame of shock-in-jet models. We discard the generation of the double peak structure in the turnover frequency-turnover flux density plane purely based on geometrical changes (variation of the Doppler factor). The detailed modeling of the spectral evolution favors a shock-shock interaction as a possible physical mechanism behind the deviations from the standard shock-in-jet model.

The environment of high-mass X-ray binaries can be characterized either by the SNR that forms these systems, or by the wind from the companion massive star. These regions should be tenuous but very hot, and surrounded by a dense and cold shocked ISM shell. The interaction between the jet and such a complex medium, also affected by the system proper motion, can lead to very different jet termination structures. The evolution of the jet termination regions during the life of a high-mass microquasar is simulated to improve the present understanding of these structures. Also, the evolving emission characteristics are modeled to inform potential observational campaigns for this class of object. We have performed 2D numerical simulations of jets propagating in different scenarios, corresponding to different epochs after the formation of the high-mass X-ray binary, using the code Ratpenat. We have also made simple estimates of the non-thermal emission that could be produced in the jet termination regions. We find that, in the way through the hot and tenuous medium of the shocked wind/SNR ejecta, the jet suffers recollimation shocks in which it loses part of its thrust and ends in a strong shock inflating a hot cocoon. The jet head propagates with a speed similar to the medium sound speed, until it eventually reaches the denser and colder shocked ISM and the unperturbed ISM later on. In these last stages of evolution, the jet is significantly slowed down and can be disrupted. For relatively old sources, the microquasar peculiar velocity becomes important, leading to complete jet destruction. Extended non-thermal radiation can be generated in the jet termination regions, being hard X-rays and TeV photons the best suited wavelengths to observe these structures.

We have found evidence for interaction between a standing and a traveling shock in the jet of the blazar CTA 102. Our result is based in the study of the spectral evolution of the turnover frequency-turnover flux density plane. The radio/mm light curves were taken during a major radio outburst in April 2006.

The radio light curve and spectral evolution of the blazar CTA 102 during its 2006 outburst can be rather well explained by the standard shock-in-jet model. The results of a pixel-to-pixel spectral analysis of multi-frequency VLBI images, together with kinematics derived from the MOJAVE survey lead to the picture of an over-pressured jet with respect to the ambient medium. The interaction of a traveling shock wave with a standing one is a possible scenario which could explain the observed spectral behaviour

The blazar CTA 102 underwent a major radio flare in April 2006. We used several 15 GHz VLBI observations from the MOJAVE program to investigate the influence of this extreme event on jet kinematics. The result of modeling and analysis lead to the suggestion of an interaction between traveling and standing shocks 0.2 mas away from the VLBI core.

Recent observations in X-rays and gamma-rays of nearby FRI radio galaxies have raised the question of the origin of the emission detected in the termination structures of their jets. The study of these structures can give information on the conditions for particle acceleration and radiation at the front shocks. In addition, an evolutionary scenario can help to disentangle the origin of the detected X-ray emission in young FRI sources, like some Gigahertz Peaked Spectrum AGNs. This work focuses on the nature and detectability of the radiation seen from the termination regions of evolving FRI jets. We use the results of a relativistic, two-dimensional numerical simulation of the propagation of an FRI jet, coupled with a radiation model, to make predictions for the spectra and lightcurves of the thermal and non-thermal emission at different stages of the FRI evolution. Our results show that under moderate magnetic fields, the synchrotron radiation would be the dominant non-thermal channel, appearing extended in radio and more compact in X-rays, with relatively small flux variations with time. The shocked jet synchrotron emission would dominate the X-ray band, although the shocked ISM/ICM thermal component alone may be significant in old sources. Inverse Compton scattering of CMB photons could yield significant fluxes in the GeV and TeV bands, with a non-negligible X-ray contribution. The IC radiation would present a bigger angular size in X-rays and GeV than in TeV, with fluxes increasing with time. We conclude that the thermal and non-thermal broadband emission from the termination regions of FRI jets could be detectable for sources located up to distances of a few 100 Mpc.

We present the first three-dimensional simulations of the evolution of a microquasar jet inside the binary-star system. The aim is to study the interaction of these jets with the stellar wind from a massive companion and the possible locations of high-energy emission sites. We have simulated two jets with different injection power in order to give a hint on the minimum power required for the jet to escape the system and become visible in larger scales. In the setup, we include a massive star wind filling the grid through which the jet evolves. We show that jets should have powers of the order of $10^{37}\rm{erg/s}$ or more in order not to be destroyed by the stellar wind. The jet-wind interaction results in regions in which high energy emission could be produced. These results imply the possible existence of a population of X-ray binaries not detected in the radio band due to jet disruption inside the region dominated by the stellar wind.

The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions (2D), it is our aim to test and extend the conclusions of such works to three dimensions. We present numerical simulations for the study of the stability properties of 3D, sheared, relativistic flows. This work uses a fully parallelized code Ratpenat that solves equations of relativistic hydrodynamics in 3D. The results of the present simulations confirm those in 2D. We conclude that the growth of resonant modes in sheared relativistic flows could be important in explaining the long-term collimation of extragalactic jets.

High-mass microquasars may produce jets that will strongly interact with surrounding stellar winds on binary system spatial scales. We study the dynamics of the collision between a mildly relativistic hydrodynamical jet of supersonic nature and the wind of an OB star. We performed numerical 3D simulations of jets that cross the stellar wind with the code Ratpenat. The jet head generates a strong shock in the wind, and strong recollimation shocks occur due to the initial overpressure of the jet with its environment. These shocks can accelerate particles up to TeV energies and produce gamma-rays. The recollimation shock also strengthens jet asymmetric Kelvin-Helmholtz instabilities produced in the wind/jet contact discontinuity. This can lead to jet disruption even for jet powers of several times $10^{36}$ erg s$^{-1}$. High-mass microquasar jets likely suffer a strong recollimation shock that can be a site of particle acceleration up to very high energies, but also eventually lead to the disruption of the jet.

The high kinetic energy outflowing in the jets of microquasars is delivered to the surrounding interstellar medium. This energy input can cause the formation of bow shocks and cocoons that may be detectable from radio to gamma-ray energies. Evidences or hints of emission from jet/medium interactions have been reported for some sources, but little has been done regarding the theoretical modeling of the resulting non-thermal emission. We have developed an analytical model based on those successfully applied to extragalactic sources or the interaction of AGN jets with their surroundings. Focusing on the adiabatic phase of the growing structures, we give estimations of the expected luminosities through synchrotron, relativistic Bremsstrahlung and inverse Compton processes. We conclude that the interaction structures may be detectable at radio wavelengths, while extreme values for the jet kinetic power, the source age and the medium density are required to make the emission at high and very high energies detectable.

The interaction of microquasar jets with their environment can produce non-thermal radiation as is the case for extragalactic outflows impacting on their surroundings. We have developed an analytical model based on those successfully applied to extragalactic sources. The jet is taken to be a supersonic and mildly relativistic hydrodynamical outflow. We focus on the jet/shocked medium structure when being in its adiabatic phase, and assume that it grows in a self-similar way. We calculate the fluxes and spectra of the radiation produced via synchrotron, Inverse Compton and relativistic Bremsstrahlung processes by electrons accelerated in strong shocks. A hydrodynamical simulation is also performed to further investigate the jet interaction with the environment and check the physical parameters used in the analytical model. We conclude that microquasar jet termination regions could be detectable at radio wavelengths for current instruments sensitive to arcminute scales while at X-rays the expected luminosities are moderate, although the emitter is more compact than the radio one. The radiation at gamma-ray energies may be within the detection limits of the next generation of satellite and ground-based instruments.

We present results on the modelling of the ejection of a superluminal component in the jet of 3C111. We propose that the component is generated by an injection of dense material followed by a decrease in the injection rate of bulk particles in the jet. Our model is supported by 1D relativistic hydrodynamics and emission simulations, and is capable of reproducing the brightness evolution of two features, as revealed by 15 GHz VLBA observations. We show that other scenarios, such as an increase of the Lorentz factor in the material of the perturbation, fails to reproduce the observed evolution of this flare.

We present a possible scenario for the ejection of a superluminal component in the jet of the Broad Line Radio Galaxy 3C111 in early 1996. VLBI observations at 15 GHz discovered the presence of two jet features on scales smaller than one parsec. The first component evolves downstream, whereas the second one fades out after 1 parsec. We propose the injection of a perturbation of dense material followed by a decrease in the injection rate of material in the jet as a plausible explanation. This scenario is supported by 1D relativistic hydrodynamics and emission simulations. The perturbation is modeled as an increase in the jet density, without modifying the original Lorentz factor in the initial conditions. We show that an increase of the Lorentz factor in the material of the perturbation fails to reproduce the observed evolution of this flare. We are able to estimate the lifetime of the ejection event in 3C111 to be 36\pm7 days.

We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources. We show that the relative brightness of the knots associated with the recollimation shocks decreases with increasing jet magnetization, suggesting that jets presenting stationary components may have a relatively weak magnetization, with magnetic fields of the order of equipartition or below.

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}$.

We present the first results on the analysis of the structures observed in the jet of the quasar 0836+710. We obtain the ridge lines of the jet at different epochs and several frequencies. We interpret the oscillatory structures obtained as waves that can be attached to the growth of instabilities. We explain how to derive information on the nature and origin of these structures by fitting together the ridge lines at different epochs and frequencies. Finally we show the predictive power of this approach: by generating an artificial wave and applying the corresponding relativistic and projection effects we show that apparent changes in the jet direction in the inner regions of jets can be attached to the transversal motion of structures.

In 1996, a major radio flux-density outburst occured in the broad-line radio galaxy 3C111. It was followed by a particularly bright plasma ejection associated with a superluminal jet component, which has shaped the parsec-scale structure of 3C111 for almost a decade. Here, we present results from 18 epochs of Very Long Baseline Array (VLBA) observations conducted since 1995 as part of the VLBA 2 cm Survey and MOJAVE monitoring programs. This major event allows us to study a variety of processes associated with outbursts of radio-loud AGN in much greater detail than has been possible in other cases: the primary perturbation gives rise to the formation of a leading and a following component, which are interpreted as a forward and a backward-shock. Both components evolve in characteristically different ways and allow us to draw conclusions about the work flow of jet-production events; the expansion, acceleration and recollimation of the ejected jet plasma in an environment with steep pressure and density gradients are revealed; trailing components are formed in the wake of the primary perturbation possibly as a result of coupling to Kelvin-Helmholtz instability pinching modes from the interaction of the jet with the external medium. The interaction of the jet with its ambient medium is further described by the linear-polarization signature of jet components traveling along the jet and passing a region of steep pressure/density gradients.

In high-mass microquasars (HMMQ), strong interactions between jets and stellar winds at binary system scales could occur. In order to explore this possibility, we have performed numerical 2-dimensional simulations of jets crossing the dense stellar material to study how the jet will be affected by these interactions. We find that the jet head generates strong shocks in the wind. These shocks reduce the jet advance speed, and compress and heat up jet and wind material. In addition, strong recollimation shocks can occur where pressure balance between the jet side and the surrounding medium is reached. All this, altogether with jet bending, could lead to the destruction of jets with power $<10^{36} \rm{erg/s}$. The conditions around the outflow shocks would be convenient for accelerating particles up to $\sim $TeV energies. These accelerated particles could emit via synchrotron and inverse Compton (IC) scattering if they were leptons, and via hadronic processes in case they were hadrons.

NRAO 150 -a compact and bright radio to mm source showing core/jet structure- has been recently identified as a quasar at redshift z=1.52 through a near-IR spectral observation. To study the jet kinematics on the smallest accessible scales and to compute the first estimates of its basic physical properties, we have analysed the ultra-high-resolution images from a new monitoring program at 86 GHz and 43 GHz with the GMVA and the VLBA, respectively. An additional archival and calibration VLBA data set, covering from 1997 to 2007, has been used. Our data shows an extreme projected counter-clock-wise jet position angle swing at an angular rate of up to ~11 deg./yr within the inner ~31 pc of the jet, which is associated with a non-ballistic superluminal motion of the jet within this region. The results suggest that the magnetic field could play an important role in the dynamics of the jet in NRAO 150, which is supported by the large values of the magnetic field strength obtained from our first estimates. The extreme characteristics of the jet swing make NRAO 150 a prime source to study the jet wobbling phenomenon.

The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, in agreement with recent X-ray observations that show the presence of bow shocks in FRI sources. The jet expands until it becomes underpressured with respect to the ambient medium, and then recollimates. Subsequent oscillations around pressure equilibrium and generation of standing shocks lead to the mass loading and disruption of the jet flow. We derive an estimate for the minimum age of the source of $t>1. 10^8 \rm{yrs}$, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak CSS sources may be the young counterparts of FRIs. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario.

Studying transversal structure in extragalactic jets is crucial for understanding their physics. The Japanese led space VLBI project VSOP has offered arguably the best opportunity for such studies, by reaching baseline lengths of up to 36,000 km and resolving structures down to an angular size of 0.3 mas at 5 GHz. VSOP observations of the jet in 0836+710 at 1.6 and 5 GHz have enabled tracing the radial structure of the flow on scales from 2 mas to 200 mas and determining the wavelengths of individual oscillatory modes responsible for the formation of the structure observed. We conclude that these modes are produced by Kelvin-Helmholtz instability in a sheared relativistic flow. Our results point towards the stratification of the jet and the growth of different modes at different jet radii. We also discuss the implications of the driving frequency on the physics of the active nucleus of the quasar.

Qualitatively new aspects of the (linear and non-linear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20; specific internal energies $\approx 60c^2$). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer ii) to confirm the result of the linear analysis with numerical simulations and, iii) more interestingly, to follow the instability development through the non-linear regime. We find that very high-order reflection modes with dominant growth rates can modify the global, long-term stability of the relativistic flow. We discuss the dependence of these resonant modes on the jet flow Lorentz factor and specific internal energy, and on the shear layer thickness. The results could have potential applications in the field of extragalactic relativistic jets.

Studying the internal structure of extragalactic jets is crucial for understanding their physics. The Japanese-led space VLBI project VSOP has presented an opportunity for such studies, by reaching baseline lengths of up to 36,000 km and resolving structures down to an angular size of $\approx 0.3$ mas at 5 GHz. VSOP observations of the jet in 0836+710 at 1.6 and 5 GHz have enabled tracing of the radial structure of the flow on scales from 2 mas to 200 mas along the jet and determination of the wavelengths of individual oscillatory modes responsible for the formation of the structure observed. We apply linear stability analysis to identify the oscillatory modes with modes of Kelvin-Helmholtz instability that match the wavelengths of the structures observed. We find that the jet structure in 0836+710 can be reproduced by the helical surface mode and a combination of the helical and elliptic body modes of Kelvin-Helmholtz instability. Our results indicate that the jet is substantially stratified and different modes of the instability grow inside the jet at different distances to the jet axis. The helical surface mode can be driven externally, and we discuss the implications of the driving frequency on the physics of the active nucleus in 0836+710.

In this paper, we present new VLBA observations of the radio jet in the quasar S5 0836+710 at 8 and 22 GHz. The identification of the ridge lines allow us to interpret the jet structure in terms of Kelvin-Helmholtz instability. Combined with previous epochs of VLBA and VSOP data at 1.6 and 5 GHz, these new observations will allow us to study the evolution of the instabilities in the jet. We have detected signatures of possible jet disruption in the jet due to the growth of instabilities, which points towards a possible morphological classification of the source as an FRI, contrary to what previously thought in terms of luminosity criteria.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well-known quasar 3C 273. We present here the results from a 3D relativistic hydrodynamics numerical simulation based on the parameters given for the jet by Lobanov & Zensus (2001), and one in which the effects of jet precession and the injection of discrete components have been taken into account. We compare the model with the structures observed in 3C 273 using very long baseline interferometry and constrain the basic properties of the flow. We find growing perturbation modes in the simulation with similar wavelengths to those observed, but with a different set of wave speeds and mode identification. If the observed longest helical structure is produced by the precession of the flow, longer precession periods should be expected. Our results show that some of the observed structures could be explained by growing Kelvin-Helmholtz instabilities in a slow moving region of the jet. However, we point towards possible errors in the mode identification that show the need of more complete linear analysis in order to interpret the observations. We conclude that, with the given viewing angle, superluminal components and jet precession cannot explain the observed structures.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well known quasar 3C 273. We compare the model with the radio structure observed in 3C 273 on parsec scales using very long baseline interferometry (VLBI) and constrain the basic properties of the flow. Our results show that Kelvin-Helmholtz instabilities are the most plausible mechanism to generate the observed structures.

The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external medium and more general perturbations. The models considered span a wide range of Lorentz factors ($2.5-20$) and internal energies ($0.08 c^2-60 c^2$) and are classified into three classes according to the main characteristics of their long-term, non-linear evolution. We observe a clear separation of these three groups in a relativistic Mach-number Lorentz-factor plane. Jets with a low Lorentz factor and small relativistic Mach number are disrupted after saturation. Those with a large Lorentz factor and large relativistic Mach number are the stablest, due to the appearance of short wavelength resonant modes which generate local mixing and heating in the shear layer around a fast, unmixed core, giving a plausible solution for the problem of the long-term stability of relativistic jets. A third group is present between them, including jets with intermediate values of Lorentz factor and relativistic Mach number, which are disrupted by a slow process of mixing favored by an efficient and continuous conversion of kinetic into internal energy. In the long term, all the models develop a distinct transversal structure (shear/transition layers) as a consequence of KH perturbation growth, depending on the class they belong to. The properties of these shear layers are analyzed in connection with the parameters of the original jet models.

We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets.

In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models undergoing qualitatively different non-linear evolution are clearly grouped in well-separated regions in a jet Lorentz factor/jet-to-ambient enthalpy diagram. Jets with a low Lorentz factor and small enthalpy ratio are disrupted by a strong shock after saturation. Those with a large Lorentz factor and enthalpy ratio are unstable although the process of mixing and momentum exchange proceeds to a longer time scale due to a steady conversion of kinetic to internal energy in the jet. In these cases, the high value of the initial Lorentz seems to prevent transversal velocity from growing far enough to generate the strong shock that breaks the slower jets. Finally, jets with either high Lorentz factors and small enthalpy ratios or low Lorentz factors and large enthalpy ratios appear as the most stable.

Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is considered to proceed in ram pressure equilibrium with the ambient medium. Finally, we also assume that the jets feeding the hot spots are relativistic and have a time dependent power. Comparison with observational data allows to interpret the CSO-FRII evolution in terms of decreasing jet power with time.

The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three different values of the jet Lorentz factor (5, 10 and 20) and a few different values of specific internal energy of the jet matter (from 0.08 to $60.0 c^2$) have been considered. Figures illustrating the evolution of the perturbations are also shown.

A set of simulations concerning the influence of internal energy on the stability of relativistic jets is presented. Results show that perturbations saturate when the amplitude of the velocity perturbation approaches the speed of light limit. Also, contrary to what predicted by linear stability theory, jets with higher specific internal energy appear to be more stable.

We present a model to determine the physical parameters of jets and hot spots of a sample of CSOs under very basic assumptions like synchrotron emission and minimum energy conditions. Based on this model we propose a simple evolutionary scenario for these sources assuming that they evolve in ram pressure equilibrium with the external medium and constant jet power. The parameters of our model are constrained from fits of observational data (radio luminosity, hot spot radius and hot spot advance speed) versus projected linear size. From these plots we conclude that CSOs evolve self-similarly and that their radio luminosity increases with linear size along the first kiloparsec. Assuming that the jets feeding CSOs are relativistic from both kinematical and thermodynamical points of view, we use the values of the pressure and particle number density within the hot spots to estimate the fluxes of momentum (thrust), energy, and particles of these relativistic jets. The mean jet power obtained in this way is within an order of magnitude that inferred for FRII sources, which is consistent with CSOs being the possible precursors of large doubles. The inferred flux of particles corresponds to, for a barionic jet, about a 10% of the mass accreted by a black hole of $10^8 {\rm M_{\odot}}$ at the Eddington limit, pointing towards a very efficient conversion of accretion flow into ejection, or to a leptonic composition of jets.

Compact symmetric objects conform a class of sources characterized by high luminosity radio emission located symmetrically on both sides of the active galactic nucleus on linear scales of less than 1 kpc. Given their small size, the hot spots of the jets in CSOs provide a unique laboratory for the study of the physics of relativistic jets and their environment close to the central engine. We present a simple model for the hot spots in CSOs assuming synchrotron emission, minimum energy and ram pressure equilibrium with the external medium. Further comparison of our model with observational data allows us to constrain the physical parameters in the hot spots and the jets feeding them, and the density profile of the external medium.