Black Hole X-ray Binary Systems Research Articles

Discovery of persistent quasi-periodic oscillations in accreting white dwarfs: a new link to X-ray binaries

ABSTRACT Almost all accreting black hole and neutron star (NS) X-ray binary systems (XRBs) exhibit prominent brightness variations on a few characteristic time-scales and their harmonics. These quasi-periodic oscillations (QPOs) are thought to be associated with the precession of a warped accretion disc, but the physical mechanism that generates the precessing warp remains uncertain. Relativistic frame dragging (Lense–Thirring precession) is one promising candidate, but a misaligned magnetic field is an alternative, especially for NS XRBs. Here, we report the discovery of five accreting white dwarf systems (AWDs) that display strong optical QPOs with characteristic frequencies and harmonic structures that suggest they are the counterpart of the QPOs seen in XRBs. Since AWDs are firmly in the classical (non-relativistic) regime, Lense–Thirring precession cannot account for these QPOs. By contrast, a weak magnetic field associated with the white dwarf can drive disc warping and precession in these systems, similar to what has been proposed for NS XRBs. Our observations confirm that magnetically driven warping is a viable mechanism for generating QPOs in disc-accreting astrophysical systems, certainly in AWDs and possibly also in NS XRBs. Additionally, they establish a new way to estimate magnetic field strengths, even in relatively weak-field systems where other methods are not available.

Implications of the Spin-Induced Accretion Disk Truncation on the X-ray Binary Broadband Emission

Black hole X-ray binary systems consist of a black hole accreting mass from its binary companion, forming an accretion disk. As a result, twin relativistic plasma ejections (jets) are launched towards opposite and perpendicular directions. Moreover, multiple broadband emission observations from X-ray binary systems range from radio to high-energy gamma rays. The emission mechanisms exhibit thermal origins from the disk, stellar companion, and non-thermal jet-related components (i.e., synchrotron emission, inverse comptonization of less energetic photons, etc.). In many attempts at fitting the emitted spectra, a static black hole is often assumed regarding the accretion disk modeling, ignoring the Kerr metric properties that significantly impact the geometry around the usually rotating black hole. In this work, we study the possible implications of the spin inclusion in predictions of the X-ray binary spectrum. We mainly focus on the most significant aspect inserted by the Kerr geometry, the innermost stable circular orbit radius dictating the disk’s inner boundary. The outcome suggests a higher-peaked and hardened X-ray spectrum from the accretion disk and a substantial increase in the inverse Compton component of disk-originated photons. Jet-photon absorption is also heavily affected at higher energy regimes dominated by hadron-induced emission mechanisms. Nevertheless, a complete investigation requires the full examination of the spin contribution and the resulting relativistic effects beyond the disk truncation.

Galactic Stellar Black Hole Binaries: Spin Effects on Jet Emissions of High-Energy Gamma-Rays

In the last few decades, galactic stellar black hole X-ray binary systems (BHXRBs) have aroused intense observational and theoretical research efforts specifically focusing on their multi-messenger emissions (radio waves, X-rays, γ-rays, neutrinos, etc.). In this work, we investigate jet emissions of high-energy neutrinos and gamma-rays created through several hadronic and leptonic processes taking place within the jets. We pay special attention to the effect of the black hole’s spin (Kerr black holes) on the differential fluxes of photons originating from synchrotron emission and inverse Compton scattering and specifically on their absorption due to the accretion disk’s black-body radiation. The black hole’s spin (dimensionless spin parameter a*) enters into the calculations through the radius of the innermost circular orbit around the black hole, the RISCO parameter, assumed to be the inner radius of the accretion disk, which determines its optical depth τdisk. In our results, the differential photon fluxes after the absorption effect are depicted as a function of the photon energy in the range 1GeV ≤E≤103GeV. It is worth noting that when the black holes’ spin (α*) increases, the differential photon flux becomes significantly lower.

NICER, NuSTAR, and Insight-HXMT Views to the Newly Discovered Black Hole X-Ray Binary Swift J1727.8-1613

Swift J1727.8–1613 is a black hole X-ray binary newly discovered in 2023. We perform spectral analysis with simultaneous Insight-HXMT, NICER, and NuSTAR observations when the source was approaching the hard intermediate state. Such a joint view reveals an additional hard component apart from the normally observed hard component with reflection in the spectrum, to be distinguished from the usual black hole X-ray binary systems. By including this extra component in the spectrum, we have measured a high spin of 0.98−0.07+0.02 and an inclination of around 40−0.8+1.2 °, which is consistent with NICER results reported before. However, we find that the additional spectral component cannot be exclusively determined due to the model degeneracy. Accordingly, a possible jet/corona configuration is adjusted to account for the spectral fitting with different model trials. The extra component may originate either from a relativistic jet or a jet base/corona underneath a slow jet.

Truncated, tilted discs as a possible source of Quasi-Periodic Oscillations

ABSTRACT Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency (≲10 Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs (≳60 Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with the black hole’s spin axis. However, existing versions of this model have not taken into account the effects of a surrounding, geometrically thin disc on the precessing, inner, geometrically thick flow. In Bollimpalli et. al 2023, using a set of GRMHD simulations of tilted, truncated accretion discs, we confirmed that the outer thin disc slows down the precession rate of the precessing torus, which has direct observational implications for type-C QPOs. In this paper, we provide a detailed analysis of those simulations and compare them with an aligned truncated disc simulation. We find that the misalignment of the disc excites additional variability in the inner hot flow, which is absent in the comparable aligned-disc simulations. This suggests that the misalignment may be a crucial requirement for producing QPOs. We attribute this variability to global vertical oscillations of the inner torus at epicyclic frequencies corresponding to the transition radius. This explanation is consistent with current observations of higher frequency QPOs in black hole X-ray binary systems.

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.

The effect of thermal winds on the outbursts evolution of LMXB systems

ABSTRACT Theoretical models of accretion discs and observational data indicate that the X-ray emission from the inner parts of an accretion disc can irradiate its outer regions and induce a thermal wind, which carries away the mass and angular momentum from the disc. Our aim is to investigate the influence of the thermal wind on the outburst light curves of black hole X-ray binary systems. We carry out numerical simulations of a non-stationary disc accretion with wind using upgraded open code freddi. We assume that the wind launches only from the ionized part of the disc and may turn off if the latter shrinks fast enough. Our estimates of the viscosity parameter α are shifted downward compared to a scenario without a wind. Generally, correction of α depends on the spectral hardness of central X-rays and the disc outer radius, but unlikely to exceed a factor of 10 in the case of a black hole low-mass X-ray binary (BH LMXB). We fit 2002 outburst of BH LMXB 4U 1543 − 47 taking into account the thermal wind. The mass-loss in the thermal wind is of order of the accretion rate on the central object at the peak of the outburst. New estimate of the viscosity parameter α for the accretion disc in this system is about two times lower than the previous one. Additionally, we calculate evolution of the number of hydrogen atoms towards 4U 1543 − 47 due to the thermal wind from the hot disc.

Probing the jet size of two black hole X-ray binaries in the hard state

ABSTRACT Using multifrequency very long baseline interferometry (VLBI) observations, we probe the jet size in the optically thick hard state jets of two black hole X-ray binary (BHXRB) systems, MAXI J1820+070 and V404 Cygni. Because of optical depth effects, the phase-referenced VLBI core positions move along the jet axis of the BHXRB in a frequency-dependent manner. We use this ‘core shift’ to constrain the physical size of the hard state jet. We place an upper limit of 0.3 au on the jet size measured between the 15 and 5 GHz emission regions of the jet in MAXI J1820+070, and an upper limit of 1.0 au between the 8.4 and 4.8 GHz emission regions of V404 Cygni. Our limit on the jet size in MAXI J1820+070 observed in the low-hard state is a factor of 5 smaller than the values previously observed in the high-luminosity hard state (using time lags between multifrequency light curves), thus showing evidence of the BHXRB jet scaling in size with jet luminosity. We also investigate whether motion of the radio-emitting region along the jet axis could affect the measured VLBI parallaxes for the two systems, leading to a mild tension with the parallax measurements of Gaia. Having mitigated the impact of any motion along the jet axis in the measured astrometry, we find the previous VLBI parallax measurements of MAXI J1820+070 and V404 Cygni to be unaffected by jet motion. With a total time baseline of 8 yr, due to having incorporated 14 new epochs in addition to the previously published ones, our updated parallax measurement of V404 Cygni is 0.450 ± 0.018 mas (2.226 ± 0.091 kpc).

Extending the theory of propagating fluctuations: the first fully relativistic treatment and analytical Fourier–Green’s functions

ABSTRACT The aperiodic variability ubiquitously observed from accreting black hole X-ray binary systems is generally analysed within the framework of the so-called ‘theory of propagating fluctuations’. In this paper we derive the Fourier transforms of the Green’s function solutions of the thin disc equations. These solutions suffice to describe all possible solutions through standard convolution techniques. Solutions are found for both Newtonian discs and general relativistic solutions with a vanishing ISCO stress. We use this new relativistic theory to highlight the Kerr black hole spin dependence of a number of observable variability properties of black hole discs. The phase lags, coherence, and power density spectra of Kerr discs are shown to be strong functions of black hole spin. Observations of the aperiodic variability of black hole accretion sources may now, at least in principle, offer a new avenue to directly constrain black hole spins.

Polarized x-rays constrain the disk-jet geometry in the black hole x-ray binary Cygnus X-1.

A black hole x-ray binary (XRB) system forms when gas is stripped from a normal star and accretes onto a black hole, which heats the gas sufficiently to emit x-rays. We report a polarimetric observation of the XRB Cygnus X-1 using the Imaging X-ray Polarimetry Explorer. The electric field position angle aligns with the outflowing jet, indicating that the jet is launched from the inner x-ray-emitting region. The polarization degree is 4.01 ± 0.20% at 2 to 8 kiloelectronvolts, implying that the accretion disk is viewed closer to edge-on than the binary orbit. These observations reveal that hot x-ray-emitting plasma is spatially extended in a plane perpendicular to, not parallel to, the jet axis.

Spin and Accretion Rate Dependence of Black Hole X-Ray Spectra

Abstract We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first-principles physical calculations. The power-law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe Kα equivalent width grows sublinearly with Fe abundance. Strikingly, the Kα profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30%–100% greater than as predicted by the Novikov–Thorne model.

The inner flow geometry in MAXI J1820+070 during hard and hard-intermediate states

Context. We present a systematic X-ray spectral-timing study of the recently discovered, exceptionally bright black hole X-ray binary system MAXI J1820+070. Our analysis focuses on the first part of the 2018 outburst, covering the rise throughout the hard state, the bright hard and hard-intermediate states, and the transition to the soft-intermediate state. Aims. We address the issue of constraining the geometry of the innermost accretion flow and its evolution throughout an outburst. Methods. We employed two independent X-ray spectral-timing methods applied to archival NICER data of MAXI J1820+070. We first identified and tracked the evolution of a characteristic frequency of soft X-ray thermal reverberation lags (lags of the thermally reprocessed disc emission after the irradiation of variable hard X-ray photons). This frequency is sensitive to intrinsic changes in the relative distance between the X-ray source and the disc. Then, we studied the spectral evolution of the quasi-thermal component responsible for the observed thermal reverberation lags. We did so by analysing high-frequency covariance spectra, which single out spectral components that vary in a linearly correlated way on the shortest sampled timescales and are thus produced in the innermost regions of the accretion flow. Results. The frequency of thermal reverberation lags steadily increases throughout most of the outburst, implying that the relative distance between the X-ray source and the disc decreases as the source softens. However, near transition this evolution breaks, showing a sudden increase (decrease) in lag amplitude (frequency). On the other hand, the temperature of the quasi-thermal component in covariance spectra, due to disc irradiation and responsible for the observed soft reverberation lags, consistently increases throughout all the analysed observations. Conclusions. This study proposes an alternative interpretation to the recently proposed contracting corona scenario. Assuming a constant height for the X-ray source, the steady increase in the reverberation lag frequency and in the irradiated disc temperature in high-frequency covariance spectra can be explained in terms of a decrease in the disc inner radius as the source softens. The behaviour of thermal reverberation lags near transition might be related to the relativistic plasma ejections detected at radio wavelengths, suggesting a causal connection between the two phenomena. Throughout most of the hard and hard-intermediate state, the disc is consistent with being truncated (with an inner radius Rin ≳ 10 Rg), reaching close to the innermost stable circular orbit only near transition.

Radiation GRMHD Simulations of the Hard State of Black Hole X-Ray Binaries and the Collapse of a Hot Accretion Flow

Abstract We present global radiation GRMHD simulations of strongly magnetized accretion onto a spinning, stellar mass black hole at sub-Eddington rates. Using a frequency-dependent Monte Carlo procedure for Compton scattering, we self-consistently evolve a two-temperature description of the ion–electron fluid and its radiation field. For an Eddington ratio L/L Edd ≳ 10−3, the emergent spectrum forms an apparent power-law shape from thermal Comptonization up to a cutoff at ≃100 keV, characteristic of that seen in the hard spectral states of black hole X-ray binary systems. At these luminosities, the radiative efficiency is high (≈24%) and results in a denser midplane region where magnetic fields are dynamically important. For L/L Edd ∼ 10−2, our hot accretion flow appears to undergo thermal runaway and collapse. Our simulations demonstrate that hot accretion flows can be radiatively efficient and provide an estimate of their maximum luminosity.

State transitions of GX 339-4 during its outburst rising phase

ABSTRACT We investigate systematically four outbursts of black hole system GX 339-4 observed by the Rossi X-ray Timing Explorer (RXTE) in both spectral and timing domains and find that these outbursts have some common properties, although they experience different ‘q’ tracks in the hardness–intensity diagram (HID). While the spectral indices are around 1.5 in the low/hard state (LHS) and 2.4 in the soft intermediate state (SIMS) and high/soft state (HSS), the spectral parameters of thermal, non-thermal, and reflection components vary significantly in transitions from the LHS to HIMS. Also, the quasi-periodic oscillation (QPO) shows a peculiar behaviour during the state transition between the LHS and HIMS: the RMS drop of the type C fundamental QPO is accompanied by the appearance of the second harmonic. Interestingly, the QPO RMS is found to have a similar linear relationship with the non-thermal fraction of emission in different outbursts. These findings provide more clues to aid our understanding of the outbursts of a black hole X-ray binary system.

Exploring higher order images with Fe Kα-lines from relativistic discs: black hole spin determination and bias

ABSTRACT We study the contributions to the relativistic Fe K α line profile from higher order images (HOIs) produced by strongly deflected rays from the disc which cross the plunging region, located between the innermost stable circular orbit (ISCO) radius and the event horizon of a Kerr black hole. We investigate the characteristics features imprinted by the HOIs in the line profile for different black hole spins, disc emissivity laws, and inclinations. We find that they extend from the red wing of the profile up to energies slightly lower than those of the blue peak, adding ∼0.4–1.3 per cent to the total line flux. The contribution to the specific flux is often in the ∼1 to 7 per cent range, with the highest values attained for low and negative spin (a ≲ 0.3) black holes surrounded by intermediate inclination angle (i ∼ 40°) discs. We simulate future observations of a black hole X-ray binary system with the Large Area Detector of the planned X-ray astronomy enhanced X-ray Timing and Polarimetry Mission (eXTP) and find that the Fe Kα line profiles of systems accreting at ≲1 per cent the Eddington rate are affected by the HOI features for a range of parameters. This would provide evidence of the extreme gravitational lensing of HOI rays. Our simulations show also that not accounting for HOI contributions to the Fe Kα line profile may systematically bias measurements of the black hole spin parameter towards values higher by up to ∼0.3 than the inputted ones.

Optical and ultraviolet monitoring of the black hole X-ray binary MAXI J1820+070/ASASSN-18ey for 18 months

ABSTRACT MAXI J1820+070 is a low-mass black hole X-ray binary system with high luminosity in both optical and X-ray bands during its outburst periods. We present extensive photometry in the X-ray, ultraviolet, and optical bands, as well as densely cadenced optical spectra, covering the phase from the beginning of the optical outburst to ∼550 days. During the rebrightening process, the optical emission precedes the X-ray by 20.80 ± 2.85 days. The spectra are characterized by blue continua and emission features of Balmer series, He i, and He ii lines and a broad Bowen blend. The pseudo equivalent widths (pEWs) of emission lines are found to show anticorrelations with the X-ray flux measured at comparable phases, due to increased suppression by the optical continuum. At around the X-ray peak, the full widths at half-maximum (FWHMs) of Hβ and He ii λ4686 tend to stabilize at 19.4 and 21.8 Å, corresponding to the line-forming regions at radii of 1.7 and 1.3 R⊙ within the disc. We further analysed the absolute fluxes of the lines and found that the fluxes of Hβ and He ii λ4686 show positive correlations with the X-ray flux, favouring the suggestion that the irradiation model is responsible for the optical emission. However, the fact that X-ray emission experiences a dramatic flux drop at t ∼ 200 days after the outburst, while the optical flux only shows a little variation, suggests that additional energy such as viscous energy may contribute to the optical radiation, in addition to X-ray irradiation.

VLBA Discovery of a Resolved Source in the Candidate Black Hole X-Ray Binary AT2019wey

Abstract AT2019wey is a Galactic low-mass X-ray binary with a candidate black hole accretor first discovered as an optical transient by ATLAS in 2019 December. It was then associated with an X-ray source discovered by SRG/eROSITA and SRG/ART-XC instruments in 2020 March. After a brightening in X-rays in 2020 August, VLA observations of the source revealed an optically thin spectrum that subsequently shifted to optically thick, as the source continued to brighten in the radio. This motivated us to observe AT2019wey with the VLBA. We found a resolved source that we interpret to be a steady compact jet, a feature associated with black hole X-ray binary systems in hard X-ray spectral states. The jet power is comparable to the accretion-disk X-ray luminosity. Here, we summarize the results from these observations.

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. II. Magnetohydrodynamic simulations

ABSTRACT Trapped inertial oscillations (r modes) provide a promising explanation for high-frequency quasi-periodic oscillations (HFQPOs) observed in the emission from black hole X-ray binary systems. An eccentricity (or warp) can excite r modes to large amplitudes, but concurrently, the oscillations are likely damped by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). We force eccentricity in global, unstratified, zero-net-flux MHD simulations of relativistic accretion discs and find that a sufficiently strong disc distortion generates trapped inertial waves despite this damping. In our simulations, eccentricities above ∼0.03 in the inner disc excite trapped waves. In addition to the competition between r-mode damping and driving, we observe that larger amplitude eccentric structures modify and in some cases suppress MRI turbulence. Given the variety of distortions (warps as well as eccentricities) capable of amplifying r modes, the robustness of trapped inertial wave excitation in the face of MRI turbulence in our simulations provides support for a discoseismic explanation for HFQPOs.

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. I. Hydrodynamic simulations

ABSTRACT High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of black hole X-ray binary systems promise insight into strongly curved spacetime. ‘Discoseismic’ oscillations with frequencies set by the intrinsic properties of the central black hole, in particular ‘trapped inertial waves’ (r modes), offer an attractive explanation for HFQPOs. To produce an observable signature, however, such oscillations must be excited to sufficiently large amplitudes. Turbulence driven by the magnetorotational instability fails to provide the necessary amplification, but r modes may still be excited via interaction with accretion disc warps or eccentricities. We present 3D global hydrodynamic simulations of relativistic accretion discs, which demonstrate for the first time the excitation of trapped inertial waves by an imposed eccentricity in the flow. While the r modes’ saturated state depends on the vertical boundary conditions used in our unstratified, cylindrical framework, their excitation is unambiguous in all runs with eccentricity ≳ 0.005 near the innermost stable circular orbit. These simulations provide a proof of concept, demonstrating the robustness of trapped inertial wave excitation in a non-magnetized context. In a companion paper, we explore the competition between this excitation, and damping by magnetohydrodynamic turbulence.

Disc and wind in black hole X-ray binary MAXI J1820+070 observed through polarized light during its 2018 outburst

ABSTRACT We describe the first complete polarimetric data set of the entire outburst of a low-mass black hole X-ray binary system and discuss the constraints for geometry and radiative mechanisms it imposes. During the decaying hard state, when the optical flux is dominated by the non-thermal component, the observed polarization is consistent with the interstellar values in all filters. During the soft state, the intrinsic polarization of the source is small, ∼0.15 per cent in B and V filters, and is likely produced in the irradiated disc. A much higher polarization, reaching ∼0.5 per cent in V and R filters, at a position angle of ∼25○ observed in the rising hard state coincides in time with the detection of winds in the system. This angle coincides with the position angle of the jet. The detected optical polarization is best explained by scattering of the non-thermal (hot flow or jet base) radiation in an equatorial wind.