Refine
Year of publication
Document Type
- Article (105)
- Other (6)
- Conference Proceeding (3)
- Doctoral Thesis (1)
- Postprint (1)
- Review (1)
Language
- English (117)
Keywords
- stars: massive (31)
- stars: early-type (30)
- X-rays: stars (26)
- X-rays: binaries (21)
- stars: winds, outflows (20)
- stars: magnetic field (17)
- stars: mass-loss (17)
- stars: atmospheres (16)
- stars: Wolf-Rayet (14)
- Magellanic Clouds (11)
Institute
- Institut für Physik und Astronomie (117) (remove)
Context. Supergiant fast X-ray (SFXT) transients are a peculiar class of supergiant X-ray binaries characterized by a remarkable variability in the X-ray domain, widely ascribed to accretion from a clumpy stellar wind. Aims. In this paper we performed a systematic and homogeneous analysis of the sufficiently bright X-ray flares observed with XMM-Newton from the supergiant fast X-ray transients to probe spectral variations on timescales as short as a few hundred seconds. Our ultimate goal is to investigate whether SFXT flares and outbursts are triggered by the presence of clumps, and to reveal whether strongly or mildly dense clumps are required. Methods. For all sources, we employ a technique developed by our group already exploited in a number of our previous papers, making use of an adaptive rebinned hardness ratio to optimally select the time intervals for the spectral extraction. A total of twelve observations performed in the direction of five SFXTs are reported, providing the largest sample of events available so far. Results. Using the original results reported here and those obtained with our technique from the analysis of two previously published XMM-Newton observations of IGR J17544-2619 and IGR J18410-0535, we show that both strongly and mildly dense clumps can trigger these events. In the former case, the local absorption column density may increase by a factor of >> 3, while in the latter case, the increase is only a factor of similar to 2-3 (or lower). An increase in the absorption column density is generally recorded during the rise of the flares/outbursts, while a drop follows when the source achieves peak flux. In a few cases, a re-increase of the absorption column density after the flare is also detected, and we discovered one absorption event related to the passage of an unaccreted clump in front of the compact object. Overall, there seems to be no obvious correlation between the dynamic ranges in the X-ray fluxes and absorption column densities in supergiant fast X-ray transients, with an indication that lower densities are recorded at the highest fluxes. Conclusions. The spectral variations measured in all sources are in agreement with the idea that the flares/outbursts are triggered by the presence of dense structures in the wind interacting with the X-rays from the compact object (leading to photoionization). The lack of correlation between the dynamic ranges in the X-ray fluxes and absorption column densities can be explained by the presence of accretion inhibition mechanism(s). Based on the knowledge acquired so far on the SFXTs, we propose a classification of the flares/outbursts from these sources in order to drive future observational investigations. We suggest that the difference between the classes of flares/outbursts is related to the fact that the mechanism(s) inhibiting accretion can be overcome more easily in some sources compared to others. We also investigate the possibility that different stellar wind structures, other than clumps, could provide the means to temporarily overcome the inhibition of accretion in supergiant fast X-ray transients.
In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting approximately 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high, we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detections of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousand seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.
3A 1954+319 has been classified for a long time as a symbiotic X-ray binary, hosting a slowly rotating neutron star and an aged M red giant. Recently, this classification has been revised thanks to the discovery that the donor star is an M supergiant. This makes 3A 1954+319 a rare type of high-mass X-ray binary consisting of a neutron star and a red supergiant donor. In this paper, we analyse two archival and still unpublished XMM-Newton and NuSTAR observations of the source. We perform a detailed hardness ratio-resolved spectral analysis to search for spectral variability that could help investigating the structures of the inhomogeneous M supergiant wind from which the neutron star is accreting. We discuss our results in the context of wind-fed supergiant X-ray binaries and show that the newest findings on 3A 1954+319 reinforce the hypothesis that the neutron star in this system is endowed with a magnetar-like magnetic field strength (greater than or similar to 10(14) G).
The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the nonstationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total timescale of several hours), the transitions of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the nonstationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favored in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.
The long-term X-ray light curves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of three-dimensional (3D) hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ 16493-4348, which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light curve are similar to the observed ones, thus providing support to our scenario. We propose that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.
Symbiotic X-ray binaries are systems hosting a neutron star accreting form the wind of a late-type companion. These are rare objects and so far only a handful of them are known. One of the most puzzling aspects of the symbiotic X-ray binaries is the possibility that they contain strongly magnetized neutron stars. These are expected to be evolutionary much younger compared to their evolved companions and could thus be formed through the (yet poorly known) accretion induced collapse of a white dwarf. In this paper, we perform a broad-band X-ray and soft gamma-ray spectroscopy of two known symbiotic binaries, Sct X-1 and 4U 1700+24, looking for the presence of cyclotron scattering features that could confirm the presence of strongly magnetized NSs. We exploited available Chandra, Swift, and NuSTAR data. We find no evidence of cyclotron resonant scattering features (CRSFs) in the case of Sct X-1 but in the case of 4U 1700+24 we suggest the presence of a possible CRSF at similar to 16 keV and its first harmonic at similar to 31 keV, although we could not exclude alternative spectral models for the broad-band fit. If confirmed by future observations, 4U 1700+24 could be the second symbiotic X-ray binary with a highly magnetized accretor. We also report about our long-term monitoring of the last discovered symbiotic X-ray binary IGR J17329-2731 performed with Swift/XRT. The monitoring revealed that, as predicted, in 2017 this object became a persistent and variable source, showing X-ray flares lasting for a few days and intriguing obscuration events that are interpreted in the context of clumpy wind accretion.
The information content of data on rotationally periodic recurrent discrete absorption components (DACs) in hot star wind emission lines is discussed. The data comprise optical depths tau(w,phi) as a function of dimensionless Doppler velocity w=(Deltalambda/lambda(0))(c/v(infinity)) and of time expressed in terms of stellar rotation angle phi. This is used to study the spatial distributions of density, radial and rotational velocities, and ionisation structures of the corotating wind streams to which recurrent DACs are conventionally attributed. The simplifying assumptions made to reduce the degrees of freedom in such structure distribution functions to match those in the DAC data are discussed and the problem then posed in terms of a bivariate relationship between tau(w, phi) and the radial velocity v(r)(r), transverse rotation rate Omega(r) and density rho(r, phi) structures of the streams. The discussion applies to cases where: the streams are equatorial; the system is seen edge on; the ionisation structure is approximated as uniform; the radial and transverse velocities are taken to be functions only of radial distance but the stream density is allowed to vary with azimuth. The last kinematic assumption essentially ignores the dynamical feedback of density on velocity and the relationship of this to fully dynamical models is discussed. The case of narrow streams is first considered, noting the result of Hamann et al. (2001) that the apparent acceleration of a narrow stream DAC is higher than the acceleration of the matter itself, so that the apparent slow acceleration of DACs cannot be attributed to the slowness of stellar rotation. Thus DACs either involve matter which accelerates slower than the general wind flow, or they are formed by structures which are not advected with the matter flow but propagate upstream (such as Abbott waves). It is then shown how, in the kinematic model approximation, the radial speed of the absorbing matter can be found by inversion of the apparent acceleration of the narrow DAC, for a given rotation law. The case of broad streams is more complex but also more informative. The observed tau(w,phi) is governed not only by v(r)(r) and Omega(r) of the absorbing stream matter but also by the density profile across the stream, determined by the azimuthal (phi(0)) distribution function F- 0(phi(0)) of mass loss rate around the stellar equator. When F-0(phi(0)) is fairly wide in phi(0), the acceleration of the DAC peak tau(w, phi) in w is generally slow compared with that of a narrow stream DAC and the information on v(r)(r), Omega(r) and F-0(phi(0)) is convoluted in the data tau(w, phi). We show that it is possible, in this kinematic model, to recover by inversion, complete information on all three distribution functions v(r)(r), Omega(r) and F- 0(phi(0)) from data on tau(w, phi) of sufficiently high precision and resolution since v(r)(r) and Omega(r) occur in combination rather than independently in the equations. This is demonstrated for simulated data, including noise effects, and is discussed in relation to real data and to fully hydrodynamic models
The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS 2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33 000 +/- 1000K and 4.0 +/- 0.1 dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H alpha shows a variable emission. The H alpha emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.
Sk 183 is the visually brightest star in the N90 nebula, a young star-forming region in the Wing of the Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope which reveals Sk 183 to be one of the most massive O-type stars in the SMC. Classified as an O3-type dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption, which suggests a later type. We propose that Sk 183 has a composite spectrum and that it is similar to another star in the SMC, MPG 324. This brings the number of rare O2- and O3-type stars known in the whole of the SMC to a mere four. We estimate physical parameters for Sk 183 from analysis of its spectrum. For a single-star model, we estimate an effective temperature of 46 +/- 2 kK, a low mass-loss rate of similar to 10(-7) M-circle dot yr(-1), and a spectroscopic mass of 46(-8)(+ 9) M-circle dot (for an adopted distance modulus of 18.7 mag to the young population in the SMC Wing). An illustrative binary model requires a slightly hotter temperature (similar to 47.5 kK) for the primary component. In either scenario, Sk 183 is the earliest-type star known in N90 and will therefore be the dominant source of hydrogen-ionizing photons. This suggests Sk 183 is the primary influence on the star formation along the inner edge of the nebula.
We have developed lists of likely B3-A0 stars (called "late B" stars) in the young cluster Trumpler 16. The following criteria were used: location within 3' of eta Car, an appropriate V and B - V combination, and proper motion (where available). Color and magnitude cuts have been made assuming an E(B - V) = 0.55 mag +/- 0.1, which is a good approximation close to the center of Trumpler 16. These lists have been cross-correlated with X-ray sources found in the Chandra Carina Complex Project. Previous studies have shown that only very rarely (if at all) do late main-sequence B stars produce X-rays. We present evidence that the X-ray-detected sources are binaries with low-mass companions, since stars less massive than 1.4 M-circle dot are strong X-ray sources at the age of the cluster. Both the median X-ray energies and X-ray luminosities of these sources are in good agreement with values for typical low-mass coronal X-ray sources. We find that 39% of the late B stars based on a list with proper motions have low-mass companions. Similarly, 32% of a sample without proper motions have low-mass companions. We discuss the X-ray detection completeness. These results on low-mass companions of intermediate-mass stars are complementary to spectroscopic and interferometric results and probe new parameter space of low-mass companions at all separations. They do not support a steeply rising distribution of mass ratios to low masses for intermediate-mass (5 M-circle dot) primaries, such as would be found by random pairing from the initial mass function.
We analyze the expansion of hydrogen-poor knots and filaments in the born-again planetary nebulae A30 and A78 based on Hubble Space Telescope ( HST) images obtained almost 20 yr apart. The proper motion of these features generally increases with distance to the central star, but the fractional expansion decreases, i.e., the expansion is not homologous. As a result, there is not a unique expansion age, which is estimated to be 610-950 yr for A30 and 600-1140 yr for A78. The knots and filaments have experienced complex dynamical processes: the current fast stellar wind is mass loaded by the material ablated from the inner knots; the ablated material is then swept up until it shocks the inner edges of the outer, hydrogen-rich nebula. The angular expansion of the outer filaments shows a clear dependence on position angle, indicating that the interaction of the stellar wind with the innermost knots channels the wind along preferred directions. The apparent angular expansion of the innermost knots seems to be dominated by the rocket effect of evaporating gas and by the propagation of the ionization front inside them. Radiation-hydrodynamical simulations show that a single ejection of material followed by a rapid onset of the stellar wind and ionizing flux can reproduce the variety of clumps and filaments at different distances from the central star found in A30 and A78.
Clumps in hot star winds can originate from shock compression due to the line driven instability. One-dimensional hydrodynamic simulations reveal a radial wind structure consisting of highly compressed shells separated by voids, and colliding with fast clouds. Two-dimensional simulations are still largely missing, despite first attempts. Clumpiness dramatically affects the radiative transfer and thus all wind diagnostics in the UV, optical, and in X-rays. The microturbulence approximation applied hitherto is currently superseded by a more sophisticated radiative transfer in stochastic media. Besides clumps, i.e. jumps in the density stratification, so-called kinks in the velocity law, i.e. jumps in dv/dr, play an eminent role in hot star winds. Kinks are a new type of radiative-acoustic shock, and propagate at super-Abbottic speed.
We discuss X-ray line formation in dense O star winds. A random distribution of wind shocks is assumed to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep decline. These box- type, blueshifted profiles resemble recent Chandra observations of the O3 star zeta Pup. For partially transparent fragments, the emission lines become similar to those from a homogeneous wind.
We discuss the connection between wind overloading and discrete absorption components in P Cygni line profiles from O stars. Overloading can create horizontal plateaus in the radial wind speed that cause the extra absorption in the line profile. The upstream propagation speed of these velocity plateaus is analyzed. The second part of the paper deals with X-ray emission from O stars. X-ray line profiles observed with Chandra and XMM are often symmetric, contrary to what is expected for lines from a homogeneous wind. We discuss the influence on line symmetry of photon escape channels in a strongly clumped wind.
Only a small fraction of massive stars seem to host a measurable structured magnetic field, whose origin is still unknown and whose implications for stellar evolution still need to be assessed. Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD 44743; B1 II/III) and epsilon CMa (HD 52089; B1.5II) in December 2013 and April 2014. For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field, approximately constant with time. We determined the physical parameters of both stars and characterise their X-ray spectrum. For the beta Cep star beta CMa, our mode identification analysis led to determining a rotation period of 13.6 +/- 1.2 days and of an inclination angle of the rotation axis of 57.6 +/- 1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of about 22 degrees and a dipolar magnetic field strength (B-d) of about 100 G (60 < B-d < 230 G within the 1 sigma level), below what is typically found for other magnetic massive stars. This conclusion is strengthened further by considerations of the star's X-ray spectrum. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Our results imply that both stars are expected to have a dynamical magnetosphere, so the magnetic field is not able to support a circumstellar disk. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.
Context.
The supergiant ionized shell SMC-SGS 1 (DEM 167), which is located in the outer Wing of the Small Magellanic Cloud (SMC), resembles structures that originate from an energetic star-formation event and later stimulate star formation as they expand into the ambient medium. However, stellar populations within and surrounding SMC-SGS 1 tell a different story. Aims. We present a photometric study of the stellar population encompassed by SMC-SGS 1 in order to trace the history of such a large structure and its potential influence on star formation within the low-density, low-metallicity environment of the SMC.
Methods.
For a stellar population that is physically associated with SMC-SGS 1, we combined near-ultraviolet (NUV) photometry from the Galaxy Evolution Explorer with archival optical (V-band) photometry from the ESO Danish 1.54 m Telescope. Given their colors and luminosities, we estimated stellar ages and masses by matching observed photometry to theoretical stellar isochrone models. Results. We find that the investigated region supports an active, extended star-formation event spanning similar to 25-40 Myr ago, as well as continued star formation into the present. Using a standard initial mass function, we infer a lower bound on the stellar mass from this period of similar to 3 x 10(4) M-circle dot, corresponding to a star-formation intensity of similar to 6 x 10(-3) M-circle dot kpc(-2) yr(-1).
Conclusions.
The spatial and temporal distributions of young stars encompassed by SMC-SGS 1 imply a slow, consistent progression of star formation over millions of years. Ongoing star formation, both along the edge and interior to SMC-SGS 1, suggests a combined stimulated and stochastic mode of star formation within the SMC Wing. We note that a slow expansion of the shell within this low-density environment may preserve molecular clouds within the volume of the shell, leaving them to form stars even after nearby stellar feedback expels local gas and dust.
The Chandra Carina Complex contains 200 known O- and B-type stars. The Chandra survey detected 68 of the 70 O stars and 61 of 127 known B0-B3 stars. We have assembled a publicly available optical/X-ray database to identify OB stars that depart from the canonical L-X/L-bol relation or whose average X-ray temperatures exceed 1 keV. Among the single O stars with high kT we identify two candidate magnetically confined wind shock sources: Tr16-22, O8.5 V, and LS 1865, O8.5 V((f)). The O4 III(fc) star HD 93250 exhibits strong, hard, variable X-rays, suggesting that it may be a massive binary with a period of > 30 days. The visual O2 If* binary HD 93129A shows soft 0.6 keV and hard 1.9 keV emission components, suggesting embedded wind shocks close to the O2 If* Aa primary and colliding wind shocks between Aa and Ab. Of the 11 known O-type spectroscopic binaries, the long orbital-period systems HD 93343, HD 93403, and QZ Car have higher shock temperatures than short-period systems such as HD 93205 and FO 15. Although the X-rays from most B stars may be produced in the coronae of unseen, low-mass pre-main-sequence companions, a dozen B stars with high L-X cannot be explained by a distribution of unseen companions. One of these, SS73 24 in the Treasure Chest cluster, is a new candidate Herbig Be star.
Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J175442619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 +/- 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e < 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (v(infinity) = 1500 km s(-1) in IGR J17544-2619 and v(infinity) = 700 km s(-1) in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.
We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeK alpha emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeK alpha in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeK alpha spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in gamma Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeK alpha is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeK alpha (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeK alpha is narrow (sigma(line) < 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V similar to 1000 km s(-1)). (6) The equivalent hydrogen column (N-H) directly correlates to the EW of FeK alpha, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object - wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor's wind in the FeK alpha emission and the absorption when the donor is a supergiant massive star.