-
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.
-
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.
-
Reverse phase transitions in binary neutron-star systems with exotic-matter cores
Authors:
Maximiliano Ujevic,
Henrique Gieg,
Federico Schianchi,
Swami Vivekanandji Chaurasia,
Ingo Tews,
Tim Dietrich
Abstract:
Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open…
▽ More
Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron star mergers, or if they appear at even higher densities that are not realized in the Cosmos. In this article, we perform numerical-relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiralling binary component. Our simulations reveal that shortly before the merger it is possible that such stars undergo a "reverse phase transition", i.e., densities decrease and the quark core inside the star disappears leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to a rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the post-merger phase transition, e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition.
△ Less
Submitted 5 January, 2023; v1 submitted 8 November, 2022;
originally announced November 2022.
-
Numerical Simulations of Dark Matter Admixed Neutron Star Binaries
Authors:
Mattia Emma,
Federico Schianchi,
Francesco Pannarale,
Violetta Sagun,
Tim Dietrich
Abstract:
Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the f…
▽ More
Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the first binary neutron star simulations of this kind. We find that the presence of dark matter reduces the lifetime of the merger remnant and favors a prompt collapse to a black hole. Furthermore, we find differences in the merger time for systems with the same total mass and mass ratio, but different amounts of dark matter. Finally, we find that electromagnetic signals produced by the merger of binary neutron stars admixed with dark matter are very unlikely to be as bright as their dark matter free counterparts. Given the increasing sensitivity of multi-messenger facilities, our analysis gives a new perspective on how to probe the presence of dark matter.
△ Less
Submitted 22 June, 2022;
originally announced June 2022.
-
Incorporating a radiative hydrodynamics scheme in the numerical-relativity code BAM
Authors:
Henrique Gieg,
Federico Schianchi,
Tim Dietrich,
Maximiliano Ujevic
Abstract:
To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, e.g., through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory…
▽ More
To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, e.g., through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory based Equations of State and a simple description of neutrino interactions through a Neutrino Leakage Scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB$Λφ$. Overall, our test simulations have good agreement with those reported in the literature.
△ Less
Submitted 10 July, 2022; v1 submitted 2 June, 2022;
originally announced June 2022.
-
Dynamical ejecta synchrotron emission as possible contributor to the changing behaviour of GRB170817A
Authors:
Vsevolod Nedora,
David Radice,
Sebastiano Bernuzzi,
Albino Perego,
Boris Daszuta,
Andrea Endrizzi,
Aviral Prakash,
Federico Schianchi
Abstract:
Over the past three years, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. (2021) indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible…
▽ More
Over the past three years, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. (2021) indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star (BNS) merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with a soft equation of state (EOS) and high mass ratio models with stiff EOS are disfavored as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in a good overall agreement with the data.
△ Less
Submitted 25 September, 2021; v1 submitted 9 April, 2021;
originally announced April 2021.
-
Predicting electromagnetic counterparts using low-latency, gravitational-wave data products
Authors:
Cosmin Stachie,
Michael W. Coughlin,
Tim Dietrich,
Sarah Antier,
Mattia Bulla,
Nelson Christensen,
Reed Essick,
Philippe Landry,
Benoit Mours,
Federico Schianchi,
Andrew Toivonen
Abstract:
Searches for gravitational-wave counterparts have been going in earnest since GW170817 and the discovery of AT2017gfo. Since then, the lack of detection of other optical counterparts connected to binary neutron star or black hole - neutron star candidates has highlighted the need for a better discrimination criterion to support this effort. At the moment, the low-latency gravitational-wave alerts…
▽ More
Searches for gravitational-wave counterparts have been going in earnest since GW170817 and the discovery of AT2017gfo. Since then, the lack of detection of other optical counterparts connected to binary neutron star or black hole - neutron star candidates has highlighted the need for a better discrimination criterion to support this effort. At the moment, the low-latency gravitational-wave alerts contain preliminary information about the binary properties and, hence, on whether a detected binary might have an electromagnetic counterpart. The current alert method is a classifier that estimates the probability that there is a debris disc outside the black hole created during the merger as well as the probability of a signal being a binary neutron star, a black hole - neutron star, a binary black hole or of terrestrial origin. In this work, we expand upon this approach to predict both the ejecta properties and provide contours of potential lightcurves for these events in order to improve follow-up observation strategy. The various sources of uncertainty are discussed, and we conclude that our ignorance about the ejecta composition and the insufficient constraint of the binary parameters, by the low-latency pipelines, represent the main limitations. To validate the method, we test our approach on real events from the second and third Advanced LIGO-Virgo observing runs.
△ Less
Submitted 28 June, 2021; v1 submitted 2 March, 2021;
originally announced March 2021.
-
Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers
Authors:
Vsevolod Nedora,
Federico Schianchi,
Sebastiano Bernuzzi,
David Radice,
Boris Daszuta,
Andrea Endrizzi,
Albino Perego,
Aviral Prakash,
Francesco Zappa
Abstract:
We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the bro…
▽ More
We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear equations of state lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.
△ Less
Submitted 11 November, 2021; v1 submitted 22 November, 2020;
originally announced November 2020.
-
Numerical Relativity Simulations of the Neutron Star Merger GW170817: Long-Term Remnant Evolutions, Winds, Remnant Disks, and Nucleosynthesis
Authors:
Vsevolod Nedora,
Sebastiano Bernuzzi,
David Radice,
Boris Daszuta,
Andrea Endrizzi,
Albino Perego,
Aviral Prakash,
Mohammadtaher Safarzadeh,
Federico Schianchi,
Domenico Logoteta
Abstract:
We present a systematic numerical-relativity study of the dynamical ejecta, winds and nucleosynthesis in neutron star merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOS) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamics turbulence up to 100 milliseconds postmer…
▽ More
We present a systematic numerical-relativity study of the dynamical ejecta, winds and nucleosynthesis in neutron star merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOS) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamics turbulence up to 100 milliseconds postmerger. Spiral density waves propagating from the neutron star remnant to the disk trigger a wind with mass flux ${\sim}0.1{-}0.5\,{\rm M_\odot/s}$ persisting for the entire simulation as long as the remnant does not collapse to black hole. This wind has average electron fraction $\gtrsim 0.3$ and average velocity ${\sim}0.1-0.17\,$c and thus is a site for the production of weak $r$-process elements (mass number $A<195$). Disks around long-lived remnants have masses ${\sim}0.1{-}0.2\,{\rm M_\odot}$, temperatures peaking at $\lesssim10\,$MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multi-component kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the comparable-mass long-lived mergers robustly accounts for all the $r$-process peaks, from mass number ${\sim}75$ to actinides in terms of solar abundances. Total abundandes are weakly dependent on the EOS, while the mass ratio affect the production of first peak elements.
△ Less
Submitted 12 August, 2020; v1 submitted 10 August, 2020;
originally announced August 2020.
-
Accretion-induced prompt black hole formation in asymmetric neutron star mergers, dynamical ejecta and kilonova signals
Authors:
Sebastiano Bernuzzi,
Matteo Breschi,
Boris Daszuta,
Andrea Endrizzi,
Domenico Logoteta,
Vsevolod Nedora,
Albino Perego,
Federico Schianchi,
David Radice,
Francesco Zappa,
Ignazio Bombaci,
Nestor Ortiz
Abstract:
We present new numerical relativity results of neutron star mergers with chirp mass $1.188M_\odot$ and mass ratios $q=1.67$ and $q=1.8$ using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model deri…
▽ More
We present new numerical relativity results of neutron star mergers with chirp mass $1.188M_\odot$ and mass ratios $q=1.67$ and $q=1.8$ using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model derived from ab-initio calculations based on the Brueckner-Hartree-Fock approach. We report for the first time evidence for accretion-induced prompt collapse in high-mass-ratio mergers, in which the tidal disruption of the companion and its accretion onto the primary star determine prompt black hole formation. As a result of the tidal disruption, an accretion disc of neutron-rich and cold matter forms with baryon masses ${\sim}0.15M_\odot$, and it is significantly heavier than the remnant discs in equal-masses prompt collapse mergers. Massive dynamical ejecta of order ${\sim}0.01M_\odot$ also originate from the tidal disruption. They are neutron rich and expand from the orbital plane with a crescent-like geometry. Consequently, bright, red and temporally extended kilonova emission is predicted from these mergers. Our results show that prompt black hole mergers can power bright electromagnetic counterparts for high-mass-ratio binaries, and that the binary mass ratio can be in principle constrained from multimessenger observations.
△ Less
Submitted 10 August, 2020; v1 submitted 12 March, 2020;
originally announced March 2020.
