Neutron Star Merger Remnants Research Articles

Binary neutron star mergers in massive scalar-tensor theory: Properties of postmerger remnants

We investigate the properties of postmerger remnants of binary neutron star mergers in the framework of Damour–Esposito-Farese-type scalar-tensor theory of gravity with a massive scalar field by numerical relativity simulation. It is found that the threshold mass for prompt collapse is raised in the presence of the excited scalar field. Our simulation results also suggest the existence of a long-lived ϕ mode in hypermassive neutron stars due to the presence of the massive scalar field that enhances the quasiradial oscillation in the remnant. We investigate the descalarization condition in hypermassive neutron stars and discover a distinctive signature in postmerger gravitational waves. Published by the American Physical Society 2024

Incidence of afterglow plateaus in gamma-ray bursts associated with binary neutron star mergers

One of the most surprising gamma-ray burst (GRB) features discovered with the Swift X-ray telescope (XRT) is a plateau phase in the early X-ray afterglow light curves. These plateaus are observed in the majority of long GRBs, while their incidence in short GRBs (SGRBs) is still uncertain due to their fainter X-ray afterglow luminosity with respect to long GRBs. An accurate estimate of the fraction of SGRBs with plateaus is of utmost relevance given the implications that the plateau may have for our understanding of the jet structure and possibly of the nature of the binary neutron star (BNS) merger remnant. This work presents the results of an extensive data analysis of the largest and most up-to-date sample of SGRBs observed with the XRT, and for which the redshift has been measured. We find a plateau incidence of 18-37<!PCT!> in SGRBs, which is a significantly lower fraction than that measured in long GRBs (>50<!PCT!>). Although still debated, the plateau phase could be explained as energy injection from the spin-down power of a newly born magnetized neutron star (NS; magnetar). We show that this scenario can nicely reproduce the observed short GRB (SGRBs) plateaus, while at the same time providing a natural explanation for the different plateau fractions between short and long GRBs. In particular, our findings may imply that only a minority of BNS mergers generating SGRBs leave behind a sufficiently stable or long-lived NS to form a plateau. From the probability distribution of the BNS remnant mass, a fraction 18-37<!PCT!> of short GRB plateaus implies a maximum NS mass in the range $ odot

Magnetic dissipation in short gamma-ray-burst jets

Aims. Short gamma-ray bursts originate when relativistic jets emerge from the remnants of binary neutron star (BNS) mergers, as observed in the first multi-messenger event GW170817–GRB 170817A, which coincided with a gravitational wave signal. Both the jet and the remnant are believed to be magnetized, and the presence of magnetic fields is known to influence the jet propagation across the surrounding post-merger environment. In the magnetic interplay between the jet and the environment itself, effects due to a finite plasma conductivity may be important, especially in the first phases of jet propagation. We aim to investigate these effects, from jet launching to its final breakout from the post-merger environment. Methods. We used the PLUTO numerical code to perform 2D axisymmetric and full 3D resistive relativistic magnetohydrodynamic (MHD) simulations, employing spherical coordinates with spatial radial stretching. We considered and compared different models for physical resistivity, which must be small but still dominating over the intrinsic numerical dissipation (which yields unwanted smearing of structures in any ideal MHD code). Stiff terms in the current density are treated with IMplicit-EXplicit Runge Kutta algorithms for time-stepping. A Synge-like gas (Taub equation of state) is also considered. All simulations were performed using an axisymmetric analytical model for both the jet propagation environment and the jet injection; we leave the case of jet propagation in a realistic environment (i.e., imported from an actual BNS merger simulation) to a future study. Results. As expected, no qualitative differences are detected due to the effect of a finite conductivity, but significant quantitative differences in the jet structure and induced turbulence are clearly seen in 2D axisymmetric simulations, and we also compare different resistivity models. We see the formation of regions with a resistive electric field parallel to the magnetic field, and nonthermal particle acceleration may be enhanced there. The level of dissipated Ohmic power is also dependent on the various recipes for resistivity. Most of the differences arise before the breakout from the inner environment, whereas once the jet enters the external weakly magnetized environment region, these differences are preserved during further propagation despite the lower grid refinement. Finally, we show and discuss the 3D evolution of the jet within the same environment in order to highlight the emergence of nonaxisymmetric features.

Energy-dependent and Energy-integrated Two-moment General-relativistic Neutrino Transport Simulations of a Hypermassive Neutron Star

We compare two-moment-based energy-dependent and three variants of energy-integrated neutrino transport general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. To study the impacts due to the choice of the neutrino transport schemes, we perform simulations with the same setups and input neutrino microphysics. We show that the main differences between energy-dependent and energy-integrated neutrino transport are found in the disk and ejecta properties, as well as in the neutrino signals. The properties of the disk surrounding the neutron star and the ejecta in energy-dependent transport are very different from the ones obtained using energy-integrated schemes. Specifically, in the energy-dependent case, the disk is more neutron-rich at early times and becomes geometrically thicker at later times. In addition, the ejecta is more massive and, on average, more neutron-rich in the energy-dependent simulations. Moreover, the average neutrino energies and luminosities are about 30% higher. Energy-dependent neutrino transport is necessary if one wants to better model the neutrino signals and matter outflows from neutron star merger remnants via numerical simulations.

Radioactive Gamma-Ray Lines from Long-lived Neutron Star Merger Remnants

The observation of a kilonova AT2017gfo associated with the gravitational wave event GW170817 provides the first strong evidence that neutron star mergers are dominant contributors to the production of heavy r-process elements. Radioactive gamma-ray lines emitted from neutron star merger remnants provide a unique probe for investigating the nuclide composition and tracking its evolution. In this work, we studied the gamma-ray line features arising from the radioactive decay of heavy nuclei in the merger remnants based on the r-process nuclear reaction network and the astrophysical inputs derived from numerical relativity simulations. The decay chain of 50126 Sn (T 1/2 = 230 kyr) → 51126 Sb (T 1/2 = 12.35 days) → 52126 Te (stable) produces several bright gamma-ray lines with energies of 415, 667, and 695 keV, making it the most promising decay chain during the remnant phase. The photon fluxes of these bright gamma-ray lines reach ∼10−5 γ cm−2 s−1 for Galactic merger remnants with ages less than 100 kyr, which can be detected by the high energy resolution MeV gamma-ray detectors like the MASS mission.

Jetlike structures in low-mass binary neutron star merger remnants

Jetlike structures in low-mass binary neutron star merger remnants

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

Symmetry breaking due to multi-angle matter-neutrino resonance in neutron star merger remnants

Neutron star merger remnants are unique sites for exploring neutrino flavor conversion in dense media. Because of the natural excess of ν̅e over νe , the neutrino-neutrino potential can cancel the matter potential, giving rise to matter-neutrino resonant flavor conversion. Under the assumption of two (anti)neutrino flavors and spatial homogeneity, we solve the neutrino quantum kinetic equations to investigate the occurrence of the matter-neutrino resonance within a multi-angle framework. We find that isotropy is broken spontaneously, regardless of the mass ordering. Relying on a hydrodynamical simulation of a binary neutron star merger remnant with a black hole of 3 M ⊙ and an accretion torus of 0.3 M ⊙, we find that complete flavor conversion caused by the matter-neutrino resonance is unlikely, although the matter and neutrino potentials cancel at various locations above the disk. Importantly, the matter-neutrino resonant flavor conversion crucially depends on the shape of the neutrino angular distributions. Our findings suggest that an accurate modeling of the neutrino angular distributions is necessary to understand flavor conversion physics in merger remnants, its implications on the disk physics and synthesis of the elements heavier than iron.

Secular Outflows from Long-Lived Neutron Star Merger Remnants

We study mass ejection from a binary neutron star merger producing a long-lived massive neutron star remnant with general-relativistic neutrino-radiation hydrodynamics simulations. In addition to outflows generated by shocks and tidal torques during and shortly after the merger, we observe the appearance of a wind driven by spiral density waves in the disk. This spiral-wave-driven outflow is predominantly located close to the disk orbital plane and have a broad distribution of electron fractions. At higher latitudes, a high electron-fraction wind is driven by neutrino radiation. The combined nucleosynthesis yields from all the ejecta components is in good agreement with Solar abundance measurements.

Thermal behavior as indicator for hyperons in binary neutron star merger remnants

Thermal behavior as indicator for hyperons in binary neutron star merger remnants

Magnetized Outflows from Short-lived Neutron Star Merger Remnants Can Produce a Blue Kilonova

We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino–matter interactions and is well suited to studying the resulting nucleosynthesis and kilonova emission. A magnetized wind is driven from the remnant and ejects neutron-rich material at a quasi-steady-state rate of 0.8 × 10−1 M ⊙s−1. We find that the ejecta in our simulations underproduce r-process abundances beyond the second r-process peak. For sufficiently long-lived remnants, these outflows alone can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.

Jets from Neutron-Star Merger Remnants and Massive Blue Kilonovae.

We perform three-dimensional general-relativistic magnetohydrodynamic simulations with weak interactions of binary neutron-star (BNS) mergers resulting in a long-lived remnant neutron star, with properties typical of galactic BNS and consistent with those inferred for the first observed BNS merger GW170817. We demonstrate self-consistently that within ≲30 ms postmerger magnetized (σ∼5-10) incipient jets emerge with asymptotic Lorentz factor Γ∼5-10, which successfully break out from the merger debris within ≲20 ms. A fast (v≲0.6c), magnetized (σ∼0.1) wind surrounds the jet core and generates a UV/blue kilonova precursor on timescales of hours, similar to the precursor signal due to free neutron decay in fast dynamical ejecta. Postmerger ejecta are quickly dominated by magnetohydrodynamically driven outflows from an accretion disk. We demonstrate that, within only 50ms postmerger, ≳2×10^{-2}M_{⊙} of lanthanide-free, quasispherical ejecta with velocities ∼0.1-0.2c is launched, yielding a kilonova signal consistent with GW170817 on timescales of ≲5 d.

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼1014.5 g cm−3, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

Bound Debris Expulsion from Neutron Star Merger Remnants

Many studies have found that neutron star mergers leave a fraction of the stars’ mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site of r-process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars’ interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star.

Role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers

Role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers

Magnetic field effects on nucleosynthesis and kilonovae from neutron star merger remnants

ABSTRACT We investigate the influence of parametric magnetic field configurations of a hypermassive neutron star (HMNS) on the outflow properties, nucleosynthesis yields, and kilonova light curves. We perform three-dimensional dynamical space–time general-relativistic magnetohydrodynamic simulations, including a neutrino leakage scheme, microphysical finite-temperature equation of state, and an initial poloidal magnetic field. We find that varying the magnetic field strength and falloff impacts the formation of magnetized winds or mildly relativistic jetted outflows, which in turn has profound effects on the outflow properties. All of the evolved configurations collapse to a black hole ∼38–40 ms after coalescence, where the ones forming jetted outflows seem more effective at redistributing angular momentum, which result in earlier collapse times. Larger mass ejecta rates and radial velocities of unbound material characterize the systems that form jetted outflows. The bolometric light curves of the kilonovae and r-process yields that are produced by the post-merger remnant system change considerably with different magnetic field parameters. We conclude that the magnetic field strength and falloff have robust effects on the outflow properties and electromagnetic observables. This can be particularly important as the total ejecta mass from our simulations (≃10−3 M⊙) makes the ejecta from HMNS a compelling source to power kilonova through radioactive decay of r-process elements.

Local magneto-shear instability in Newtonian gravity

ABSTRACT The magneto-rotational instability (MRI), which is due to an interplay between a sheared background and the magnetic field, is commonly considered a key ingredient for developing and sustaining turbulence in the outer envelope of binary neutron star merger remnants. To assess whether (or not) the instability is active and resolved, criteria originally derived in the accretion disc literature, thus exploiting the symmetries of such systems, are often used. In this paper, we discuss the magneto-shear instability as a truly local phenomenon, relaxing common symmetry assumptions on the background on top of which the instability grows. This makes the discussion well suited for highly dynamical environments such as binary mergers. We find that, although this is somewhat hidden in the usual derivation of the MRI dispersion relation, the instability crucially depends on the assumed symmetries. Relaxing the symmetry assumptions in the background, we find that the role of the magnetic field is significantly diminished, as it affects the modes’ growth but does not drive it. We conclude by making contact with a suitable filtering operation, as this is key to separating background and fluctuations in highly dynamical systems.

Flares from Merged Magnetars: Their Prospects as a New Population of Gamma-Ray Counterparts of Binary Neutron Star Mergers

Long-lived massive magnetars are expected to be remnants of some binary neutron star (BNS) mergers. In this paper, we argue that the magnetic powered flaring activities of these merged magnetars would occur dominantly in their early millisecond-period-spin phase, which is in the timescale of days. Such flares endure significant absorption by the ejecta from the BNS collision, and their detectable energy range is from 0.1 to 10 MeV, in a time lag of approximately days after the merger events indicated by the gravitational wave chirps. We estimate the rate of such flares in different energy ranges, and find that there could have been 0.1–10 cases detected by Fermi/GBM. A careful search for ∼10 ms spin-period modulation in weak short gamma-ray bursts (GRBs) may identify them from the archival data. The next-generation megaelectronvolt detectors could detect them at a mildly higher rate. The recent report on the Quasi-Period-Oscillation found in two BASTE GRBs should not be considered as cases of such flares, for they were detected in a lower energy range and with a much shorter period spin modulation.

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.