-
Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes
Authors:
Liwei Ji,
Vassilios Mewes,
Yosef Zlochower,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Manuela Campanelli,
Federico Cipolletta,
Zachariah B. Etienne
Abstract:
Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of flui…
▽ More
Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition timestep limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double FFT filter and implementing it within the fully MPI-parallelized SphericalNR framework in the Einstein Toolkit. We demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the timestep significantly compared to unfiltered simulations.
△ Less
Submitted 9 November, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD
Authors:
Leonardo R. Werneck,
Zachariah B. Etienne,
Ariadna Murguia-Berthier,
Roland Haas,
Federico Cipolletta,
Scott C. Noble,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Bruno Giacomazzo,
Thiago Assumpção,
Joshua Faber,
Tanmayee Gupte,
Bernard J. Kelly,
Julian H. Krolik
Abstract:
We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp…
▽ More
We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit.
△ Less
Submitted 14 December, 2023; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations
Authors:
Federico G. Lopez Armengol,
Zachariah B. Etienne,
Scott C. Noble,
Bernard J. Kelly,
Leonardo R. Werneck,
Brendan Drachler,
Manuela Campanelli,
Federico Cipolletta,
Yosef Zlochower,
Ariadna Murguia-Berthier,
Lorenzo Ennoggi,
Mark Avara,
Riccardo Ciolfi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Julian H. Krolik,
Vassilios Mewes,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril,
Jeremy Schnittman
Abstract:
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg…
▽ More
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach.
△ Less
Submitted 31 October, 2022; v1 submitted 17 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.
-
HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage
Authors:
Ariadna Murguia-Berthier,
Scott C. Noble,
Luke F. Roberts,
Enrico Ramirez-Ruiz,
Leonardo R. Werneck,
Michael Kolacki,
Zachariah B. Etienne,
Mark Avara,
Manuela Campanelli,
Riccardo Ciolfi,
Federico Cipolletta,
Brendan Drachler,
Lorenzo Ennoggi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Bernard J. Kelly,
Julian H. Krolik,
Federico G. Lopez Armengol,
Ben Margalit,
Tim Moon,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril
, et al. (3 additional authors not shown)
Abstract:
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge…
▽ More
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.
△ Less
Submitted 31 January, 2022; v1 submitted 9 June, 2021;
originally announced June 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.
-
Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
Authors:
Vassilios Mewes,
Yosef Zlochower,
Manuela Campanelli,
Thomas W. Baumgarte,
Zachariah B. Etienne,
Federico G. Lopez Armengol,
Federico Cipolletta
Abstract:
We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-S…
▽ More
We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.
△ Less
Submitted 5 May, 2020; v1 submitted 14 February, 2020;
originally announced February 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.
-
On the last stable orbit around rapidly rotating neutron stars
Authors:
Federico Cipolletta,
Christian Cherubini,
Simonetta Filippi,
Jorge A. Rueda,
Remo Ruffini
Abstract:
We compute the binding energy and angular momentum of a test-particle at the last stable circular orbit (LSO) on the equatorial plane around a general relativistic, rotating neutron star (NS). We present simple, analytic, but accurate formulas for these quantities that fit the numerical results and which can be used in several astrophysical applications. We demonstrate the accuracy of these formul…
▽ More
We compute the binding energy and angular momentum of a test-particle at the last stable circular orbit (LSO) on the equatorial plane around a general relativistic, rotating neutron star (NS). We present simple, analytic, but accurate formulas for these quantities that fit the numerical results and which can be used in several astrophysical applications. We demonstrate the accuracy of these formulas for three different equations of state (EOS) based on nuclear relativistic mean-field theory models and argue that they should remain still valid for any NS EOS that satisfy current astrophysical constraints. We compare and contrast our numerical results with the corresponding ones for the Kerr metric characterized by the same mass and angular momentum.
△ Less
Submitted 14 June, 2017; v1 submitted 7 December, 2016;
originally announced December 2016.
-
Fast Rotating Neutron Stars with Realistic Nuclear Matter Equation of State
Authors:
Federico Cipolletta,
Christian Cherubini,
Simonetta Filippi,
Jorge A. Rueda,
Remo Ruffini
Abstract:
We construct equilibrium configurations of uniformly rotating neutron stars for selected relativistic mean-field nuclear matter equations of state (EOS). We compute in particular the gravitational mass ($M$), equatorial ($R_{\rm eq}$) and polar ($R_{\rm pol}$) radii, eccentricity, angular momentum ($J$), moment of inertia ($I$) and quadrupole moment ($M_2$) of neutron stars stable against mass-she…
▽ More
We construct equilibrium configurations of uniformly rotating neutron stars for selected relativistic mean-field nuclear matter equations of state (EOS). We compute in particular the gravitational mass ($M$), equatorial ($R_{\rm eq}$) and polar ($R_{\rm pol}$) radii, eccentricity, angular momentum ($J$), moment of inertia ($I$) and quadrupole moment ($M_2$) of neutron stars stable against mass-shedding and secular axisymmetric instability. By constructing the constant frequency sequence $f=716$ Hz of the fastest observed pulsar, PSR J1748-2446ad, and constraining it to be within the stability region, we obtain a lower mass bound for the pulsar, $M_{\rm min}=[1.2$-$1.4] M_\odot$, for the EOS employed. Moreover we give a fitting formula relating the baryonic mass ($M_b$) and gravitational mass of non-rotating neutron stars, $M_b/M_\odot=M/M_\odot+(13/200)(M/M_\odot)^2$ [or $M/M_\odot=M_b/M_\odot-(1/20)(M_b/M_\odot)^2$], which is independent on the EOS. We also obtain a fitting formula, although not EOS independent, relating the gravitational mass and the angular momentum of neutron stars along the secular axisymmetric instability line for each EOS. We compute the maximum value of the dimensionless angular momentum, $a/M\equiv c J/(G M^2)$ (or "Kerr parameter"), $(a/M)_{\rm max}\approx 0.7$, found to be also independent on the EOS. We compare and contrast then the quadrupole moment of rotating neutron stars with the one predicted by the Kerr exterior solution for the same values of mass and angular momentum. Finally we show that, although the mass quadrupole moment of realistic neutron stars never reaches the Kerr value, the latter is closely approached from above at the maximum mass value, as physically expected from the no-hair theorem. In particular the stiffer the EOS is, the closer the Kerr solution is approached.
△ Less
Submitted 22 June, 2015; v1 submitted 19 June, 2015;
originally announced June 2015.
-
Angular Momentum Role in the Hypercritical Accretion of Binary-Driven Hypernovae
Authors:
L. Becerra,
F. Cipolletta,
C. L. Fryer,
Jorge A. Rueda,
R. Ruffini
Abstract:
The induced gravitational collapse (IGC) paradigm explains a class of energetic, $E_{\rm iso}\gtrsim 10^{52}$~erg, long-duration gamma-ray bursts (GRBs) associated with Ic supernovae, recently named binary-driven hypernovae (BdHNe). The progenitor is a tight binary system formed of a carbon-oxygen (CO) core and a neutron star companion. The supernova ejecta of the exploding CO core triggers a hype…
▽ More
The induced gravitational collapse (IGC) paradigm explains a class of energetic, $E_{\rm iso}\gtrsim 10^{52}$~erg, long-duration gamma-ray bursts (GRBs) associated with Ic supernovae, recently named binary-driven hypernovae (BdHNe). The progenitor is a tight binary system formed of a carbon-oxygen (CO) core and a neutron star companion. The supernova ejecta of the exploding CO core triggers a hypercritical accretion process onto the neutron star, which reaches in a few seconds the critical mass, and gravitationally collapses to a black hole emitting a GRB. In our previous simulations of this process we adopted a spherically symmetric approximation to compute the features of the hypercritical accretion process. We here present the first estimates of the angular momentum transported by the supernova ejecta, $L_{\rm acc}$, and perform numerical simulations of the angular momentum transfer to the neutron star during the hyperaccretion process in full general relativity. We show that the neutron star: i) reaches in a few seconds either mass-shedding limit or the secular axisymmetric instability depending on its initial mass; ii) reaches a maximum dimensionless angular momentum value, $[c J/(G M^2)]_{\rm max}\approx 0.7$; iii) can support less angular momentum than the one transported by supernova ejecta, $L_{\rm acc} > J_{\rm NS,max}$, hence there is an angular momentum excess which necessarily leads to jetted emission.
△ Less
Submitted 7 September, 2015; v1 submitted 28 May, 2015;
originally announced May 2015.