We investigate with relativistic MHD simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically-dominated plasma column. The nonlinear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = J d/B is above some predefined threshold (where B is the field strength, J the current density and d the grid scale), and assume that the particle injection efficiency scales as proportional to the square of the current. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical-emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multi-wavelength emission from BL Lacs.

We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces and are evolved using the constrained transport formalism. Direct application of Stokes' theorem to Faraday's and Ampere's laws ensures that the resulting discretization is divergence-free for the magnetic field and charge-conserving for the electric field. Hydrodynamic variables retain, instead, the usual zone-centred representation commonly adopted in finite-volume schemes. Temporal discretization is based on Runge-Kutta implicit-explicit (IMEX) schemes in order to resolve the temporal scale disparity introduced by the stiff source term in Ampere's law. The implicit step is accomplished by means of an improved and more efficient Newton-Broyden multidimensional root-finding algorithm. The explicit step relies on a multidimensional Riemann solver to compute the line-averaged electric and magnetic fields at zone edges and it employs a one-dimensional Riemann solver at zone interfaces to update zone-centred hydrodynamic quantities. For the latter, we introduce a five-wave solver based on the frozen limit of the relaxation system whereby the solution to the Riemann problem can be decomposed into an outer Maxwell solver and an inner hydrodynamic solver. A number of numerical benchmarks demonstrate that our method is superior in stability and robustness to the more popular charge-conserving divergence cleaning approach where both primary electric and magnetic fields are zone-centered. In addition, the employment of a less diffusive Riemann solver noticeably improves the accuracy of the computations.

We carry out a linear stability analysis of a magnetized relativistic rotating cylindrical jet flow using the approximation of zero thermal pressure. We identify several modes of instability in the jet: Kelvin-Helmholtz, current driven and two kinds of centrifugal-buoyancy modes -- toroidal and poloidal. The Kelvin-Helmholtz mode is found at low magnetization and its growth rate depends very weakly on the pitch parameter of the background magnetic field and on rotation. The current driven mode is found at high magnetization, the values of its growth rate and the wavenumber, corresponding to the maximum growth, increase as we decrease the pitch parameter of the background magnetic field. This mode is stabilized by rotation, especially, at high magnetization. The centrifugal-buoyancy modes, arising due to rotation, tend also to be more stable when magnetization is increased. Overall, relativistic jet flows appear to be more stable with respect to their non-relativistic counterpart.

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 present a linear mode analysis of the relativistic MHD equations in the presence of finite electrical conductivity. Starting from the fully relativistic covariant formulation, we derive the dispersion relation in the limit of small linear perturbations. It is found that the system supports ten wave modes which can be easily identified in the limits of small or large conductivities. In the resistive limit, matter and electromagnetic fields decouple and solution modes approach pairs of light and acoustic waves as well as a number of purely damped (non-propagating) modes. In the opposite (ideal) limit, the frozen-in condition applies and the modes of propagation coincide with a pair of fast magnetosonic, a pair of slow and Alfv\'en modes, as expected. In addition, the contact mode is always present and it is unaffected by the conductivity. For finite values of the conductivity, the dispersion relation gives rise to either pairs of opposite complex conjugate roots or purely imaginary (damped) modes. In all cases, the system is dissipative and also dispersive as the phase velocity depends nonlineary on the wavenumber. Occasionally, the group velocity may exceed the speed of light although this does not lead to superluminal signal propagation.

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.

We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged particles, or cosmic rays (CR). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques. The module can be used to describe either test particles motion in the fluid electromagnetic field or to solve the fully coupled MHD-PIC system of equations with particle backreaction on the fluid as originally introduced by \cite{Bai_etal.2015}. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm's law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided the Larmor gyration scale is properly resolved. When applicable, this formulation avoids to resolve microscopic scales offering a substantial computational saving with respect to PIC simulations. We present a fully-conservative formulation which is second-order accurate in time and space and extends to either Runge-Kutta (RK) or corner-transport-upwind (CTU) time-stepping schemes (for the fluid) while a standard Boris integrator is employed for the particles. For highly-energetic relativistic CRs and in order to overcome the time step restriction a novel sub-cycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration and test particle acceleration in reconnecting layers are discussed.

We present the results of state-of-the-art simulations of recollimation shocks induced by the interaction of a relativistic jet with an external medium, including the effect of radiative losses of the shocked gas. Our simulations confirm that -- as suggested by earlier semi-analytical models -- the post-shock pressure loss induced by radiative losses may lead to a stationary equilibrium state characterized by a very strong focusing of the flow, with the formation of quite narrow nozzles, with cross-sectional radii as small as $10^{-3}$ times the length scale of the jet. We also study the time-dependent evolution of the jet structure induced of a density perturbation injected at the flow base. The set-up and the results of the simulations are particularly relevant for the interpretation of the observed rapid variability of the $\gamma$-ray emission associated to flat spectrum radio quasars. In particular, the combined effects of jet focusing and Doppler beaming of the observed radiation make it possible to explain the sub-hour flaring events such as that observed in the FSRQ PKS 1222+216 by MAGIC.

We investigate magnetohydrodynamic turbulence driven by the magnetorotational instability (MRI) in Keplerian disks with a nonzero net azimuthal magnetic field using shearing box simulations. As distinct from most previous studies, we analyze turbulence dynamics in Fourier (${\bf k}$-) space to understand its sustenance. The linear growth of MRI with azimuthal field has a transient character and is anisotropic in Fourier space, leading to anisotropy of nonlinear processes in Fourier space. As a result, the main nonlinear process appears to be a new type of angular redistribution of modes in Fourier space -- the \emph{nonlinear transverse cascade} -- rather than usual direct/inverse cascade. We demonstrate that the turbulence is sustained by interplay of the linear transient growth of MRI (which is the only energy supply for the turbulence) and the transverse cascade. These two processes operate at large length scales, comparable to box size and the corresponding small wavenumber area, called \emph{vital area} in Fourier space is crucial for the sustenance, while outside the vital area direct cascade dominates. The interplay of the linear and nonlinear processes in Fourier space is generally too intertwined for a vivid schematization. Nevertheless, we reveal the \emph{basic subcycle} of the sustenance that clearly shows synergy of these processes in the self-organization of the magnetized flow system. This synergy is quite robust and persists for the considered different aspect ratios of the simulation boxes. The spectral characteristics of the dynamical processes in these boxes are qualitatively similar, indicating the universality of the sustenance mechanism of the MRI-turbulence.

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.

We examine the relationship between magnetic flux generation, taken as an indicator of large-scale dynamo action, and magnetic helicity, computed as an integral over the dynamo volume, in a simple dynamo. We consider dynamo action driven by Magneto-Rotational Turbulence (MRT) within the shearing-box approximation. We consider magnetically open boundary conditions that allow a flux of helicity in or out of the computational domain. We circumvent the problem of the lack of gauge invariance in open domains by choosing a particular gauge -- the winding gauge -- that provides a natural interpretation in terms of average winding number of pairwise field lines. We use this gauge precisely to define and measure the helicity and helicity flux for several realizations of dynamo action. We find in these cases, that the system as a whole does not break reflectional symmetry and the total helicity remains small even in cases when substantial magnetic flux is generated. We find no particular connection between the generation of magnetic flux and the helicity or the helicity flux through the boundaries. We suggest that this result may be due to the essentially nonlinear nature of the dynamo processes in MRT.

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.

Magnetic reconnection is a plasma phenomenon where a topological rearrangement of magnetic field lines with opposite polarity results in dissipation of magnetic energy into heat, kinetic energy and particle acceleration. Such a phenomenon is considered as an efficient mechanism for energy release in laboratory and astrophysical plasmas. An important question is how to make the process fast enough to account for observed explosive energy releases. The classical model for steady state magnetic reconnection predicts reconnection times scaling as $S^{1/2}$ (where $S$ is the Lundquist number) and yields times scales several order of magnitude larger than the observed ones. Earlier two-dimensional MHD simulations showed that for large Lundquist number the reconnection time becomes independent of $S$ ("fast reconnection" regime) due to the presence of the secondary tearing instability that takes place for $S \gtrsim 1 \times 10^4$. We report on our 3D MHD simulations of magnetic reconnection in a magnetically confined cylindrical plasma column under either a pressure balanced or a force-free equilibrium and compare the results with 2D simulations of a circular current sheet. We find that the 3D instabilities acting on these configurations result in a fragmentation of the initial current sheet in small filaments, leading to enhanced dissipation rate that becomes independent of the Lundquist number already at $S \simeq 1\times 10^3$.

We perform a linear stability analysis of magnetized rotating cylindrical jet flows in the approximation of zero thermal pressure. We focus our analysis on the effect of rotation on the current driven mode and on the unstable modes introduced by rotation. We find that rotation has a stabilizing effect on the current driven mode only for rotation velocities of the order of the Alfv\'en velocity. Rotation introduces also a new unstable centrifugal buoyancy mode and the "cold" magnetorotational instability. The first mode is analogous to the Parker instability with the centrifugal force playing the role of effective gravity. The magnetorotational instability can be present, but only in a very limited region of the parameter space and is never dominant. The current driven mode is characterized by large wavelenghts and is dominant at small values of the rotational velocity, while the buoyancy mode becomes dominant as rotation is increased and is characterized by small wavelenghts.

We present an approach to deriving global properties of accretion disks from the knowledge of local solutions derived from numerical simulations based on the shearing box approximation. The approach consists of a two-step procedure. First a local solution valid for all values of the disk height is constructed by piecing together an interior solution obtained numerically with an analytical exterior radiative solution. The matching is obtained by assuming hydrostatic balance and radiative equilibrium. Although in principle the procedure can be carried out in general, it simplifies considerably when the interior solution is fully convective. In these cases, the construction is analogous to the derivation of the Hayashi tracks for protostars. The second step consists of piecing together the local solutions at different radii to obtain a global solution. Here we use the symmetry of the solutions with respect to the defining dimensionless numbers--in a way similar to the use of homology relations in stellar structure theory--to obtain the scaling properties of the various disk quantities with radius.

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

We investigate the linear and nonlinear evolution of current-carrying jets in a periodic configuration by means of high resolution three-dimensional numerical simulations. The jets under consideration are strongly magnetized with a variable pitch profile and initially in equilibrium under the action of a force-free magnetic field. The growth of current-driven (CDI) and Kelvin-Helmholtz (KHI) instabilities is quantified using three selected cases corresponding to static, Alfvenic and super-Alfvenic jets. During the early stages, we observe large-scale helical deformations of the jet corresponding to the growth of the initially excited CDI mode. A direct comparison between our simulation results and the analytical growth rates obtained from linear theory reveals good agreement on condition that high-resolution and accurate discretization algorithms are employed. After the initial linear phase, the jet structure is significantly altered and, while slowly-moving jets show increasing helical deformations, larger velocity shear are violently disrupted on a few Alfven crossing time leaving a turbulent flow structure. Overall, kinetic and magnetic energies are quickly dissipated into heat and during the saturated regime the jet momentum is redistributed on a larger surface area with most of the jet mass travelling at smaller velocities. The effectiveness of this process is regulated by the onset of KHI instabilities taking place at the jet/ambient interface and can be held responsible for vigorous jet braking and entrainment.

We consider the problem of convergence in stratified isothermal shearing boxes with zero net magnetic flux. We present results with the highest resolution to-date--up to 200 grid-point per pressure scale height--that show no clear evidence of convergence. Rather, the Maxwell stresses continue to decrease with increasing resolution. We propose some possible scenarios to explain the lack of convergence based on multi-layer dynamo systems.

We investigate linear dynamics of non-axisymmetric perturbations in incompressible, vertically stratified Keplerian discs with a weak vertical magnetic field in the shearing box approximation. Perturbations are decomposed into shearing waves whose evolution is followed via numerical integration of the linearized ideal MHD equations. There are two basic modes in the system -- inertia-gravity waves and magnetic mode, which displays the magnetorotational instability (MRI). As distinct from previous studies, we introduce `eigenvariables' characterizing each (counter-propagating) component of the inertia-gravity and magnetic modes, which are governed by a set of four first order coupled ordinary differential equations. This allowed us to identify a new process of linear coupling of the two above non-axisymmetric modes due to the disc's differential rotation. We did a comparative analysis of the dynamics of non-axisymmetric and axisymmetric magnetic mode perturbations. It is shown that the growth of optimal and close-to-optimal non-axisymmetric harmonics of this mode, having transient nature, can prevail over the exponential growth of axisymmetric ones (i.e., over the axisymmetric MRI) during dynamical time. A possible implication of this result for axisymmetric channel solutions is discussed. Specifically, the formation of the channel may be affected/impeded by non-axisymmetric modes already at the early linear stage leading to its untimely disruption -- the outcome strongly depends on the amplitude and spectrum of initial perturbation. So, this competition may result in an uncertainty in the magnetic mode's non-linear dynamics. It is also shown that a maximum non-axisymmetric growth is at vertical wavelengths close to the scale-height for which compressibility effects are important. This indirectly suggests that compressibility plays a role in the dynamics of the non-axisymmetric MRI.

We perform a linear analysis of the stability of a magnetized relativistic non-rotating cylindrical flow in the aproximation of zero thermal pressure, considering only the m = 1 mode. We find that there are two modes of instability: Kelvin-Helmholtz and current driven. The Kelvin-Helmholtz mode is found at low magnetizations and its growth rate depends very weakly on the pitch parameter. The current driven modes are found at high magnetizations and the value of the growth rate and the wavenumber of the maximum increase as we decrease the pitch parameter. In the relativistic regime the current driven mode is splitted in two branches, the branch at high wavenumbers is characterized by the eigenfunction concentrated in the jet core, the branch at low wavenumbers is instead characterized by the eigenfunction that extends outside the jet velocity shear region.

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. We choose radiative boundary conditions at the vertical boundaries in which the heat flux is propor- tional to the fourth power of the temperature. We compare the results with the corresponding cases in which fixed temperature boundary conditions are applied. The most notable result is that the formation of a fully convective state in which the density is nearly constant as a function of height and the heat is transported to the upper and lower boundaries by overturning motions is robust and persists even in cases with radiative boundary conditions. Interestingly, in the convective regime, although the diffusive transport is negligible the mean stratification does not relax to an adiabatic state.

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero mean magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. The calculations begin from an isothermal state spanning three scale heights above and below the mid-plane. After a long transient the layers settle to a stationary state in which thermal losses out of the boundaries are balanced by dissipative heating. We identify two regimes. A conductive regime in which the heat is transported mostly by conduction and the density decreases with height. In the limit of large thermal diffusivity this regime resembles the more familiar isothermal case. Another, the convective regime, observed at smaller values of the thermal diffusivity, in which the layer becomes unstable to overturning motions, the heat is carried mostly by advection and the density becomes nearly constant throughout the layer. In this latter constant-density regime we observe evidence for large-scale dynamo action leading to a substantial increase in transport efficiency relative to the conductive cases.

We present a description of the adaptive mesh refinement (AMR) implementation of the PLUTO code for solving the equations of classical and special relativistic magnetohydrodynamics (MHD and RMHD). The current release exploits, in addition to the static grid version of the code, the distributed infrastructure of the CHOMBO library for multidimensional parallel computations over block-structured, adaptively refined grids. We employ a conservative finite-volume approach where primary flow quantities are discretized at the cell-center in a dimensionally unsplit fashion using the Corner Transport Upwind (CTU) method. Time stepping relies on a characteristic tracing step where piecewise parabolic method (PPM), weighted essentially non-oscillatory (WENO) or slope-limited linear interpolation schemes can be handily adopted. A characteristic decomposition-free version of the scheme is also illustrated. The solenoidal condition of the magnetic field is enforced by augmenting the equations with a generalized Lagrange multiplier (GLM) providing propagation and damping of divergence errors through a mixed hyperbolic/parabolic explicit cleaning step. Among the novel features, we describe an extension of the scheme to include non-ideal dissipative processes such as viscosity, resistivity and anisotropic thermal conduction without operator splitting. Finally, we illustrate an efficient treatment of point-local, potentially stiff source terms over hierarchical nested grids by taking advantage of the adaptivity in time. Several multidimensional benchmarks and applications to problems of astrophysical relevance assess the potentiality of the AMR version of PLUTO in resolving flow features separated by large spatial and temporal disparities.

We consider the problem of convergence in homogeneous shearing box simula- tions of magneto-rotationally driven turbulence. When there is no mean magnetic flux, if the equations are non dimensionalized with respect to the diffusive scale, the only free parameter in the problem is the size of the computational domain. The problem of convergence then relates to the asymptotic form of the solutions as the computational box size becomes large. By using a numerical code with a high order of accuracy we show that the solutions become asymptotically inde- pendent of domain size. We also show that cases with weak magnetic flux join smoothly to the zero flux cases as the flux vanishes. These results are consistent with the operation of a subcritical small-scale dynamo driving the turbulence. We conclude that for this type of turbulence the angular momentum transport is a proportional to the diffusive flux and therefore has limited relevance in as- trophysical situations.

We present and compare third- as well as fifth-order accurate finite difference schemes for the numerical solution of the compressible ideal MHD equations in multiple spatial dimensions. The selected methods lean on four different reconstruction techniques based on recently improved versions of the weighted essentially non-oscillatory (WENO) schemes, monotonicity preserving (MP) schemes as well as slope-limited polynomial reconstruction. The proposed numerical methods are highly accurate in smooth regions of the flow, avoid loss of accuracy in proximity of smooth extrema and provide sharp non-oscillatory transitions at discontinuities. We suggest a numerical formulation based on a cell-centered approach where all of the primary flow variables are discretized at the zone center. The divergence-free condition is enforced by augmenting the MHD equations with a generalized Lagrange multiplier yielding a mixed hyperbolic/parabolic correction, as in Dedner et al. (J. Comput. Phys. 175 (2002) 645-673). The resulting family of schemes is robust, cost-effective and straightforward to implement. Compared to previous existing approaches, it completely avoids the CPU intensive workload associated with an elliptic divergence cleaning step and the additional complexities required by staggered mesh algorithms. Extensive numerical testing demonstrate the robustness and reliability of the proposed framework for computations involving both smooth and discontinuous features.

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.

We investigate mode coupling in a two dimensional compressible disc with radial stratification and differential rotation. We employ the global radial scaling of linear perturbations and study the linear modes in the local shearing sheet approximation. We employ a three-mode formalism and study the vorticity (W), entropy (S) and compressional (P) modes and their coupling properties. The system exhibits asymmetric three-mode coupling: these include mutual coupling of S and P-modes, S and W-modes, and asymmetric coupling between the W and P-modes. P-mode perturbations are able to generate potential vorticity through indirect three-mode coupling. This process indicates that compressional perturbations can lead to the development of vortical structures and influence the dynamics of radially stratified hydrodynamic accretion and protoplanetary discs.

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.

We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for the equations of ideal MHD. The solution to the Riemann problem is approximated by a five wave pattern, comprised of two outermost fast shocks, two rotational discontinuities and a contact surface in the middle. The proposed scheme is considerably more elaborate than in the classical case since the normal velocity is no longer constant across the rotational modes. Still, proper closure to the Rankine-Hugoniot jump conditions can be attained by solving a nonlinear scalar equation in the total pressure variable which, for the chosen configuration, has to be constant over the whole Riemann fan. The accuracy of the new Riemann solver is validated against one dimensional tests and multidimensional applications. It is shown that our new solver considerably improves over the popular HLL solver or the recently proposed HLLC schemes.

Recent Cerenkov observations of the two BL Lac objects PKS 2155-304 and Mkn 501 revealed TeV flux variability by a factor ~2 in just 3-5 minutes. Even accounting for the effects of relativistic beaming, such short timescales are challenging simple and conventional emitting models, and call for alternative ideas. We explore the possibility that extremely fast variable emission might be produced by particles streaming at ultra-relativistic speeds along magnetic field lines and inverse Compton scattering any radiation field already present. This would produce extremely collimated beams of TeV photons. While the probability for the line of sight to be within such a narrow cone of emission would be negligibly small, one would expect that the process is not confined to a single site, but can take place in many very localised regions, along almost straight magnetic lines. A possible astrophysical setting realising these conditions is magneto-centrifugal acceleration of beams of particles. In this scenario, the variability timescale would not be related to the physical dimension of the emitting volume, but might be determined by either the typical duration of the process responsible for the production of these high energy particle beams or by the coherence length of the magnetic field. It is predicted that even faster TeV variability - with no X-ray counterpart - should be observed by the foreseen more sensitive Cerenkov telescopes.

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.

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.

Aims: We study the changes in the properties of turbulence driven by the magnetorotational instability in a shearing box, as the computational domain size in the radial direction is varied relative to the height Methods: We perform 3D simulations in the shearing box approximation, with a net magnetic flux, and we consider computational domains with different aspect ratios Results: We find that in boxes of aspect ratio unity the transport of angular momentum is strongly intermittent and dominated by channel solutions in agreement with previous work. In contrast, in boxes with larger aspect ratio, the channel solutions and the associated intermittent behavior disappear. Conclusions: There is strong evidence that, as the aspect ratio becomes larger, the characteristics of the solution become aspect ratio independent. We conclude that shearing box calculations with aspect ratio unity or near unity may introduce spurious effects.

We investigate the stability, nonlinear development and equilibrium structure of vortices in a background shearing Keplerian flow. We make use of high-resolution global two-dimensional compressible hydrodynamic simulations. We introduce the concept of nonlinear adjustment to describe the transition of unbalanced vortical fields to a long-lived configuration. We discuss the conditions under which vortical perturbations evolve into long-lived persistent structures and we describe the properties of these equilibrium vortices. The properties of equilibrium vortices appear to be independent from the initial conditions and depend only on the local disk parameters. In particular we find that the ratio of the vortex size to the local disk scale height increases with the decrease of the sound speed, reaching values well above the unity. The process of spiral density wave generation by the vortex, discussed in our previous work, appear to maintain its efficiency also at nonlinear amplitudes and we observe the formation of spiral shocks attached to the vortex. The shocks may have important consequences on the long term vortex evolution and possibly on the global disk dynamics. Our study strengthens the arguments in favor of anticyclonic vortices as the candidates for the promotion of planetary formation. Hydrodynamic shocks that are an intrinsic property of persistent vortices in compressible Keplerian flows are an important contributor to the overall balance. These shocks support vortices against viscous dissipation by generating local potential vorticity and should be responsible for the eventual fate of the persistent anticyclonic vortices. Numerical codes have be able to resolve shock waves to describe the vortex dynamics correctly.

To investigate the efficiency of centrifugal acceleration of particles as a possible mechanism for the generation of ultra-high $\gamma$-ray nonthermal emission from TeV blazars, we study the centrifugal acceleration of electrons by rotating magnetic field lines, for an extended range of inclination angles and determine the maximum Lorentz factors $\gamma_{max}$ attainable by the electrons via this process. {Two principal limiting mechanisms for the particle acceleration, inverse Compton scattering and breakdown of the bead-on-the-wire approximation, are examined. Particles may be centrifugally accelerated up to $\gamma_{max} \simeq 10^8$ and the main limiting mechanism for the $\gamma_{max}$ is the inverse Compton scattering. The energy of centrifugally accelerated particles can be amply sufficient for the generation (via inverse Compton scattering) of the ultra-high energy (up to $20TeV$) gamma emission in TeV blazars.

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.

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The HLLC solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper (Mignone & Bodo 2004) to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave. Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the contrained tranport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one and two dimensional numerical test problems.

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.

We present an extension of the HLLC approximate Riemann solver by Toro, Spruce and Speares to the relativistic equations of fluid dynamics. The solver retains the simplicity of the original two-wave formulation proposed by Harten, Lax and van Leer (HLL) but it restores the missing contact wave in the solution of the Riemann problem. The resulting numerical scheme is computationally efficient, robust and positively conservative. The performance of the new solver is evaluated through numerical testing in one and two dimensions.

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

We present an extension of the Piecewise Parabolic Method to special relativistic fluid dynamics in multidimensions. The scheme is conservative, dimensionally unsplit, and suitable for a general equation of state. Temporal evolution is second-order accurate and employs characteristic projection operators; spatial interpolation is piece-wise parabolic making the scheme third-order accurate in smooth regions of the flow away from discontinuities. The algorithm is written for a general system of orthogonal curvilinear coordinates and can be used for computations in non-cartesian geometries. A non-linear iterative Riemann solver based on the two-shock approximation is used in flux calculation. In this approximation, an initial discontinuity decays into a set of discontinuous waves only implying that, in particular, rarefaction waves are treated as flow discontinuities. We also present a new and simple equation of state which approximates the exact result for the relativistic perfect gas with high accuracy. The strength of the new method is demonstrated in a series of numerical tests and more complex simulations in one, two and three dimensions.

We perform a detailed analytical and numerical study of the dynamics of perturbations (vortex/aperiodic mode, Rossby and spiral-density waves) in 2D compressible disks with a Keplerian law of rotation. We draw attention to the process of spiral-density wave generation from vortices, discussing, in particular, the initial, most peculiar stages of wave emission. We show that the linear phenomenon of wave generation by vortices in smooth (without inflection points) shear flows found by using the so-called non-modal approach, is directly applicable to the present case. After an analytical non-modal description of the physics and characteristics of the spiral-density wave generation/propagation in the local shearing-sheet model, we follow the process of wave generation by small amplitude coherent circular vortex structures, by direct global numerical simulation, describing the main features of the generated waves.

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.

Recently it was found that helical magnetized flows efficiently amplify Alfv\'en waves (Rogava et al. 2003, A&A, v.399, p.421). This robust and manifold nonmodal effect was found to involve regimes of transient algebraic growth (for purely ejectional flows), and exponential instabilities of both usual and parametric nature. However the study was made in the incompressible limit and an important question remained open - whether this amplification is inherent to swirling MHD flows per se and what is the degree of its dependence on the incompressibility condition. In this paper, in order to clear up this important question, we consider full compressible spectrum of MHD modes: Alfv\'en waves (AW), slow magnetosonic waves (SMW) and fast magnetosonic waves (FMW). We find that helical flows inseparably blend these waves with each other and make them unstable, creating the efficient energy transfer from the mean flow to the waves. The possible role of these instabilities for the onset of the MHD turbulence, self-heating of the flow and the overall dynamics of astrophysical flows are discussed.

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.

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 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 minimum value of the power in kinetic bulk motion of the galactic superluminal source GRS 1915+105. This value far exceeds the Eddington luminosity for accretion onto a black hole of 10 solar masses. This large value severely limits the possible carriers of the kinetic luminosity at the base of the jet, and favours a jet production and acceleration controlled by a magnetic field whose value, at the base of the jet, exceeds $10^8$ Gauss. The Blandford and Znajek process can be responsible of the extraction of the rotational energy of a Kerr black hole, if lasting long enough to provide the required kinetic energy. This time, of the order of a day, implies that the process must operate in a stationary, not impulsive, mode.

We present the results of the new optical polarimetric observations of the superluminal source GRO J1655--40, carried out in July 1997, with the multichannel photopolarimeter of the Torino Observatory, using the 2.15-m telescope of Complejo Astronomico El Leoncito (Argentina). The observed amount of polarization shows (mainly in the $V$ and $I$ bands) oscillations which are consistent with the orbital period of the system. This gives information about physical and geometrical properties of the binary system and confirms that the optical continuum is intrinsically polarized. The origin of the polarized flux is likely to be related to electron scattering of photons from the accretion disk.

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 calculate the flux and size ratio of radio emitting features moving in antiparallel directions at relativistic speeds, taking into account that the pattern we observe may move at a velocity different from the one of the emitting plasma. Both velocities can be determined in sources in which both the approaching and the receding features are observed. This is the case of the galactic superluminal source GRS 1915+105, for which the pattern speed (responsible for the apparent superluminal velocity) was found to be 0.92$c$. For this source we find that the velocity of the plasma is 0.73$c$. The found velocity helps decreasing the estimate of the associated kinetic power. Since intrinsically identical knots will be observed to have different sizes, high resolution radio observations of GRS 1915+105 will be a test for the proposed difference of the pattern and plasma velocities.