Stars In Binary Systems Research Articles

On the Formation of the Double Neutron Star Binary PSR J1846-0513

Abstract The Double Neutron Star binary PSR J1846-0513 is discovered by the Five hundred-meter Aperture Spherical radio Telescope (FAST) in Commensal Radio Astronomy FAST Survey. The pulsar is revealed to be harbored in an eccentric orbit with e = 0.208 and orbital period of 0.613 day. The total mass of the binary is constrained to be 2.6287(35)M⊙, with the upper mass limit of the pulsar at 1.3455 M⊙ and the lower mass limit of the com panion at 1.2845 M⊙. To reproduce its evolution history, we simulated the 1D model of its progenitor from neutron star - helium (He) star binary system via MESA code. Since the large eccentricity is widely believed to originate from an asymmetric supernova explosion, we examined the dynamical effects of the supernova. Our results show that the binary may
stem from a helium star with initial mass of 3.3 − 4.0 M⊙ accompanied by a neutron star in a circle orbit with initial period of ∼ 0.5 day.

Binary progenitor systems for Type Ic supernovae

Core-collapse supernovae are explosions of massive stars at the end of their evolution. They are responsible for metal production and for halting star formation, having a significant impact on galaxy evolution. The details of these processes depend on the nature of supernova progenitors, but it is unclear if Type Ic supernovae (without hydrogen or helium lines in their spectra) originate from core-collapses of very massive stars (>30 M⊙) or from less massive stars in binary systems. Here we show that Type II (with hydrogen lines) and Ic supernovae are located in environments with similar molecular gas densities, therefore their progenitors have comparable lifetimes and initial masses. This supports a binary interaction for most Type Ic supernova progenitors, which explains the lack of hydrogen and helium lines. This finding can be implemented in sub-grid prescriptions in numerical cosmological simulations to improve the feedback and chemical mixing.

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.

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.

A Channel to Form Fast-spinning Black Hole–Neutron Star Binary Mergers as Multimessenger Sources. II. Accretion-induced Spin-up

In this work, we investigate an alternative channel for the formation of fast-spinning black hole–neutron star (BHNS) binaries, in which super-Eddington accretion is expected to occur in accreting BHs during the stable mass transfer phase within BH-stripped helium (BH–He-rich) star binary systems. We evolve intensive MESA grids of close-orbit BH–He-rich star systems to systematically explore the projected aligned spins of BHs in BHNS binaries, as well as the impact of different accretion limits on the tidal disruption probability and electromagnetic (EM) signature of BHNS mergers. Most of the BHs in BHNS mergers cannot be effectively spun up through accretion if the accretion rate is limited to ≲10ṀEdd , where ṀEdd is the standard Eddington accretion limit. In order to reach high spins (e.g., χ BH ≳ 0.5), the BHs are required to be born less massive (e.g., ≲3.0 M ⊙) in binary systems with initial periods of ≲0.2–0.3 days and accrete material at ∼100ṀEdd . However, even under this high accretion limit, ≳6 M ⊙ BHs are typically challenging to significantly spin up and generate detectable associated EM signals. Our population simulations suggest that different accretion limits have a slight impact on the ratio of tidal disruption events. However, as the accretion limit increases, the EM counterparts from the cosmological BHNS population can become bright overall.

Stellar Cycle and Evolution of Polar Spots in an M+WD Binary

Stellar activity cycles reveal continuous relaxation and induction of magnetic fields. The activity cycle is typically traced through the observation of cyclic variations in total brightness or Ca H&K emission flux of stars, as well as cyclic variations in the orbital periods of binary systems. In this work, we report the identification of a semidetached binary system (TIC 16320250) consisting of a white dwarf (0.67 M ⊙) and an active M dwarf (0.56 M ⊙). The long-term multiband optical light curves spanning twenty years revealed three repeated patterns, suggestive of a possible activity cycle of about 10 years of the M dwarf. Light-curve fitting indicates the repeated variation is caused by the evolution, particularly the motion, of polar spots. The significant Ca H&K, Hα, ultra-violet, and X-ray emissions imply that the M dwarf is one of the most magnetically active stars. We propose that in the era of large time-domain photometric sky surveys (e.g., ASAS-SN, Zwicky Transient Facility, LSST, Sitian), long-term light-curve modeling can be a valuable tool for tracing and revealing stellar activity cycle, especially for stars in binary systems.

Underestimation of the tidal force and apsidal motion in close binary systems by the perturbative approach: Comparisons with non-perturbative models

Context. Stellar deformations play a significant role in the dynamical evolution of stars in binary systems, impacting the tidal dissipation and the outcomes of mass transfer processes. The prevalent method for modelling the deformations and tidal interactions of celestial bodies solely relies on the perturbative approach, which assumes that stellar deformations are minor perturbations to the spherical symmetry. An observable consequence of stellar deformations is the apsidal motion in eccentric systems, which has be observationally determined across numerous binary systems. Aims. Our objective is to assert the reliability of the perturbative approach when applied to close and strongly deformed binary systems. Methods. We have developed a non-perturbative 3D modelling method designed to account for high stellar deformations. We focus on comparing the properties of perturbatively deformed stellar models with our 3D models, particularly in terms of apsidal motion. Results. Our research highlights that the perturbative model becomes imprecise and underestimates the tidal force and rate of apsidal motion at a short orbital separation. This discrepancy primarily results from the first-order treatment in the perturbative approach, and cannot be rectified using straightforward mathematical corrections due to the strong non-linearity and numerous parameters of the problem. We have determined that our methodology affects the modelling of approximately 42% of observed binary systems with measured apsidal motion, introducing a discrepancy greater than 2% when the normalised orbital separation verifies q−1/5a(1 − e2)/R1 ≲ 6.5 (q is the mass ratio of the system, a is its semi-major axis, e is its orbital eccentricity and R1 is the radius of the primary star). Conclusions. The perturbative approach underestimates tidal interactions between bodies up to ∼40% for close low-mass binaries. All the subsequent modelling is impacted by our findings, in particular, the tidal dissipation is significantly underestimated. As a result, all binary stellar models are imprecise when applied to systems with a low orbital separation, and the outcomes of these models are also affected by these inaccuracies.

Explosions in Roche-lobe Distorted Stars: Relativistic Bullets in Binaries

State-of-the-art surveys reveal that most massive stars in the Universe evolve in close binaries. Massive stars in such systems are expected to develop aspherical envelopes due to tidal interactions and/or rotational effects. Recently, it was shown that point explosions in oblate stars can produce relativistic equatorial ring-like outflows. Moreover, since stripped-envelope stars in binaries can expand enough to fill their Roche lobes anew, it is likely that these stars die with a greater degree of asphericity than the oblate spheroid geometry previously studied. We investigate the effects of this asymmetry by studying the gas dynamics of axisymmetric point explosions in stars in various stages of filling their Roche lobes. We find that point explosions in these pear-shaped stars produce transrelativistic ejecta that coalesce into bullets pointed both toward and away from the binary companion. We present this result and comment on key morphological differences between core-collapse explosions in spherical versus distorted stars in binary systems, effects on gravitational wave sources, and observational signatures that could be used to glean these explosion geometries from current and future surveys.

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.

Integral Fluxes of Neutrinos and Gamma-Rays Emitted from Neighboring X-ray Binaries

Astrophysical plasma ejections (jets) are formed and powered by black holes that accrete material from their companion star in binary systems. Black hole X-ray binary systems constitute potential powerful galactic and extragalactic neutrino and gamma-ray sources. After being accelerated to highly relativistic velocities and subjected to various energy-consuming interactions, the lepto-hadronic content of the jets produces secondary particles such as pions and muons that decay to gamma-ray photons and neutrinos heading towards the Earth. In this work, we employ a jet emission model in order to predict the neutrino and gamma-ray integral fluxes emanating from some of the most investigated and prominent stellar black hole X-ray binary systems in the Milky Way, such as GRO J1655-40, Cygnus X-1, SS 433, and GRS 1915+105. For the sake of comparison, we also include an extragalactic system, namely, LMC X-1, located in the Large Magellanic Cloud. For the case of gamma-ray emissions, we also include absorption effects due to X-ray emission from the accretion disk and the black hole corona, as well as ultraviolet (UV) emission from the binary system’s companion star.

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)

Properties of Binary Systems in a One-dimensional Approximation

Evolutionary calculations for stars in close binary systems are in high demand to obtain better constraints on gravitational-wave source progenitors, understand transient events from stellar interactions, and more. Modern one-dimensional (1D) stellar codes make use of the Roche lobe radius R L concept in order to treat stars in binary systems. If the stellar companion is approaching its R L, mass transfer treatment is initiated. However, the effective acceleration also affects the evolution of a star in a close binary system. This is different from the gravity inside a single star, whether that single star is rotating or not. Here, we present numerically obtained tables of properties of stars in a binary system as a function of the effective potential: volume-equivalent radii of the equipotential surfaces, effective accelerations and the inverse effective accelerations averaged over the same equipotential surfaces, and the properties of the L 1-plane cross sections. The tables are obtained for binaries where the ratios of the primary star mass to the companion star mass are from 10−6 to 105 and include equipotential surfaces up to the star’s outer Lagrangian point. We describe the numerical methods used to obtain these quantities and report how we verified the numerical results. We also describe and verify the method to obtain the effective acceleration for nonpoint mass distributions. We supply a sample code showing how to use our tables to get the average effective accelerations in 1D stellar codes.

New constraints on the kinematic, relativistic, and evolutionary properties of the PSR J1757−1854 double neutron star system

ABSTRACT PSR J1757−1854 is one of the most relativistic double neutron star binary systems known in our Galaxy, with an orbital period of $P_\text{b}=4.4\, \text{h}$ and an orbital eccentricity of e = 0.61. As such, it has promised to be an outstanding laboratory for conducting tests of relativistic gravity. We present the results of a 6-yr campaign with the 100-m Green Bank and 64-m Parkes radio telescopes, designed to capitalize on this potential. We identify secular changes in the profile morphology and polarization of PSR J1757−1854, confirming the presence of geodetic precession and allowing the constraint of viewing geometry solutions consistent with General Relativity. We also update PSR J1757−1854’s timing, including new constraints of the pulsar’s proper motion, post-Keplerian parameters, and component masses. We conclude that the radiative test of gravity provided by PSR J1757−1854 is fundamentally limited to a precision of 0.3 per cent due to the pulsar’s unknown distance. A search for pulsations from the companion neutron star is also described, with negative results. We provide an updated evaluation of the system’s evolutionary history, finding strong support for a large kick velocity of $w\ge 280\, \rm{km\,s}^{-1}$ following the second progenitor supernova. Finally, we reassess PSR J1757−1854’s potential to provide new relativistic tests of gravity. We conclude that a 3-σ constraint of the change in the projected semimajor axis ($\dot{x}$) associated with Lense–Thirring precession is expected no earlier than 2031. Meanwhile, we anticipate a 3-σ measurement of the relativistic orbital deformation parameter δθ as soon as 2026.

A small survey for massive stars in the LMC/SMC: binaries and pulsations

ABSTRACT Binarity and pulsations are two notable properties that frequently appear together in massive stars. To place constraints on how binarity affects stellar evolution, the exact parameters of massive stars in binary systems must be known. Their exact and accurate masses and radii, which were calculated from binary dynamics, present model-independent restrictions and challenge existing theories of stellar evolution. In this investigation, our goal is to determine the atmospheric characteristics of nine double-lined spectroscopic binaries that are located in the Large Magellanic Cloud and Small Magellanic Cloud. Nine newly discovered double-lined eclipsing binaries with B-type massive components have been studied; these binaries were found using the OGLE variable star catalogues. The absolute physical parameters of each component have been measured with an accuracy of less than or equal to 3 per cent. By analysing through Transiting Exoplanet Survey Satellite data of the systems, we were able to identify five candidates for slowly pulsating B stars in massive binary systems with masses ranging from 2 to 25 M⊙. According to the results of the preliminary research, the orbital periods of these binary systems span anywhere from 1.8 to 6.3 d, whereas the pulsating periods range anywhere from 0.3 to 1.5 d. The achievement in identifying a significant number of pulsating eclipsing systems opens the door to the initial sensible approach of constraining the internal physics of a group of massive stars via asteroseismology of eclipsing binaries including such stars.

Expectations for fast radio bursts in neutron star–massive star binaries

Context. Recent observations of a small sample of repeating fast radio bursts (FRBs) have revealed a periodicity in their bursting activity that suggests a binary origin for the modulation. Aims. We set out to explore the scenario where a subset of repeating FRBs originates in binary systems that host a highly energetic neutron star and a massive companion star, akin to γ-ray binaries and young high-mass X-ray binaries. Methods. In this scenario, we infer observables, compare them with current observational constraints, and make predictions for future observations. Firstly, we specifically focused on the host galaxy properties and binary formation rates. Subsequently, we investigated the expected evolution of the rotation and dispersion measure in this scenario, the predicted birth site offsets, and the origin of the persistent radio emission observed in a subset of these systems. Results. The host galaxies for repeating FRBs favour the formation of neutron star–massive star binary systems, but any conclusive evidence will require future discoveries and localisations of FRBs. The birth rate of high-mass X-ray binaries, used as a proxy for all considered binaries, significantly exceeds the estimated rate of FRBs, which can be explained if only a small subset of these systems produce FRBs. We show that, under simple assumptions, we can reproduce the dispersion measure and rotation measure evolution that is seen in a subset of repeating FRBs. We also discuss the possibility of detecting a persistent radio source associated with the FRB due to an intra-binary shock between the companion star wind and either the pulsar wind or giant magnetar flares. The observed long-term luminosity stability of the persistent radio sources is most consistent with a giant flare-powered scenario. However, this explanation is highly dependent on the magnetic field properties of the neutron star. Conclusions. With these explorations, we provide a framework to discuss future FRB observations in the context of neutron star–massive star binary scenarios. We conclude that more localisations and observations of repeaters will be necessary to conclusively determine or rule out a connection between (repeating) FRBs and such binaries.

Combining HIPPARCOS and Gaia data for the study of binaries: The BINARYS tool

Context. Orbital motion in binary and planetary systems is the main source of precise stellar and planetary mass measurements, and the joint analysis of data obtained using multiple observational methods can both lift degeneracies and improve precision. Aims. We set out to measure the masses of individual stars in binary systems using all the information brought by the HIPPARCOS and Gaia absolute astrometric missions. Methods. We present BINARYS, a tool that uses the HIPPARCOS and Gaia absolute astrometric data and combines them with relative astrometry and/or radial velocity measurements to determine the orbit of a binary system. This tool rigorously combines the HIPPARCOS and Gaia data (here EDR3) and can use the HIPPARCOS Transit Data as needed for binaries where HIPPARCOS detects significant flux from the secondary component. It also supports the case where Gaia has resolved the system, giving an astrometric solution for both components. Results. We determine model-independent individual masses for the first time for three systems: the two mature binaries Gl 494 (M1 = 0.584 ± 0.003 M⊙ and M2 = 87 ± 1 MJup) and HIP 88745 (M1 = 0.96 ± 0.02 M⊙ and 0.62-0.008+0.009 ), and the younger AB Dor member GJ 2060 (1926-6+7 and 0.882 -0.005+0.004). The latter provides a rare test of evolutionary model predictions at young ages in the low-stellar-mass range and sets a lower age limit of 100 Myr for the moving group.

Constraining the anomalous coupling of gravitational waves with the double pulsar

We revisit the decay of the orbital period in binary systems that occurs due to the emission of gravitational waves in the context of modifie- gravity models where the coupling ${G}_{gw}$ between matter and on shell gravitons is allowed to differ from the Newton constant ${G}_{N}$. Using the most precisely-timed double neutron star binary system, those of the double pulsar, we constrain the ratio ${G}_{gw}/{G}_{N}$ to the level of ${10}^{\ensuremath{-}4}$, improving by two orders of magnitude the present bound on this quantity.

An observationally derived kick distribution for neutron stars in binary systems

ABSTRACT Understanding the natal kicks received by neutron stars (NSs) during formation is a critical component of modelling the evolution of massive binaries. Natal kicks are an integral input parameter for population synthesis codes, and have implications for the formation of double NS systems and their subsequent merger rates. However, many of the standard observational kick distributions that are used are obtained from samples created only from isolated NSs. Kick distributions derived in this way overestimate the intrinsic NS kick distribution. For NSs in binaries, we can only directly estimate the effect of the natal kick on the binary system, instead of the natal kick received by the NS itself. Here, for the first time, we present a binary kick distribution for NSs with low-mass companions. We compile a catalogue of 145 NSs in low-mass binaries with the best available constraints on proper motion, distance, and systemic radial velocity. For each binary, we use a three-dimensional approach to estimate its binary kick. We discuss the implications of these kicks on system formation, and provide a parametric model for the overall binary kick distribution, for use in future theoretical modelling work. We compare our results with other work on isolated NSs and NSs in binaries, finding that the NS kick distributions fit using only isolated pulsars underestimate the fraction of NSs that receive low kicks. We discuss the implications of our results on modelling double NS systems, and provide suggestions on how to use our results in future theoretical works.

A spectroscopic modelling method for the detached eclipsing binaries to derive atmospheric parameters

Based on luminosity contributions, we develop aspectroscopic modelling method to derive atmospheric parameters of component stars in binary systems. The method is designed for those spectra of binaries that show double-lined features due to radial-velocity differences between the component stars. We first derive the orbital parameters and the stellar radii by solving the light and radial-velocity curves. The luminosity contributions in different phases can then be calculated. We construct a synthesised double-lined spectra model by superposing theoretical single-star spectra according to the luminosity contributions. Finally, we derive the atmospheric parameters of each component star using the model-fitting method. For multi-epoch double-lined spectra observed by the Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) Medium Resolution Survey (R ~ 7500), our method gives robust results for detached eclipsing binary systems observed in different orbital phases. Furthermore, this method can also be applied to other spectroscopic data with different resolutions as long as the systems are detached eclipsing binaries with nearly spherical stars.

Comparison of different models for gravitational and electromagnetic radiation of binary neutron star mergers

The neutron star binary merger is one of the most energetic phenomena in our Universe. Based on the calculation of gravitational radiation (GR) with the separate consideration of revolution and stellar rotation, and electromagnetic radiation (ER) with unipolar induction DC circuit model and magnetic dipole model, the results are compared to analyse the extent to which the stellar rotation will affect the total gravitational radiation power. Besides, the relationships between radiation power and the related parameters (e.g., orbital radius and stellar mass) are investigated. Furthermore, the feasibility of different types of binary star merger as the progenitor of fast radio bursts and gamma-ray bursts were studied based on the two models of ER, obtaining that binary neutron star and neutron star-white dwarf systems are among the possible progenitors. Finally, the radiation power of both GR and ER are compared under the same conditions. According to the results, the energy dissipation of the system is dominated by gravitational radiation. These results shed lights on further studies on the radiation processes of binary neutron star mergers and binary star systems.