Hypermassive Neutron Star Research Articles

Modeling mergers of known galactic systems of binary neutron stars

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between 0.75 and 0.99. Specifically, these systems are J1756-2251, J0737-3039A, J1906 + 0746, B1534 + 12, J0453 + 1559 and B1913 + 16. We follow the dynamics of the merger from the late stage of the inspiral process up to ∼20ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems (q = 0.75, J0453 + 1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapses within the simulated time.

Equation of state effects and one-arm spiral instability in hypermassive neutron stars formed in eccentric neutron star mergers

We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for all equations of state in our survey. We find that mergers with stiffer equations of state tend to produce HMNSs with stronger m = 2 azimuthal mode density deformations, and weaker m = 1 components, relative to softer equations of state. We also find that for equations of state that can give rise to double-core HMNSs, large m = 1 density modes can already be present due to asymmetries in the two cores. This results in the generation of l = 2, m = 1 gravitational wave modes even before the dominance of a one-arm mode that ultimately arises following merger of the two cores. Our results further suggest that stiffer equations of state give rise to HMNSs generating lower m = 1 gravitational wave frequencies. Thus, if gravitational waves from the one-arm instability are detected, they could in principle constrain the neutron star equation of state. We estimate that, depending on the equation of state, the one-arm mode could potentially be detectable by second generation gravitational wave detectors at ∼10 Mpc and by third generation ones at ∼100 Mpc. Finally, we provide estimates of the properties of dynamical ejecta, as well as the accompanying kilonovae signatures.

Binary neutron star mergers and short gamma-ray bursts: Effects of magnetic field orientation, equation of state, and mass ratio

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Unequal mass binary neutron star mergers and multimessenger signals

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state.

How loud are neutron star mergers?

We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $\kappa^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB).

BINARY NEUTRON STAR MERGERS: A JET ENGINE FOR SHORT GAMMA-RAY BURSTS.

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, n = 1 polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter a/MBH ~ 0.74) immersed in a magnetized accretion disk. About 4000M ~ 60(MNS/1.625M⊙) ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow-an incipient jet-is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is Δ t ~ 0.1 (MNS/1.625M⊙) s. In contrast to black hole-neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars.

Modeling equal and unequal mass binary neutron star mergers using public codes

We present three-dimensional simulations of the dynamics of binary neutron star (BNS) mergers from the late stage of the inspiral process up to $\sim 20$ ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH). We investigate five equal-mass models of total gravitational mass $2.207$, $2.373$, $2.537$, $2.697$ and $2.854 M_\odot$, respectively, and four unequal mass models with $M_{\mathrm{ADM}}\simeq 2.53\ M_\odot$ and $q\simeq 0.94$, $0.88$, $0.82$, and $0.77$ (where $q = M^{(1)}/M^{(2)}$ is the mass ratio). We use a semi-realistic equation of state (EOS) namely, the seven-segment piece-wise polytropic SLyPP with a thermal component given by $\Gamma_{th} = 1.8$. We have also compared the resulting dynamics (for one model) using both, the BSSN-NOK and CCZ4 methods for the evolution of the gravitational sector, and also different reconstruction methods for the matter sector, namely PPM, WENO and MP5. Our results show agreement and high resolution, but superiority of BSSN-NOK supplemented by WENO reconstruction at lower resolutions. One of the important characteristics of the present investigation is that, for the first time, this has been done using only publicly available open source software, in particular, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. All of the source code and parameters used for the simulations have been made publicly available. This not only makes it possible to re-run and re-analyze our data; it also enables others to directly build upon this work for future research.

Relativistic simulations of eccentric binary neutron star mergers: One-arm spiral instability and effects of neutron star spin

We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star spin. For cases where the spin is aligned with the orbital angular momentum, the additional centrifugal support in the remnant hypermassive neutron star can prevent the prompt collapse to a black hole, while for antialigned cases the decreased total angular momentum can facilitate the collapse to a black hole. We show that even moderate spins can significantly increase the amount of ejected material, including the amount unbound with velocities greater than half the speed of light, leading to brighter electromagnetic signatures associated with kilonovae and interaction of the ejecta with the interstellar medium. Furthermore, we find that the initial neutron star spin can strongly affect the already rich phenomenology in the postmerger gravitational wave signatures that arise from the oscillation modes of the hypermassive neutron star. In several of our simulations, the resulting hypermassive neutron star develops the one-arm ($m=1$) spiral instability, the most pronounced cases being those with small but non-negligible neutron star spins. For long-lived hypermassive neutron stars, the presence of this instability leads to improved prospects for detecting these events through gravitational waves, and thus may give information about the neutron star equation of state.

One-arm spiral instability in hypermassive neutron stars formed by dynamical-capture binary neutron star mergers

Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular momentum. We find that this instability develops when post-merger fluid vortices lead to the generation of a toroidal remnant - a configuration whose maximum density occurs in a ring around the center-of-mass - with high vorticity along its rotation axis. The instability quickly saturates on a timescale of $\sim 10$ ms, with the $m=1$ azimuthal density multipole mode dominating over higher modes. The instability also leaves a characteristic imprint on the post-merger gravitational wave signal that could be detectable if the instability persists in long-lived remnants.

MONTE CARLO NEUTRINO TRANSPORT THROUGH REMNANT DISKS FROM NEUTRON STAR MERGERS

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.

Modeling the Complete Gravitational Wave Spectrum of Neutron Star Mergers.

In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency f2 related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that f(2) is correlated with the tidal coupling constant κ(2)^T that characterizes the binary tidal interactions during the late-inspiral merger. The relation f(2)(κ(2)^T) depends very weakly on the binary total mass, mass ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.

Effects of the microphysical equation of state in the mergers of magnetized neutron stars with neutrino cooling

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB 130603B and GRB 060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft equation of state (EOS) might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a subgrid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature postmerger within a time scale of $\ensuremath{\approx}10\text{ }\text{ }\mathrm{ms}$ but can have subtle effects on the postmerger dynamics.

Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

We study equal and unequal-mass neutron star mergers by means of new numerical relativity simulations in which the general relativistic hydrodynamics solver employs an algorithm that guarantees mass conservation across the refinement levels of the computational mesh. We consider eight binary configurations with total mass $M=2.7\,M_\odot$, mass-ratios $q=1$ and $q=1.16$, and four different equation of states (EOSs), and one configuration with a stiff EOS, $M=2.5M_\odot$ and $q=1.5$. We focus on the post-merger dynamics and study the merger remnant, dynamical ejecta and the postmerger gravitational wave spectrum. Although most of the merger remnants form a hypermassive neutron star collapsing to a black hole+disk system on dynamical timescales, stiff EOSs can eventually produce a stable massive neutron star. Ejecta are mostly emitted around the orbital plane; favored by large mass ratios and softer EOS. The postmerger wave spectrum is mainly characterized by non-axisymmetric oscillations of the remnant. The stiff EOS configuration consisting of a $1.5M_\odot$ and a $1.0M_\odot$ neutron star shows a rather peculiar dynamics. During merger the companion star is very deformed; about~$\sim0.03M_\odot$ of rest-mass becomes unbound from the tidal tail due torque; and the merger remnant forms stable neutron star surrounded by a massive accretion disk $\sim0.3M_\odot$. Similar configurations might be particularly interesting for electromagnetic counterparts. Comparing results obtained with and without the conservative mesh refinement algorithm, we find that post-merger simulations can be affected by systematic errors if mass conservation is not enforced in the mesh refinement strategy. However, mass conservation also depends on grid details and on the artificial atmosphere setup. [abridged]

Properties of hypermassive neutron stars formed in mergers of spinning binaries

We present numerical simulations of binary neutron star mergers, comparing irrotational binaries to binaries of NSs rotating aligned to the orbital angular momentum. For the first time, we study spinning BNSs employing nuclear physics equations of state, namely the ones of Lattimer and Swesty as well as Shen, Horowitz, and Teige. We study mainly equal mass systems leading to a hypermassive neutron star (HMNS), and analyze in detail its structure and dynamics. In order to exclude gauge artifacts, we introduce a novel coordinate system used for post-processing. The results for our equal mass models show that the strong radial oscillations of the HMNS modulate the instantaneous frequency of the gravitational wave (GW) signal to an extend that leads to separate peaks in the corresponding Fourier spectrum. In particular, the high frequency peaks which are often attributed to combination frequencies can also be caused by the modulation of the $m=2$ mode frequency in the merger phase. As a consequence for GW data analysis, the offset of the high frequency peak does not necessarily carry information about the radial oscillation frequency. Further, the low frequency peak in our simulations is dominated by the contribution of the plunge and the first 1--2 bounces. The amplitude of the radial oscillations depends on the initial NS spin, which therefore has a complicated influence on the spectrum. Another important result is that HMNSs can consist of a slowly rotating core with an extended, massive envelope rotating close to Keplerian velocity, contrary to the common notion that a rapidly rotating core is necessary to prevent a prompt collapse. Finally, our estimates on the amount of unbound matter show a dependency on the initial NS spin, explained by the influence of the latter on the amplitude of radial oscillations, which in turn cause shock waves.

Spectral properties of the post-merger gravitational-wave signal from binary neutron stars

Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.

Constraining the Equation of State of Neutron Stars from Binary Mergers

Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve this riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we have found that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the compactness of the progenitor stars. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We have tested our approach with simulated data and verified it works well for all the equations of state considered.

High resolution numerical relativity simulations for the merger of binary magnetized neutron stars

We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70 m, which yields the highest-resolution results among those derived so far. By an in-depth resolution study, we clarify several amplification mechanisms of magnetic fields during the binary neutron star merger for the first time. First, the Kelvin-Helmholtz instability developed in the shear layer at the onset of the merger significantly amplifies the magnetic fields. A hypermassive neutron star (HMNS) formed after the merger is then subject to the nonaxisymmetric magnetorotational instability, which amplifies the magnetic field in the HMNS. These two amplification mechanisms cannot be found with insufficient-resolution runs. We also show that the HMNS eventually collapses to a black hole surrounded by an accretion torus which is strongly magnetized at birth.

Constraining the equation of state of neutron stars from binary mergers.

Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve this riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we find that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the total compactness of the stars in the binary. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We test our approach with simulated data and verify it works well for all the equations of state considered.

Red or blue? A potential kilonova imprint of the delay until black hole formation following a neutron star merger

Mergers of binary neutron stars (NSs) usually result in the formation of a hypermassive neutron star (HMNS). Whether- and when this remnant collapses to a black hole (BH) depends primarily on the equation of state and on angular momentum transport processes, both of which are uncertain. Here we show that the lifetime of the merger remnant may be directly imprinted in the radioactively powered kilonova emission following the merger. We employ axisymmetric, time-dependent hydrodynamic simulations of remnant accretion disks orbiting a HMNS of variable lifetime, and characterize the effect of this delay to BH formation on the disk wind ejecta. When BH formation is relatively prompt (~ 100 ms), outflows from the disk are sufficiently neutron rich to form heavy r-process elements, resulting in ~ week-long emission with a spectral peak in the near-infrared (NIR), similar to that produced by the dynamical ejecta. In contrast, delayed BH formation allows neutrinos from the HMNS to raise the electron fraction in the polar direction to values such that potentially Lanthanide-free outflows are generated. The lower opacity would produce a brighter, bluer, and shorter-lived ~ day-long emission (a `blue bump') prior to the late NIR peak from the dynamical ejecta and equatorial wind. This new diagnostic of BH formation should be useful for events with a signal to noise lower than that required for direct detection of gravitational waveform signatures.

NECESSARY CONDITIONS FOR SHORT GAMMA-RAY BURST PRODUCTION IN BINARY NEUTRON STAR MERGERS

The central engine of short gamma-ray bursts (sGRBs) is hidden from direct view, operating at a scale much smaller than that probed by the emitted radiation. Thus we must infer its origin not only with respect to the formation of the {\it trigger} - the actual astrophysical configuration that is capable of powering a sGRB - but also from the consequences that follow from the various evolutionary pathways that may be involved in producing it. Considering binary neutron star mergers we critically evaluate, analytically and through numerical simulations, whether the neutrino-driven wind produced by the newly formed hyper-massive neutron star can allow the collimated relativistic outflow that follows its collapse to actually produce a sGRB or not. Upon comparison with the observed sGRB duration distribution, we find that collapse cannot be significantly delayed (<=100 ms) before the outflow is choked, thus limiting the possibility that long-lived hyper-massive remnants can account for these events. In the case of successful breakthrough of the jet through the neutrino-driven wind, the energy stored in the cocoon could contribute to the precursor and extended emission observed in sGRBs.