Binary Neutron Star Systems Research Articles

Improved all-sky search method for continuous gravitational waves from unknown neutron stars in binary systems

Continuous gravitational waves from spinning deformed neutron stars have not been detected yet and are one of the most promising signals for future detection. All-sky searches for continuous gravitational waves from unknown neutron stars in binary systems are the most computationally challenging search type. Consequently, very few search algorithms and implementations exist for these sources, and only a handful of such searches have been performed so far. In this paper, we present a new all-sky binary search method, BinarySkyHou$\mathcal{F}$, which extends and improves upon the earlier BinarySkyHough method and which was the basis for a recent search [Covas et al., Astrophys. J. Lett. 929, L19 (2022)]. We compare the sensitivity and computational cost to the previous method, showing that it is both more sensitive and computationally efficient, which allows for broader and more sensitive searches.

Black hole mass and spin measurements through the relativistic precession model: XTE J1859+226

ABSTRACT The X-ray light curves of accreting black holes and neutron stars in binary systems show various types of quasi-periodic oscillations (QPOs), the origin of which is still debated. The relativistic precession model identifies the QPO frequencies with fundamental time-scales from General Relativity, and has been proposed as a possible explanation of certain types of such oscillations. Under specific conditions (i.e. the detection of a particular QPOs triplet) such a model can be used to obtain self-consistent measurements of the mass and spin of the compact object. So far this has been possible only in the black hole binary GRO J1655-40. In the RXTE/PCA data from the 1999–2000 outburst of the black hole transient XTE J1859+226 we found a QPO triplet, and used the relativistic precession model to obtain high-precision measurements of the black hole mass and spin – M = (7.85 ± 0.46) M⊙, a* = 0.149 ± 0.005 – the former being consistent with the most recent dynamical mass determination from optical measurements. Similarly to what has been already observed in other black hole systems, the frequencies of the QPOs and broad-band noise components match the general relativistic frequencies of particle motion close to the compact object predicted by the model. Our findings confirm previous results and further support the validity of the relativistic precession model, which is the only electromagnetic-measurement-based method that so far has consistently yielded spins close to those from the gravitational waves produced by merging binary black holes.

Bare collapse, formation of neutron star binaries and fast optical transients

ABSTRACT ‘Bare collapse’, the collapse of a bare stellar core into a neutron star with a very small mass ejection links two seemingly unrelated phenomena: the formation of binary neutron star (BNS) systems and the observations of fast and luminous optical transients. We carried out calculations of the collapse due to electron-capture of both evolutionary and synthetic isentropic bare stellar cores. We find that the collapse results in the formation of a light ∼1.3M⊙ neutron star and an ejection of ∼0.1M⊙ at ∼0.1c. The outer shell of the ejecta is composed of 56Ni that can power an ultra-stripped supernova. The models we explored can explain most of the observed fast optical flares but not the brightest ones. Collapse of cores surrounded by somewhat more massive envelopes can produce larger amounts of 56Ni and explain brighter flares. Alternatively, those events can arise due to interaction of the very energetic ejecta with winds that were ejected from the progenitor a few days before the collapse.There are no data bases used.

Nonlinear mode-tide coupling in coalescing binary neutron stars with relativistic corrections

We compute the internal modes of a nonspinning neutron star and its tidal metric perturbation in general relativity, and determine the effect of relativistic corrections to the modes on mode coupling. Claims have been made that a new hydrodynamic instability can occur in a neutron star in a binary neutron star system triggered by the nonlinear coupling of the companion's tidal field to pairs of p-modes and g-modes in it as the binary inspirals toward merger. This "PG" instability may be significant since it can influence the binary's inspiral phase by extracting orbital energy, thereby potentially causing large deviations in their gravitational waveforms from those predicted by theoretical models that do not account for it. This can result in incorrect parameter estimation, at best, or mergers going undetected, at worst, owing to the use of deficient waveform models. On the other hand, better modeling of this instability and its effect on binary orbits can unravel a new phenomenon and shed light on stellar instabilities, via gravitational wave observations. So far, all mode-tide coupling instability studies have been formulated in Newtonian perturbation theory. Neutron stars are compact objects, so relativistic corrections might be important. We present and test a new code to calculate the relativistic eigenmodes of nonrotating relativistic stars. We use these relativistic tide and neutron star eigenmodes to compute the mode-tide coupling strength (MTCS) for a few selected equations of state. The MTCS thus calculated can be at most tens of percent different from its purely Newtonian value, but we confirm the dependencies on orbital separation and equation of state found by Newtonian calculations. For some equations of state, the MTCS is very sensitive to the neutron star crust region, demonstrating the importance of treating this region accurately.

Numerical relativity for Horndeski gravity

In this paper, we present an overview of recent developments in the numerical solution of Horndeski gravity theories, which are the class of all scalar–tensor theories of gravity that have second-order equations of motion. We review several methods that have been used to establish well-posed initial value problems for these theories, and discuss well-posed formulations of the constraint equations. We also discuss global aspects of exact, strongly coupled solutions to some of the Horndeski gravity theories: the formation of shocks, the loss of hyperbolicity, and the formation of naked curvature singularities. Finally, we discuss numerical solutions to binary black hole and neutron star systems for several Horndeski theories.

Modelling the host galaxies of binary compact object mergers with observational scaling relations

ABSTRACT The merger rate density evolution of binary compact objects and the properties of their host galaxies carry crucial information to understand the sources of gravitational waves. Here, we present galaxy$\mathcal {R}$ate, a new code that estimates the merger rate density of binary compact objects and the properties of their host galaxies, based on observational scaling relations. We generate our synthetic galaxies according to the galaxy stellar mass function. We estimate the metallicity according to both the mass–metallicity relation (MZR) and the fundamental metallicity relation (FMR). Also, we take into account galaxy–galaxy mergers and the evolution of the galaxy properties from the formation to the merger of the binary compact object. We find that the merger rate density changes dramatically depending on the choice of the star-forming galaxy main sequence, especially in the case of binary black holes (BBHs) and black hole neutron star systems (BHNSs). The slope of the merger rate density of BBHs and BHNSs is steeper if we assume the MZR with respect to the FMR, because the latter predicts a shallower decrease of metallicity with redshift. In contrast, binary neutron stars (BNSs) are only mildly affected by both the galaxy main sequence and metallicity relation. Overall, BBHs and BHNSs tend to form in low-mass metal-poor galaxies and merge in high-mass metal-rich galaxies, while BNSs form and merge in massive galaxies. We predict that passive galaxies host at least ∼5–10 per cent, ∼15–25 per cent, and ∼15–35 per cent of all BNS, BHNS, and BBH mergers in the local Universe.

Investigating the mass-ratio dependence of the prompt-collapse threshold with numerical-relativity simulations

The next observing runs of advanced gravitational-wave detectors will lead to a variety of binary neutron star detections and numerous possibilities for multimessenger observations of binary neutron star systems. In this context a clear understanding of the merger process and the possibility of prompt black hole formation after merger is important, as the amount of ejected material strongly depends on the merger dynamics. These dynamics are primarily affected by the total mass of the binary, however, the mass ratio also influences the postmerger evolution. To determine the effect of the mass ratio, we investigate the parameter space around the prompt-collapse threshold with a new set of fully relativistic simulations. The simulations cover three equations of state and seven mass ratios in the range of $1.0\ensuremath{\le}q\ensuremath{\le}1.75$, with five to seven simulations of binary systems of different total mass in each case. The threshold mass is determined through an empirical relation based on the collapse time, which allows us to investigate effects of the mass ratio on the threshold mass and also on the properties of the remnant system. Furthermore, we model effects of mass ratio and equation of state on tidal parameters of threshold configurations.

Thermodynamics of Hot Neutron Stars and Universal Relations

Over the last few years, the detection of gravitational waves from binary neutron star systems has rekindled our hopes for a deeper understanding of the unknown nature of ultradense matter. In particular, gravitational wave constraints on the tidal deformability of a neutron star can be translated into constraints on several neutron star properties using a set of universal relations. Apart from binary neutron star mergers, supernova explosions are also important candidates for the detection of multimessenger signals. Such observations may allow us to impose significant constraints on the binding energy of neutron stars. The purpose of the present study is twofold. Firstly, we investigate the agreement of finite temperature equations of state with established universal relations. Secondly, we examine the possible existence of a universal relation between the binding energy and the dimensionless tidal deformability, which are the bulk properties connected to the most promising sources for multimessenger signals. We find that hot equations of state are not always compatible with accepted universal relations. Therefore, the use of such expressions for probing general relativity or imposing constraints on the structure of neutron stars would be inconclusive (when thermal effects are present). Additionally, we show that the binding energy and the dimensionless tidal deformability exhibit a universal trend at least for moderate neutron star masses. The latter allows us to set bounds on the binding energy of a 1.4 M⊙ neutron star using data from the GW170817 event. Finally, we provide a relation between the compactness, the binding energy and the dimensionless tidal deformability of a neutron star that is accurate for cold and hot isentropic equations of state.

High-accuracy high-mass-ratio simulations for binary neutron stars and their comparison to existing waveform models

The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical relativity simulations for four different physical setups with mass ratios $q=1.25$, 1.50, 1.75, 2.00, and total gravitational mass $M=2.7\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately second-order converging results for the dominant (2,2) mode, but also the subdominant (2,1), (3,3), and (4,4) modes, while generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current nrtidal model used to describe tidal effects is a valid description for high-mass-ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.

Incorporating a Radiative Hydrodynamics Scheme in the Numerical-Relativity Code BAM

To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHBΛϕ. Overall, our test simulations have good agreement with those reported in the literature.

Thermal effects on tidal deformability in the last orbits of an inspiraling binary neutron star system

The study of binary neutron stars mergers by the detection of the emitted gravitational waves is one of the most promised tools to study the properties of dense nuclear matter at high densities. It is worth claiming that strong evidence that the temperature of the stars during the last orbits before coalescing is very low, T≪1 MeV (hereafter cold case or cold star; kB=1), does not exist. Nevertheless, theoretical studies suggest that the temperature concerning the inspiral phase, could reach even a few MeV. The heating process of the interior of the neutron stars is as follows; tides transfer mechanical energy and angular momentum to the star at the expense of the orbit, where friction within the star converts the mechanical energy into heat. During the inspiral, these effects are potentially detectable. Different treatments have been used to estimate the transfer of the mechanical energy and the size of the tidal friction, leading to different conclusions about the importance of pre-merger tidal effects. The present work is dedicated to the study of the effect of temperature on the tidal deformability of neutron stars during the inspiral of a neutron star system just before the merger. We applied a class of hot equations of state, both isothermal and adiabatic, originated from various nuclear models. We found that even for low values of temperature (T<1 MeV), the effects on the basic ingredients of tidal deformability are not negligible. On the other hand, in the case of the adiabatic star, the thermal effects on tidal deformability remain imperceptible, up to the entropy per baryon value S=0.2 (in units of kB). According to the main finding, the effects of the temperature on tidal deformability, at least for low temperatures, are almost negligible. This is a surprising effect since both radius and Love number, the basic ingredients of the tidal deformability, are sensitive to temperature. The consequences of the above result are discussed and analyzed.

Systematic bias due to eccentricity in parameter estimation for merging binary neutron stars

We study the impact of eccentricity on gravitational-wave parameter estimation for binary neutron star systems. For signals with small eccentricity injected into the advanced LIGO sensitivity, we perform Bayesian parameter estimation using the circular waveform model and show how the recovered parameters can be biased from their true values, focusing on the intrinsic parameters the chirp mass ($M_c$), the symmetric mass ratio ($\eta$), and the tidal deformability ($\tilde{\lambda}$). By comparing the results between the Bayesian and the analytic Fisher-Cutler-Vallisneri (FCV) methods, we obtain the valid criteria for the FCV approach. Employing the FCV method and using the realistic population of binary neutron star sources distributed in the $m_1$-$m_2$-$e_0$ space, where $e_0$ indicates the eccentricity at 10Hz, we calculate the measurement errors ($\sigma_{\theta}$) and the systematic biases ($\Delta \theta/\sigma_{\theta}$) and obtain their generalized distributions in the range of $0 \leq e_0 \leq 0.025$. We find that for all of the three parameters, the biases increase with increasing $e_0$, and this increase is faster for larger $e_0$. The bias is mainly dependent on the value of $e_0$ and weakly dependent on the component masses, and thus the distribution shows a narrow band in the $e_0$-$\Delta \theta/\sigma_{\theta}$ plane. We present various posterior examples to illustrate our new findings, such as the bimodality of posteriors. In particular, we give a specific injection-recovery example to demonstrate the importance of including eccentricity in parameter estimation to avoid incorrect predictions of the neutron star equation of state.

Constraints on r-modes and Mountains on Millisecond Neutron Stars in Binary Systems

Continuous gravitational waves are nearly monochromatic signals emitted by asymmetries in rotating neutron stars. These signals have not yet been detected. Deep all-sky searches for continuous gravitational waves from isolated neutron stars require significant computational expense. Deep searches for neutron stars in binary systems are even more expensive, but these targets are potentially more promising emitters, especially in the hundreds of Hertz region, where ground-based gravitational-wave detectors are most sensitive. We present here an all-sky search for continuous signals with frequency between 300 and 500 Hz, from neutron stars in binary systems with orbital periods between 15 and 60 days and projected semimajor axes between 10 and 40 lt-s. This is the only binary search on Advanced LIGO data that probes this frequency range. Compared to previous results, our search is over an order of magnitude more sensitive. We do not detect any signals, but our results exclude plausible and unexplored neutron star configurations, for example, neutron stars with relative deformations greater than 3 × 10−6 within 1 kpc from Earth and r-mode emission at the level of α ∼ a few 10−4 within the same distance.

NEUTRON STAR BINARYAS GRAVITATIONAL QUANTUM SCATTERING SYSTEM

Using ensemble of N gravitons, we evaluate the integrated density of states for gravitons arising from the spiral-in motion of gravitationally bound neutron star (NS) binaries, under the Dirichlet boundary conditions. Consequently, the two massive objects making up of the binary system are immersed in a gravito-bosonic environment. The corrections to the density of states for gravitons give expression to the quantum fluctuations of gravitational fields, i.e., the gravitational Casimir energy. By considering the extended scales of the star components, we find the gravitational Casimir effects are larger than that obtained by treating them as point particles and that the gravitational Casimir energy depends both on the scales of objects and on the separation of the binary, which increases for lager scales of star components in a typical NS binary system, with separation of , i.e. low-frequency gravitational wave sources.

Discriminating same-mass neutron stars and black holes gravitational waveforms

Gravitational wave-forms from coalescences of binary black hole systems and binary neutron star systems with low tidal effects can hardly be distinguished if the two systems have similar masses. In the absence of discriminating power based on the gravitational wave-forms, the classification of sources into binary neutron stars, binary black holes and mixed systems containing a black hole and a neutron star can only be unambiguous when assuming the standard model of stellar evolution and using the fact that there exists a mass gap between neutron stars and black holes. This approach is however limited by its own assumptions: for instance the 2.6 solar mass object detected in the GW190814 event remains unclassified, and models of new physics can introduce new compact objects, like primordial black holes, which may have masses in the same range as neutron stars. Then, without an electromagnetic counterpart (kilonova), classifying mergers of compact objects without mass gap criteria remains a difficult task, unless the source is close enough. In what follows we investigate a procedure to discriminate a model between binary neutron star merger and primordial binary black hole merger by using a Bayes factor in simulated wave-forms that we superimpose to realistic detector noise.

High-accuracy simulations of highly spinning binary neutron star systems

With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $\chi \sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections are combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyse them with a set of four different waveform models. Based on our analysis with about 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to $\pm 250\rm m$ ($2\%$ at $90\%$ credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.

Pulsar Star Detection: A Comparative Analysis of Classification Algorithms using SMOTE

A Pulsar is a highly magnetized rotating compact star whose magnetic poles emit beams of radiation. The application of pulsar stars has a great application in the field of astronomical study. Applications like the existence of gravitational radiation can be indirectly confirmed from the observation of pulsars in a binary neutron star system. Therefore, the identification of pulsars is necessary for the study of gravitational waves and general relativity. Detection of pulsars in the universe can help research in the field of astrophysics. At present, there are millions of pulsar candidates present to be searched. Machine learning techniques can help detect pulsars from such a large number of candidates. The paper discusses nine common classification algorithms for the prediction of pulsar stars and then compares their performances using various classification metrics such as classification accuracy, precision and recall value, ROC score and f-score on both balanced and unbalanced data. SMOTE-technique is used to balance the data for better results. Among the nine algorithms, XGBoosting algorithm achieved the best results. The paper is concluded with prospects of machine learning for pulsar detection in the field of astronomy.

Multimessenger Constraints on Magnetic Fields in Merging Black Hole–Neutron Star Binaries

Abstract The LIGO–Virgo–KAGRA Collaboration recently detected gravitational waves (GWs) from the merger of black hole–neutron star (BHNS) binary systems GW200105 and GW200115. No coincident electromagnetic (EM) counterparts were detected. While the mass ratio and BH spin in both systems were not sufficient to tidally disrupt the NS outside the BH event horizon, other, magnetospheric mechanisms for EM emission exist in this regime and depend sensitively on the NS magnetic field strength. Combining GW measurements with EM flux upper limits, we place upper limits on the NS surface magnetic field strength above which magnetospheric emission models would have generated an observable EM counterpart. We consider fireball models powered by the black hole battery mechanism, where energy is output in gamma rays over ≲1 s. Consistency with no detection by Fermi-GBM or INTEGRAL SPI-ACS constrains the NS surface magnetic field to ≲1015 G. Hence, joint GW detection and EM upper limits rule out the theoretical possibility that the NSs in GW200105 and GW200115, and the putative NS in GW190814, retain dipolar magnetic fields ≳1015 G until merger. They also rule out formation scenarios where strongly magnetized magnetars quickly merge with BHs. We alternatively rule out operation of the BH-battery-powered fireball mechanism in these systems. This is the first multimessenger constraint on NS magnetic fields in BHNS systems and a novel approach to probe fields at this point in NS evolution. This demonstrates the constraining power that multimessenger analyses of BHNS mergers have on BHNS formation scenarios, NS magnetic field evolution, and the physics of BHNS magnetospheric interactions.

Periodic Repeating Fast Radio Bursts: Interaction between a Magnetized Neutron Star and Its Planet in an Eccentric Orbit

Abstract Fast radio bursts (FRBs) are mysterious transient phenomena. The study of repeating FRBs may provide useful information about their nature due to their redetectability. The two most famous repeating sources are FRBs 121102 and 180916, with a period of 157 days and 16.35 days, respectively. Previous studies suggest that the periodicity of FRBs is likely associated with neutron star (NS) binary systems. Here we introduce a new model which proposes that periodic repeating FRBs are due to the interaction of a NS with its planet in a highly elliptical orbit. The periastron of the planet is very close to the NS so that it would be partially disrupted by tidal force every time it passes through the periastron. Fragments generated in the process could interact with the compact star through the Alfvén wing mechanism and produce FRBs. The model can naturally explain the repeatability of FRBs, with a period ranging from a few days to several hundred days, but it generally requires that the eccentricity of the planet’s orbit should be large enough. Taking FRBs 121102 and 180916 as examples, it is shown that the main features of the observed repeating behaviors can be satisfactorily accounted for.