-
General Relativistic Magneto-Hydrodynamic Simulations with BAM: Implementation and Code Comparison
Authors:
Anna Neuweiler,
Tim Dietrich,
Bernd Brügmann,
Edoardo Giangrandi,
Kenta Kiuchi,
Federico Schianchi,
Philipp Mösta,
Swapnil Shankar,
Bruno Giacomazzo,
Masaru Shibata
Abstract:
Binary neutron star mergers are among the most energetic events in our Universe, with magnetic fields significantly impacting their dynamics, particularly after the merger. While numerical-relativity simulations that correctly describe the physics are essential to model their rich phenomenology, the inclusion of magnetic fields is crucial for realistic simulations. For this reason, we have extende…
▽ More
Binary neutron star mergers are among the most energetic events in our Universe, with magnetic fields significantly impacting their dynamics, particularly after the merger. While numerical-relativity simulations that correctly describe the physics are essential to model their rich phenomenology, the inclusion of magnetic fields is crucial for realistic simulations. For this reason, we have extended the BAM code to enable general relativistic magneto-hydrodynamic (GRMHD) simulations employing a hyperbolic `divergence cleaning' scheme. We present a large set of standard GRMHD tests and compare the BAM code to other GRMHD codes, SPRITZ, GRaM-X, and SACRA$_{\rm KK22}$, which employ different schemes for the evolution of the magnetic fields. Overall, we find that the BAM code shows a good performance in simple special-relativistic tests. In addition, we find good agreement and consistent results when comparing GRMHD simulation results between BAM and SACRA$_{\rm KK22}$.
△ Less
Submitted 30 July, 2024;
originally announced July 2024.
-
Black-hole formation in binary neutron star mergers: The impact of spin on the prompt-collapse scenario
Authors:
Federico Schianchi,
Maximiliano Ujevic,
Anna Neuweiler,
Henrique Gieg,
Ivan Markin,
Tim Dietrich
Abstract:
Accurate modeling of the multi-messenger signatures connected to binary neutron star mergers requires proper knowledge on the final remnant's fate and the conditions under which black holes (BHs) can form in such mergers. In this article, we use a suite of 84 numerical-relativity simulations in 28 different physical setups to explore the impact of the individual stars' spin on the merger outcome a…
▽ More
Accurate modeling of the multi-messenger signatures connected to binary neutron star mergers requires proper knowledge on the final remnant's fate and the conditions under which black holes (BHs) can form in such mergers. In this article, we use a suite of 84 numerical-relativity simulations in 28 different physical setups to explore the impact of the individual stars' spin on the merger outcome and on the early postmerger dynamics. We find that for setups close to the prompt-collapse threshold, the stars' intrinsic spin significantly changes the lifespan of the remnant before collapse and that the mass of the debris disk surrounding the BH is also altered. To enable a better understanding of BH formation, we check if there is at least a theoretical chance of observing densities that are above the maximum density allowed in a stable isolated neutron star, and we investigate the importance of different pressure contributions on the evolution of the postmerger remnant and BH formation.
△ Less
Submitted 1 March, 2024; v1 submitted 26 February, 2024;
originally announced February 2024.
-
Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817
Authors:
Shreya Anand,
Peter T. H. Pang,
Mattia Bulla,
Michael W. Coughlin,
Tim Dietrich,
Brian Healy,
Thomas Hussenot-Desenonges,
Theophile Jegou du Laz,
Mansi M. Kasliwal,
Nina Kunert,
Ivan Markin,
Kunal Mooley,
Vsevolod Nedora,
Anna Neuweiler
Abstract:
GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the i…
▽ More
GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta.
△ Less
Submitted 20 July, 2023;
originally announced July 2023.
-
M1 neutrino transport within the numerical-relativistic code BAM with application to low mass binary neutron star mergers
Authors:
Federico Schianchi,
Henrique Gieg,
Vsevolod Nedora,
Anna Neuweiler,
Maximiliano Ujevic,
Mattia Bulla,
Tim Dietrich
Abstract:
Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakag…
▽ More
Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme (M1) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes $\texttt{POSSIS}$ and $\texttt{Skynet}$, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
General-Relativistic Hydrodynamics Simulation of a Neutron Star - Sub-Solar-Mass Black Hole Merger
Authors:
Ivan Markin,
Anna Neuweiler,
Adrian Abac,
Swami Vivekanandji Chaurasia,
Maximiliano Ujevic,
Mattia Bulla,
Tim Dietrich
Abstract:
Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary sys…
▽ More
Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary system composed of a black hole with mass $0.5M_\odot$ and a neutron star with mass $1.4 M_\odot$. We investigate the merger dynamics of this exotic system and provide information about the connected gravitational-wave and kilonova signals. Our study indicates that current gravitational-waveform models are unable to adequately describe such systems and that phenomenological relations connecting the binary parameters with the ejecta and remnant properties are not applicable to our system. Furthermore, we find a dependence of the kilonova signal on the azimuthal viewing angle due to the asymmetric mass ejection. This first-of-its-kind simulation opens the door for the study of sub-solar mass black hole - neutron star mergers and could serve as a testing ground for future model development.
△ Less
Submitted 23 April, 2023;
originally announced April 2023.
-
Long-Term Simulations of Dynamical Ejecta: Homologous Expansion and Kilonova Properties
Authors:
Anna Neuweiler,
Tim Dietrich,
Mattia Bulla,
Swami Vivekanandji Chaurasia,
Stephan Rosswog,
Maximiliano Ujevic
Abstract:
Accurate numerical-relativity simulations are essential to study the rich phenomenology of binary neutron star systems. In this work, we focus on the material that is dynamically ejected during the merger process and on the kilonova transient it produces. Typically, radiative transfer simulations of kilonova light curves from ejecta make the assumption of homologous expansion, but this condition m…
▽ More
Accurate numerical-relativity simulations are essential to study the rich phenomenology of binary neutron star systems. In this work, we focus on the material that is dynamically ejected during the merger process and on the kilonova transient it produces. Typically, radiative transfer simulations of kilonova light curves from ejecta make the assumption of homologous expansion, but this condition might not always be met at the end of usually very short numerical-relativity simulations. In this article, we adjust the infrastructure of the BAM code to enable longer simulations of the dynamical ejecta with the aim of investigating when the condition of homologous expansion is satisfied. In fact, we observe that the deviations from a perfect homologous expansion are about 30% at roughly 100ms after the merger. To determine how these deviations might affect the calculation of kilonova light curves, we extract the ejecta data for different reference times and use them as input for radiative transfer simulations. Our results show that the light curves for extraction times later than 80ms after the merger deviate by less than 0.4mag and are mostly consistent with numerical noise. Accordingly, deviations from the homologous expansion for the dynamical ejecta component are negligible for the purpose of kilonova modelling.
△ Less
Submitted 31 August, 2022; v1 submitted 29 August, 2022;
originally announced August 2022.