Neutron Star Remnant Research Articles

Numerical Relativity Simulations of the Neutron Star Merger GW170817: Long-term Remnant Evolutions, Winds, Remnant Disks, and Nucleosynthesis

Abstract We present a systematic numerical relativity study of the dynamical ejecta, winds, and nucleosynthesis in neutron star (NS) merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOSs) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamic turbulence up to 100 ms postmerger. Spiral density waves propagating from the NS remnant to the disk trigger a wind with mass flux ∼0.1–0.5 M ⊙ s−1, which persists for the entire simulation as long as the remnant does not collapse to a black hole. This wind has average electron fraction ≳0.3 and average velocity ∼0.1–0.17 c and thus is a site for the production of weak r-process elements (mass number A < 195). Disks around long-lived remnants have masses ∼0.1–0.2 M ⊙, temperatures peaking at ≲10 MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multicomponent kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the long-lived mergers of comparable mass robustly accounts for all the r-process peaks, from mass number ∼75 to actinides in terms of solar abundances. Total abundances are weakly dependent on the EOS, while the mass ratio affects the production of first-peak elements.

On rare core collapse supernovae inside planetary nebulae

ABSTRACT We conduct simulations using mesa of the reverse formation of a white dwarf (WD)–neutron star (NS) binary system in which the WD forms before the NS. We conclude that a core collapse supernova (CCSN) explosion might occur inside a planetary nebula (PN) only if a third star forms the PN. In this WD–NS reverse binary evolution, the primary star evolves and transfers mass to the secondary star, forms a PN, and leaves a WD remnant. If the mass-transfer brings the secondary star to have a mass of $\gtrsim 8\, \mathrm{ M}_\odot$ before it develops a helium core, and if the secondary does not suffer an enhanced mass-loss before it develops a massive helium core, e.g. by mass-transfer, it explodes as a CCSN and leaves an NS remnant. The time period from the formation of the PN by the primary to the explosion of the secondary is $\gtrsim 10^6 {~\rm yr}$. By that time, the PN has long dispersed into the interstellar medium. In a binary system with nearly equal-mass components, the first mass-transfer episode takes place after the secondary star has developed a helium core and it ends its life forming a PN and a WD. The formation of a CCSN inside a PN (CCSNIP) requires the presence of a third star. The third star should be less massive than the secondary star but by no more than few ×0.01 M⊙. We estimate that the rate of CCSNIP is ≈10−4 times the rate of all CCSNe.

Magnetically Driven Baryon Winds from Binary Neutron Star Merger Remnants and the Blue Kilonova of 2017 August

Abstract The observation of a radioactively powered kilonova associated with the first binary neutron star (BNS) merger detected in gravitational waves proved that these events are ideal sites for the production of heavy r-process elements. However, the physical origin of the ejected material responsible for the early (“blue”) and late (“red”) components of this kilonova is still debated. Here, we investigate the possibility that the early/blue kilonova originated from the magnetically driven baryon wind launched after merger by the metastable neutron star remnant. Exploiting a magnetized BNS merger simulation with over 250 ms of post-merger evolution, we can follow for the first time the full mass-ejection process up to its final decline. We find that the baryon wind carries ≃0.010–0.028 M ⊙ of unbound material, proving that the high mass estimated for the blue kilonova can be achieved. We also find expansion velocities of up to ∼0.2c, consistent with the lower end of the observational estimates, and we discuss possible effects neglected here that could further increase the final ejecta velocity. Overall, our results show that the magnetically driven baryon wind represents a viable channel to explain the blue kilonova.

Intrinsic Properties of the Engine and Jet that Powered the Short Gamma-Ray Burst Associated with GW170817

Abstract GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15°–30° away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet–wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant.

The Pre-explosion Mass Distribution of Hydrogen-poor Superluminous Supernova Progenitors and New Evidence for a Mass–Spin Correlation

Abstract Despite indications that superluminous supernovae (SLSNe) originate from massive progenitors, the lack of a uniformly analyzed statistical sample has so far prevented a detailed view of the progenitor mass distribution. Here we present and analyze the pre-explosion mass distribution of hydrogen-poor SLSN progenitors as determined from uniformly modeled light curves of 62 events. We construct the distribution by summing the ejecta mass posteriors of each event, using magnetar light-curve models presented in our previous works (and using a nominal neutron star remnant mass). The resulting distribution spans 3.6–40 M ⊙, with a sharp decline at lower masses, and is best fit by a broken power law described by at 3.6–8.6 M ⊙ and at 8.6–40 M ⊙. We find that observational selection effects cannot account for the shape of the distribution. Relative to Type Ib/c SNe, the SLSN mass distribution extends to much larger masses and has a different power-law shape, likely indicating that the formation of a magnetar allows more massive stars to explode as some of the rotational energy accelerates the ejecta. Comparing the SLSN distribution with predictions from single and binary star evolution models, we find that binary models for a metallicity of Z ≲ 1/3 Z ⊙ are best able to reproduce its broad shape, in agreement with the preference of SLSNe for low metallicity environments. Finally, we uncover a correlation between the pre-explosion mass and the magnetar initial spin period, where SLSNe with low masses have slower spins, a trend broadly consistent with the effects of angular momentum transport evident in models of rapidly rotating carbon–oxygen stars.

Late X-ray flares from the interaction of a reverse shock with a stratified ejecta in GRB afterglows: simulations on a moving mesh

ABSTRACT Late activity of the central engine is often invoked in order to explain the flares observed in the early X-ray afterglow of gamma-ray bursts, either in the form of an active neutron star remnant or (fall-back) accretion on to a black hole. However, these scenarios are not always plausible, in particular when flares are delayed to very late times after the burst. Recently, a new scenario was proposed that suggests X-ray flares can be the result of the passing of a long-lived reverse shock through a stratified ejecta, with the advantage that it does not require late-time engine activity. In this work, we numerically demonstrate this scenario to be physically plausible, by performing one-dimensional simulations of ejecta dynamics and emission using our novel moving-mesh relativistic hydrodynamics code. Improved efficiency and precision over previous work enables the exploration of a broader range of set-ups. We can introduce a more physically realistic description of the circumburst medium mass density. We can also locally trace the cooling of electrons when computing the broad-band emission from these set-ups. We show that the synchrotron cooling time-scale can dominate the flare decay time if the stratification in the ejecta is constrained to a localized angular region inside the jet, with size corresponding to the relativistic causal connection angle, and that it corresponds to values reported in observations. We demonstrate that this scenario can produce a large range of observed flare times, suggesting a connection between flares and initial ejection dynamics rather than with late-time remnant activity.

Low-energy core-collapse supernovae in the frame of the jittering jets explosion mechanism

ABSTRACT We relate the pre-explosion binding energy of the ejecta of core-collapse supernovae (CCSNe) of stars with masses in the lower range of CCSNe and the location of the convection zones in the pre-collapse core of these stars, to explosion properties in the frame of the jittering jets explosion mechanism. Our main conclusion is that in the frame of the jittering jets explosion mechanism the remnant of a pulsar in these low-energy CCSNe has some significance, in that the launching of jets by the newly born neutron star (NS) spins-up the NS and create a pulsar. We crudely estimated the period of the pulsars to be tens of milliseconds in these cases. The convective zones seed perturbations that lead to accretion of stochastic angular momentum that in turn is assumed to launch jittering jets in this explosion mechanism. We calculate the binding energy and the location of the convective zones with the stellar evolution code mesa. For the lowest stellar masses, we study, MZAMS ≃ 8.5–11 M⊙, the binding energy above the convective zones is low, and so is the expected explosion energy in the jittering jets explosion mechanism that works in a negative feedback cycle. The expected mass of the NS remnant is MNS ≈ 1.25–1.6 M⊙, even for these low-energy CCSNe.

Collimated outflows from long-lived binary neutron star merger remnants

ABSTRACT The connection between short gamma-ray bursts (SGRBs) and binary neutron star (BNS) mergers was recently confirmed by the association of GRB 170817A with the merger event GW170817. However, no conclusive indications were obtained on whether the merger remnant that powered the SGRB jet was an accreting black hole (BH) or a long-lived massive neutron star (NS). Here, we explore the latter case via BNS merger simulations covering up to 250 ms after merger. We report, for the first time in a full merger simulation, the formation of a magnetically driven collimated outflow along the spin axis of the NS remnant. For the system at hand, the properties of such an outflow are found largely incompatible with an SGRB jet. With due consideration of the limitations and caveats of our present investigation, our results favour a BH origin for GRB 170817A and SGRBs in general. Even though this conclusion needs to be confirmed by exploring a larger variety of physical conditions, we briefly discuss possible consequences of all SGRB jets being powered by accreting BHs.

Magnetohydrodynamic simulations of binary neutron star mergers in general relativity: Effects of magnetic field orientation on jet launching.

Binary neutron star mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission across the spectrum. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of highly spinning binary neutrons stars initially on quasicircular orbits that merge and undergo delayed collapse to a black hole. The binaries consist of two identical stars modeled as Γ = 2 polytropes with spin χ NS = 0.36 aligned along the direction of the total orbital angular momentum L. Each star is initially threaded by a dynamical unimportant interior dipole magnetic field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment μ is either aligned or perpendicular to L and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where μ in both stars is aligned along L. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta (which can give rise to a detectable kilonova) are very sensitive to the magnetic field orientation. By contrast, the physical properties of the black hole + disk remnant, such as the mass and spin of the black hole, the accretion rate, and the electromagnetic (Poynting) luminosity, are roughly independent of the initial magnetic field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results suggest that, as in the black hole-neutron star merger scenario, an incipient jet emerges only when the seed magnetic field has a sufficiently large-scale poloidal component aligned to the initial orbital angular momentum. The lifetime [Δt ≳ 140(M NS/1.625 M ⊙) ms] and Poynting luminosities [L EM ≃ 1052 erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford-Znajek mechanism for launching jets.

Effects of dark matter on the upper bound mass of neutron stars

Abstract Observations have indicated that we do not see neutron stars (NS) of mass near the theoretical upper limit as predicted. Here we invoke the role of dark matter (DM) particles in star formation, and their role in lowering the mass of remnants eventually formed from these stars. Massive stars can capture DM particles more effectively than the lower mass stars, thus further softening the equation of state of the remnant neutron stars. We also look at the capture of DM particles by the NS, which could further soften the upper mass limit of NS. The admixture of DM particles would be higher at earlier epochs (high z).

Constraining the Long-lived Magnetar Remnants in Short Gamma-Ray Bursts from Late-time Radio Observations

Abstract The joint detection of GW170817 and GRB 170817A indicated that at least a fraction of short gamma-ray bursts (SGRBs) originate from binary neutron star (BNS) mergers. One possible remnant of a BNS merger is a rapidly rotating, strongly magnetized neutron star, which has been discussed as one possible central engine for gamma-ray bursts. For a rapidly rotating magnetar central engine, the deposition of the rotation energy into the ejecta launched from the merger could lead to bright radio emission. The brightness of radio emission years after an SGRB would provide an estimate of the kinetic energy of ejecta and, hence, a possible constraint on the BNS merger product. We perform a more detailed calculation on the brightness of radio emission from the interaction between the merger ejecta and circumburst medium in the magnetar scenario, invoking several important physical processes such as generic hydrodynamics, relativistic effects, and the deep Newtonian phase. We use the model to constrain the allowed parameter space for 15 SGRBs that have late radio observations. Our results show that an injection energy of E inj ∼ 1052 erg is allowed for all the cases, which suggests that the possibility of a supramassive or hypermassive neutron star remnant is not disfavored by the available radio data.

Constraining the gravitational-wave afterglow from a binary neutron star coalescence

ABSTRACT Binary neutron star mergers are rich laboratories for physics, accessible with ground-based interferometric gravitational-wave detectors such as the Advanced LIGO and Advanced Virgo. If a neutron star remnant survives the merger, it can emit gravitational waves that might be detectable with the current or next generation detectors. The physics of the long-lived post-merger phase is not well understood and makes modelling difficult. In particular the phase of the gravitational-wave signal is not well modelled. In this paper, we explore methods for using long duration post-merger gravitational-wave signals to constrain the parameters and the properties of the remnant. We develop a phase-agnostic likelihood model that uses only the spectral content for parameter estimation and demonstrate the calculation of a Bayesian upper limit in the absence of a signal. With the millisecond magnetar model, we show that for an event like GW170817, the ellipticity of a long-lived remnant can be constrained to less than about 0.5 in the parameter space used.

Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

Abstract Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2–4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a Neutron Star Extreme Matter Observatory (NEMO): a gravitational-wave interferometer optimised to study nuclear physics with merging neutron stars. The concept uses high-circulating laser power, quantum squeezing, and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above 1 kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year and potentially allow for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

A Deep Targeted Search for Fast Radio Bursts from the Sites of Low-redshift Short Gamma-Ray Bursts

Abstract Some short gamma-ray bursts (SGRBs) are thought to be caused by the mergers of binary neutron stars which may sometimes produce massive neutron star remnants capable of producing extragalactic fast radio bursts (FRBs). We conducted a deep search for FRBs from the sites of six low-redshift SGRBs. We collected high time- and frequency-resolution data from each of the sites for 10 hr using the 2 GHz receiver of the Green Bank Telescope (GBT). Two of the SGRB sites we targeted were visible with the Arecibo Radio Telescope with which we conducted an additional 10 hr of 1.4 GHz observations for each. We searched our data for FRBs using the GPU-optimized dedispersion algorithm heimdall and the machine-learning-based package Fast Extragalactic Transient Candidate Hunter. We did not discover any FRBs, but would have detected any with peak flux densities in excess of 87 mJy at the GBT or 21 mJy at Arecibo with a signal-to-noise ratio of at least 10. The isotropic-equivalent energy of any FRBs emitted from these sites in our bands during our observations must not have exceeded a few times 1038 erg, comparable to some of the lowest energy bursts yet seen from the first known repeating FRB 121102.

J0453+1559: A Neutron Star–White Dwarf Binary from a Thermonuclear Electron-capture Supernova?

Abstract The compact binary radio pulsar system J0453+1559 consists of a recycled pulsar as primary component of 1.559(5) M ⊙ and an unseen companion star of 1.174(4) M ⊙. Because of the relatively large orbital eccentricity of e = 0.1125, it was argued that the companion is a neutron star (NS), making it the NS with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictory concerning the possibility of producing such a low-mass NS remnant. Here we challenge the NS interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen–neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several 0.1 M ⊙ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR J0453+1559 by a NS–white dwarf system. More work on tECSNe is needed to assess the viability of this scenario.

Kilohertz gravitational waves from binary neutron star remnants: Time-domain model and constraints on extreme matter

The remnant star of a neutron star merger is an anticipated loud source of kiloHertz gravitational waves that conveys unique information on the equation of state of hot matter at extreme densities. Observations of such signals are hampered by the photon shot noise of ground-based interferometers and pose a challenge for gravitational-wave astronomy. We develop an analytical time-domain waveform model for postmerger signals informed by numerical relativity simulations. The model completes effective-one-body waveforms for quasi-circular nonspinning binaries in the kiloHertz regime. We show that a template-based analysis can detect postmerger signals with a minimal signal-to-noise ratios (SNR) of 8, corresponding to GW170817-like events for third-generation interferometers. Using Bayesian model selection and the complete inspiral-merger-postmerger waveform model it is possible to infer whether the merger outcome is a prompt collapse to a black hole or a remnant star. In the latter case, the radius of the maximum mass (most compact) nonrotating neutron star can be determined to kilometer precision. We demonstrate the feasibility of inferring the stiffness of the equation of state at extreme densities using the quasiuniversal relations deduced from numerical-relativity simulations.

Searching for Hypermassive Neutron Stars with Short Gamma-Ray Bursts

Abstract Neutron star mergers can form a hypermassive neutron star (HMNS) remnant, which may be the engine of a short gamma-ray burst (SGRB) before it collapses to a black hole, possibly several hundred milliseconds after the merger. During the lifetime of an HMNS, numerical relativity simulations indicate that it will undergo strong oscillations and emit gravitational waves with frequencies of a few kilohertz, which are unfortunately too high for detection to be probable with the Advanced Laser Interferometer Gravitational-Wave Observatory. Here we discuss the current and future prospects for detecting these frequencies as modulation of the SGRB. The understanding of the physical mechanism responsible for the HMNS oscillations will provide information on the equation of state of the hot HMNS, and the observation of these frequencies in the SGRB data would give us insight into the emission mechanism of the SGRB.

No Pulsar Wind Nebula in the Southern Blowout Region of the Cygnus Loop

Abstract We report on optical observations of the Katsuda et al. candidate X-ray pulsar and pulsar wind nebula in the Cygnus Loop supernova remnant. We determine that the point source suggested to be a pulsar is actually the nucleus of a Seyfert 1 galaxy at redshift z = 0.2080, while the diffuse X-ray source, which is displaced by 2.′6 from the point source, is a cluster of galaxies at z = 0.223. We also analyze an archival follow-up XMM-Newton observation of this field, the results of which support our extragalactic identifications. Thus, a long expected neutron star remnant of the Cygnus Loop explosion remains elusive.

Thermodynamics conditions of matter in neutron star mergers

Matter in neutron star collisions can reach densities up to few times the nuclear saturation threshold and temperatures up to one hundred MeV. Understanding the structure and composition of such matter requires many-body nonperturbative calculations that are currently highly uncertain.Unique constraints on the neutron star matter are provided by gravitational-wave observations aided by numerical relativity simulations. In this work, we explore the thermodynamical conditions of matter and radiation along the merger dynamics. We consider 3 microphysical equation of state models and numerical relativity simulations including an approximate neutrino transport scheme. The neutron star cores collision and their multiple centrifugal bounces heat the initially cold matter to several tens of MeV. Streams of hot matter with initial densities $\sim1-2\rho_0$ move outwards and cool due to decompression and neutrino emission. The merger can result in a neutron star remnant with densities up to $3-5\rho_0$ and temperatures $\sim 50$~MeV. The highest temperatures are confined in an approximately spherical annulus at densities $\sim\rho_0$. Such temperatures favour positron-neutron capture at densities $\sim\rho_0$, thus leading to a neutrino emission dominated by electron antineutrinos. We study the impact of trapped neutrinos on the remnant matter's pressure, electron fraction and temperature and find that it has a negligible effect. Disks around neutron star or black hole remnant are neutron rich and not isentropic, but they differ in size, entropy and lepton fraction depending on the nature of the central object. In presence of a black hole, disks are smaller and mostly transparent to neutrinos; in presence of a massive neutron star, they are more massive, geometrically and optically thick.

Constraint on the maximum mass of neutron stars using GW170817 event

We revisit the constraint on the maximum mass of cold spherical neutron stars coming from the observational results of GW170817. We develop a new framework for the analysis by employing both energy and angular momentum conservation laws as well as solid results of latest numerical-relativity simulations and of neutron stars in equilibrium. The new analysis shows that the maximum mass of cold spherical neutron stars can be only weakly constrained as $M_{\rm max} \alt 2.3M_\odot$. Our present result illustrates that the merger remnant neutron star at the onset of collapse to a black hole is not necessarily rapidly rotating and shows that we have to take into account the angular momentum conservation law to impose the constraint on the maximum mass of neutron stars.