-
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.
-
Nucleosynthesis in magnetorotational supernovae: impact of the magnetic field configuration
Authors:
M. Reichert,
M. Bugli,
J. Guilet,
M. Obergaulinger,
M. Á. Aloy,
A. Arcones
Abstract:
The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for \emph{extreme} supernova explosions, most likely h…
▽ More
The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for \emph{extreme} supernova explosions, most likely hosted by magnetorotational effects. We investigate the nucleosynthesis of five state-of-the-art magnetohydrodynamic models with fast rotation that have been previously calculated in full 3D and that involve an accurate neutrino transport (M1). One of the models does not contain any magnetic field and synthesizes elements around the iron group, in agreement with other CC-SNe models in literature. All other models host a strong magnetic field of the same intensity, but with different topology. For the first time, we investigate the nucleosynthesis of MR-SNe models with a quadrupolar magnetic field and a 90 degree tilted dipole. We obtain a large variety of ejecta compositions reaching from iron nuclei to nuclei up to the third r-process peak. We assess the robustness of our results by considering the impact of different nuclear physics uncertainties such as different nuclear masses, $β^{-}$-decays and $β^{-}$-delayed neutron emission probabilities, neutrino reactions, fission, and a feedback of nuclear energy on the temperature. We find that the qualitative results do not change with different nuclear physics input. The properties of the explosion dynamics and the magnetic field configuration are the dominant factors determining the ejecta composition.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
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.
-
Prospects for realtime characterization of core-collapse supernova and neutrino properties
Authors:
Meriem Bendahman,
Isabel Goos,
Joao Coelho,
Matteo Bugli,
Alexis Coleiro,
Sonia El Hedri,
Thierry Foglizzo,
Davide Franco,
Jérôme Guilet,
Antoine Kouchner,
Raphaël Raynaud,
Yahya Tayalati
Abstract:
Core-collapse supernovae (CCSNe) offer extremely valuable insights into the dynamics of galaxies. Neutrino time profiles from CCSNe, in particular, could reveal unique details about collapsing stars and particle behavior in dense environments. However, CCSNe in our galaxy and the Large Magellanic Cloud are rare and only one supernova neutrino observation has been made so far. To maximize the infor…
▽ More
Core-collapse supernovae (CCSNe) offer extremely valuable insights into the dynamics of galaxies. Neutrino time profiles from CCSNe, in particular, could reveal unique details about collapsing stars and particle behavior in dense environments. However, CCSNe in our galaxy and the Large Magellanic Cloud are rare and only one supernova neutrino observation has been made so far. To maximize the information obtained from the next Galactic CCSN, it is essential to combine analyses from multiple neutrino experiments in real time and transmit any relevant information to electromagnetic facilities within minutes. Locating the CCSN, in particular, is challenging, requiring disentangling CCSN localization information from observational features associated with the properties of the supernova progenitor and the physics of the neutrinos. Yet, being able to estimate the progenitor distance from the neutrino signal would be of great help for the optimisation of the electromagnetic follow-up campaign that will start soon after the propagation of the neutrino alert. Existing CCSN distance measurement algorithms based on neutrino observations hence rely on the assumption that neutrino properties can be described by the Standard Model. This paper presents a swift and robust approach to extract CCSN and neutrino physics information, leveraging diverse next-generation neutrino detectors to counteract potential measurement biases from Beyond the Standard Model effects.
△ Less
Submitted 23 November, 2023; v1 submitted 10 November, 2023;
originally announced November 2023.
-
A $4^{\rm th}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
Authors:
Vittoria Berta,
Andrea Mignone,
Matteo Bugli,
Giancarlo Mattia
Abstract:
We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exp…
▽ More
We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than one-dimensional finite volume ones), employs the generic upwind constrained transport averaging and sophisticated limiting strategies that include both a discontinuity detector and an order reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness of the method.
△ Less
Submitted 18 December, 2023; v1 submitted 18 October, 2023;
originally announced October 2023.
-
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.
-
Three-dimensional core-collapse supernovae with complex magnetic structures: II. Rotational instabilities and multi-messenger signatures
Authors:
Matteo Bugli,
Jérôme Guilet,
Thierry Foglizzo,
Martin Obergaulinger
Abstract:
The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low $T/\|W\|$ instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on…
▽ More
The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low $T/\|W\|$ instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on the growth of such instability has not yet been cleared. We analyze a series of three-dimensional magnetohydrodynamic models to characterize the effects of different magnetic configurations on the development of the low $T/\|W\|$ and the related multi-messenger features. In the absence of magnetic fields, we observe the growth on dynamical time scales of the low $T/\|W\|$, associated with a strong burst of GW and a correlated modulation of the neutrino emission. However, models with a strong magnetic field show a quenching of the low $T/\|W\|$, due to a flattening of the rotation profile in the first $\sim100$ ms after shock formation caused by the magnetic transport of angular momentum. The associated GW emission is weakened by an order of magnitude, exhibits a broader spectral shape, and has no dominant feature associated with the PNS large-scale oscillation modes. Neutrino luminosities are damped along the equatorial plane due to a more oblate PNS, and the only clear modulation in the signal is due to SASI activity. Finally, magnetized models produce lower luminosities for $ν_e$ than for $\barν_e$, which is connected to a higher concentration of neutron-rich material in the PNS surroundings.
△ Less
Submitted 9 February, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
MRI-driven dynamo at very high magnetic Prandtl numbers
Authors:
Jerome Guilet,
Alexis Reboul-Salze,
Raphael Raynaud,
Matteo Bugli,
Basile Gallet
Abstract:
The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to re…
▽ More
The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to resistivity). Despite its relevance to proto-neutron stars and neutron star merger remnants, the numerically challenging regime of high Pm is still largely unknown. Using zero-net flux shearing box simulations in the incompressible approximation, we studied MRI-driven dynamos at unprecedentedly high values of Pm reaching 256. The simulations show that the stress and turbulent energies are proportional to Pm up to moderately high values ($\mathrm{Pm} \sim 50$). At higher Pm, they transition to a new regime consistent with a plateau independent of Pm for $\rm Pm \gtrsim 100$. This trend is independent of the Reynolds number, which may suggest an asymptotic regime where the energy injection and dissipation are independent of the diffusive processes. Interestingly, large values of Pm not only lead to intense small-scale magnetic fields but also to a more efficient dynamo at the largest scales of the box.
△ Less
Submitted 1 September, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
MRI-driven $α-Ω$ dynamos in protoneutron stars
Authors:
Alexis Reboul-Salze,
Jérôme Guilet,
Raphaël Raynaud,
Matteo Bugli
Abstract:
Magnetars are highly magnetized neutron stars that can produce X-ray and soft gamma-ray emissions and that have a dipole of $10^{14}$ G to $10^{15}$ G. A promising mechanism for explaining magnetar formation is magnetic field amplification by the MRI in fast-rotating protoneutron stars (PNS). This scenario is supported by recent global models, which showed that small-scale turbulence can generate…
▽ More
Magnetars are highly magnetized neutron stars that can produce X-ray and soft gamma-ray emissions and that have a dipole of $10^{14}$ G to $10^{15}$ G. A promising mechanism for explaining magnetar formation is magnetic field amplification by the MRI in fast-rotating protoneutron stars (PNS). This scenario is supported by recent global models, which showed that small-scale turbulence can generate a dipole with magnetar-like intensity. However, the impact of buoyancy and density stratification on the efficiency of the MRI in generating a dipole is still unknown. We assess the impact of the density and entropy profiles on the MRI dynamo in a global model of a fast-rotating PNS, which focuses on its outer stratified region that is stable to convection. Using the pseudo-spectral code MagIC, we performed 3D Boussinesq and anelastic MHD simulations in spherical geometry with explicit diffusivities. We performed a parameter study in which we investigate the effect of different approximations and of thermal diffusion. We obtain a self-sustained turbulent MRI-driven dynamo. This confirms most of our previous incompressible results once rescaled for density. The MRI also generates a nondominant equatorial dipole, which represents about 4.3% of the averaged magnetic field strength. Interestingly, in the presence of a density gradient, an axisymmetric magnetic field at large scales oscillates with time, which can be described as a mean-field $αΩ$ dynamo. Buoyancy damps turbulence in the equatorial plane but it has overall a relatively weak influence with a realistic high thermal diffusion. Our results support the ability of the MRI to generate magnetar-like large-scale magnetic fields. They furthermore predict the presence of an $αΩ$ dynamo in the protoneutron star, which could be important to model in-situ magnetic field amplification in core-collapse supernovae. [abridged]
△ Less
Submitted 21 September, 2022; v1 submitted 3 November, 2021;
originally announced November 2021.
-
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.
-
Three-dimensional core-collapse supernovae with complex magnetic structures: I. Explosion dynamics
Authors:
Matteo Bugli,
Jérôme Guilet,
Martin Obergaulinger
Abstract:
Magnetic fields can play a major role in the dynamics of outstanding explosions associated to violent events such as GRBs and hypernovae, since they provide a natural mechanism to harness the rotational energy of the central proto-neutron star and power relativistic jets through the stellar progenitor. As the structure of such fields is quite uncertain, most numerical models of MHD-driven core-col…
▽ More
Magnetic fields can play a major role in the dynamics of outstanding explosions associated to violent events such as GRBs and hypernovae, since they provide a natural mechanism to harness the rotational energy of the central proto-neutron star and power relativistic jets through the stellar progenitor. As the structure of such fields is quite uncertain, most numerical models of MHD-driven core-collapse supernovae consider an aligned dipole as initial magnetic field, while the field's morphology can actually be much more complex. We present three-dimensional simulations of core-collapse supernovae with more realistic magnetic structures, such as quadrupolar fields and, for the first time, an equatorial dipolar field. Configurations other than an aligned dipole produce weaker explosions and less collimated outflows, but can at the same time be more efficient in extracting the rotational energy from the PNS. This energy is then stored in the surroundings of the PNS, rather than powering the polar jets. A significant axial dipolar component is also produced by models starting with a quadrupolar field, pointing to an effective dynamo mechanism operating in proximity of the PNS surface.
△ Less
Submitted 23 July, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
A global model of the magnetorotational instability in protoneutron stars
Authors:
Alexis Reboul-Salze,
Jérôme Guilet,
Raphaël Raynaud,
Matteo Bugli
Abstract:
Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MR…
▽ More
Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MRI at generating a dipole field is still unknown. To answer this question, we study the MRI dynamo in an idealized global model of a fast rotating protoneutron star with differential rotation. We perform 3D incompressible MHD simulations in spherical geometry with explicit diffusivities where the differential rotation is forced at the outer boundary. We vary the initial magnetic field and investigated different magnetic boundary conditions. These simulations were compared to local shearing box simulations. We obtain a self-sustained turbulent MRI-driven dynamo, whose saturated state is independent of the initial magnetic field. The MRI generates a strong turbulent magnetic field of $B \geq 2\times 10^{15}$ G and a non-dominant magnetic dipole, which represents systematically about $5\%$ of the averaged magnetic field strength. Interestingly, this dipole is tilted towards the equatorial plane. We find that local shearing box models can reproduce fairly well several characteristics of global MRI turbulence such as the kinetic and magnetic spectra. The turbulence is nonetheless more vigorous in the local models than in the global ones. Overall, our results support the ability of the MRI to form magnetar-like large-scale magnetic fields. They furthermore predict the presence of a stronger small-scale magnetic field. The resulting magnetic field could be important to power outstanding stellar explosions, such as superluminous supernovae and GRBs.
△ Less
Submitted 27 November, 2020; v1 submitted 7 May, 2020;
originally announced May 2020.
-
Honing and proofing Astrophysical codes on the road to Exascale. Experiences from code modernization on many-core systems
Authors:
Salvatore Cielo,
Luigi Iapichino,
Fabio Baruffa,
Matteo Bugli,
Christoph Federrath
Abstract:
The complexity of modern and upcoming computing architectures poses severe challenges for code developers and application specialists, and forces them to expose the highest possible degree of parallelism, in order to make the best use of the available hardware. The Intel$^{(R)}$ Xeon Phi$^{(TM)}$ of second generation (code-named Knights Landing, henceforth KNL) is the latest many-core system, whic…
▽ More
The complexity of modern and upcoming computing architectures poses severe challenges for code developers and application specialists, and forces them to expose the highest possible degree of parallelism, in order to make the best use of the available hardware. The Intel$^{(R)}$ Xeon Phi$^{(TM)}$ of second generation (code-named Knights Landing, henceforth KNL) is the latest many-core system, which implements several interesting hardware features like for example a large number of cores per node (up to 72), the 512 bits-wide vector registers and the high-bandwidth memory. The unique features of KNL make this platform a powerful testbed for modern HPC applications. The performance of codes on KNL is therefore a useful proxy of their readiness for future architectures. In this work we describe the lessons learnt during the optimisation of the widely used codes for computational astrophysics P-Gadget-3, Flash and Echo. Moreover, we present results for the visualisation and analysis tools VisIt and yt. These examples show that modern architectures benefit from code optimisation at different levels, even more than traditional multi-core systems. However, the level of modernisation of typical community codes still needs improvements, for them to fully utilise resources of novel architectures.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
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.
-
The impact of non-dipolar magnetic fields in core-collapse supernovae
Authors:
M. Bugli,
J. Guilet,
M. Obergaulinger,
P. Cerdá-Durán,
M. Á. Aloy
Abstract:
The magnetic field is believed to play an important role in at least some core-collapse supernovae if its magnitude reaches $10^{15}\,\rm{G}$, which is a typical value for a magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much mor…
▽ More
The magnetic field is believed to play an important role in at least some core-collapse supernovae if its magnitude reaches $10^{15}\,\rm{G}$, which is a typical value for a magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much more complex with strong multipolar and small-scale components and the consequences for the explosion are so far unclear. We investigate the effects of the magnetic field topology on the dynamics of core-collapse supernovae and the properties of forming proto-neutron star (PNS) by comparing pre-collapse fields of different multipolar orders and radial profiles. Using axisymmetric special relativistic MHD simulations and a two-moment neutrino transport, we find that higher multipolar magnetic configurations lead to generally less energetic explosions, slower expanding shocks and less collimated outflows. Models with a low order multipolar configuration tend to produce more oblate PNS, surrounded in some cases by a rotationally supported toroidal structure of neutron-rich material. Moreover, magnetic fields which are distributed on smaller angular scales produce more massive and faster rotating central PNS, suggesting that higher-order multipolar configurations tend to decrease the efficiency of the magnetorotational launching mechanism. Even if our dipolar models systematically display a far more efficient extraction of the rotational energy of the PNS, fields distributed on smaller angular scales are still capable of powering magnetorotational explosions and shape the evolution of the central compact object.
△ Less
Submitted 16 January, 2020; v1 submitted 6 September, 2019;
originally announced September 2019.
-
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.
-
ECHO-3DHPC: Advance the performance of astrophysics simulations with code modernization
Authors:
Matteo Bugli,
Luigi Iapichino,
Fabio Baruffa
Abstract:
We present recent developments in the parallelization scheme of ECHO-3DHPC, an efficient astrophysical code used in the modelling of relativistic plasmas. With the help of the Intel Software Development Tools, like Fortran compiler and Profile-Guided Optimization (PGO), Intel MPI library, VTune Amplifier and Inspector we have investigated the performance issues and improved the application scalabi…
▽ More
We present recent developments in the parallelization scheme of ECHO-3DHPC, an efficient astrophysical code used in the modelling of relativistic plasmas. With the help of the Intel Software Development Tools, like Fortran compiler and Profile-Guided Optimization (PGO), Intel MPI library, VTune Amplifier and Inspector we have investigated the performance issues and improved the application scalability and the time to solution. The node-level performance is improved by $2.3 \times$ and, thanks to the improved threading parallelisation, the hybrid MPI-OpenMP version of the code outperforms the MPI-only, thus lowering the MPI communication overhead.
△ Less
Submitted 10 October, 2018;
originally announced October 2018.
-
ECHO-3DHPC: relativistic accretion disks onto black holes
Authors:
Matteo Bugli
Abstract:
Current state-of-the-art simulations of accretion flows onto black holes require a significant level of numerical sophistication, in order to allow the three-dimensional modelling of relativistic magnetised plasma in a regime of strong gravity. We present here a new version of the GRMHD code \texttt{ECHO} developed in collaboration with the Max Planck Computing and Data Facility (MPCDF) and the Le…
▽ More
Current state-of-the-art simulations of accretion flows onto black holes require a significant level of numerical sophistication, in order to allow the three-dimensional modelling of relativistic magnetised plasma in a regime of strong gravity. We present here a new version of the GRMHD code \texttt{ECHO} developed in collaboration with the Max Planck Computing and Data Facility (MPCDF) and the Leibniz Rechenzentrum (LRZ), which employs a hybrid multidimensional MPI-OpenMP coupled with the production of MPI-HDF5 output files. The code's high degree of parallelisation has been crucial for the study of some fundamental properties of thick accretion disks around black holes, in particular the excitation of non-axisymmetric modes in presence of both hydrodynamic and magnetohydrodynamic instabilities.
△ Less
Submitted 8 October, 2018;
originally announced October 2018.
-
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.
-
rHARM: Accretion and Ejection in Resistive GR-MHD
Authors:
Qian Qian,
Christian Fendt,
Scott Noble,
Matteo Bugli
Abstract:
Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code…
▽ More
Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code rHARM are compared with an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magneto-rotational instability (MRI) in tori around black holes for a range of magnetic diffusivities. We find indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the black hole is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild black holes that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units.
△ Less
Submitted 14 October, 2016;
originally announced October 2016.
-
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.
-
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.
-
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.