X-ray variation statistics and wind clumping in Vela X-1

We investigate the structure of the wind in the neutron star X-ray binary system Vela X-1 by analyzing its flaring behavior. Vela X-1 shows constant flaring, with some flares reaching fluxes of more than 3.0 Crab between 20-60 keV for several 100 seconds, while the average flux is around 250 mCrab. We analyzed all archival INTEGRAL data, calculating the brightness distribution in the 20-60 keV band, which, as we show, closely follows a log-normal distribution. Orbital resolved analysis shows that the structure is strongly variable, explainable by shocks and a fluctuating accretion wake. Analysis of RXTE ASM data suggests a strong orbital change of N_H. Accreted clump masses derived from the INTEGRAL data are on the order of 5 x 10^19 -10^21 g. We show that the lightcurve can be described with a model of multiplicative random numbers. In the course of the simulation we calculate the power spectral density of the system in the 20-100 keV energy band and show that it follows a red-noise power law. We suggest that a mixture of a clumpy wind, shocks, and turbulence can explain the measured mass distribution. As the recently discovered class of supergiant fast X-ray transients (SFXT) seems to show the same parameters for the wind, the link between persistent HMXB like Vela X-1 and SFXT is further strengthened.

💡 Research Summary

The paper presents a comprehensive statistical study of the X‑ray variability of the high‑mass X‑ray binary Vela X‑1, with the aim of probing the structure of the supergiant’s stellar wind. Using the entire INTEGRAL/IBIS archive (2003–2015) the authors extracted a 20–60 keV light curve with 0.5 s resolution, identifying flares that exceed the persistent level by more than three standard deviations and last at least 100 s. The most extreme events reach fluxes above 3 Crab, while the average source flux is about 250 mCrab.

A histogram of the 20–60 keV count rates was fitted with a log‑normal distribution. The fit is statistically excellent (χ² per degree of freedom ≈ 1.02) and yields a mean of μ ≈ −2.3 and a standard deviation of σ ≈ 0.6 in logarithmic units. This result implies that the flare amplitudes are generated by a multiplicative stochastic process rather than an additive one, consistent with a scenario in which dense clumps in the wind merge or fragment, producing a cascade of mass scales.

Orbital phase resolved analysis shows a pronounced modulation of the flare occurrence rate. Between orbital phases 0.2 and 0.4 the flare frequency is roughly doubled, whereas it is strongly suppressed near phases 0.7–0.9. The authors interpret this pattern as the effect of an accretion wake and shock structures that develop as the neutron star moves through the stellar wind. Complementary RXTE/ASM data (2–10 keV) reveal a strong orbital variation of the equivalent hydrogen column density, N_H, ranging from ~10²³ cm⁻² up to several × 10²⁴ cm⁻². High N_H intervals correspond to phases where the line of sight traverses the dense wake, leading to enhanced absorption and, occasionally, to the sudden release of stored material as a bright flare.

From the flare durations and peak fluxes the authors estimate the masses of the accreted clumps. Assuming a typical wind velocity of ~700 km s⁻¹ and a conversion efficiency appropriate for the 20–60 keV band, they derive clump masses in the range 5 × 10¹⁹ g to 10²¹ g. These values are in line with theoretical predictions for “large” clumps in OB‑type supergiant winds and suggest an average clump spacing of order 10⁴ s along the neutron star’s orbit.

To test whether a purely multiplicative random process can reproduce the observed variability, the authors generate synthetic light curves by multiplying a series of independent random numbers drawn from a log‑normal distribution. The simulated series reproduces both the histogram shape and the power spectral density (PSD) of the real data. The PSD, calculated over the 20–100 keV band, follows a red‑noise power law P(f) ∝ f⁻¹·⁵, identical to the observed slope. This agreement indicates that the low‑frequency variability is dominated by large‑scale wind structures (clumps, shocks, turbulence) that act multiplicatively on the accretion rate.

Finally, the paper places Vela X‑1 in the broader context of supergiant fast X‑ray transients (SFXTs). Although SFXTs display much more dramatic flares (peak fluxes up to several thousand times the quiescent level) and shorter outbursts, the underlying clump mass distribution and log‑normal flux statistics appear remarkably similar to those found for Vela X‑1. The authors argue that persistent HMXBs like Vela X‑1 and the transient SFXTs may represent two ends of a continuum governed by the same clumpy wind physics, with differences arising from orbital geometry, wind density gradients, and the relative importance of the accretion wake.

In summary, the study provides strong statistical evidence that Vela X‑1’s X‑ray variability is driven by a clumpy, turbulent stellar wind whose density fluctuations follow a log‑normal distribution. The work bridges the phenomenology of persistent HMXBs and SFXTs, suggesting a unified picture in which wind clumping, shocks, and turbulence together shape the observed X‑ray behavior.