Characterisation of the stellar wind in Cyg X-1 via modelling of colour-colour diagrams
E. V. Lai,B. De Marco,Y. Cavecchi,I. El Mellah,M. Cinus,C. M. Diez,V. Grinberg,A. A. Zdziarski,P. Uttley,M. Bachetti,J. José,G. Sala,A. Różańska,J. Wilms
2024-08-12
Abstract:Cygnus X-1 is a high mass X-ray binary where accretion onto the black hole is mediated by the stellar wind from the blue supergiant companion star HDE 226868. Depending on the position of the black hole along the orbit, X-ray observations can probe different layers of the stellar wind. Deeper wind layers can be investigated at superior conjunction (i.e. null orbital phases). We aim at characterising the stellar wind in the Cyg X-1/HDE 226868 system analysing one passage at superior conjunction covered by XMM-Newton during the CHOCBOX campaign via modelling of colour-colour diagrams. Since X-ray absorption is energy-dependent, colour indices provide information on the parameters of the stellar wind, such as the column density $N_{H,w}$ and the covering factor $f_c$. We fitted colour-colour diagrams with models that include both a continuum and a stellar wind component. We used the KDE method to infer the unknown probability distribution of the data points in the colour-colour diagram, and selected the model corresponding to the highest likelihood. In order to study the temporal evolution of the wind around superior conjunction, we extracted and fitted time-resolved colour-colour diagrams. We found that the model that best describes the shape of the colour-colour diagram of Cyg X-1 at superior conjunction requires the wind to be partially ionised. The shape of the colour-colour diagram strongly varies during the analysed observation, as due to concurrent changes of the mean $N_{H,w}$ and the $f_c$ of the wind. Our results suggest the existence of a linear scaling between the rapid variability amplitude of $N_{H,w}$ (on time scales between 10 s and 11 ks) and its long term variations (on time scales 11>ks). Using the inferred best-fit values, we estimated the stellar mass loss rate to be $\sim 7\times10^{-6} {\rm M_{\odot}yr^{-1}}$ and the clumps to have a mass of $\sim10^{17}$ g.
High Energy Astrophysical Phenomena