No Arabic abstract
We present self-consistent global, steady-state MHD models and synthetic optically thin synchrotron emission maps for the jet of M87. The model consist of two distinct zones: an inner relativistic outflow, which we identify with the observed jet, and an outer cold disk-wind. While the former does not self-collimate efficiently due to its high effective inertia, the latter fulfills all the conditions for efficient collimation by the magneto-centrifugal mechanism. Given the right balance between the effective inertia of the inner flow and the collimation efficiency of the outer disk wind, the relativistic flow is magnetically confined into a well collimated beam and matches the measurements of the opening angle of M87 over several orders of magnitude in spatial extent. The synthetic synchrotron maps reproduce the morphological structure of the jet of M87, i.e. center-bright profiles near the core and limb-bright profiles away from the core. At the same time, they also show a local increase of brightness at some distance along the axis associated to a recollimation shock in the MHD model. Its location coincides with the position of the optical knot HST-1. In addition our best fitting model is consistent with a number of observational constraints such as the magnetic field in the knot HST-1, and the jet-to-counterjet brightness ratio.
(abridged) Significant progress has been made in the last years in the understanding of the jet formation mechanism through a combination of numerical simulations and analytical MHD models for outflows characterized by the symmetry of self-similarity. In a previous article we introduced models of truncated jets from disks, i.e. evolved in time numerical simulations based on a radially self-similar MHD solution, but including the effects of a finite radius of the jet-emitting disk and thus the outflow. These models need now to be compared with available observational data. A direct comparison of the results of combined analytical theoretical models and numerical simulations with observations has not been performed as yet. In order to compare our models with observed jet widths inferred from recent optical images taken with HST and AO observations, we use a new set of tools to create emission maps in different forbidden lines, from which we determine the jet width as the FWHM of the emission. It is shown that the untruncated analytical disk outflow solution considered here cannot fit the small jet widths inferred by observations of several jets. Various truncated disk-wind models are examined, whose extracted jet widths range from higher to lower values compared to the observations. Thus we can fit the observed range of jet widths by tuning our models. We conclude that truncation is necessary to reproduce the observed jet widths and our simulations limit the possible range of truncation radii. We infer that the truncation radius, which is the radius on the disk mid-plane where the jet-emitting disk switches to a standard disk, must be between around 0.1 up to about 1 AU in the observed sample for the considered disk-wind solution. One disk-wind simulation with an inner truncation radius at about 0.11 AU also shows potential for reproducing the observations, but a parameter study is needed.
(Abridged) We present maps for various Galactic longitudes and latitudes at 1.4 GHz, which is the frequency where deep SKA surveys are proposed. The maps are about 1.5 deg in size and have an angular resolution of about 1.6 arcsec. We analyse the maps in terms of their probability density functions (PDFs) and structure functions. Total intensity emission is more smooth in the plane than at high latitudes due to the different contributions from the regular and random magnetic field. The high latitude fields show more extended polarized emission and RM structures than those in the plane, where patchy emission structures on very small scales dominate. The RM PDFs in the plane are close to Gaussians, but clearly deviate from that at high latitudes. The RM structure functions show smaller amplitudes and steeper slopes towards high latitudes. These results emerge from the fact that much more turbulent cells are passed through by the line-of-sights in the plane. Although the simulated random magnetic field components distribute in 3D, the magnetic field spectrum extracted from the structure functions of RMs conforms to 2D in the plane and approaches 3D at high latitudes. This is partly related to the outer scale of the turbulent magnetic field, but mainly to the different lengths of the line-of-sights.
Analytical radially self-similar models are the best available solutions describing disk-winds but need several improvements. In a previous article, we introduced models of jets from truncated disks, i.e. evolved in time numerical simulations based on a radially self-similar MHD solution but including the effects of a finite radius of the jet-emitting disk and thus the outflow. In paper I of this series, we compared these models with available observational data varying the jet density and velocity, the mass of the protostar and the radius of the aforementioned truncation. In paper I, we assumed that the jet lies in the plane of the sky. In this paper, we investigate the effect of different inclinations of the jet. In order to compare our models with observed jet widths inferred from recent optical images taken with HST and AO, we create again emission maps in different forbidden lines and from such emission maps, we determine the jet width as the full-width half-maximum of the emission. We can reproduce the jet width of DG Tau and its variations very well and the derived inclination of 40$^circ$ is in excellent agreement with literature values of 32--52$^circ$. In CW Tau we overestimate the inclination in our best-fit model. In the other objects, we cannot find appropriate models which reproduce the variations of the observed jet widths, only the average jet width itself is well modeled as in paper I. We conclude that truncation -- i.e. taking into account the finite radius of the jet launching region -- is necessary to reproduce the observed jet widths and our simulations limit the possible range of truncation radii. The effects of inclination are important for modeling the intrinsic variations seen in observed jet widths. Our models can be used to infer independently the inclinations in the observed sample, however, a parameter study with a finer grid of parameters is needed.
M87 is one of the nearest radio galaxies with a prominent jet extending from sub-pc to kpc-scales. Because of its proximity and large mass of the central black hole, it is one of the best radio sources to study jet formation. We aim at studying the physical conditions near the jet base at projected separations from the BH of $sim7-100$ Schwarzschild radii ($R_{rm sch}$). Global mm-VLBI Array (GMVA) observations at 86 GHz ($lambda=3.5,$mm) provide an angular resolution of $sim50mu$as, which corresponds to a spatial resolution of only $7~R_{rm sch}$ and reach the small spatial scale. We use five GMVA data sets of M87 obtained during 2004--2015 and present new high angular resolution VLBI maps at 86GHz. In particular, we focus on the analysis of the brightness temperature, the jet ridge lines, and the jet to counter-jet ratio. The imaging reveals a parabolically expanding limb-brightened jet which emanates from a resolved VLBI core of $sim(8-13) R_{rm sch}$ size. The observed brightness temperature of the core at any epoch is $sim(1-3)times10^{10},$K, which is below the equipartition brightness temperature and suggests magnetic energy dominance at the jet base. We estimate the diameter of the jet at its base to be $sim5 R_{rm sch}$ assuming a self-similar jet structure. This suggests that the sheath of the jet may be anchored in the very inner portion of the accretion disk. The image stacking reveals faint emission at the center of the edge-brightened jet on sub-pc scales. We discuss its physical implication within the context of the spine-sheath structure of the jet.
In this paper three dimensional relativistic hydrodynamic simulations of AGN jets are presented to investigate the FR I/FR II dichotomy. Three simulations are presented which illustrates the difference in morphology for high/low Lorentz factor injection as well as a stratified background medium. Lorentz factors of 10 and 1.0014 were used for the high and low Lorentz factor cases respectively. The hydrodynamic simulations show a division in the morphology of jets based on their initial injection luminosity. An additional simulation was set up to investigate the evolution of the low Lorentz factor jet if the mass injection was lowered after a certain time. A synchrotron emission model was applied to these simulations to reproduce intensity maps at radio frequencies (1.5GHz) which were compared to the observed emission structures of FR I/FR II radio galaxies. The effect of Doppler boosting on the intensity maps was also investigated for different polar angles. The intensity maps of both the high and low Lorentz factor cases reproduced emission structures that resemble those of FR II type radio galaxies with a dominant cocoon region containing time dependent hot spots and filaments. An FR I like structure was, however, produced for the low Lorentz factor case if the mass injection rate was lowered after a set time period.