Showing 1–50 of 70 results for author: Del Zanna, L

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

A GPU-Accelerated Modern Fortran Version of the ECHO Code for Relativistic Magnetohydrodynamics

Authors: Luca Del Zanna, Simone Landi, Lorenzo Serafini, Matteo Bugli, Emanuele Papini

Abstract: The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme magnetization, curved spacetimes near compact objects) and the large variety of spatial scales needed to resolve turbulent motions. A great benefit comes from the porti… ▽ More The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme magnetization, curved spacetimes near compact objects) and the large variety of spatial scales needed to resolve turbulent motions. A great benefit comes from the porting of existing codes running on standard processors to GPU-based platforms. However, this usually requires a drastic rewriting of the original code, the use of specific languages like CUDA, and a complex analysis of data management and optimization of parallel processes. Here we describe the porting of the ECHO code for special and general relativistic MHD to accelerated devices, simply based on native Fortran language built-in constructs, especially 'do concurrent' loops, few OpenACC directives, and the straightforward data management provided by the Unified Memory option of NVIDIA compilers.Thanks to these very minor modifications to the original code, the new version of ECHO runs at least 16 times faster on GPU platforms compared to CPU-based ones. The chosen benchmark is the 3D propagation of a relativistic MHD Alfvén wave, for which strong and weak scaling tests performed on the LEONARDO pre-exascale supercomputer at CINECA are provided (using up to 256 nodes corresponding to 1024 GPUs, and over 14 billion cells). Finally, an example of high-resolution relativistic MHD Alfvénic turbulence simulation is shown, demonstrating the potential for astrophysical plasmas of the new GPU-based version of ECHO. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: Accepted for publication on Fluids, MDPI, 17 pages

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)

The essential role of multi-point measurements in investigations of turbulence, three-dimensional structure, and dynamics: the solar wind beyond single scale and the Taylor Hypothesis

Authors: W. H. Matthaeus, S. Adhikari, R. Bandyopadhyay, M. R. Brown, R. Bruno, J. Borovsky, V. Carbone, D. Caprioli, A. Chasapis, R. Chhiber, S. Dasso, P. Dmitruk, L. Del Zanna, P. A. Dmitruk, Luca Franci, S. P. Gary, M. L. Goldstein, D. Gomez, A. Greco, T. S. Horbury, Hantao Ji, J. C. Kasper, K. G. Klein, S. Landi, Hui Li , et al. (27 additional authors not shown)

Abstract: Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer h… ▽ More Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere. △ Less

Submitted 26 November, 2022; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

Quasi-universality of the magnetic deformation for neutron stars in general relativity and beyond

Authors: J. Soldateschi, N. Bucciantini, L. Del Zanna

Abstract: Neutron stars harbour extremely powerful magnetic fields, leading to their shape being deformed. Their magnetic deformation depends both on the geometry - and strength - of their internal magnetic field and on their composition, encoded by the equation of state. However, both the details of the internal magnetic structure and the equation of state of the innermost part of neutron stars are mostly… ▽ More Neutron stars harbour extremely powerful magnetic fields, leading to their shape being deformed. Their magnetic deformation depends both on the geometry - and strength - of their internal magnetic field and on their composition, encoded by the equation of state. However, both the details of the internal magnetic structure and the equation of state of the innermost part of neutron stars are mostly unkown. We performed a study of numerical models of magnetised, static, axisymmetric neutron stars in general relativity and in one of its most promising extensions, scalar-tensor theories. We did so by using several realistic equations of state currently allowed by observational and nuclear physics constraints, considering also those for strange quark stars. We show that it is possible to find simple relations among the magnetic deformation of a neutron star, its Komar mass, and its circumferential radius in the case of purely poloidal and purely toroidal magnetic configurations satisfying the equilibrium criterion in the Bernoulli formalism. These relations are quasi-universal, in the sense that they mostly do not depend on the equation of state. Our results, being formulated in terms of potentially observable quantities, could help to understand the magnetic properties of neutron stars interiors and the detectability of continuous gravitational waves by isolated neutron stars, independently of their equation of state. In the case of scalar-tensor theories, these relations depend also on the scalar charge of the neutron stars, thus potentially providing a new way to set constraints on the theory of gravitation. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 19 pages, 5 figures, submitted as proceeding to the Sixteenth Marcel Grossmann Meeting, this article draws heavily from arXiv:2106.00603

Are GRMHD Mean-Field Dynamo Models of Thick Accretion Disks SANE?

Authors: Niccolò Tomei, Luca Del Zanna, Matteo Bugli, Niccolò Bucciantini

Abstract: The remarkable results by the Event Horizon Telescope collaboration concerning the emission from M87* and, more recently, its polarization properties, require an increasingly accurate modeling of the plasma flows around the accreting black hole. Radiatively inefficient sources such as M87* and Sgr A* are typically modeled with the SANE (standard and normal evolution) paradigm, if the accretion dyn… ▽ More The remarkable results by the Event Horizon Telescope collaboration concerning the emission from M87* and, more recently, its polarization properties, require an increasingly accurate modeling of the plasma flows around the accreting black hole. Radiatively inefficient sources such as M87* and Sgr A* are typically modeled with the SANE (standard and normal evolution) paradigm, if the accretion dynamics is smooth, or with the MAD (magnetically arrested disk) paradigm, if the black hole's magnetosphere reacts by halting the accretion sporadically, resulting in a highly dynamical process. While the recent polarization studies seem to favor MAD models, this may not be true for all sources, and SANE accretion surely still deserves attention. In this work, we investigate the possibility of reaching the typical degree of magnetization and other accretion properties expected for SANE disks by resorting to the mean-field dynamo process in axisymmetric GRMHD simulations, which are supposed to mimic the amplifying action of an unresolved magnetorotational instability-driven turbulence. We show that it is possible to reproduce the main diagnostics present in the literature by starting from very unfavorable initial configurations, such as a purely toroidal magnetic field with negligible magnetization. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Quasi-universality of the magnetic deformation of neutron stars in general relativity and beyond

Authors: J. Soldateschi, N. Bucciantini, L. Del Zanna

Abstract: Neutron stars are known to host extremely powerful magnetic fields. Among other effects, one of the consequences of harbouring such fields is the deformation of the neutron star structure, leading, together with rotation, to the emission of continuous gravitational waves. On the one hand, the details of their internal magnetic fields are mostly unknown. Likewise, their internal structure, encoded… ▽ More Neutron stars are known to host extremely powerful magnetic fields. Among other effects, one of the consequences of harbouring such fields is the deformation of the neutron star structure, leading, together with rotation, to the emission of continuous gravitational waves. On the one hand, the details of their internal magnetic fields are mostly unknown. Likewise, their internal structure, encoded by the equation of state, is highly uncertain. Here we present a study of axisymmetric models of isolated magnetised neutron stars, for various realistic equations of state considered viable by observations and nuclear physics constraints. We show that it is possible to find simple relations between the magnetic deformation of a neutron star, its Komar mass and its circumferential radius. Such relations are quasi-universal, meaning that they are mostly independent on the equation of state of the neutron star and only slightly dependent on the magnetic field configuration. Being formulated in terms of potentially observable quantities, as we discuss, our results could help to constrain the magnetic properties of the neutron star interior and to better assess the detectability of continuous gravitational waves by isolated neutron stars, without knowing their equation of state. Our results are derived both in general relativity and in scalar-tensor theories - one of the most promising extensions of general relativity - in this case by considering also the scalar charge. We show that even in this case general relations hold that account for deviations from general relativity, that could potentially be used to set constraints on the gravitational theory. △ Less

Submitted 11 October, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A; minor language corrections

Journal ref: A&A 654, A162 (2021)

Axisymmetric equilibrium models of magnetised neutron stars in scalar-tensor theories

Authors: Jacopo Soldateschi, Niccolò Bucciantini, Luca Del Zanna

Abstract: General relativity probably is not the definitive theory of gravity, due a number or issues, both from the theoretical and from the observational point of view. Alternative theories of gravity were conceived to extend general relativity and account for such issues. Among the most promising ones are scalar-tensor theories, which predict an enrichment of the phenomenology of compact objects, like ne… ▽ More General relativity probably is not the definitive theory of gravity, due a number or issues, both from the theoretical and from the observational point of view. Alternative theories of gravity were conceived to extend general relativity and account for such issues. Among the most promising ones are scalar-tensor theories, which predict an enrichment of the phenomenology of compact objects, like neutron stars. We updated the well-tested XNS code to numerically solve the Einstein-Maxwell equations for a stationary, magnetised neutron star in a class of scalar-tensor theories containing the spontaneous scalarisation phenomenon. We found that there exist "quasi-universal relations" among the mass, radius, scalar charge and magnetic deformation of a neutron star that are true independently of the equation of state, both in general relativity and in scalar-tensor theories. This result could potentially provide new tools to test scalar-tensor theories and the magnetic field geometry inside neutron stars. △ Less

Submitted 12 May, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: Contribution to the 2021 Gravitation session of the 55th Rencontres de Moriond; corrected a minor mistake

Magnetic deformation of neutron stars in scalar-tensor theories

Authors: J. Soldateschi, N. Bucciantini, L. Del Zanna

Abstract: Scalar-tensor theories are among the most promising alternatives to general relativity that have been developed to account for some long standing issues in our understanding of gravity. Some of these theories predict the existence of a non-linear phenomenon, spontaneous scalarisation, which can lead to the appearance of sizeable modifications to general relativity in the presence of compact matter… ▽ More Scalar-tensor theories are among the most promising alternatives to general relativity that have been developed to account for some long standing issues in our understanding of gravity. Some of these theories predict the existence of a non-linear phenomenon, spontaneous scalarisation, which can lead to the appearance of sizeable modifications to general relativity in the presence of compact matter distributions, namely neutron stars. On the one hand, one of the effects of the scalar field is to modify the emission of gravitational waves, both due to variations in the quadrupolar deformation of the star and to the presence of additional modes of emission. On the other hand, neutron stars are known to harbour extremely powerful magnetic fields which can affect their structure and shape, leading in turn to the emission of gravitational waves, this time due to a magnetic quadrupolar deformation. We investigate here how the presence of spontaneous scalarisation can affect the magnetic deformation of neutron stars and their emission of quadrupolar gravitational waves, both of tensor and scalar nature. We will show that it is possible to provide simple parameterisations of the magnetic deformation and gravitational wave power of neutron stars in terms of their baryonic mass, circumferential radius and scalar charge, and that a universal scaling exists independently of the magnetic field geometry and of the parameters of the scalar-tensor theory. Finally, we comment on the observability of the deviations of the gravitational waves strain from general relativity by current and future observatories. △ Less

Submitted 9 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 9 pages, 2 figures, accepted for publication in A&A

Journal ref: A&A 645, A39 (2021)

General Relativistic Mean-field Dynamo Model for Proto-neutron Stars

Authors: K. Franceschetti, L. Del Zanna

Abstract: Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo d… ▽ More Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo due to sub-scale, unresolved turbulence. The dynamo process, combined with the differential rotation of the (proto-)star, is able to produce an exponential growth of any initial magnetic seed field up to the values required to explain the observations. By varying the dynamo coefficient we obtain different growth rates. We find a quasi-linear dependence of the growth rates on the intensity of the dynamo. Furthermore, the time interval in which exponential growth occurs and the growth rates also seems to depend on the initial configuration of the magnetic field. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 12 pages, 4 figures, accepted for publication on Universe

3D Relativistic MHD Simulations of Pulsar Bow Shock Nebulae

Authors: N. Bucciantini, B. Olmi, L. Del Zanna

Abstract: Pulsars out of their parent SNR directly interact with the ISM producing so called Bow-Shock Pulsar Wind Nebulae, the relativistic equivalents of the heliosphere/heliotail system. These have been directly observed from Radio to X-ray, and are found also associated to TeV halos, with a large variety of morphologies. They offer a unique environment where the pulsar wind can be studied by modelling i… ▽ More Pulsars out of their parent SNR directly interact with the ISM producing so called Bow-Shock Pulsar Wind Nebulae, the relativistic equivalents of the heliosphere/heliotail system. These have been directly observed from Radio to X-ray, and are found also associated to TeV halos, with a large variety of morphologies. They offer a unique environment where the pulsar wind can be studied by modelling its interaction with the surrounding ambient medium, in a fashion that is different/complementary from the canonical Plerions. These systems have also been suggested as the possible origin of the positron excess detected by AMS and PAMELA, in contrast to dark matter. I will present results from 3D Relativistic MHD simulations of such nebulae. On top of these simulations we computed the expected emission signatures, the properties of high energy particle escape, the role of current sheets in channeling cosmic rays, the level of turbulence and magnetic amplification, and how they depend on the wind structure and magnetisation. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 10 pages, 3 figures, Proceedings of ASTRONUM-2019, Accepted in Journal of Physics Conference Series

Axisymmetric equilibrium models for magnetised neutron stars in Scalar-Tensor Theories

Authors: J. Soldateschi, N. Bucciantini, L. Del Zanna

Abstract: Among the possible extensions of General Relativity that have been put forward in order to address some long standing issues in our understanding of the Universe, Scalar-Tensor Theories have received a lot of attention for their simplicity. Interestingly, some of these predict a potentially observable non-linear phenomenon, known as \textit{spontaneous scalarisation}, in the presence of highly com… ▽ More Among the possible extensions of General Relativity that have been put forward in order to address some long standing issues in our understanding of the Universe, Scalar-Tensor Theories have received a lot of attention for their simplicity. Interestingly, some of these predict a potentially observable non-linear phenomenon, known as \textit{spontaneous scalarisation}, in the presence of highly compact matter distributions, like the case of neutron stars. Neutron stars are ideal laboratories to investigate the properties of matter under extreme conditions, and in particular they are known to harbour the strongest magnetic fields in the Universe. Here, for the first time, we present a detailed study of magnetised neutron stars in Scalar-Tensor Theories. First, we show that the formalism developed for the study of magnetised neutron stars in General Relativity, based on the \textit{eXtended Conformally Flat Condition}, can easily be extended in the presence of a non-minimally coupled scalar field, retaining many of its numerical advantages. We then carry out a study of the parameter space considering the two extreme geometries of purely toroidal and purely poloidal magnetic fields, varying both the strength of the magnetic field and the intensity of scalarisation. We compare our results with magnetised general-relativistic solutions and un-magnetised scalarised solutions, showing how the mutual interplay between magnetic and scalar fields affect the magnetic and the scalarisation properties of neutron stars. In particular, we focus our discussion on magnetic deformability, maximum mass and range of scalarisation. △ Less

Submitted 16 July, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Comments: accepted for publication by A&A; minor language corrections; minor typos correction

Journal ref: A&A 640, A44 (2020)

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

General relativistic magnetohydrodynamic dynamo in thick accretion disks: fully nonlinear simulations

Authors: N. Tomei, L. Del Zanna, M. Bugli, N. Bucciantini

Abstract: The recent imaging of the M87 black hole at millimeter wavelengths by the Event Horizon Telescope (EHT) collaboration has triggered a renewed interest in numerical models for the accretion of magnetized plasma in the regime of general relativistic magnetohydrodynamics (GRMHD). Here non-ideal simulations, including both the resistive effects and, above all, the mean-field dynamo action due to sub-s… ▽ More The recent imaging of the M87 black hole at millimeter wavelengths by the Event Horizon Telescope (EHT) collaboration has triggered a renewed interest in numerical models for the accretion of magnetized plasma in the regime of general relativistic magnetohydrodynamics (GRMHD). Here non-ideal simulations, including both the resistive effects and, above all, the mean-field dynamo action due to sub-scale, unresolved turbulence, are applied for the first time to such systems in the fully nonlinear regime. Combined with the differential rotation of the disk, the dynamo process is able to produce an exponential growth of any initial seed magnetic field up to the values required to explain the observations, when the instability tends to saturate even in the absence of artificial quenching effects. Before reaching the final saturation stage we observe a secondary regime of exponential growing, where the magnetic field increases more slowly due to accretion, which is modifying the underlying equilibrium. By varying the dynamo coefficient we obtain different growth rates, though the field seems to saturate at approximately the same level, at least for the limited range of parameters explored here, providing substantial values for the MAD parameter for magnetized accretion. For reasonable values of the central mass density and the commonly employed recipes for synchrotron emission by relativistically hot electrons, our model is able to reproduce naturally the observed flux of Sgr A*, the next target for EHT. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Comments: 15 pages, 13 figures

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

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

[Plasma 2020 Decadal] The essential role of multi-point measurements in turbulence investigations: the solar wind beyond single scale and beyond the Taylor Hypothesis

Authors: W. H. Matthaeus, R. Bandyopadhyay, M. R. Brown, J. Borovsky, V. Carbone, D. Caprioli, A. Chasapis, R. Chhiber, S. Dasso, P. Dmitruk, L. Del Zanna, P. A. Dmitruk, Luca Franci, S. P. Gary, M. L. Goldstein, D. Gomez, A. Greco, T. S. Horbury, Hantao Ji, J. C. Kasper, K. G. Klein, S. Landi, Hui Li, F. Malara, B. A. Maruca , et al. (24 additional authors not shown)

Abstract: This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for su… ▽ More This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence. △ Less

Submitted 16 March, 2019; originally announced March 2019.

Comments: White paper submitted to the PLASMA 2020 Decadal Survey Committee

Covariant and 3+1 Equations for Dynamo-Chiral General Relativistic Magnetohydrodynamics

Authors: Luca Del Zanna, Niccolò Bucciatini

Abstract: The exponential amplification of initial seed magnetic fields in relativistic plasmas is a very important topic in astrophysics, from the conditions in the early Universe to the interior of neutron stars. While dynamo action in a turbulent plasma is often invoked, in the last years a novel mechanism of quantum origin has gained increasingly more attention, namely the Chiral Magnetic Effect (CME).… ▽ More The exponential amplification of initial seed magnetic fields in relativistic plasmas is a very important topic in astrophysics, from the conditions in the early Universe to the interior of neutron stars. While dynamo action in a turbulent plasma is often invoked, in the last years a novel mechanism of quantum origin has gained increasingly more attention, namely the Chiral Magnetic Effect (CME). This has been recognized in semi-metals and it is most likely at work in the quark-gluon plasma formed in heavy-ion collision experiments, where the highest magnetic fields in nature, up to B~10^18 G, are produced. This effect is expected to survive even at large hydrodynamical/MHD scales and it is based on the chiral anomaly due to an imbalance between left- and right-handed relativistic fermions in the constituent plasma. Such imbalance leads to an electric current parallel to an external magnetic field, which is precisely the same mechanism of an alpha-dynamo action in classical MHD. Here we extend the close parallelism between the chiral and the dynamo effects to relativistic plasmas and we propose a unified, fully covariant formulation of the generalized Ohm's law. Moreover, we derive for the first time the 3+1 general relativistic MHD equations for a chiral plasma both in flat and curved spacetimes, in view of numerical investigation of the CME in compact objects, especially magnetars, or of the interplay among the non-ideal magnetic effects of dynamo, the CME and reconnection. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: 11 pages, 3 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Fast magnetic reconnection: The "ideal" tearing instability in classic, Hall, and relativistic plasmas

Authors: Emanuele Papini, Simone Landi, Luca Del Zanna

Abstract: Magnetic reconnection is believed to be the driver of many explosive phenomena in Astrophysics, from solar to gamma-ray flares in magnetars and in the Crab nebula. However, reconnection rates from classic MHD models are far too slow to explain such observations. Recently, it was realized that when a current sheet gets sufficiently thin, the reconnection rate of the tearing instability becomes "ide… ▽ More Magnetic reconnection is believed to be the driver of many explosive phenomena in Astrophysics, from solar to gamma-ray flares in magnetars and in the Crab nebula. However, reconnection rates from classic MHD models are far too slow to explain such observations. Recently, it was realized that when a current sheet gets sufficiently thin, the reconnection rate of the tearing instability becomes "ideal", in the sense that the current sheet destabilizes on the "macroscopic" Alfvénic timescales, regardless of the Lundquist number of the plasma. Here we present 2D compressible MHD simulations in the classical, Hall, and relativistic regimes. In particular, the onset of secondary tearing instabilities is investigated within Hall-MHD for the first time. In the frame of relativistic MHD, we summarize the main results from Del Zanna et al. [1]: the relativistic tearing instability is found to be extremely fast, with reconnection rates of the order of the inverse of the light crossing time, as required to explain the high-energy explosive phenomena. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 12 pages, 6 figures, submitted to the Proceedings of the 12th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2017)

Relativistic MHD modeling of magnetized neutron stars, pulsar winds, and their nebulae

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

Abstract: Neutron stars are among the most fascinating astrophysical sources, being characterized by strong gravity, densities about the nuclear one or even above, and huge magnetic fields. Their observational signatures can be extremely diverse across the electromagnetic spectrum, ranging from the periodic and low-frequency signals of radio pulsars, up to the abrupt high-energy gamma-ray flares of magnetar… ▽ More Neutron stars are among the most fascinating astrophysical sources, being characterized by strong gravity, densities about the nuclear one or even above, and huge magnetic fields. Their observational signatures can be extremely diverse across the electromagnetic spectrum, ranging from the periodic and low-frequency signals of radio pulsars, up to the abrupt high-energy gamma-ray flares of magnetars, where energies of ~10^46 erg are released in a few seconds. Fast-rotating and highly magnetized neutron stars are expected to launch powerful relativistic winds, whose interaction with the supernova remnants gives rise to the non-thermal emission of pulsar wind nebulae, which are known cosmic accelerators of electrons and positrons up to PeV energies. In the extreme cases of proto-magnetars (magnetic fields of ~10^15 G and millisecond periods), a similar mechanism is likely to provide a viable engine for the still mysterious gamma-ray bursts. The key ingredient in all these spectacular manifestations of neutron stars is the presence of strong magnetic fields in their constituent plasma. Here we will present recent updates of a couple of state-of-the-art numerical investigations by the high-energy astrophysics group in Arcetri: a comprehensive modeling of the steady-state axisymmetric structure of rotating magnetized neutron stars in general relativity, and dynamical 3-D MHD simulations of relativistic pulsar winds and their associated nebulae. △ Less

Submitted 3 November, 2017; originally announced November 2017.

Comments: EPS 44th Conference on Plasma Physics (June 2017, Belfast), paper accepted for publication on Plasma Physics and Controlled Fusion

Journal ref: Plasma Phys. Control. Fusion 60 (2018) 014027

Sub-structure formation in starless cores

Authors: C. Toci, D. Galli, A. Verdini, L. Del Zanna, S. Landi

Abstract: Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as… ▽ More Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: 8 pages, 4 figures

Papaloizou-Pringle instability suppression by the magnetorotational instability in relativistic accretion discs

Authors: Matteo Bugli, Jérôme Guilet, Ewald Müller, Luca Del Zanna, Niccolò Bucciantini, Pedro J. Montero

Abstract: Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, how… ▽ More Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events. △ Less

Submitted 18 October, 2017; v1 submitted 6 July, 2017; originally announced July 2017.

Comments: 15 pages, 13 figures, submitted to MNRAS

General relativistic models for rotating magnetized neutron stars in conformally flat spacetime

Authors: A. G. Pili, N. Bucciantini, L. Del Zanna

Abstract: The extraordinary energetic activity of magnetars is usually explained in terms of dissipation of a huge internal magnetic field of the order of $10^{15-16}$G. How such a strong magnetic field can originate during the formation of a neutron star is still subject of active research. An important role can be played by fast rotation: if magnetars are born as millisecond rotators dynamo mechanisms may… ▽ More The extraordinary energetic activity of magnetars is usually explained in terms of dissipation of a huge internal magnetic field of the order of $10^{15-16}$G. How such a strong magnetic field can originate during the formation of a neutron star is still subject of active research. An important role can be played by fast rotation: if magnetars are born as millisecond rotators dynamo mechanisms may efficiently amplify the magnetic field inherited from the progenitor star during the collapse. In this case, the combination of rapid rotation and strong magnetic field determine the right physical condition not only for the development of a powerful jet driven explosion, manifesting as a gamma ray burst, but also for a copious gravitational waves emission. Strong magnetic fields are indeed able to induce substantial quadrupolar deformations in the star. In this paper we analyze the joint effect of rotation and magnetization on the structure of a polytropic and axisymmetric neutron star, within the ideal magneto-hydrodynamic regime. We will consider either purely toroidal or purely poloidal magnetic field geometries. Through the sampling of a large parameter space, we generalize previous results in literature, inferring new quantitative relations that allow for a parametrization of the induced deformation, that takes into account also the effects due to the stellar compactness and the current distribution. Finally, in the case of purely poloidal field, we also discuss how different prescriptions on the surface charge distribution (a gauge freedom) modify the properties of the surrounding electrosphere and its physical implications. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Comments: 25 pages, 17 figures, 6 tables, accepted for publication in MNRAS

Modeling the effect of small-scale magnetic turbulence on the X-ray properties of Pulsar Wind Nebulae

Authors: N. Bucciantini, R. Bandiera, B. Olmi, L. Del Zanna

Abstract: Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to test our current understanding of relativistic plasma processes. It is well known that magnetic fields play a crucial role in their dynamics and emission properties. At present, one of the main issues concerns the level of magnetic turbulence present in these systems, which in the absence of space resolved X-ray polarizati… ▽ More Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to test our current understanding of relativistic plasma processes. It is well known that magnetic fields play a crucial role in their dynamics and emission properties. At present, one of the main issues concerns the level of magnetic turbulence present in these systems, which in the absence of space resolved X-ray polarization measures cannot be directly constrained. In this work we investigate, for the first time using simulated synchrotron maps, the effect of a small scale fluctuating component of the magnetic field on the emission properties in X-ray. We illustrate how to include the effects of a turbulent component in standard emission models for PWNe, and which consequences are expected in terms of net emissivity and depolarization, showing that the X-ray surface brightness maps can provide already some rough constraints. We then apply our analysis to the Crab and Vela nebulae and, by comparing our model with Chandra and Vela data, we found that the typical energies in the turbulent component of the magnetic field are about 1.5 to 3 times the one in the ordered field. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS

Multidimensional relativistic MHD simulations of Pulsar Wind Nebulae: dynamics and emission

Authors: Luca Del Zanna, Barbara Olmi

Abstract: Pulsar Wind Nebulae, and the Crab nebula in particular, are the best cosmic laboratories to investigate the dynamics of magnetized relativistic outflows and particle acceleration up to PeV energies. Multidimensional MHD modeling by means of numerical simulations has been very successful at reproducing, to the very finest details, the innermost structure of these synchrotron emitting nebulae, as ob… ▽ More Pulsar Wind Nebulae, and the Crab nebula in particular, are the best cosmic laboratories to investigate the dynamics of magnetized relativistic outflows and particle acceleration up to PeV energies. Multidimensional MHD modeling by means of numerical simulations has been very successful at reproducing, to the very finest details, the innermost structure of these synchrotron emitting nebulae, as observed in the X-rays. Therefore, the comparison between the simulated source and observations can be used as a powerful diagnostic tool to probe the physical conditions in pulsar winds, like their composition, magnetization, and degree of anisotropy. However, in spite of the wealth of observations and of the accuracy of current MHD models, the precise mechanisms for magnetic field dissipation and for the acceleration of the non-thermal emitting particles are mysteries still puzzling theorists to date. Here we review the methodologies of the computational approach to the modeling of Pulsar Wind Nebulae, discussing the most relevant results and the recent progresses achieved in this fascinating field of high-energy astrophysics. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: 29 pages review, preliminary version. To appear in the book "Modelling Nebulae" edited by D. Torres for Springer, based on the invited contributions to the workshop held in Sant Cugat (Barcelona), June 14-17, 2016

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

Activation of MHD reconnection on ideal timescales

Authors: S. Landi, E. Papini, L. Del Zanna, A. Tenerani, F. Pucci

Abstract: Magnetic reconnection in laboratory, space and astrophysical plasmas is often invoked to explain explosive energy release and particle acceleration. However, the timescales involved in classical models within the macroscopic MHD regime are far too slow to match the observations. Here we revisit the tearing instability by performing visco-resistive two-dimensional numerical simulations of the evolu… ▽ More Magnetic reconnection in laboratory, space and astrophysical plasmas is often invoked to explain explosive energy release and particle acceleration. However, the timescales involved in classical models within the macroscopic MHD regime are far too slow to match the observations. Here we revisit the tearing instability by performing visco-resistive two-dimensional numerical simulations of the evolution of thin current sheets, for a variety of initial configurations and of values of the Lunquist number $S$, up to $10^7$. Results confirm that when the critical aspect ratio of $S^{1/3}$ is reached in the reconnecting current sheets, the instability proceeds on ideal (Alfvénic) macroscopic timescales, as required to explain observations. Moreover, the same scaling is seen to apply also to the local, secondary reconnection events triggered during the nonlinear phase of the tearing instability, thus accelerating the cascading process to increasingly smaller spatial and temporal scales. The process appears to be robust, as the predicted scaling is measured both in inviscid simulations and when using a Prandtl number $P=1$ in the viscous regime. △ Less

Submitted 12 October, 2016; originally announced October 2016.

Comments: Accepted for publication in Plasma Physics and Controlled Fusion

Quark deconfinement in the proto-magnetar model of Long Gamma-Ray Bursts

Authors: A. G. Pili, N. Bucciantini, A. Drago, G. Pagliara, L. Del Zanna

Abstract: We investigate the possible implications of quark deconfinement on the phenomenology of Long Gamma-Ray Bursts focusing, in particular, on the possibility to describe multiple prompt emission phases in the context of the proto-magnetar model. Starting from numerical models of rotating Hadron Stars and Quark Stars in full general relativity we track the electromagnetic spin-down evolution in both th… ▽ More We investigate the possible implications of quark deconfinement on the phenomenology of Long Gamma-Ray Bursts focusing, in particular, on the possibility to describe multiple prompt emission phases in the context of the proto-magnetar model. Starting from numerical models of rotating Hadron Stars and Quark Stars in full general relativity we track the electromagnetic spin-down evolution in both the hadronic and quark phase, linking the two families through conservation of baryon number and angular momentum. We give estimates of the timescales and the energetics involved in the spin-down process deriving, in the relevant spin range, the relation between the initial and the final masses and rotational energies, whenever hadron-quark conversion is possible. We show how the results can be used in relevant astrophysical cases such as the double burst GRB 110709B. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Comments: 5 pages, 3 figures, 2 tables

Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited

Authors: L. Del Zanna, E. Papini, S. Landi, M. Bugli, N. Bucciantini

Abstract: Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequ… ▽ More Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequent nonlinear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics, as appropriate in situations where the Alfven velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S^-1/2, with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S^-1/3 are considered, the so-called "ideal" tearing regime is retrieved, with modes growing independently on S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfven velocity. In the fully nonlinear stage we observe an inverse cascade towards the fundamental mode, with Petschek-type supersonic jets propagating at the external Alfven speed from the X-point, and a fast reconnection rate at the predicted value R~(ln S)^-1. △ Less

Submitted 20 May, 2016; originally announced May 2016.

Comments: 14 pages, 9 figures, accepted for publication (MNRAS)

The "ideal" tearing mode: theory and resistive MHD simulations

Authors: L. Del Zanna, S. Landi, E. Papini, F. Pucci, M. Velli

Abstract: Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the L… ▽ More Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L~S^(-1/2) and thus can be up to ~100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 10^5 to 10^7, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: Proceedings of ASTRONUM 2015, Avignon, France (11 pages, 6 figures)

Modeling the structure of magnetic fields in Neutron Stars: from the interior to the magnetosphere

Authors: N. Bucciantini, A. G. Pili, L. Del Zanna

Abstract: The phenomenology of the emission of pulsars and magnetars depends dramatically on the structure and properties of their magnetic field. In particular it is believed that the outbursting and flaring activity observed in AXPs and SRGs is strongly related to their internal magnetic field. Recent observations have moreover shown that charges are present in their magnetospheres supporting the idea tha… ▽ More The phenomenology of the emission of pulsars and magnetars depends dramatically on the structure and properties of their magnetic field. In particular it is believed that the outbursting and flaring activity observed in AXPs and SRGs is strongly related to their internal magnetic field. Recent observations have moreover shown that charges are present in their magnetospheres supporting the idea that their magnetic field is tightly twisted in the vicinity of the star. In principle these objects offer a unique opportunity to investigate physics in a regime beyond what can be obtained in the laboratory. We will discuss the properties of equilibrium models of magnetized neutron stars, and we will show how internal and external currents can be related. These magnetic field configurations will be discussed considering also their stability, relevant for their origin and possibly connected to events like SNe and GRBs. We will also show what kind of deformations they induce in the star, that could lead to emission of gravitational waves. In the case of a twisted magnetosphere we will show how the amount of twist regulates their general topology. A general formalism based on the simultaneous numerical solution of the general relativistic Grad-Shafranov equation and Einstein equations will be presented. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 9 pages, 3 figures, Proceedings of the 10th International Conference on Numerical Modeling of Space Plasma Flows, 8-12 June 2015, Avignon, France

Resistive magnetohydrodynamic simulations of the ideal tearing mode

Authors: Simone Landi, Luca Del Zanna, Emanuele Papini, Fulvia Pucci, Marco Velli

Abstract: We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield… ▽ More We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $S\to\infty$. Our simulations confirm that the tearing instability growth rate can be as fast as $γ\approx 0.6\,{τ_A}^{-1}$, where $τ_A$ is the ideal Alfvénic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the \emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvénic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas. △ Less

Submitted 27 April, 2015; originally announced April 2015.

Comments: 9 pages, 5 figures, to be published in Astrophysical Journal

Constraints on particle acceleration sites in the Crab Nebula from relativistic MHD simulations

Authors: Barbara Olmi, Luca Del Zanna, Elena Amato, Niccolò Bucciantini

Abstract: The Crab Nebula is one of the most efficient accelerators in the Galaxy and the only galactic source showing direct evidence of PeV particles. In spite of this, the physical process behind such effective acceleration is still a deep mystery. While particle acceleration, at least at the highest energies, is commonly thought to occur at the pulsar wind termination shock, the properties of the upstre… ▽ More The Crab Nebula is one of the most efficient accelerators in the Galaxy and the only galactic source showing direct evidence of PeV particles. In spite of this, the physical process behind such effective acceleration is still a deep mystery. While particle acceleration, at least at the highest energies, is commonly thought to occur at the pulsar wind termination shock, the properties of the upstream flow are thought to be non-uniform along the shock surface, and important constraints on the mechanism at work come from exact knowledge of where along this surface particles are being accelerated. Here we use axisymmetric relativistic MHD simulations to obtain constraints on the acceleration site(s) of particles of different energies in the Crab Nebula. Various scenarios are considered for the injection of particles responsible for synchrotron radiation in the different frequency bands, radio, optical and X-rays. The resulting emission properties are compared with available data on the multi wavelength time variability of the inner nebula. Our main result is that the X-ray emitting particles are accelerated in the equatorial region of the pulsar wind. Possible implications on the nature of the acceleration mechanism are discussed. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: 12 pages, 7 figures, 2 tables

The role of currents distribution in general relativistic equilibria of magnetized neutron stars

Authors: N. Bucciantini, A. G. Pili, L. Del Zanna

Abstract: Magnetic fields play a critical role in the phenomenology of neutron stars. There is virtually no observable aspect which is not governed by them. Despite this, only recently efforts have been done to model magnetic fields in the correct general relativistic regime, characteristic of these compact objects. In this work we present, for the first time a comprehensive and detailed parameter study, in… ▽ More Magnetic fields play a critical role in the phenomenology of neutron stars. There is virtually no observable aspect which is not governed by them. Despite this, only recently efforts have been done to model magnetic fields in the correct general relativistic regime, characteristic of these compact objects. In this work we present, for the first time a comprehensive and detailed parameter study, in general relativity, of the role that the current distribution, and the related magnetic field structure, have in determining the precise structure of neutron stars. In particular, we show how the presence of localized currents can modify the field strength at the stellar surface, and we look for general trends, both in terms of energetic properties, and magnetic field configurations. Here we verify that, among other things, for a large class of different current distributions the resulting magnetic configurations are always dominated by the poloidal component of the current. △ Less

Submitted 17 December, 2014; originally announced December 2014.

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

General relativistic neutron stars with twisted magnetosphere

Authors: A. G. Pili, N. Bucciantini, L. Del Zanna

Abstract: Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars are extreme manifestations of the most magnetized neutron stars: magnetars. The phenomenology of their emission and spectral properties strongly support the idea that the magnetospheres of these astrophysical objects are tightly twisted in the vicinity of the star. Previous studies on equilibrium configurations have so far focused on either the… ▽ More Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars are extreme manifestations of the most magnetized neutron stars: magnetars. The phenomenology of their emission and spectral properties strongly support the idea that the magnetospheres of these astrophysical objects are tightly twisted in the vicinity of the star. Previous studies on equilibrium configurations have so far focused on either the internal or the external magnetic field configuration, without considering a real coupling between the two fields. Here we investigate numerical equilibrium models of magnetized neutron stars endowed with a confined twisted magnetosphere, solving the general relativistic Grad-Shafranov equation both in the interior and in the exterior of the compact object. A comprehensive study of the parameters space is provided to investigate the effects of different current distributions on the overall magnetic field structure. △ Less

Submitted 12 December, 2014; originally announced December 2014.

Comments: 15 pages, 14 figures. Accepted for publication in MNRAS

Parametric decay of parallel and oblique Alfven waves in the expanding solar wind

Authors: L. Del Zanna, L. Matteini, S. Landi, A. Verdini, M. Velli

Abstract: The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case i… ▽ More The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case is also considered. In the oblique case, the direct excitation of daughter modes transverse to the local background field is found for the first time in an expanding environment, and this transverse cascade seems to be favored for monochromatic mother waves. The expansion effect reduces the instability growth rate, and it can even suppress its onset for the lowest frequency modes considered here, possibly explaining the persistence of these outgoing waves in the solar wind. △ Less

Submitted 22 July, 2014; originally announced July 2014.

Comments: 20 pages, 9 figures, accepted for publication in Joural of Plasma Physics

Journal ref: J. Plasma Phys. 81 (2015) 325810102

Numerical simulations of stellar jets and comparison between synthetic and observed maps: clues to the launch mechanism

Authors: Francesco Rubini, Lorenzo Maurri, Gabriele Inghirami, Francesca Bacciotti, Luca Del Zanna

Abstract: High angular resolution spectra obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS) provide rich morphological and kinematical information about the stellar jet phenomenon, which allows us to test theoretical models efficiently. In this work, numerical simulations of stellar jets in the propagation region are executed with the PLUTO code, by adopting inflow conditions that ari… ▽ More High angular resolution spectra obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS) provide rich morphological and kinematical information about the stellar jet phenomenon, which allows us to test theoretical models efficiently. In this work, numerical simulations of stellar jets in the propagation region are executed with the PLUTO code, by adopting inflow conditions that arise from former numerical simulations of magnetized outflows, accelerated by the disk-wind mechanism in the launching region. By matching the two regions, information about the magneto-centrifugal accelerating mechanism underlying a given astrophysical object can be extrapolated by comparing synthetic and observed position-velocity diagrams (PVDs). We show that quite different jets, like those from the young T Tauri stars DG-Tau and RW-Aur, may originate from the same disk-wind model for different configurations of the magnetic field at the disk surface. This result supports the idea that all the observed jets may be generated by the same mechanism. △ Less

Submitted 29 May, 2014; originally announced May 2014.

Comments: 15 pages, 18 figures, accepted for publication by A&A

Axisymmetric equilibrium models for magnetized neutron stars in General Relativity under the Conformally Flat Condition

Authors: A. G. Pili, N. Bucciantini, L. Del Zanna

Abstract: Extremely magnetized neutron stars with magnetic fields as strong as $\sim 10^{15-16}$ G, or magnetars, have received considerable attention in the last decade due to their identification as a plausible source for Soft Gamma Repeaters and Anomalous X-ray Pulsars. Moreover, this class of compact objects has been proposed as a possible engine capable of powering both Long and Short Gamma-Ray Bursts,… ▽ More Extremely magnetized neutron stars with magnetic fields as strong as $\sim 10^{15-16}$ G, or magnetars, have received considerable attention in the last decade due to their identification as a plausible source for Soft Gamma Repeaters and Anomalous X-ray Pulsars. Moreover, this class of compact objects has been proposed as a possible engine capable of powering both Long and Short Gamma-Ray Bursts, if the rotation period in their formation stage is short enough (~1 ms). Such strong fields are expected to induce substantial deformations of the star and thus to produce the emission of gravitational waves. Here we investigate, by means of numerical modeling, axisymmetric static equilibria of polytropic and strongly magnetized stars in full general relativity, within the ideal magneto-hydrodynamic regime. The eXtended Conformally Flat Condition (XCFC) for the metric is assumed, allowing us to employ the techniques introduced for the X-ECHO code [Bucciantini & Del Zanna, 2011, Astron. Astrophys. 528, A101], proven to be accurate, efficient, and stable. The updated XNS code for magnetized neutron star equilibria is made publicly available for the community (see www.arcetri.astro.it/science/ahead/XNS). Several sequences of models are here retrieved, from the purely toroidal (resolving a controversy in the literature) or poloidal cases, to the so-called twisted torus mixed configurations, expected to be dynamically stable, which are solved for the first time in the non-perturbative regime. △ Less

Submitted 6 February, 2014; v1 submitted 17 January, 2014; originally announced January 2014.

Comments: 24 pages, 22 figures, 5 tables

Dynamo action in thick disks around Kerr black holes: high-order resistive GRMHD simulations

Authors: M. Bugli, L. Del Zanna, N. Bucciantini

Abstract: We present the first kinematic study of an $αΩ$-dynamo in the General Relativistic Magneto-HydroDynamics (GRMHD) regime, applied to thick disks orbiting around Kerr black holes and using a fully covariant mean field dynamo closure for the Ohm law. We show that the $αΩ$-dynamo mechanism leads to a continuous exponential growth of the magnetic field within the disk and to the formation of dynamo wav… ▽ More We present the first kinematic study of an $αΩ$-dynamo in the General Relativistic Magneto-HydroDynamics (GRMHD) regime, applied to thick disks orbiting around Kerr black holes and using a fully covariant mean field dynamo closure for the Ohm law. We show that the $αΩ$-dynamo mechanism leads to a continuous exponential growth of the magnetic field within the disk and to the formation of dynamo waves drifting away or toward the equatorial plane. Since the evolution of the magnetic field occurs qualitatively in the same fashion as in the Sun, we present also butterfly diagrams that characterize our models and show the establishment of an additional timescale, which depends on the microscopic properties of the turbulent motions, possibly providing an alternative explanation to periodicities observed in many high-energy astrophysical sources where accretion onto a rotating black hole is believed to operate. △ Less

Submitted 21 January, 2014; v1 submitted 16 January, 2014; originally announced January 2014.

Comments: 5 pages, 4 figures. Accepted for publication in MNRAS

High-order schemes for non-ideal $3+1$ GRMHD: a study of the kinematic dynamo process in accretion tori

Authors: L. Del Zanna, M. Bugli, N. Bucciantini

Abstract: We present the first astrophysical application of the ECHO code in its recent version supplemented by a generalized Ohm law, namely a kinematic study of dynamo effects in thick accretion disks. High-order \emph{implicit-explicit} Runge-Kutta time-stepping routines are implemented and validated within $3+1$ \emph{General Relativistic MagnetoHydroDynamics} (GRMHD). The scheme is applied to a differe… ▽ More We present the first astrophysical application of the ECHO code in its recent version supplemented by a generalized Ohm law, namely a kinematic study of dynamo effects in thick accretion disks. High-order \emph{implicit-explicit} Runge-Kutta time-stepping routines are implemented and validated within $3+1$ \emph{General Relativistic MagnetoHydroDynamics} (GRMHD). The scheme is applied to a differentially rotating torus orbiting a Kerr black hole, where the mean-field dynamo process leads to strong amplification of seed magnetic fields. We show that the interplay between the toroidal and poloidal components occurs qualitatively in the same fashion as in the Sun, butterfly diagrams are reproduced, and a typical time-scale for the field evolution is found, depending on the dynamo and resistivity numbers, which could explain periodicities as observed in several accreting systems. △ Less

Submitted 14 January, 2014; originally announced January 2014.

Comments: 6 pages, Proceedings of ASTRONUM 2013, Biarritz. To appear as ASP Conference Series

Solving the 3+1 GRMHD equations in the eXtended Conformally Flat Condition: the XNS code for magnetized neutron stars

Authors: N. Bucciantini, A. G. Pili, L. Del Zanna

Abstract: High-energy phenomena in astrophysics involve quite generally a combination of relativistic motions and strong gravity. The simultaneous solution of Einstein equations and General Relativistic MHD equations is thus necessary to model with accuracy such phenomena. The so-called Conformally Flat Condition (CFC) allows a simplified treatment of Einstein equations, that can be particularly efficient i… ▽ More High-energy phenomena in astrophysics involve quite generally a combination of relativistic motions and strong gravity. The simultaneous solution of Einstein equations and General Relativistic MHD equations is thus necessary to model with accuracy such phenomena. The so-called Conformally Flat Condition (CFC) allows a simplified treatment of Einstein equations, that can be particularly efficient in those contexts where gravitational wave emission is negligible, like core-collapse, or the formation/evolution of neutron stars. We have developed a set of codes to model axisymmetric MHD flows, in General Relativity, where the solution of Einstein equations is achieved with a semi-spectral scheme. Here, we will show how this framework is particularly well suited to investigate neutron star equilibrium models in the presence of strong magnetic fields and we will present the XNS code, that has been recently developed and here updated to treat poloidal and mixed configurations. △ Less

Submitted 15 January, 2014; v1 submitted 14 January, 2014; originally announced January 2014.

Comments: 6 pages, 2 figures, Proceedings of the 8th International Conference on Numerical Modeling of Space Plasma Flows, 1-5 July 2013, Biarritz, France

General Relativistic Equilibrium Models of Magnetized Neutron Stars

Authors: A. G. Pili, N. Bucciantini, L. Del Zanna

Abstract: Magnetic fields play a crucial role in many astrophysical scenarios and, in particular, are of paramount importance in the emission mechanism and evolution of Neutron Stars (NSs). To understand the role of the magnetic field in compact objects it is important to obtain, as a first step, accurate equilibrium models for magnetized NSs. Using the conformally flat approximation we solve the Einstein's… ▽ More Magnetic fields play a crucial role in many astrophysical scenarios and, in particular, are of paramount importance in the emission mechanism and evolution of Neutron Stars (NSs). To understand the role of the magnetic field in compact objects it is important to obtain, as a first step, accurate equilibrium models for magnetized NSs. Using the conformally flat approximation we solve the Einstein's equations together with the GRMHD equations in the case of a static axisymmetryc NS taking into account different types of magnetic configuration. This allows us to investigate the effect of the magnetic field on global properties of NSs such as their deformation. △ Less

Submitted 19 December, 2013; v1 submitted 29 November, 2013; originally announced November 2013.

Comments: 4 pages, 6 figures, Proceedings of the 4th International Conference High Energy Phenomena in Relativistic Outflows "HEPRO IV", 23-26 July 2013, Heidelberg, Germany. To be publised in International Journal of Modern Physics: Conference Series

On the MHD modeling of the Crab Nebula radio emission

Authors: Barbara Olmi, Luca Del Zanna, Elena Amato, Rino Bandiera, Niccolò Bucciantini

Abstract: In recent years, it has become a well-established paradigm that many aspects of the physics of Pulsar Wind Nebulae (PWNe) can be fully accounted for within a relativistic MHD description. Numerical simulations have proven extremely successful in reproducing the X-ray morphology of the Crab Nebula, down to very fine detail. Radio emission, instead, is currently one of the most obscure aspects of th… ▽ More In recent years, it has become a well-established paradigm that many aspects of the physics of Pulsar Wind Nebulae (PWNe) can be fully accounted for within a relativistic MHD description. Numerical simulations have proven extremely successful in reproducing the X-ray morphology of the Crab Nebula, down to very fine detail. Radio emission, instead, is currently one of the most obscure aspects of the physics of these objects, and one that holds important information about pulsar properties and their role as antimatter factories. Here we address the question of radio emission morphology and integrated spectrum from the Crab Nebula, by using for the first time an axisymmetric dynamical model with parameters chosen to best reproduce its X-ray morphology. Based on our findings we discuss constraints on the origin of the radio emitting particles. △ Less

Submitted 31 October, 2013; originally announced October 2013.

Comments: 8 pages, 5 figures, 1 table

Relativistic viscous hydrodynamics for heavy-ion collisions with ECHO-QGP

Authors: L. Del Zanna, V. Chandra, G. Inghirami, V. Rolando, A. Beraudo, A. De Pace, G. Pagliara, A. Drago, F. Becattini

Abstract: We present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top of the \emph{Eulerian Conservative High-Order} astrophysical code for general relativistic magneto-hydrodynamics [\emph{Del Zanna et al., Astron. Astrophys. 473,… ▽ More We present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top of the \emph{Eulerian Conservative High-Order} astrophysical code for general relativistic magneto-hydrodynamics [\emph{Del Zanna et al., Astron. Astrophys. 473, 11, 2007}] and here it has been upgraded to handle the physics of the Quark-Gluon Plasma. ECHO-QGP features second-order treatment of causal relativistic viscosity effects in both Minkowskian or Bjorken coordinates; partial or complete chemical equilibrium of hadronic species before kinetic freeze-out; initial conditions based on the optical Glauber model, including a Monte-Carlo routine for event-by-event fluctuating initial conditions; a freeze-out procedure based on the Cooper-Frye prescription. The code is extensively validated against several test problems and results always appear accurate, as guaranteed by the combination of the conservative (shock-capturing) approach and the high-order methods employed. ECHO-QGP can be extended to include evolution of the electromagnetic fields coupled to the plasma. △ Less

Submitted 29 July, 2013; v1 submitted 30 May, 2013; originally announced May 2013.

Comments: 25 pages, two column, Final version: accepted for publication in European Physical Journal C

Journal ref: Eur. Phys. J. C (2013) 73: 2524

A fully covariant mean-field dynamo closure for numerical 3+1 resistive GRMHD

Authors: N. Bucciantini, L. Del Zanna

Abstract: The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong gravitational fields, and relativistic motions. In recent years numerical schemes for General Relativistic MHD (GRMHD) have been developed to model the multidimensional dy… ▽ More The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong gravitational fields, and relativistic motions. In recent years numerical schemes for General Relativistic MHD (GRMHD) have been developed to model the multidimensional dynamics of such systems, including the possibility of an evolving spacetime. Such schemes have been also extended beyond the ideal limit including the effects of resistivity, in an attempt to model dissipative physical processes acting on small scales (sub-grid effects) over the global dynamics. Along the same lines, magnetic fields could be amplified by the presence of turbulent dynamo processes, as often invoked to explain the high values of magnetization required in accretion disks and neutron stars. Here we present, for the first time, a further extension to include the possibility of a mean-field dynamo action within the framework of numerical 3+1 (resistive) GRMHD. A fully covariant dynamo closure is proposed, in analogy with the classical theory, assuming a simple alpha-effect in the comoving frame. Its implementation into a finite-difference scheme for GRMHD in dynamical spacetimes [the X-ECHO code: (Bucciantini and Del Zanna 2011)] is described, and a set of numerical test is presented and compared with analytical solutions wherever possible. △ Less

Submitted 14 September, 2012; v1 submitted 14 May, 2012; originally announced May 2012.

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

General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion

Authors: Olindo Zanotti, Constanze Roedig, Luciano Rezzolla, Luca Del Zanna

Abstract: We present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows wi… ▽ More We present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find significant differences with respect to purely hydrodynamical evolutions. In particular, once the system relaxes to a radiation-pressure dominated regime, the accretion rates become about two orders of magnitude smaller than in the purely hydrodynamical case, remaining however super-Eddington as are the luminosities. Furthermore, when increasing the Mach number of the inflowing gas, the accretion rates become smaller because of the smaller cross section of the black hole, but the luminosities increase as a result a stronger emission in the shocked regions. Overall, our approach provides the first self-consistent calculation of the Bondi-Hoyle luminosity, most of which is emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7, and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of computing luminosities self-consistently has also allowed us to compare with the bremsstrahlung luminosity often used in modelling the electromagnetic counterparts to supermassive black-hole binaries, to find that in the optically-thick regime these more crude estimates are about 20 times larger than our radiation-hydrodynamics results. △ Less

Submitted 10 March, 2015; v1 submitted 27 May, 2011; originally announced May 2011.

Comments: With updated bibliographyc information

Journal ref: Mon. Not. R. Astron. Soc. 417, 2899-2915 (2011)

GRMHD in axisymmetric dynamical spacetimes: the X-ECHO code

Authors: N. Bucciantini, L. Del Zanna

Abstract: We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an important role, such as for the evolution of magnetized neutron stars or for the gravitational collapse of the magnetized rotating cores of massive stars, whic… ▽ More We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an important role, such as for the evolution of magnetized neutron stars or for the gravitational collapse of the magnetized rotating cores of massive stars, which is the astrophysical scenario believed to eventually lead to (long) GRB events. The code is based on the extension of the Eulerian conservative high-order (ECHO) scheme [Del Zanna et al., A&A 473, 11 (2007)] for GRMHD, here coupled to a novel solver for the Einstein equations in the extended conformally flat condition (XCFC). We fully exploit the 3+1 Eulerian formalism, so that all the equations are written in terms of familiar 3D vectors and tensors alone, we adopt spherical coordinates for the conformal background metric, and we consider axisymmetric spacetimes and fluid configurations. The GRMHD conservation laws are solved by means of shock-capturing methods within a finite-difference discretization, whereas, on the same numerical grid, the Einstein elliptic equations are treated by resorting to spherical harmonics decomposition and solved, for each harmonic, by inverting band diagonal matrices. As a side product, we build and make available to the community a code to produce GRMHD axisymmetric equilibria for polytropic relativistic stars in the presence of differential rotation and a purely toroidal magnetic field. This uses the same XCFC metric solver of the main code and has been named XNS. Both XNS and the full X-ECHO codes are validated through several tests of astrophysical interest. △ Less

Submitted 26 January, 2011; v1 submitted 18 October, 2010; originally announced October 2010.

Comments: 18 pages, 9 figures, accepted for publication in A&A

Report number: NORDITA-2010-88

Journal ref: Astron. Astrophys. 528, A101 (2011)

EM counterparts of recoiling black holes: general relativistic simulations of non-Keplerian discs

Authors: Olindo Zanotti, Luciano Rezzolla, Luca Del Zanna, Carlos Palenzuela

Abstract: We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determini… ▽ More We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. We argue that the estimates of the bremsstrahlung luminosity computed without properly taking into account the radiation transfer yield cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic isothermal evolution, that the luminosity produced can reach a peak value above $L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system. △ Less

Submitted 12 September, 2010; v1 submitted 22 February, 2010; originally announced February 2010.

Comments: 17 pages, 11 figures, accepted by A&A, movies available at http://numrel.aei.mpg.de/Visualisations/Archive/BinaryBlackHoles/EMCounterparts/EMCounterparts.html

Journal ref: A&A 523, A8 (2010)

Polarization in Pulsar Wind Nebulae

Authors: D. Volpi, L. Del Zanna, E. Amato, N. Bucciantini

Abstract: The main goal of our present work is to provide, for the first time, a simple computational tool that can be used to compute the brightness, the spectral index, the polarization, the time variability and the spectrum of the non-thermal light (both synchrotron and inverse Compton, IC) associated with the plasma dynamics resulting from given relativistic magnetohydrodynamics (RMHD) simulations. Th… ▽ More The main goal of our present work is to provide, for the first time, a simple computational tool that can be used to compute the brightness, the spectral index, the polarization, the time variability and the spectrum of the non-thermal light (both synchrotron and inverse Compton, IC) associated with the plasma dynamics resulting from given relativistic magnetohydrodynamics (RMHD) simulations. The proposed method is quite general, and can be applied to any scheme for RMHD and to all non-thermal emitting sources, e.g. pulsar wind nebulae (PWNe), and in particular to the Crab Nebula (CN) as in the present proceeding. Here only the linear optical and X-ray polarization that characterizes the PWNe synchrotron emission is analyzed in order to infer information on the inner bulk flow structure, to provide a direct investigation of the magnetic field configuration, in particular the presence and the strength of a poloidal component, and to understand the origin of some emitting features, such as the knot, whose origins are still uncertain. The inverse Compton radiation is examined to disentangle the different contributions to radiation from the magnetic field and the particle energy distribution function, and to search for a possible hadronic component in the emitting PWN, and thus for the presence of ions in the wind. If hadronic radiation was found in a PWN, young supernova remnants would provide a natural birth-place of the cosmic-rays (CRs) up to the so-called knee in the CR spectrum. △ Less

Submitted 24 March, 2009; originally announced March 2009.

Comments: 6 pages, 2 figures, workshop: Polarimetry days in Rome: Crab status, theory and prospects, October 16-17 2008, Rome, Italy

Magnetized Relativistic Jets and Long-Duration GRBs from Magnetar Spindown during Core-Collapse Supernovae

Authors: N. Bucciantini, E. Quataert, B. D. Metzger, T. A. Thompson, J. Arons, L. Del Zanna

Abstract: We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of… ▽ More We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum, and energy are very similar to the corresponding free wind values (<~ 30% differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (~5-10 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN. △ Less

Submitted 23 January, 2009; originally announced January 2009.

Comments: 15 pages, 7 figures, submitted to MNRAS