S. Massaglia

PLUTO: A Numerical Code for Computational Astrophysics

We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2, and 3 spatial dimensions and different systems of coordinates. The code provides a multiphysics, multialgorithm modular environment particularly oriented toward the treatment of astrophysical flows in presence of discontinuities. Different hydrodynamic modules and algorithms may be independently selected to properly describe Newtonian, relativistic, MHD, or relativistic MHD fluids. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Although a plethora of numerical methods has been successfully developed over the past two decades, the vast majority shares a common discretization recipe, involving three general steps: a piecewise polynomial reconstruction followed by the solution of Riemann problems at zone interfaces and a final evolution stage. We have checked and validated the code against several benchmarks available in literature. Test problems in 1, 2, and 3 dimensions are discussed.

MHD simulations of jet acceleration from Keplerian accretion disks. The effects of disk resistivity

Context: Accretion disks and astrophysical jets are used to model many active astrophysical objects, such as young stars, relativistic stars, and active galactic nuclei. However, existing proposals for how these structures may transfer angular momentum and energy from disks to jets through viscous or magnetic torques do not yet provide a full understanding of the physical mechanisms involved. Thus, global stationary solutions have not explained the stability of these structures; and global numerical simulations that include both the disk and jet physics have so far been limited to relatively short time scales and narrow (and possibly astrophysically unlikely) ranges of viscosity and resistivity parameters that may be crucial to defining the coupling of the inflow-outflow dynamics. Aims: We present self-consistent, time-dependent simulations of supersonic jets launched from magnetized accretion disks, using high-resolution numerical techniques. In particular we study the effects of the disks magnetic resistivity, parametrized through an α-prescription, in determining the properties of the inflow-outflow system. Moreover we analyze under which conditions steady state solutions of the type proposed in the self-similar models of Blandford & Payne can be reached and maintained in a self-consistent nonlinear stage. Methods: We used the resistive MHD FLASH code with adaptive mesh refinement (AMR), allowing us to follow the evolution of the structure on a long enough time scale to reach steady state. A detailed analysis of the initial configuration state is given. Results: We obtain the expected solutions within the axisymmetric (2.5 D) limit. Assuming a magnetic field around equipartition with the thermal pressure of the disk, we show how the characteristics of the disk-jet system, such as the ejection efficiency and the energetics, are affected by the anomalous resistivity acting inside the disk.

High‐resolution 3D relativistic MHD simulations of jets

Relativistic magnetized jets are key elements in active galactic nuclei and in other astrophysical environments. Their structure and evolution involve a complex non-linear physics that can be approached by numerical studies only. Still, owing to a number of challenging computational aspects, only a few numerical investigations have been undertaken so far. In this paper, we present high-resolution three-dimensional numerical simulations of relativistic magnetized jets carrying an initially toroidal magnetic field. The presence of a substantial toroidal component of the field is nowadays commonly invoked and held responsible for the process of jet acceleration and collimation. We find that the typical nose cone structures, commonly observed in axisymmetric two-dimensional simulations, are not produced in the three-dimensional case. Rather, the toroidal field gives rise to strong current-driven kink instabilities leading to jet wiggling. However, it appears to be able to maintain a highly relativistic spine along its full length. By comparing low- and high-resolution simulations, we emphasize the impact of resolution on the jet dynamical properties.

Formation of dynamical structures in relativistic jets: the FRI case

Context. Strong observational evidence indicates that all extragal actic jets associated with AGNs move at relativistic speed up to 100 pc - 1 kpc scales from the nucleus. At larger distances, reflecting the Fanaro ff-Riley radio source classification, we observe an abrupt dec eleration in FR-I jets while relativistic motions persist up to Mpc scale in FR -II. Moreover, VLBI observations of some object like B2 1144+35, Mrk501 and M87 show limb brightening of the jet radio emission at the parsec scale. This effect is interpreted kinematically as due to the presence of a d eboosted central spine at high Lorentz factor and of a weakly relativi stic external layer. Aims. In this paper we investigate whether these effects can be interpreted by a breaking of the collimated flow by external medium entrainment favored by sh ear instabilities, namely KelvinHelmholtz instabilities. We examine in details the physical conditions under which significant deceleration of a relativistic flow is produced. Methods. We investigate the phenomenon by means of high-resolution three-dimensional relativistic hydrodynamic simulations using the PLUTO code for computational astrophysics. Results. We find that the parameter of utmost importance in determinin g the instability evolution and the entrainment properties is the ambient/jet density contrast. We show that lighter jets suffer stronger slowing down in the external layer than in the central part and conserve a central spine at high Lorentz factor. Conclusions. Our model is verified by constructing synthetic emission map s from the numerical simulations that compare reasonably well with VLBI observations of the inner part of FR-I sources.

Heating groups and clusters of galaxies: The role of AGN jets

X-ray observations of groups and clusters of galaxies show that the Intra-Cluster Medium (ICM) in their cores is hotter than expected from cosmological numerical simulations of cluster formation which include star formation, radiative cooling and SN feedback. We investigate the effect of the injection of supersonic AGN jets into the ICM using axisymmetric hydrodynamical numerical simulations. A simple model for the ICM, describing the radial properties of gas and the gravitational potential in cosmological N-Body+SPH simulations of one cluster and three groups of galaxies at redshift z = 0, is obtained and used to set the environment in which the jets are injected. We varied the kinetic power of the jet and the emission-weighted X-ray temperature of the ICM. The jets transfer their energy to the ICM mainly by the effects of their terminal shocks. A high fraction of the injected energy can be deposited through irreversible processes in the cluster gas, up to 75% in our simulations. We show how one single, powerful jet can reconcile the predicted X-Ray properties of small groups, e.g. the LX − TX relation, with observations. We argue that the interaction between AGN jets and galaxy groups and cluster atmospheres is a viable feedback mechanism.

On the origin of X-shaped radio-sources: New insights from the properties of their host galaxies

A significant fraction of extended radio sources presents a peculiar X-shaped radio morphology: in addition to the classical double lobed structure, radio emission is also observed along a second axis of symmetry in the form of diffuse wings or tails. We re-examine the origin of these extensions relating the radio morphology to the properties of their host galaxies. The orientation of the wings shows a striking connection with the structure of the host galaxy as they are preferentially aligned with its minor axis. Furthermore, wings are only observed in galaxies of high projected ellipticity. Hydrodynamical simulations of the radio-source evolution show that X-shaped radio-sources naturally form in this geometrical situation: as a jet propagates in a non-spherical gas distribution, the cocoon surrounding the radio-jets expands laterally at a high rate producing wings of radio emission, in a way that is reminiscent of the twin-exhaust model for radio-sources.

The curved X-ray spectrum of PSR B1509 58 observed with BeppoSAX

We report timing and spectral results for PSR B1509 58 observed by BeppoSAX in February 1998. We obtained pulse proles with high statistical signicance from 0.1 to 300 keV that exhibit the well-known broad asymmetric single pulse. The shape of this pulse does not change across this energy window and can be described as the sum of a narrow and a broader Gaussian component separated0:13 in phase. The spectral distribution can be accurately represented, over the entire BeppoSAX energy range, by a curved function rather than a simple power-law shape. The extrapolation of this model to higher energies is also consistent with the COMPTEL fluxes in the 0.75{30 MeV range, with the maximum luminosity for the broad pulse of PSR B1509 58 reached at 5:4 MeV. The comparison with the Crab pulsar spectrum suggests a possible origin of the X-ray emission in terms of synchrotron radiation from secondary pair particles near the neutron star.

Two-component jet simulations I. Topological stability of analytical MHD outflow solutions

Context. Observations of collimated outflows in young stellar objects indicate that several features of the jets can be understood by adopting the picture of a two-component outflow, wherein a central stellar component around the jet axis is surrounded by an extended disk wind. The precise contribution of each component may depend on the intrinsic physical properties of the YSO-disk system as well as its evolutionary stage. Aims. This article reports a systematic separate investigation of these jet components via time-dependent simulations of two prototypical and complementary analytical solutions, each closely related to the properties of stellar outflows and disk winds. These models describe a meridionally and a radially self-similar exact solution of the steady-state, ideal hydromagnetic equations, respectively. Methods. Using the PLUTO code to carry out the simulations, the study focuses on the topological stability of each of the two analytical solutions, which are successfully extended to all space by removing their singularities. In addition, their behavior and robustness over several physical and numerical modifications is extensively examined. Therefore, this work serves as the starting point for the analysis of the two-component jet simulations. Results. It is found that radially self-similar solutions (disk winds) always reach a final steady-state while maintaining all their welldefined properties. The different ways to replace the singular part of the solution around the symmetry axis, being a first approximation towards a two-component outflow, lead to the appearance of a shock at the super-fast domain corresponding to the fast magnetosonic separatrix surface. These conclusions hold true independently of the numerical modifications and/or evolutionary constraints that the models have undergone, such as starting with a sub-modified-fast initial solution or different types of heating/cooling assumptions. Furthermore, the final outcome of the simulations remains close enough to the initial analytical configurations, thus showing their topological stability. Conversely, the asymptotic configuration and the stability of meridionally self-similar models (stellar winds) is related to the heating processes at the base of the wind. If the heating is modified by assuming a polytropic relation between density and pressure, a turbulent evolution is found. On the other hand, adiabatic conditions lead to the replacement of the outflow by an almost static atmosphere.

X-ray emission from expanding cocoons

X-ray observations of extragalactic radiosources show strong evidences of interaction between the radio emitting plasma and the X-ray emitting ambient gas. In this paper we perform a detailed study of this interaction by numerical sim- ulations. We study the propagation of an axisymmetric supersonic jet in an isothermal King atmosphere and we analyze the evolution of the resulting X-ray properties and their dependence on the jet physical parameters. We show the existence of two distinct and observationally subsequent dierent regimes of interaction, with strong and weak shocks. In the first case shells of enhanced X-ray emission are to be expected, while in the second case we expect deficit of X-ray emission coincident with the cocoon. By a comparison between analytical models and the results of our numerical simulations, we discuss the dependence of the transition between these two regimes on the jet parameters and we find that the mean controlling quantity results to be the jet kinetic power. We then discuss how the observed jets can be used to constrain the jet properties.

Two-component jet simulations - II. Combining analytical disk and stellar MHD outflow solutions

Context. Theoretical arguments along with observational data of YSO jets suggest the presence of two steady components: a disk wind type outflow needed to explain the observed high mass loss rates and a stellar wind type outflow probably accounting for the observed stellar spin down. Each component’s contribution depends on the intrinsic physical properties of the YSO-disk system and its evolutionary stage. Aims. The main goal of this paper is to understand some of the basic features of the evolution, interaction and co-existence of the two jet components over a parameter space and when time variability is enforced. Methods. Having studied separately the numerical evolution of each type of the complementary disk and stellar analytical wind solutions in Paper I of this series, we proceed here to mix together the two models inside the computational box. The evolution in time is performed with the PLUTO code, investigating the dynamics of the two-component jets, the modifications each solution undergoes and the potential steady state reached. Results. The co-evolution of the two components, indeed, results in final steady state configurations with the disk wind effectively collimating the inner stellar component. The final outcome stays close to the initial solutions, supporting the validity of the analytical studies. Moreover, a weak shock forms, disconnecting the launching region of both outflows with the propagation domain of the twocomponent jet. On the other hand, several cases are being investigated to identify the role of each two-component jet parameter. Time variability is not found to considerably affect the dynamics, thus making all the conclusions robust. However, the flow fluctuations generate shocks, whose large scale structures have a strong resemblance to observed YSO jet knots. Conclusions. Analytical disk and stellar solutions, even sub modified fast ones, provide a solid foundation to construct twocomponent jet models. Tuning their physical properties along with the two-component jet parameters allows a broad class of realistic scenarios to be addressed. The applied flow variability provides very promising perspectives for the comparison of the models with observations.