Showing 1–50 of 101 results for author: Mignone, A

A Fourth-Order Finite Volume Scheme for Resistive Relativistic Magnetohydrodynamics

Authors: Andrea Mignone, Vittoria Berta, Marco Rossazza, Matteo Bugli, Giancarlo Mattia, Luca Del Zanna, Lorenzo Pareschi

Abstract: We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advance… ▽ More We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advanced as cell volume-averages. Spatial accuracy employs 5th-order accurate WENO-Z reconstruction from point values (as described in a companion paper) to obtain left and right states at zone interfaces. Explicit flux evaluation is carried out by solving a Riemann problem at cell interfaces, using the Maxwell-Harten-Lax-van Leer with contact wave resolution (MHLLC). Time stepping is based on the implicit-explicit (IMEX) Runge-Kutta (RK) methods, of which we consider both the 3rd-order strong stability preserving SSP3(4,3,3) and a recent 4th-order additive RK scheme, to cope with the stiffness introduced by the source term in Ampere's law. Numerical benchmarks are presented in order to assess the accuracy and robustness of our implementation. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 17 pages, 8 figures - Accepted 2024 July 6. Received 2024 July 3; in original form 2024 March 19

Probing the role of self-gravity in clouds impacted by AGN-driven winds

Authors: Ankush Mandal, Dipanjan Mukherjee, Christoph Federrath, Geoffrey V. Bicknell, Nicole P. H. Nesvadba, Andrea Mignone

Abstract: The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of… ▽ More The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of self-gravity. Our results demonstrate that the cloudlets fragmented by the wind can become gravitationally bound, significantly increasing their survival time. While external pressurisation leads to a global collapse of the clouds in cases of weaker winds ($10^{42}-10^{43}~{\rm erg~s^{-1}}$), higher-power winds ($10^{44}-10^{45}~{\rm erg~s^{-1}}$) disperse the gas and cause localised collapse of the cloudlets. We also demonstrate that a kinetic energy-dominated wind is more efficient in accelerating and dispersing the gas than a thermal wind with the same power. The interaction can give rise to multi-phase outflows with velocities ranging from a few 100 to several 1000~${\rm km\,s^{-1}}$. The mass outflow rates are tightly correlated with the wind power, which we explain by an ablation-based mass-loss model. Moreover, the velocity dispersion and the virial parameter of the cloud material can increase by up to one order of magnitude through the effect of the wind. Even though the wind can suppress or quench star formation for about 1 Myr during the initial interaction, a substantial number of gravitationally bound dense cloudlets manage to shield themselves from the wind's influence and subsequently undergo rapid gravitational collapse, leading to an enhanced star formation rate (SFR). △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. Accepted 2024 May 16. Received 2024 May 16; in original form 2024 January 18

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Authors: Asmita Bhandare, Benoît Commerçon, Guillaume Laibe, Mario Flock, Rolf Kuiper, Thomas Henning, Andrea Mignone, Gabriel-Dominique Marleau

Abstract: Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. W… ▽ More Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. We explore if thermal and dynamical conditions developed during disc formation can generate gas flows that efficiently mix and transport well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a 2D axisymmetric geometry and follow the azimuthal component of velocity. Dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10 and 100 micron-sized neutral spherical dust. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 Msun and explore effects of initial rotation rates and cloud core sizes. Our study underlines mechanisms for early transport of dust from inner hot disc regions via the occurrence of meridional flows and outflow. The vortical flow fosters dynamical mixing and retention of dust while thermal pressure driven outflow replenishes dust in the outer disc. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse. [Abridged] △ Less

Submitted 12 May, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: Language edited version

A $4^{\rm th}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions

Authors: Vittoria Berta, Andrea Mignone, Matteo Bugli, Giancarlo Mattia

Abstract: We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exp… ▽ More We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than one-dimensional finite volume ones), employs the generic upwind constrained transport averaging and sophisticated limiting strategies that include both a discontinuity detector and an order reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness of the method. △ Less

Submitted 18 December, 2023; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 30 pages, 14 figures

Journal ref: Journal of Computational Physics, 499, 112701 (2024)

Resistive relativistic MHD simulations of astrophysical jets

Authors: Giancarlo Mattia, Luca Del Zanna, Matteo Bugli, Andrea Pavan, Riccardo Ciolfi, Gianluigi Bodo, Andrea Mignone

Abstract: Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics (RRMHD) simulations. We aim at investigating different values and models for the plasma resistivity coefficient, and at assessing their impact on the level of… ▽ More Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics (RRMHD) simulations. We aim at investigating different values and models for the plasma resistivity coefficient, and at assessing their impact on the level of turbulence, the formation of current sheets and reconnection plasmoids, the electromagnetic energy content, and the dissipated power. Methods. We use the PLUTO code for simulations and we assume an axisymmetric setup for jets, endowed with both poloidal and toroidal magnetic fields, and propagating in a uniform magnetized medium. The gas is assumed to be characterized by a realistic Synge-like equation of state (Taub equation), appropriate for such type of astrophysical jets. The Taub equation is combined here for the first time with the Implicit-Explicit Runge-Kutta time-stepping procedure, as required in RRMHD simulations. Results. The main result is that turbulence is clearly suppressed for the highest values of resistivity (low Lundquist numbers), current sheets are broader, and plasmoids are barely present, while for low values of resistivity results are very similar to ideal runs, where dissipation is purely numerical. We find that recipes employing a variable resistivity based on the advection of a jet tracer or on the assumption of a uniform Lundquist number improve on the use of a constant coefficient and are probably more realistic, preserving the development of turbulence and of sharp current sheets, possible sites for the acceleration of the non-thermal particles producing the observed high-energy emission. △ Less

Submitted 12 September, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 17 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 679, A49 (2023)

A self-gravity module for the PLUTO code

Authors: Ankush Mandal, Dipanjan Mukherjee, Andrea Mignone

Abstract: We present a novel implementation of an iterative solver for the solution of the Poisson equation in the PLUTO code for astrophysical fluid dynamics. Our solver relies on a relaxation method in which convergence is sought as the steady-state solution of a parabolic equation, whose time-discretization is governed by the \textit{Runge-Kutta-Legendre} (RKL) method. Our findings indicate that the RKL-… ▽ More We present a novel implementation of an iterative solver for the solution of the Poisson equation in the PLUTO code for astrophysical fluid dynamics. Our solver relies on a relaxation method in which convergence is sought as the steady-state solution of a parabolic equation, whose time-discretization is governed by the \textit{Runge-Kutta-Legendre} (RKL) method. Our findings indicate that the RKL-based Poisson solver, which is both fully parallel and rapidly convergent, has the potential to serve as a practical alternative to conventional iterative solvers such as the \textit{Gauss-Seidel} (GS) and \textit{successive over-relaxation} (SOR) methods. Additionally, it can mitigate some of the drawbacks of these traditional techniques. We incorporate our algorithm into a multigrid solver to provide a simple and efficient gravity solver that can be used to obtain the gravitational potentials in self-gravitational hydrodynamics. We test our implementation against a broad range of standard self-gravitating astrophysical problems designed to examine different aspects of the code. We demonstrate that the results match excellently with the analytical predictions (when available), and the findings of similar previous studies. △ Less

Submitted 14 September, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Published in ApJS

Journal ref: ApJS 268 40 (2023)

Particle acceleration with Magnetic Reconnection in large scale RMHD simulations: I. Current sheet identification and characterization

Authors: Matteo Nurisso, Annalisa Celotti, Andrea Mignone, Gianluigi Bodo

Abstract: We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites in 2D and 3D large scale relativisticmagnetohydrodynamic numerical simulations. This has been implemented in the PLUTO code and tested in the cases of a single current sheet, a 2D jet and a 3D unstable plasma column. Its main features are: a) a computational cost which allows it… ▽ More We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites in 2D and 3D large scale relativisticmagnetohydrodynamic numerical simulations. This has been implemented in the PLUTO code and tested in the cases of a single current sheet, a 2D jet and a 3D unstable plasma column. Its main features are: a) a computational cost which allows its use in large scale simulations; b) the capability to deal with complex 2D and 3D structures of the reconnection sites. In the performed simulations, we identify the computational cells that are part of a current sheet by a measure of the gradient of the magnetic field along different directions. Lagrangian particles, which follow the fluid, are used to sample plasma parameters before entering the reconnection sites that form during the evolution of the different configurations considered. Specifically, we track the distributions of the magnetization parameter $σ$ and the thermal to magnetic pressure ratio $β$ that - according to particle-in-cell simulation results - control the properties of particle acceleration in magnetic reconnection regions. Despite the initial conditions of the simulations were not chosen "ad hoc", the 3D simulation returns results suitable for efficient particle acceleration and realistic non-thermal particle distributions. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 10 figures; MNRAS, in press

A guiding center implementation for relativistic particle dynamics in the PLUTO code

Authors: A. Mignone, H. Haudemand, E. Puzzoni

Abstract: We present a numerical implementation of the guiding center approximation to describe the relativistic motion of charged test particles in the PLUTO code for astrophysical plasma dynamics. The guiding center approximation (GCA) removes the time step constraint due to particle gyration around magnetic field lines by following the particle center of motion rather than its full trajectory. The gyrati… ▽ More We present a numerical implementation of the guiding center approximation to describe the relativistic motion of charged test particles in the PLUTO code for astrophysical plasma dynamics. The guiding center approximation (GCA) removes the time step constraint due to particle gyration around magnetic field lines by following the particle center of motion rather than its full trajectory. The gyration can be detached from the guiding center motion if electromagnetic fields vary sufficiently slow compared to the particle gyration radius and period. Our implementation employs a variable step-size linear multistep method, more efficient when compared to traditional one-step Runge Kutta schemes. A number of numerical benchmarks is presented in order to assess the validity of our implementation. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Jet-environment interplay in magnetized binary neutron star mergers

Authors: Andrea Pavan, Riccardo Ciolfi, Jay Vijay Kalinani, Andrea Mignone

Abstract: GRB 170817A, the first short gamma-ray burst (sGRB) to be detected in coincidence with a gravitational wave signal, demonstrated that merging binary neutron star (BNS) systems can power collimated ultra-relativistic jets and, in turn, produce sGRBs. Moreover, it revealed that sGRB jets possess an intrinsic angular structure that is imprinted in the observable prompt and afterglow emission. Advance… ▽ More GRB 170817A, the first short gamma-ray burst (sGRB) to be detected in coincidence with a gravitational wave signal, demonstrated that merging binary neutron star (BNS) systems can power collimated ultra-relativistic jets and, in turn, produce sGRBs. Moreover, it revealed that sGRB jets possess an intrinsic angular structure that is imprinted in the observable prompt and afterglow emission. Advanced numerical simulations represent the leading approach to investigate the physical processes underlying the evolution of sGRB jets breaking out of post-merger environments, and thus connect the final angular structure and energetics with specific jet launching conditions. In a previous paper, we carried out the first three-dimensional (3D) special-relativistic hydrodynamic simulations of incipient (top-hat) sGRB jets propagating across the realistic environment resulting from a general-relativistic (GR) hydrodynamic BNS merger simulation. While the above work marked an important step toward a consistent end-to-end description of sGRB jets from BNS mergers, those simulations did not account for the presence of magnetic fields, which are expected to play a key role. Here, we overcome this limitation, reporting the first 3D special-relativistic magnetohydrodynamic (MHD) simulation of a magnetized (structured and rotating) sGRB jet piercing through a realistic magnetized post-merger environment, wherein the initial conditions of the latter are directly imported from the outcome of a previous GRMHD BNS merger simulation. △ Less

Submitted 14 June, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

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

The impact of resistive electric fields on particle acceleration in reconnection layers

Authors: E. Puzzoni, A. Mignone, G. Bodo

Abstract: In the context of particle acceleration in high-energy astrophysical environments featuring magnetic reconnection, the importance of the resistive term of the electric field compared to the convective one is still under debate. In this work, we present a quantitative analysis through 2D magnetohydrodynamic numerical simulations of tearing-unstable current sheets coupled to a test-particles approac… ▽ More In the context of particle acceleration in high-energy astrophysical environments featuring magnetic reconnection, the importance of the resistive term of the electric field compared to the convective one is still under debate. In this work, we present a quantitative analysis through 2D magnetohydrodynamic numerical simulations of tearing-unstable current sheets coupled to a test-particles approach, performed with the PLUTO code. We find that the resistive field plays a significant role in the early-stage energization of high-energy particles. Indeed, these particles are firstly accelerated due to the resistive electric field when they cross an X-point, created during the fragmentation of the current sheet. If this preliminary particle acceleration mechanism dominated by the resistive field is neglected, particles cannot reach the same high energies. Our results support therefore the conclusion that the resistive field is not only non-negligible but it does actually play an important role in the particle acceleration mechanism. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: 7 pages, 5 figures

A numerical study of the interplay between Fermi acceleration mechanisms in radio lobes of FR-II radio galaxies

Authors: Sayan Kundu, Bhargav Vaidya, Andrea Mignone, Martin J. Hardcastle

Abstract: Context: Radio-loud AGNs are thought to possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Stochastic turbulent acceleration (STA) and diffusive shock acceleration (DSA) are commonly cited as potential sources of high-energy particles in weakly magnetized environments. Together, these acceleration processes and various radiative losses determine t… ▽ More Context: Radio-loud AGNs are thought to possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Stochastic turbulent acceleration (STA) and diffusive shock acceleration (DSA) are commonly cited as potential sources of high-energy particles in weakly magnetized environments. Together, these acceleration processes and various radiative losses determine the emission characteristics of these extra-galactic radio sources. Aims: The purpose of this research is to investigate the dynamical interplay between the STA and DSA in the radio lobes of FR-II radio galaxies, as well as the manner in which these acceleration mechanisms, along with a variety of radiative losses, collectively shape the emission features seen in these extra-galactic sources. Methods: A phenomenologically motivated model of STA is considered and subsequently employed on a magneto-hydrodynamically simulated radio lobe through a novel hybrid Eulerian-Lagrangian framework. Results: STA gives rise to a curved particle spectrum that is morphologically different from the usual shock-accelerated spectrum. As a consequence of this structural difference in the underlying particle energy spectrum, various multi-wavelength features arise in the spectral energy distribution of the radio lobe. Additionally, we observe enhanced diffuse X-ray emission from radio lobes for cases where STA is taken into account in addition to DSA. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 15 pages, 10 figures, Accepted for publication in A&A. Comments are welcome

Journal ref: A&A 667, A138 (2022)

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

A comparison of approximate non-linear Riemann solvers for Relativistic MHD

Authors: Giancarlo Mattia, Andrea Mignone

Abstract: We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate mo… ▽ More We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate modes that are included. Five different approaches - namely the HLL, HLLC, HLLD, HLLEM and GFORCE schemes - are compared in terms of accuracy and robustness against one- and multi-dimensional standard numerical benchmarks. Our results demonstrate that - for weak or moderate magnetizations - the HLLD Riemann solver yields the most accurate results, followed by HLLC solver(s). The GFORCE approach provides a valid alternative to the HLL solver being less dissipative and equally robust for strongly magnetized environments. Finally, our tests show that the HLLEM Riemann solver is not cost-effective in improving the accuracy of the solution and reducing the numerical dissipation. △ Less

Submitted 23 November, 2021; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Journal ref: MNRAS 500,1 (2022) 481-499

Acceleration of Solar Energetic Particles by the shock of Interplanetary Coronal Mass Ejection

Authors: Shanwlee Sow Mondal, Aveek Sarkar, Bhargav Vaidya, Andrea Mignone

Abstract: Interplanetary Coronal Mass Ejection (ICME) shocks are known to accelerate particles and contribute significantly to Solar Energetic Particle (SEP) events. We have performed Magnetohydrodynamic-Particle in Cell (MHD-PIC) simulations of ICME shocks to understand the acceleration mechanism. These shocks vary in Alfvénic Mach numbers as well as in magnetic field orientations (parallel \& quasi-perpen… ▽ More Interplanetary Coronal Mass Ejection (ICME) shocks are known to accelerate particles and contribute significantly to Solar Energetic Particle (SEP) events. We have performed Magnetohydrodynamic-Particle in Cell (MHD-PIC) simulations of ICME shocks to understand the acceleration mechanism. These shocks vary in Alfvénic Mach numbers as well as in magnetic field orientations (parallel \& quasi-perpendicular). We find that Diffusive Shock Acceleration (DSA) plays a significant role in accelerating particles in a parallel ICME shock. In contrast, Shock Drift Acceleration (SDA) plays a pivotal role in a quasi-perpendicular shock. High-Mach shocks are seen to accelerate particles more efficiently. Our simulations suggest that background turbulence and local particle velocity distribution around the shock can indirectly hint at the acceleration mechanism. Our results also point towards a few possible \textit{in situ} observations that could validate our understanding of the topic. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the Astrophysical Journal, 16 pages, 20 figures

On the Impact of the Numerical Method on Magnetic Reconnection and Particle Acceleration -- I. The MHD case

Authors: Eleonora Puzzoni, Andrea Mignone, Gianluigi Bodo

Abstract: We present 2D MHD numerical simulations of tearing-unstable current sheets coupled to a population of non-thermal test-particles, in order to address the problem of numerical convergence with respect to grid resolution, numerical method and physical resistivity. Numerical simulations are performed with the PLUTO code for astrophysical fluid dynamics through different combinations of Riemann solver… ▽ More We present 2D MHD numerical simulations of tearing-unstable current sheets coupled to a population of non-thermal test-particles, in order to address the problem of numerical convergence with respect to grid resolution, numerical method and physical resistivity. Numerical simulations are performed with the PLUTO code for astrophysical fluid dynamics through different combinations of Riemann solvers, reconstruction methods, grid resolutions at various Lundquist numbers. The constrained transport method is employed to control the divergence-free condition of magnetic field. Our results indicate that the reconnection rate of the background tearing-unstable plasma converges only for finite values of the Lundquist number and for sufficiently large grid resolutions. In general, it is found that (for a 2nd-order scheme) the minimum threshold for numerical convergence during the linear phases requires the number of computational zones covering the initial current sheet width to scale roughly as $\sim \sqrt{\bar{S}}$, where $\bar{S}$ is the Lundquist number defined on the current sheet width. On the other hand, the process of particle acceleration is found to be nearly independent of the underlying numerical details inasmuch as the system becomes tearing-unstable and enters in its nonlinear stages. In the limit of large $\bar{S}$, the ensuing power-law index quickly converge to $p \approx 1.7$, consistently with the fast reconnection regime. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 14 pages, 11 figures

Journal ref: MNRAS 508 (2021) 2771-2783

Morphology of radio relics II: Properties of polarised emission

Authors: P. Domínguez-Fernández, M. Brüggen, F. Vazza, M. Hoeft, W. E. Banda-Barragán, K. Rajpurohit, D. Wittor, A. Mignone, D. Mukherjee, B. Vaidya

Abstract: Radio relics are diffuse radio sources in galaxy clusters that are associated with merger shock waves. Detailed observations of radio relics in total intensity and in polarisation show complex structures on kiloparsec scales. The relation between the observed features and the underlying morphology of the magnetic field is not clear. Using three-dimensional magneto-hydrodynamical simulations, we st… ▽ More Radio relics are diffuse radio sources in galaxy clusters that are associated with merger shock waves. Detailed observations of radio relics in total intensity and in polarisation show complex structures on kiloparsec scales. The relation between the observed features and the underlying morphology of the magnetic field is not clear. Using three-dimensional magneto-hydrodynamical simulations, we study the polarised emission produced by a shock wave that propagates through a turbulent medium that resembles the intracluster medium. We model the polarised synchrotron emission on the basis of diffusive shock acceleration of cosmic-ray electrons. We find that the synchrotron emission produced in a shocked turbulent medium can reproduce some of the observed features in radio relics. Shock compression can give rise to a high polarisation fraction at the shock front and a partial alignment of the polarisation $E$-vectors with the shock normal. Our work confirms that radio relics can also be formed in an environment with a tangled magnetic field. We also discuss the effect of Faraday Rotation intrinsic to the source, and how our results depend on the angular resolution of observations. △ Less

Submitted 13 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS, 22 pages, 23 figures

Numerical modeling and physical interplay of stochastic turbulent acceleration for non-thermal emission processes

Authors: Sayan Kundu, Bhargav Vaidya, Andrea Mignone

Abstract: Particle acceleration is an ubiquitous phenomenon in astrophysical and space plasma. Diffusive shock acceleration (DSA) and stochastic turbulent acceleration are known to be the possible mechanisms for producing very high energetic particles, particularly in weakly magnetized regions. An interplay of different acceleration processes along with various radiation losses is typically observed in astr… ▽ More Particle acceleration is an ubiquitous phenomenon in astrophysical and space plasma. Diffusive shock acceleration (DSA) and stochastic turbulent acceleration are known to be the possible mechanisms for producing very high energetic particles, particularly in weakly magnetized regions. An interplay of different acceleration processes along with various radiation losses is typically observed in astrophysical sources. While DSA is a systematic acceleration process that energizes particles in the vicinity of shocks, stochastic turbulent acceleration (STA) is a random energizing process, where the interaction between cosmic ray particles and electromagnetic fluctuations results in particle acceleration. This process is usually interpreted as a biased random walk in energy space, modelled through a Fokker-Planck equation. In the present work, we describe a novel Eulerian algorithm, adopted to incorporate turbulent acceleration in the presence of DSA and radiative processes like synchrotron and Inverse-Compton emission. The developed framework extends the hybrid Eulerian-Lagrangian module in a full-fledged relativistic Magneto-hydrodynamic (RMHD) code PLUTO. From our validation tests and case studies, we showcase the competing and complementary nature of both acceleration processes. Axisymmetric simulations of an RMHD jet with this extended hybrid framework clearly demonstrate that emission due to shocks is localized while that due to turbulent acceleration originates in the backflow and is more diffuse, particularly in the high energy X-ray band. △ Less

Submitted 6 August, 2021; originally announced August 2021.

Comments: 16 Pages, 15 Figures, Accepted for publication in The Astrophysical Journal

3D MHD astrospheres: applications to IRC-10414 and Betelgeuse

Authors: D. M. -A. Meyer, A. Mignone, M. Petrov, K. Scherer, P. F. Velazquez, P. Boumis

Abstract: A significative fraction of all massive stars in the Milky Way move supersonically through their local interstellar medium (ISM), producing bow shock nebulae by wind-ISM interaction. The stability of these observed astrospheres around cool massive stars challenges precedent two-dimensional (magneto-)hydrodynamical simulations of their surroundings. We present three-dimensional magneto-hydrodynamic… ▽ More A significative fraction of all massive stars in the Milky Way move supersonically through their local interstellar medium (ISM), producing bow shock nebulae by wind-ISM interaction. The stability of these observed astrospheres around cool massive stars challenges precedent two-dimensional (magneto-)hydrodynamical simulations of their surroundings. We present three-dimensional magneto-hydrodynamical (3D MHD) simulations of the circumstellar medium of runaway M-type red supergiant stars moving with velocity v_star= 50 km/s. We treat the stellar wind with a Parker spiral and assume a 7 microG magnetisation of the ISM. Our free parameter is the angle theta_mag between ISM flow and magnetisation, taken to 0, 45 and 90 degrees. It is found that simulation dimension, coordinate systems and grid effects can greatly affect the development of the modelled astrospheres. Nevertheless, as soon as the ISM flow and magnetisation directions differs by more than a few degrees (theta_mag>5 degree), the bow shock is stabilised, most clumpiness and ragged structures vanishing. The complex shape of the bowshocks induce important projection effects, e.g. at optical Ha line, producing complex of astrospheric morphologies. We speculate that those effects are also at work around earlier-type massive stars, which would explain their diversity of their observed arc-like nebula around runaway OB stars. Our 3D MHD models are fitting well observations of the astrospheres of several runaway red supergiant stars. The results interpret the smoothed astrosphere of IRC-10414 and Betelgeuse aOri) are stabilised by an organised, non-parallel ambient magnetic field. Our findings suggest that IRC-10414 is currently in a steady state of its evolution, and that Betelgeuse's bar is of interstellar origin. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: Accepted to MNRAS

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

Short gamma-ray burst jet propagation in binary neutron star merger environments

Authors: Andrea Pavan, Riccardo Ciolfi, Jay Vijay Kalinani, Andrea Mignone

Abstract: The multimessenger event GW170817/GRB 170817A confirmed that binary neutron star (BNS) mergers can produce short gamma-ray burst (SGRB) jets. This evidence promoted new investigations on the mechanisms through which a BNS merger remnant can launch such a powerful relativistic outflow and on the propagation of the latter across the surrounding post-merger environment. In particular, great strides h… ▽ More The multimessenger event GW170817/GRB 170817A confirmed that binary neutron star (BNS) mergers can produce short gamma-ray burst (SGRB) jets. This evidence promoted new investigations on the mechanisms through which a BNS merger remnant can launch such a powerful relativistic outflow and on the propagation of the latter across the surrounding post-merger environment. In particular, great strides have been made in jet propagation models, establishing connections between the initial jet launching conditions, including the incipient jet launching time (with respect to merger) and the injection parameters, and the observable SGRB prompt and afterglow emission. However, present semi-analytical models and numerical simulations (with one notable exception) adopt simple hand-made prescriptions to account for the post-merger environment, lacking a direct association with any specific merging BNS system. Here, we present the first three-dimensional relativistic hydrodynamics simulations of incipient SGRB jets propagating through a post-merger environment that is directly imported from the outcome of a previous general relativistic BNS merger simulation. Our results show that the evolution and final properties of the jet can be largely affected by the anisotropies and the deviations from axisymmetry and homologous expansion characterizing more realistic BNS merger environments. In addition, we find that the inclusion of the gravitational pull from the central compact object, often overlooked, can have a major impact. Finally, we consider different jet launching times referred to the same BNS merger model and discuss the consequences for the ultimate jet properties. △ Less

Submitted 25 June, 2021; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: 16 pages, 20 figures

Streaming instability in a global patch simulation of protoplanetary disks

Authors: Mario Flock, Andrea Mignone

Abstract: In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two comp… ▽ More In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two components there. We address this problem by means of numerical 2.5D simulations in order to study the gas and dust interaction in fully global stratified disks. To this purpose, we employ the recently developed dust grain module in the PLUTO code. Our model focuses on a typical T Tauri disk model, simulating a short patch of the disk at 10 au which includes grains of constant Stokes number of $St=0.01$ and $St=0.1$, corresponding to grains with sizes of 0.9 cm and 0.9 mm, respectively, for the given disk model. By injecting a constant pebble flux at the outer domain, the system reaches a quasi steady state of turbulence and dust concentrations driven by the streaming instability. For our given setup and using resolutions up to 2500 cells per scale height we resolve the streaming instability, leading to local dust clumping and concentrations. Our results show dust density values of around 10-100 times the gas density with a steady state pebble flux between $3.5 \times 10^{-4}$ and $2.5 \times 10^{-3} M_{\rm Earth}/\mathit{year}$ for the models with $\mathit{St}=0.01$ and $\mathit{St}=0.1$. The grain size and pebble flux for model $\mathit{St}=0.01$ compares well with dust evolution models of the first million years of disk evolution. For those grains the scatter opacity dominates the extinction coefficient at mm wavelengths. These types of global dust and gas simulations are a promising tool for studies of the gas and dust evolution at pressure bumps in protoplanetary disks. △ Less

Submitted 28 March, 2021; originally announced March 2021.

Comments: accepted at A&A

SWELTO -- Space WEather Laboratory in Turin Observatory

Authors: A. Bemporad, L. Abbo, D. Barghini, C. Benna, R. Biondo, D. Bonino, G. Capobianco, F. Carella, A. Cora, S. Fineschi, F. Frassati, D. Gardiol, S. Giordano, A. Liberatore, S. Mancuso, A. Mignone, S. Rasetti, F. Reale, A. Riva, F. Salvati, R. Susino, A. Volpicelli, L. Zangrilli

Abstract: SWELTO -- Space WEather Laboratory in Turin Observatory is a conceptual framework where new ideas for the analysis of space-based and ground-based data are developed and tested. The input data are (but not limited to) remote sensing observations (EUV images of the solar disk, Visible Light coronagraphic images, radio dynamic spectra, etc...), in situ plasma measurements (interplanetary plasma dens… ▽ More SWELTO -- Space WEather Laboratory in Turin Observatory is a conceptual framework where new ideas for the analysis of space-based and ground-based data are developed and tested. The input data are (but not limited to) remote sensing observations (EUV images of the solar disk, Visible Light coronagraphic images, radio dynamic spectra, etc...), in situ plasma measurements (interplanetary plasma density, velocity, magnetic field, etc...), as well as measurements acquired by local sensors and detectors (radio antenna, fluxgate magnetometer, full-sky cameras, located in OATo). The output products are automatic identification, tracking, and monitoring of solar stationary and dynamic features near the Sun (coronal holes, active regions, coronal mass ejections, etc...), and in the interplanetary medium (shocks, plasmoids, corotating interaction regions, etc...), as well as reconstructions of the interplanetary medium where solar disturbances may propagate from the Sun to the Earth and beyond. These are based both on empirical models and numerical MHD simulations. The aim of SWELTO is not only to test new data analysis methods for future application for Space Weather monitoring and prediction purposes, but also to procure, test and deploy new ground-based instrumentation to monitor the ionospheric and geomagnetic responses to solar activity. Moreover, people involved in SWELTO are active in outreach to disseminate the topics related with Space Weather to students and the general public. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Journal ref: Published in the INAF Technical Reports series, num. 40, 2020 (https://openaccess.inaf.it/handle/20.500.12386/27715)

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)

Morphology of radio relics I: What causes the substructure of synchrotron emission?

Authors: P. Domínguez-Fernández, M. Brüggen, F. Vazza, W. E. Banda-Barragán, K. Rajpurohit, A. Mignone, D. Mukherjee, B. Vaidya

Abstract: High-resolution radio observations of cluster radio relics often show complex spatial and spectral features. However, it is not clear what these features reveal about the underlying magnetic field properties. We performed three-dimensional magneto-hydrodynamical simulations of merger shock waves propagating through a magnetised, turbulent intracluster medium. Our model includes the diffusive shock… ▽ More High-resolution radio observations of cluster radio relics often show complex spatial and spectral features. However, it is not clear what these features reveal about the underlying magnetic field properties. We performed three-dimensional magneto-hydrodynamical simulations of merger shock waves propagating through a magnetised, turbulent intracluster medium. Our model includes the diffusive shock acceleration of cosmic-ray electrons, their spatial advection and energy losses at run-time. With this set-up we can investigate the relation between radio substructure and pre-shock plasma conditions in the host cluster. We find that upstream turbulence plays a major role in shaping the properties of radio relics produced downstream. Within the assumption of diffusive shock acceleration, we can reproduce the observed discrepancy between the X-ray derived Mach number of shocks, and the Mach number inferred from radio spectra. Our simulated spectral index maps and profiles across the radio relic also suggest that the standard deviation of the upstream magnetic field must be relatively small ($σ_B\leq 1 \, μ$G) in order to reproduce observations and therefore, radio relics can potentially constrain the distribution of magnetic fields in galaxy clusters outskirts. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS, 22 pages, 25 figures

Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics

Authors: Rolf Kuiper, Harold W. Yorke, Andrea Mignone

Abstract: Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a new… ▽ More Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a newly developed module for photoionization, SEDNA, coupled to the magnetohydrodynamics code PLUTO. These extensions are currently not publicly available; access can be granted on a case-by-case basis. We explain the theoretical background of the equations solved, elaborate on the numerical layout, and present a comprehensive test suite for radiation-ionization hydrodynamics. The grid-based radiation and ionization modules support static one-dimensional, two-dimensional, and three-dimensional grids in Cartesian, cylindrical, and spherical coordinates. Each module splits the radiation field into two components, one originating directly from a point source - solved using a ray-tracing scheme - and a diffuse component - solved with a three-dimensional flux-limited diffusion (FLD) solver. The FLD solver for the continuum radiation transport makes use of either the equilibrium one-temperature approach or the linearization two-temperature approach. The FLD solver for the photoionization module enables accounting for the temporal evolution of the radiation field from direct recombination of free electrons into hydrogen's ground state as an alternative to the on-the-spot approximation. A brief overview of completed and ongoing scientific studies is given to explicitly illustrate the multipurpose nature of the numerical framework presented. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal Supplement Series, 2020, Volume 250, Issue 1, id.13, 28 pp

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

A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks

Authors: Julio David Melon Fuksman, Hubert Klahr, Mario Flock, Andrea Mignone

Abstract: We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly… ▽ More We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly anisotropic radiation transport as can be expected in the vicinity of an accreting planet in a protoplanetary disk. To reduce the computational cost caused by the timescale disparity between radiation and matter fields, we integrate their evolution equations separately in an operator-split way, using substepping to evolve the radiation equations. We further increase the code's efficiency by adopting the reduced speed of light approximation (RSLA). Our integration scheme for the evolution equations of radiation fields relies on implicit-explicit schemes, in which radiation-matter interaction terms are integrated implicitly while fluxes are integrated via Godunov-type solvers. The module is suitable for general astrophysical computations in 1, 2, and 3 dimensions in Cartesian, spherical and cylindrical coordinates, and can be implemented on rotating frames. We demonstrate the algorithm performance on different numerical benchmarks, paying particular attention to the applicability of the RSLA for computations of physical processes in protoplanetary disks. We show 2D simulations of vertical convection in disks and 3D simulations of gas accretion by planetary cores, which are the first of their kind to be solved with a two-moment approach. △ Less

Submitted 18 January, 2021; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 27 pages, 16 figures

Journal ref: The Astrophysical Journal, 906:78, 2021

Systematic construction of upwind constrained transport schemes for MHD

Authors: Andrea Mignone, Luca Del Zanna

Abstract: The constrained transport (CT) method reflects the state of the art numerical technique for preserving the divergence-free condition of magnetic field to machine accuracy in multi-dimensional MHD simulations performed with Godunov-type, or upwind, conservative codes. The evolution of the different magnetic field components, located at zone interfaces using a staggered representation, is achieved b… ▽ More The constrained transport (CT) method reflects the state of the art numerical technique for preserving the divergence-free condition of magnetic field to machine accuracy in multi-dimensional MHD simulations performed with Godunov-type, or upwind, conservative codes. The evolution of the different magnetic field components, located at zone interfaces using a staggered representation, is achieved by calculating the electric field components at cell edges, in a way that has to be consistent with the Riemann solver used for the update of cell-centered fluid quantities at interfaces. Albeit several approaches have been undertaken, the purpose of this work is, on the one hand, to compare existing methods in terms of robustness and accuracy and, on the other, to extend the \emph{upwind contrained transport} (UCT) method by Londrillo \& Del Zanna (2004) and Del Zanna et al. (2007) for the systematic construction of new averaging schemes. In particular, we propose a general formula for the upwind fluxes of the induction equation which simply involves the information available from the base Riemann solver employed for the fluid part, provided it does not require full spectral decomposition, and 1D reconstructions of velocity and magnetic field components from nearby intercell faces to cell edges. Our results are presented here in the context of second-order schemes for classical MHD, but they can be easily generalized to higher than second order schemes, either based on finite volumes or finite differences, and to other physical systems retaining the same structure of the equations, such as that of relativistic or general relativistic MHD. △ Less

Submitted 31 July, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 29 pages, 16 figures

Laboratory disruption of scaled astrophysical outflows by a misaligned magnetic field

Authors: G. Revet, B. Khiar, E. Filippov, C. Argiroffi, J. Béard, R. Bonito, M. Cerchez, S. N. Chen, T. Gangolf, D. P. Higginson, A. Mignone, B. Olmi, M. Ouillé, S. N. Ryazantsev, I. Yu. Skobelev, M. I. Safronova, M. Starodubtsev, T. Vinci, O. Willi, S. Pikuz, S. Orlando, A. Ciardi, J. Fuchs

Abstract: The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growin… ▽ More The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growing misalignments, this nozzle becomes increasingly asymmetric, disrupting jet formation. Our results thus suggest outflow/magnetic field misalignment to be a plausible key process regulating jet collimation in a variety of objects from our Sun's outflows to extragalatic jets. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field. △ Less

Submitted 20 December, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

A Particle Module for the PLUTO Code: III -- Dust

Authors: A. Mignone, M. Flock, B. Vaidya

Abstract: The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust p… ▽ More The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust particles which overcomes the stiffness introduced by the linear coupling term. Besides being time-reversible and globally second-order accurate in time, the exponential integrator provides energy errors which are always bounded and it remains stable in the limit of arbitrarily small particle stopping times yielding the correct asymptotic solution. Such properties make this method preferable to the more widely used semi-implicit or fully implicit schemes at a very modest increase in computational cost. Coupling between particles and grid quantities is achieved through particle deposition and field-weighting techniques borrowed from Particle-In-Cell simulation methods. In this respect, we derive new weight factors in curvilinear coordinates that are more accurate than traditional volume- or area-weighting. A comprehensive suite of numerical benchmarks is presented to assess the accuracy and robustness of the algorithm in Cartesian, cylindrical and spherical coordinates. Particular attention is devoted to the streaming instability which is analyzed in both local and global disk models. The module is part of the PLUTO code for astrophysical gas-dynamics and it is mainly intended for the numerical modeling of protoplanetary disks in which solid and gas interact via aerodynamic drag. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Comments: 22 pages, 13 figures

A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks

Authors: Siddhartha Gupta, Prateek Sharma, Andrea Mignone

Abstract: Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form.… ▽ More Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modeling, and future directions are discussed. △ Less

Submitted 1 February, 2021; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: 18 pages, 12 figures, 2 tables; MNRAS accepted version

A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations

Authors: A. Mignone, G. Mattia, G. Bodo, L. Del Zanna

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

Submitted 2 April, 2019; originally announced April 2019.

Comments: 25 pages, 14 figures

A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code

Authors: Julio David Melon Fuksman, Andrea Mignone

Abstract: We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer equations are solved under the grey approximation and imposing the M1 closure, which allows the radiation transport to be handled in both the free-streaming and diffusi… ▽ More We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer equations are solved under the grey approximation and imposing the M1 closure, which allows the radiation transport to be handled in both the free-streaming and diffusion limits. Equations are integrated following an implicit-explicit scheme, where radiation-matter interaction terms are integrated implicitly, whereas transport and all of the remaining source terms are solved explicitly by means of the same Godunov-type solvers included in PLUTO. Among these, we introduce a new Harten-van Leer-contact (HLLC) solver for optically thin radiation transport. The code is suitable for multidimensional computations in Cartesian, spherical and cylindrical coordinates, using either a single processor or parallel architectures. Adaptive grid computations are also made possible, by means of the CHOMBO library. The algorithm performance is demonstrated through a series of numerical benchmarks by investigating various different configurations with a particular emphasis on the behavior of the solutions in the free-streaming and diffusion limits. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: 25 pages, 17 figures

Journal ref: ApJS 242, 20 (2019)

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

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)

Linear Wave Propagation for Resistive Relativistic Magnetohydrodynamics

Authors: A. Mignone, G. Mattia, G. Bodo

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

Submitted 4 September, 2018; originally announced September 2018.

Comments: 20 pages, 13 figures

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

A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations

Authors: A. Mignone, G. Bodo, B. Vaidya, G. Mattia

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

Submitted 5 April, 2018; originally announced April 2018.

Comments: 27 pages, 16 figures. Accepted for publication in ApJ Supplement series

HESS J0632+057: hydrodynamics and nonthermal emission

Authors: Valentí Bosch-Ramon, Maxim V. Barkov, Andrea Mignone, Pol Bordas

Abstract: HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the… ▽ More HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the only gamma-ray emitting high-mass binary known to host a non-accreting pulsar. In this Letter, we present hydrodynamical simulations of HESS J0632+057 in the context of a pulsar and a stellar wind interacting in an eccentric binary, and propose a scenario for the non-thermal phenomenology of the source. In this scenario, the non-thermal activity before and around apastron is linked to the accumulation of non-thermal particles in the vicinity of the binary, and the sudden drop of the emission before apastron is produced by the disruption of the two-wind interaction structure, allowing these particles to efficiently escape. In addition to providing a framework to explain the non-thermal phenomenology of the source, this scenario predicts extended, moving X-ray emitting structures similar to those observed in PSR B1259-63. △ Less

Submitted 31 July, 2017; originally announced August 2017.

Comments: 6 pages, 5 figures, accepted by MNRAS Letters

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

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

Authors: Bhargav Vaidya, Deovrat Prasad, Andrea Mignone, Prateek Sharma, Luca Rickler

Abstract: An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat… ▽ More An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat parabolic terms arising from the diffusion operator. An explicit treatment of parabolic terms will considerably reduce the simulation time step due to its dependence on the square of the grid resolution ($Δx$) for stability. Although an implicit scheme relaxes the constraint on stability, it is difficult to distribute efficiently on a parallel architecture. Treating parabolic terms with accelerated super-time stepping (STS) methods has been discussed in literature but these methods suffer from poor accuracy (first order in time) and also have difficult-to-choose tuneable stability parameters. In this work we highlight a second order (in time) Runge Kutta Legendre (RKL) scheme (first described by Meyer et. al. 2012) that is robust, fast and accurate in treating parabolic terms alongside the hyperbolic conversation laws. We demonstrate its superiority over the first order super time stepping schemes with standard tests and astrophysical applications. We also show that explicit conduction is particularly robust in handling saturated thermal conduction. Parallel scaling of explicit conduction using RKL scheme is demonstrated up to more than $10^4$ processors. △ Less

Submitted 22 August, 2017; v1 submitted 17 February, 2017; originally announced February 2017.

Comments: 15 pages, 9 figures, incorporated comments from the referee. This version is now accepted for publication in MNRAS

Numerical Simulations of Torsional Alfvén Waves in Axisymmetric Solar Magnetic Flux Tubes

Authors: D. Wójcik, K. Murawski, Z. E. Musielak, P. Konkol, A. Mignone

Abstract: We investigate numerically Alfvén waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azim… ▽ More We investigate numerically Alfvén waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azimuthal direction, behave differently as a result of gradients of physical parameters along the tube. To explain these differences in the wave behavior, the time evolution of the wave variables and the resulting cutoff period for each wave variable are calculated, and used to determine regions in the solar chromosphere where strong wave reflection may occur. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: Submitted to Solar Physics (accepted)

Multi-D magnetohydrodynamic modelling of pulsar wind nebulae: recent progress and open questions

Authors: B. Olmi, L. Del Zanna, E. Amato, N. Bucciantini, A. Mignone

Abstract: In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting pa… ▽ More In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting particles and the quest for the acceleration sites of particles of different energies along the termination shock, by using wisps motion as a diagnostic tool; the study of the magnetic dissipation process in high magnetization nebulae by means of new long-term three-dimensional simulations of the pulsar wind nebula evolution; the investigation of the relativistic tearing instability in thinning current sheets, leading to fast reconnection events that might be at the origin of the Crab nebula gamma-ray flares. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 30 pages, 12 figures

Bow shock nebulae of hot massive stars in a magnetized medium

Authors: D. M. -A. Meyer, A. Mignone, R. Kuiper, A. Raga, W. Kley

Abstract: A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic field… ▽ More A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic fields, aligned with the motion of the flow, on the formation and emission properties of these circumstellar structures. Our axisymmetric, magneto-hydrodynamical simulations with optically-thin radiative cooling, heating and anisotropic thermal conduction show that the presence of the background ISM magnetic field affects the projected optical emission our bow shocks at Ha and [OIII] lambda 5007 which become fainter by about 1-2 orders of magnitude, respectively. Radiative transfer calculations against dust opacity indicate that the magnetic field slightly diminishes their projected infrared emission and that our bow shocks emit brightly at 60 micron. This may explain why the bow shocks generated by ionizing runaway massive stars are often difficult to identify. Finally, we discuss our results in the context of the bow shock of Zeta Ophiuchi and we support the interpretation of its imperfect morphology as an evidence of the presence of an ISM magnetic field not aligned with the motion of its driving star. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: Accepted by MNRAS, 21 pages, 13 figures

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

MHD simulations of three-dimensional Resistive Reconnection in a cylindrical plasma column

Authors: Edoardo Striani, Andrea Mignone, Bhargav Vaidya, Gianluigi Bodo, Attilio Ferrari

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

Submitted 25 July, 2016; originally announced July 2016.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

3D MHD modeling of twisted coronal loops

Authors: F. Reale, S. Orlando, M. Guarrasi, A. Mignone, G. Peres, A. W. Hood, E. R. Priest

Abstract: We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tu… ▽ More We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 29 pages, 14 figures, accepted for publication, movies available upon request to the first author

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