EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 626 (June 2019) A&A, 626 (2019) A64 Full HTML
Free Access
Issue
A&A
Volume 626, June 2019
Article Number A64
Number of page(s) 18
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201935074
Published online 13 June 2019

© ESO 2019

1 Introduction

Discovered during a balloon flight in 1964 (Bowyer et al. 1965), Cygnus X-1 is one of the best studied black hole X-ray binaries. Based on radio parallax data, the distance of the system was measured to be at 1.86−0.11+0.12$1.86^{+0.12}_{-0.11}$ kpc (Xiang et al. 2011; Reid et al. 2011)1. Cyg X-1 consists of a (14.8 ± 1.0) M⊙ black hole that accretes from the strong stellar wind of the (19.2 ± 1.9) M⊙ supergiant O9.7 Iab star HDE 226868 (Murdin & Webster 1971; Webster & Murdin 1972; Walborn 1973; Herrero et al. 1995; Caballero-Nieves et al. 2009; Orosz et al. 2011). The star and the black hole are in a quasi-circular (eccentricity e = 0.018 ± 0.002; Orosz et al. 2011) 5.599829(16) d orbit (Webster & Murdin 1972; Brocksopp et al. 1999; Gies et al. 2003) with an inclination of i = 27.1° ± 0.8° (Orosz et al. 2011). This corresponds to a separation between the center of mass of the star and the black hole of only 42.2 R⊙, or 2.5 stellar radii. Combined with the high mass-loss rate of HDE 226868 (~10−6 M⊙ yr−1; Puls et al. 2006; Herrero et al. 1995), this small separation means that the black hole is continuously accreting material from the stellar wind, making Cyg X-1 one of the few persistent black holes in our Galaxy.

With its high equivalent hydrogen column density, NH, of around 6 × 1021 cm−2 (Dotani et al. 1997; Schulz et al. 2002a; Miller et al. 2002; Hanke et al. 2008), Cyg X-1 is associated with an X-ray dust scattering halo (see Bode et al. 1985, for the first analysis of this halo). Xiang et al. (2011) found a dust containing cloud in the interstellar medium (ISM) at ~0.885 D, where D is the distance of the system, to be the region responsible for the scattering halo.

Already soon after the identification of the optical counterpart of Cyg X-1, X-ray light curves were found to show a strong orbital modulation of the X-ray absorption column, NH, due to absorption of X-rays from the black hole in the stellar wind (Li & Clark 1974; Remillard & Canizares 1984; Bałucińska-Church et al. 2000; Poutanen et al. 2008; Miškovičová et al. 2016; Grinberg et al. 2015). Together with observations of the orbital modulation of optical lines from HDE 226868 (Gies & Bolton 1986a,b; Gies et al. 2003), these phenomena led to the picture of the stellar wind of HDE 226868 as a line driven wind or CAK wind after Castor et al. (1975, see also Friend & Castor 1982 and Morton 1967) with an asymptotic velocity of v∞ ~ 2000 km s−1 (Muijres et al. 2012). This wind is disturbed by the gravitational potential of the black hole, which leads to a focusing of the wind toward the black hole; see Friend & Castor (1982) and Gies & Bolton (1986b) for early models and El Mellah et al. (2019) for a modern treatment ofthis process. In addition, the wind is also affected by the X-rays from the compact object: strong orbital modulation of NH is seen during the canonical hard state of the black hole, where the X-ray spectrum is dominated by a Comptonized power law (Parker et al. 2015; Nowak et al. 2011; Wilms et al. 2006, and references therein). The typical bolometric luminosity of Cyg X-1 in this state is around 2 × 1037 erg s−1 (e.g., Wilms et al. 2006; Nowak et al. 1999), albeit with a large uncertainty due to our lack of knowledge of the UV spectral shape. Only very little NH modulation is seen during the thermally dominated X-ray soft state (Wen et al. 1999; Boroson & Vrtilek 2010), in which the typical bolometric luminosity is at most a factor two higher than in the hard state (e.g., Tomsick et al. 2014; Zhang et al. 1997). Optical spectra show the stellar wind to be strongly photoionized during the latter state (Gies et al. 2003, 2008).

Line driven winds are not expected to be smooth flows, but show strong density perturbations or “clumps” (Owocki et al. 1988; Feldmeier et al. 1997; Puls et al. 2006, 2008; Oskinova et al. 2012; Sundqvist & Owocki 2013). In X-ray binaries, thedensity contrast could even be further enhanced by the interaction between the wind and the strong X-rays from the compact object (Blondin 1994; Blondin & Woo 1995; Manousakis & Walter 2011, 2015, and references therein). For Vela X-1 and Cyg X-1 it has been estimated that more than 90% of the wind mass is contained in less than 10% of the wind volume (Sako et al. 1999; Rahoui et al. 2011). When the line of sight to the compact object passes through one of these clumps, X-rays are absorbed by the moderately ionized material in the clump, leading to a so-called dipping event. This is also observed for other sources (Hemphill et al. 2014; Grinberg et al. 2017). During the hard state of Cyg X-1, such short-term dipping events are observed predominantly during the upper conjunction of the black hole, i.e., when the line of sight passes through the densest region of the stellar wind and is most likely to pass through a clump (Li & Clark 1974; Mason et al. 1974; Parsignault et al. 1976; Pravdo et al. 1980; Remillard & Canizares 1984; Kitamoto et al. 1984; Bałucińska-Church et al. 2000; Feng & Cui 2002; Poutanen et al. 2008; Hanke et al. 2009; Miškovičová et al. 2016; Grinberg et al. 2015). The precise structure of the clumps, i.e., their density and ionization structure, is unknown. Most recent 2D simulations of such a stellar wind show a very complex evolution of velocity and density structures with the formation of characteristic small-scale clumps of various shapes embedded in areas with lower density (Sundqvist et al. 2018). Sundqvist et al. (2018) have found a typical clump mass of 1017 g and an average clump size of 1% of the stellar radius at a distance of two stellar radii. These results qualitatively confirm earlier theoretical models (e.g., Oskinova et al. 2012; Sundqvist & Owocki 2013) and observations (e.g., Grinberg et al. 2015). See also the review paper by Martínez-Núñez et al. (2017).

In this paper we present the first detailed high-resolution study of absorption dips in Cyg X-1 using time-resolved X-ray spectroscopy during dips observed with the High Energy Transmission Grating Spectrometer (HETGS) on board Chandra. The analysis of a series of Chandra-HETGS observations of the source taken during the low/hard state allows us to investigate the dipping mechanism and to probe the ionization state of the absorber directly. This is the third paper in a series devoted to a study of Cyg X-1 using Chandra high-resolution spectroscopy. Our previous analyses addressed the stellar wind in Cyg X-1 in the hard state as seen outside the absorption dips, i.e., we analyzed the hot tenuous phase of the wind. In the first paper of the series we investigated the wind at ϕorb ~ 0 (Hanke et al. 2009, hereafter Paper I). We extended this work to the orbital modulation of the wind using non-dip observations at ϕorb ~ 0.0, ~0.2, ~0.5, and ~0.75 (Miškovičová et al. 2016, hereafter Paper II). These papers provide us with a reference for the “normal” (non-dip) hard-state spectrum of Cyg X-1.

In this paper, we address the spectral changes during dipping episodes. We discuss our data reduction and identification of the dipping events in Sect. 2. In Sect. 3 we study the spectral evolution of the silicon and sulfur line regions from the non-dip to the deepest dip and we show that absorption lines of lower ionized species appear as the line of sight crossesthrough denser regions of the absorbing clump. During the deepest dips the line of sight is fully blocked, revealing line emission from the photoionized plasma around the black hole. We summarize our results in Sect. 4.

2 Data reduction

Out of the Chandra-HETGS observations discussed in Paper II for their non-dip properties, we selected a subsample of three observations(Table 1), namely ObsIDs 3814, 8525, and 9847, which show distinct dipping episodes in their light curves (Fig. 1, Sect. 2.1). As discussed in Paper II, the other ObsIDs available were taken in the intermediate and soft states and thus the material in the system was in a different radiative environment. All observations selected for further analysis of the dips were performed in timed exposure (TE) mode. In this mode the usual frame time is 3.2 s exposure before the data are transferred into a frame store for readout. Cyg X-1 is a bright source even in the low/hard state. Therefore, for the discussed observations only a half array of the chips (512 CCD rows) was read out to minimize pileup by reducing the frame time to 1.7 s. We re-extracted the data using the Chandra Interactive Analysis of Observations (CIAO) software version 4.6 with all parameters set to default except for the HETG cross dispersion extraction width (width_factor_hetg parameter in CIAO), which we reduced from 35 to 10 for better coverage at the short wavelength end.

For the light curves and spectral analysis we used the first order spectra of the high and medium energy gratings (HEG, MEG; Canizares et al. 2005). Because of the excellent background discrimination of the data extraction process (order sorting), any remaining background is negligible compared to the bright source. We therefore did not subtract any background from the final spectra. All further data analysis was performed with the Interactive Spectral Interpretation System version 1.6.2 (ISIS; Houck & Denicola 2000; Houck 2002; Noble & Nowak 2008).

2.1 Light curves

The light curves of ObsIDs 3814 (~48 ks) and 8525 (~30 ks) slightly overlap in orbital phase (ϕorb = 0.93–0.03 and ϕorb = 0.02–0.08) and show avery similar morphology of strong absorption dips. ObsID 9847 (~19 ks) at ϕorb = 0.17–0.21 was taken within the same binary orbit as ObsID 8525; it contains absorption dips as well, although they are less pronounced than at ϕorb ~ 0. Dips can last from several seconds to more than 10 min (Kitamoto et al. 1984). We extracted light curves with a 25.5 s resolution to uncover the shorter dips as well. For a detailed description of the data set and a detailed analysis of the non-dip spectrum of individual observations, we refer to Papers I and II. Paper II also discusses ObsID 11044 (~30 ks), which at ϕorb = 0.48–0.54 for the first time provides high-resolution spectroscopy of Cyg X-1 during lower conjunction, i.e., when our line of sight passes outside of the focused wind. As this light curve appears to be virtually free of dipping, we do not discuss this light curve in detail in this work. Cyg X-1 was in a comparable hard state during ObsIDs 3814, 8525, and 9847 according to the Rossi X-ray Timing Explorer (RXTE) All Sky Monitor classification by Grinberg et al. (2013); see also the more detailed discussion of the source state in Paper II. We do not use ObsID 3815 as it shows much less dipping than the other hard state observations, which does not allow for a distinction between different dip stages (Paper II, see also Miškovičová et al. 2011).

thumbnail Fig. 1

Top panel: light curves with a 25.5 s time resolution of all three observations as a function of orbital phase. Bottom panel: corresponding hardness ratios. All observations were in the low/hard state. Dips are strongest at ϕorb ~ 0.0, become weaker at ϕorb ~ 0.2 and ϕorb ~ 0.75, and completely fade at ϕorb ~ 0.5 (see also Paper II). Colors indicate the selected dipping stages (see right panel of Fig. 3 for the color denotation and Sect. 2.3 for details on the data selection).
Table 1

Chandra-HETGS observations used in this paper.

2.2 Color–color diagrams

As we discuss in more detail below, the dips are due to absorption events caused by material in the line of sight to the primary source of X-rays. Dipping events such as thos seen in Cyg X-1 are transient events characterized by quickly varying spectral shape. This strong variability complicates the spectral analysis. Ideally, we would want to study how the spectral shape – photon index, absorption/emission lines, and continuum absorption – varies with time, but the limited signal to noise of our observations renders this impossible. We therefore have to resort to some kind of averaging technique in which we extract spectra from time intervals where we believe that the spectral shape is at least representative for a given part of a dip. We find these time intervals by looking at the time resolved spectral behavior of the source as represented in so-called color-color diagrams.

Figure 2 shows how absorption affects an observed primary power law continuum with Γ = 1.73 (typical for Cyg X-1 in the hard state, e.g., Grinberg et al. 2013) that is fully covered by material with a fixed column density NH(0)$N_{\mathrm{H}}^{(0)}$ (e.g., the outer parts of the stellar wind or absorption in the ISM). For this paper, we set NH(0)$N_{\mathrm{H}}^{(0)}$ to the fit result of Paper I, NH(0) = 5.4 × 1021 cm−2$N_{\mathrm{H}}^{(0)}\,{=}\,5.4\,{\times}\,10^{21}\,\mathrm{cm}^{-2}$. In most astrophysical sources, the structures responsible for dipping do not cover the whole primary X-ray source, rather a partial coverer is present in these systems that covers a fraction fc of the source with a column NH (i.e., 1 − fc remains uncovered). A possible physical picture for such a partial coverer would be an optically thick cloud that is smaller in (angular) size than the X-ray source, or a cloud that passes very quickly over the X-ray source, covering the source only for part of the integration time. The observer therefore sees the sum of the uncovered spectrum (dash-dotted line in Fig. 2) and the covered spectrum (dashed lines). The summed spectrum is shown as solid lines in Fig. 2. The left panel of Fig. 2 shows the observed spectral shapes for a constant covering factor fc = 90% and for varying optical depths of the covering medium (up to 128NH(0)$128N_{\mathrm{H}}^{(0)}$). As NH is increased, the partial coverer removes most of the flux at soft energies and only the direct component remains visible.

Because absorption events are typically of rather short duration, the detailed spectral shape is often not directly observable. It is, however, possible to characterize the spectral shape using X-ray colors or hardness ratios (see, e.g., Hanke et al. 2008; Nowak et al. 2011, for similar approaches). In this work, we define the X-ray hardness ratio as the ratio of the count rates in two energy bands. In order to be consistent with Paper I we define the ratio such that its value increases as the spectrum softens, so technically this ratio is a “softness ratio”. We calculate hardness ratios using the fluxes measured in three energy bands – denoted as A (0.5–1.5 keV ≡ 24.8–8.27 Å), B (1.5–3.0 keV ≡ 8.27–4.13 Å), and C (3.0–10.0 keV ≡ 4.13–1.24 Å) with the ranges consistent with Paper II and the previous work of Nowak et al. (2011) – to characterize the spectral shape in this way.

The right-hand panel of Fig. 2 shows the locus of such colors for the spectral shape discussed above and several different covering fractions. For a constant fc, a characteristic track in the color–color diagram becomes apparent: At low NH~NH(0)$N_{\mathrm{H}}\sim N_{\mathrm{H}}^{(0)}$, the covered component dominates the spectrum and the source is found in the top right of the diagram. Since photoabsorption first influences the soft bands, increasing NH decreases both, A∕B and B∕C, and the source moves toward the bottom left of the diagram. At intermediate NH, the A band is dominated by the (constant) uncovered fraction, but increasing NH still decreases B and C. For this reason, A∕B starts to increase again, while B∕C continues to decrease, such that the track of the source in the color-color-diagram starts to turn toward the right. When the contribution of the covered component to both A and B bands is almost negligible, A∕B remains constant, while B∕C increases with increasing NH. This behavior leads to a horizontal track in which the source moves to the right in the color–color diagram. Finally, for the largest NH the entire covered fraction is removed from the observable spectrum such that we expect the X-ray color to move asymptotically back to its unabsorbed value.

thumbnail Fig. 2

Left panel: effect of partial absorption with increasing NH at constant covering fraction with NH(0) = 5.4 × 1021 cm−2$N_{\mathrm{H}}^{(0)}\,{=}\,5.4\,{\times}\,10^{21}\,\mathrm{cm}^{-2}$ after Paper I. Right panel: color–color diagram showing the tracks for different covering fractions and column densities of the second absorber. See text for details.

2.3 Dip selection

We now apply the ideas discussed above to Cyg X-1. Figure 3 shows the hardness ratios measured during the three observations considered in this work. The data show a behavior similar to that expected from partial covering: data outside of dips (orange data points) are barely absorbed; as the source enters deeper dips (bluer colors), it follows a track that resembles the theoretical tracks shown in Fig. 2, including a turning point where the soft color becomes softer, while the hard color barely changes.

To study the spectral shapes during different phases of the dips in more detail, we extract data from time intervals that correspond to various stages of dipping, i.e., phases of similar spectral shape. A direct comparison of the theoretical tracks derived from the simple partial covering model, however, shows deviations between the data and the tracks. These deviations could be indicative of changes in both the covering fraction and NH, ionization effects and because the underlying X-ray continuum is more complex than a simple absorbed power law. A detailed treatment of these effects, however, would require us to fully understand the ionization structure of the absorber, which is beyond the scope of this paper. We therefore determine the shape of the track empirically by fitting an empirical curve to the scatter plot. This curve is described through a parameterized polynomial of second degree for each of the two colors in the diagram. To find the polynomial coefficients, a χ2 minimization algorithm is used to optimize the shortest distance of each data point to the curve.

To select data from the dips, we first remove the non-dip data as defined in Paper II2. This ensures consistency with the earlier results from Paper II. For the remaining data, corresponding to various degrees of dipping, the goal is to find the highest possible resolution in terms of number of dipping stages, while maintaining a good enough signal-to-noise ratio (S/N) to be able to constrain the spectral fits well. To find this balance, we start out with a large number of 12 segments, which are chosen such that each slice contains roughly the same number of counts. Then, starting at the deepest dipping stage (bottom right-hand corner of the color-color diagram) we successively combine these small segments until the signal-to-noise ratio of the resulting spectra is sufficient to detect all possible lines of the Si and S series, but keep the number of segments for this as low as possible to be able to distinguish a greater number of dip stages. These combined segments constitute the selection for the deepest dipping stage for further analysis. Subsequently, the selection of the next dipping stages follow the same approach until all 12 segments are sorted.

For all observations we obtain the best results by defining three dipping stages in addition to the non-dip phase, each consisting of four of the smaller segments. Consequently, the three stages each have roughly the same number of counts within an observation. The four dipping stages for each observation are classified as “non-dip”, “weak dip”, “dip”, and “strong dip” (see Fig. 3 and Table 2 for the count rates and exposure of each dip stage).

thumbnail Fig. 3

Color–color diagrams for all three observations, showing hardness ratios calculated for soft, intermediate, and hard X-ray bands. As discussed in the text, the observed tracks are determined by absorption: following the track from the top right corner,representing the persistent flux, to the bottom left corner, showing dips, absorption increases and therefore both colors become harder. During the deepest dips, the softer color becomes softer again, while the hard color does not change. Thisbehavior can be explained by partial covering (Sect. 2.2). The spectrum hardens considerably during thesedips. The colors denote the four different dip stages (see right panel). The black line indicates our polynomial fit to the track; see Sect. 2.3 for details.
Table 2

Exposure time, texp, and dip stage averaged count rates for the total observations and silicon (1.605–2.045 keV ≡ 7.725–6.063 Å) and sulfur (2.295–2.7 keV ≡ 5.402–4.6 Å) bands for all observations.

3 Spectral evolution from non-dip to dip

Figure 4 shows the evolution of the spectrum during dipping for ObsID 8525. In addition to the general change in spectral shape due to photoelectric absorption, we see strong changes between the individual dipping stages in the region between 1.6 and 2.7 keV (7.7–4.6 Å), where absorption lines of silicon and sulfur ions are the most prominent spectral features (box in Fig. 4, see also Miškovičová et al. 2011). We therefore concentrate on these Si and S lines.

3.1 Si and S regions

For a detailed analysis of the Si and S absorption lines, we take into account the ± 1 order HEG and MEG spectra. We combine both spectra for the HEG and the MEG, and fit the spectra without rebinning them. For clarity reasons, all figures show the combined HEG and MEG spectra, where the higher resolution of the HEG is rebinned to the lower resolution of the MEG. This combination of HEG and MEG, however, only applies to the display of the data, not to the actual fitting.

As the emphasis of this paper is on the behavior of the Si and S lines, we do not attempt to model the broadband Chandra continuum or the continuum absorption, but rather describe the local spectra in the Si (1.605–2.045keV ≡ 7.725–6.063 Å) and the S (2.295–2.7 keV ≡ 5.402–4.6 Å) regions. These regions are narrow enough in energy that the curvature due to the absorption is negligible. We can therefore describe the local absorbed continuum by a simple power law; obviously, the photon indices of the continuum are different in both bands.

A further complication is that as a consequence of the source brightness the non-dip, weak dip, and to a lesser extent also the dip spectra are affected by slight pileup. The major effect that pileup has on our narrowband data is a change in count rate, which is seen as a change in the relative flux normalization of both spectra, as pileup affects the MEG more strongly than the HEG because of the significantly lower spectral resolution3. We can compensate for this effect by including a multiplicative, detector dependent constant in the spectral modeling. Because of the narrow energy bands considered here, a further pileup correction is not necessary. The pileup fraction is a slowly changing function of energy such that it does not appreciably affect the equivalent width of the lines considered in this work; because of the nature of X-ray gratings, line energies are not affected. See Hanke et al. (2009) for a discussion of pileup effects in the non-dip spectrum of Cyg X-1.

The local continuum being accounted for (see Table 3 for parameters), we model the absorption lines in both regions using additive Gaussian line profiles (Fig. 5). These profiles give a slightly better description of the line shapes than Voigt profiles. For the absorption lines, the line width σ was frozen at 1.2 eV (silicon) and 1.8 eV (sulfur), well below the detector resolution. Lines were identified using energies from Hell et al. (2016), Porquet et al. (2010), and from AtomDB, version 2.0.04. We find blended line complexes for L-shell ions, He w lines for He-like ions, and Ly α lines for H-like ions of both silicon and sulfur. In addition to the K α absorption line series, there is the forbidden He-like Si XIII z line in emission (1.8394 keV ≡ 6.7405 Å, 1s 2s 3S1 → 1s2 1S0; Porquet et al. 2010); this line is sometimes also called the He-like Si XIII f line. We note that there is a degeneracy between the Be-like Si XI and Li-like Si XII absorption lines at 1.8275 keV ≡ 6.7844 Å and 1.8450 keV ≡ 6.7200 Å (both energies from Hell et al. 2016) and the forbidden He-like Si XIII z emission line. This increases the uncertainty of the respective fit values. Finally, in the strong dip spectrum of ObsID 9847, there are indications for a O-like S IX K α line, which appears once the HEG data are rebinned to match the (unbinned) MEG grid and both are combined for display. This line is too weak for us to be able to constrain its parameters such that we only claim a tentative detection of O-like S IX. The valuesfor central line energies and equivalent widths of all detected silicon and sulfur absorption lines can be found in Tables 4 and 5. We note that the lower ionization stages (N-like and O-like ions) only appear deeper in the dips. This hints toward a shell-like ionization structure of the clumps with a hot, highly ionized surface and cooler and less ionized parts in the core.

In addition to these lines from S and Si we also find an Al XIII Ly α5 absorption line (1.7285 keV ≡ 7.1729 Å), Mg He γ (1.6591 keV ≡ 7.4730 Å, 1s 2 1 S0–1s4p 1P1) and He δ (1.6961 keV ≡ 7.3100 Å, 1s 2 1 S0–1s5p 1P1) absorption lines, and we detect Mg XII Ly β (1.7447 keV ≡ 7.1063 Å) and, in some cases, a hint of Mg Ly δ (1.8843 keV ≡ 6.5799 Å) blended with an Fe XXIV transition (1.8851 keV ≡ 5.5771 Å, 1s 2 2s 2S1∕2–1s27p 2P3∕2). These lines imply the presence of Mg Ly γ at 1.840 keV ≡ 6.738 Å, i.e., blended with He-like Si z. We expect Mg Ly γ to be very weak but it may slightly contaminate our measurements of the Si He z and the Li-like Si XII line, thus further contributing to the uncertainty of the fit values.

thumbnail Fig. 4

Rebinned HEG spectra of the different dipping stages of ObsID 8525. For clarity, the non-dip spectrum (orange, S/N ≥ 13) is scaled by afactor of 8, the weak dip spectrum (red, S/N ≥ 15) by a factor of 4, and the dip spectrum (violet, S/N ≥ 15) by a factor of 2. The strong dip spectrum (blue, S/N ≥ 17) is not scaled. The box denotes the region containing the Si and S lines. The color scheme is the same as in Figs. 1 and 3.

3.1.1 Equivalent widths or line strengths and column densities

The equivalent widths of the silicon and sulfur lines vary around 1 eV and have rather large error bars for the weaker lines (Table 5). We note that the He-like Si XIII z emission line additionally increases the uncertainties of the Si XI and Si XII lines. Comparing the lines in the different dipping stages reveals that the line strengths of the low charge states – and thus the respective column densities – increase with dipping, while those of the higher charge states decrease (Fig. 6). This shift in charge balance alone is already strong evidence supporting the theory of a highly structured medium where clumps of colder, denser material are embedded in a highly ionized plasma (Oskinova et al. 2012; Sundqvist & Owocki 2013, and references therein).

For optically thin absorption lines, the equivalent width directly translates into the column density of the parent ion, Ni. In this section and in the following we use the convention that the equivalent width of absorption lines is negative and that of emission lines is positive. We derive the column densities for all Si and S ions visible in our spectra from (Paper I, Eq. (11)) Ni = mec2|Wλ|πe2fijλ02 = 9.11 × 1015 cm−2fij |WE|eV,\begin{equation*} N_i\,{=}\,\frac{m_{\mathrm{e}}c^2 |W_{\lambda}|}{\pi e^2 f_{ij} \lambda_0^2}\,{=}\,\frac{9.11\,{\times}\,10^{15}\,\mathrm{cm}^{-2}}{f_{ij}}\,\frac{|W_{\textrm{E}}|}{\mathrm{eV}}, \vspace*{-3.6pt}\end{equation*}(1)

where Wλ and WE are the equivalent width in wavelength and in energy space, and where fij is the oscillator strength, taken as the sum of all transitions blending into the absorption line of the respective ion. We use the compilation of Verner et al. (1996) for H-like ions, and the compilation of Palmeri et al. (2008) for all other ionization stages. These oscillator strengths include Auger damping, which is important at the densities expected. All other symbols have their usual meanings.

Summing the column densities for all ions of an element then yields a lower limit for the total column density of the element. For comparison with the continuum fitting values, the abundance of the element can then be utilized to convert this column density into a corresponding NH value. The values listed in Tables 6 and 7 list two NH values. One value is based on the solar abundances as summarized by Wilms et al. (2000). In our second conversion we consider that in theiranalysis of the optical spectrum of HDE 226868 Herrero et al. (1995) found that helium is overabundant by a factor of 2.53 with respect to the solar value. Assuming that metals scale with the same factor as helium with respect to hydrogen, we scale the solar abundances of Wilms et al. (2000) by this factor and use these corrected abundances to derive the corresponding hydrogen column.

Figure 6 shows how the single ion column densities vary with depth of the dip. As we enter the dip, for the lowest ionization states (Fig. 6, bottom panels) the columns increase toward the deepest dip stages. This trend reverses for the highest ionization states (Fig. 6, top panels). This behavior is also indicative of an absorber with a layered ionization structure, i.e., a medium whose outer regions are more highly ionized than the central core region.

For the non-dip data, column densities of H-like and He-like silicon have already been presented in Paper I. Despite the fact that our dip selection criterion for the non-dip data has not changed, there is a slight difference in the equivalent width determined for He-like Si XIII, where Paper I finds an equivalent width of − 1.73 eV, while we find − 1.92 eV, resulting in a ~10% difference inthe column. Both values are still in agreement within their uncertainties. The difference in equivalent width for H-like Si XIV, on the other hand, is a surprising factor 1.3. A possible explanation is that the non-dip data allowed a global fit, including interstellar absorption, more than a hundred lines, and a pileup correction. The three dip stages, however, do not have enough signal to constrain such a model, and consequently for consistency we also model the non-dip data with local fits. Differences at the 15% level would therefore be expected owing to these different analysis approaches. Taking into account these additional 15%, the H-like Si XIV columns are also consistent within their uncertainties.

Figure 7 shows the column densities as a function of the ionization potential of the respective ions. The closed shell ions He-like Si XIII and S XV exhibit enhanced column densities as expected. However, whereas the appearance of lower ionization stages in the strong dip points to a lower temperature in the core of a clump, it is also possible that we see a constant temperature clump with a higher density in its core. As the absorber is ionized, the column density measured from continuum absorption is only a lower limit; we cannot see the fully ionized material. The figure also clearly illustrates the large column of the hydrogenic lines. It is therefore very likely that a significant amount of the S and Si in the system are fully ionized and thus cannot be detected (see also Sect. 3.1.3).

Table 3

Detector constants and continuum parameters for all observations.

thumbnail Fig. 5

Evolution of silicon (left column) and sulfur (right column) lines with dipping stages for the three observations considered (from top to bottom). The panels show combined HEG and MEG ± 1st spectra relative to their respective continua. The solid black line shows the line fits.
Table 4

Observed silicon and sulfur line centers for all observations.

Table 5

Evolution of the equivalent widths of the silicon and sulfur lines during the dipping stages for all observations.

thumbnail Fig. 6

Column densities for each dip stage derived from the respective equivalent width for all Si and S ion absorption lines detected in ObsIDs 3814 (yellow diamonds), 8525 (green triangles), and 9847 (blue squares). Negative values of the equivalent width correspond to absorption lines. The column densities of lower charged ions increase with increasing dip depth, whereas column densities of highly charged ions decrease with dip level. Be-like Si XI and Li-like Si XII lines are contaminated by the Si He z emission line.

3.1.2 Doppler shifts

To analyze the morphology of the absorbing material, we next take a look at the bulk motion, measuring the Doppler shifts of the lines in the Si and S regions for each dipping stage. Since the uncertainty of the theoretical rest-wavelengths is on the order of the expected line shifts (Hell et al. 2013), we measured the centroids of the line blends for each ion in the laboratory using an electron beam ion trap and a microcalorimeter (Hell et al. 2016). For the shifts of the He- and H-like 1s → 2p transitions, we used the tables of Drake (1988) and Garcia & Mack (1965), respectively, as reference. These are the same reference publications that we also used for the calibration of the laboratory data.

Figure 8 shows the resulting velocities for each line in the different dipping stages of each observation. Overall, for the samespectral lines the Doppler shifts are consistent throughout the different dipping stages. This is more obvious for the uncontaminated lines in the Si spectra, where statistics are much better than for the S lines and line centers can be determined more accurately. The O-like S IX line detected in ObsID 9847 is the only outlier. As the width of this line is narrow (only one energy bin), and since we only claim a tentative detection of this line (see Sect. 3.1), the derived Doppler velocity is probably inaccurate and we do not include the line in Fig. 8.

Figure 9 shows the distribution of the Doppler shifts in the three observations. For ObsID 3814, the velocities scatter around 0 km s−1. For ObsID 8525, the velocities increase toward negative values and scatter around − 60 km s−1, and for ObsID 9847, the scatter moves down to around − 230 km s−1.

This Doppler shift cannot be due to the free fall velocity of donor material onto the black hole. As we look at the system near ϕorb = 0, the free fall velocity would cause a redshift, not a blueshift of lines as seen in the spectra. The Keplerian velocity of the black hole cannot be the cause for the energy shift either, as we would expect the maximum amplitude of the shift to be on the order of 100 km s−1 and not much larger than that; the projected semi-amplitude of the black hole is KBH ~ 91 km s−1 (e.g., Brocksopp et al. 1999; Gies et al. 2003; Orosz et al. 2011).

Table 6

Column densities as calculated from the silicon line equivalent widths.

Table 7

Column densities as calculated from the sulfur line equivalent widths.

3.1.3 A clumpy absorber

The consistency of the Doppler shifts of the lines during dipping provides evidence for inhomogeneities in the absorber that are forming a unified structure that has some kind of density stratification. We call these structures clumps, even though we do not want to imply that the structures are static phenomena as opposed to transient structures. Either owing to self shielding, or because the material is in some kind of approximate pressure equilibrium similar to the clouds in active galaxies (Krolik et al. 1981), the part of the clump facing away from the black hole is less ionized. As the clump moves through the line of sight (Fig. 10), we first look through its outer, more ionized regions. We only see the lower ionization stages during the deeper dip where the middle of the clump is located in our line of sight.

Regarding the orbital phase dependency, from the three ObsIDs analyzed we can already see a trend toward higher Doppler velocities further away from ϕorb ~ 0. This is consistent with Paper II. Figure 11 shows the median, mean, first, and third quartile of the velocity distribution of the Si and S lines in comparison to the non-dip results shown in Fig. 11 of Paper II. The phase dependence of the Doppler shifts measured within the dips as well as the amplitude are roughly consistent with Paper II. The material producing the dips thus seems to be at the same distance from the black hole as what can be seen during non-dip.

The projected wind velocities in this figure are calculated using a CAK-model with a focused wind based on Gies & Bolton (1986b, see Paper II for more details6). The model consists of a radially symmetric wind that is focused onto the black hole in a cone of 20°. Gies & Bolton (1986b) stated that their model is valid only out to 0.9d away from the stellar center, where d is the distance between the star and the black hole, as the density at this distance is already too small to produce a significant contribution to the absorption line profiles. It is possible to argue that this is only the case for a smooth wind and clumps have to be treated differently. Even for ϕorb = 0 the regions we are interested in (distances of 0–0.75d from the black hole, in the direction of the line of sight at different phases) are mostly outside of the valid region of the model of Gies & Bolton (1986b), but as our observations are near ϕorb = 0, where the distance to the star is not much more than 0.9d, we still take the model for a rough estimate. However, both measured non-dip and dip Doppler shifts fit neither the projected wind velocity (difference in phase) nor the black hole (difference in amplitude). Thus, we cannot make an exact statement on the location of the clumps in the system from our data, although a comparison of the amplitude implies a distance of less than 0.25 d from the black hole.

In orderto investigate the origin of the observed ionization, we use XSTAR (Kallman & Bautista 2001), which models the photoionization structure of spherically symmetric gas clouds that are irradiated by a central source of arbitrary spectral energy distribution. We assume gas clouds with the default elemental composition used in XSTAR, which is based on Grevesse et al. (1996). The clouds have a covering fraction of 100%, a temperature of 106 K, a constant total hydrogen particle density of n = 1010 cm−3, and are irradiated by an incident 1–1000 Ry luminosity of L = 1.4 × 1037 erg s−1. In this configuration, the definition of the ionization parameter is ξ = L/(nR2)$\xi\,{=}\,L/\left(n R^2\right)$, where R is the inner radius of the shell. Our incident spectral energy density (SED) is adapted from Pepe et al. (2015) and assumes that the source of X-rays is point like. Any mechanical input from the jet in the system is ignored. In order to gauge the influence of the strong UV photon field of the donor star, we perform simulations with and without a stellar contribution to the total SED (see Fig. 12).

We calculate model grids by varying ξ and the total hydrogen column, NH = nℓ, where ℓ is the geometric size of the cloud. Since L and n are fixed, this approach effectively models photoionized clouds of varying sizes and distances from the source of ionizing photons. The left panel of Fig. 13 shows the columns for the ions of Si and S as a function of log ξ for NH = 1022 cm−2. There is no single value of logξ that produces the measured columns for all ions under consideration. This result holds for all values of NH. Based on a grid of XSTAR simulations in which we vary logξ and NH, in the right-hand panel of Fig. 13, we show the contours in the log ξ-NH-plane in which the theoretical column of each ion of Si and S ion equals that measured during the strong dip of observation 8525. In the ideal case ofa simple medium of constant density and temperature, all of these contours would intersect in a single point; the modulo measurement errors are, however, small compared to the large dynamic range shown in Fig. 13. We note that the inclusion of the optical companion has only a marginal effect on the observed ion column densities of Si and S (Fig. 13, right, panels c and d). This is important because the binary geometry of the system with two sources of radiation cannot be accounted for by the XSTAR simulation. Neglecting the optical companion therefore only introduces a small systematic effect on ionization structure modeling of these ions.

In practice, assuming that the dip is produced by a single absorbing cloud, we would expect some kind of a density variation in the cloud, as discussed above. We can use our simple XSTAR simulations to identify regions of particular interest. From Fig. 13, columns of moderately ionized ions can be roughly reproduced with hydrogen column densities around 3 × 1020 cm−2 and ionization parameters of log ξ ~1–2. The high ionization states of silicon and sulfur, on the other hand, require column densities and an ionization parameter that are roughly an order of magnitude higher. This result points at a more complex origin for the ionization structure than photoionization alone. One possible explanation would be ionization due to the strong shocks that are present in hydrodynamic simulations for massive X-ray binary systems (Blondin et al. 1990, 1991; Blondin 1994; Manousakis & Walter 2011, 2015; Sundqvist et al. 2018, and references therein) and in winds from early-type stars (Owocki et al. 1988; Sundqvist & Owocki 2013, and references therein).

thumbnail Fig. 7

Column densities vs. ionization potential for the respective silicon and sulfur ions in all dip stages.
thumbnail Fig. 8

Doppler velocities for the various ionization states of Si VIII–Si XIV (empty) and S X–S XVI (filled) in all dipping stages (non-dip: orange circles; weak dip: red triangles; dip: violet squares; and strong dip: blue stars) over all observations in the different orbital phases.
thumbnail Fig. 9

Distribution of velocity shifts of all Si and S lines detected in the three ObsIDs with their medians. The peak velocity clearly shifts between the observations.
thumbnail Fig. 10

Sketch of a clump passing our line of sight and causing the different dip stages. Red color indicates more highly ionized regions; blue color regions indicate lower ionization balance.
thumbnail Fig. 11

Radial velocities of the black hole (red dashed line) and the companion (orange dashed line), and wind velocities projected onto the line of sight as a function of distance from the black hole. Distances are given in units of the distance between the black hole and its donor, d. Colored boxes indicate the range of measured non-dip Doppler shifts from the Chandra observations (red: ObsID 3815; yellow: ObsID 3814; green: ObsID 8525; and blue: ObsID 9847, modified and corrected from Paper II, Fig. 11). Symbols show the median (star), mean (circle), first (point-up triangle), and third quartile (point-down triangle) of the distribution of Doppler shifts for the Si and S lines.
thumbnail Fig. 12

Incident SED used for the XSTAR simulations. The dotted and dashed lines show the SED with and without the optical companion star, respectively. We note that the SED is always renormalized to the assumed model luminosity.

3.1.4 Line symmetry

Although a detailed analysis of the complex photoionization line profiles observed in Cyg X-1 is beyond the scope ofthis paper, we can at least take a quantitative look at the line shapes in the region of ± 8 bins around the respective rest energy, adjusted for the most probable Doppler shift for each ObsID (3814: 0 km s−1, 8525: − 80 km s−1, 9847: − 320 km s−1; see also Fig. 9). Figure 14 shows the original data (colored) and the same data flipped (black) around the center of the bin containing the adjusted rest energy (gray dashed line). For a symmetric line shape, the original and flipped data should match, as is the case, for example, for 8525 He-like Si XIII in the dip data, or 3814 H-like Si XIV non-dip, weak dip and dip data. However, often we see the lines shifted by one or even two bins (e.g., 3814 C-like Si IX strong dip or 9847 B-like S XII weak dip, dip and strong dip). There are also signs of asymmetry (e.g., 9847 C-like Si IX strong dip or 3814 C-like S XI dip), as well as lines thatshow both asymmetry and a shift (e.g., 3814 C-like Si IX strong dip or 8525 B-like Si X non-dip and dip).

Clear P Cygni-like profiles for Cyg X-1 were reported only at ϕorb = 0.5 during ObsID 11044 (Miškovičová et al. 2011). Schulz et al. (2002a) also found indications of P Cygni profiles at ϕorb = 0.74. In the UV, P Cygni profiles were found by Vrtilek et al. (2008), also at ϕorb = 0.5, where the stellar wind is focused toward the observer by the black hole. At around ϕorb = 0, the focused wind is moving away from the observer toward the black hole, and the redshift and ionization are highest. Asymmetries in the lines arethus expected at that phase. This was observed, for example, by Feng et al. (2003) and discussed in Paper II. Our observationsfit this picture.

3.2 Emission in the deepest dips

While the absorption lines discussed so far show a clear dependency on the depth of the dip, this is not the case for the He-like Si XIII z emission line, whose total line flux is constant throughout the dip stages. This result could either indicate that the emitting material originates in a geometrically much larger region than the absorber, such that the emission line is not affected by the absorber, or, alternatively, that the emission line originates in an area where the line of sight does not pass through the absorber.

A possible origin for the Si XIII emission line could be the photoionization region in the stellar wind that surrounds the X-ray emission region. As discussed in Paper I, the high ionization region around the X-ray source, similar to the Strömgren sphere around stars, is comparable in size to the separation of the black hole and its donor star, i.e., it is much larger than the absorbing cloud. Given the low optical depth of this region we expect a region whose X-ray emission is dominated by emission lines, similar to that seen, for example, in Vela X-1 (Schulz et al. 2002b; Watanabe et al. 2006; Grinberg et al. 2017) or in Cen X-3 (Wojdowski et al. 2003). The flux from the photoionization region, however, is much fainter than the X-rays from the accretion flow. The region has therefore been mainly studied for edge-on systems, where our line of sight onto the bright X-ray source is blocked by the donor star during eclipses. While the Cyg X-1/HDE 226868 system does not show eclipses, we can utilize the strong dipping to at least partially block the line of sight to the black hole. The best opportunity for this is presented by the long, deep dip at the beginning ObsID 8525, which is the most pronounced of all the dips in our observations. The deep part of this dip (according to our dip stage identification, i.e., the blue data points in Fig. 1 belonging to that dip) lasted for 5.3 ks (1.5 h), which is long enough to accumulate a decent gratings spectrum of one single dip. Because of the strong absorption, pileup is no concern in the resulting spectrum and the HEG and MEG spectra agree with each other.

As shown in Fig. 15, during the deepest parts of the dip in ObsID 8525, the photoionized zone around the black hole indeed dominates the spectrum, resulting in the presence of emission lines. The most prominent lines detected are Si XIII He z (1.8394 keV ≡ 6.7405; Porquet et al. 2010), Mg Ly α (1.4723 keV ≡ 8.4211 Å), the Mg He α triplet (~1.3522 keV ≡ 9.1691 Å), as well as Ne Ly β (1.2109 keV ≡ 10.2390 Å), Ne Ly α (1.0280 keV ≡ 12.0601 Å), and the Ne IX He α triplet (~0.9220 keV ≡ 13.4473 Å). Unfortunately, however, the signal to noise in the lines is too low to allow a further characterization of their parameters.

thumbnail Fig. 13

Left panel: predicted columns for the ions of Si and S for a spherical cloud with NH = 1022 cm−2 that is irradiated by a X-ray source with the spectrum defined by Fig. 12. Right panel: lines of constant ion column density color-coded for each different ion of Si (panels a and b) and S (panels c and d), for XSTAR simulated values of equivalent hydrogen column density, NH, and ionization parameter, logξ. The assumed ion columns are the values measured during the strong dip of observation 8525 (Table 5). If a simple XSTAR model described the observation, we would expect these lines to cross in one point (within their error bars). The left column of the figure shows calculations without the effect of the optical companion (OC); the right column includes the UV radiation from the companion.
thumbnail Fig. 14

Original line profiles (colored) for the non-dip (top, orange), weak dip (second from top, red), dip (third from top, purple), and strong dip (bottom, blue) lines for selected ions, as well as the same data flipped (black, see also text for further information) at the energy corresponding to adding a constant Doppler shift for each observation (vertical dashed line) to the line’s laboratory energy. For an ideal, symmetric line, the colored and black curves should match. Most of the lines observed have complex line profiles.
thumbnail Fig. 15

Deepdip spectrum from ObsID 8525 from a subset of the strong dip data in the 0.9–2.3 keV (13.8–5.4 Å) range. The data are rebinned to S∕N ≥ 5, revealing emission lines from the photoionized zone around the black hole.

4 Summary

We have analyzed three Chandra-HETGS observations of Cyg X-1 around superior conjunction of the system, where the structure of the stellar wind causes absorption features, so-called dips, in the X-ray light curve. With the help of color–color diagrams, we divide the data into four different dip stages containing a comparable total number of counts. For each dip stage, we fit the silicon and sulfur regions of the spectrum by applying a local power law continuum and Gaussian lines to model the absorption lines for different ionization stages of silicon and sulfur. We find lower charge ions appearing in the deeper dip stages, equivalent widths/absorption columns rising for deeper dip stages, phase dependent Doppler shifts for the lines, asymmetric line shapes, and emission lines emerging in the deepest dip.

The lower charge states of both elements only appear in the deeper dip stages, pointing toward colder material shielded from the irradiation by the black hole. The total columns in the dip stages derived from the equivalent widths of the lines represent lower limits to the real values, as we cannot see the fully ionized material adding to the column, which is likely to carry most of the mass. The derived NH values of NH ~ 2–4 × 1021 cm−2 are in agreement with previous results from Paper I, although there are slight differences due to differences in the continuum.

The Doppler shifts of the lines with respect to measured laboratory values (Hell et al. 2016) show a scatter between the different dipstages, but the same trend within a single observation, pointing toward a single structure containing different ionization stages, which moves as a whole with a certain speed. The Doppler velocities show a clear modulation with orbital phase and are in agreement with the results found in Paper II. However, these velocities match neither the expected wind velocities calculated from a toy model after Gies & Bolton (1986a) nor free fall nor the Kepler velocity around the black hole. Thus, from the line shifts it is not possible to tell where the material causing the absorption dips in the light curve is located.

For the three observations around ϕorb = 0, the lines show asymmetries. This is expected from previous observations near superior conjunction (see Paper I and references therein). Clear P Cygni-like profiles are reported only near ϕorb = 0.5.

Within the deepest dip, there is an emission line spectrum emerging that is outshone otherwise. Together with the Si XIII emission line that does not change its total line flux during the dip stages, this emission line spectrum indicates a much larger area responsible for the emission than what is absorbed by the wind structures, or an area distant from the absorber such that the line of sight toward it is not affected. A possibility is the almost fully ionized photoionization region in the stellar wind surrounding the X-ray emission region.

5 Conclusions

Our analysis gives us a glimpse at the properties of the structures causing the dips. From the light curve, we can constrain the duration of a dip, and thus of a passage of the clump through our line of sight, from 0.5 ks for the shorter dips up to 5 ks for the longer dips. This assumption is only valid if a dip is caused by a single, bigger clump with a cold or dense core instead of many small clumps in the line of sight.

Assuming a distance of 0.25 d between the clump and the black hole (d being the distance between black hole and donor; see Sect. 3.1.2), a simple Keplerian approach (i.e., the clump orbiting the black hole) gives a velocity of ~460 km s−1 perpendicular to the line of sight. From the velocity of the clump and the duration of the dip, we can estimate the diameter of the clumps to (2.3–23) × 108 m, equivalent to (0.33–3.3) R⊙ or (0.02–0.2) R*. As an advancing clump grows in size (Sundqvist et al. 2012), we can calculate its ejection size (Grinberg et al. 2015) to about a third of the size that we measure at the distance of 0.25 d.

With the current data we can only speculate whether this picture of single, big clumps is the correct one. The morphology of the light curve in ObsID 8525 supports this idea, whereas ObsID 3814 with its rapid variability appears to be more in favor of a large number of small clumps in our line of sight. For a sufficiently high number of smaller clumps, those closer to the black hole are partly shielding the outer ones from the ionizing radiation, which is in agreement with our observation of lower ionization stages during strong dipping. A possible explanation would then be a mixture of some larger clumps, causing the distinct dips in the light curve, and many smaller clumps.

This idea is consistent with the picture of structured, radiation pressure driven O-star winds (Castor et al. 1975). As shown, for example, by Owocki & Rybicki (1984), the stationary solution for a line-driven wind is unstable. These instabilities grow quickly and result in strong shocks. Nonstationary hydrodynamic simulations show that dense cool shells of gas already form in deep wind regions close to the photosphere of the O star (Feldmeier et al. 1997; Dessart & Owocki 2003; Oskinova et al. 2012; Sundqvist et al. 2018). Large density, velocity, and temperature variations due to the de-shadowing instability further compress the gas in these shells and fragment these shells into clumps. Theory predicts thus that the wind is actually a two-phase medium consisting of tenuous, hot gas (covering most of the volume) and embedded cool, dense clumps (containing most of the mass). Observational evidence for this clumping has been found for isolated O-type supergiants (e.g., Eversberg et al. 1998; Markova et al. 2005; Bouret et al. 2005; Fullerton et al. 2006; Oskinova et al. 2006), but also for high mass X-ray binaries (e.g., Torrejón et al. 2015).

Speculatively, the picture of the dip behavior would then be that observations at phase 0, which sample the region of the stellar wind closer to the donor star, should show a wide spectrum of clump sizes. Observations at phase 0.75 sample regions that are farther away from the wind, where the clump size distribution has changed, probably favoring larger clumps, resulting in less short-term NH variation. Clearly, however, using only three observations from different orbits and because of the need to perform joint spectroscopy of multiple dips, we do not have sufficient statistics to make a firm statement about the distribution and size of the clumps. Our estimates of (0.02–0.2) R* for the diameter of a clump are in agreement with the high end of the recent 2D simulation results of Sundqvist et al. (2018), however, who have found a typical length scale for clumps of  0.01 R* at a distance of two stellar radii from the star. We cannot distinguish between quasi-spherical clumps and, for instance, pancake-shaped clumps as proposed by Oskinova et al. (2012), which would have different optical depths for different lines of sight. Sundqvist et al. (2018) have even found various different clump shapes coexisting in their 2D simulations. Athena with its larger effective area will provide more insight as it will be able to look at single dips.

Acknowledgements

We thank the referee for their very constructive report which improved the paper. The research leading to these resultswas funded by the Bundesministerium für Wirtschaft und Technologie under grant numbers DLR 50 OR 0701 and DLR 50 OR1113, by LLNL under Contract DE-AC52-07NA27344, and is supported by NASA grants to LLNL. V.G. is supported through the Margarete von Wrangell fellowship by the ESF and the Ministry of Science, Research and the Arts Baden-Württemberg. This research has made use of ISIS functions provided by ECAP/Remeis observatory and MIT (http:// www.sternwarte.uni-erlangen.de/isis/). We thank John E. Davis for providing the slxfig module used for creating the presented plots. We thank Ivica Miškovičová for her initial reduction of the data sets used here.

References

  1. Bałucińska-Church, M., Church, M. J., Charles, P. A., et al. 2000, MNRAS, 311, 861 [NASA ADS] [CrossRef] [Google Scholar]
  2. Blondin, J. M. 1994, ApJ, 435, 756 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Blondin, J. M., & Woo, J. W. 1995, ApJ, 445, 889 [NASA ADS] [CrossRef] [Google Scholar]
  4. Blondin, J. M., Kallman, T. R., Fryxell, B. A., & Taam, R. E. 1990, ApJ, 356, 591 [NASA ADS] [CrossRef] [Google Scholar]
  5. Blondin, J. M., Stevens, I. R., & Kallman, T. R. 1991, ApJ, 371, 684 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bode, M. F., Priedhorsky, W. C., Norwell, G. A., & Evans, A. 1985, ApJ, 299, 845 [NASA ADS] [CrossRef] [Google Scholar]
  7. Boroson, B., & Vrtilek, S. D. 2010, ApJ, 710, 197 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bouret, J., Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bowyer, S., Byram, E. T., Chubb, T. A., & Friedman, H. 1965, Ann. Astrophys., 28, 791 [NASA ADS] [Google Scholar]
  10. Brocksopp, C., Fender, R. P., Larionov, V., et al. 1999, MNRAS, 309, 1063 [NASA ADS] [CrossRef] [Google Scholar]
  11. Caballero-Nieves, S. M., Gies, D. R., Bolton, C. T., et al. 2009, ApJ, 701, 1895 [NASA ADS] [CrossRef] [Google Scholar]
  12. Canizares, C. R., Davis, J. E., Dewey, D., et al. 2005, PASP, 117, 1144 [NASA ADS] [CrossRef] [Google Scholar]
  13. Castor, J. I., Abbott, D. C., Klein, R. I. 1975, ApJ, 195, 157 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dessart, L., & Owocki, S. P. 2003, A&A, 406, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Dotani, T., Inoue, H., Mitsuda, K., et al. 1997, ApJ, 485, L87 [NASA ADS] [CrossRef] [Google Scholar]
  16. Drake, G. W. 1988, Can. J. Phys., 66, 586 [NASA ADS] [CrossRef] [Google Scholar]
  17. El Mellah, I., Sander, A. A. C., Sundqvist, J. O., & Keppens, R. 2019, A&A, 622, A189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Eversberg, T., Lepine, S., & Moffat, A. F. J. 1998, ApJ, 494, 799 [NASA ADS] [CrossRef] [Google Scholar]
  19. Feldmeier, A., Puls, J., Pauldrach, A. W. A. 1997, A&A, 322, 878 [NASA ADS] [Google Scholar]
  20. Feng, Y. X., & Cui, W. 2002, ApJ, 564, 953 [NASA ADS] [CrossRef] [Google Scholar]
  21. Feng, Y. X., Tennant, A. F., & Zhang, S. N. 2003, ApJ, 597, 1017 [NASA ADS] [CrossRef] [Google Scholar]
  22. Friend, D. B., & Castor, J. I. 1982, ApJ, 261, 293 [NASA ADS] [CrossRef] [Google Scholar]
  23. Fullerton, A. W., Massa D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gandhi, P., Rao, A., Johnson, M. A. C., et al. 2019, MNRAS, 485, 2642 [NASA ADS] [CrossRef] [Google Scholar]
  25. Garcia, J. D., & Mack, J. E. 1965, J. Opt. Soc. Am., 55, 654 [NASA ADS] [CrossRef] [Google Scholar]
  26. Gies, D. R., & Bolton, C. T. 1986a, ApJ, 304, 371 [NASA ADS] [CrossRef] [Google Scholar]
  27. Gies, D. R., & Bolton, C. T. 1986b, ApJ, 304, 389 [NASA ADS] [CrossRef] [Google Scholar]
  28. Gies, D. R., Bolton, C. T., Thomson, J. R., et al. 2003, ApJ, 583, 424 [NASA ADS] [CrossRef] [Google Scholar]
  29. Gies, D. R., Bolton, C. T., Blake, R. M., et al. 2008, ApJ, 678, 1237 [NASA ADS] [CrossRef] [Google Scholar]
  30. Grevesse, N., Noels, A., & Sauval, A. J. 1996, in Cosmic Abundances, eds. S. S. Holt & G. Sonneborn (San Francisco: ASP Conf. Ser.), 99, 117 [NASA ADS] [Google Scholar]
  31. Grinberg, V., Hell, N., Pottschmidt, K., et al. 2013, A&A, 554, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Grinberg, V., Leutenegger, M. A., Hell, N., et al. 2015, A&A, 576, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Grinberg, V., Hell, N., El Mellah, I., et al. 2017, A&A, 608, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Hanke, M., Wilms, J., Nowak, M. A., et al. 2008, in VII Microquasar Workshop: Microquasars and Beyond. PoS (MQW7), 29 [Google Scholar]
  35. Hanke, M., Wilms, J., Nowak, M. A., et al. 2009, ApJ, 690, 330 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hell, N., Miškovičová, I., Brown, G. V., et al. 2013, Phys. Scr., T156, 014008 [NASA ADS] [CrossRef] [Google Scholar]
  37. Hell, N., Brown, G. V., Wilms, J., et al. 2016, ApJ, 830, 26 [NASA ADS] [CrossRef] [Google Scholar]
  38. Hemphill, P. B., Rothschild, R. E., Markowitz, A., et al. 2014, ApJ, 792, 14 [NASA ADS] [CrossRef] [Google Scholar]
  39. Herrero, A., Kudritzki, R. P., Gabler, R., et al. 1995, A&A, 297, 556 [NASA ADS] [Google Scholar]
  40. Houck, J. C. 2002, in High Resolution X-ray Spectroscopy with XMM-Newton and Chandra, Mullard Space Science Laboratory, Holmbury St Mary [Google Scholar]
  41. Houck, J. C., & Denicola, L. A. 2000, in Astronomical Data Analysis Software and Systems IX, eds. N. Manset, C. Veillet, & D. Crabtree, ASP Conf. Ser., 216, 591 [NASA ADS] [Google Scholar]
  42. Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kitamoto, S., Miyamoto, S., Tanaka, Y., et al. 1984, PASJ, 36, 731 [NASA ADS] [Google Scholar]
  44. Krolik, J. H., McKee, C. F., & Tarter, C. B. 1981, ApJ, 249, 422 [NASA ADS] [CrossRef] [Google Scholar]
  45. Li, K. F., & Clark, G. W. 1974, ApJ, 191, L27 [NASA ADS] [CrossRef] [Google Scholar]
  46. Manousakis, A., & Walter, R. 2011, A&A, 526, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Manousakis, A., & Walter, R. 2015, A&A, 575, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Markova, N., Puls, J., Scuderi, S., & Markov, H. 2005, A&A, 440, 1133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Martínez-Núñez, S., Kretschmar, P., Bozzo, E., et al. 2017, Space Sci. Rev., 212, 59 [NASA ADS] [CrossRef] [Google Scholar]
  50. Mason, K. O., Hawkins, F. J., Sanford, P. W., et al. 1974, ApJ, 192, L65 [NASA ADS] [CrossRef] [Google Scholar]
  51. Miller, J. M., Fabian, A. C., Wijnands, R., et al. 2002, ApJ, 578, 348 [NASA ADS] [CrossRef] [Google Scholar]
  52. Miškovičová, I., Hanke, M., Wilms, J., et al. 2011, Acta Polytech., 51, 4 [Google Scholar]
  53. Miškovičová, I., Hell, N., Hanke, M., et al. 2016, A&A, 590, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Morton, D. C. 1967, ApJ, 150, 535 [NASA ADS] [CrossRef] [Google Scholar]
  55. Muijres, L. E., Vink, J. S., de Koter, A., et al. 2012, A&A, 537, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Murdin, P., & Webster, B. L. 1971, Nature, 233, 110 [NASA ADS] [CrossRef] [Google Scholar]
  57. Noble, M. S., & Nowak, M. A. 2008, PASP, 120, 821 [NASA ADS] [CrossRef] [Google Scholar]
  58. Nowak, M. A., Vaughan, B. A., Wilms, J., et al. 1999, ApJ, 510, 874 [NASA ADS] [CrossRef] [Google Scholar]
  59. Nowak, M. A., Hanke, M., Trowbridge, S. N., et al. 2011, ApJ, 728, 13 [NASA ADS] [CrossRef] [Google Scholar]
  60. Orosz, J. A., McClintock, J. E., Aufdenberg, J. P., et al. 2011, ApJ, 742, 84 [NASA ADS] [CrossRef] [Google Scholar]
  61. Oskinova, L. M., Feldmeier, A., & Hamann, W. 2006, MNRAS, 372, 313 [NASA ADS] [CrossRef] [Google Scholar]
  62. Oskinova, L. M., Feldmeier, A., & Kretschmar, P. 2012, MNRAS, 421, 2820 [NASA ADS] [CrossRef] [Google Scholar]
  63. Owocki, S. P., & Rybicki, G. B. 1984, ApJ, 284, 337 [NASA ADS] [CrossRef] [Google Scholar]
  64. Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 [NASA ADS] [CrossRef] [Google Scholar]
  65. Palmeri, P., Quinet, P., Mendoza, C., et al. 2008, ApJS, 177, 408 [NASA ADS] [CrossRef] [Google Scholar]
  66. Parker, M. L., Tomsick, J. A., Miller, J. M., et al. 2015, ApJ, 808, 9 [NASA ADS] [CrossRef] [Google Scholar]
  67. Parsignault, D. R., Epstein, A., Grindlay, J., et al. 1976, Ap&SS, 42, 175 [NASA ADS] [CrossRef] [Google Scholar]
  68. Pepe, C., Vila, G. S., & Romero, G. E. 2015, A&A, 584, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Porquet, D., Dubau, J., & Grosso, N. 2010, Space Sci. Rev., 157, 103 [NASA ADS] [CrossRef] [Google Scholar]
  70. Poutanen, J., Zdziarski, A. A., Ibragimov, A. 2008, MNRAS, 389, 1427 [NASA ADS] [CrossRef] [Google Scholar]
  71. Pravdo, S. H., White, N. E., Becker, R. H., et al. 1980, ApJ, 237, L71 [NASA ADS] [CrossRef] [Google Scholar]
  72. Puls, J., Markova, N., Scuderi, S., et al. 2006, A&A, 454, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Puls, J., Vink,J. S., & Najarro, F. 2008, A&ARv, 16, 209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Rahoui, F., Lee, J. C., Heinz, S., et al. 2011, ApJ, 736, 63 [NASA ADS] [CrossRef] [Google Scholar]
  75. Reid, M. J., McClintock, J. E., Narayan, R., et al. 2011, ApJ, 742, 83 [NASA ADS] [CrossRef] [Google Scholar]
  76. Remillard, R. A., & Canizares, C. R. 1984, ApJ, 278, 761 [NASA ADS] [CrossRef] [Google Scholar]
  77. Sako, M., Liedahl, D. A., Kahn, S. M., & Paerels, F. 1999, ApJ, 525, 921 [NASA ADS] [CrossRef] [Google Scholar]
  78. Schulz, N. S., Cui, W., Canizares, C. R., et al. 2002a, ApJ, 565, 1141 [NASA ADS] [CrossRef] [Google Scholar]
  79. Schulz, N. S., Canizares, C. R., Lee, J. C., & Sako, M. 2002b, ApJ, 564, L21 [NASA ADS] [CrossRef] [Google Scholar]
  80. Sundqvist, J. O., & Owocki, S. P. 2013, MNRAS, 428, 1837 [NASA ADS] [CrossRef] [Google Scholar]
  81. Sundqvist, J. O., Owocki, S. P., Cohen, D. H., et al. 2012, MNRAS, 420, 1553 [NASA ADS] [CrossRef] [Google Scholar]
  82. Sundqvist, J. O., Owocki, S. P., Puls, J. 2018, A&A, 611, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Tomsick, J. A., Nowak, M. A., Parker, M., et al. 2014, ApJ, 780, 78 [NASA ADS] [CrossRef] [Google Scholar]
  84. Torrejón, J. M., Schulz, N. S., Nowak, M. A., et al. 2015, ApJ, 810, 102 [NASA ADS] [CrossRef] [Google Scholar]
  85. Verner, D. A., Verner, E. M., & Ferland, G. J. 1996, At. Data Nucl. Data Tables, 64, 1 [NASA ADS] [CrossRef] [Google Scholar]
  86. Vrtilek, S. D., Boroson, B. S., Hunacek, A., et al. 2008, ApJ, 678, 1248 [NASA ADS] [CrossRef] [Google Scholar]
  87. Walborn, N. R. 1973, ApJ, 179, p.L123 [NASA ADS] [CrossRef] [Google Scholar]
  88. Watanabe, S., Sako, M., Ishida, M., et al. 2006, ApJ, 651, 421 [NASA ADS] [CrossRef] [Google Scholar]
  89. Webster, B. L., & Murdin, P. 1972, Nature, 235, 37 [NASA ADS] [CrossRef] [Google Scholar]
  90. Wen, L., Cui, W., Levine, A. M., & Bradt, H. V. 1999, ApJ, 525, 968 [NASA ADS] [CrossRef] [Google Scholar]
  91. Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 [NASA ADS] [CrossRef] [Google Scholar]
  92. Wilms, J., Nowak, M. A., Pottschmidt, K., et al. 2006, A&A, 447, 245 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Wojdowski, P. S., Liedahl, D. A., Sako, M., et al. 2003, ApJ, 582, 959 [NASA ADS] [CrossRef] [Google Scholar]
  94. Xiang, J., Lee, J. C., Nowak, M. A., & Wilms, J. 2011, ApJ, 738, 78 [NASA ADS] [CrossRef] [Google Scholar]
  95. Zhang, S. N., Cui, W., Harmon, B. A., et al. 1997, ApJ, 477, L95 [NASA ADS] [CrossRef] [Google Scholar]
1

New Gaia measurements appear to favor a slightly larger distance of 2.38−0.17+0.20$2.38^{+0.20}_{-0.17}$ kpc, possibly due to systematic error caused by the high optical brightness of the system and/or an orbital wobble (Gandhi et al. 2019).

2

The simpler non-dip selection in Paper II was done using only the ratio A∕C.

3

See “The Chandra ABC Guide to Pileup”, CXC, 2010, http://cxc.harvard.edu/ciao/download/doc/pileup_abc.pdf

4

See http://www.atomdb.org; line energies without a reference refer to AtomDB.

5

Lines denoted as part of the Lyman series typically have two unresolved components with the configuration of 1s1∕2 2S1∕2–np3∕2 2P3∕2 and 1s1∕2 2S1∕2–np1∕2 2P1∕2, respectively. The labels α, β, γ etc. represent the principal quantum number n = 2, 3, 4, …. The line energy given in the text corresponds to the mean energy of the components, weighted with the respective statistical weight (2:1).

6

There was a sign error in the code that produced the figure in Paper II. This led to a flip of sign for the velocity curves and a slight phase shift. These changes do not have an influence on the interpretation of the figure in Paper II. Figure 11 shows the corrected velocity curves.

All Tables

Table 1

Chandra-HETGS observations used in this paper.

In the text
Table 2

Exposure time, texp, and dip stage averaged count rates for the total observations and silicon (1.605–2.045 keV ≡ 7.725–6.063 Å) and sulfur (2.295–2.7 keV ≡ 5.402–4.6 Å) bands for all observations.

In the text
Table 3

Detector constants and continuum parameters for all observations.

In the text
Table 4

Observed silicon and sulfur line centers for all observations.

In the text
Table 5

Evolution of the equivalent widths of the silicon and sulfur lines during the dipping stages for all observations.

In the text
Table 6

Column densities as calculated from the silicon line equivalent widths.

In the text
Table 7

Column densities as calculated from the sulfur line equivalent widths.

In the text

All Figures

thumbnail Fig. 1

Top panel: light curves with a 25.5 s time resolution of all three observations as a function of orbital phase. Bottom panel: corresponding hardness ratios. All observations were in the low/hard state. Dips are strongest at ϕorb ~ 0.0, become weaker at ϕorb ~ 0.2 and ϕorb ~ 0.75, and completely fade at ϕorb ~ 0.5 (see also Paper II). Colors indicate the selected dipping stages (see right panel of Fig. 3 for the color denotation and Sect. 2.3 for details on the data selection).
In the text
thumbnail Fig. 2

Left panel: effect of partial absorption with increasing NH at constant covering fraction with NH(0) = 5.4 × 1021 cm−2$N_{\mathrm{H}}^{(0)}\,{=}\,5.4\,{\times}\,10^{21}\,\mathrm{cm}^{-2}$ after Paper I. Right panel: color–color diagram showing the tracks for different covering fractions and column densities of the second absorber. See text for details.
In the text
thumbnail Fig. 3

Color–color diagrams for all three observations, showing hardness ratios calculated for soft, intermediate, and hard X-ray bands. As discussed in the text, the observed tracks are determined by absorption: following the track from the top right corner,representing the persistent flux, to the bottom left corner, showing dips, absorption increases and therefore both colors become harder. During the deepest dips, the softer color becomes softer again, while the hard color does not change. Thisbehavior can be explained by partial covering (Sect. 2.2). The spectrum hardens considerably during thesedips. The colors denote the four different dip stages (see right panel). The black line indicates our polynomial fit to the track; see Sect. 2.3 for details.
In the text
thumbnail Fig. 4

Rebinned HEG spectra of the different dipping stages of ObsID 8525. For clarity, the non-dip spectrum (orange, S/N ≥ 13) is scaled by afactor of 8, the weak dip spectrum (red, S/N ≥ 15) by a factor of 4, and the dip spectrum (violet, S/N ≥ 15) by a factor of 2. The strong dip spectrum (blue, S/N ≥ 17) is not scaled. The box denotes the region containing the Si and S lines. The color scheme is the same as in Figs. 1 and 3.
In the text
thumbnail Fig. 5

Evolution of silicon (left column) and sulfur (right column) lines with dipping stages for the three observations considered (from top to bottom). The panels show combined HEG and MEG ± 1st spectra relative to their respective continua. The solid black line shows the line fits.
In the text
thumbnail Fig. 6

Column densities for each dip stage derived from the respective equivalent width for all Si and S ion absorption lines detected in ObsIDs 3814 (yellow diamonds), 8525 (green triangles), and 9847 (blue squares). Negative values of the equivalent width correspond to absorption lines. The column densities of lower charged ions increase with increasing dip depth, whereas column densities of highly charged ions decrease with dip level. Be-like Si XI and Li-like Si XII lines are contaminated by the Si He z emission line.
In the text
thumbnail Fig. 7

Column densities vs. ionization potential for the respective silicon and sulfur ions in all dip stages.
In the text
thumbnail Fig. 8

Doppler velocities for the various ionization states of Si VIII–Si XIV (empty) and S X–S XVI (filled) in all dipping stages (non-dip: orange circles; weak dip: red triangles; dip: violet squares; and strong dip: blue stars) over all observations in the different orbital phases.
In the text
thumbnail Fig. 9

Distribution of velocity shifts of all Si and S lines detected in the three ObsIDs with their medians. The peak velocity clearly shifts between the observations.
In the text
thumbnail Fig. 10

Sketch of a clump passing our line of sight and causing the different dip stages. Red color indicates more highly ionized regions; blue color regions indicate lower ionization balance.
In the text
thumbnail Fig. 11

Radial velocities of the black hole (red dashed line) and the companion (orange dashed line), and wind velocities projected onto the line of sight as a function of distance from the black hole. Distances are given in units of the distance between the black hole and its donor, d. Colored boxes indicate the range of measured non-dip Doppler shifts from the Chandra observations (red: ObsID 3815; yellow: ObsID 3814; green: ObsID 8525; and blue: ObsID 9847, modified and corrected from Paper II, Fig. 11). Symbols show the median (star), mean (circle), first (point-up triangle), and third quartile (point-down triangle) of the distribution of Doppler shifts for the Si and S lines.
In the text
thumbnail Fig. 12

Incident SED used for the XSTAR simulations. The dotted and dashed lines show the SED with and without the optical companion star, respectively. We note that the SED is always renormalized to the assumed model luminosity.
In the text
thumbnail Fig. 13

Left panel: predicted columns for the ions of Si and S for a spherical cloud with NH = 1022 cm−2 that is irradiated by a X-ray source with the spectrum defined by Fig. 12. Right panel: lines of constant ion column density color-coded for each different ion of Si (panels a and b) and S (panels c and d), for XSTAR simulated values of equivalent hydrogen column density, NH, and ionization parameter, logξ. The assumed ion columns are the values measured during the strong dip of observation 8525 (Table 5). If a simple XSTAR model described the observation, we would expect these lines to cross in one point (within their error bars). The left column of the figure shows calculations without the effect of the optical companion (OC); the right column includes the UV radiation from the companion.
In the text
thumbnail Fig. 14

Original line profiles (colored) for the non-dip (top, orange), weak dip (second from top, red), dip (third from top, purple), and strong dip (bottom, blue) lines for selected ions, as well as the same data flipped (black, see also text for further information) at the energy corresponding to adding a constant Doppler shift for each observation (vertical dashed line) to the line’s laboratory energy. For an ideal, symmetric line, the colored and black curves should match. Most of the lines observed have complex line profiles.
In the text
thumbnail Fig. 15

Deepdip spectrum from ObsID 8525 from a subset of the strong dip data in the 0.9–2.3 keV (13.8–5.4 Å) range. The data are rebinned to S∕N ≥ 5, revealing emission lines from the photoionized zone around the black hole.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.