Binary Neutron Star Systems Research Articles

The Completeness of Accreting Neutron Star Binary Candidates from the Chinese Space Station Telescope

Abstract A neutron star (NS) has many extreme physical conditions, and one may obtain some important information about an NS via accreting neutron star binary (ANSB) systems. The upcoming Chinese Space Station Telescope (CSST) provides an opportunity to search for a large sample of ANSB candidates. Our goal is to check the completeness of the potential ANSB samples from CSST data. In this paper, we generate some ANSBs and normal binaries under the CSST photometric system by binary evolution and binary population synthesis method and use a machine learning method to train a classification model. Although the Precision (94.56%) of our machine learning model is as high as before study, the Recall is only about 63.29%. The Precision/Recall is mainly determined by the mass transfer rate between the NSs and their companions. In addition, we also find that the completeness of ANSB samples from CSST photometric data by the machine learning method also depends on the companion mass and the age of the system. ANSB candidates with a low initial mass companion star (0.1 M ⊙ to 1 M ⊙) have a relatively high Precision (94.94%) and high Recall (86.32%), whereas ANSB candidates with a higher initial mass companion star (1.1 M ⊙ to 3 M ⊙) have similar Precision (93.88%) and quite low Recall (42.67%). Our results indicate that although the machine learning method may obtain a relatively pure sample of ANSBs, a completeness correction is necessary for one to obtain a complete sample.

Skynet Astromancer Suite: Gravitymancer

The Gravitymancer processing tool within Skynet’s Astromancer suite opens a unique window for students of introductory astronomy courses to explore the most powerful events observed in our Universe – the merger of binary neutron-star systems and black-hole systems. These merger events produce short-lived bursts of gravitational radiation, as observed by one or more of the following gravitational-wave observatories: LIGO, Virgo, and KAGRA. Students can use Gravitymancer to load and interpret the archival event data, finding useful properties like the total mass and mass ratio of the binary system, distance, and inclination angle. While on the surface, students can interpret gravitational-wave events with an easy-to-use GUI, we present here the complex processing utilized to make this tool work.

The Peculiar Precursor of a Gamma-Ray Burst from a Binary Merger Involving a Magnetar

The milestone discovery of GW170817-GRB 170817A-AT 2017gfo has shown that gravitational waves (GWs) could be produced during the merger of a neutron star–neutron star/black hole and that in electromagnetic (EM) waves, a gamma-ray burst (GRB) and a kilonova (KN) are generated in sequence after the merger. Observationally, however, EM properties before the merger phase are still unclear. Here we report a peculiar precursor in a KN-associated long-duration GRB 211211A, providing evidence of the EM before the merger. This precursor lasts ∼0.2 s, and the waiting time between the precursor and the main burst is ∼1 s, comparable to that between GW170817 and GRB 170817A. The spectrum of the precursor could be well fit with a nonthermal cutoff power-law model instead of a blackbody. In particular, a ∼22 Hz quasiperiodic oscillation candidate (∼3σ) is detected in the precursor. These temporal and spectral properties indicate that this precursor is probably produced by a catastrophic flare accompanied with magnetoelastic or crustal oscillations of a magnetar in a binary compact merger. The strong magnetic field of the magnetar can also account for the prolonged duration of GRB 211211A. However, it poses a challenge to reconcile the rather short lifetime of a magnetar with the rather long spiraling time of a binary neutron star system only by the GW radiation before the merger.

Thermally emitting isolated neutron star candidates from the SRG/eROSITA All-Sky Survey

The SRG/eROSITA All-Sky Survey (eRASS) allows for the creation of a complete sample of X-ray dim isolated neutron stars (XDINSs), which will significantly facilitate the study of their population properties, evolution, and connection to other families of isolated neutron stars (INSs). In this work, we conduct a systematic search for XDINSs on the eROSITA western Galactic hemisphere down to an X-ray flux limit of 10$^ and discuss the resulting candidate sample. Consistently with the properties of the known XDINSs, we selected all eRASS sources possessing a soft X-ray spectral distribution and that are unlikely to be associated with optical or infrared sources. Our selection criteria allowed us to recover all known XDINSs and previously proposed candidates. In addition, we put forward 33 new candidate members for dedicated follow-up identification campaigns. We found the resulting candidate sample to be about 30 -- 50<!PCT!> complete, mainly due to source confusion and the stringent cross-matching criteria adopted to select the most promising candidates for immediate follow-up investigation. The candidates of the sample presented here can be broadly divided into two groups: 13 rather soft and 20 hot and somewhat hard X-ray emitters. Interestingly, the remarkably thermal nature of the candidates in the first group as well as their spatial distribution, lack of known counterparts, and absence of significant flux variability agree well with the properties of other confirmed thermally emitting INSs. For the candidates in the second group, the current observational data do not allow one to discern between rotation-powered or recycled pulsars, cataclysmic variables, or quiescent neutron stars in binary systems or even to rule out an extragalactic nature. On the basis of population synthesis and the estimated source completeness of the search, we expect that between one and three new XDINSs are among the already singled-out list of XDINS candidates -- a long-sought increase in the proposed number of members of this elusive class of X-ray emitters.

Pulsar as a Weber detector of gravitational waves and a probe to its internal phase transitions

It is believed that cores of neutron stars provide a natural laboratory where exotic high baryon density phases of quantum chromo dynamics (QCD) may exist. In fact, the theoretically well-established neutron superfluid phase is also believed to be found only inside neutron stars. Focus on neutrons stars has tremendously intensified in recent years with the direct detection of gravitational waves (GWs) by LIGO/Virgo from binary neutron star (BNS) merger events which has allowed the possibility of directly probing the properties of the interior of a neutron star. A truly remarkable phenomenon manifested by rapidly rotating neutron stars is in their avatar as Pulsars. The accuracy of pulsar timing can reach the level of one part in 10[Formula: see text], comparable to that of atomic clocks. Indeed, it was such a great accuracy which had allowed the first indirect detection of GWs from a BNS system. Such an incredible accuracy of pulse timings points to a very interesting possibility. Any deformation of the pulsar, even if it is extremely tiny, has the potential of leaving its imprints on the pulses through introduction of tiny perturbations in the entire moment of inertia (MI) tensor. While, the diagonal components of perturbed MI tensor affect the pulse timings, the off-diagonal components lead to wobbling of pulsar, directly affecting the pulse profile. This opens up a new window of opportunity for exploring various phase transitions occurring inside a pulsar core, through induced density fluctuations, which may be observable as perturbations in the pulse timing as well as its profile. Such perturbations also naturally induce a rapidly changing quadrupole moment of the star, thereby providing a new source of GW emission. Another remarkable possibility arises when we consider the effect of an external GW on neutron star. With the possibility of detecting any minute changes in its configuration through pulse observations, the neutron star has the potential of performing as a Weber detector of GW. This brief review will focus on these specific aspects of a pulsar. Specifically, the focus will be on the type of physics which can be probed by utilizing the effect of changes in the MI tensor of the pulsar on pulse properties.

Cost of Inferred Nuclear Parameters toward the f-mode Dynamical Tide in Binary Neutron Stars

Gravitational-wave (GW) observations from neutron stars (NSs) in a binary system provide an excellent scenario to constrain the nuclear parameters. The investigation of Pratten et al. has shown that the ignorance of f-mode dynamical tidal correction in the GW waveform model of the binary NS system can lead to substantial bias in the measurement of NS properties and NS equations of state. In this work, we investigate the bias in the nuclear parameters resulting from the ignorance of dynamical tidal correction. In addition, this work demonstrates the sensitivity of the nuclear parameters and the estimated constraints on nuclear parameters and NS properties from future GW observations. We infer the nuclear parameters from GW observations by describing the NS matter within the relativistic mean field model. For a population of GW events, we notice that the ignorance of dynamical tide predicts a lower median for nucleon effective mass (m*) by ∼6% compared to the scenario when dynamical tidal correction is considered. Whereas, at a 90% credible interval, m* gets constrained up to ∼5% and ∼3% in A+ (the LIGO-Virgo detectors with a sensitivity of the fifth observing run) and Cosmic Explorer, respectively. We also discuss the resulting constraints on all other nuclear parameters, including compressibility, symmetry energy, and slope of symmetry energy, considering an ensemble of GW events. We do not notice any significant impact in analyzing nuclear parameters other than m* due to the ignorance of f-mode dynamical tides.

GLADEnet: A progressive web app for multi-messenger cosmology and electromagnetic follow-ups of gravitational-wave sources

Multi-messenger astronomy is an emerging field of research aimed at unravelling the physics governing astrophysical transients. GW170817 stands out as the first multi-messenger observation of the coalescence of a binary system of neutron stars, detected by the LIGO and Virgo gravitational-wave interferometers, along with space- and ground-based electromagnetic telescopes. It is a striking example of how multi-messenger observations significantly enhance our understanding of the physics of compact objects, relativistic outflows, and nucleosynthesis. It shows a new way of making cosmology and has the potential to resolve the tension between different measurements of the expansion rate of the Universe. To optimise multi-messenger observational strategies, to evaluate the efficiency of the searches for counterparts, and to identify the host galaxy of the source in a large sky localisation, information about the volumes of galaxies within the gravitational-wave localisation is of paramount importance. This requires the use of galaxy catalogues and appropriate knowledge of their completeness. Here, we describe a new interactive web tool named GLADEnet that allows us to identify catalogued galaxies and to assess the incompleteness of the catalogue of galaxies in real time across the gravitational-wave sky localisation. This measure is of particular importance when using catalogues such as the GLADE catalogue (Galaxy List for the Advanced Detector Era), which includes a collection of various catalogues that make completeness differ across different regions of the sky. We discuss the analysis steps to defining a completeness coefficient and provide a comprehensive guide on how to use the web app, detailing its functionalities. The app is geared towards managing the vast collection of over 22 million objects in GLADE. The completeness coefficient and the GLADE galaxy list will be disseminated in real time via GLADEnet, powered by the Virtual Observatory (VO) standard and tools.

Total r-process Yields of Milky Way Neutron Star Mergers

While it is now known that the mergers of double neutron star binary systems (NSMs) are copious producers of heavy elements, there remains much speculation about whether they are the sole or even principal site of rapid neutron-capture (r-process) nucleosynthesis, one of the primary ways in which heavy elements are produced. The occurrence rates, delay times, and galactic environments of NSMs hold sway over estimating their total contribution to the elemental abundances in the solar system and the Galaxy. Furthermore, the expected elemental yields of NSMs may depend on the merger parameters themselves—such as their stellar masses and radii—which are not currently considered in many galactic chemical evolution models. Using the characteristics of the observed sample of double neutron star (DNS) systems in the Milky Way as a guide, we predict the expected nucleosynthetic yields that a population of DNSs would produce upon merger, and we compare that nucleosynthetic signature to the heavy-element abundance pattern of solar system elements. We find that with our current models, the present DNS population favors the production of lighter r-process elements, while underproducing the heaviest elements relative to the solar system. This inconsistency could imply an additional site for the heaviest elements or a population of DNSs much different from that observed today.

A Theory for Neutron Star and Black Hole Kicks and Induced Spins

Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s−1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s−1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M ⊙, is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M ⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

New and robust gravitational-waveform model for high-mass-ratio binary neutron star systems with dynamical tidal effects

For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the model, which contains a closed-form expression that describes tidal effects, focusing on the description of binary neutron star systems. The model improves upon previous versions by employing a larger set of numerical-relativity data for its calibration, by including high-mass ratio systems covering also a wider range of equations of state. It also takes into account dynamical tidal effects and the known post-Newtonian mass-ratio dependence of individual calibration parameters. We implemented the model in the publicly available uite software library by augmenting different binary black hole waveform models (, , and ). We test the validity of by comparing it with numerical-relativity waveforms, as well as other tidal models. Finally, we perform parameter estimation for GW170817 and GW190425 with the new tidal approximant and find overall consistent results with respect to previous studies. Published by the American Physical Society 2024

Anisotropic energy injection from magnetar central engines in short GRBs

ABSTRACT A long-lived magnetar, potentially originating from a binary neutron star system, has been proposed to explain the extended emission observed in certain short-duration gamma-ray bursts (sGRBs), and is posited as a potential central engine to power the engine-fed kilonovae. Previously, the process by which energy is injected into the surrounding ejecta/jet was widely believed to be nearly isotropic. In this study, we employ special relativity magnetohydrodynamic (SRMHD) simulations to investigate the wind injection process from a magnetar central engine. We explore the dynamics and energy distribution within the system and found that the parameter α = uA/uMWN can be used to indicate the collimation of the magnetar wind energy injection, where uA is the local Alfven four-speed and uMWN is the four-speed of the magnetar wind nebular (MWN) formed from wind-ejecta collision. A significant portion of the injected energy from the magnetar spin-down wind will be channeled to the jet axis due to collimation within the MWN. Achieving isotropic energy injection requires a significantly small α that necessitates either an ultra-relativistic expanding MWN or an extremely low magnetization MWN, both of which are challenging to attain in sGRBs. Consequently, a considerably reduced energy budget (i.e. energy per solid angle reduced by a factor of up to 10 with respect to the value under isotropic assumption) is anticipated to be injected into the ejecta for engine-fed kilonovae. Engine-fed kilonovae would appear fainter than originally anticipated.

Quasiequilibrium configurations of binary systems of dark matter admixed neutron stars

Quasiequilibrium configurations of binary systems of dark matter admixed neutron stars

Effects of nuclear matter properties in neutron star mergers

ABSTRACT The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multimessenger astronomy. Today, a variety of nuclear EOSs are available with various underlying microphysical models. This calls for a study to focus on EOS effects from different physical nuclear matter properties and their influence on BNS mergers. We perform simulations of equal-mass BNS mergers with a set of nine different EOSs based on Skyrme density functionals. In the models, we systematically vary the effective nucleon mass, incompressibility, and symmetry energy at saturation density. This allows us to investigate the influence of specific nuclear matter properties on the dynamics of BNS mergers. We analyse the impact of these properties on the merger dynamics, the fate of the remnant, disc formation, ejection of matter, and GW emission. Our results indicate that some aspects of the merger, such as the frequencies of the post-merger GW spectrum and the shock-heated ejecta mass, are sensitive to the EOS around saturation density while others, such as the contraction of the remnant and the tidal ejecta mass, are sensitive to the behaviour towards higher densities, e.g. characterized by the slope of the pressure versus density. The detailed density dependence of the EOS thus needs to be taken into account to describe its influence on BNS mergers.

General relativistic simulations of high-mass binary neutron star mergers: rapid formation of low-mass stellar black holes

ABSTRACT Almost a hundred compact binary mergers have been detected via gravitational waves by the LIGO-Virgo-KAGRA collaboration in the past few years providing us with a significant amount of new information on black holes and neutron stars. In addition to observations, numerical simulations using newly developed modern codes in the field of gravitational wave physics will guide us to understand the nature of single and binary degenerate systems and highly energetic astrophysical processes. We here presented a set of new fully general relativistic hydrodynamic simulations of high-mass binary neutron star systems using the open-source EINSTEIN TOOLKIT and LORENE codes. We considered systems with total baryonic masses ranging from 2.8$\, \mathrm{M}_\odot$ to 4.0$\, \mathrm{M}_\odot$ and used the SLy equation of state. We analysed the gravitational wave signal for all models and report potential indicators of systems undergoing rapid collapse into a black hole that could be observed by future detectors like the Einstein Telescope and the Cosmic Explorer. The properties of the post-merger black hole, the disc and ejecta masses, and their dependence on the binary parameters were also extracted. We also compared our numerical results with recent analytical fits presented in the literature and provided parameter-dependent semi-analytical relations between the total mass and mass ratio of the systems and the resulting black hole masses and spins, merger frequency, BH formation time, ejected mass, disc mass, and radiated gravitational wave energy.

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.

Erratum: All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems [Phys. Rev. D 103 , 064017 (2021)

Erratum: All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems [Phys. Rev. D <b>103</b> , 064017 (2021)

Event Rate of Fast Radio Bursts from Binary Neutron Star Mergers

It is proposed that a one-off fast radio burst (FRB) with periodic structure may be produced during the inspiral phase of a binary neutron star (BNS) merger. In this paper, we study the event rate of such kind of FRB. We first investigate the properties of two one-off FRBs with periodic structure (i.e., FRB 20191221A and FRB 20210213A) in this scenario, by assuming a fast magnetosonic wave is responsible for their radio emission. For the luminosities and periods of these bursts, it is found that for the NSs in the premerger BNS, magnetic field strengths of B ≳ 1012 G are required. This is relatively high compared with those of most of the BNSs observed in our Galaxy, of which their magnetic fields are around 109 G. Since the observed BNSs in our Galaxy are binaries that have not suffered a merger, a credible event rate of BNS-merger-originated FRBs should be estimated by considering the evolution of both the BNS systems and their magnetic fields. Based on population synthesis and adopting decaying magnetic fields of the NSs, we estimate the event rate of BNS mergers relative to their final magnetic fields. We find that rapidly merging BNSs tend to merge with high magnetization, and the event rate of BNS-merger-originated FRBs, i.e., BNS mergers with both NSs’ magnetic fields being higher than 1012 G, is ∼8 × 104 yr−1 (19% of all BNS mergers) for redshifts z < 1.

Using grey sirens to resolve the Hubble–Lemaître tension

ABSTRACT The measurement of the Hubble–Lemaître constant (H0) from the cosmic microwave background and the Type IA supernovae are at odds with each other. One way to resolve this tension is to use an independent way to measure H0. This can be accomplished by using gravitational-wave (GW) observations. Previous works have shown that with the onset of the next generation of GW detector networks, it will be possible to constrain H0 to better than 2 per cent precision (which is enough to resolve the tension) with binary black hole systems that are extremely well localized in the sky, also called golden dark sirens. Bright sirens like binary neutron star systems can also help resolve the tension if both the GW and the following electromagnetic counterpart are detected. In this work, we show that neutron star-black hole (NSBH) mergers can act both as golden dark sirens as well as bright sirens, thus, assigning them the term grey sirens. We assess the potential of using NSBH mergers to measure H0 and find that the Voyager network might be able to resolve the tension in an observation span of 5 yr. The next generation networks, which include the Cosmic Explorer detectors and the Einstein Telescope will be able to measure H0 to sub-per cent level just by using NSBH mergers.

Searching For Black Holes Using Auto Differentiation

This study presents GravAD, a novel approach for detecting gravitational waves using automatic differentiation and JAX. GravAD demonstrates comparable signal-to-noise ratio and mass values to established LIGO pipelines with a significant reduction in the number of templates. Limitations include the inability to handle binary neutron star systems and some lower-mass black holes. Leveraging JAX’s acceleration, GravAD offers potential as a rapid preliminary tool for gravitational wave detection. Future work includes further optimisation of functions, exploration of alternative optimisation algorithms, real-time data analysis adaptation, and expanding the scope to handle a broader range of astrophysical sources.

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.