Radio sources of low power are the most common in the universe. Their jets typically move at nonrelativistic velocity and show plume-like morphologies that in many instances appear distorted and bent. We investigate the role of magnetic field on the propagation and evolution of low-power jets and the connection between the field intensity and the resulting morphology. The problem is addressed by means of three-dimensional magnetohydrodynamic (MHD) simulations. We consider supersonic jets that propagate in a stratified medium. The ambient temperature increases with distance from the jet origin maintaining constant pressure. Jets with low magnetization show an enhanced collimation at small distances with respect to hydrodynamic (HD) cases studied in a previous paper. These jets eventually evolve in a way similar to the HD cases. Jets with higher magnetization are affected by strong nonaxisymmetric modes that lead to the sudden jet energy release. From there on, distorted plumes of jet material move at subsonic velocities. This transition is associated with the formation of structures reminiscent of the `warm spots' observed in wide-angle-tail (WAT) sources.

We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through an Eulerian grid where the (relativistic) MHD equations are solved concurrently. Lagrangian particles follow fluid streamlines and represent ensembles of (real) relativistic particles with a finite energy distribution. The spectral distribution of each particle is updated in time by solving the relativistic cosmic ray transport equation based on local fluid conditions. This enables us to account for a number of physical processes, such as adiabatic expansion, synchrotron and inverse Compton emission. An accurate semi-analytically numerical scheme that combines the method of characteristics with a Lagrangian discretization in the energy coordinate is described. In presence of (relativistic) magnetized shocks, a novel approach to consistently model particle energization due to diffusive shock acceleration has been presented. Our approach relies on a refined shock-detection algorithm and updates the particle energy distribution based on the shock compression ratio, magnetic field orientation and amount of (parameterized) turbulence. The evolved distribution from each \textit{Lagrangian} particle is further used to produce observational signatures like emission maps and polarization signals accounting for proper relativistic corrections. We further demonstrate the validity of this hybrid framework using standard numerical benchmarks and evaluate the applicability of such a tool to study high energy emission from extra-galactic jets.

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 simmetry in the form of diffuse wings or tails. In a previous investigation we showed the existence of a connection between the radio morphology and the properties of the host galaxies. Motivated by this connection we performed two-dimensional numerical simulations showing that X-shaped radio sources may naturally form as a jet propagates along the major axis a highly elliptical density distribution, because of the fast expansion of the cocoon along the minor axis of the distribution. We intend to extend our analysis by performing three-dimensional numerical simulations and investigating the role of different parameters is determining the formation of the X-shaped morphology. The problem is addressed by numerical means, carrying out three-dimensional relativistic magnetohydrodynamic simulations of bidirectional jets propagating in a triaxial density distribution. We show that only jets with power $\lesssim 10^{44}$ erg s$^{-1}$ can give origin to an X-shaped morphology and that a misalignment of $30^o$ between the jet axis and the major axis of the density distribution is still favourable to the formation of this kind of morphology. In addition we compute synthetic radio emission maps and polarization maps. In our scenario for the formation of X-shaped radio sources only low power FRII can give origin to such kind of morphology. Our synthetic emission maps show that the different observed morphologies of X-shaped sources can be the result of similar structures viewed under different perspectives.

Extragalactic radio sources have been classified into two classes, Fanaroff-Riley I and II, which differ in morphology and radio power. Strongly emitting sources belong to the edge-brightened FR II class, and weakly emitting sources to the edge-darkened FR I class. The origin of this dichotomy is not yet fully understood. Numerical simulations are successful in generating FR II morphologies, but they fail to reproduce the diffuse structure of FR Is. By means of hydro-dynamical 3D simulations of supersonic jets, we investigate how the displayed morphologies depend on the jet parameters. Bow shocks and Mach disks at the jet head, which are probably responsible for the hot spots in the FR II sources, disappear for a jet kinetic power L_kin < 10^43 erg/s. This threshold compares favorably with the luminosity at which the FR I/FR II transition is observed. The problem is addressed by numerical means carrying out 3D HD simulations of supersonic jets that propagate in a non-homogeneous medium with the ambient temperature that increases with distance from the jet origin, which maintains constant pressure. The jet energy in the lower power sources, instead of being deposited at the terminal shock, is gradually dissipated by the turbulence. The jets spread out while propagating, and they smoothly decelerate while mixing with the ambient medium and produce the plumes characteristic of FR I objects. Three-dimensionality is an essential ingredient to explore the FR I evolution because the properties of turbulence in two and three dimensions are very different, since there is no energy cascade to small scales in two dimensions, and two-dimensional simulations with the same parameters lead to FRII-like behavior.

An Equation of State (\textit{EoS}) closes the set of fluid equations. Although an ideal EoS with a constant \textit{adiabatic index} $\Gamma$ is the preferred choice due to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Here, we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation as well as temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation- equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes.We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice-versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS.The first one employs root-finder methods and it is best suited for EoS in analytical form. The second one leans on lookup table and interpolation and results in a more computationally efficient approach although care must be taken to ensure thermodynamic consistency. A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is $500-10^4$ K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although using lookup table methods can significantly reduce the overhead by a factor $3\sim 4$.

Disc-winds originating from the inner parts of accretion discs are considered as the basic component of magnetically collimated outflows. The only available analytical MHD solutions to describe disc-driven jets are those characterized by the symmetry of radial self-similarity. However, radially self-similar MHD jet models, in general, have three geometrical shortcomings, (i) a singularity at the jet axis, (ii) the necessary assumption of axisymmetry, and (iii) the non-existence of an intrinsic radial scale, i.e. the jets formally extend to radial infinity. Hence, numerical simulations are necessary to extend the analytical solutions towards the axis, by solving the full three-dimensional equations of MHD and impose a termination radius at finite radial distance. We focus here on studying the effects of relaxing the (ii) assumption of axisymmetry, i.e. of performing full 3D numerical simulations of a disc-wind crossing all magnetohydrodynamic critical surfaces. We compare the results of these runs with previous axisymmetric 2.5D simulations. The structure of the flow in all simulations shows strong similarities. The 3D runs reach a steady state and stay close to axisymmetry for most of the physical quantities, except for the poloidal magnetic field and the toroidal velocity which slightly deviate from axisymmetry. The latter quantities show signs of instabilities, which, however, are confined to the region inside the fast magnetosonic separatrix surface. The forces present in the flow, both of collimating and accelerating nature, are in good agreement in both the 2.5D and the 3D runs. We conclude that the analytical solution behaves well also after relaxing the basic assumption of axisymmetry.

Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.

It is a well established fact that some YSO jets (e.g. RW Aur) display different propagation speeds between their blue and red shifted parts, a feature possibly associated with the central engine or the environment in which the jet propagates. In order to understand the origin of asymmetric YSO jet velocities, we investigate the efficiency of two candidate mechanisms, one based on the intrinsic properties of the system and one based on the role of the external medium. In particular, a parallel or anti-parallel configuration between the protostellar magnetosphere and the disk magnetic field is considered and the resulting dynamics are examined both in an ideal and a resistive magneto-hydrodynamical (MHD) regime. Moreover, we explore the effects of a potential difference in the pressure of the environment, as a consequence of the non-uniform density distribution of molecular clouds. Ideal and resistive axisymmetric numerical simulations are carried out for a variety of models, all of which are based on a combination of two analytical solutions, a disk wind and a stellar outflow. We find that jet velocity asymmetries can indeed occur both when multipolar magnetic moments are present in the star-disk system as well as when non-uniform environments are considered. The latter case is an external mechanism that can easily explain the large time scale of the phenomenon, whereas the former one naturally relates it to the YSO intrinsic properties. [abridged]

We investigate supersonic, axisymmetric magnetohydrodynamic (MHD) jets with a time-dependent injection velocity by numerical simulations with the PLUTO code. Using a comprehensive set of parameters, we explore different jet configurations in the attempt to construct models that can be directly compared to observational data of microjets. In particular, we focus our attention on the emitting properties of traveling knots and construct, at the same time, accurate line intensity ratios and surface brightness maps. Direct comparison of the resulting brightness and line intensity ratios distributions with observational data of microjets shows that a closer match can be obtained only when the jet material is pre-ionized to some degree. A very likely source for a pre-ionized medium is photoionization by X-ray flux coming from the central object.

Stars form in regions of the galaxy that are denser and cooler than the mean interstellar medium. These regions are called Giant Molecular Clouds. At the beginning of their life, up to $10^5-10^6$ years, stars accrete matter from their rich surrounding environment and are origin of a peculiar phenomenon that is the jet emission. Jets from Young Stellar Objects (YSOs) are intensively studied by the astrophysical community by observations at different wavelengths, analytical and numerical modeling and laboratory experiments. Indications about the jet propagation and its resulting morphologies are here obtained by means of a combined study of hypersonic jets carried out both in the laboratory and by numerical simulations.

The coupling between time-dependent, multidimensional MHD numerical codes and radiative line emission is of utmost importance in the studies of the interplay between dynamical and radiative processes in many astrophysical environments, with particular interest for problems involving radiative shocks. There is a widespread consensus that line emitting knots observed in Herbig-Haro jets can be interpreted as radiative shocks. In this paper we address two different aspects relevant to the time-dependent calculations of the line intensity ratios of forbidden transitions, resulting from the excitation by planar, time-dependent radiative shocks traveling in a stratified medium. The first one concerns the impact of the radiation and ionization processes included in the cooling model, and the second one the effects of the numerical grid resolution. In this paper we apply the AMR methodology to the treatment of radiating shocks and show how this method is able to vastly reduce the integration time. The technique is applied to the knots of the HH 30 jet to obtain the observed line intensity ratios and derive the physical parameters, such as density, temperature and ionization fraction. We consider the impact of two different cooling functions and different grid resolutions on the results. We conclude that the use of different cooling routines has effects on results whose weight depends upon the line ratio considered. Moreover, we find the minimum numerical resolution of the simulation grid behind the shock to achieve convergence in the results. This is crucial for the forthcoming 2D calculations of radiative shocks.

Relativistic magnetized jets are key elements in Active Galactic Nuclei and in other astrophysical environments. Their structure and evolution involves a complex nonlinear 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 3D 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 an 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.

Recent observations as well as theoretical studies 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. In this framework, we construct numerical two-component jet models by properly mixing an analytical disk wind solution with a complementary analytically derived stellar outflow. Their combination is controlled by both spatial and temporal parameters, in order to address different physical conditions and time variable features. We study the temporal evolution and the interaction of the two jet components on both small and large scales. The simulations reach steady state configurations close to the initial solutions. Although time variability is not found to considerably affect the dynamics, flow fluctuations generate condensations, whose large scale structures have a strong resemblance to observed YSO jet knots.

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

We present the general properties of the Active Galactic Nuclei (AGNs) and discuss the origin and structure of jets that are associated to a fraction of these objects. We then we address the problems of particle acceleration at highly relativistic energies and set limits on the luminosity of AGN jets for being origin of UHECRs.

We investigate the linear stability properties of the plane interface separating two relativistic magnetized flows in relative motion. The two flows are governed by the (special) relativistic equations for a magnetized perfect gas in the infinite conductivity approximation. By adopting the vortex-sheet approximation, the relativistic magnetohydrodynamics equations are linearized around the equilibrium state and the corresponding dispersion relation is derived and discussed. The behavior of the configuration and the regimes of instability are investigated following the effects of four physical parameters, namely: the flow velocity, the relativistic and Alfv\'enic Mach numbers and the inclination of the wave vector on the plane of the interface. From the numerical solution of the dispersion relation, we find in general two separate regions of instability, associated respectively with the slow and fast magnetosonic modes. Modes parallel to the flow velocity are destabilized only for sufficiently low magnetization. For the latter case, stabilization is attained, additionally, at sufficiently large relativistic velocities between the two flows in relative motion. The relevance of these results to the study of the stability of astrophysical jets is briefly commented.

Time-dependent cooling processes are of paramount importance in the evolution of astrophysical gaseous nebulae and, in particular, when radiative shocks are present. The present work introduces a necessary set of tools that can be used to model radiative astrophysical flows in the optically-thin plasma limit. We aim to provide reliable and accurate predictions of emission line ratios and radiative cooling losses in astrophysical simulations of shocked flows. Moreover, we discuss numerical implementation aspects to ease future improvements and implementation in other MHD numerical codes. The most important source of radiative cooling for our plasma conditions comes from the collisionally-excited line radiation. We evolve a chemical network, including 29 ion species, to compute the ionization balance in non-equilibrium conditions. After a series of validations and tests, typical astrophysical setups are simulated in 1D and 2D, employing both the present cooling model and a simplified one. The influence of the cooling model on structure morphologies can become important, especially for emission line diagnostic purposes.

Strong observational evidence indicates that all extragalactic jets associated with AGNs move at relativistic speed up to 100 pc - 1 kpc scales from the nucleus. At larger distances, reflecting the Fanaroff-Riley radio source classification, we observe an abrupt deceleration 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 deboosted central spine at high Lorentz factor and of a weakly relativistic external layer. In this paper we investigate whether these effects can be interpreted by a breaking of the collimated flow by external medium entrainment favored by shear instabilities, namely Kelvin-Helmholtz instabilities. We examine in details the physical conditions under which significant deceleration of a relativistic flow is produced. We investigate the phenomenon by means of high-resolution three-dimensional relativistic hydrodynamic simulations using the PLUTO code for computational astrophysics. We find that the parameter of utmost importance in determining 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. Our model is verified by constructing synthetic emission maps from the numerical simulations that compare reasonably well with VLBI observations of the inner part of FR-I sources.

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. In this context, the present article starts a systematic investigation of two-component jet models 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. By 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. It is found that radially self-similar solutions (disk-winds) always reach a final steady-state while maintaining all their well-defined 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. 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.

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 multi-physics, multi-algorithm modular environment particularly oriented towards 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.

Accretion disks and astrophysical jets are used to model many active astrophysical objects, viz., young stars, relativistic stars, and active galactic nuclei. In this paper 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 disk magnetic resistivity, parametrized through an alpha-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 and Payne can be reached and maintained in a self consistent nonlinear stage. We use the resistive MHD FLASH code with adaptive mesh refinement, allowing us to follow the evolution of the structure for a time scale long enough to reach steady state. A detailed analysis of the initial configuration state is given. We obtain the expected solutions in the axisymmetric (2.5D) limit. Assuming a magnetic field around equipartition with the thermal pressure of the disk, we show how the characteristics of the disk jet system, as the ejection efficiency and the energetics, are affected by the anomalous resistivity acting inside the disk.

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.

The line emission from a growing number of Herbig-Haro jets can be observed and resolved at angular distances smaller than a few arcseconds from the central source. The interpretation of this emission is problematic, since the simplest model of a cooling jet cannot sustain it. It has been suggested that what one actually observes are shocked regions with a filling factor of $\sim 1%$. In this framework, up to now, comparisons with observations have been based on stationary shock models. Here we introduce for the first time the self-consistent dynamics of such shocks and we show that considering their properties at different times, i.e. locations, we can reproduce observational data of the DG Tau microjet. In particular, we can interpret the spatial behavior of the [SII]6716/6731 and [NII]/[OI]6583/6300 line intensity ratios adopting a set of physical parameters that yield values of mass loss rates and magnetic fields consistent with previous estimates. We also obtain the values of the mean ionization fraction and electron density along the jet, compare these values with the ones derived from observations using the sulfur doublet to constrain the electron density (e.g. Bacciotti et al. 1995).

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.

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

Comparing the properties AGN and Herbig-Haro jets can be a useful exercise for understanding the physical mechanisms at work in collimated outflows that propagate in such different environments. In the case of Herbig-Haro jets, the presence of emission lines in the spectra and the continuous evolution of the observation techniques greatly favor our knowledge of the physical parameters of the jets instead, for AGN jets, the process of constraining the jet parameters is hampered by the nature of the emission from these objects that is non-thermal. I will discuss how one cannot directly constrain the basic parameters of extragalactic jets by observations but must treat and interpret the data either by statistical means or by comparing observed and simulated morphologies in order to gain some indications on the values of these parameters.

We show that a helical shear flow of a magnetized plasma may serve as an efficient amplifier of Alfven waves. We find that even when the flow is purely ejectional (i.e., when no rotation is present) Alfven waves are amplified through the transient, shear-induced, algebraic amplification process. Series of transient amplifications, taking place sequentially along the flow, may result in a cascade amplification of these waves. However, when a flow is swirling or helical (i.e., some rotation is imposed on the plasma motion), Alfven waves become subject to new, much more powerful shear instabilities. In this case, depending on the type of differential rotation, both usual and parametric instabilities may appear. We claim that these phenomena may lead to the generation of large amplitude Alfven waves and the mechanism may account for the appearance of such waves in the solar atmosphere, in accretion-ejecion flows and in accretion columns. These processes may also serve as an important initial (linear and nonmodal) phase in the ultimate subcritical transition to MHD Alfvenic turbulence in various kinds of astrophysical shear flows.

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.

We present the results of a BeppoSAX observation of the Supernova Remnant MSH 15-52, associated with the pulsar PSR B1509-58, and discuss its main morphological and spectroscopic properties in the 1.6--200 keV energy range (MECS and PDS instruments). The two main structures of the remnant, the Southern Nebula, the plerion centered on the pulsar, and the Northern Nebula, are clearly visible in the MECS, with the former showing a much a harder spectrum. Furthermore, a diffuse extended emission surrounds the whole remnant up to ~ 17' from the center. Non-thermal flux is detected in the PDS up to 200 keV as well, and it appears that also in this energy range the emission is not concentrated in the central region around the pulsar. These data imply that the plerion extends up to a few tens of parsecs from the pulsar.

We present BeppoSAX observations of the two FR I type radio galaxies Hercules A (3C 348) and MRC 0625-536 in the energy range $0.2 - 200$ keV. Data analysis shows that the X-ray flux from Hercules A is consistent with a diffuse thermal plasma emitting at $T \approx 4-5$ keV with a possible, but somewhat uncertain, contribution of a softer component at $T \approx 3 $ keV. The non thermal emission from the active nucleus must be significantly smaller than the thermal one, and no indication of relevant core obscuration by a surrounding torus was detected. The flux from MRC 0625-536 originates form an extended region and has been fitted to a thermal law with $T \approx 5.7$ keV and with a column density consistent with the galactic absorption. A spatially resolved spectral analysis does not show a relevant variation of the temperature and the metallicity across the diffuse emission zone. A non thermal spectral component, related to the nuclear activity, may be present in the innermost region with some possible amount of local obscuration, contributing $\lapp 10%$ to the total luminosity. Hard X-ray emission from MRC 0625-536 has been detected in the PDS (15 - 200 keV) that may be related either to its galactic core or to the intracluster region.

We obtained EVN observations of the X-ray transient source GRS1915+105 in its radio-loud state. In 2 images obtained before a major flare, the source is resolved, showing two opposite asymmetric jets. This indicates that steady jets exist also before a major flare, with an estimated velocity of 0.2-0.6c.

We study, by numerical simulations, the entrainment process in a supersonic, radiative jet flow, during the evolution of Kelvin-Helmholtz instabilities, in the context of the the acceleration problem of molecular bipolar outflows, observed in Giant Molecular Clouds. Our results show that a large fraction of the initial jet momentum can be transferred to the ambient medium by this process. We therefore analyze in detail the instability evolution and compare some of the main observational properties of molecular outflows with those of the entrained material that we observe in our simulations. In particular, we find a good agreement for the mass vs. velocity distribution and for the outflow collimation structure, especially when a light jet is moving into a denser ambient medium. This is probably the case for (obscured) optical jets driving powerful molecular outflows in the denser environment of the inner regions of molecular clouds.

The analysis of the stability properties of astrophysical jets against Kelvin-Helmholtz (or shear-layer) instabilities plays a basic role in the understanding the origin and physical characteristics of these objects. Numerical simulations by Bodo et al. (1998) have shown that the three-dimensional non-linear evolution of KH instabilities in supersonic jets is substantially faster than in the two-dimensional case, leading to a cascade of modes towards smaller scales and a very effective mixing and momentum transfer to the ambient medium. On the other hand, Rossi et al. (1997) and Micono et al. (1998) found, in two dimensions, that radiative losses tend to reduce and delay mixing effects and momentum transfer to the ambient medium. In this paper, as a logical next step, we investigate the effects of radiative losses on the stability of 3D supersonic jets, assuming that the internal jet density is initially lower, equal and higher than the ambient medium, respectively. We find that light and equal-density radiative jets evolve in a qualitatively similar fashion with respect to the corresponding adiabatic ones. Conversely, we note substantial differences in the evolution of heavy jets: they remain more collimated and do not spread out, while the momentum gained by the ambient medium stays within ~ 5 jet radii.

Recent HST observations suggest that the NLR in Seyfert Galaxies can be the result of interaction between jet and external inhomogeneous medium; following this suggestion we perform numerical simulations considering the impact of a radiative jet on a dense cloud. We approach the problem adopting a hydrodynamical code, that consents us to study in detail the jet hydrodynamics, while we choose a more simplified treatment of radiative processes, in order to give a qualitatively good interpretation of the emission processes. Our three main purposes are: i) to reproduce in our simulations the physical conditions observed in the NLR of Seyfert Galaxies, ii) to obtain physical constraints of the jet parameters and iii) to study the jet capability to photoionize the surrounding medium.

We investigate analytically the propagation of linear waves in a three-dimensional, nonmagnetic, isothermal atmosphere stratified in plane-parallel layers. The motivation is to study oscillations in the non-magnetic chromosphere and to assess the limitations of one-dimensional simulations of the Kv bright point phenomenon. We consider an impulsively excited acoustic disturbance, emanating from a point source, and propagating outward as a spherical acoustic wave accompanied by an internal gravity wave. The waves amplify exponentially in the upward direction. A significant wave amplitude is therefore found only in a relatively narrow cone about the vertical. The amplitude of the wave decreases with time. Because of the lateral spread, the wave amplitude decays faster in 2D and 3D simulations than in 1D. The initial pulse, which travels at the sound speed, carries most of the energy injected into the medium. Subsequent wave crests leave the source region at ever-increasing phase speed, but slow to the sound speed as they approach the head of the wave. Important conclusions from the 3D solution that were not anticipated from the plane-wave solution are: 1. The bulk of the energy is emitted in the upward (and downward) direction; much less goes into the horizontal direction. 2. The wave profile narrows from the initial pulse through the amplitude maxima in the wake of the pulse. As a consequence of both points, the shock-heated regions in the wake of the initial pulse would weaken in strength and shrink in size.

We calculate the temporal evolution of distributions of relativistic electrons subject to synchrotron and adiabatic processes and Fermi-like acceleration in shocks. The shocks result from Kelvin-Helmholtz instabilities in the jet. Shock formation and particle acceleration are treated in a self-consistent way by means of a numerical hydrocode. We show that in our model the number of relativistic particles is conserved during the evolution, with no need of further injections of supra-thermal particles after the initial one. From our calculations, we derive predictions for values and trends of quantities like the spectral index and the cutoff frequency that can be compared with observations.

We present the first results of BeppoSAX observations of a small sample of low brightness FRI radio galaxies. The flux of all the targets is consistent with a thermal spectrum, as due to the presence of hot intracluster gas or galactic corona. Moreover in three sources a non thermal absorbed spectrum can be present in the MECS spectrum at energies larger than 7 keV, while for a fourth object a high energy flux has been detected in the PDS instrument at energies larger than 15 keV. This component could be related to the inner AGN surrounded by an obscuring torus.

We present the results of a numerical analysis of the propagation and interaction of a supersonic jet with the external medium. We discuss the motion of the head of the jet into the ambient in different physical conditions, carrying out calculations with different Mach numbers and density ratios of the jet to the exteriors. Performing the calculation in a reference frame in motion with the jet head, we can follow in detail its long term dynamics. This numerical scheme allows us also to study the morphology of the cocoon for different physical parameters. We find that the propagation velocity of the jet head into the ambient medium strongly influences the morphology of the cocoon, and this result can be relevant in connection to the origin and structure of lobes in extragalactic radiosources.

We report X-ray observations of the radio galaxies B2 1615+35 (NGC 6109) and B2 1621+38 (NGC 6137), which belong to poor groups of galaxies and are characterized by a "head-tail" structure. We also analyze the X-ray behaviour of the radio galaxy NGC 6107, present within the field of view and showing a symmetric extended radio structure. The X-ray emission from the group of 1615+35 is complex, including both emission from a hot intergroup gas, and emission from the individual radio galaxies (1615+35 itself, which is pointlike, and NGC 6107, which is extended). In the case of 1621+38, only extended X-ray emission from the galaxy itself is detected. The extended X-ray emission of the galaxies 1621+38 and NGC 6107 is probably related to a gaseous galactic atmosphere. The pointlike emission of the radio galaxy 1615+35, instead, is likely to originate in the nucleus. We discuss the confinement of radio components by the ambient gas, and find that the external thermal pressure exceeds the internal equipartition pressure. The X-ray emitting plasma in the 1615+35 group is characterized by luminosity, temperature and metallitcity typical for groups, i.e. lower than for rich clusters.