Binary Neutron Star Merger Research Articles

Impact of pions on binary neutron star mergers

Impact of pions on binary neutron star mergers

Quenched Galaxies are Important Host Candidates of Binary Black Hole and Binary Neutron Star Mergers

In this work, we present the probabilities of mergers of binary black holes (BBHs) and binary neutron stars (BNSs) as functions of stellar mass, metallicity, specific star formation rate (sSFR), and age for galaxies with redshift z ≤ 0.1. Using the binary-star evolution (BSE) code and some fitting formulae, we construct a phenomenological model of cosmic gravitational wave (GW) merger events. By using the Bayesian analysis method and the observations from Advanced LIGO and Virgo, we obtain the relevant parameters of the phenomenological model (such as the maximum black hole mass is ). Combining the above model results with the galaxy catalog given by the EMERGE empirical galaxy model, we find the normalized probability of occurrence of a merger event varying with for galaxies with z ≤ 0.1 is different from that in previous studies, that is, two peaks exist in this work while there is only one peak (log10(sSFR/yr −1) = −10) in the previous work. The sSFR value corresponding to the new peak is log10(sSFR/yr −1) = −12 and in line with the value () of NGC 4493, the host galaxy of BNS merger event GW170817. The new peak is caused by today’s quenched galaxies, which give a large contribution to the total SFR at high redshift in the EMERGE empirical galaxy model. Moreover, we find that the BNS mergers are most likely detected in galaxies with age ∼11 Gyr, which is greater than previous results (6−8 Gyr) and close to the age of NGC 4993, Gyr.

Hydrodynamic Mixing of Accretion Disk Outflows in Collapsars: Implications for r-process Signatures

The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement—or even dominate—r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities of r-process elements, has not been observed. Recent work identified the distribution of r-process material within the SN ejecta as a key predictor of the ease with which signals associated with r-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds—the alleged locus of r-production in rare SNe—mix with initially r-process-free ejecta. We study mixing as a function of the wind mass, wind duration, and the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs of r-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence of r-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-driven r-process-enriched SN model is unlikely to explain standard energetic SNe.

Pre/Post-Merger Consistency Test for Gravitational Signals from Binary Neutron Star Mergers

Gravitational waves from binary neutron star (BNS) mergers can constrain nuclear models, predicting their equation of state (EOS). Matter effects on the inspiral-merger signal are encoded in the multipolar tidal polarizability parameters, whose leading order combination is sufficient to capture, with high accuracy, the key features of the merger waveform. Similar EOS-insensitive relations exist for the post-merger signal and can be used to model the emissions from the remnant. Several works suggested that the appearance of new degrees of freedom in high-density post-merger matter can be inferred by observing a violation of these EOS-insensitive relations. Here, we demonstrate a Bayesian method to test such an EOS-insensitive relation between the tidal polarizability parameters (or any other equivalent parameter) and the dominant post-merger frequency using information from the pre-and-post-merger signal. Technically, the method is similar to the inspiral-merger-ringdown consistency tests of General Relativity with binary black holes. However, differently from the latter, BNS pre/post-merger consistency tests are conceptually less informative and they only address the consistency of the assumed EOS-insensitive relation. Specifically, we discuss how such tests cannot conclusively discriminate between an EOS without respecting such a relation and the appearance of new degrees of freedom (or phase transitions) in high-density matter.

Opacities of singly and doubly ionized neodymium and uranium for kilonova emission modeling

ABSTRACT Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium (Z = 60) as well as the corresponding actinide uranium (Z = 92). Using two different codes [flexible atomic code (fac) and hartree–fock-relativistic (hfr)] for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the fac calculations, we optimize the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimization and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed. We find that the opacity of uranium is almost double the neodymium opacity.

Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation.

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017 s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Jetted and Turbulent Stellar Deaths: New LVK-detectable Gravitational-wave Sources

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ∼10–100 Hz, over a characteristic jet activity timescale ∼10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ∼10 events yr−1. The detection rate drops to ∼1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

Panning for gold, but finding helium: Discovery of the ultra-stripped supernova SN 2019wxt from gravitational-wave follow-up observations

We present the results from multi-wavelength observations of a transient discovered during an intensive follow-up campaign of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN 2019wxt, a young transient in a galaxy whose sky position (in the 80% GW contour) and distance (∼150 Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient’s tightly constrained age, its relatively faint peak magnitude (Mi ∼ −16.7 mag), and the r-band decline rate of ∼1 mag per 5 days appeared suggestive of a compact binary merger. However, SN 2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of only ∼0.1 M⊙, with 56Ni comprising ∼20% of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitor channels that could give rise to the observed properties of SN 2019wxt and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling genuine electromagnetic counterparts to GW events from transients such as SN 2019wxt soon after discovery is challenging: in a bid to characterise this level of contamination, we estimated the rate of events with a volumetric rate density comparable to that of SN 2019wxt and found that around one such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500 Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns.

Search for Ultra-high-energy Photons from Gravitational Wave Sources with the Pierre Auger Observatory

A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at ∼7 MeV cm−2 (time period 1000 s) and ∼35 MeV cm−2 (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total GW energy. These are the first limits on UHE photons from GW sources.

Superheavy Elements in Kilonovae

As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first “kilonova” associated with a binary neutron star merger (NSM) suggests that these sites are hosts of the rapid neutron capture (“r”) process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron-rich nuclei for which superheavy elements (Z ≥ 104) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy-element synthesis in nature. Favorable NSM conditions yield a mass fraction of superheavy elements X Z≥104 ≈ 3 × 10−2 at 7.5 hr post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.

Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations

Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of equations closely resembling the hydrodynamic ones. Solving the M1 scheme is still mathematically challenging, since it is necessary to model the radiation-matter interaction in regimes where the evolution equations become stiff and behave as an advection-diffusion problem. Here, we present different global, high-order time integration schemes based on Implicit-Explicit Runge-Kutta methods designed to overcome the time-step restriction caused by such behavior while allowing us to use the explicit Runge-Kutta commonly employed for the magneto-hydrodynamics and Einstein equations. Finally, we analyze their performance in several numerical tests.

Characterizing quasisteady states of fast neutrino-flavor conversion by stability and conservation laws

The question of what ingredients characterize the quasisteady state of fast neutrino-flavor conversion (FFC) is one of the longstanding riddles in neutrino oscillation. Addressing this issue is necessary for accurate modeling of neutrino transport in core-collapse supernova and binary neutron star merger. Recent numerical simulations of FFC have shown, however, that the quasisteady state is sensitively dependent on boundary conditions in space, and the physical reason for the dependence is not clear at present. In this study, we provide a physical interpretation of this issue based on arguments with stability and conservation laws. The stability can be determined by the disappearance of electron neutrino-lepton number--heavy-leptonic one (ELN-XLN) angular crossings, and we also highlight two conserved quantities characterizing the quasisteady state of FFC: (1) lepton number conservation along each neutrino trajectory and (2) conservation law associated with angular moments, depending on boundary conditions, for each flavor of neutrinos. We present an analytic prescription that matches the results of the nonlinear simulations presented in this work. This study represents a major step forward to a unified picture determining asymptotic states of FFCs.

Jet-environment interplay in magnetized binary neutron star mergers

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. 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 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 earlier 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.

Binary boson stars: Merger dynamics and formation of rotating remnant stars

Scalar boson stars have attracted attention as simple models for exploring the nonlinear dynamics of a large class of ultra compact and black hole mimicking objects. Here, we study the impact of interactions in the scalar matter making up these stars. In particular, we show the pivotal role the scalar phase and vortex structure play during the late inspiral, merger, and post-merger oscillations of a binary boson star, as well as their impact on the properties of the merger remnant. To that end, we construct constraint satisfying binary boson star initial data and numerically evolve the nonlinear set of Einstein-Klein-Gordon equations. We demonstrate that the scalar interactions can significantly affect the inspiral gravitational wave amplitude and phase, and the length of a potential hypermassive phase shortly after merger. If a black hole is formed after merger, we find its spin angular momentum to be consistent with similar binary black hole and binary neutron star merger remnants. Furthermore, we formulate a mapping that approximately predicts the remnant properties of any given binary boson star merger. Guided by this mapping, we use numerical evolutions to explicitly demonstrate, for the first time, that rotating boson stars can form as remnants from the merger of two non-spinning boson stars. We characterize this new formation mechanism and discuss its robustness. Finally, we comment on the implications for rotating Proca stars.

Joint analysis of gravitational-wave and electromagnetic data of mergers: breaking an afterglow model degeneracy in GW170817 and in future events

ABSTRACT On 2017 August 17, Advanced LIGO and Virgo observed GW170817, the first gravitational-wave (GW) signal from a binary neutron star merger. It was followed by a short-duration gamma-ray burst, GRB 170817A, and by a non-thermal afterglow emission. In this work, a combined simultaneous fit of the electromagnetic (EM, specifically, afterglow) and GW domains is implemented, both using the posterior distribution of a GW standalone analysis as prior distribution to separately process the EM data, and fitting the EM and GW domains simultaneously. These approaches coincide mathematically, as long as the actual posterior of the GW analysis, and not an approximation, is used as prior for the EM analysis. We treat the viewing angle, θv, as shared parameter across the two domains. In the afterglow modelling with a Gaussian structured jet this parameter and the jet core angle, θc, are correlated, leading to high uncertainties on their values. The joint EM + GW analysis relaxes this degeneracy, reducing the uncertainty compared to an EM-only fit. We also apply our methodology to hypothetical GW170817-like events occurring in the next GW observing run at ∼140 and 70 Mpc. At 70 Mpc the existing EM degeneracy is broken, thanks to the inclusion of the GW domain in the analysis. At 140 Mpc, the EM-only fit cannot constrain θv nor θc because of the lack of detections in the afterglow rising phase. Folding the GW data into the analysis leads to tighter constraints on θv, still leaving θc unconstrained, requiring instruments with higher sensitivities, such as Athena.

Compact object mergers: exploring uncertainties from stellar and binary evolution with sevn

ABSTRACT Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code sevn, which overcomes this issue by interpolating the main stellar properties from a set of pre-computed evolutionary tracks. We describe the new interpolation and adaptive time-step algorithms of sevn, and the main upgrades on single and binary evolution. With sevn, we evolved 1.2 × 109 binaries in the metallicity range 0.0001 ≤ Z ≤ 0.03, exploring a number of models for electron-capture, core-collapse and pair-instability supernovae, different assumptions for common envelope, stability of mass transfer, quasi-homogeneous evolution, and stellar tides. We find that stellar evolution has a dramatic impact on the formation of single and BCOs. Just by slightly changing the overshooting parameter (λov = 0.4, 0.5) and the pair-instability model, the maximum mass of a black hole can vary from ≈60 to ≈100 M⊙. Furthermore, the formation channels of BCOs and the merger efficiency we obtain with sevn show significant differences with respect to the results of other population-synthesis codes, even when the same binary-evolution parameters are used. For example, the main traditional formation channel of BCOs is strongly suppressed in our models: at high metallicity (Z ≳ 0.01) only <20 per cent of the merging binary black holes and binary neutron stars form via this channel, while other authors found fractions >70 per cent.

High-energy Neutrino Emission Associated with Gravitational-wave Signals: Effects of Cocoon Photons and Constraints on Late-time Emission

We investigate prospects for the detection of high-energy neutrinos produced in the prolonged jets of short gamma-ray bursts (sGRBs). The X-ray light curves of sGRBs show extended emission components lasting for 100–1000 s, which are considered to be produced by prolonged engine activity. Jets produced by such activity should interact with photons in the cocoon formed by the propagation of the jet inside the ejecta of neutron star mergers. We calculate neutrino emission from jets produced by prolonged engine activity, taking account of the interaction between photons provided from the cocoon and cosmic rays accelerated in the jets. We find that IceCube-Gen2, a future neutrino telescope, with second-generation gravitational-wave detectors will probably be able to observe neutrino signals associated with gravitational waves with around 10 years of operation, regardless of the assumed value of the Lorentz factor of the jets. Neutrino observations may enable us to constrain the dissipation region of the jets. We apply this model to GRB 211211A, a peculiar long GRB whose origin may be a binary neutron star merger. Our model predicts that IceCube is unlikely to detect any associated neutrinos, but a few similar events will be able to put a meaningful constraint on the physical quantities of the prolonged engine activities.

Anomalous coupling in radiation mediated shocks

We summarize recent attempts to unravel the role of plasma kinetic effects in radiation mediated shocks. Such shocks form in all strong stellar explosions and are responsible for the early electromagnetic emission released from these events. A key issue that has been overlooked in all previous works is the nature of the coupling between the charged leptons, that mediate the radiation force, and the ions, which are the dominant carriers of the shock energy. Our preliminary investigation indicates that in the case of relativistic shocks, as well as Newtonian shocks in multi-ion plasma, this coupling is driven by either, transverse magnetic fields of a sufficiently magnetized upstream medium, or plasma microturbulence if strong enough magnetic fields are absent. We discuss the implications for the shock breakout signal, as well as abundance evolution and kilonova emission in binary neutron star mergers.

Revisiting Vainshtein screening for fast N-body simulations

We revisit a method to incorporate the Vainshtein screening mechanism in N-body simulations proposed by R. Scoccimarro in [1]. We further extend this method to cover a subset of Horndeski theories that evade the bound on the speed of gravitational waves set by the binary neutron star merger GW170817. The procedure consists of the computation of an effective gravitational coupling that is time and scale dependent, G eff (k,z), where the scale dependence will incorporate the screening of the fifth-force. This is a fast procedure that when contrasted to the alternative of solving the full equation of motion for the scalar field inside N-body codes, reduces considerably the computational time and complexity required to run simulations. To test the validity of this approach in the non-linear regime, we have implemented it in a COmoving Lagrangian Approximation (COLA) N-body code, and ran simulations for two gravity models that have full N-body simulation outputs available in the literature, nDGP and Cubic Galileon. We validate the combination of the COLA method with this implementation of the Vainshtein mechanism with full N-body simulations for predicting the boost function: the ratio between the modified gravity non-linear matter power spectrum and its General Relativity counterpart. This quantity is of great importance for building emulators in beyond-ΛCDM models, and we find that the method described in this work has an agreement of below 2% for scales down to k ≈ 3h/Mpc with respect to full N-body simulations.

Nonparametric Representation of Neutron Star Equation of State Using Variational Autoencoder

We introduce a new nonparametric representation of the neutron star (NS) equation of state (EOS) by using the variational autoencoder (VAE). As a deep neural network, the VAE is frequently used for dimensionality reduction since it can compress input data to a low-dimensional latent space using the encoder component and then reconstruct the data using the decoder component. Once a VAE is trained, one can take the decoder of the VAE as a generator. We employ 100,000 EOSs that are generated using the nonparametric representation method based on Han et al. as the training set and try different settings of the neural network, then we get an EOS generator (the trained VAE’s decoder) with four parameters. We use the mass–tidal-deformability data of binary NS merger event GW170817, the mass–radius data of PSR J0030+0451, PSR J0740+6620, PSR J0437-4715, and 4U 1702-429, and the nuclear constraints to perform the Bayesian inference. The overall results of the analysis that includes all the observations are , , and (90% credible levels), where R 1.4/Λ1.4 are the radius/tidal deformability of a canonical 1.4 M ⊙ NS, and is the maximum mass of a nonrotating NS. The results indicate that the implementation of these VAE techniques can obtain reasonable results, while accelerating calculation by a factor of ∼3–10 or more, compared with the original method.