Star Low-mass X-ray Binaries Research Articles

Kilohertz quasi-periodic oscillations as probes of the X-ray color-color diagram and neutron star accretion-disk structure for Z sources

Based on the detected kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low-mass X-ray binaries (NS-LMXBs), we investigate the evolution of the NS magnetosphere-disk structure along the Z track in the X-ray color-color diagram (CCD) for luminous Z sources, such as Cyg X-2, GX 5-1, GX 17+2, and Sco X-1. We find that the magnetosphere-disk radiusrinferred by kHz QPOs for all the sources shows a monotonically decreasing trend along the Z track from the horizontal branch (HB) to the normal branch (NB), implying that the dominated radiation components may dramatically change as the accretion disk moves toward the NS surface. In addition, the specific radius that corresponds to the HB/NB vertex is found to be aroundr ∼ 20 km, implying a potential characteristic position of transiting for the X-ray radiation mode. Furthermore, we find that the NBs that occur near the NS surface have a radius ofr ∼ 16−20 km, which is systematically smaller than those of HBs that have radii ofr ∼ 20−29 km. To interpret the relation between the CCD properties and the special magnetosphere-disk radii of Z sources, we suggest that the magnetic field lines corresponding to NB are “frozen-in” to the plasma, and move further inward with the shrinking of the NS magnetosphere-disk radius and pile up near the NS surface. They then form a strong magnetic field region aroundr ∼ 16−20 km, where the high magnetic energy density and high plasma mass density may dominate the radiation process in NB.

A MeerKAT survey of nearby nova-like cataclysmic variables

ABSTRACT We present the results of MeerKAT radio observations of 11 nearby nova-like cataclysmic variables (CVs). We have detected radio emission from IM Eri, RW Sex, V3885 Sgr, and V603 Aql. While RW Sex, V3885 Sgr, and V603 Aql had been previously detected, this is the first reported radio detection of IM Eri. Our observations have doubled the sample of non-magnetic CVs with sensitive radio data. We observe that at our radio detection limits, a specific optical luminosity ${\gtrsim}2.2\times 10^{18}\,$ erg s−1 Hz−1 (corresponding to MV ≲ 6.0) is required to produce a radio detection. We also observe that the X-ray and radio luminosities of our detected nova-like CVs are on an extension of the $L_X\propto L_R^{\sim 0.7}$ power law originally proposed for non-pulsating neutron star low-mass X-ray binaries. We find no other correlations between the radio emission and emission in other wavebands or any other system parameters for the existing sample of radio-detected non-magnetic CVs. We measure in-band (0.9–1.7 GHz) radio spectral indices that are consistent with reports from earlier work. Finally, we constructed broad spectral energy distributions for our sample from published multiwavelength data, and use them to place constraints on the mass transfer rates of these systems.

The advection-dominated accretion flow for the anticorrelation between the X-ray photon index and the X-ray luminosity in neutron star low-mass X-ray binaries

ABSTRACT Observationally, an anticorrelation between the X-ray photon index Γ (obtained by fitting the X-ray spectrum between 0.5 and 10 keV with a single power law) and the X-ray luminosity L0.5-10 keV, i.e. a softening of the X-ray spectrum with decreasing L0.5-10 keV, is found in neutron star low-mass X-ray binaries (NS-LMXBs) in the range of $L_{\rm 0.5\!-\!10\,keV}\sim 10^{34}\!-\!10^{36}\ \rm erg\ s^{-1}$. In this paper, we explain the observed anticorrelation between Γ and L0.5–10 keV within the framework of the self-similar solution of the advection-dominated accretion flow (ADAF) around a weakly magnetized NS. The ADAF model intrinsically predicts an anticorrelation between Γ and L0.5–10 keV. In the ADAF model, there is a key parameter, fth, which describes the fraction of the ADAF energy released at the surface of the NS as thermal emission to be scattered in the ADAF. We test the effect of fth on the anticorrelation between Γ and L0.5–10 keV. It is found that the value of fth can significantly affect the anticorrelation between Γ and L0.5–10 keV. Specifically, the anticorrelation between Γ and L0.5–10 keV becomes flatter with decreasing fth as taking fth = 0.1, 0.03, 0.01, 0.005, 0.003, and 0, respectively. By comparing with a sample of non-pulsating NS-LMXBs with well measured Γ and L0.5–10 keV, we find that indeed only a small value of 0.003 ≲ fth ≲ 0.1 is needed to match the observed anticorrelation between Γ and L0.5–10 keV. Finally, we argue that the small value of fth ≲ 0.1 derived in this paper further confirms our previous conclusion that the radiative efficiency of NSs with an ADAF accretion may not be as high as $\epsilon \sim {\dot{M} GM\over R_{*}}/{\dot{M} c^2}\sim 0.2$.

Silicon K-edge Dust Properties of Neutron Star Low-mass X-ray Binaries

AbstractThe dust properties of the line-of-sight materials in neutron star low-mass X-ray binaries (LMXBs) can be probed by X-ray observations and laboratory experiments. We use a Markov chain Monte Carlo (MCMC) method to conduct a spectral analysis of Chandra ACIS-S/HETG archival data of a sample of LMXBs, including GX 5-1 and GX 13+1. Our MCMC-based analysis puts constraints on the Si K-edge dust properties of the outflowing disk winds in this sample. Further X-ray observations of other LMXBs will help us better understand the grain features of dense outflows and accretion flows in neutron star binary systems.

The connection between the UV/optical and X-ray emission in the neutron star low-mass X-ray binary Aql X-1

ABSTRACT Accreting neutron stars and black holes in low-mass X-ray binaries (LMXBs) radiate across the electromagnetic spectrum. Linking the emission produced at different wavelengths can provide valuable information about the accretion process and any associated outflows. In this work, we study simultaneous X-ray and ultraviolet (UV)/optical observations of the neutron star LMXB Aql X-1, obtained with the Neil Gehrels Swift Observatory during its 2013, 2014, and 2016 accretion outbursts. We find that the UV/optical and X-ray emission are strongly correlated during all three outbursts. For the 2013 and 2014 episodes, which had the best Swift sampling, we find that the correlation between the UV/optical and X-ray fluxes is significantly steeper during the decay (soft state) of the outburst than during the rise (hard-to-soft state). We observe a UV/optical hysteresis behaviour that is likely linked to the commonly known X-ray spectral hysteresis pattern. For the decays of the three outbursts, we obtain a correlation index that cannot be directly explained by any single model. We suspect that this is a result of multiple emission processes contributing to the UV/optical emission, but we discuss alternative explanations. Based on these correlations, we discuss which mechanisms are likely dominating the UV/optical emission of Aql X-1.

A strange star scenario for the formation of isolated millisecond pulsars

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about 20% Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of ∼60 km s−1 can produce ∼6 × 103 isolated MSPs and birth rate of ∼6.6 × 10−7 yr−1 in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

The Comptonizing medium of the neutron star in 4U 1636 − 53 through its lower kilohertz quasi-periodic oscillations

ABSTRACT Inverse Compton scattering dominates the high-energy part of the spectra of neutron star (NS) low-mass X-ray binaries (LMXBs). It has been proposed that inverse Compton scattering also drives the radiative properties of kilohertz quasi-periodic oscillations (kHz QPOs). In this work, we construct a model that predicts the energy dependence of the rms amplitude and time lag of the kHz QPOs. Using this model, we fit the rms amplitude and time lag energy spectra of the lower kHz QPO in the NS LMXB 4U 1636 − 53 over 11 frequency intervals of the QPO and report three important findings: (i) A medium that extends 1–8 km above the NS surface is required to fit the data; this medium can be sustained by the balance between gravity and radiation pressure, without forcing any equilibrium condition. (ii) We predict a time delay between the oscillating NS temperature, due to feedback, and the oscillating electron temperature of the medium, which, with the help of phase resolved spectroscopy, can be used as a probe of the geometry and the feedback mechanism. (iii) We show that the observed variability as a function of QPO frequency is mainly driven by the oscillating electron temperature of the medium. This provides strong evidence that the Comptonizing medium in LMXBs significantly affects, if not completely drives, the radiative properties of the lower kHz QPOs regardless of the nature of the dynamical mechanism that produces the QPO frequencies.

Quasi-simultaneous radio and X-ray observations of Aql X-1 : probing low luminosities

ABSTRACT Aql X-1 is one of the best-studied neutron star low-mass X-ray binaries. It was previously targeted using quasi-simultaneous radio and X-ray observations during at least seven different accretion outbursts. Such observations allow us to probe the interplay between accretion inflow (X-ray) and jet outflow (radio). Thus far, these combined observations have only covered one order of magnitude in radio and X-ray luminosity range; this means that any potential radio–X-ray luminosity correlation, LR ∝ LXβ, is not well constrained (β ≈ 0.4–0.9, based on various studies) or understood. Here we present quasi-simultaneous Very Large Array and Swift-XRT observations of Aql X-1’s 2016 outburst, with which we probe one order of magnitude fainter in radio and X-ray luminosity compared to previous studies (6 × 1034 erg s−1 < LX <3 × 1035 erg s−1, i.e. the intermediate to low-luminosity regime between outburst peak and quiescence). The resulting radio non-detections indicate that Aql X-1’s radio emission decays more rapidly at low X-ray luminosities than previously assumed – at least during the 2016 outburst. Assuming similar behaviour between outbursts, and combining all available data in the hard X-ray state, this can be modelled as a steep β =$1.17^{+0.30}_{-0.21}$ power-law index or as a sharp radio cut-off at LX ≲ 5 × 1035 erg s−1 (given our deep radio upper limits at X-ray luminosities below this value). We discuss these results in the context of other similar studies.

Neutron star QPOs from oscillating, precessing hot, thick flow

ABSTRACT Across black hole (BH) and neutron star (NS) low-mass X-ray binaries (LMXBs), there appears to be some correlation between certain high- and low-frequency quasi-periodic oscillations (QPOs). In a previous paper, we showed that for BH LMXBs, this could be explained by the simultaneous oscillation and precession of a hot, thick, torus-like corona. In the current work, we extend this idea to NS LMXBs by associating the horizontal branch oscillations (HBOs) with precession and the upper-kiloHertz (ukHz) QPO with vertical epicyclic motion. For the Atoll source 4U 1608-52, the model can match many distinct, simultaneous observations of the HBO and ukHz QPO by varying the inner and outer radius of the torus, while maintaining fixed values for the mass (MNS) and spin (a*) of the NS. The best-fitting values are MNS = 1.38 ± 0.03 M⊙ and a* = 0.325 ± 0.005. By combining these constraints with the measured spin frequency, we are able to obtain an estimate for the moment of inertia of INS = 1.40 ± 0.02 × 1045 g cm2, which places constraints on the equation of state. The model is unable to fit the lower-kHz QPO, but evidence suggests that QPO may be associated with the boundary layer between the accretion flow and the NS surface, which is not treated in this work.

R-modes in stratified neutron stars with entrainment

We calculate the temperature-dependent r-mode spectrum of superfluid neutron stars with npeµ (neutron, proton, electron, muon) core composition. This study is an extension of the previous work by Kantor, Gusakov [MNRAS 469, 3928 (2017)], where such spectrum was calculated under simplifying assumption of vanishing entrainment between superfluid neutrons and superconducting protons. We show that accounting for the entrainment leads to non-analytic behavior of the spectrum at small rotation rates. Namely, we find that in the leading order in rotation accounting for any non-zero value of entrainment eliminates superfluid r-modes. We show that next-to-leading order in rotation restores the superfluid r-modes in the spectrum. We calculate this spectrum and show that for certain neutron star models normal r-mode experiences stabilizing resonances with superfluid r-modes. This confirms the scenario of Gusakov, Chugunov, and Kantor [PRL 112, 151101 (2014)] that explains neutron stars in low-mass X-ray binaries.

Crust of accreting neutron stars within simplified reaction network

ABSTRACT Transiently accreting neutron stars in low-mass X-ray binaries are generally believed to be heated up by nuclear reactions in accreted matter during hydrostatic compression. Detailed modelling of these reactions is required for the correct interpretation of observations. In this paper, we construct a simplified reaction network, which can be easily implemented and depends mainly on atomic mass tables as nuclear physics input. We show that it reproduces results of the detailed network by Lau et al. very well, if one applies the same mass model. However, the composition and the heating power are shown to be sensitive to the mass table used and treatment of mass tables boundary, if one applies several of them in one simulation. In particular, the impurity parameter Qimp at density ρ = 2 × 1012 g cm−3 can differ for a factor of few, and even increase with density increase. The profile of integrated heat release shown to be well confined between results by Fantina et al. and Lau et al.

Insight-HXMT observation on 4U 1608–52: Evolving spectral properties of a bright type-I X-ray burst

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by ∼ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

On the frequency correlations of low-frequency QPOs with kilohertz QPOs in accreting millisecond X-ray pulsars

ABSTRACT We investigate frequency correlations of low frequency (LF, <80 Hz) and kHz quasi-periodic oscillations (QPOs) using the complete RXTE data sets on six accreting millisecond X-ray pulsars (AMXPs) and compare them to those of non-pulsating neutron star (NS) low-mass X-ray binaries with known spin. For the AMXPs SAX J1808.4−3658 and XTE J1807−294, we find frequency-correlation power-law indices that, surprisingly, are significantly lower than in the non-pulsars, and consistent with the relativistic precession model (RPM) prediction of 2.0 appropriate to test-particle orbital and Lense–Thirring precession frequencies. As previously reported, power-law normalizations are significantly higher in these AMXPs than in the non-pulsating sources, leading to requirements on the NS specific moment of inertia in this model that cannot be satisfied with realistic equations of state. At least two other AMXPs show frequency correlations inconsistent with those of SAX J1808.4−3658 and XTE J1807−294, and possibly similar to those of the non-pulsating sources; for two AMXPs no conclusions could be drawn. We discuss these results in the context of a model that has had success in black hole (BH) systems involving a torus-like hot inner flow precessing due to (prograde) frame dragging, and a scenario in which additional (retrograde) magnetic and classical precession torques not present in BH systems are also considered. We show that a combination of these interpretations may accommodate our results.

Long-term temperature evolution of neutron stars undergoing episodic accretion outbursts

Context. Transient neutron star low-mass X-ray binaries undergo episodes of accretion, alternated with quiescent periods. During an accretion outburst, the neutron star heats up due to exothermic accretion-induced processes taking place in the crust. Besides the long-known deep crustal heating of nuclear origin, a likely non-nuclear source of heat, dubbed “shallow heating”, is present at lower densities. Most of the accretion-induced heat slowly diffuses into the core on a timescale of years. Over many outburst cycles, a state of equilibrium is reached when the core temperature is high enough that the heating and cooling (photon and neutrino emission) processes are in balance.Aims. We investigate how stellar characteristics and outburst properties affect the long-term temperature evolution of a transiently accreting neutron star. For the first time the effects of crustal properties are considered, particularly that of shallow heating.Methods. Using our codeNSCool, we tracked the thermal evolution of a neutron star undergoing outbursts over a period of 105yr. The outburst sequence is based on the regular outbursts observed from the neutron star transient Aql X-1. For each model we calculated the timescale over which equilibrium was reached and we present these timescales along with the temperature and luminosity parameters of the equilibrium state.Results. We performed several simulations with scaled outburst accretion rates, to vary the amount of heating over the outburst cycles. The results of these models show that the equilibrium core temperature follows a logarithmic decay function with the equilibrium timescale. Secondly, we find that shallow heating significantly contributes to the equilibrium state. Increasing its strength raises the equilibrium core temperature. We find that if deep crustal heating is replaced by shallow heating alone, the core would still heat up, reaching only a 2% lower equilibrium core temperature. Deep crustal heating may therefore not be vital to the heating of the core. Additionally, shallow heating can increase the quiescent luminosity to values higher than previously expected. The thermal conductivity in the envelope and crust, including the potentially low-conductivity pasta layer at the bottom of the crust, is unable to significantly alter the long-term internal temperature evolution. Stellar compactness and nucleon pairing in the core change the specific heat and the total neutrino emission rate as a function of temperature, with the consequences for the properties of the equilibrium state depending on the exact details of the assumed pairing models. The presence of direct Urca emission leads to the lowest equilibrium core temperature and the shortest equilibrium timescale.

Population synthesis of accreting neutron stars emitting gravitational waves

ABSTRACT The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable rmodes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains, and unstable r modes. However, based on the resultant distributions alone, it is difficult to distinguish between the competing mechanisms.

A model for the radio/X-ray correlation in three neutron star low-mass X-ray binaries 4U 1728-34, Aql X-1ad, and EXO 1745-248

Observationally, for neutron star low-mass X-ray binaries, so far, the correlation between the radio luminosity $L_{\rm R}$ and the X-ray luminosity $L_{\rm X}$, i.e., $L_{\rm R}\propto L_{\rm X}^{\beta}$, has been reasonably well-established only in three sources 4U 1728-34, Aql X-1 and EXO 1745-248 in their hard state. The slope $\beta$ of the radio/X-ray correlation of the three sources is different, i.e., $\beta \sim 1.4$ for 4U 1728-34, $\beta \sim 0.4$ for Aql X-1, and $\beta \sim 1.6$ for EXO 1745-248. In this paper, for the first time we explain the different radio/X-ray correlation of 4U 1728-34, Aql X-1 and EXO 1745-248 with the coupled advection-dominated accretion (ADAF)-jet model respectively. We calculate the emergent spectrum of the ADAF-jet model for $L_{\rm X}$ and $L_{\rm R}$ at different $\dot m$ ($\dot m=\dot M/\dot M_{\rm Edd}$), adjusting $\eta$ ($\eta \equiv \dot M_{\rm jet}/\dot M$, describing the fraction of the accreted matter in the ADAF transfered vertically forming the jet) to fit the observed radio/X-ray correlations. Then we derive a fitting formula of $\eta$ as a function of $\dot m$ for 4U 1728-34, Aql X-1 and EXO 1745-248 respectively. If the relation between $\eta$ and $\dot m$ can be extrapolated down to a lower value of $\dot m$, we find that in a wide range of $\dot m$, the value of $\eta$ in Aql X-1 is greater than that of in 4U 1728-34 and EXO 1745-248, implying that Aql X-1 may have a relatively stronger large-scale magnetic field, which is supported by the discovery of the coherent millisecond X-ray pulsation in Aql X-1.

R-mode constraints from neutron star equation of state

The gravitational radiation has been proposed a long time before, as an explana- tion for the observed relatively low spin frequencies of young neutron stars and of accreting neutron stars in low-mass X-ray binaries as well. In the present work we studied the effects of the neutron star equation of state on the r-mode instability window of rotating neutron stars.

NuSTAR Observations of the Accreting Atolls GX 3+1, 4U 1702-429, 4U 0614+091, and 4U 1746-371

Abstract Atoll sources are accreting neutron star (NS) low-mass X-ray binaries. We present a spectral analysis of four persistent atoll sources (GX 3+1, 4U 1702−429, 4U 0614+091, and 4U 1746−371) observed for ∼20 ks each with NuSTAR to determine the extent of the inner accretion disk. These sources range from an apparent luminosity of 0.006–0.11 of the Eddington limit (assuming the empirical limit of 3.8 × 1038 erg s−1). Broad Fe emission features shaped by Doppler and relativistic effects close to the NS were firmly detected in three of these sources. The position of the disk appears to be close to the innermost stable circular orbit (ISCO) in each case. For GX 3+1, we determine (90% confidence level) and an inclination of 27°–31°. For 4U 1702−429, we find a and inclination of 53°–64°. For 4U 0614+091, the disk has a position of and inclination of 50°–62°. If the disk does not extend to the innermost stable circular orbit, we can place conservative limits on the magnetic field strength in these systems in the event that the disk is truncated at the Alfvén radius. This provides the limit at the poles of B ≤ 6.7 × 108 G, 3.3 × 108 G, and 14.5 × 108 G for GX 3+1, 4U 1702−429, and 4U 0614+091, respectively. For 4U 1746−371, we argue that the most plausible explanation for the lack of reflection features is a combination of source geometry and strong Comptonization. We place these sources among the larger sample of NSs that have been observed with NuSTAR.

Analytical model of compact star in low-mass X-ray binary with de Sitter spacetime

In this article, we suggest a Sitter model in favor of compact stars in low-mass X-ray binaries. Here, we have considered the presence of a cosmological constant (on a small scale) to investigate the stellar structure. We conclude that this approach is very suitable for the familiar physicalmodel of compact stars in low-mass X-ray binaries. We calculate the probable radius, compactness (u) and surface redshift (Zs) of six compact stars in low-mass X-ray binaries, namely Cyg X-2, V395 Carinae/2S 0921−630, XTE J2123−058, X1822−371 (V691 CrA), 4U 1820−30 and GR Mus (XB 1254−690). We also offer a possible equation of state (EOS) for the stellar objects.

A NuSTAR Observation of the Low-mass X-Ray Binary GX 349+2 throughout the Z-track

Abstract Although the most luminous class of neutron star (NS) low-mass X-ray binaries, known as Z sources, have been well studied, their behavior is not fully understood. In particular, what causes these sources to trace out the characteristic Z-shaped pattern on color–color or hardness–intensity diagrams (HIDs) is not well known. By studying the physical properties of the different spectral states of these sources, we may better understand such variability. With that goal in mind, we present a recent NuSTAR observation of the Z source GX 349+2, which spans approximately 2 days and covers all its spectral states. By creating an HID we were able to extract five spectra and trace the change in spectral parameters throughout the Z-track. GX 349+2 shows a strong, broad Fe Kα line in all states, regardless of the continuum model used. Through modeling of the reflection spectrum and Fe Kα line we find that in most states the inner disk radius is consistent with remaining unchanged at an average radius of 17.5 R g or 36.4 km for a canonical 1.4 M ⊙ NS. During the brightest flaring branch, however, the inner disk radius from reflection is not well constrained.