Issue 
A&A
Volume 609, January 2018



Article Number  A80  
Number of page(s)  16  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201731048  
Published online  12 January 2018 
Jettorus connection in radio galaxies
Relativistic hydrodynamics and synthetic emission
^{1} Institut für Theoretische Physik, Goethe Universität, MaxvonLaueStr. 1, 60438 Frankfurt, Germany
email: cfromm@th.physik.unifrankfurt.de
^{2} MaxPlanckInstitut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
^{3} Departament d’Astronomia i Astrofísica, Universitat de València, Dr. Moliner 50, 46100 Burjassot, València, Spain
^{4} Observatori Astronòmic, Parc Científic, Universitat de València, C/ Catedràtic José Beltrán 2, 46980 Paterna, València, Spain
^{5} Frankfurt Institute for Advanced Studies, RuthMoufangStrasse 1, 60438 Frankfurt, Germany
Received: 26 April 2017
Accepted: 28 August 2017
Context. 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 jets or they may be induced by the different conditions in the surrounding ambient medium, including the obscuring torus, or a combination of the two.
Aims. In this paper we investigate the influence of the ambient medium, including the obscuring torus, on the observed properties of jets from radio galaxies.
Methods. We performed specialrelativistic hydrodynamic (SRHD) simulations of overpressured and pressurematched jets using the specialrelativistic hydrodynamics code Ratpenat, which is based on a secondorder accurate finitevolume method and an approximate Riemann solver. Using a newly developed radiative transfer code to compute the electromagnetic radiation, we modelled several jets embedded in various ambient medium and torus configurations and subsequently computed the nonthermal emission produced by the jet and thermal absorption from the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory.
Results. The detailed analysis of our simulations shows that the observed properties such as core shift could be used to distinguish between overpressured and pressure matched jets. In addition to the properties of the jets, insights into the extent and density of the obscuring torus can be obtained from analyses of the singledish spectrum and spectral index maps.
Key words: galaxies: active / galaxies: jets / radiation mechanisms: nonthermal / radio continuum: galaxies
© ESO, 2018
1. Introduction
Active galactic nuclei (AGNs) are generally assumed to consist of a central, compact massive object (i.e. a supermassive black hole), an accretion disk surrounded by a molecular torus, and a jet launched perpendicular to the blackhole disk system. Within this standard model, AGNs may be classified according to their orientation with respect to the line of sight. If the AGN is seen at a viewing angle of ϑ = 90° (i.e. edgeon), the innermost regions are obscured by the torus and more internal features become observable as the viewing angle decreases (Antonucci 1993; Urry & Padovani 1995).
In addition to the classification of AGNs by their viewing angle, Fanaroff & Riley (1974) introduced a division based on the morphology of the jets: FanaroffRiley I (FR I) jets are characterised by a bright core and FR II jets are dominated by their bright lobes.
The largescale structure of AGNs observed at high viewing angles can be resolved through radio interferometric observations (see e.g. Kharb et al. 2008) and high resolution very long baseline interferometry (VLBI) observations enable us to probe the morphological evolution of AGNs on parsec scales (see e.g. Lister et al. 2009). These jets remain collimated over long distances, up to kiloparsec scales, at which they can be modelled using specialrelativistic (magneto)hydrodynamical [R(M)HD] simulations and the computed emission and polarisation properties provide insights into the underlying fluid properties and magnetic field configurations (see e.g. Perucho & Martí 2007; Laing & Bridle 2002; Hardcastle & Krause 2013).
In the infrared band, the electromagnetic emission from these AGNs is assumed to be generated by the reemission by the obscuring torus of the radiation initially produced in the accretion disk. Based on observations, two kinds of torus models were constructed, namely smooth torus models and clumpy torus models (for a review of torus models, see e.g. Hoenig 2013).
So far both models have been applied successfully to explain observational results either at kiloparsec scales (via jet simulations) or at subparsec scales (via torus simulations). However, to explain the VLBI observations of AGNs viewed almost edgeon , jet and torus models need to be combined. In this paper we combine both approaches, addressing the question of the impact that the combination of the torus and the ambient medium have on the observed properties of the jet at parsec scales.
Additionally, we aim to make the connection between numerical simulations of relativistic outflows and their emission with observations in a deeper way than has been carried out so far. Jet acceleration, recollimation, and other phenomena have previously been analysed with the help of onedimensional (1D), 2D, and 3D numerical simulations, including various characteristics such as magnetic fields and special relativity (see e.g. Mizuno et al. 2015). Some emission codes have even been applied (Gómez et al. 1997; Mimica et al. 2009) to explain the emission. Here we go one step further and consider additional components in the AGN picture, such as the obscuring torus, and including time delays, emission and absorption in a code, which will enable us to compare with real observations of twosided sources in future works.
The organisation of the paper is as follows. In Sect. 2 we introduce our numerical setup for the jet and torus. The results of the simulations and emission calculations are presented in Sects. 3.1 and 3.2. The discussion of our results is provided in Sect. 4. Throughout the paper we assume an idealfluid equation of state , where p is the pressure, ρ the restmass density, ϵ the specific internal energy, and the adiabatic index (Rezzolla & Zanotti 2013).
2. Specialrelativistic hydrodynamic simulations
We performed several 2D axisymmetric simulations of supersonic specialrelativistic hydrodynamical jets via the finitevolume code Ratpenat (for more details see Perucho et al. 2010, and references therein).
2.1. Jet model
The numerical setup is similar to that used in Fromm et al. (2016) and for completeness we provide below a short summary.
The numerical grid consisted of 320 cells in the radial direction and 400 cells in the axial direction. Using a numerical resolution of 4 cells per jet radius^{1} (R_{j}), the grid covered 80 R_{j} × 100 R_{j}. We assumed the z axis to be in the direction of propagation of the jet and xaxis to represent the radial cylindrical coordinate.
We used outflow conditions at z = z_{max} and R = R_{max} and, owing to the axisymmetry, we applied reflection conditions at R = 0. At z = 0 we used inflow conditions for r ≤ 1 R_{j}, i.e. we injected a velocity v_{b} parallel to the zaxis and outflow conditions for R>R_{j}. We initialised a preexisting jet along the simulation box surrounded by a dense ambient medium with decreasing pressure gradient. The jet is characterised by the following parameters at injection (z = 0): the velocity of the jet v_{b}, the classical Mach number of the jet M, the restmass density of the jet ρ_{b}, and the adiabatic index . The bulk Lorentz factor Γ was computed from the jet velocity and pressure, p_{b} at the jet nozzle from the Mach number using an idealfluid equation of state. The ratio d_{k} = p_{b}/p_{a} between the pressure in the jet, p_{b}, and the pressure in the ambient medium, p_{a}, led to either a pressure matched (PM) jet (i.e. with d_{k} = 1) or to an overpressured (OP) jet (i.e. with d_{k}> 1).
Depending on this ratio, internal structures in the jet, the socalled recollimation shocks were formed in the case of an OP jet or, in the case of a PM jet, the jet appeared featureless (Mizuno et al. 2015). The properties of the ambient medium (e.g. the gradient in the ambient pressure and ambient restmass density) played a crucial role in the shape of the jet and the properties of the created recollimation shocks in the case of OP jets. We modelled the decrease in the ambient medium pressure with a pressure profile presented in Gómez et al. (1997) as (1)where z_{c} can be considered as the core radius and the exponents n and m control the steepening of the ambient pressure. The initial conditions, given in code units, where speed of light c = 1, jet radius is R_{j}, and ambient medium restmass density, ρ_{a} = 1 g cm^{3}, are listed in Table 1.
Fig. 1
Distribution of ambient pressure and ambient restmass density in code units for the gradients in the ambient medium. 
Initial parameters for the jet.
The distributions of the ambient pressure and ambient restmass density for the parameters for the jet models considered here are shown in Fig. 1. The external pressure gradient is kept constant during the simulation by adding an external force, which preserves the pressure gradient without affecting the dynamics of the jet (see also Gómez et al. 1997). Throughout the paper we used a single value of n = 1.5, while setting m = 2, 3, 4 and excluding gravity and cooling processes during the SRHD simulations.
2.2. Torus model
Once a stationary state for the jet was reached, typically after five grid longitudinal crossing times, we inserted a steadystate torus. The included torus acts only as an absorber for the synchrotron radiation generated in the jet and has no influence on the dynamics of the jet. Therefore the torus should be regarded as a phenomenological model; Stalevski et al. (2016) provide a detailed discussion of the torus modelling.
Fig. 2
Geometry of the assumed torus model (adapted from Stalevski et al. 2012). 
More specifically, our torus model is adapted from Stalevski et al. (2012, 2016) and the geometry of the torus can be seen in Fig. 2. The geometry of the model resembles a flared disk and is characterised by three parameters: the inner radius of the torus, R_{in}, outer radius of the torus, R_{out}, and the halfopening angle Θ. For the inner radius of the torus we used (2)where L_{AGN} is the bolometric luminosity and T_{sub} is the sublimation temperature of the dust grains. We assumed a power law distribution for the restmass density in the torus, given by (3)where the cosΘdependence in the Eq. (3) leads to a concentration of restmass density towards the equatorial plane. Since we did not perform a radiative modelling of the torus, its temperature distribution is simply assumed to follow a certain prescription. In particular, following Schartmann et al. (2005), the temperature decrease moving outward in the radial direction; furthermore, because of the direct illumination of the inner torus walls, the temperature should decrease when going from high to low latitudes. Such behaviour may be modelled with the following temperature profile (4)Figure 3 illustrates the distribution of the temperature and restmass density in the torus with L_{AGN} = 10^{43} erg s^{1}, T_{sub} = 1500 K, which leads to an inner radius of R_{in} = 4 R_{j}. The outer radius is set to R_{out} = 30 R_{j} and the angular thickness of the torus is θ = 50°. For the scaling of the restmass density we set , while the exponents for the distribution of the temperature and restmass density are k = 1 and l = 2.
Fig. 3
Left panel: distribution of temperature in the meridional plane. Right panel: corresponding restmass density distribution. The Θdependence in the distribution of the temperature and the restmass density is clearly visible. 
Fig. 4
Example of the 3D geometry of the jettorus model. 
Fig. 5
Stationary results for the jet simulations. The panels show the 2D distribution of the restmass density for different ambient medium configurations as indicated by the exponent m. In each panel the left side (−30 <r/R_{j}< 0) corresponds to the OP jet and the right side (0 <y/R_{j}< 30) to the PM jet. The top row spans the entire simulation grid whereas the bottom row shows a magnified view of the nozzle region (−7 <z/R_{j}< 7). The white lines in the bottom row show stream lines visualising the direction of the flow and the bold dashed lines correspond to the inward travelling and reflected shock wave. 
The combined system consisting of jet, counterjet, and torus is shown in Fig. 4. The observed image is the projection of the jettorus system onto the x−y plane. In Fig. 4 a distant observer would see the jettorus system at a viewing angle of ϑ = 0°.
3. Results
3.1. Specialrelativistic hydrodynamics
In Fig. 5 we show the 2D distribution of the restmass density of the stationary jets from our (SRHD) simulations. Each panel corresponds to a different ambient medium configuration as characterised by the exponent m (see Eq. (1)). The left half of each panel (−30 <r/R_{j}< 0) shows the distribution for the OP jet and the right half of each panel (0 <r/R_{j}< 30) for the PM jet. The top row presents the entire simulation grid and the bottom row gives a magnified view of the nozzle region of the jets (0 <z/R_{j}< 7).
In the case of OP jets, the initial pressure mismatch at the nozzle led to the generation of two shock waves: one travelling radially outward and the other propagating inward in the direction of the jet axis. Between these two shocks the jet radially expanded until a pressure equilibrium with the ambient medium was reached. This equilibrium is first obtained at the jetambient medium boundary and leads to an inward travelling wave, which is indicated by the bold dashed lines in the bottom row of Fig. 5. The radial expansion of the jet is stopped as soon as the expansion regions are crossed by the inward travelling sound waves. Since these waves propagate with the local sound speed, the outer regions in the jet already stop expanding while the inner regions continue to expand (further decreasing the pressure and restmass density). The innermost region in the jet stops expanding as soon as the sound wave reaches the jet axis. At this position, a recollimation shock is formed, which can be regarded as a new jet nozzle (local maximum in pressure and restmass density) and so the flow expands again (see e.g. Daly & Marscher 1988; Falle 1991; Mizuno et al. 2015). The recollimation shocks are clearly visible at z ~ 6 R_{j} in the top half of the panels in the left column of Fig. 5. At the position of the recollimation shock, a new radially outward travelling wave is formed and the jet expands again.
If the pressure gradient in the ambient medium at the position of the first recollimation shock (z ~ 6 R_{j}) is steep enough, the sound waves reflected from the boundary between the jet and ambient medium (generated after reaching radial pressure equilibrium) cannot reach the jet axis and therefore no additional recollimation shock is formed and the jet continues to expand (for detailed discussion see Porth & Komissarov 2015). With our settings for the ambient medium, this behaviour was obtained for the models with m> 2. In these cases, even for the OP jets no additional recollimation shocks were formed beyond z_{c} ~ 6 R_{j} and a conical jet shape was obtained further downstream.
In contrast to the OP jets, in PM jets there is a radialpressure equilibrium at the nozzle. Therefore, only minor variations are created within the jet and the pressure and restmass density are continuously decreasing. If the pressure in the jet decreases faster than the pressure in the ambient medium, there is a position along the jet at which the pressure in the ambient is higher than in the jet. At this position a reconfinement shock is created that is travelling towards the jet axis (Komissarov & Falle 1997). This behaviour can be seen in case of PM jets embedded in an ambient medium with pressure gradient m = 4 at z ~ 60 R_{j} in Figs. 5 and 6. The opening angle of these jets depends on the gradient in the ambient medium: the steeper the gradient the larger the opening angle of the jet. This behaviour is visible in the first column of Fig. 5.
Fig. 7
Variation of the jet parameters along the jet axis. The left column corresponds to OP jets and the right column to PM jets. From top to bottom: the jet radius, axial pressure, and axial velocity in the zdirection are shown (see text for details). 
A detailed picture of the influence of the pressure ratio at the jet nozzle, d_{k}, and the configuration of the ambient medium on the shape of the jet and its fluid parameters can be obtained by analysing their evolution along the jet axis. The variation of the normalised jet radius, r/R_{j}, of the axial pressure, p_{axial}, and of the axial velocity, v_{axial}, for the different jet and ambient medium models is shown in Fig. 7. The largest variations in the jet parameters occur at small scales (i.e. for z ≲ 20 R_{j}) and for different jet models, i.e. pressure matched or overpressured. Within a jet model the modification of the gradient in the ambient medium leads to variations, but does not lead to a different evolution of the parameters along the jet axis.
In the case of PM jets, stronger gradients in the ambient medium lead to larger jet radii, r, faster velocities, v_{axial}, and smaller axial pressures, p_{axial}. This behaviour can be understood in the following way: the radial position where a pressure equilibrium between the jet and the ambient medium is reached (given the same initial pressure at the jet nozzle) shifts downstream with the gradient in the ambient medium. Since the jet is in pressure equilibrium with the ambient medium, the axial pressure in the jet is equal to the pressure in the ambient medium, so that the steeper the gradient in the ambient medium the stronger the decrease in the jet pressure. According to the conservation equations of hydrodynamics, a hot supersonic flow accelerates if its cross section increases (Rezzolla & Zanotti 2013); thus a PM jet that is embedded in an ambient medium with larger gradients (hence larger jet radii) develops higher velocities than the same jet that is surrounded by an ambient medium with smaller gradients (see bottom panel of Fig. 7).
In contrast to the smooth evolution of the jet parameters (r, p_{axial}, and v_{z,axial}) for the PM jet, large variations are present in the parameters for OP jets (see Fig. 7). As already mentioned the initial pressure mismatch at the jet nozzle led to an opening of the jet and as soon as a radial pressure equilibrium was reached a radially inward travelling wave was formed. The larger the gradient in the ambient medium, the larger the jet radii where the pressure equilibrium was obtained became. Because of the finite sound speed, the crossing time of this wave increased with jet radius. The inner parts of the jet stopped expanding as soon as the sound wave crossed the jet axis. Therefore, the pressure decreased for larger jet radii, i.e. a steeper gradient in the ambient medium (compare solid red and blue lines in the middle panel of Fig. 7).
After the sound wave crossed the jet axis, the pressure increases to values comparable with the pressure in the ambient medium (see middle panel of Fig. 7 at z ~ 8 R_{j}). The difference in position of the pressure jump and its magnitude reflects the difference in the sound crossing time and the gradient in the ambient medium. Since the OP jets exhibit larger jet radii than PM jets, the fluid was accelerated to higher velocities, v_{z,axial} at z ~ 5 R_{j}. The formation of the recollimation shock at z ~ 6 R_{j} leads to a sharp decrease in the axial velocity v_{z,axial}, and the expansion of the fluid after crossing the recollimation shock leads to an increase in the axial velocity. After z> 10 R_{j}, the rapid decrease in the ambientmedium pressure induces an expansion of the jet, which is accompanied by a decrease in the jet pressure, p_{axial}, and a slight increase in the axial jet velocity. In the case of the smallest gradients in the ambient medium (m = 2) additional recollimation shocks are formed at z ~ 20 and at z ~ 40 (see the solid red curve in Fig. 7).
3.2. Emission
For the calculation of the emission we followed the methods described in Gómez et al. (1997) and Mimica et al. (2009). For completeness we provide here the basic equations and underlying assumptions for the calculation of nonthermal emission. We assumed a power law distribution of relativistic electrons and included adiabatic, but neglected synchrotron cooling during the calculation of the nonthermal emission. Since the magnetic fields used in the computation of the emission were well below the equipartition ratio (ϵ_{b}< 0.2), we expected no or only a minor influence of radiative cooling on the nonthermal emission in the radio regime , (5)where n_{0} is a normalisation coefficient, γ is the electron Lorentz factor, γ_{min} and γ_{max} are the lower and upper electron Lorentz factors and s is the spectral slope. In order to construct the distribution of nonthermal particles from the thermal population we assumed that the number density of relativistic particles is a fraction ζ_{e} of the thermal particles (6)where m_{p} is the proton mass. The second assumption relates the energy in the relativistic particles to the energy in the thermal particles via the parameter ϵ_{e} as follows: (7)where m_{e} is the electron mass. After performing the integration in Eqs. (6) and (7) a relation for lower electron Lorentz factor, γ_{e, min}, may be derived as (8)For the upper electron Lorentz factor we used a constant ratio of the lower electron Lorentz factor, namely (9)The normalisation coefficient of the nonthermal particle distribution can be obtained by performing the integral in Eq. (6) within the boundaries given by γ_{e, min} and γ_{e, max}, yielding (10)Since the information on the magnetic field cannot be obtained from our purely hydrodynamical numerical simulations, we introduced it through its fraction ϵ_{B} of the equipartition magnetic field, namely, we set (11)Following Mimica et al. (2009), the emission and absorption coefficients for the nonthermal emission, respectively ϵ_{ν,nt} and α_{ν,nt}, may be written as^{2}where the characteristic frequency, ν_{0} = 3eBsinθ/ (4πm_{e}c) with θ the angle between the direction to the observer and the direction of the magnetic field, , and e is the electron charge.
The value of H(ξ) depends on the orientation of the magnetic field. If the magnetic field is ordered, H(ξ) = F(ξ), i.e. (14)where K_{5/3} is the modified Bessel function of the second kind of order 5/3; in the case of a random magnetic field H(ξ) = R(ξ), so that (15)with α the angle between the comoving magnetic field and the line of sight. For the functions F(ξ) and R(ξ) we used the approximations given by Crusius & Schlickeiser (1986) and Joshi & Böttcher (2011), namely In order to compute the observed emission, we transformed the coefficients of emission ϵ_{ν,nt} and absorption α_{ν,nt}, into the observer’s frame and corrected for cosmological effects where the Doppler factor is defined as D: = Γ^{1}^{(}1−v_{b}cosϑ^{)}^{1} and z denotes the redshift.
In addition to the nonthermal component, the emission also includes the contribution from a dusty torus, such that the thermal absorption coefficient α_{ν,th} is given by (20)where Z = μ_{i}/μ_{e}, with μ_{e} = 2/(1 + X); and μ_{i} = 4/(1 + 3X), where X corresponds to the abundance of hydrogen (set here to be X = 0.71); and is the velocityaveraged gaunt factor (see e.g. Rybicki & Lightman 1986). For the gaunt factor we used the values tabulated by van Hoof et al. (2014) and applied a bilinear interpolation in log u – log γ^{2} space, where u = (hν)/(k_{b}T) and γ^{2} = (Z^{2}Ry)/(k_{b}T), and Ry is the Rydberg energy.
The total intensity, I_{ν}, was a mixture of the nonthermal emission/absorption and thermal absorption (if the ray propagated through the dusty torus) and was computed via the transport equation as (21)where ds is the path length along a ray. For the calculation of the emission we computed the starting point of each ray at the position where the optical depth τ = τ_{limit} and followed the ray to the observer. We used a passive tracer advected during the SRHD simulations to distinguish between ambient medium (tracer = 0), jet (tracer ≈ 1) and torus (tracer = 2). The optical depth is defined as (22)The observed flux density can be computed from the intensity according to (23)where D_{L} is the luminosity distance and Δx and Δy represent the resolution of the detector’s frame. The luminosity distance is computed via the equation provided by Pen (1999), where we took Ω_{m} = 0.27 and a Hubble constant of H_{0} = 71 km s^{1} Mpc^{1}. The final jettorus geometry at a viewing angle of ϑ = 0° used for the calculation of the emission is shown in Fig. 4; furthermore, hereafter we assumed (D_{L} = 21.1 Mpc), which corresponds to a redshift of z = 0.005.
3.3. Parameter space study
In this section we analyse the influence of the emission and torus parameters on the observed spectrum. In Table 2 we provide an overview of the parameters used and of the values for our reference model. Generically, such parameters can be divided into three classes: i) scaling parameters; ii) emission parameters; and iii) torus parameters; we discuss these separately below.
Parameters for the emission simulations and values for the reference model.
In order to cover a wide frequency range, i.e. 10^{9} Hz <ν< 10^{12} Hz, we used a logarithmic grid with 100 frequency bins. The ray tracing was performed with a Delaunay triangulation on the axisymmetric RHD simulations with (r,z) coordinates to a 3D Cartesian grid (x,y,z) with a typical resolution of 300^{3} cells (see Fig. 4 for a 3D representation of the jettorus system).
3.3.1. Singledish spectra
To demonstrate the influence of the torus on the spectrum, we computed the emission for the reference model with and without a torus. The result of these calculations is presented in Fig. 8, where the black lines indicate simulations including a torus, while red lines are used for simulations without a torus. Furthermore, solid lines refer to the total flux density, while dashed and dotted lines correspond to the emission from the jet and counterjet, respectively.
Not surprisingly, in the lowfrequency regime, i.e. for ν < 10^{11} Hz, simulations with a torus show a smaller flux than simulations without a torus, simply as a result of the increased absorption along the line of sight. At higher frequencies, the flux densities of both simulations are very similar since the torus becomes optically thin at these frequencies. The comparison between the emission from the jet (dashed lines) and counterjet (dotted lines) shows that the counterjet is more affected by the absorption of the torus than the jet. Similar to the case of the total flux density, the emission from jet and counterjet including a torus converges to their counter parts without torus (at ν ~ 10^{11} Hz for the jet and at ν ~ 3 × 10^{11} Hz for the counterjet).
Fig. 8
Singledish spectra for a simulation including a torus (black) and without a torus (red). 
Fig. 9
Singledish spectra for the different jet and ambient medium models. 
Fig. 10
Influence of the scaling parameters R_{j}, ϑ, and L_{kin} on the shape of the observed singledish spectrum. 
Fig. 11
Influence of the emission parameters ρ_{a},ϵ_{e},ϵ_{b},ζ_{e},ϵ_{γ}, and s on the shape of the observed singledish spectrum. 
Figure 9 shows the singledish spectrum for the different jettorus and ambient medium models (see Table 1). Within the same ambient medium, the PM jets (d_{k} = 1; dashed lines in Fig. 9) exhibit higher flux densities than the OP jets (d_{k} = 2.5, solid lines in Fig. 9). Furthermore, OP jets show a flatter spectrum than PM jets within 10^{9} Hz ≤ ν ≤ 10^{11} Hz. The general trend could be described as follows: the steeper the gradient in the ambient medium, the flatter the spectrum. The change in the spectral slope at ν ~ 10^{11} Hz is indicative of the turnover frequency of the torus, i.e. the frequency at which the torus becomes optically thin (see also Fig. 8).
Figure 10 is meant to summarise the impact of the different scaling parameters, R_{j},ϑ, and,L_{kin}, on the singledish spectrum. More specifically, we computed the emission when only one of the scaling parameters is varied while the others are kept fixed to the values listed in Table 2. The variation in the observed flux can then be summarised as follows:

Jet radius, R_{j}: the jet radius acts as the scaling length of thesimulations and has the largest impact on the observedsingledish spectrum (see panel a in Fig. 10). Theincrease in flux with R_{j} follows a simple scaling , as expected from Eq. (23).

Viewing angle, ϑ: owing to Doppler boosting, the observed spectrum is shifted towards higher flux densities and higher turnover frequencies with decreasing viewing angle. Since we use large viewing angles for our simulations the variations in the spectrum are minor.

Bolometric luminosity, L_{kin}: changes in the bolometric luminosity lead to a variation in the inner radius of the torus (cf. Eq. (2)). Lowering the kinetic luminosity decreases the inner radius of the torus. Since we keep the outer radius of the torus fixed, the total size of the torus increases with decreasing kinetic luminosity and lower values for restmass density and temperature are obtained. Since the opacity of the torus scales as τ_{torus} ∝ T^{− 1/2}ρ^{2} (when assuming T and ρ to decay as r^{− k}), lower densities lead to lower opacities and a large emission from the jet. As a result, the observed flux density increases with decreasing kinetic luminosity.
In Fig. 11 we report the effect of the emission parameters ϵ_{e},ϵ_{b},ζ_{e},ϵ_{γ}, and s on the observed singledish spectrum while keeping the other parameters fixed. The emission simulations were performed using our reference model (see Table 2).

Ambient medium restmass density, ρ_{a}: a variation in the ambientmedium restmass densityρ_{a} leads to large changes in the shape and total flux density of the spectrum (see panel a in Fig. 11). The singledish spectrum steepens and is shifted to higher turnover flux densities with increasing ambient medium restmass density.
Fig. 12 Influence of the torus parameter on the observed singledish spectrum. Top row from left to right: the torus outer radius, torus thickness, and torus restmass density are shown. Bottom row from left to right: the torus temperature, distribution of the torus restmass density, and distribution of the torus temperature are shown.

Thermal to nonthermal energy ratio,ϵ_{e}: if more energy is stored in the nonthermal particles via an increase in ϵ_{e}, the turnover frequency increases, while the turnover flux density rises only slightly. The spectral index within ν ~ 10^{11} Hz decreases in absolute values, thus leading to a flatter spectrum (see panel b in Fig. 11).

Equipartition ratio, ϵ_{b}: an increase of the equipartition ratio ϵ_{b} leads to higher observed flux densities on the singledish spectrum and to a shift of the turnover position towards lower frequencies and higher flux densities (see panel c of Fig. 11). Additionally, the spectrum steepens for ν< 10^{11} Hz.

Ratio between thermal to nonthermal number densities, ζ_{e}: increasing the number of nonthermal particles leads to higher observed flux densities and shifted the turnover frequency towards lower frequencies. Additionally, the spectral slope increases (see panel d of Fig. 11).

Ratio between upper and lower e^{−}Lorentz factors, ϵ_{γ}: the variation of ϵ_{γ} had only a minor effect on the singledish spectrum (see panel e in Fig. 11). An increase in γ_{max} leads to a flattening of the spectrum at high frequencies (ν > 10^{11} Hz) and shifted the highfrequency cutoff in the spectrum towards higher frequencies, which are clearly visible for γ_{e} = 10^{2} (see blue line in panel e).

Spectral slope, s: changing the spectral slope of the nonthermal particle distribution induced only minor variations in the observed flux densities. However, the slope of the spectrum for ν ≤ 10^{11} Hz flattens, while for ν > 10^{11} Hz the spectral slope steepens with increasing s (see panel f in Fig. 11).
Since the properties of the obscuring torus strongly affect the observed emission, we also performed a parameterspace study varying the properties of the torus reported in Table 2. In our model the torus is characterised by its geometry (R_{out} and Θ) and by the distribution of the restmass density and temperature (see Eqs. (3), (4)). The first row in Fig. 12 summarises the changes in the singledish spectrum that result from the torus geometry, while the second row shows the modification in the spectral shape depending on the distribution of the restmass density and temperature. Overall, the variations in the spectrum can be summarised as

Torus geometry: the variation of the torus geometry and its effecton the singledish spectrum is presented in the top row ofFig. 12. From left to right we show the effects of thetorus outer radius, R_{out}, the torus thickness, Θ, and the torus restmass density, ρ_{t}. Increasing the outer radius of the torus leads to a drop in the observed emission between 10^{9} Hz ≤ ν ≤ 10^{11} Hz. The influence of the torus thickness is presented in panel b of Fig. 12, which shows that the observed flux density decreases as the thickness of the torus is increased.

Density distribution: the observed flux density decreases with increasing torus restmass density. At the same time, the turnover frequency is shifted towards higher frequencies (see panel d in Fig. 12).

Temperature distribution: variations in the temperature of the torus only have minor impacts on the observed flux density, so that the hotter the torus the higher the observed flux density (see panel g in Fig. 12).

Radial and angular distribution of the torus restmass density and torus temperature: the gradient in the radial direction mainly modifies the total observed flux density with a steeper gradient in the radial direction leading to an increase in the flux density (see panel e in Fig. 12). On the other hand, a larger effect on the shape of the spectrum comes from changing the gradient in the polar direction (see panel f). More specifically, if the temperature distribution is more concentrated at the edges and the restmass density falls off more rapidly towards the edges of the torus (i.e. if higher values for l are used), the spectrum steepens and the observed flux density is higher than for smoother restmass density and temperature distributions (i.e. for lower values for l).
All things considered, the parameter study on the singledish spectrum can be summarised as follows: the scaling parameter R_{j} and the ambient medium restmass density ρ_{a} have the largest impact on the observed singledish flux density (see panel a in Figs. 10 and 11). Changes in the emission parameters, ϵ_{e},ϵ_{b},ζ_{e},ϵ_{γ}, and s mainly modify the shape of the spectrum and lead only to minor variations in the flux density (see Fig. 11). The torus parameters show the smallest effect and mainly alter the slope in the optically thin regime of the spectrum (see Fig. 12).
VLBA properties used for the convolution of the radio maps.
Fig. 13
Convolved radio maps at 15 GHz. In each panel the range from −3 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 3 × 10^{18} cm corresponds to the counterjet. The panels correspond to different gradients used in the ambient medium indicated by n and m. In each panel the the top half corresponds to an OP jet and the bottom half to a PM jet. The convolving beam is plotted in bottom left corner of each panel. 
Fig. 14
Same as Fig. 13 at a frequency of 86 GHz. In each panel the range from −1.5 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 1.5 × 10^{18} cm corresponds to the counterjet. The panels correspond to different gradients used in the ambient medium indicated by n and m. In each panel the top half corresponds to an OP jet and the bottom half to a PM jet. The convolving beam is plotted in the bottom left corner of each panel. 
3.3.2. Synthetic radio maps
To facilitate comparison with VLBI observations, we computed radio maps at different frequencies and convolved these frequencies with a 2D Gaussian beam. The size of the beam depends on the properties of the array and on the observing frequency. In Table 3 we present the typical beam sizes (uniform weighting, fulltrack (u,v) coverage for a source at −8° declination) for the Very Long Baseline Array (VLBA) and its sensitivity, i.e. the lowest detectable thermal noise (see the EVN calculator^{3}, for details). If not stated otherwise, a fluxdensity cutoff of 5 × S_{limit} was used for the synthetic radio maps.
Fig. 15
Influence of the viewing angle (left column), the thermal to nonthermal energy ratio (middle column) and the torus restmass density (right column) on the 15 GHz radio maps. In each panel the range from −3 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 
In Fig. 13 we show the convolved radio maps at 15 GHz for different ambient medium and jet models (similar to Fig. 5). The panels show, from top to bottom, various models for the ambient medium and in each panel the top half, i.e. (0–1) × 10^{18} cm corresponds to the OP jet, and bottom half, i.e. (− 1−0) × 10^{18} cm, to the PM jet. The width of the jet increases with the gradient in the ambient medium, from top to bottom, from 0.5 × 10^{18} cm to 0.8 × 10^{18} cm, while the observed flux density decreases. Within a common atmosphere, the OP jets (top half of each panel) show a less elongated core region as compared to the PM jets (bottom half of each panel) and the axial flux density decreases faster in OP jets than in PM jets. These effects become more visible as the restmass density gradient in the ambient medium is increased.
Analogous to the parameter study performed on the singledish spectra (see Sect. 3.3.1), we investigated the impact of the emission parameters presented in Table 2 on the synthetic radio maps. Figure 15 summarises the effects of the viewing angle, ϑ, the thermal to nonthermal energy ratio ϵ_{e}, and the restmass density of the torus, ρ_{t}, on the 15 GHz radio maps. For a better comparison, the reference model is included in the second row of Fig. 15.
At smaller viewing angles (i.e. for ϑ = 60°), the jets appear more asymmetric because of the increased Doppler boosting and larger opacity along the line of sight, especially for the counterjet (see the top panel in the first row of Fig. 15). However, with increasing viewing angle the asymmetries between jet and counterjet tend to disappear, while the total flux density in the radio map decreases. Upon varying the nonthermal to thermal energy ratio ϵ_{e}, the flux density increases and affects both jets in the same way (see second row in Fig. 15). On the other hand, a change in ϵ_{e} does not alter the jet structure or the emission gap between the jet and counterjet.
In the third row of Fig. 15 we show the changes in the 15 GHz radio maps due to a variation in the torus restmass density, ρ_{t}. Since we keep the torus size constant, a change in the torus restmass density does not have an impact on the observed emission at x > 1 × 10^{18} cm. However, the observed gap between jet and counterjet does change depending on the torus restmass density. In particular, the gap is nearly invisible at low densities (top panel in the third row) while it is clearly noticeable at higher densities (see last panel in the third row of Fig. 15). Besides the changes in the gap between the jet and counterjet, a variation in the torus restmass densities leads to changes in the observed flux densities within x < 1 × 10^{18} cm. This can be easily explained as due to the increased opacity, which reduces the flux reaching the observer (see the decreasing flux density at x = 0.5 × 10^{18} cm).
In the analysis of VLBI observations it is common to compute spectral index maps, where the spectral index, α, between to adjacent frequencies ν_{1} and ν_{2} is given by (24)To avoid resolutioninduced artefacts in the spectral index maps, a common beam size was used for the convolution. Furthermore, as usual in observational VLBI to avoid overresolution from the highfrequency image, the beam size was set on the lowfrequency map.
In Fig. 16 we present the spectral index computed between 22 GHz and 43 GHz for different torus densities indicated in the top left corner. The spectral index (SPIX) is found to be negative (i.e. α ~ −0.6) for regions  x  > 0.5 × 10^{18} cm to reflect optically thin regions, and positive (i.e. α ~ 2.5) for  x  < 0.5 × 10^{18} cm to reflect optically thick jet regions.
These inner regions absorb a large amount of the emission because of the high restmass density in the surrounding torus. Indeed, as the torus restmass density is increased, more emission is found to be absorbed and the optically thick region expands towards the edges of the torus (compare the top to the bottom panels in Fig. 16). Besides the expansion of the optically thick region, the spectral index also generally increases for larger torus restmass densities. Thus, the extent of the optically thick region and the value of the spectral index can be used to deduce the size and restmass density of the obscuring torus.
Fig. 16
Spectral index maps computed between 22 GHz and 43 GHz for various torus densities. In each panel the range from −3 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 
Fig. 17
Synthetic multifrequency radio maps from 8 GHz to 86 GHz. In each panel the range from −3 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 
Multifrequency VLBI observations enabled us to investigate the opacity and therefore the frequencydependent position of the socalled “radio core”, i.e. of the first detectable fluxdensity maximum. The variation of this position with frequency provides details of the restmass density and the magnetic field in the source. This effect is called core shift (see Lobanov 1998, for details). In VLBI observations the initial position of the jet is lost during the imaging process, and various techniques are used to reconstruct the core shift between adjacent frequencies (see Fromm et al. 2013, for a summary). However, this is not needed when working with synthetic radio maps, since the initial position is known. Here, we defined the core as the position of the first fluxdensity maximum along an axial slice through the convolved radio map.
In Fig. 17 we present multifrequency VLBI radio maps from 8 GHz to 86 GHz based on our reference model (see Table 2 for the parameters used). The total flux density of the radio maps and the transversal size of the jets decreased with frequency, where the latter could be explained by the convolution with the observing beam (see Table 3). Given the parameters used for our reference model, the turnover frequency is at ~ 8 GHz, so that the flux density decreases as S_{ν} ∝ ν^{− (s−1)/2} for ν> 8 GHz. Since the flux density is lower at higher frequencies and the sensitivity of the array also decreases with frequency (see third column of Table 3), the jets appear shorter in size at higher frequencies. Furthermore, besides a drop in the observed emission with frequency, the gap between the jet and counterjet is reduced with frequency and the position of the fluxdensity maximum in the jet and the counterjet is shifted upstream.
In Fig. 18 we present the core shift obtained from the synthetic multifrequency radio maps shown in Fig. 17. The red stars correspond to the position of the radio core in the counterjet and the blue stars to the position of the radio core in the jet. The inset panel in Fig. 18 shows a magnified view of the nozzle region of the jets and the symbols indicate the position of the core. As mentioned previously, the core positions in the jet and counterjet shift upstream with frequency; at the same time, however, the amount of the shift differs between the jet and counterjet. This behaviour could be explained in terms of the various amounts of absorption along the line of sight for rays crossing the obscuring torus. The opacity is clearly larger for rays emerging from the counterjet since they have to cross the entire torus; this is in contrast to rays from the jet, which need to cross only a fraction of the torus. At ν ~ 80 GHz, the torus becomes optically thin and both cores are shifted towards the footpoint of the jets (here at x = 0).
Figure 18 also shows the core shifts obtained from an OP jet (diamond markers). The computed core shifts show the same behaviour as the PM jet for ν< 20 GHz. For higher frequencies, both the jet and counterjet show a plateau occurring between ν ~ 20 GHz and ν ~ 70 GHz. In this frequency range, the recollimation shocks in the jets mimic the radio core (i.e. they produce a local fluxdensity maximum) and since their position is frequency independent, the core shift does not decrease with frequency. On the other hand, if the observing frequency is increased, the upstream emission of the recollimation shocks becomes optically thin and a new fluxdensity maximum is produced further upstream. Similar to the PM jet, the torus becomes optically thin around ν ~ 80 GHz and the core positions converge towards the footpoint of the jet.
Fig. 18
Variation of the core shift with frequency as obtained from our reference model. The inset panel shows a contour plot of the synthetic radio maps overplotted with corresponding core position for different observing frequencies. In the inset panel the range from −1.5 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 1.5 × 10^{18} cm corresponds to the counterjet. 
4. Summary and outlook
We investigated the connection between jets in an ambient medium of decreasing pressure surrounded by a stationary torus. We then simulated jets using the stateofthe art SRHD code Ratpenat (Perucho et al. 2010) and included a stationary torus. The properties of the torus were motivated by both observations and theoretical modelling. In the postprocessing we used a newly developed emission code and computed the thermal and nonthermal emission. To investigate the effect of the scaling parameters on the observed emission (singledish spectra and multifrequency radio maps) we performed a detailed parameter study and produced observable signatures, which can be used to distinguish between different jet, ambient medium, and torus configurations based on singledish and VLBI observations. We applied a coreshift analysis to the synthetic radio maps produced and demonstrated how one can distinguish between OP jets and PM jets depending on the slope of the core shift. Our simulations and the extensive analysis of the parameter space study have shown that the observed emission is the result of a delicate interplay between hydrodynamical effects and emission and absorption processes.
A particularly robust result of this analysis is that the observed spectra of PM jets exhibit higher flux densities than the spectra of OP jets. This is because the initial pressure mismatch at the jet nozzle causes the OP jets to expand faster and convert internal energy into kinetic energy more efficiently than PM jets. At large viewing angles, the emission is strongly deboosted for faster jets and, as a result, less emission is received by the observer. Also, while the spectra alone make it difficult to distinguish these two types of jets, highfrequency radio maps show a clear difference in both the morphology of the jet and evolution of the flux density along the jet. More specifically, apart from the core region, OP jets exhibit local maxima in the flux density at the position of the recollimation shocks (see Fig. 14). On the other hand, in PM jets the flux density decreases continuously along the jet axis outside the core region Gómez et al. (1997), Mimica et al. (2009), Fromm et al. (2016).
Another important result of our analysis is that the inclusion of a stationary torus in the synthetic emission consistently leads to a strong absorption in the singledish spectrum. This behaviour is mostly due to the enhanced restmass density between the jet and observer, which in turn increases the opacity along the line of sight. Additionally, depending on properties of the torus, the spectrum flattens within a certain frequency range. Overall, the geometry of the torus has the greatest influence on the observed spectrum. In particular, increasing the geometrical size of the torus leads to a longer path for the rays generated in the jets, which then pass through a high restmass density field, thereby increasing the opacity along the line of sight. Indeed, if the torus is large enough and its restmass density sufficiently high, it can block a considerable fraction of the jet, which then becomes invisible to an observer. This effect is reflected in the flattening of the spectrum, in the appearance of an emission gap between the jet and counterjet in the synthetic radio maps, and in optically thick spectralindex maps. Since the thermal absorption coefficient scales as ν^{3}, the aforementioned effects are frequency dependent, i.e. the spectra including a torus should converge to that without a torus for ν> 10^{11} Hz, with the emission gap between the jet and counterjet vanishing with increasing frequency. Observationally, therefore, the spectral break could be used to constrain the restmass density within the torus. Given a dense frequency sampling, such as that used in the FGamma programme (Fuhrmann et al. 2016), the restmass density of the torus, together with additional emission parameters, can in principle be fitted.
This behaviour and in particular the frequency dependence of the synthetic radio maps can be used to distinguish between OP and PM jets. In particular, since the gap is frequency dependent, the observed onset of the jet and counterjet varies with frequency. At low frequencies the torus should block the observer’s view and the onset of the jets are separated by an emission gap; at such frequencies, OP and PM jets are indistinguishable. However, at higher observing frequencies the gap between the jets decreases, exposing the differences in the two types of jets. Since one of the key features of OP jets are recollimation shocks, which lead to local maxima in the pressure and restmass density along the jet, their presence becomes visible as a local enhancement of the emission. At a certain frequency, the onset of the jet coincides with the position of a recollimation shock, which itself is frequency independent and does not change until the upstream region becomes optically thinner and brighter than the recollimation shock. This effect leads to the formation of plateaus in the coreshift plot. Because the behaviour described above is absent in the case of PM jets, for which no strong recollimation shocks are formed and no plateaus in the core shift occur, observations should be able to provide information to discriminate between these two basic classes of jets. More specifically, using multifrequency VLBI observations, the size of the torus can be estimated by the dimensions of the optically thick region in the spectral index maps. The frequencydependent size of the emission gap can be used to obtain estimates of the restmass density of the torus and analysis of the core shift would enable one to differentiate between OP and PM jets.
The simulations performed in this work were able reproduce some of the observed properties of jets from radio galaxies. However, a more detailed treatment of the torus, in particular relaxing the assumption of homogeneity, i.e. using a clumpy torus could lead to more complex singledish spectra and radio maps. Including an independent treatment of the magnetic field, i.e. switching from SRHD to SRMHD simulations would enable us to test the influence of higher magnetisation on the nonthermal emission. Owing to the 2D axisymmetric SRHD setup, the jets are very stable against hydrodynamical instabilities. Therefore, future studies should include 3D SRMHD in order to test the stability of the jets against hydrodynamical and magnetohydrodynamical instabilities.
In a followup paper we will address some of the abovementioned simplifications and apply the techniques presented here to actual observational data and include the propagation of shock waves using the socalled “slowlight” approach, which will enable us to investigate the temporal variations in the singledish spectra and multifrequency radio maps.
Acknowledgments
Support comes from the ERC Synergy Grant “BlackHoleCam – Imaging the Event Horizon of Black Holes” (Grant 610058). M.P. acknowledges support by the Spanish “Ministerio de Economía y Competitividad” grants AYA201348226C32P and AYA201566899P. Z.Y. acknowledges support from an Alexander von Humboldt Fellowship. E.R. acknowledges support from the grants AYA201238491C0201, AYA201563939C22P, and PROMETEOII/2014/057. C.M.F. wants to thank Walter Alef and Helge Rottmann for support during the computation on the MPIfR cluster and Mariafelicia De Laurentis and Hector Olivares for fruitful discussion and helpful comments to the manuscript.
References
 Antonucci, R. 1993, ARA&A, 31, 473 [Google Scholar]
 Crusius, A., & Schlickeiser, R. 1986, A&A, 164, L16 [NASA ADS] [Google Scholar]
 Daly, R. A., & Marscher, A. P. 1988, ApJ, 334, 539 [NASA ADS] [CrossRef] [Google Scholar]
 Falle, S. A. E. G. 1991, MNRAS, 250, 581 [NASA ADS] [CrossRef] [Google Scholar]
 Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31 [NASA ADS] [CrossRef] [Google Scholar]
 Fromm, C. M., Ros, E., Perucho, M., et al. 2013, A&A, 557, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fromm, C. M., Perucho, M., Mimica, P., & Ros, E. 2016, A&A, 588, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fuhrmann, L., Angelakis, E., Zensus, J. A., et al. 2016, A&A, 596, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gómez, J. L., Marti, J. M., Marscher, A. P., Ibanez, J. M., & Alberdi, A. 1997, ApJ, 482, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Hardcastle, M. J., & Krause, M. G. H. 2013, MNRAS, 430, 174 [NASA ADS] [CrossRef] [Google Scholar]
 Hoenig, S. F. 2013, ArXiv eprints [arXiv:1301.1349] [Google Scholar]
 Joshi, M., & Böttcher, M. 2011, ApJ, 727, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Kharb, P., O’Dea, C. P., Baum, S. A., et al. 2008, ApJS, 174, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. S., & Falle, S. A. E. G. 1997, MNRAS, 288, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Laing, R. A., & Bridle, A. H. 2002, MNRAS, 336, 328 [NASA ADS] [CrossRef] [Google Scholar]
 Lister, M. L., Aller, H. D., Aller, M. F., et al. 2009, AJ, 137, 3718 [NASA ADS] [CrossRef] [Google Scholar]
 Lobanov, A. P. 1998, A&A, 330, 79 [NASA ADS] [Google Scholar]
 Mimica, P., Aloy, M. A., Agudo, I., et al. 2009, ApJ, 696, 1142 [NASA ADS] [CrossRef] [Google Scholar]
 Mizuno, Y., Gómez, J. L., Nishikawa, K.I., et al. 2015, ApJ, 809, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Pen, U.L. 1999, ApJS, 120, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Perucho, M., & Martí, J. M. 2007, MNRAS, 382, 526 [NASA ADS] [CrossRef] [Google Scholar]
 Perucho, M., Martí, J. M., Cela, J. M., et al. 2010, A&A, 519, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Porth, O., & Komissarov, S. S. 2015, MNRAS, 452, 1089 [NASA ADS] [CrossRef] [Google Scholar]
 Rezzolla, L., & Zanotti, O. 2013, Relativistic Hydrodynamics (Oxford University Press) [Google Scholar]
 Rybicki, G. B., & Lightman, A. P. 1986, Rad. Process. Astrophys., 400 [Google Scholar]
 Schartmann, M., Meisenheimer, K., Camenzind, M., Wolf, S., & Henning, T. 2005, A&A, 437, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stalevski, M., Fritz, J., Baes, M., Nakos, T., & Popović, L. Č. 2012, MNRAS, 420, 2756 [NASA ADS] [CrossRef] [Google Scholar]
 Stalevski, M., Ricci, C., Ueda, Y., et al. 2016, MNRAS, 458, 2288 [NASA ADS] [CrossRef] [Google Scholar]
 Sutherland, R. S., & Bicknell, G. V. 2007, ApJS, 173, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef] [Google Scholar]
 van Hoof, P. A. M., Williams, R. J. R., Volk, K., et al. 2014, MNRAS, 444, 420 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A:
We performed a convergence test for our emission calculations and the results are presented in Fig. A.1. For the test we used the reference model (see Table 2) and increased the number of cells used in the computation in 100^{3} steps. In Fig. A.1 we show the singledish spectrum between 1 GHz and 1 THz for different resolutions and the inset panel shows the variation of the total flux density with respect to the number of simulation grid cells. The calculated flux density converges to a common value and no significant differences was obtained above 400^{3} cells. At the lowest resolution (100^{3}) the highfrequency flux density is significantly different compared to the high resolution runs. This difference occurred since the highfrequency emission was created close to the jet nozzle and therefore emanates from a region of small transversal jet size. If the resolution was not sufficiently high to resolve these regions, the triangulation and interpolation of these cells onto the 3D grid failed (overestimating the contribution from the ambient medium) and led to unphysically low flux densities.
Fig. A.1
Convergence test for the emission calculations. Singledish spectrum for seven different resolutions. 
Similar to the test on emission calculations we conducted a convergence study of the SRHD simulations. For this study we employed the procedure described in Sect. 2, increasing the numerical resolution by a factor 2, resulting in 4, 8, and 16 cells per jet radii. We simulated both, a PM and OP jet using the initial conditions presented in Table 2. A steady state was obtained and after roughly after 5 to 6 longitudinal grid crossing times, we computed the emission with a fixed detector frame resolution of 300^{3} cells. In Fig. A.2 we present the result of the convergence study on the SRHD simulations. The solid lines correspond to the PM jet (d_{k} = 1) and the dashed lines to the OP jet (d_{k} = 2.5). The result of this study shows that our initial resolution of 4 cells per jet radii can reproduce 95% of the total flux density obtained from the simulation with the highest resolution in the case of PM jets and 90% in the case of OP jets. Since we are mainly interested in emission from the largescale structure, 4 cells per jet radii represents a sufficient resolution for our study. However, a higher resolution would be needed to sufficiently resolve the structure induced by travelling shocks and study the time variability of nonthermal emission from the jets.
Fig. A.2
Convergence test for the SRHD calculations. Singledish spectrum for three different resolutions for both PM jet (solid lines) and OP jet (dashed lines). 
All Tables
All Figures
Fig. 1
Distribution of ambient pressure and ambient restmass density in code units for the gradients in the ambient medium. 

In the text 
Fig. 2
Geometry of the assumed torus model (adapted from Stalevski et al. 2012). 

In the text 
Fig. 3
Left panel: distribution of temperature in the meridional plane. Right panel: corresponding restmass density distribution. The Θdependence in the distribution of the temperature and the restmass density is clearly visible. 

In the text 
Fig. 4
Example of the 3D geometry of the jettorus model. 

In the text 
Fig. 5
Stationary results for the jet simulations. The panels show the 2D distribution of the restmass density for different ambient medium configurations as indicated by the exponent m. In each panel the left side (−30 <r/R_{j}< 0) corresponds to the OP jet and the right side (0 <y/R_{j}< 30) to the PM jet. The top row spans the entire simulation grid whereas the bottom row shows a magnified view of the nozzle region (−7 <z/R_{j}< 7). The white lines in the bottom row show stream lines visualising the direction of the flow and the bold dashed lines correspond to the inward travelling and reflected shock wave. 

In the text 
Fig. 6
Same as Fig. 5 for the pressure. 

In the text 
Fig. 7
Variation of the jet parameters along the jet axis. The left column corresponds to OP jets and the right column to PM jets. From top to bottom: the jet radius, axial pressure, and axial velocity in the zdirection are shown (see text for details). 

In the text 
Fig. 8
Singledish spectra for a simulation including a torus (black) and without a torus (red). 

In the text 
Fig. 9
Singledish spectra for the different jet and ambient medium models. 

In the text 
Fig. 10
Influence of the scaling parameters R_{j}, ϑ, and L_{kin} on the shape of the observed singledish spectrum. 

In the text 
Fig. 11
Influence of the emission parameters ρ_{a},ϵ_{e},ϵ_{b},ζ_{e},ϵ_{γ}, and s on the shape of the observed singledish spectrum. 

In the text 
Fig. 12
Influence of the torus parameter on the observed singledish spectrum. Top row from left to right: the torus outer radius, torus thickness, and torus restmass density are shown. Bottom row from left to right: the torus temperature, distribution of the torus restmass density, and distribution of the torus temperature are shown. 

In the text 
Fig. 13
Convolved radio maps at 15 GHz. In each panel the range from −3 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 3 × 10^{18} cm corresponds to the counterjet. The panels correspond to different gradients used in the ambient medium indicated by n and m. In each panel the the top half corresponds to an OP jet and the bottom half to a PM jet. The convolving beam is plotted in bottom left corner of each panel. 

In the text 
Fig. 14
Same as Fig. 13 at a frequency of 86 GHz. In each panel the range from −1.5 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 1.5 × 10^{18} cm corresponds to the counterjet. The panels correspond to different gradients used in the ambient medium indicated by n and m. In each panel the top half corresponds to an OP jet and the bottom half to a PM jet. The convolving beam is plotted in the bottom left corner of each panel. 

In the text 
Fig. 15
Influence of the viewing angle (left column), the thermal to nonthermal energy ratio (middle column) and the torus restmass density (right column) on the 15 GHz radio maps. In each panel the range from −3 × 10^{18} cm <x< 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas <x< 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 

In the text 
Fig. 16
Spectral index maps computed between 22 GHz and 43 GHz for various torus densities. In each panel the range from −3 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 

In the text 
Fig. 17
Synthetic multifrequency radio maps from 8 GHz to 86 GHz. In each panel the range from −3 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 3 × 10^{18} cm corresponds to the counterjet. The convolving beam is plotted in the bottom left corner of each panel. 

In the text 
Fig. 18
Variation of the core shift with frequency as obtained from our reference model. The inset panel shows a contour plot of the synthetic radio maps overplotted with corresponding core position for different observing frequencies. In the inset panel the range from −1.5 × 10^{18} cm < x < 0 cm indicates the jet (pointing towards the observer) and the range from 0 mas < x < 1.5 × 10^{18} cm corresponds to the counterjet. 

In the text 
Fig. A.1
Convergence test for the emission calculations. Singledish spectrum for seven different resolutions. 

In the text 
Fig. A.2
Convergence test for the SRHD calculations. Singledish spectrum for three different resolutions for both PM jet (solid lines) and OP jet (dashed lines). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.