Properties of the propagating shock wave in the accretion flow around GX 339-4 in the 2010 outburst

Context. The black hole candidate GX 339-4 exhibited an X-ray outburst in January 2010, which is still continuing. We here discuss the timing and the spectral properties of the outburst using RXTE data. Aims. Our goal is to study the timing and spectral properties of GX 339-4 using its recent outburst data and extract information about the nature of the accretion flow. Methods. We use RXTE archival data of the recent GX 339-4 outburst and analyze them with the NASA HEAsoft package, version 6.8. We then compare the observed quasi-periodic oscillation (QPO) frequencies with those from existing shock oscillation model and obtain the nature of evolution of the shock locations during the outburst. Results. We found that the QPO frequencies are monotonically increasing from 0.102 Hz to 5.69 Hz within a period of ~ 26 days. We explain this evolution with the propagating oscillatory shock (POS) solution and find the variation of the initial and final shock locations and strengths. The model fits also give the velocity of the propagating shock wave, which is responsible for the generation of QPOs and their evolutions, at ~ 10 m/s. We observe from the spectra that up to 2010 April 10, the object was in a hard state. After that, it went to the hard-intermediate state. On April 18, it had a state transition and went to the soft-intermediate state. On May 15, another state transition was observed and the source moved to the soft state. Conclusions. As in the previously fitted outburst sources, this source also showed the tendency of a rapidly increasing QPO frequency ($\nu_{QPO}$) in a viscous time scale, which can be modeled quite accurately. In this case, the shock seems to have disappeared at about ~ 172 Schwarzschild radii, unlike in the 2005 outburst of GRO J1655-40, where the shock disappeared behind the horizon.

💡 Research Summary

The paper presents a comprehensive analysis of the 2010 outburst of the black‑hole candidate GX 339‑4 using archival Rossi X‑ray Timing Explorer (RXTE) data. The authors focus on two complementary aspects: timing properties, particularly the evolution of low‑frequency quasi‑periodic oscillations (QPOs), and spectral properties, including state transitions as traced by hardness‑intensity diagrams (HIDs).

Data and Methodology
The study employs both the All‑Sky Monitor (ASM) and the Proportional Counter Array (PCA) instruments. For timing analysis, the authors extract 0.01 s binned light curves from PCA Event mode data, generate power density spectra (PDS) with the XRONOS “powspec” tool, and fit the QPO peaks with Lorentzian profiles to obtain centroid frequencies, quality factors, and rms amplitudes. Spectral analysis uses PCA Standard‑2 data in the 3–25 keV band, modeled with a multicolour disc blackbody (diskbb) plus a power‑law component, supplemented by a Gaussian Fe Kα line at ~6.5 keV. Absorption is fixed at N_H = 5 × 10²¹ cm⁻², and a 1 % systematic error is added.

QPO Evolution and the Propagating Oscillatory Shock (POS) Model
QPOs first appear on 22 March 2010 (MJD 55277) at 0.102 Hz and increase monotonically to 5.69 Hz by 17 April 2010 (MJD 55303), a span of ~26 days. The authors interpret this monotonic rise with the POS solution, originally proposed by Chakrabarti and collaborators. In this framework, the QPO frequency is the inverse of the infall time from the post‑shock region; as the shock front propagates inward, the infall time shortens and the QPO frequency rises. By fitting the observed frequencies, they infer an initial shock location of ≈1500 r_g (Schwarzschild radii) moving inward at a constant speed of ≈10 m s⁻¹. The shock weakens during the propagation: the compression ratio R = ρ₋/ρ₊ declines from an initial value R₀≈4.2 to ≈1 (the weakest possible shock) following R⁻¹ = R₀⁻¹ + α t_d², with α≈6.8 × 10⁻⁴ day⁻². The shock disappears at ≈172 r_g, significantly farther out than the event horizon, indicating that the QPOs are generated well outside the black‑hole’s immediate vicinity.

Spectral Evolution and State Transitions
The HID constructed from PCA count rates (2–20 keV) and hardness ratios (6–20 keV / 2–6 keV) shows the classic “q‑track” with four distinct branches:

1. Hard State (A, up to MJD 55296) – Dominated by a hard power‑law (photon index Γ < 2) and strong 15–30 keV emission.
2. Hard‑Intermediate State (B, MJD 55296–55304) – Rapid softening; hard flux drops sharply while the soft (2–4 keV) disc component rises.
3. Soft‑Intermediate State (C, MJD 55304–55331) – Both disc and power‑law components become comparable; QPOs persist but at a roughly constant high frequency.
4. Soft State (D, after MJD 55331) – Disc dominates, overall count rate declines, and the photon index steepens to Γ ≈ 2.5.

These transitions are corroborated by the timing behavior: QPOs are present and type‑C (hard and hard‑intermediate states) until the shock weakens; sporadic QPOs of type‑B appear in the soft‑intermediate state, and disappear in the soft state.

Physical Interpretation
The authors adopt a two‑component accretion flow model: a sub‑Keplerian, low‑angular‑momentum halo that forms a standing shock, and a Keplerian, high‑angular‑momentum disc that moves inward as viscosity increases. During the rising phase, the halo dominates, producing a strong shock and the observed low‑frequency QPOs. As the Keplerian disc penetrates deeper, the shock weakens (compression ratio declines) and eventually vanishes, coinciding with the hard‑to‑soft state transition. The inferred shock speed (∼10 m s⁻¹) is slower than that derived for GRO J1655‑40 and XTE J1550‑564 (∼20 m s⁻¹), suggesting source‑specific viscosity or angular‑momentum transport properties.

Conclusions
The study demonstrates that the POS model successfully reproduces the rapid QPO frequency increase observed in GX 339‑4, extending the applicability of this framework beyond previously studied outbursts (e.g., GRO J1655‑40, XTE J1550‑564). The shock’s disappearance at ~172 r_g, rather than plunging into the horizon, provides a clear observational signature linking shock dynamics to spectral state evolution. The work reinforces the view that QPOs are a manifestation of oscillating shocks in the inner accretion flow and that state transitions are driven by the relative dominance of Keplerian versus sub‑Keplerian components, modulated by changes in disk viscosity. Future high‑resolution timing missions and magnetohydrodynamic simulations will be essential to refine the shock‑oscillation paradigm and to explore its connection with jet ejection events often observed in GX 339‑4.