-
AsterX: a new open-source GPU-accelerated GRMHD code for dynamical spacetimes
Authors:
Jay V. Kalinani,
Liwei Ji,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Lucas Timotheo Sanches,
Bing-Jyun Tsao,
Steven R Brandt,
Manuela Campanelli,
Riccardo Ciolfi,
Bruno Giacomazzo,
Roland Haas,
Erik Schnetter,
Yosef Zlochower
Abstract:
We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel…
▽ More
We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel applications. AsterX employs the Valencia formulation for GRMHD, coupled with the 'Z4c' formalism for spacetime evolution, while incorporating high resolution shock capturing schemes to accurately handle the hydrodynamics. AsterX has undergone rigorous testing in both static and dynamic spacetime, demonstrating remarkable accuracy and agreement with other codes in literature. Benchmarking the code through scaling tests on OLCF's Frontier supercomputer, we demonstrate a weak scaling efficiency of about 67%-77% on 4096 nodes compared to an 8-node performance.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
Jet-environment interplay in magnetized binary neutron star mergers
Authors:
Andrea Pavan,
Riccardo Ciolfi,
Jay Vijay Kalinani,
Andrea Mignone
Abstract:
GRB 170817A, the first short gamma-ray burst (sGRB) to be detected in coincidence with a gravitational wave signal, demonstrated that merging binary neutron star (BNS) systems can power collimated ultra-relativistic jets and, in turn, produce sGRBs. Moreover, it revealed that sGRB jets possess an intrinsic angular structure that is imprinted in the observable prompt and afterglow emission. Advance…
▽ More
GRB 170817A, the first short gamma-ray burst (sGRB) to be detected in coincidence with a gravitational wave signal, demonstrated that merging binary neutron star (BNS) systems can power collimated ultra-relativistic jets and, in turn, produce sGRBs. Moreover, it revealed that sGRB jets possess an intrinsic angular structure that is imprinted in the observable prompt and afterglow emission. Advanced numerical simulations represent the leading approach to investigate the physical processes underlying the evolution of sGRB jets breaking out of post-merger environments, and thus connect the final angular structure and energetics with specific jet launching conditions. In a previous paper, we carried out the first three-dimensional (3D) special-relativistic hydrodynamic simulations of incipient (top-hat) sGRB jets propagating across the realistic environment resulting from a general-relativistic (GR) hydrodynamic BNS merger simulation. While the above work marked an important step toward a consistent end-to-end description of sGRB jets from BNS mergers, those simulations did not account for the presence of magnetic fields, which are expected to play a key role. Here, we overcome this limitation, reporting the first 3D special-relativistic magnetohydrodynamic (MHD) simulation of a magnetized (structured and rotating) sGRB jet piercing through a realistic magnetized post-merger environment, wherein the initial conditions of the latter are directly imported from the outcome of a previous GRMHD BNS merger simulation.
△ Less
Submitted 14 June, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Turbulent magnetic field amplification in binary neutron star mergers
Authors:
C. Palenzuela,
R. Aguilera-Miret,
F. Carrasco,
R. Ciolfi,
J. V. Kalinani,
W. Kastaun,
B. Miñano,
D. Viganò
Abstract:
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insu…
▽ More
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.
△ Less
Submitted 3 July, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz
Authors:
Jay V. Kalinani,
Riccardo Ciolfi,
Wolfgang Kastaun,
Bruno Giacomazzo,
Federico Cipolletta,
Lorenzo Ennoggi
Abstract:
General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on `conserved' variables which then need to be converted back into the fundamenta…
▽ More
General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on `conserved' variables which then need to be converted back into the fundamental (`primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the early evolution (~50 ms) of a Fishbone-Moncrief BH-accrection disk system, show that RePrimAnd is able to support magnetized, low density environments with magnetic-to-fluid pressure ratios as high as 10^4, in situations where the previously used recovery scheme fails.
△ Less
Submitted 12 May, 2022; v1 submitted 22 July, 2021;
originally announced July 2021.
-
Short gamma-ray burst jet propagation in binary neutron star merger environments
Authors:
Andrea Pavan,
Riccardo Ciolfi,
Jay Vijay Kalinani,
Andrea Mignone
Abstract:
The multimessenger event GW170817/GRB 170817A confirmed that binary neutron star (BNS) mergers can produce short gamma-ray burst (SGRB) jets. This evidence promoted new investigations on the mechanisms through which a BNS merger remnant can launch such a powerful relativistic outflow and on the propagation of the latter across the surrounding post-merger environment. In particular, great strides h…
▽ More
The multimessenger event GW170817/GRB 170817A confirmed that binary neutron star (BNS) mergers can produce short gamma-ray burst (SGRB) jets. This evidence promoted new investigations on the mechanisms through which a BNS merger remnant can launch such a powerful relativistic outflow and on the propagation of the latter across the surrounding post-merger environment. In particular, great strides have been made in jet propagation models, establishing connections between the initial jet launching conditions, including the incipient jet launching time (with respect to merger) and the injection parameters, and the observable SGRB prompt and afterglow emission. However, present semi-analytical models and numerical simulations (with one notable exception) adopt simple hand-made prescriptions to account for the post-merger environment, lacking a direct association with any specific merging BNS system. Here, we present the first three-dimensional relativistic hydrodynamics simulations of incipient SGRB jets propagating through a post-merger environment that is directly imported from the outcome of a previous general relativistic BNS merger simulation. Our results show that the evolution and final properties of the jet can be largely affected by the anisotropies and the deviations from axisymmetry and homologous expansion characterizing more realistic BNS merger environments. In addition, we find that the inclusion of the gravitational pull from the central compact object, often overlooked, can have a major impact. Finally, we consider different jet launching times referred to the same BNS merger model and discuss the consequences for the ultimate jet properties.
△ Less
Submitted 25 June, 2021; v1 submitted 26 April, 2021;
originally announced April 2021.
-
Spritz: General Relativistic Magnetohydrodynamics with Neutrinos
Authors:
Federico Cipolletta,
Jay Vijay Kalinani,
Edoardo Giangrandi,
Bruno Giacomazzo,
Riccardo Ciolfi,
Lorenzo Sala,
Beatrice Giudici
Abstract:
We here present a new version of the publicly available general relativistic magnetohydrodynamic (GRMHD) code $\texttt{Spritz}$, which now includes an approximate neutrino leakage scheme able to handle neutrino cooling and heating. The leakage scheme is based on the publicly available $\texttt{ZelmaniLeak}$ code, with a few modifications in order to properly work with $\texttt{Spritz}$. We discuss…
▽ More
We here present a new version of the publicly available general relativistic magnetohydrodynamic (GRMHD) code $\texttt{Spritz}$, which now includes an approximate neutrino leakage scheme able to handle neutrino cooling and heating. The leakage scheme is based on the publicly available $\texttt{ZelmaniLeak}$ code, with a few modifications in order to properly work with $\texttt{Spritz}$. We discuss the involved equations, physical assumptions, and implemented numerical methods, along with a large battery of general relativistic tests performed with and without magnetic fields. Our tests demonstrate the correct implementation of the neutrino leakage scheme, paving the way for further improvements of our neutrino treatment and the first application to magnetized binary neutron star mergers. We also discuss the implementation in the $\texttt{Spritz}$ code of high-order methods for a more accurate evolution of hydrodynamical quantities.
△ Less
Submitted 12 March, 2021; v1 submitted 18 December, 2020;
originally announced December 2020.
-
Magnetically driven baryon winds from binary neutron star merger remnants and the blue kilonova of August 2017
Authors:
Riccardo Ciolfi,
Jay Vijay Kalinani
Abstract:
The observation of a radioactively-powered kilonova associated with the first binary neutron star (BNS) merger detected in gravitational waves proved that these events are ideal sites for the production of heavy r-process elements. However, the physical origin of the ejected material responsible for the early ("blue") and late ("red") components of this kilonova is still debated. Here, we investig…
▽ More
The observation of a radioactively-powered kilonova associated with the first binary neutron star (BNS) merger detected in gravitational waves proved that these events are ideal sites for the production of heavy r-process elements. However, the physical origin of the ejected material responsible for the early ("blue") and late ("red") components of this kilonova is still debated. Here, we investigate the possibility that the early/blue kilonova originated from the magnetically driven baryon wind launched after merger by the metastable neutron star remnant. Exploiting a magnetized BNS merger simulation with over 250 ms of post-merger evolution, we can follow for the first time the full mass ejection process up to its final decline. We find that the baryon wind carries $\simeq\!0.010\!-\!0.028\,M_\odot$ of unbound material, proving that the high mass estimated for the blue kilonova can be achieved. We also find expansion velocities of up to $\sim\!0.2\,c$, consistent with the lower end of the observational estimates, and we discuss possible effects neglected here that could further increase the final ejecta velocity. Overall, our results show that the magnetically driven baryon wind represents a viable channel to explain the blue kilonova.
△ Less
Submitted 25 August, 2020; v1 submitted 23 April, 2020;
originally announced April 2020.
-
Spritz: a new fully general-relativistic magnetohydrodynamic code
Authors:
Federico Cipolletta,
Jay Vijay Kalinani,
Bruno Giacomazzo,
Riccardo Ciolfi
Abstract:
The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and electromagnetic radiation, from radio to $γ$-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical…
▽ More
The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and electromagnetic radiation, from radio to $γ$-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical mechanisms involved in such systems, it is necessary to employ fully general relativistic magnetohydrodynamic (GRMHD) simulations able to include the effects of a composition and temperature dependent equation of state describing neutron star matter as well as neutrino emission and reabsorption. Here, we present our new code named $\texttt{Spritz}$ that solves the GRMHD equations in 3D Cartesian coordinates and on a dynamical spacetime. The code can support tabulated equations of state, taking into account finite temperature effects and allowing for the inclusion of neutrino radiation. In this first paper, we present the general features of the code and a series of tests performed in special and general relativity to assess the robustness of the basic GRMHD algorithms implemented. Among these tests, we also present the first comparison between a non-staggered and a staggered formulation of the vector potential evolution, which is used to guarantee the divergence-less character of the magnetic field. With respect to other publicly available GRMHD codes, $\texttt{Spritz}$ combines the robust approach of a staggered formulation of the vector potential together with the use of an equation of state driver ($\texttt{EOS_Omni}$) that can allow the code to use finite temperature equations of state. A next version of the code will fully test the $\texttt{EOS_Omni}$ driver by coupling it with a neutrino leakage scheme.
△ Less
Submitted 15 June, 2020; v1 submitted 10 December, 2019;
originally announced December 2019.
-
First 100 ms of a long-lived magnetized neutron star formed in a binary neutron star merger
Authors:
Riccardo Ciolfi,
Wolfgang Kastaun,
Jay Vijay Kalinani,
Bruno Giacomazzo
Abstract:
The recent multimessenger observation of the short gamma-ray burst (SGRB) GRB 170817A together with the gravitational wave (GW) event GW170817 provides evidence for the long-standing hypothesis associating SGRBs with binary neutron star (BNS) mergers. The nature of the remnant object powering the SGRB, which could have been either an accreting black hole (BH) or a long-lived magnetized neutron sta…
▽ More
The recent multimessenger observation of the short gamma-ray burst (SGRB) GRB 170817A together with the gravitational wave (GW) event GW170817 provides evidence for the long-standing hypothesis associating SGRBs with binary neutron star (BNS) mergers. The nature of the remnant object powering the SGRB, which could have been either an accreting black hole (BH) or a long-lived magnetized neutron star (NS), is, however, still uncertain. General relativistic magnetohydrodynamic (GRMHD) simulations of the merger process represent a powerful tool to unravel the jet launching mechanism, but so far most simulations focused the attention on a BH as the central engine, while the long-lived NS scenario remains poorly investigated. Here, we explore the latter by performing a GRMHD BNS merger simulation extending up to ~100 ms after merger, much longer than any previous simulation of this kind. This allows us to (i) study the emerging structure and amplification of the magnetic field and observe a clear saturation at magnetic energy $E_\mathrm{mag} \sim 10^{51}$ erg, (ii) follow the magnetically supported expansion of the outer layers of the remnant NS and its evolution into an ellipsoidal shape without any surrounding torus, and (iii) monitor density, magnetization, and velocity along the axis, observing no signs of jet formation. We also argue that the conditions at the end of the simulation disfavor later jet formation on subsecond timescales if no BH is formed. Furthermore, we examine the rotation profile of the remnant, the conversion of rotational energy associated with differential rotation, the overall energy budget of the system, and the evolution of the GW frequency spectrum. Finally, we perform an additional simulation where we induce the collapse to a BH ~70 ms after merger, in order to gain insights on the prospects for massive accretion tori in case of a late collapse. We find that...
△ Less
Submitted 18 July, 2019; v1 submitted 23 April, 2019;
originally announced April 2019.