🔭 astrophysics

State-Dependent X-ray Variability in Cygnus X-1: A 12-Year NuSTAR Timing Study of Accretion Flow Geometry

This 12-year NuSTAR timing study of Cygnus X-1 reveals energy-dependent bimodal flux distributions and systematic state-dependent changes in variability frequencies, time lags, and disk truncation radii that map the geometric evolution of the accretion flow, while identifying a unique failed state transition that challenges standard models of wind-fed black hole accretion.

Original authors: Kshitij Duraphe, Kartik Mandar, Chooda Khanal, Abha Pareek, Tejaswi Kondhiya, V Sree Suswara, Deeksha Dinesh, Vidyasagar Bhat, Gopal Bhatta

Published 2026-03-24

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kshitij Duraphe, Kartik Mandar, Chooda Khanal, Abha Pareek, Tejaswi Kondhiya, V Sree Suswara, Deeksha Dinesh, Vidyasagar Bhat, Gopal Bhatta

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a cosmic dance floor where a massive, invisible partner (a black hole) is spinning wildly while trying to eat a meal from a giant, glowing star next to it. This cosmic couple is called Cygnus X-1, and for the last 12 years, astronomers have been watching it with a very special, high-powered camera called NuSTAR.

This paper is like a 12-year diary of that dance, revealing how the black hole eats, how it moves, and how it sometimes gets confused. Here is the story in simple terms:

1. The Two Ways of Eating (The States)

Think of the black hole's "meal" (accretion flow) as having two distinct moods:

The Hard State (The Grumpy Eater): The black hole is surrounded by a hot, chaotic cloud of gas (a corona). It eats slowly but steadily, and the food gets super-heated, glowing with high-energy X-rays. It's like a campfire: hot, bright, and crackling with noise.
The Soft State (The Smooth Eater): The hot cloud collapses, and the food (gas) forms a smooth, cool, flat disk right up against the black hole. It glows with a softer, steadier light. It's like a calm, flowing river.

Usually, the black hole switches between these two moods. Sometimes it's grumpy, sometimes it's smooth.

2. The "Bimodal" Mystery (The Two Peaks)

The researchers noticed something fascinating about the light coming from the black hole.

If you look at the low-energy light (like the warm glow of the river), the brightness changes smoothly. It's like a dimmer switch; you can see every step between bright and dark.
But if you look at the high-energy light (the crackling fire), the brightness is bimodal. This means it only really exists in two modes: either "Very Bright" or "Very Dim." It rarely stays in the middle. It's like a light switch that is either fully ON or fully OFF, with almost no "halfway" setting.

The Analogy: Imagine a thermostat. The temperature of the room (the soft light) changes gradually. But the heater itself (the hard light) only has two settings: "Full Blast" or "Off." The paper proves that this "On/Off" switch happens specifically in the high-energy part of the spectrum.

3. The Rhythm of the Dance (Timing and Frequency)

Just like a drummer has a beat, the black hole has a rhythm.

In the Hard State, the "beat" is slow and heavy (about 0.05 beats per second). This means the gas is swirling further away from the black hole.
In the Intermediate State (the transition), the beat speeds up (0.074 beats per second). This is the "aha!" moment: as the beat speeds up, it proves the gas disk is physically moving closer to the black hole, like a figure skater pulling their arms in to spin faster.
In the Soft State, the beat disappears into a smooth hum (red noise), because the gas is now right up against the black hole's edge.

4. The "Echo" Effect (Time Lags)

The researchers also listened to the "echoes" of the light.

In the Hard State, the high-energy light arrives slightly after the low-energy light (about 50 milliseconds later). Think of this like a shout in a canyon: the sound (soft light) leaves first, and the echo (hard light, bounced off the hot cloud) comes back a split second later.
In the Soft State, the echo vanishes. The shout and the echo happen at the same time. This tells us the "canyon" (the hot cloud) has collapsed, and the sound is traveling directly from the source.

5. The Glitch: The "Failed" Transition

The most exciting discovery is about one specific observation (Observation 30302019006).

What happened? The black hole started to switch from "Grumpy" to "Smooth." The light got softer, suggesting the gas was moving inward.
The Glitch: But then, the dance froze. The light stopped fluctuating. It became incredibly calm, almost too calm.
The Analogy: Imagine a car trying to shift gears. It revs up, the engine noise changes, but then the car just... stalls. It doesn't go into the new gear, and it doesn't go back to the old one either.
Why? The researchers think the black hole's companion star blew a strong "wind" that messed up the gas flow. The disk tried to move closer to the black hole but got pushed back or disrupted before it could settle. It was a failed state transition.

Why Does This Matter?

This study is like finding the "Rosetta Stone" for black holes.

It proves the geometry: We now have direct proof that when a black hole changes its mood, the physical shape of the gas around it actually changes (moving closer or further away).
It explains the "On/Off" switch: We learned that the "bimodal" behavior is a property of the hot cloud, not the whole system.
It found a new path: The "failed transition" suggests that black holes don't always follow a perfect script. Sometimes, the environment (like stellar winds) can interrupt the process, creating weird, unstable states we haven't seen before.

In a nutshell: This paper used 12 years of data to show us exactly how a black hole eats, how its "room" changes size, and how sometimes, even a black hole can get stuck in a traffic jam.

1. Problem Statement

Cygnus X-1, the prototypical high-mass X-ray binary (HMXB) containing a stellar-mass black hole, exhibits complex state transitions between hard and soft spectral states. While the coupling between accretion flow geometry, spectral states, and jet ejection is well-established, the precise physical mechanisms governing the state transitions and the geometry of the accretion flow (specifically the truncation radius of the disk and the extent of the Comptonizing corona) remain subjects of active debate.

Previous studies have often relied on transient black hole binaries or lower-energy instruments. There is a need for a comprehensive, long-term timing analysis of a persistent source like Cygnus X-1 using high-sensitivity hard X-ray data (3–79 keV) to:

Determine the energy dependence of flux bimodality.
Quantify the evolution of characteristic frequencies and time lags across states.
Investigate anomalies in state transitions, such as "failed" transitions where spectral softening occurs without expected timing variability.

2. Methodology

The authors conducted a comprehensive timing analysis of 26 NuSTAR observations of Cygnus X-1 spanning 2012–2024 (a 12-year baseline).

Data Reduction: Raw Level 1 data were processed using NuSTARDAS (HEASoft) and the June 2025 CALDB. Standard screening removed bad time intervals and high particle background. Light curves were extracted for both focal plane modules (FPMA/FPMB) with a 0.1s time resolution to capture rapid aperiodic variability.
Energy Band Definition: A data-driven approach was used to define optimal energy bands. The authors analyzed count rate distributions across nine energy bands to identify where the well-known bimodal state distribution emerges. They selected 3–8 keV (Soft) and 8–79 keV (Hard) to separate thermal disk emission from Comptonized emission.
Statistical Analysis:
- Bimodality Testing: Count rate distributions were fitted with Normal, Lognormal, and Gaussian Mixture Models (GMM) using the Bayesian Information Criterion (BIC).
- Variability Metrics: Fractional variability ( $F_{var}$ ) and modulation factors were calculated.
- Power Spectral Density (PSD): PSDs were computed using Welch's method (256s segments, 50% overlap) in the 0.01–5 Hz range. The spectra were modeled with Lorentzian functions to identify break frequencies.
- Time Lag Analysis: Frequency-dependent time lags were derived via cross-spectral analysis between soft and hard bands to probe causal connections and propagation delays.
- RMS-Flux Relations: Linear correlations between RMS variability and mean flux were computed for 10s segments to test the propagating fluctuation model.

3. Key Contributions

Energy-Dependent Bimodality: The study definitively establishes that raw count rate bimodality is intrinsic to the Comptonized component above ~10 keV, while thermal emission below 8 keV remains unimodal. This refines the understanding of how spectral states manifest in different energy regimes.
Geometric Mapping: By converting characteristic break frequencies to physical scales, the authors provide direct observational evidence of the inward progression of the accretion disk during state transitions, mapping the truncation radius evolution from ~5.5 $R_g$ (hard) to ~2 $R_g$ (intermediate).
Discovery of a "Failed" Transition: The identification of a specific observation (ObsID 30302019006) exhibiting anomalous behavior challenges standard transition models, suggesting a new pathway for accretion flow evolution in wind-fed systems.

4. Key Results

A. Spectral State and Bimodality

Hard vs. Soft Bands: Below 8 keV, the count rate distribution is unimodal, reflecting a continuous range of thermal disk temperatures. Above 10 keV, the distribution becomes distinctly bimodal, reflecting the bistable nature of the Comptonizing corona (extended in hard states, suppressed in soft states).
Hardness-Intensity Diagram (HID): The source spends ~~50% of its time in the hard state. Unlike transient sources, Cygnus X-1 does not show a clear q-shaped hysteresis loop, consistent with its persistent nature and stable mass accretion rate (~~0.01–0.1 $L_{Edd}$ ).

B. Timing and Power Spectra

Characteristic Frequencies: The break frequency in the PSD shifts systematically:
- Hard State: Median $\sim$ 0.050 Hz.
- Intermediate State: Median $\sim$ 0.074 Hz.
- Soft State: Featureless red noise (no discernible break).
Truncation Radius: Using the empirical relation $R_g \sim 6 (50/m \cdot f_{LB})^{2/3}$ , the break frequency shift corresponds to the disk truncation radius moving inward from ~5.5 $R_g$ (hard) to ~2 $R_g$ (intermediate).
Time Lags:
- Hard State: Shows systematic positive (hard) lags of ~50 ms at low frequencies, decreasing with frequency. This indicates propagating fluctuations from the outer corona to the inner region.
- Soft State: Lags are suppressed to near-zero, indicating the collapse of the Comptonizing corona and the disk extending close to the Innermost Stable Circular Orbit (ISCO).
- Intermediate State: Shows complex behavior with a peak lag of ~75 ms near 0.2 Hz.

C. RMS-Flux Relations

A linear relationship between RMS variability and mean flux persists across all states, confirming the multiplicative coupling of accretion fluctuations.
Parameter Evolution: The slope of the relation decreases from 0.609 (soft) to 0.581 (hard), while the y-intercept increases. This tracks the changing relative contributions of thermal (disk) vs. non-thermal (corona/jet) emission.

D. The "Failed" State Transition (ObsID 30302019006)

Anomalies: This observation showed an intermediate hardness ratio but:
- Suppressed Variability: $F_{var} < 1.38\%$ (far below the typical >10% for intermediate states).
- Anticorrelated Flux: Soft and hard band fluxes were negatively correlated, unlike the positive correlation seen in other states.
- Sub-ISCO Truncation: The break frequency implied a truncation radius $< 1 R_g$ , seemingly inside the ISCO.
Interpretation: The authors propose this represents a failed transition where the disk began to cool and move inward but failed to establish a stable soft-state configuration, possibly due to disruption by the stellar wind of the O9.7Iab companion.

5. Significance

Validation of Propagating Fluctuation Model: The persistence of linear rms-flux relations and the systematic evolution of time lags strongly support the model where accretion rate fluctuations propagate inward through the flow.
Refinement of State Transition Physics: The discovery of the "failed transition" suggests that state transitions in wind-fed HMXBs (like Cygnus X-1) may follow different evolutionary paths compared to Roche-lobe overflow systems, potentially explaining anomalous behaviors in other sources and Changing-Look AGNs.
Benchmark for General Relativity: The precise constraints on the disk truncation radius and corona geometry provide critical inputs for testing GR models and understanding the interaction between the accretion flow and the black hole spin.
Universal Scaling: The tight correlations between spectral hardness, timescales, and variability amplitudes suggest universal scaling relations that may apply from stellar-mass to supermassive black holes.

In conclusion, this 12-year NuSTAR study provides the most detailed characterization to date of Cygnus X-1's accretion flow, linking spectral state changes directly to geometric reconfigurations of the disk and corona, while uncovering new complexities in how these transitions can fail.