X-ray Binaries Research Articles

Nebular emission from young stellar populations including binary stars

ABSTRACT We investigate the nebular emission produced by young stellar populations using the new galsevn model based on the combination of the sevn population-synthesis code including binary-star processes and the galaxev code for the spectral evolution of stellar populations. Photoionization calculations performed with the cloudy code confirm that accounting for binary-star processes strongly influences the predicted emission-line properties of young galaxies. In particular, we find that our model naturally reproduces the strong He ii λ4686/H β ratios commonly observed at high H β equivalent widths in metal-poor, actively star-forming galaxies, which have proven challenging to reproduce using previous models. Including bursty star formation histories broadens the agreement with observations, while the most extreme He ii λ1640 equivalent widths can be reproduced by models dominated by massive stars. galsevn also enables us to compute, for the first time in a way physically consistent with stellar emission, the emission from accretion discs of X-ray binaries (XRBs) and radiative shocks driven by stellar winds and supernova explosions. We find that these contributions are unlikely to prominently affect the predicted He ii λ4686/H β ratio, and that previous claims of a significant contribution by XRBs to the luminosities of high-ionization lines are based on models predicting improbably high ratios of X-ray luminosity to star formation rate, inconsistent with the observed average luminosity function of XRBs in nearby galaxies. The results presented here provide a solid basis for a more comprehensive investigation of the physical properties of observed galaxies with galsevn using Bayesian inference.

Unveiling hidden variability components in accreting X-ray binaries using both the Fourier power and cross-spectra

ABSTRACT We present a novel method for measuring the lags of (weak) variability components in neutron-star and black-hole low-mass X-ray binaries (LMXBs). For this we assume that the power and cross-spectra of these sources consists of a number of components that are coherent in different energy bands, but are incoherent with one another. The technique is based on fitting simultaneously the power spectrum (PS) and the Real and Imaginary parts of the cross-spectrum (CS) with a combination of Lorentzian functions. We show that, because the PS of LMXBs is insensitive to signals with a large Imaginary part and a small Real part in the CS, this approach allows us to uncover new variability components that are only detected in the CS. We also demonstrate that, contrary to earlier claims, the frequency of the type-C quasi-periodic oscillation (QPO) in the black-hole binary GRS 1915+105 does not depend on energy. Rather, the apparent energy dependence of the QPO frequency can be explained by the presence of a separate QPO component with a slightly higher frequency than that of the QPO, whose rms amplitude increases faster with energy than the rms amplitude of the QPO. From all the above we conclude that, as in the case of the PS, the CS of black-hole and neutron-star binaries can be fitted by a combination of Lorentzian components. Our findings provide evidence that the frequency-dependent part of the transfer function of these systems can be described by a combination of responses, each of them acting over relatively well-defined time-scales. This conclusion challenges models that assume that the main contribution to the lags comes from a global, broadband, transfer function of the accreting system.

Comprehensive X-ray and multiwavelength study of ULXs in NGC 1566

ABSTRACT This paper presents a comprehensive X-ray and multiwavelength study of ultraluminous X-ray sources (ULXs) in NGC 1566 using archival Chandra, Swift/XRT, JWST, and HST observations. The main results are, first, from the hardness ratio diagram, where spectral state transitions from hard to soft as seen in typical Galactic high-mass X-ray binaries for ULX-3 were observed. Secondly, a new transient ULX candidate (ULX-4) was identified, reaching a peak luminosity of ∼1040 erg s−1. Thirdly, the optical and NIR (near-infrared) counterparts of the ULXs were searched from the precise astrometric calculations. For ULX-1 and ULX-2, evidence was found that the observed NIR emission is due to the circumbinary disc/dust disrupted by X-rays. Lastly, the optical observations suggest that the possible donor star of ULX-3 is a B-type supergiant. In the case of ULX-4, the multiwavelength properties are not clear since many sources are detected within the astrometric error radius.

The influence of black holes on the binary population of the globular cluster Palomar 5

ABSTRACT The discovery of stellar-mass black holes (BHs) in globular clusters (GCs) raises the possibility of long-term retention of BHs within GCs. These BHs influence various astrophysical processes, including merger-driven gravitational waves and the formation of X-ray binaries. They also impact cluster dynamics by heating and creating low-density cores. Previous N-body models suggested that Palomar 5, a low-density GC with long tidal tails, may contain more than 100 BHs. To test this scenario, we conduct N-body simulations of Palomar 5 with primordial binaries to explore the influence of BHs on binary populations and the stellar mass function. Our results show that primordial binaries have minimal effect on the long-term evolution. In dense clusters with BHs, the fraction of wide binaries with periods >105 d decreases, and the disruption rate is independent of the initial period distribution. Multi-epoch spectroscopic observations of line-of-sight velocity changes can detect most bright binaries with periods below 104 d, significantly improving velocity dispersion measurements. Four BH-MS binaries in the model with BHs suggests their possible detection through the same observation method. Including primordial binaries leads to a flatter inferred mass function because of spatially unresolved binaries, leading to a better match of the observations than models without binaries, particularly in Palomar 5’s inner region. Future observations should focus on the cluster velocity dispersion and binaries with periods of 104–105 d in Palomar 5’s inner and tail regions to constrain BH existence.

A systematic study of the high-frequency bump in the black-hole low-mass X-ray binary GX 339 − 4

ABSTRACT The high-frequency bump, characterized by a frequency exceeding ∼30 Hz, represents a seldom-explored time-variability feature in the power density spectrum (PDS) of black-hole X-ray binaries. In the 2002, 2004, 2007, and 2010 outbursts of GX 339 − 4, the bump has been occasionally observed in conjunction with type-C quasi-periodic oscillations (QPOs). We systematically study the properties of the bump during these four outbursts observed by Rossi X-ray Timing Explorer (RXTE) in the 2–60 keV bands and detect the bump in 39 observations. While the frequencies of the type-C QPOs are in the range of ∼0.1–9 Hz, the root-mean-square (rms) amplitude of the bump shows an evolution in the hardness ratio versus the type-C QPO frequency plot. By comparing the rms amplitude of the bump with the corona temperature and simultaneous radio jet flux of the source, as previously studied in GRS 1915 + 105, we establish that in the hard state of GX 339 − 4, the bump is always strong, with the measurements of the rms amplitude in the range of 4–10 per cent. At the same time, the corona temperature is high and the radio flux is low. These findings indicate that, using the bump as a proxy, the majority of the accretion energy is directed towards the hot corona rather than being channeled into the radio jet. We discuss this phenomenon in terms of an inefficient energy transfer mechanism between the corona and jet in GX 339 − 4.

Synthetic X-ray emission from white dwarf accreting planetary material

ABSTRACT The emission of hard X-rays associated with white dwarfs (WDs) can be generated by the presence of a stellar companion, either by the companion’s coronal emission or by an accretion disc formed by the material stripped from the companion. Recent studies have suggested that a Jupiter-like planet can also be a donor of the material whose accretion on to the WD can generate hard X-rays. We use the guacho code to reproduce the conditions of this WD–planet scenario. With the example of the hard X-ray WD KPD 0005+5106, we explore different terminal wind velocities and mass-loss rates of a donor planet for a future network of simulations to investigate the luminosity and the spectral and temporal properties of the hard X-ray emission in the WD–planet systems. Our simulations show that the material stripped from the planet forms a disc and accretes on to the WD to reach temperatures high enough to generate hard X-rays as usually seen in X-ray binaries with low-mass companions. For high terminal wind velocities, the planet material does not form a disc, but it rather accretes directly on to the WD surface. The simulations reproduce the X-ray luminosity of another X-ray accreting WD (G 29-38), and only sometimes reaches the hard X-ray luminosity of KPD 0005+5106. The X-ray variability is stochastic and does not reproduce the period of KPD 0005+5106, suggesting that the additional physical processes (e.g. hot spots resulting from magnetic channelling of the accreting material) need to be explored.

Delayed Massive-star Mechanical Feedback at Low Metallicity

The classical model of massive-star mechanical feedback is based on effects at solar metallicity (Z ⊙), yet feedback parameters are very different at low metallicity. Metal-poor stellar winds are much weaker, and more massive supernova progenitors likely collapse directly to black holes without exploding. Thus, for ∼0.4 Z ⊙ we find reductions in the total integrated mechanical energy and momentum of ∼40% and 75%, respectively, compared to values classically expected at solar metallicity. But in particular, these changes effectively delay the onset of mechanical feedback until ages of ∼10 Myr. Feedback from high-mass X-ray binaries could slightly increase mechanical luminosity between ages 5 and 10 Myr, but it is stochastic and unlikely to be significant on this timescale. Stellar dynamical mechanisms remove most massive stars from clusters well before 10 Myr, which would further promote this effect; this process is exacerbated by gas retention implied by weak feedback. Delayed mechanical feedback implies that radiation feedback therefore dominates at early ages, which is consistent with the observed absence of superwinds in some extreme starbursts. This scenario may lead to higher star formation efficiencies, multiple stellar populations in clusters, and higher Lyman continuum escape. This could explain the giant star-forming complexes in metal-poor galaxies and the small sizes of OB superbubble shells relative to their inferred ages. It could also drive modest effects on galactic chemical evolution, including on oxygen abundances. Thus, delayed low-metallicity mechanical feedback may have broad implications, including for early cosmic epochs.

Variability as a Predictor for the Hard-to-soft State Transition in GX 339−4

During the outbursts of black hole X-ray binaries (BHXRBs), their accretion flows transition through several states. The source luminosity rises in the hard state, dominated by nonthermal emission, before transitioning to the blackbody-dominated soft state. As the luminosity decreases, the source transitions back into the hard state and fades to quiescence. This picture does not always hold, as ≈40% of the outbursts never leave the hard state. Identifying the physics that govern state transitions remains one of the outstanding open questions in black hole astrophysics. In this paper we present an analysis of archival RXTE data of multiple outbursts of GX 339−4. We compare the properties of the X-ray variability and time-averaged energy spectrum and demonstrate that the variability (quantified by the power spectral hue) systematically evolves ≈10–40 days ahead of the canonical state transition (quantified by a change in spectral hardness); no such evolution is found in hard-state-only outbursts. This indicates that the X-ray variability can be used to predict if and when the hard-to-soft state transition will occur. Finally, we find a similar behavior in 10 outbursts of four additional BHXRBs with more sparse observational coverage. Based on these findings, we suggest that state transitions in BHXRBs might be driven by a change in the turbulence in the outer regions of the disk, leading to a dramatic change in variability. This change is only seen in the spectrum days to weeks later, as the fluctuations propagate inwards toward the corona.

The Impact of Angular Momentum Loss on the Outcomes of Binary Mass Transfer

We use the rapid binary population synthesis code COMPAS to investigate commonly used prescriptions for the determination of mass transfer stability in close binaries and the orbital separations after stable mass transfer. The degree of orbital tightening during nonconservative mass transfer episodes is governed by the poorly constrained angular momentum carried away by the ejected material. Increased orbital tightening drives systems toward unstable mass transfer leading to a common envelope. We find that the fraction of interacting binaries that will undergo only stable mass transfer throughout their lives fluctuates between a few and ∼20% due to uncertainty in the angular momentum loss alone. If mass transfer is significantly nonconservative, stability prescriptions that rely on the assumption of conservative mass transfer underpredict the number of systems which experience unstable mass transfer and stellar mergers. This may substantially impact predictions about the rates of various transients, including luminous red novae, stripped-envelope supernovae, X-ray binaries, and the progenitors of coalescing compact binaries.

High-soft to Low-hard State Transition in Black Hole X-Ray Binaries with GRMHD Simulations

To understand the decaying phase of outbursts in the black hole (BH) X-ray binaries (BH-XRBs), we performed very long general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a Kerr BH with slowly rotating matter injected from outside. We thoroughly studied the flow properties, dynamical behavior of the accretion rate, magnetic flux rate, and jet properties during the temporal evolution. Due to the interaction between the thin disk and injected matter, the accretion flow near the BH goes through different phases. The sequence of phases is: soft state → soft-intermediate state → hard-intermediate state → hard state → quiescent state. For the accretion rate (and hence the luminosity) to decrease (as observed) in our model, the mass injection should not decay slower than the angular momentum injection. We also observed quasiperiodic oscillations (QPOs) in the accretion flow. Throughout the evolution, we observed low-frequency QPOs (∼10 Hz) and high-frequency QPOs (∼200 Hz). Our simple unified accretion flow model for state transitions is able to describe outbursts in BH-XRBs.

The Fundamental Plane of GRBs

Gamma-ray bursts (GRBs) exhibit powerful radiation and relativistic jets similar to blazars. However, the central engine of GRBs remains unknown. In this paper, we use the fundamental plane to analyze a sample of GRBs with measured mass. We extend, over ∼12 orders of magnitude, the correlation analysis and fundamental plane with a sample of X-ray binaries, active galactic nuclei (AGNs) including blazars, and brightest cluster galaxies. The fundamental plane of our de-beamed sample, with a measured mass ( = (0.60 ± 0.03)+(0.78 ± 0.02)+7.23 ± 0.95), closely aligns with the findings of previous work on AGNs. This finding suggests that GRBs adhere to the fundamental plane of AGNs and supports the theory proposed in previous work that the central engine of GRBs may be black holes. This observation provides a plausible explanation for the striking similarities between GRBs and AGNs.

Testing theories of accretion and gravity with super-extremal Kerr discs

ABSTRACT Fitting the thermal continuum emission of accreting black holes observed across X-ray bands represents one of the principle means of constraining the properties (mass and spin) of astrophysical black holes. Recent ‘continuum fitting’ studies of Galactic X-ray binaries in the soft state have found best-fitting dimensionless spin values which run into the prior bounds placed on traditional models (a⋆ = 0.9999). It is of critical importance that these results are robust, and not a result solely of the presence of these prior bounds and deficiencies in conventional models of accretion. Motivated by these results we derive and present superkerr, an xspec model comprising of a thin accretion disc solution valid in the Kerr geometry for arbitrary spin parameter a⋆, extending previous models valid only for black holes (|a⋆| < 1). This extension into ‘super-extremal’ space–times with |a⋆| > 1 includes solutions which describe discs evolving around naked singularities, not black holes. While being valid solutions of Einstein’s field equations these naked singularities are not expected to be present in nature. We discuss how the ‘measurement’ of a Kerr spin parameter 1 < a⋆ < 5/3 would present compelling evidence for the requirement of a rethink in either standard accretion theory, or our theories of gravity.

A relativistic outflow model of the X-ray polarization in Cyg X-1

ABSTRACT We propose that the polarization of the emission from the X-ray binary Cygnus X-1, measured using the Imaging X-ray Polarimetry Explorer, is imprinted by bulk Comptonization of coronal emission in a mildly relativistic wind or jet with a hollow-cone geometry. Models based on scattering in a static corona overlying a thin accretion disc have difficulty reproducing the relatively high-polarization degree (PD $\sim 4~{{\ \rm per\ cent}}$) concurrently with the low inclination (∼30°) of the binary orbit. We show that bulk outflow with a Lorentz factor ≳1.5 is adequate to reproduce the observed PD, with position angle parallel to the large-scale jet, provided that the scattering occurs in a conical sheath offset from the jet axis and our line of sight aligns roughly with the opening angle of the cone. Physically, this flow geometry could represent the entrainment of dense material near the base of an accelerating jet as it passes through the disc corona, or a slow (but still relativistic) sheath around a fast jet. If similar outflows are present in other X-ray binaries at higher inclination, we might expect to see still higher degrees of linear polarization ${\lesssim} 10~{{\ \rm per\ cent}}$ with an orientation perpendicular to the jet direction.

Population synthesis of Be X-ray binaries: metallicity dependence of total X-ray outputs

ABSTRACT X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionization histories during the Epoch of Reionization and Cosmic Dawn (z ∼ 5–30). Theoretical predictions of the X-ray emission from XRBs are important for modelling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion discs around rapidly rotating O/B stars, has not been investigated systematically. Be-XRBs are the largest class of high-mass XRBs (HMXBs) identified in local observations and are expected to play even more important roles in metal-poor environments at high redshifts. In light of this, we build a physically motivated model for Be-XRBs based on recent hydrodynamic simulations and observations of decretion discs. Our model is able to reproduce the observed population of Be-XRBs in the Small Magellanic Cloud with appropriate initial conditions and binary stellar evolution parameters. We derive the X-ray output from Be-XRBs as a function of metallicity in the (absolute) metallicity range Z ∈ [10−4, 0.03] with a large suite of binary population synthesis (BPS) simulations. The simulated Be-XRBs can explain a non-negligible fraction ($\gtrsim 30{{\ \rm per\ cent}}$) of the total X-ray output from HMXBs observed in nearby galaxies for Z ∼ 0.0003–0.02. The X-ray luminosity per unit star formation rate from Be-XRBs in our fiducial model increases by a factor of ∼8 from Z = 0.02 to Z = 0.0003, which is similar to the trend seen in observations of all types of HMXBs. We conclude that Be-XRBs are potentially important X-ray sources that deserve greater attention in BPS of XRBs.

4U 1210−64: a new member of the rare intermediate-mass X-ray binary subclass

ABSTRACT 4U 1210−64 is a peculiar X-ray binary that was first discovered in 1978 by the Uhuru satellite. The analysis of the X-ray data revealed a 6.7-d orbital period and an additional long-term modulation that is manifested as low- and high-flux states. Based on the previous classification of the donor star from the analysis of its optical spectra, the system has been suggested to be a high-mass X-ray binary. We revisit the optical classification where, based on the spectra from the Southern African Large Telescope (SALT), we conclude that the donor star is of spectral class A8 III–IV, making it a member of the rare intermediate-mass X-ray binaries. We perform radial velocity analysis using the SALT spectra where we consider circular and eccentric orbits. From the mass function derived and the mass constraints of the donor star, we demonstrate that a neutron star is favoured as the compact object in the binary system. We show, for the first time, the folded optical light curves, whose shape is interpreted to be due to a combination of ellipsoidal variations, irradiation of the donor star, and mutual eclipses of the star and accretion disc.

Infrared spectroscopy of SWIFT J0850.8−4219: identification of the second red supergiant X-ray binary in the Milky Way

ABSTRACT High-mass X-ray binaries hosting red supergiant (RSG) donors are a rare but crucial phase in massive stellar evolution, with only one source previously known in the Milky Way. In this letter, we present the identification of the second Galactic RSG X-ray binary SWIFT J0850.8−4219. We identify the source 2MASS 08504008−4211514 as the likely infrared counterpart with a chance coincidence probability ≈5 × 10−6. We present a $1.0{\!-\!}2.5\, \mu$m spectrum of the counterpart, exhibiting features characteristic of late-type stars and an exceptionally strong He i emission line, corroborating the identification. Based on i) the strength of the 12CO(2,0) band, ii) strong CN bandheads and absent TiO bandheads at ≈1.1 µm and iii) equivalent width of the Mg i$1.71\, \mu$m line, we classify the counterpart to be a K3–K5 type RSG with an effective temperature of 3820 ± 100 K, located at a distance of ≈12 kpc. We estimate the source X-ray luminosity to be (4 ± 1) × 1035 erg s−1, with a hard photon index (Γ < 1), arguing against a white dwarf accretor but consistent with a magnetized neutron star in the propeller phase. Our results highlight the potential of systematic near-infrared spectroscopy of Galactic hard X-ray sources in completing our census of the local X-ray binary population.

VLA monitoring of LS V +44 17 reveals scatter in the X-ray–radio correlation of Be/X-ray binaries

ABSTRACT LS V +44 17 is a persistent Be/X-ray binary (BeXRB) that displayed a bright, double-peaked period of X-ray activity in late 2022/early 2023. We present a radio monitoring campaign of this outburst using the Very Large Array. Radio emission was detected, but only during the second, X-ray brightest, peak, where the radio emission followed the rise and decay of the X-ray outburst. LS V +44 17 is therefore the third neutron star BeXRB with a radio counterpart. Similar to the other two systems (Swift J0243.6+6124 and 1A 0535+262), its X-ray and radio luminosity are correlated: we measure a power-law slope $\beta = 1.25^{+0.64}_{-0.30}$ and a radio luminosity of LR = (1.6 ± 0.2) × 1026 erg s−1 at a 0.5–10 keV X-ray luminosity of 2 × 1036 erg s−1 (i.e. $\sim 1~{{\ \rm per\ cent}}$LEdd). This correlation index is slightly steeper than measured for the other two sources, while its radio luminosity is higher. We discuss the origin of the radio emission, specifically in the context of jet launching. The enhanced radio brightness compared to the other two BeXRBs is the first evidence of scatter in the giant BeXRB outburst X-ray–radio correlation, similar to the scatter observed in subclasses of low-mass X-ray binaries. While a universal explanation for such scatter is not known, we explore several options: we conclude that the three sources do not follow proposed scalings between jet power and neutron star spin or magnetic field, and instead briefly explore the effects that ambient stellar wind density may have on BeXRB jet luminosity.

SRG/eROSITA catalogue of X-ray active SDSS dwarf galaxies

ABSTRACT We present a sample of 99 dwarf galaxies (M* < 109.5 M⊙) with X-ray activity in their central regions. The sample was obtained from a match of the SRG/eROSITA X-ray catalogue in the eastern galactic hemisphere with the MPA-JHU SDSS catalogue. The obtained matches were cleaned rigorously with the help of external optical catalogues to increase the purity of the sample. This work is the largest study of this kind – X-ray activity in ≈85 per cent of matched dwarfs was not reported before. The majority of X-ray active dwarfs are identified as star-forming galaxies. However, the X-ray luminosity of 82 objects cannot be explained by the collective emission of X-ray binaries, rendering them strong candidates for dwarf galaxies with an active accreting black hole in their centre. We find that the fraction of AGN among dwarf galaxies drops from ∼2 · 10−2 at LX ∼ 1039 erg s–1 to ∼(2 − 4) · 10−4 at LX ∼ 1041 erg s–1 and increases with the stellar mass of the host galaxy. We serendipitously discovered sources with unexpected properties. We report on a tidal disruption event (TDE) candidate in a dwarf galaxy, a massive black hole in a dwarf galaxy with a soft thermal spectrum, a luminous dwarf galaxy with an obscured X-ray spectrum and a few other peculiar sources. We found three Ultra-luminous X-ray (ULX) source candidates and a sample of X-ray bright galaxy pairs, in four of which both members shine in X-rays.

Insights into Density and Location Diagnostics of Photoionized Outflows in X-Ray Binaries

The population of metastable levels is key to high-precision density diagnostics of astrophysical plasmas. In photoionized plasmas, density is used to infer the distance from the ionizing source, which is otherwise difficult to obtain. Perfecting models that compute these populations is thus crucial. The present paper presents a semianalytic hydrogenic approximation for assessing the relative importance of different processes in populating atomic levels. This approximation shows that in the presence of a radiation source, photoexcitations and collisional excitations are both important over a wide range of plasma temperatures and ionizing spectra, while radiative recombination is orders of magnitude weaker. The interesting case of Fe+21 with a collisional radiative model with photoexcitation demonstrates this effect. The population of the first excited metastable level in Fe+21 is sensitive to the electron number density in the critical range of n e = 1012–1015 cm−3; it was observed to be significantly populated in the X-ray spectrum of the 2005 outburst of the X-ray binary GRO J1655-40. The present model shows that photoexcitation is the predominant process indirectly populating the metastable level. For the photoionized plasma in the GRO J1655-40 outflow, the model indicates a measured value of n e = (2.6 ± 0.5) × 1013 cm−3, implying a distance from the source of r = (4.4 ± 0.4) × 1010 cm. Finally, we show how the computed critical density and distance of Fe+21 yield the correct ionization parameter of the ion, independent of ionization balance calculations.

Discovery of X-Ray Polarization from the Black Hole Transient Swift J1727.8−1613

We report the first detection of the X-ray polarization of the bright transient Swift J1727.8−1613 with the Imaging X-ray Polarimetry Explorer. The observation was performed at the beginning of the 2023 discovery outburst, when the source resided in the bright hard state. We find a time- and energy-averaged polarization degree of 4.1% ± 0.2% and a polarization angle of 2.°2 ± 1.°3 (errors at 68% confidence level; this translates to ∼20σ significance of the polarization detection). This finding suggests that the hot corona emitting the bulk of the detected X-rays is elongated, rather than spherical. The X-ray polarization angle is consistent with that found in submillimeter wavelengths. Since the submillimeter polarization was found to be aligned with the jet direction in other X-ray binaries, this indicates that the corona is elongated orthogonal to the jet.