Neutron Star Remnant Research Articles

Features and prospects for kilonova remnant detection with current and future surveys

We study the observable spectral and temporal properties of kilonova remnants (KNRs) analytically, and point out quantitative differences with respect to supernova remnants. We provide detection prospects of KNRs in the context of ongoing radio surveys. We find that there is a good chance to expect tens of these objects in future surveys with a flux threshold of $ 0.1$ mJy. Kilonova remnants from a postulated population of long-lived supermassive neutron star remnants of neutron star mergers are even more likely to be detected, as they are extremely bright and peak earlier. For an ongoing survey with a threshold of sim 1 mJy, we expect to find tens to hundreds of such objects if they are a significant fraction of the total kilonova (KN) population. Considering that there are no such promising KN candidates in presently ongoing surveys, we constrain the fraction of these extreme KN to be no more than 30 percent of the overall KN population. This constraint depends sensitively on the details of ejecta mass and external density distribution.

TIC 290061484: A Triply Eclipsing Triple System with the Shortest Known Outer Period of 24.5 Days

We have discovered a triply eclipsing triple-star system, TIC 290061484, with the shortest known outer period, P out, of only 24.5 days. This “eclipses” the previous record set by λ Tauri at 33.02 days, which held for 68 yr. The inner binary, with an orbital period of P in = 1.8 days, produces primary and secondary eclipses and exhibits prominent eclipse timing variations with the same periodicity as the outer orbit. The tertiary star eclipses, and is eclipsed by, the inner binary with pronounced asymmetric profiles. The inclinations of both orbits evolve on observable timescales such that the third-body eclipses exhibit dramatic depth variations in TESS data. A photodynamical model provides a complete solution for all orbital and physical parameters of the triple system, showing that the three stars have masses of 6.85, 6.11, and 7.90 M ⊙, radii near those corresponding to the main sequence, and T eff in the range of 21,000–23,700 K. Remarkably, the model shows that the triple is in fact a subsystem of a hierarchical 2+1+1 quadruple with a distant fourth star. The outermost star has a period of ∼3200 days and a mass comparable to the stars in the inner triple. In ∼20 Myr, all three components of the triple subsystem will merge, undergo a Type II supernova explosion, and leave a single remnant neutron star. At the time of writing, TIC 290061484 is the most compact triple system and one of the tighter known compact triples (i.e., P out/P in = 13.7).

The Jittering Jets Explosion Mechanism in Electron Capture Supernovae

We conduct one-dimensional stellar-evolution simulations of stars with zero-age main-sequence (ZAMS) masses of M ZAMS = 8.8 − 9.45 M ⊙ toward core collapse by electron capture and find that the convective zone of the precollapse core can supply the required stochastic angular momentum fluctuations to set a jet-driven electron capture supernova explosion in the frame of the jittering jets explosion mechanism. By our assumed criteria of a minimum convective specific angular momentum and an accreted mass during jet launching of M acc ≃ 0.001−0.01 M ⊙, the layer in the convective zone that when accreted launches the exploding jittering jets resides in the helium-rich zone. Depending on the model, this exploding layer is accreted at about a minute to a few hours after core collapse occurs, much shorter than the time the exploding shock crosses the star. The final (gravitational) mass of the neutron star (NS) remnant is in the range of M NS = 1.25−1.43 M ⊙.

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.

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.

Morphology-independent characterization method of postmerger gravitational wave emission from binary neutron star coalescences

Gravitational waves (GWs) emitted during the coalescence of binary neutron star (BNS) systems carry information about the equation of state (EoS) describing the extremely dense matter inside neutron stars (NSs). In particular, the EoS determines the fate of the binary after the merger: a prompt collapse to black hole (BH), or the formation of a NS remnant that is either stable or survives up to a few seconds before collapsing to a BH. Determining the evolution of a BNS system will therefore place strong constraints on the EoS. We present a morphology-independent method, developed in the framework of the coherentWaveBurst analysis of signals from ground-based interferometric detectors of GWs. The method characterizes the time-frequency postmerger GW emission from a BNS system, and determines whether, after the merger, it formed a remnant NS or promptly collapsed to a BH. We measure the following quantities to characterize the postmerger emission: ratio of signal energies and match of luminosity profile in different frequency bands, weighted central frequency and bandwidth. From these quantities, based on the study of signals simulated through injections of numerical relativity waveforms, we build a statistics to discriminate between the different scenarios after the merger. Finally, we test our method on a set of signals simulated with new models, to estimate its efficiency as a function of the source distance.

Gravitational Wave Eigenfrequencies from Neutrino-driven Core-collapse Supernovae

Core-collapse supernovae (CCSNe) are predicted to produce gravitational waves (GWs) that may be detectable by Advanced LIGO/Virgo. These GW signals carry information from the heart of these cataclysmic events, where matter reaches nuclear densities. Recent studies have shown that it may be possible to infer the properties of the proto-neutron star (PNS) via GWs generated by hydrodynamic perturbations of the PNS. However, we lack a comprehensive understanding of how these relationships may change with the properties of CCSNe. In this work, we build a self-consistent suite of over 1000 exploding CCSNe from a grid of progenitor masses and metallicities combined with six different nuclear equations of state (EOS). Performing a linear perturbation analysis on each model, we compute the resonant GW frequencies of the PNS, and we motivate a time-agnostic method for identifying characteristic frequencies of the dominant GW emission. From this, we identify two characteristic frequencies, of the early- and late-time signal, that measure the surface gravity of the cold remnant neutron star, and simultaneously constrain the hot nuclear EOS. However, we find that the details of the CCSN model, such as the treatment of gravity or the neutrino transport, and whether it explodes, noticeably change the magnitude and evolution of the PNS eigenfrequencies.

Kilonovae of binary neutron star mergers leading to short-lived remnant neutron star formation

ABSTRACT We study kilonova emission from binary neutron star (BNS) mergers for the case that a remnant massive neutron star (MNS) forms and collapses to a black hole within 20 ms after the onset of the merger (which we refer to as ‘a short-lived case’) by consistently employing numerical relativity and nucleosynthesis results. We find that such kilonovae are fainter and last shorter than those for BNSs resulting in the formation of long-lived (${\gg} 1\, {\rm s}$) MNSs, in particular in the optical band. The resulting light curves are too faint and last for a too short duration to explain the kilonova observation for the BNS associated with GW170817, indicating that the merger remnant formed in GW170817 is unlikely to have collapsed to a black hole within a short period of time (∼20 ms) after the onset of the merger. Our present result implies that early observation is necessary to detect kilonovae associated with BNSs leading to short-lived MNS formation in particular for the optical blue band as well as that kilonovae could be hidden by the gamma-ray burst afterglow for nearly face-on observation. We provide a possible approximate scaling law for near-infrared light curves with the given reference time and magnitude when the decline power of the z-band magnitude, dMz/dlog10t, reaches 2.5. This scaling law suggests that the HK-band follow-up observation should be at least 1 mag deeper than that for the z-band reference magnitude and earlier than 4 times the reference time.

The Neutron Star to Black Hole Mass Gap in the Frame of the Jittering Jets Explosion Mechanism (JJEM)

I build a toy model in the frame of the jittering jets explosion mechanism (JJEM) of core collapse supernovae that incorporates both the stochastically varying angular momentum component of the material that the newly born neutron star (NS) accretes and the constant angular momentum component, and show that the JJEM can account for the ≃2.5–5M ⊙ mass gap between NSs and black holes (BHs). The random component of the angular momentum results from pre-collapse core convection fluctuations that are amplified by post-collapse instabilities. The fixed angular momentum component results from pre-collapse core rotation. For slowly rotating pre-collapse cores the stochastic angular momentum fluctuations form intermittent accretion disks (or belts) around the NS with varying angular momentum axes in all directions. The intermittent accretion disk/belt launches jets in all directions that expel the core material in all directions early on, hence leaving an NS remnant. Rapidly rotating pre-collapse cores form an accretion disk with angular momentum axis that is about the same as the pre-collapse core rotation. The NS launches jets along this axis and hence the jets avoid the equatorial plane region. Inflowing core material continues to feed the central object from the equatorial plane increasing the NS mass to form a BH. The narrow transition from slow to rapid pre-collapse core rotation, i.e., from an efficient to inefficient jet feedback mechanism, accounts for the sparsely populated mass gap.

Constraining the ellipticity and frequency of binary neutron star remnant via its gravitational-wave and electromagnetic radiations

ABSTRACT The nature of the merger remnant of binary neutron star remains an open question. From the theoretical point of view, one possible outcome is a supra-massive neutron star (SMNS), which is supported by rigid rotation and through its survival of hundreds of seconds before collapsing into a black hole. If this is the case, the SMNS can emit continuous gravitational waves (GW) and electromagnetic radiation, particularly in the X-ray band. In this work, the ellipticity and initial frequency of SMNS are constrained with a Bayesian framework using simulated X-ray and GW signals, which could be detected by The Transient High Energy Sky and Early Universe Surveyor and Einstein Telescope, respectively. We found that only considering the X-ray emission cannot completely constrain the initial frequency and ellipticity of the SMNS, but it can reduce the ranges of the parameters. Afterwards, we can use the posterior distribution of the X-ray parameter estimates as a prior for the GW parameter estimates. It was found that the 95 per cent credible region of the joint X-ray–GW analysis was about 105 times smaller than that of the X-ray analysis alone.

Plausible presence of new state in neutron stars with masses above [formula omitted

We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030 + 0451, PSR J0740 + 6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory (χEFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forward neural network model embedded in the Bayesian nonparametric inference, the structure of NS matter’s sound speed cs is explored in a model-agnostic way. It is found that a peak structure is common in the cs2 posterior, locating at (2.4-4.8)ρsat (nuclear saturation density) and cs2 exceeds c2/3 at 90% credibility. The non-monotonic behavior suggests evidence of the state deviating from the hadronic matter inside the very massive NSs. Assuming the new/exotic state is featured as it is softer than typical hadronic models or even with hyperons, we find that a sizable (⩾10-3M⊙) exotic core, likely made of quark matter, is plausible for the NS with a gravitational mass above about 0.98MTOV, where MTOV represents the maximum gravitational mass of a non-rotating cold NS. The inferred MTOV=(2.18-0.13+0.27)M⊙ (90% credibility) is well consistent with the value of (2.17-0.12+0.15)M⊙ estimated independently with GW170817/GRB 170817A/AT2017gfo assuming a temporary supramassive NS remnant formed after the merger. PSR J0740 + 6620, the most massive NS detected so far, may host a sizable exotic core with a probability of ≈0.36.

The implications of large binding energies of massive stripped core collapse supernova progenitors on the explosion mechanism

ABSTRACT We examine the binding energies of massive stripped-envelope core collapse supernova (SECCSN) progenitors with the stellar evolution code mesa, and find that the jittering jets explosion mechanism is preferred for explosions where carbon–oxygen cores with masses of ${\gtrsim} 20 \, \mathrm{M}_\odot$ collapse to leave a neutron star (NS) remnant. We calculate the binding energy at core collapse under the assumption that the remnant is an NS. Namely, stellar gas above mass coordinate of ${\simeq} 1.5\text{{--}}2.5 \, \mathrm{M}_\odot$ is ejected in the explosion. We find that the typical binding energy of the ejecta of stripped-envelope (SE) progenitors with carbon–oxygen core masses of $M_{\rm CO} \gtrsim 20 \, \mathrm{M}_\odot$ is $E_{\rm bind} \gtrsim 2 \times 10^{51} {~\rm erg}$. We claim that jets are most likely to explode such cores as jet-driven explosion mechanisms can supply high energies to the explosion. We apply our results to SN 2020qlb, which is an SECCSN with a claimed core mass of ${\simeq} 30\!-\!50 \, \mathrm{M}_\odot$, and conclude that the jittering jets explosion mechanism best accounts for such an explosion that leaves an NS.

Central engine of GRB170817A: Neutron star versus Kerr black hole based on multimessenger calorimetry and event timing

Context.LIGO–Virgo–KAGRA observations may identify the remnant of compact binary coalescence and core-collapse supernovae associated with gamma-ray bursts. The multimessenger event GW170817–GRB170817A appears ripe for this purpose thanks to its fortuitous close proximity at 40 Mpc. Its post-merger emission, ℰGW, in a descending chirp can potentially break the degeneracy in spin-down of a neutron star or black hole remnant by the relatively large energy reservoir in the angular momentum,EJ, of the latter according to the Kerr metric.Aims.The complex merger sequence of GW170817 is probed for the central engine of GRB170817A by multimessenger calorimetry and event timing.Methods.We used model-agnostic spectrograms with equal sensitivity to ascending and descending chirps generated by time-symmetric butterfly matched filtering. The sensitivity was calibrated by response curves generated by software injection experiments, covering a broad range in energies and timescales. The statistical significance for candidate emission from the central engine of GRB170817A is expressed by probabilities of false alarm (PFA; type I errors) derived from an event-timing analysis. Probability density functions (PDF) were derived for start-timets, identified via high-resolution image analyses of the available spectrograms. For merged (H1,L1)-spectrograms of the LIGO detectors, a PFAp1derives from causality intsgiven GW170817–GRB17081A (contextual). A statistically independent confirmation is presented in individual H1 and L1 analyses, quantified by a second PFAp2of consistency in their respective observations ofts(acontextual). A combined PFA derives from their product since the mean and (respectively) the difference in timing are statistically independent.Results.Applied to GW170817–GRB170817A, PFAs of event timing intsproducep1 = 8.3 × 10−4andp2 = 4.9 × 10−5of a post-merger output ℰGW≃ 3.5%M⊙c2(p1p2 = 4.1 × 10−8, equivalentZ-score 5.48). ℰGWexceedsEJof the hyper-massive neutron star in the immediate aftermath of GW170817, yet it is consistent withEJrejuvenated in gravitational collapse to a Kerr black hole. Similar emission may be expected from energetic core-collapse supernovae producing black holes of interest to upcoming observational runs by LIGO–Virgo–KAGRA.

Magnetar Wind-Driven Shock Breakout Emission after Double Neutron Star Mergers: The Effect of the Anisotropy of the Merger Ejecta

A rapidly rotating and highly magnetized remnant neutron star (NS; magnetar) could survive from a merger of double NSs and drive a powerful relativistic wind. The early interaction of this wind with the previous merger ejecta can lead to shock breakout (SBO) emission mainly in ultraviolet and soft X-ray bands, which provides an observational signature for the existence of the remnant magnetar. Here, we investigate the effect of an anisotropic structure of the merger ejecta on the SBO emission. It is found that the bolometric light curve of the SBO emission can be broadened, since the SBO can occur at different times for different directions. In more detail, the profile of the SBO light curve can be highly dependent on the ejecta structure and, thus, we can in principle use the SBO light curves to probe the structure of the merger ejecta in future.

R-process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants

Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS.

On the diversity of magnetar-driven kilonovae

ABSTRACT A non-negligible fraction of binary neutron star mergers are expected to form long-lived neutron star remnants, dramatically altering the multimessenger signatures of a merger. Here, we extend existing models for magnetar-driven kilonovae and explore the diversity of kilonovae and kilonova afterglows. Focusing on the role of the (uncertain) magnetic field strength, we study the resulting electromagnetic signatures as a function of the external dipolar and internal toroidal fields. These two parameters govern, respectively, the competition between magnetic-dipole spin-down and gravitational-wave spin-down (due to magnetic-field deformation) of the rapidly rotating remnant. We find that even in the parameter space where gravitational-wave emission is dominant, a kilonova with a magnetar central engine will be significantly brighter than one without an engine, as this parameter space is where more of the spin-down luminosity is thermalized. In contrast, a system with minimal gravitational-wave emission will produce a kilonova that may be difficult to distinguish from ordinary kilonovae unless early epoch observations are available. However, as the bulk of the energy in this parameter space goes into accelerating the ejecta, such a system will produce a brighter kilonova afterglow that will peak in shorter times. To effectively hide the presence of the magnetar from the kilonova and kilonova afterglow, the rotational energy inputted into the ejecta must be ≲10−3to 10−2Erot. We discuss the different diagnostics available to identify magnetar-driven kilonovae in serendipitous observations and draw parallels to other potential magnetar-driven explosions, such as superluminous supernovae and broad-line supernovae Ic.

Correction to: Models of binary neutron star remnants with tabulated equations of state

Correction to: Models of binary neutron star remnants with tabulated equations of state

Finite-temperature effects in dynamical spacetime binary neutron star merger simulations: validation of the parametric approach

ABSTRACT Parametric equations of state (EoSs) provide an important tool for systematically studying EoS effects in neutron star merger simulations. In this work, we perform a numerical validation of the M*-framework for parametrically calculating finite-temperature EoS tables. The framework, introduced by Raithel et al., provides a model for generically extending any cold, β-equilibrium EoS to finite temperatures and arbitrary electron fractions. In this work, we perform numerical evolutions of a binary neutron star merger with the SFHo finite-temperature EoS, as well as with the M*-approximation of this same EoS, where the approximation uses the zero-temperature, β-equilibrium slice of SFHo and replaces the finite-temperature and composition-dependent parts with the M*-model. We find that the approximate version of the EoS is able to accurately recreate the temperature and thermal pressure profiles of the binary neutron star remnant, when compared to the results found using the full version of SFHo. We additionally find that the merger dynamics and gravitational wave signals agree well between both cases, with differences of $\lesssim 1\!-\!2\,{\textrm{per cent}}$ introduced into the post-merger gravitational wave peak frequencies by the approximations of the EoS. We conclude that the M*-framework can be reliably used to probe neutron star merger properties in numerical simulations.