Neutron Star Low-mass X-ray Binaries Research Articles

High-Energy Spectra of Black Hole and Neutron Star Low-Mass X-Ray Binaries

In the case of low-mass X-ray binaries, the companion star is often too faint for detection; therefore, there is no chance for dynamical studies to independently determine the mass of the compact object. In the absence of a mass estimate, one cannot make a distinction as to whether the binary hosts a black hole or neutron star. Therefore, the question arises whether this distinction can be made based on the X-ray data alone, even when there are no bursts or pulsations. These would automatically imply a neutron star, but they are not always present. Black hole systems are known to emit radiation with an unbroken power–law shape up to several hundred keV energies in their high/soft states. If the non-thermal Comptonization processes that are responsible for this are somehow related to the lack of a solid surface, and to the fact that more gravitational potential energy can be released for a black hole, then there would be a definite method to reliably distinguish between the two sources. This work intends to review the available observations and studies to compare how these two populations behave during their different spectral states. A conclusion can be made that high/soft-state spectra are really different for black holes and neutron stars, for the low/hard state; however, the same conclusion cannot be safely made.

Orbital-period Changes of Low-mass X-Ray Binaries Driven by Magnetic Braking

Magnetic braking (MB) plays an important role in driving the evolution of low-mass X-ray binaries (LMXBs). The modified MB prescription, the convection and rotation boosted (CARB) model, is very successful in reproducing the detected mass-transfer rates of persistent neutron star (NS) LMXBs. In this work, we investigate whether the CARB MB prescription could account for the formation and evolution of some NS and black hole (BH) LMXBs with an observed orbital-period derivative. Using the MESA code, we perform a detailed binary evolution model for six NS and three BH LMXBs. Our simulations find that the CARB MB prescription can successfully reproduce the observed donor-star masses, orbital periods, and period derivatives of four NS LMXBs and one BH LMXB. Our calculated effective temperatures are in good agreement with the detected spectral types of two NS LMXBs and one BH LMXB. However, the standard MB model makes it difficult to produce the observed period derivatives of those LMXBs experiencing a rapid orbital shrinkage or expansion.

Probing the origin of the extended flaring branch of Z-type X-ray binaries GX 340+0 and GX 5-1 using AstroSat

Abstract ‘Z’ type neutron star low-mass X-ray binaries typically show a ‘Z’-like three-branched track in their hardness intensity diagram. However, a few such ‘Z’ sources show an additional branch known as the extended flaring branch (EFB). EFB has been poorly studied, and its origin is not known. It is thought to be an extension of the flaring branch (FB) or associated with Fe Kα complex or an additional continuum due to the radiative recombination continuum (RRC) process. Using AstroSat observations, we have detected the EFB from two ‘Z’ sources, GX 340+0 and GX 5–1, and performed a broadband spectral analysis in the 0.5–22 keV energy range. During EFB, both sources show the presence of a significant RRC component with absorption edges at $7.91^{+0.16}_{-0.15}$ keV and $8.10^{+0.16}_{-0.17}$ keV, respectively along with blackbody radiation and thermal Comptonisation. No signature of RRC was detected during the FB, which is adjoint to the EFB. No Fe Kα complex is detected. Interestingly, inside EFB dips of GX 5-1, for the first time, we have detected flaring events of 30–60s, which can be modelled with a single blackbody radiation. During the FB to EFB transition, an increase in the blackbody radius by a factor of 1.5–2 is observed in both sources. Our analysis strongly suggests that EFB is not an extension of FB or caused by the Fe Kα complex. Rather, it is caused by a sudden expansion of the hot, thermalised boundary layer and subsequent rapid cooling.

Radiation mechanism of twin kilohertz quasi-periodic oscillations in neutron star low-mass X-ray binaries

Context. The connection between quasi-periodic oscillations (QPOs) and magnetic fields has been investigated in various celestial bodies. Magnetohydrodynamic (MHD) waves have been employed to explain the simultaneous upper and lower kilohertz (kHz) QPOs. Nevertheless, the intricate and undefined formation pathways of twin kHz QPOs present a compelling avenue for exploration. This study area holds great interest as it provides an opportunity for deriving crucial parameters related to compact stars. Aims. We strive to develop a self-consistent model elucidating the radiation mechanism of twin kHz QPOs, which we then compare it with observations. Methods. A sample of 28 twin kHz QPOs detected from the X-ray binary 4U 1636–53 was used for a comparison with the results of the Markov chain Monte Carlo calculations based on our model of the radiation mechanism of twin kHz QPOs, which is related to twin MHD waves. Results. We obtained 28 groups of parameters of 4U 1636–53 and a tight exponential fit between the flux and the temperature of seed photons to Compton up-scattering and find that the electron temperature in the corona around the neutron star decreases with increasing seed photon temperature. Conclusions. The origin of twin kHz QPOs are dual disturbances that arise from twin MHD waves that are generated at the innermost radius of an accretion disc. The seed photons can be transported through a high temperature corona and are Compton up-scattered. The photons that vary the frequencies of twin MHD waves lead to the observed twin kHz QPOs.

Long-term evolution of spin and other properties of neutron star low-mass X-ray binaries: implications for millisecond X-ray pulsars

ABSTRACT A neutron star (NS) accreting matter from a companion star in a low-mass X-ray binary (LMXB) system can spin up to become a millisecond pulsar (MSP). Properties of many such MSP systems are known, which is excellent for probing fundamental aspects of NS physics when modelled using the theoretical computation of NS LMXB evolution. Here, we systematically compute the long-term evolution of NS, binary, and companion parameters for NS LMXBs using the stellar evolution code mesa. We consider the baryonic to gravitational mass conversion to calculate the NS mass evolution and show its cruciality for the realistic computation of some parameters. With computations using many combinations of parameter values, we find the general nature of the complex NS spin frequency ($\nu$) evolution, which depends on various parameters, including accretion rate, fractional mass-loss from the system, and companion star magnetic braking. Further, we utilize our results to precisely match some main observed parameters, such as $\nu$, orbital period ($P_{\rm orb}$), etc., of four accreting millisecond X-ray pulsars (AMXPs). By providing the $\nu$, $P_{\rm orb}$, and the companion mass spaces for NS LMXB evolution, we indicate the distribution and plausible evolution of a few other AMXPs. We also discuss the current challenges in explaining the parameters of AMXP sources with brown dwarf companions and indicate the importance of modelling the transient accretion in LMXBs as a possible solution.

Search for thermonuclear burst oscillations in the Swift/BAT data set

This study comprehensively analyzes type I X-ray bursts observed by Swift/BAT from 2005 to April 2024 to search for X-ray burst oscillations (XBOs) in neutron star low-mass X-ray binaries. XBOs, periodic signals detected within type I X-ray bursts, typically range from 11 to 620 Hz and are often observed in the soft X-ray data of these bursts. Using the high-sensitivity and precise timing capabilities of the Swift/BAT, we found 50 type I X-ray bursts from 37 neutron star low-mass X-ray binaries. We conducted a detailed timing analysis of these bursts. For sources with known burst oscillation frequencies, our findings largely corroborate previous studies. However, many sources displayed low confidence levels in the oscillation signals, with Z12 values between 10 and 20. For sources without known oscillation/spin frequencies, we utilized FFT analysis to search for signals across a broad frequency range. This approach revealed potential oscillation signals, with several bursts showing significance levels exceeding 3σ, including those from MAXI J1421–613, XTE J1701–407, XMM J174457–2850.3, Swift J1734.5–3027, IGR J17473–2721, Swift J174805.3–244637, Swift J181723.1–164300, and X 1832–330.

Accretion tori around rotating neutron stars II. Oscillations and precessions

The four characteristic oscillation frequencies of accretion flows (in addition to the Keplerian orbital frequency) are often discussed in the context of the time variability of black hole and neutron star (NS) low-mass X-ray binaries (LMXBs). These four frequencies are the frequencies of the axisymmetric radial and vertical epicyclic oscillations, and the frequencies of non-axisymmetric oscillations corresponding to the periastron (radial) and Lense-Thirring (vertical) precessions. In this context, we investigated the effect of the quadrupole moment of a slowly rotating NS and provide complete formulae for calculating these oscillation and precession frequencies, as well as convenient approximations. Simple formulae corresponding to the geodesic limit of a slender torus (and test-particle motion) and the limit of a marginally overflowing torus (a torus exhibiting a critical cusp) are presented, and more general approximate formulae are included to allow calculations for arbitrarily thick tori. We provide the Wolfram Mathematica code used for our calculations together with the C++ and PYTHON codes for calculating the frequencies. Our formulae can be used for various calculations regarding the astrophysical signatures of the NS super-dense matter equation of state. For instance, we demonstrate that even for a given fixed number of free parameters, a model that accounts for fluid flow precession matches the frequencies of twin-peak quasiperiodic oscillations observed in NS LMXBs better than a model that uses geodesic precession.

The Influence of the Magnetic Braking Laws on the Evolution of Persistent and Transient Low-mass X-Ray Binaries

Swift J1858.6−0814 (hereafter J1858) is a transient neutron star (NS) low-mass X-ray binary (LMXB). There is controversy regarding its donor mass derived from observations and theoretical calculations. In this paper, we adopt seven magnetic braking (MB) prescriptions suggested in the literature and different metallicity Z to simulate the evolution of the LMXB. Our results show that, employing the MB model proposed by A. Reiners & S. Mohanty (“rm12”), the convection- and rotation-boosted (“carb”) model, and the intermediate (“inter”) and convection-boosted (“cboost”) models in K. X. Van et al. can match (part of) the observational parameters of J1858 well. We then apply our method to other observed LMXBs and find that the “rm12” and “inter” MB laws are most promising in explaining transient LMXBs. In comparison, the simulations with the “cboost” and “carb” MB laws are more inclined to reproduce persistent LMXBs and ultracompact X-ray binaries, respectively. Our results, though subject to computational and/or observational bias, show that it is challenging to find a unified MB law that applies to the NS LMXB subpopulations simultaneously, indicating our lack of understanding of the true MB law. In addition, we explore the influence of various MB laws on the magnitude of the bifurcation periods in LMXBs.

The IXPE View of Neutron Star Low-Mass X-ray Binaries

Low-mass X-ray binaries hosting weakly magnetized neutron stars (NS-LMXBs) are among the brightest sources in the X-ray sky. Since 2021, the Imaging X-ray Polarimetry Explorer (IXPE) has provided new measurements of the X-ray polarization of these sources. IXPE observations have revealed that most NS-LMXBs are significantly polarized in the X-rays, providing unprecedented insight into the geometry of their accretion flow. In this review paper, we summarize the first results obtained by IXPE on NS-LMXBs, the emerging trends within each class of sources (atoll/Z), and possible physical interpretations.

On the Magnetic Braking Law in Black Hole Low-mass X-Ray Binaries

Magnetic braking (MB) plays an important role in the evolution of close low-mass X-ray binaries (LMXBs). It is also essential to the formation of ultracompact X-ray binaries (UCXBs). There have been lively investigations on the MB mechanism(s) in both single stars and close binaries including cataclysmic variables and neutron star (NS) LMXBs but with diverse conclusions. In this paper we explore the effect of MB on the black hole (BH) LMXB evolution. We combine binary population synthesis with detailed binary evolution to obtain the expected properties of the Galactic BH LMXB population. The simulated results are compared with the observational data including the BH mass, companion mass, companion temperature, orbital period, and mean accretion rate. Our results reveal that the MB laws with relatively low efficiency (i.e., RM12 and RVJ83) exhibit better agreement with observations, contrary to what was found for NS LMXBs. This raises the interesting question about whether MB really follows the same unified law in different types of binaries. We also predict that only a very small fraction (≲2.5%) of BH LMXBs can evolve to be UCXBs. This explains why no BH UCXB has been discovered so far.

A Study of Interstellar Medium in the Line of Sight of Transient Neutron Star Low-Mass X-ray Binary, MXB 1659-298, by Timing and Spectral Analysis

This work is dedicated to the study of interstellar medium (ISM) along the line of sight (LOS) of the transient low-mass X-ray binary, MXB 1659-298, capitalizing the high resolving power of XMM-Newton in the soft energy range. We emphasized the analysis of reflection grating spectrometer (RGS) data in the energy range 0.5–2.15 keV, suitable for the study of ISM. The paper includes an explanation of why, in the soft X-ray energy range, only two observations (out of seven) were deemed eligible for analysis. Three absorption lines associated with highly ionized Fe XX (1s22p2-2p2 (3p) 4d), Si XIV (1s2-1s2p), and Mg XI (1s2-1s6p) were identified in the observations, with IDs of 8620701(2001) and 748391601(2015). These new absorption lines and the absorption edge due to the neutral oxygen K edge seen in the spectra validate the multiphase structure of ISM. The predominance of interstellar medium over the ionized absorber is established along the direction of the source. The equivalent hydrogen column density measured is nearly equal to the galactic HI value derived previously. The small value of the ionic column density of Fe, Si, and Mg in the site of the high-temperature region resembles previous findings.

Probing outbursts of the transient neutron star low-mass X-ray binary Aql X-1 with NICER: a study of spectral evolution

ABSTRACT X-ray observations of neutron star (NS) low-mass X-ray binaries (LMXBs) are useful for probing the physical processes close to the NS and for constraining source parameters. Aql X-1 is a transient NS LMXB that frequently undergoes outbursts and provides an excellent opportunity to study source properties and accretion mechanisms in a strong-gravity regime over a wide range of accretion rates. In this work, we systematically investigate the spectral evolution of Aql X-1 usingNICER observations during the source outbursts in 2019 and 2020. The NICER observations cover the complete transition of the source from its canonical hard state to a soft state and back. The spectra extracted from most observations can be explained by a partially Comptonized accretion disc. We find that the system can be described by an accretion disc with an inner temperature of $\sim 0.7$ keV and a Comptonizing medium of thermal electrons at $\sim 2$ keV, while the photon index is strongly degenerate with the covering fraction of the medium. We also find evidence of Fe K$\alpha$ fluorescence emission in the spectra, indicating reprocessing of the Comptonized photons. We observe an absorption column density higher than the Galactic column density for most of the observations, indicating a significant local absorption. For some of the observations in the 2020 outburst, however, the local absorption is negligible.

On the lack of X-ray pulsation in most neutron star low-mass X-ray binaries

ABSTRACT We have investigated whether the lack of X-ray pulsations from most neutron star (NS) low-mass X-ray binaries (LMXBs) could be due to the extension of their inner disc to the NS surface. To estimate the inner disc radii, we have employed the model, recently proposed to account for the torque reversals of LMXBs. In this model, the inner disc radius depends on the spin period as well as the dipole moment and the mass inflow rate of the disc. Our model results indicate that most LMXBs have mass accretion rates above the minimum critical rates required for the inner disc to reach down to the NS surface and thereby quench the pulsed X-ray emission. For most sources X-ray pulsations are allowed when the period decreases below a certain critical value. For the same parameters, the model is also consistent with the observed X-ray luminosity ranges of the individual accreting millisecond X-ray pulsars (AMXPs). The paucity of AMXPs compared to the majority population of non-pulsing LMXBs is explained, as well as the fact that AMXPs are transient sources.

Simultaneous Chandra and HST observations of the quiescent neutron star low-mass X-ray binaries in 47 Tucanae

ABSTRACT We present simultaneous Chandra X-ray Observatory and Hubble Space Telescope observations of three certain (X5, X7, W37) and two likely (X4, W17) quiescent neutron star low-mass X-ray binaries (qLMXBs) in the globular cluster 47 Tuc. We study these systems in the X-ray, optical, and near-ultraviolet (NUV) using the simultaneous data and additional non-contemporaneous HST data. We have discovered a blue and variable NUV counterpart to W17. We have not securely identified the eclipsing qLMXB W37 in the optical or NUV. Deeper high-resolution imaging is needed to further investigate the faint NUV excess near the centre of the W37 error circle. We suggest that a previously identified optical astrometric match to X7 is likely the true counterpart. The H α emission and the location of the counterpart in the colour–magnitude diagram, indicate that the secondary is probably a non-degenerate, H-rich star. This is consistent with previous results from fitting X7’s X-ray spectrum. In X4, the simultaneous X-ray and optical behaviour supports the earlier suggestion that the X-ray variability is driven by changes in accretion rate. The X-ray eclipses in X5 coincide with minima in the optical/NUV light curves. Comparison of the 47 Tuc qLMXBs with the cataclysmic variables (CVs) in the cluster confirms that overall the qLMXBs have larger X-ray to optical flux ratios. Based on their optical/NUV colours, we conclude that the accretion discs in the qLMXBs are less prominent than in CVs. This makes the ratio of X-ray flux to excess blue-optical flux a powerful discriminator between CVs and qLMXBs.

A large population of neutron star low-mass X-ray binaries with long outburst recurrence time?

ABSTRACT Low-mass X-ray binaries (LMXBs) with neutron stars show quite different features that depend on the rate of mass transfer from the donor star. With a high transfer rate, the Z sources are in a persistent soft spectral state, and with a moderate transfer rate the transient Atoll sources have outburst cycles like the black hole X-ray binaries. The observations document very long outburst recurrence times for quite a number of sources. We follow with our computations the evolution of the accretion disc until the onset of the ionization instability. For sources with a low mass transfer rate, the accumulation of matter in the disc is essentially reduced due to the continuous evaporation of matter from the disc to the coronal flow. Different mass transfer rates result in nearly the same amount of matter accumulated for the outburst, which means that the outburst properties are similar for sources with short outburst cycles and sources with long outburst cycles, contrary to some expectations. Then, for systems with long recurrence time, less sources will be detected and the total population of LMXBs could be larger than it appears. This would relieve the apparent problem that the observed number of LMXBs as progenitors of millisecond pulsars (MSPs) is too small compared to the number of MSPs. Concerning the few quasi-persistent sources with year-long soft states, we argue that these states are not outbursts, but quasi-stationary hot states as in Z sources.

Discovery of a variable energy-dependent X-ray polarization in the accreting neutron star GX 5−1

We report on the coordinated observations of the neutron star low-mass X-ray binary (NS-LMXB) GX 5−1 in X-rays (IXPE, NICER, NuSTAR, and INTEGRAL), optical (REM and LCO), near-infrared (REM), mid-infrared (VLT VISIR), and radio (ATCA). This Z-source was observed by IXPE twice in March–April 2023 (Obs. 1 and 2). In the radio band the source was detected, but only upper limits to the linear polarization were obtained at a 3σ level of 6.1% at 5.5 GHz and 5.9% at 9 GHz in Obs. 1 and 12.5% at 5.5 GHz and 20% at 9 GHz in Obs. 2. The mid-IR, near-IR, and optical observations suggest the presence of a compact jet that peaks in the mid- or far-IR. The X-ray polarization degree was found to be 3.7%±0.4% (at 90% confidence level) during Obs. 1 when the source was in the horizontal branch of the Z-track and 1.8%±0.4% during Obs. 2 when the source was in the normal-flaring branch. These results confirm the variation in polarization degree as a function of the position of the source in the color-color diagram, as for previously observed Z-track sources (Cyg X-2 and XTE 1701−462). Evidence of a variation in the polarization angle of ∼20° with energy is found in both observations, likely related to the different, nonorthogonal polarization angles of the disk and Comptonization components, which peak at different energies.

Analysis of Bright Source Hardness Ratios in the 4 yr Insight-HXMT Galactic Plane Scanning Survey Catalog

We conduct a statistical analysis of the hardness ratio (HR) for bright sources in the 4 yr Galactic Plane Scanning Survey catalog of Insight-HXMT. Depending on the stable (variable) flux F s (F v) or spectrum S s(S v) of each source, the bright sources are classified into three groups: F v&S v, F v&S s, and F s&S s. Our study of the HR characteristics in different types of sources reveals that accretion-powered neutron star (NS) low-mass X-ray binaries (LMXBs) exhibit softer energy spectra than NS high-mass X-ray binaries (HMXBs), but harder energy spectra than black hole binaries in most cases. This difference is probably due to their different magnetic field strengths. Additionally, F v&S v LMXBs tend to be harder than F v&S s LMXBs below 7 keV, while the opposite is true for HMXBs. Our results suggest that LMXBs may dominate unclassified sources, and NS binaries are likely to be the primary type of X-ray binaries with ambiguous compact stars. By comparing the HR of transient sources in their outburst and low-flux states, it is found that the averaged HR of four sources in the two states are roughly comparable within uncertainties. We also investigate the spatial properties of the three groups and find that F v&S v sources are mainly located in the longitude of −20° < l < 9°, F v&S s sources cross the Galactic Plane, and F s&S s sources are predominantly concentrated in 19° < l < 42°. In addition, analyzing the HR spatial distributions shows the absorption of soft X-rays (primarily below 2 keV) in the Galactic Plane.

Insight-HXMT observations of thermonuclear X-ray bursts in 4U 1636−53

ABSTRACT We conducted an analysis of 45 bursts observed from 4U 1636−53. To investigate the mechanism behind the light-curve profiles and the impact of thermonuclear X-ray bursts on the accretion environment in accreting neutron star low-mass X-ray binaries. This analysis employed both light-curve and time-resolved spectroscopy methodologies, with data collected by the Insight-Hard X-ray Modulation Telescope instrument. We found that 30 bursts exhibited similar light-curve profiles and were predominantly in the hard state, and two photospheric radius expansion (PRE) bursts were in the soft state. The light curves of most bursts did not follow a single exponential decay but displayed a dual-exponential behaviour. The initial exponent had a duration of approximately 6 s. We utilized both the standard method and the ‘fa’ method to fit the burst spectra. The majority of the ‘fa’ values exceeded 1, indicating an enhancement of the persistent emission during the burst. Under the two Comptonization components assumption, we suggest that the scattering of burst photons by the inner corona may mainly contribute to the persistent emission enhancement. We also observed an inverse correlation between the maximum fa and the persistent emission flux in the non-PRE burst. This anticorrelation suggests that when the accretion rate is lower, there is a greater enhancement of persistent emission during the burst peak. The prediction based on Poynting–Robertson drag (P–R drag) aligns with this observed anticorrelation.

Reflecting on accretion in neutron star low-mass X-ray binaries

Reflecting on accretion in neutron star low-mass X-ray binaries

AstroSat View of the Neutron Star Low-mass X-Ray Binary GX 340+0

Understanding the spectral evolution along the “Z”-shaped track in the hardness–intensity diagram of Z sources, which are a class of luminous neutron star low-mass X-ray binaries, is crucial to probe accretion processes close to the neutron star. Here, we study the horizontal branch (HB) and the normal branch (NB) of the Z source GX 340+0 using AstroSat data. We find that the HB and the NB appear as two different types of X-ray intensity dips, which can appear in any sequence and with various depths. Our 0.8–25 keV spectra of dips and the hard apex can be modeled by the emissions from an accretion disk, a Comptonizing corona covering the inner disk, and the neutron star surface. We find that as the source moves onto the HB, the corona is replenished and energized by the disk and a reduced amount of disk matter reaches the neutron star surface. We also conclude that quasiperiodic oscillations during HB/NB are strongly associated with the corona and explain the evolution of strength and hard lag of this timing feature using the estimated coronal optical depth evolution.