Showing 1–49 of 49 results for author: Rossi, P

X-Shaped Radio Galaxies: Probing Jet Evolution, Ambient Medium Dynamics, and Their Intricate Interconnection

Authors: Gourab Giri, Christian Fendt, Kshitij Thorat, Gianluigi Bodo, Paola Rossi

Abstract: This review explores the field of X-shaped radio galaxies (XRGs), a distinctive subset of winged radio sources that are identified by two pairs of jetted lobes which aligned by a significant angle, resulting in an inversion-symmetric structure. These lobes, encompassing active (primary) and passive (secondary) phases, exhibit a diverse range of properties across the multiple frequency bands, posin… ▽ More This review explores the field of X-shaped radio galaxies (XRGs), a distinctive subset of winged radio sources that are identified by two pairs of jetted lobes which aligned by a significant angle, resulting in an inversion-symmetric structure. These lobes, encompassing active (primary) and passive (secondary) phases, exhibit a diverse range of properties across the multiple frequency bands, posing challenges in discerning their formation mechanism. The proposed mechanisms can broadly be categorized into those related either to a triaxial ambient medium, into which the jet propagates, or to a complex, central AGN mechanism, where the jet is generated. The observed characteristics of XRGs as discovered in the most substantial sample to date, challenge the idea that there is universal process at work that produces the individual sources of XRGs. Instead, the observational and numerical results rather imply the absence of an universal model and infer that distinct mechanisms may be at play for the specific sources. By scrutinizing salient and confounding properties, this review intends to propose the potential direction for future research to constrain and constrict individual models applicable to XRGs. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Accepted for publication in "Frontiers in Astronomy and Space Sciences"

The different flavors of extragalactic jets: Magnetized relativistic flows

Authors: Paola Rossi, Gianluigi Bodo, Silvano Massaglia, Alessandro Capetti

Abstract: We perform three-dimensional numerical simulations of magnetized relativistic jets propagating in a uniform density environment in order to study the effect of the entrainment and the consequent deceleration, extending a previous work in which magnetic effects were not present. As in previous papers, our aim is to understand the connection between the jet properties and the resulting Fanaroff-Rile… ▽ More We perform three-dimensional numerical simulations of magnetized relativistic jets propagating in a uniform density environment in order to study the effect of the entrainment and the consequent deceleration, extending a previous work in which magnetic effects were not present. As in previous papers, our aim is to understand the connection between the jet properties and the resulting Fanaroff-Riley classification. We consider jets with different low densities, and therefore low power, and different magnetizations. We find that lower magnetization jets effectively decelerate to sub-relativistic velocities and may then result in an FR~I morphology on larger scales. At the opposite, in the higher magnetization cases, the entrainment and consequent deceleration are substantially reduced. } △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 13 pages, 11 figures, accepted in A&A

FR0 jets and recollimation-induced instabilities

Authors: A. Costa, G. Bodo, F. Tavecchio, P. Rossi, A. Capetti, S. Massaglia, A. Sciaccaluga, R. D. Baldi, G. Giovannini

Abstract: The recently discovered population of faint FR0 radiogalaxies has been interpreted as the extension to low power of the classical FRI sources. Their radio emission appears to be concentrated in very compact (pc-scale) cores, any extended emission is very weak or absent and VLBI observations show that jets are already mildly or sub-relativistic at pc scales. Based on these observational properties… ▽ More The recently discovered population of faint FR0 radiogalaxies has been interpreted as the extension to low power of the classical FRI sources. Their radio emission appears to be concentrated in very compact (pc-scale) cores, any extended emission is very weak or absent and VLBI observations show that jets are already mildly or sub-relativistic at pc scales. Based on these observational properties we propose here that the jets of FR0s are strongly decelerated and disturbed at pc scale by hydrodynamical instabilities.}{With the above scenario in mind, we study the dynamics of a low-power relativistic jet propagating into a confining external medium, focusing on the effects of entrainment and mixing promoted by the instabilities developing at the jet-environment interface downstream of a recollimation shock. We perform a 3D relativistic hydrodynamical simulation of a recollimated jet by means of the state-of-the-art code PLUTO. The jet is initially conical, relativistic (with initial Lorentz Factor $Γ$=5), cold and light with respect to the confining medium, whose pressure is assumed to slowly decline with distance. The magnetic field is assumed to be dynamically unimportant. The 3D simulation shows that, after the first recollimation/reflection shock system, a rapidlygrowing instability develops, as a result of the interplay between recollimation-induced instabilities and Richtmyer-Meshkov modes. In turn, the instabilities promote strong mixing and entrainment that rapidly lead to the deceleration of the jet and spread its momentum to slowly moving, highly turbulent external gas. We argue that this mechanism could account for the peculiarities of the low-power FR0 jets. For outflows with higher power, Lorentz factor or magnetic field, we expect that the destabilizing effects are less effective, allowing the survival of the jet up to the kpc scale, as observed in FRIs. △ Less

Submitted 23 January, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 8 pages, 6 figures, accepted for publication in A&A

A polarization study of jets interacting with turbulent magnetic fields

Authors: Moun Meenakshi, Dipanjan Mukherjee, Gianluigi Bodo, Paola Rossi

Abstract: We investigate the effect of the jet's immediate surroundings on the non-thermal synchrotron emission and its polarization properties. The ambient medium is equipped with a turbulent magnetic field, which is compressed and amplified by the jets as they progress. This leads to high polarization at the forward shock surface. The randomness in the magnetic polarities of the external fields in the sho… ▽ More We investigate the effect of the jet's immediate surroundings on the non-thermal synchrotron emission and its polarization properties. The ambient medium is equipped with a turbulent magnetic field, which is compressed and amplified by the jets as they progress. This leads to high polarization at the forward shock surface. The randomness in the magnetic polarities of the external fields in the shocked ambient medium (SAM) results in vector cancellation of the polarized components from the jet, thereby causing depolarization of the radiation from the cocoon. We find that due to the slow decay of the fields in the SAM, such depolarization by the fields with large correlation lengths is more prominent when compared to the small-scale fields. Also, the low-power jets, which have magnetic fields comparable in strength to those in the SAM, are more severely affected by the SAM's depolarizing effect, than the high-power ones. The turbulent backflows in the cocoon, as well as the shearing of fields near the contact discontinuity, strengthen the poloidal component in the jet. This causes internal depolarization due to the cancellation of the orthogonally polarized components along the Line of Sight as the field transitions from ordered toroidal to poloidal. The synchrotron maps display high-emission filaments in the cocoon with magnetic fields aligned along them. The kink instability leads to the wiggling motion of the jet's spine, resulting in hotspot complexes in low-power sources. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS

Magnetocentrifugal mechanism of pair creation in AGN

Authors: Zaza N. Osmanov, Gianluigi Bodo, Paola Rossi

Abstract: In the manuscript, we study the efficiency of pair creation by means of the centrifugal mechanism. The strong magnetic field and the effects of rotation, which always take place in Kerr-type black holes, guarantee the frozen-in condition, leading to the generation of an exponentially amplifying electrostatic field. This field, when reaching the Schwinger threshold, leads to efficient pair producti… ▽ More In the manuscript, we study the efficiency of pair creation by means of the centrifugal mechanism. The strong magnetic field and the effects of rotation, which always take place in Kerr-type black holes, guarantee the frozen-in condition, leading to the generation of an exponentially amplifying electrostatic field. This field, when reaching the Schwinger threshold, leads to efficient pair production. The process has been studied for a wide range of AGN luminosities and black hole masses, and it was found that the mechanism is very efficient, indicating that for AGNs where centrifugal effects are significant, the annihilation lines in the MeV range will be very strong. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 15 pages, 5 figures

Modelling X-shaped Radio Galaxies: Dynamical and Emission Signatures from the Back-flow model

Authors: Gourab Giri, Bhargav Vaidya, Paola Rossi, Gianluigi Bodo, Dipanjan Mukherjee, Andrea Mignone

Abstract: Some of the radio galaxies show distortion in their jets, forming tailed or winged sources. X-shaped radio galaxies are a subclass of winged sources formation mechanism of which is still unclear. The focus of this work is to understand hydro-dynamical back-flows and their role in dynamics and non-thermal emission signatures (in presence of radiative losses and diffusive shock acceleration) during… ▽ More Some of the radio galaxies show distortion in their jets, forming tailed or winged sources. X-shaped radio galaxies are a subclass of winged sources formation mechanism of which is still unclear. The focus of this work is to understand hydro-dynamical back-flows and their role in dynamics and non-thermal emission signatures (in presence of radiative losses and diffusive shock acceleration) during the initial phase of these galaxies. We have performed relativistic MHD simulations of an under-dense jet travelling in a tri-axial ambient using a hybrid Eulerian-Lagrangian framework to incorporate effects of micro-physical processes. We have demonstrated the dominant role played by pressure gradient in shaping XRGs in thermally dominated cases. We show that the prominence of the formed structure decreases as the jet deviates from the major axis of the ambient. The wing evolution is mainly governed by re-energized particles due to shocks that keep the structure active during the evolution time. The synthetic intensity maps of the radio galaxy show similarities with morphologies that are typically found in observed XRGs. This includes the cases with wider wings than the active lobes. The characteristic emission signatures in terms of its synchrotron spectra and implication of equipartition condition in age estimation are also discussed here. Additionally, we show that discrepancy of age can be attributed to mixing of different aged particle populations. Further, the effect of viewing angle on the difference of spectral index of the active lobes and the wings $(Δα)$ shows a large variation and degenerate behaviour. We have demonstrated the role of diffusive shocks in the obtained variation and have concluded that the spread of $(Δα)$ is not a dependable characteristic in determining the formation model of XRGs. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 15 pages, 17 figures, 3 tables; Accepted for publication in Astronomy & Astrophysics journal

Journal ref: A&A 662, A5 (2022)

Making Fanaroff-Riley I radio sources III. The effects of the magnetic field on relativistic jets' propagation and source morphologies

Authors: S. Massaglia, G. Bodo, P. Rossi, A. Capetti, A. Mignone

Abstract: Extragalactic radio sources appear under different morphologies, the most frequent ones are classified as Fanaroff-Riley type I (FR I), typically with lower luminosities, and Fanaroff-Riley type II, (FR II), typically more luminous. This simple classification, however, has many exceptions that we intend to investigate. Following previous analyses in the three-dimensional Hydrodynamic and Magneto-H… ▽ More Extragalactic radio sources appear under different morphologies, the most frequent ones are classified as Fanaroff-Riley type I (FR I), typically with lower luminosities, and Fanaroff-Riley type II, (FR II), typically more luminous. This simple classification, however, has many exceptions that we intend to investigate. Following previous analyses in the three-dimensional Hydrodynamic and Magneto-Hydrodynamic limits, we extend the numerical investigation to the Relativistic Magneto-Hydrodynamic regime, to include sources whose jet kinetic power sets in the range that separates FR Is from FR IIs. We consider weakly and mildly relativistic, underdense, supersonic jets that propagate in a stratified medium. In the model, the ambient temperature increases with distance from the jet origin maintaining constant pressure. We present three cases with low, high and intermediate kinetic luminosity that evolve into different morphologies. We find that the resulting morphology can be highly time dependent and that, apart from the jet power, the jet-to-ambient density ratio and the magnetization parameter play a crucial role in the jet evolution as well. △ Less

Submitted 24 December, 2021; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted for publication on Astronomy & Astrophysics

Report number: Number DF22

Journal ref: A&A 659, A139 (2022)

Current driven kink instabilities in relativistic jets: dissipation properties

Authors: G. Bodo, G. Mamatsashvili, P. Rossi, A. Mignone

Abstract: We analyze the evolution of current driven kink instabilities of a highly magnetized relativistic plasma column, focusing in particular on its dissipation properties. The instability evolution leads to the formation of thin current sheets where the magnetic energy is dissipated. We find that the total amount of dissipated magnetic energy is independent of the dissipation properties. Dissipation oc… ▽ More We analyze the evolution of current driven kink instabilities of a highly magnetized relativistic plasma column, focusing in particular on its dissipation properties. The instability evolution leads to the formation of thin current sheets where the magnetic energy is dissipated. We find that the total amount of dissipated magnetic energy is independent of the dissipation properties. Dissipation occurs in two stages: a peak when the instability saturates, which is characterized by the formation of a helicoidal current sheet at the boundary of the deformed plasma column, followed by a weaker almost flat phase, in which turbulence develops. The detailed properties of these two phases depend on the equilibrium configuration and other parameters, in particular on the steepness of the pitch radial profile, on the presence of an external axial magnetic field and on the amount of magnetization. These results are relevant for high energy astrophysical sources, since current sheets can be the sites of magnetic reconnection where particles can be accelerated to relativistic energies and give rise to the observed radiation. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 16 pages, 22 figures, MNRAS accepted

Simulating the dynamics and synchrotron emission from relativistic jets II. Evolution of non-thermal electrons

Authors: Dipanjan Mukherjee, Gianluigi Bodo, Paola Rossi, Andrea Mignone, Bhargav Vaidya

Abstract: We have simulated the evolution of non-thermal cosmic ray electrons (CREs) in 3D relativistic magneto hydrodynamic (MHD) jets evolved up to a height of 9 kpc. The CREs have been evolved in space and in energy concurrently with the relativistic jet fluid, duly accounting for radiative losses and acceleration at shocks. We show that jets stable to MHD instabilities show expected trends of regular fl… ▽ More We have simulated the evolution of non-thermal cosmic ray electrons (CREs) in 3D relativistic magneto hydrodynamic (MHD) jets evolved up to a height of 9 kpc. The CREs have been evolved in space and in energy concurrently with the relativistic jet fluid, duly accounting for radiative losses and acceleration at shocks. We show that jets stable to MHD instabilities show expected trends of regular flow of CREs in the jet spine and acceleration at a hotspot followed by a settling backflow. However, unstable jets create complex shock structures at the jet-head (kink instability), the jet spine-cocoon interface and the cocoon itself (Kelvin-Helmholtz modes). CREs after exiting jet-head undergo further shock crossings in such scenarios and are re-accelerated in the cocoon. CREs with different trajectories in turbulent cocoons have different evolutionary history with different spectral parameters. Thus at the same spatial location, there is mixing of different CRE populations, resulting in a complex total CRE spectrum when averaged over a given area. Cocoons of unstable jets can have an excess build up of energetic electrons due to re-acceleration at turbulence driven shocks and slowed expansion of the decelerated jet. This will add to the non-thermal energy budget of the cocoon. △ Less

Submitted 22 October, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Comments: Published in MNRAS. arxiv V1 has 2 updated references

Multiple subglacial water bodies below the south pole of Mars unveiled by new MARSIS data

Authors: Sebastian Emanuel Lauro, Elena Pettinelli, Graziella Caprarelli, Luca Guallini, Angelo Pio Rossi, Elisabetta Mattei, Barbara Cosciotti, Andrea Cicchetti, Francesco Soldovieri, Marco Cartacci, Federico Di Paolo, Raffaella Noschese, Roberto Orosei

Abstract: The detection of liquid water by the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) at the base of the south polar layered deposits in Ultimi Scopuli has reinvigorated the debate about the origin and stability of liquid water under present-day Martian conditions. To establish the extent of subglacial water in this region, we acquired new data, achieving extended radar coverage… ▽ More The detection of liquid water by the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) at the base of the south polar layered deposits in Ultimi Scopuli has reinvigorated the debate about the origin and stability of liquid water under present-day Martian conditions. To establish the extent of subglacial water in this region, we acquired new data, achieving extended radar coverage over the study area. Here, we present and discuss the results obtained by a new method of analysis of the complete MARSIS dataset, based on signal processing procedures usually applied to terrestrial polar ice sheets. Our results strengthen the claim of the detection of a liquid water body at Ultimi Scopuli and indicate the presence of other wet areas nearby. We suggest that the waters are hypersaline perchlorate brines, known to form at Martian polar regions and thought to survive for an extended period of time on a geological scale at below-eutectic temperatures. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Numerical study of Kelvin-Helmholtz instability and its impact on synthetic emission from magnetized jets

Authors: Nikhil Borse, Sriyasriti Acharya, Bhargav Vaidya, Dipanjan Mukherjee, Gianluigi Bodo, Paola Rossi, Andrea Mignone

Abstract: Non-thermal emission from Active Galactic Nuclei (AGN) jets extends up-to large scales in-spite of them being prone to a slew of magneto-hydrodynamic instabilities. The main focus of this study is to understand the impact of MHD instabilities on the non-thermal emission from large-scale AGN jets. We perform high-resolution three-dimensional numerical magneto-hydrodynamic simulations of a plasma co… ▽ More Non-thermal emission from Active Galactic Nuclei (AGN) jets extends up-to large scales in-spite of them being prone to a slew of magneto-hydrodynamic instabilities. The main focus of this study is to understand the impact of MHD instabilities on the non-thermal emission from large-scale AGN jets. We perform high-resolution three-dimensional numerical magneto-hydrodynamic simulations of a plasma column to investigate the dynamical and emission properties of jet configurations at kilo-parsec scales with different magnetic field profiles, jet speeds, and density contrast. We also obtain synthetic non-thermal emission signatures for different viewing angles using an approach that assumes static particle spectra and that obtained by evolving the particle spectra using Lagrangian macro-particles incorporating the effects of shock acceleration and radiative losses. We find that the shocks due to Kelvin-Helmholtz (KH) instability in the axial magnetic field configurations can strongly affect the jet dynamics. Additionally, we also find the presence of weak biconical shocks in the under-dense jet columns. The inclusion of a helical magnetic field hinders the vortex growth at the shear surface thereby stabilizing the jet column. With the evolving particle spectra approach, the synthetic SEDs obtained for cases with strong KH instability show the presence of multiple humps ranging from radio to TeV gamma-ray band. We conclude that the high-energy electrons accelerated in the vicinity of freshly formed shocks due to KH instability, result in high X-ray emission. △ Less

Submitted 11 March, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: 18 pages, 12 figures, Accepted for publication in A&A

Journal ref: A&A 649, A150 (2021)

Simulating the dynamics and non-thermal emission of relativistic magnetised jets I. Dynamics

Authors: Dipanjan Mukherjee, Gianluigi Bodo, Andrea Mignone, Paola Rossi, Bhargav Vaidya

Abstract: We have performed magneto-hydrodynamic simulations of relativistic jets from supermassive blackholes over a few tens of kpc for a range of jet parameters. One of the primary aims were to investigate the effect of different MHD instabilities on the jet dynamics and their dependence on the choice of jet parameters. We find that two dominant MHD instabilities affect the dynamics of the jet, small sca… ▽ More We have performed magneto-hydrodynamic simulations of relativistic jets from supermassive blackholes over a few tens of kpc for a range of jet parameters. One of the primary aims were to investigate the effect of different MHD instabilities on the jet dynamics and their dependence on the choice of jet parameters. We find that two dominant MHD instabilities affect the dynamics of the jet, small scale Kelvin- Helmholtz (KH) modes and large scale kink modes, whose evolution depend on internal jet parameters like the Lorentz factor, the ratio of the density and pressure to the external medium and the magnetisation and hence consequently on the jet power. Low power jets are susceptible to both instabilities, kink modes for jets with higher central magnetic field and KH modes for lower magnetisation. Moderate power jets do not show appreciable growth of kink modes, but KH modes develop for lower magnetisation. Higher power jets are generally stable to both instabilities. Such instabilities decelerate and decollimate the jet while inducing turbulence in the cocoon, with consequences on the magnetic field structure. We model the dynamics of the jets following a generalised treatment of the Begelman-Cioffi relations which we present here. We find that the dynamics of stable jets match well with simplified analytic models of expansion of non self-similar FRII jets, whereas jets with prominent MHD instabilities show a nearly self-similar evolution of the morphology as the energy is more evenly distributed between the jet head and the cocoon. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS

The different flavors of extragalactic jets: the role of relativistic flow deceleration

Authors: P. Rossi, G. Bodo, S. Massaglia, A. Capetti

Abstract: We perform three-dimensional numerical simulations of relativistic (with a Lorentz factor of 10), non magnetized jets propagating in a uniform density environment, in order to study the effect of the entrainment and the consequent deceleration. Our simulations investigate the jet propagation inside the galaxy core, where, most likely, the deceleration occurs more efficiently. We compare cases with… ▽ More We perform three-dimensional numerical simulations of relativistic (with a Lorentz factor of 10), non magnetized jets propagating in a uniform density environment, in order to study the effect of the entrainment and the consequent deceleration. Our simulations investigate the jet propagation inside the galaxy core, where, most likely, the deceleration occurs more efficiently. We compare cases with different density and pressure ratios with respect to the ambient medium finding that low density jets are efficiently decelerated and reach a quasi steady state in which, over a length of 600 jet radii, slow down from highly to sub-relativistic velocities. At the opposite, denser jets keep highly relativistic velocity over the same length. We discuss these results in relation to the Faranoff Riley (FR) radio-sources classification. We infer that lower density jets can give rise to FR 0 and FR I radio-sources, while higher density jets may be connected to FR II radio-sources. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 11 pages, 16 figures, A&A accepted

Journal ref: A&A 642, A69 (2020)

Virtual European Solar & Planetary Access (VESPA): a Planetary Science Virtual Observatory cornerstone

Authors: S. Erard, B. Cecconi, P. Le Sidaner, C. Chauvin, A. P. Rossi, M. Minin, T. Capria, S. Ivanovski, B. Schmitt, V. Genot, N. Andre, C. Marmo, A. C. Vandaele, L. Trompet, M. Scherf, R. Hueso, A. Maattanen, B. Carry, N. Achilleos, J. Soucek, D. Pisa, K. Benson, P. Fernique, E. Millour

Abstract: The Europlanet-2020 programme, which ended on Aug 31st, 2019, included an activity called VESPA (Virtual European Solar and Planetary Access), which focused on adapting Virtual Observatory (VO) techniques to handle Planetary Science data. This paper describes some aspects of VESPA at the end of this 4-years development phase and at the onset of the newly selected Europlanet-2024 programme starting… ▽ More The Europlanet-2020 programme, which ended on Aug 31st, 2019, included an activity called VESPA (Virtual European Solar and Planetary Access), which focused on adapting Virtual Observatory (VO) techniques to handle Planetary Science data. This paper describes some aspects of VESPA at the end of this 4-years development phase and at the onset of the newly selected Europlanet-2024 programme starting in 2020. The main objectives of VESPA are to facilitate searches both in big archives and in small databases, to enable data analysis by providing simple data access and online visualization functions, and to allow research teams to publish derived data in an interoperable environment as easily as possible. VESPA encompasses a wide scope, including surfaces, atmospheres, magnetospheres and planetary plasmas, small bodies, helio-physics, exoplanets, and spectroscopy in solid phase. This system relies in particular on standards and tools developed for the Astronomy VO (IVOA) and extends them where required to handle specificities of Solar System studies. It also aims at making the VO compatible with tools and protocols developed in different contexts, for instance GIS for planetary surfaces, or time series tools for plasma-related measurements. An essential part of the activity is to publish a significant amount of high-quality data in this system, with a focus on derived products resulting from data analysis or simulations. △ Less

Submitted 27 January, 2020; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: Submitted to Data Science Journal

Linear stability analysis of magnetized relativistic rotating jets

Authors: Gianluigi Bodo, George Mamatsashvili, Paola Rossi, Andrea Mignone

Abstract: 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 weakl… ▽ More 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. △ Less

Submitted 27 February, 2019; originally announced February 2019.

Comments: 15 pages, 15 figures, accepted for pubblication in MNRAS

Fluids mobilization in Arabia Terra, Mars: depth of pressurized reservoir from mounds self-similar clustering

Authors: Riccardo Pozzobon, Francesco Mazzarini, Matteo Massironi, Angelo Pio Rossi, Monica Pondrelli, Gabriele Cremonese, Lucia Marinangeli

Abstract: Arabia Terra is a region of Mars where signs of past-water occurrence are recorded in several landforms. Broad and local scale geomorphological, compositional and hydrological analyses point towards pervasive fluid circulation through time. In this work we focus on mound fields located in the interior of three casters larger than 40 km (Firsoff, Kotido and unnamed crater 20 km to the east) and sho… ▽ More Arabia Terra is a region of Mars where signs of past-water occurrence are recorded in several landforms. Broad and local scale geomorphological, compositional and hydrological analyses point towards pervasive fluid circulation through time. In this work we focus on mound fields located in the interior of three casters larger than 40 km (Firsoff, Kotido and unnamed crater 20 km to the east) and showing strong morphological and textural resemblance to terrestrial mud volcanoes and spring-related features. We infer that these landforms likely testify the presence of a pressurized fluid reservoir at depth and past fluid upwelling. We have performed morphometric analyses to characterize the mound morphologies and consequently retrieve an accurate automated mapping of the mounds within the craters for spatial distribution and fractal clustering analysis. The outcome of the fractal clustering yields information about the possible extent of the percolating fracture network at depth below the craters. We have been able to constrain the depth of the pressurized fluid reservoir between ~2.5 and 3.2 km of depth and hence, we propose that mounds and mounds alignments are most likely associated to the presence of fissure ridges and fluid outflow. Their process of formation is genetically linked to the formation of large intra-crater bulges previously interpreted as large scale spring deposits. The overburden removal caused by the impact crater formation is the inferred triggering mechanism for fluid pressurization and upwelling, that through time led to the formation of the intra-crater bulges and, after compaction and sealing, to the widespread mound fields in their surroundings. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Making Faranoff-Riley I radio sources II. The effects of jet magnetization

Authors: S. Massaglia, G. Bodo, P. Rossi, S. Capetti, A. Mignone

Abstract: 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 m… ▽ More 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. △ Less

Submitted 7 January, 2019; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: In press on Astronomy & Astrophysics

Journal ref: A&A 621, A132 (2019)

A Particle Module for the PLUTO code: II - Hybrid Framework for Modeling Non-thermal emission from Relativistic Magnetized flows

Authors: Bhargav Vaidya, Andrea Mignone, Gianluigi Bodo, Paola Rossi, Silvano Massaglia

Abstract: 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)… ▽ More 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. △ Less

Submitted 27 August, 2018; originally announced August 2018.

Comments: 23 pages, 14 figures, Accepted for publication in The Astrophysical Journal

Magnetic Helicities and Dynamo Action in Magneto-rotationally Driven Turbulence

Authors: G. Bodo, F. Cattaneo, A. Mignone, P. Rossi

Abstract: 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 o… ▽ More 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. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 26 pages, 10 figures, ApJ accepted

Interoperability in Planetary Research for Geospatial Data Analysis

Authors: Trent M. Hare, Angelo P. Rossi, Alessandro Frigeri, Chiara Marmo

Abstract: For more than a decade there has been a push in the planetary science community to support interoperable methods for accessing and working with geospatial data. Common geospatial data products for planetary research include image mosaics, digital elevation or terrain models, geologic maps, geographic location databases (e.g., craters, volcanoes) or any data that can be tied to the surface of a pla… ▽ More For more than a decade there has been a push in the planetary science community to support interoperable methods for accessing and working with geospatial data. Common geospatial data products for planetary research include image mosaics, digital elevation or terrain models, geologic maps, geographic location databases (e.g., craters, volcanoes) or any data that can be tied to the surface of a planetary body (including moons, comets or asteroids). Several U.S. and international cartographic research institutions have converged on mapping standards that embrace standardized geospatial image formats, geologic mapping conventions, U.S. Federal Geographic Data Committee (FGDC) cartographic and metadata standards, and notably on-line mapping services as defined by the Open Geospatial Consortium (OGC). △ Less

Submitted 8 June, 2017; originally announced June 2017.

Comments: 15 pages, 3 figures, from American Geophysical Union 2015 invited talk, In Press

VESPA: a community-driven Virtual Observatory in Planetary Science

Authors: S. Erard, B. Cecconi, P. Le Sidaner, A. P. Rossi, T. Capria, B. Schmitt, V. Génot, N. André, A. C. Vandaele, M. Scherf, R. Hueso, A. Määttänen, W. Thuillot, B. Carry, N. Achilleos, C. Marmo, O. Santolik, K. Benson, P. Fernique, L. Beigbeder, E. Millour, B. Rousseau, F. Andrieu, C. Chauvin, M. Minin , et al. (17 additional authors not shown)

Abstract: The VESPA data access system focuses on applying Virtual Observatory (VO) standards and tools to Planetary Science. Building on a previous EC-funded Europlanet program, it has reached maturity during the first year of a new Europlanet 2020 program (started in 2015 for 4 years). The infrastructure has been upgraded to handle many fields of Solar System studies, with a focus both on users and data p… ▽ More The VESPA data access system focuses on applying Virtual Observatory (VO) standards and tools to Planetary Science. Building on a previous EC-funded Europlanet program, it has reached maturity during the first year of a new Europlanet 2020 program (started in 2015 for 4 years). The infrastructure has been upgraded to handle many fields of Solar System studies, with a focus both on users and data providers. This paper describes the broad lines of the current VESPA infrastructure as seen by a potential user, and provides examples of real use cases in several thematic areas. These use cases are also intended to identify hints for future developments and adaptations of VO tools to Planetary Science. △ Less

Submitted 14 June, 2017; v1 submitted 26 May, 2017; originally announced May 2017.

Comments: Planetary and Space Sciences (in press), Special Issue "Enabling Open and Interoperable Access to Planetary Science and Heliophysics Databases and Tools". 43 pages, 14 figures, 1 table

Journal ref: Planetary and Space Sciences 150, 65-85, 2018

Radar sounding of Lucus Planum, Mars, by MARSIS

Authors: Roberto Orosei, Angelo Pio Rossi, Federico Cantini, Graziella Caprarelli, Lynn M. Carter, Irene Papiano, Marco Cartacci, Andrea Cicchetti, Raffaella Noschese

Abstract: Lucus Planum, extending for a radius of approximately 500 km around 181° E, 5° S, is part of the Medusae Fossae Formation (MFF), a set of several discontinuous deposits of fine-grained, friable material straddling across the Martian highland-lowland boundary. The MFF has been variously hypothesized to consist of pyroclastic flows, pyroclastic airfall, paleopolar deposits, or atmospherically-deposi… ▽ More Lucus Planum, extending for a radius of approximately 500 km around 181° E, 5° S, is part of the Medusae Fossae Formation (MFF), a set of several discontinuous deposits of fine-grained, friable material straddling across the Martian highland-lowland boundary. The MFF has been variously hypothesized to consist of pyroclastic flows, pyroclastic airfall, paleopolar deposits, or atmospherically-deposited icy dust driven by climate cycles. MARSIS, a low-frequency subsurface-sounding radar carried by ESA's Mars Express, acquired 238 radar swaths across Lucus Planum, providing sufficient coverage for the study of its internal structure and dielectric properties. Subsurface reflections were found only in three areas, marked by a distinctive surface morphology, while the central part of Lucus Planum appears to be made of radar-attenuating material preventing the detection of basal echoes. The bulk dielectric properties of these areas were estimated and compared with those of volcanic rocks and ice-dust mixtures. Previous interpretations that east Lucus Planum and the deposits on the north-western flanks of Apollinaris Patera consist of high-porosity pyroclastic material are strongly supported by the new results. The north-western part of Lucus Planum is likely to be much less porous, although interpretations about the nature of the subsurface materials are not conclusive. The exact origin of the deposits cannot be constrained by radar data alone, but our results for east Lucus Planum are consistent with an overall pyroclastic origin, likely linked to Tharsis Hesperian and Amazonian activity. △ Less

Submitted 25 May, 2017; originally announced May 2017.

Comments: 39 pages, 8 figures, 2 tables

3D Relativistic MHD numerical simulations of X-shaped radio sources

Authors: P. Rossi, G. Bodo, A. Capetti, S. Massaglia

Abstract: 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… ▽ More 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. △ Less

Submitted 28 June, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 13 pages, 16 figures, Astronomy and Astrophysics accepted

Journal ref: A&A 606, A57 (2017)

Online characterization of planetary surfaces: PlanetServer, an open-source analysis and visualization tool

Authors: R. Marco Figuera, B. Pham Huu, A. P. Rossi, M. Minin, J. Flahaut, A. Halder

Abstract: The lack of open-source tools for hyperspectral data visualization and analysiscreates a demand for new tools. In this paper we present the new PlanetServer,a set of tools comprising a web Geographic Information System (GIS) and arecently developed Python Application Programming Interface (API) capableof visualizing and analyzing a wide variety of hyperspectral data from differentplanetary bodies.… ▽ More The lack of open-source tools for hyperspectral data visualization and analysiscreates a demand for new tools. In this paper we present the new PlanetServer,a set of tools comprising a web Geographic Information System (GIS) and arecently developed Python Application Programming Interface (API) capableof visualizing and analyzing a wide variety of hyperspectral data from differentplanetary bodies. Current WebGIS open-source tools are evaluated in orderto give an overview and contextualize how PlanetServer can help in this mat-ters. The web client is thoroughly described as well as the datasets availablein PlanetServer. Also, the Python API is described and exposed the reason ofits development. Two different examples of mineral characterization of differenthydrosilicates such as chlorites, prehnites and kaolinites in the Nili Fossae areaon Mars are presented. As the obtained results show positive outcome in hyper-spectral analysis and visualization compared to previous literature, we suggestusing the PlanetServer approach for such investigations. △ Less

Submitted 29 April, 2017; v1 submitted 6 January, 2017; originally announced January 2017.

A strategy for an accurate estimation of the basal permittivity in the Martian North Polar Layered Deposits

Authors: S. E. Lauro, G. Gennarelli, E. Pettinelli, F. Soldovieri, F. Cantini, A. P. Rossi, R. Orosei

Abstract: The paper deals with the investigation of the Mars subsurface by means of data collected by the Mars Advanced Radar for Subsurface and Ionosphere Sounding working at few MHz frequencies. A data processing strategy, which combines a simple inversion model and an accurate procedure for data selection is presented. This strategy permits to mitigate the theoretical and practical difficulties of the in… ▽ More The paper deals with the investigation of the Mars subsurface by means of data collected by the Mars Advanced Radar for Subsurface and Ionosphere Sounding working at few MHz frequencies. A data processing strategy, which combines a simple inversion model and an accurate procedure for data selection is presented. This strategy permits to mitigate the theoretical and practical difficulties of the inverse problem arising because of the inaccurate knowledge of the parameters regarding both the scenario under investigation and the radiated electromagnetic field impinging on the Mars surface. The results presented in this paper show that, it is possible to reliably retrieve the electromagnetic properties of deeper structures, if such strategy is accurately applied. An example is given here, where the analysis of the data collected on Gemina Lingula, a region of the North Polar layer deposits, allowed us to retrieve permittivity values for the basal unit in agreement with those usually associated to the Earth basaltic rocks. △ Less

Submitted 27 September, 2016; originally announced September 2016.

Journal ref: Geophysical Prospecting 2016

Making Fanaroff-Riley I radio sources. Numerical Hydrodynamic 3D Simulations of Low Power Jets

Authors: S. Massaglia, G. Bodo, P. Rossi, S. Capetti, A. Mignone

Abstract: 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 t… ▽ More 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. △ Less

Submitted 18 October, 2016; v1 submitted 8 September, 2016; originally announced September 2016.

Comments: 11 pages, 12 figures, to appear on A&A

Linear stability analysis of magnetized jets: the rotating case

Authors: G. Bodo, G. Mamatsashvili, P. Rossi, A. Mignone

Abstract: 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én velocity. Rota… ▽ More 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én 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. △ Less

Submitted 6 July, 2016; originally announced July 2016.

Comments: 25 pages, 19 figures, accepted in MNRAS

Global Properties of Fully Convective Accretion Disks from Local Simulations

Authors: G. Bodo, F. Cattaneo, A. Mignone, F. Ponzo, P. Rossi

Abstract: 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… ▽ More 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. △ Less

Submitted 19 June, 2015; originally announced February 2016.

Comments: ApJ accepted, 12 pages, 5 figures

Origin of X-shaped radio-sources: further insights from the properties of their host galaxies

Authors: Melissa Gillone, Alessandro Capetti, Paola Rossi

Abstract: We analyze the properties of a sample of X-shaped radio-sources (XRSs). These objects show, in addition to the main lobes, a pair of wings producing their peculiar radio morphology. We obtain our sample by selecting from the initial list of Cheung (2007, AJ, 133, 2097) the 53 galaxies with the better defined wings and with available SDSS images. We identified the host galaxies and measured their o… ▽ More We analyze the properties of a sample of X-shaped radio-sources (XRSs). These objects show, in addition to the main lobes, a pair of wings producing their peculiar radio morphology. We obtain our sample by selecting from the initial list of Cheung (2007, AJ, 133, 2097) the 53 galaxies with the better defined wings and with available SDSS images. We identified the host galaxies and measured their optical position angle, obtaining a positive result in 22 cases. The orientation of the secondary radio structures shows a strong connection with the optical axis, with all (but one) wing forming a angle larger than 40 degrees with the host major axis. The probability that this is compatible with a uniform distribution is P = 0.9 10E-4. Spectra are available from the SDSS for 28 XRSs. We modeled them to extract information on their emission lines and stellar population properties. The sample is formed by approximately the same number of high and low excitation galaxies (HEG and LEG); this classification is essential for a proper comparison with non-winged radio-galaxies. XRSs follow the same relations between radio and line luminosity defined by radio-galaxies in the 3C sample. While in HEGs a young stellar population is often present, this is not detected in the 13 LEGs, again in agreement with the properties of non XRSs. The lack of young stars in LEGs support the idea that they did not experiences a recent gas rich merger. The connection between the optical axis and the wings orientation, as well as the stellar population and emission lines properties, provide further support for an hydro-dynamic origin of the radio-wings (for example associated with the expansion of the radio cocoon in an asymmetric external medium) rather than with a change of orientation of the jet axis. △ Less

Submitted 18 December, 2015; originally announced December 2015.

Comments: Accepted for publication in A&A

Linear and nonlinear evolution of current-carrying highly magnetized jets

Authors: M. Anjiri, A. Mignone, G. Bodo, P. Rossi

Abstract: 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) insta… ▽ More 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. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Comments: 14 pages, 11 figures

On the convergence of Magnetorotational turbulence in stratified isothermal shearing boxes

Authors: G. Bodo, F. Cattaneo, A. Mignone, P. Rossi

Abstract: 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 converge… ▽ More 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. △ Less

Submitted 24 April, 2014; originally announced April 2014.

Comments: 10 pages, 4 figures, accepted for publication in ApJ Letters

Revisiting linear dynamics of non-axisymmetric perturbations in weakly magnetized accretion discs

Authors: G. R. Mamatsashvili, G. D. Chagelishvili, G. Bodo, P. Rossi

Abstract: 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… ▽ More 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. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Linear stability analysis of magnetized relativistic jets: the nonrotating case

Authors: G. Bodo, G. Mamatsashvili, P. Rossi, A. Mignone

Abstract: 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 curre… ▽ More 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. △ Less

Submitted 24 July, 2013; originally announced July 2013.

Comments: 22 pages, 13 figures, MNRAS in press

Fully Convective Magnetorotational Turbulence in Stratified Shearing Boxes

Authors: G. Bodo, F. Cattaneo, A. Mignone, P. Rossi

Abstract: 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 corresp… ▽ More 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. △ Less

Submitted 14 June, 2013; originally announced June 2013.

Comments: 11 pages, 4 figures, accepted for publication in ApJ Letters

Magnetorotational turbulence in stratified shearing boxes with perfect gas equation of state and finite thermal diffusivity

Authors: G. Bodo, F. Cattaneo, A. Mignone, P. Rossi

Abstract: 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 th… ▽ More 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. △ Less

Submitted 24 October, 2012; originally announced October 2012.

Comments: 18 pages, 10 figures, ApJ accepted

Symmetries, scaling laws and convergence in shearing-box simulations of MRI driven turbulence

Authors: Gianluigi Bodo, Fausto Cattaneo, Attilio Ferrari, Andrea Mignone, Paola Rossi

Abstract: 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 a… ▽ More 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. △ Less

Submitted 28 June, 2011; originally announced June 2011.

Comments: 28 pages, 7 figures, submitted to ApJ

Linear coupling of modes in 2D radially stratified astrophysical discs

Authors: A. G. Tevzadze, G. D. Chagelishvili, G. Bodo, P. Rossi

Abstract: 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 ex… ▽ More 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. △ Less

Submitted 15 September, 2009; originally announced September 2009.

Comments: 10 pages, 10 figures, MNRAS (accepted)

High Resolution 3D Relativistic MHD Simulations of Jets

Authors: A. Mignone, P. Rossi, G. Bodo, A. Ferrari, S. Massaglia

Abstract: 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-resolutio… ▽ More 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. △ Less

Submitted 31 August, 2009; originally announced August 2009.

Comments: 7 pages, 5 figures

Formation of dynamical structures in relativistic jets: the FRI case

Authors: P. Rossi, A. Mignone, G. Bodo, S. Massaglia, A. Ferrari

Abstract: 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… ▽ More 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. △ Less

Submitted 10 June, 2008; originally announced June 2008.

Comments: 20 pages, 17 figures, A&A accepted

Aspect ratio dependence in magnetorotational instability shearing box simulations

Authors: G. Bodo, A. Mignone, F. Cattaneo, P. Rossi, A. Ferrari

Abstract: 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 i… ▽ More 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. △ Less

Submitted 8 May, 2008; originally announced May 2008.

Comments: 5 pages, 6 figures, Astronomy and Astrophysics accepted

Stability and nonlinear adjustment of vortices in Keplerian flows

Authors: G. Bodo, A. Tevzadze, G. Chagelishvili, A. Mignone, P. Rossi, A. Ferrari

Abstract: 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… ▽ More 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. △ Less

Submitted 25 August, 2007; originally announced August 2007.

Comments: 12 pages, 10 figures

Spiral density wave generation by vortices in Keplerian flows

Authors: G. Bodo, G. Chagelishvili, G. Murante, A. Tevzadze, P. Rossi, A. Ferrari

Abstract: 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… ▽ More 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. △ Less

Submitted 22 March, 2005; originally announced March 2005.

Comments: 18 pages, 16 figures, Astronomy & Astrophysics in press

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

Authors: C. Zanni, G. Murante, G. Bodo, S. Massaglia, P. Rossi, A. Ferrari

Abstract: 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… ▽ More 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. △ Less

Submitted 30 August, 2004; originally announced August 2004.

Comments: 18 pages, 14 figures, to appear on Astronomy and Astrophysics

Journal ref: Astron.Astrophys. 429 (2005) 399-415

X-ray emission from expanding cocoons

Authors: C. Zanni, G. Bodo, P. Rossi, S. Massaglia, A. Durbala, A. Ferrari

Abstract: 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… ▽ More 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. △ Less

Submitted 14 February, 2003; originally announced February 2003.

Comments: 15 pages, 12 figures, to appear on Astronomy and Astrophysics

Journal ref: Astron.Astrophys. 402 (2003) 949-962

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

Authors: A. Capetti, S. Zamfir, P. Rossi, G. Bodo, C. Zanni, S. Massaglia

Abstract: 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… ▽ More 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. △ Less

Submitted 16 July, 2002; originally announced July 2002.

Comments: 7 pages, 6 figures, 1 table, accepted for publication in A&A

Journal ref: Astron.Astrophys. 394 (2002) 39-46

On the matter entrainment by stellar jets and the acceleration of molecular outflows

Authors: M. Micono, G. Bodo, S. Massaglia, P. Rossi, A. Ferrari

Abstract: 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 the… ▽ More 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. △ Less

Submitted 6 September, 2000; originally announced September 2000.

Comments: 9 pages containing 10 postscript figures; uses A&A macros. To appear in Astronomy and Astrophysics

Kelvin-Helmholtz instability in three dimensional radiative jets

Authors: M. Micono, G. Bodo, S. Massaglia, P. Rossi, A. Ferrari, R. Rosner

Abstract: 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… ▽ More 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. △ Less

Submitted 15 June, 2000; originally announced June 2000.

Comments: 15 pages containing 13 postscript figures; uses A&A macros. To appear in Astronomy and Astrophysics

Understanding NLR in Seyfert Galaxies: numerical simulation of jet-cloud interaction

Authors: P. Rossi, A. Capetti, G. Bodo, S. Massaglia, A. Ferrari

Abstract: 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 m… ▽ More 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. △ Less

Submitted 1 February, 2000; originally announced February 2000.

Comments: 11 pages, accepted A&A

Acoustic Waves in a Stratified Atmosphere: II. Three-Dimensional Hydrodynamics

Authors: G. Bodo, W. Kalkofen, S. Massaglia, P. Rossi

Abstract: 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… ▽ More 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. △ Less

Submitted 17 November, 1999; originally announced November 1999.

Comments: 10 pages, 11 figures to be published in A&A