X-RAYING THE INTERGALACTIC O vi ABSORBERS

, , , , , , , and

Published 2009 May 14 © 2009. The American Astronomical Society. All rights reserved.
, , Citation Y. Yao et al 2009 ApJ 697 1784 DOI 10.1088/0004-637X/697/2/1784

0004-637X/697/2/1784

ABSTRACT

The observed intergalactic O vi absorbers at z>0 have been regarded as a significant reservoir of the "missing baryons." However, to fully understand how these absorbers contribute to the baryon inventory, it is crucial to determine whether the systems are collisionally ionized or photoionized (or both). Using the identified intergalactic O vi absorbers as tracers, we search for the corresponding X-ray absorption lines, which are useful for finding the missing baryons and for revealing the nature of the O vi absorbers. Stacking the Chandra grating spectra along six AGN sight lines, we obtain three spectra with signal-to-noise ratios of 32, 28, and 10 per 12.5 mÅ spectral bin around the expected O vii Kα wavelength. These spectra correspond to O vi absorbers with various dynamic properties. We find no detectable Ne ix, O vii, O viii, N vii, or C vi absorption lines in the spectra, but the high counting statistics allows us to obtain firm upper limits on the corresponding ionic column densities (in particular $N_{{\rm O}{\mathsc{vii}}}\lesssim 10 N_{{\rm O}{\mathsc{vi}}}$ on average at the 95% confidence level). Jointly analyzing these nondetected X-ray lines with the averaged O vi column density, we further limit the average temperature of the O vi-bearing gas to be T ≲ 105.7 K in collisional ionization equilibrium. We discuss the implications of these results for physical properties of the putative warm–hot intergalactic medium and its detection in future X-ray observations.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Hydrodynamic simulations of cosmic large-scale structure indicate that, in the local universe (z ≲ 1), the shock-heated gas in the tenuous warm–hot intergalactic medium (WHIM) at temperatures T ∼ 105−7 K contains a substantial amount of baryonic matter (Cen & Ostriker 1999; Davé et al. 2001), providing a possible solution for the so-called "missing baryons" problem (e.g., Persic & Salucci 1992; Fukugita et al. 1998). In a collisionally ionized gas at these temperatures, the most abundant heavy elements (e.g., C, N, O, and Ne) are in their high-ionization (e.g., Li-, He-, and H-like) states (Sutherland & Dopita 1993), whose K- and L-shell transitions are in the X-ray and far-ultraviolet (far-UV) wavelength bands, respectively. The existence of the WHIM has been suggested through detections of X-ray emission from high-density regions near galaxies, groups, and clusters (e.g., Wang et al. 1997, 2004; Finoguenov et al. 2003; Sołtan 2007; Mannucci et al. 2007; see also Durret et al. 2008 for a review), but these regions are expected to contain only a small portion of the WHIM gas. The majority is located in low-density cosmic-web filaments (Davé et al. 2001), whose emissions are hard to detect with current X-ray and far-UV telescopes. These filaments are most easily probed through absorption lines imprinted on spectra of background sources.

With the high-sensitivity spectrographs aboard the Hubble Space Telescope (HST) and the Far Ultraviolet Spectroscopic Explorer (FUSE), the intergalactic O vi absorption lines (at rest-frame wavelengths 1031.93 Å and 1037.62 Å) have routinely been detected in many background active galactic nucleus (AGN) spectra (e.g., Savage et al. 1998; Tripp et al. 2000; Tripp & Savage 2000; Shull et al. 2003; Prochaska et al. 2004; Danforth & Shull 2005; Cooksey et al. 2008). However, the nature of these O vi absorbers is still under debate. In a survey of the intergalactic medium (IGM) absorption lines, Danforth & Shull (2008) found good correlations between O vi and N v in their column densities and power-law slopes of the column density distributions, which are distinct from those of H i, C iii, and Si iii. These features, along with the theoretical expectations of high post-shock temperatures and long cooling times of the WHIM, motivated them to argue for a multiphase nature of the IGM. In this picture, O vi and N v trace the canonical WHIM at temperatures of 105–106 K, while the low-ionization ions such as C ii–C iv and Si ii–Si iv are predominantly photoionized. But in two other independent surveys that include a number of common sight lines to those used by Danforth & Shull (2008), Tripp et al. (2008), and Thom & Chen (2008) found that 30%–40% of the O vi absorbers have velocity centroids that are well aligned with those of the associated H i absorbers and that some absorber temperatures inferred from the line widths, T < 105 K, are well below that expected in the canonical WHIM. Tripp et al. (2008) also showed that these velocity-aligned absorbers can be naturally explained by photoionization models. These latter findings are consistent with previous investigations along individual sight lines based on the kinematic and chemical properties of the high- and low-ionization absorbers (e.g., Tripp & Savage 2000; Savage et al. 2002; Sembach et al. 2004; Prochaska et al. 2004; Lehner et al. 2006). However, the UV data usually do not have high enough quality to uniquely constrain the ionization mechanism; Tripp et al. (2008) also showed that the O vi lines could alternatively arise in hot interface layers on the surface of low-ionization clouds if the O vi systems are always multiphase media. Current data cannot rule out this possibility, and many O vi absorbers are fully consistent with this hypothesis (see Section 4.2 in Tripp et al. 2008). So far, the strongest indication of the collisional ionization origin of the O vi absorbers is the detection of Ne viii absorption lines (at rest-frame wavelengths 770.4 Å and 780.3 Å) toward HE 0226-4110 at za = 0.2070 and toward 3C 263 at za = 0.3257 (Savage et al. 2005; Narayanan et al. 2008). However, a systematic search toward other sight lines has failed to find similar systems (Lehner et al. 2006).

The claimed detections of the WHIM in the X-ray band are still controversial. In the spectrum of PKS 2155–304 observed with the Chandra Advanced CCD Imaging Spectrometer (ACIS), Fang et al. (2002) reported an O viii Kα absorption line at z = 0.0553, claiming the first detection of the WHIM in the X-ray. This detection has not been confirmed, although it cannot be ruled out either, by the Chandra High-Resolution Camera (HRC) observations and the subsequent XMM-Newton Reflection Grating Spectrometer (RGS) observations (Cagnoni et al. 2004; Williams et al. 2007; Fang et al. 2007). Nicastro et al. (2005) reported detections of two z>0 O vii WHIM systems in the Chandra spectra of Mrk 421, which were not confirmed from the XMM-Newton RGS observations of the same source (Kaastra et al. 2006; Rasmussen et al. 2007). Kaastra et al. (2006) also questioned the statistical significance of the Mrk 421 detections in the Chandra spectra. The other claimed z>0 X-ray WHIM absorptions are mostly detected at marginal significance (e.g., Mathur et al. 2003). Very recently, Buote et al. (2009) have reported a 3σ detection of O vii Kα absorption affiliated with a large-scale structure (the Sculptor Wall). The O vii and O viii absorption lines at z ≃ 0, which may be partially due to the WHIM in the Local Group (e.g., Nicastro et al. 2002; Fang et al. 2003), have been relatively well detected, but these absorptions are severely confused by the contributions from the Galactic hot gas close to (or within) the Milky Way disk (Sembach et al. 2003; Yao & Wang 2005, 2007; Wang et al. 2005; Fang et al. 2006; Bregman & Lloyd-Davies 2007; Yao et al. 2008).

With the controversial X-ray detections of the WHIM and the debatable nature of the observed far-UV absorbers, two key questions still remain open: (1) Does the WHIM exist? (2) If it does, what are its physical properties?

The highly ionized X-ray absorption lines are believed to be useful for finding the WHIM gas and for probing the nature of the O vi absorbers. The O vii ion can trace gas over a broad temperature range, and its column density is expected to be ≳10 times higher than that of O vi in a shock-heated gas with shock temperatures T ≳ 105.7 K (for overdensity δ∼10–100; e.g., Furlanetto et al. 2005). Moreover, cosmological simulations predict that hot O vi absorbers should have comparably strong (or even stronger) affiliated O vii absorption lines (see, e.g., Figures 8–10 in Cen & Fang 2006 and Figure 14 in Chen et al. 2003). However, because of the limited sensitivity and spectral resolution of the current X-ray observatories like Chandra and XMM-Newton, searching for X-ray lines in the WHIM is currently a difficult task. For instance, the column densities of the X-ray absorbing ions, even for the most promising one (O vii), are expected to be small with NOvii ≈ 1015 cm−2 (e.g., Chen et al. 2003). For the most part, the background AGNs are also relatively faint in the X-ray band. Consequently, very long exposures (see Section 5) are required to collect enough photons to conduct a blind search for X-ray absorption lines from the WHIM along a random line of sight. These difficulties contribute to the controversies and frustrations surrounding the reported X-ray detections.

In this work, we search for the X-ray absorption lines of O vii, as well as O viii, Ne ix, N vii, and C vi, by extensively exploring the archived Chandra grating observations. We use the identified O vi absorbers as tracers to avoid uncertainties of a blind search. We also develop a technique to stack all the applicable X-ray observations and obtain spectra with unprecedented counting statistics.

The paper is organized as follows. In Section 2, we list the identified O vi absorbers and the corresponding Chandra observations utilized in this work, and describe the data reduction process. In Section 3, we present our search method and stack the X-ray observations in the rest frame of the observed O vi absorbers. We search for the X-ray absorption lines in the stacked spectra, present the results of the data analysis in Section 4, and discuss the implications of our results in Section 5.

2. O vi ABSORBERS, CHANDRA OBSERVATIONS, AND DATA REDUCTION

There are ∼80 identified IGM O vi systems at redshifts of 0 < zabs < 0.5 toward ∼30 AGN sight lines (Savage et al. 2002; Sembach et al. 2004; Danforth et al. 2006; Danforth & Shull 2005, 2008; Tripp et al. 2008; Thom & Chen 2008). Chandra has observed 12 of these AGNs with high spectral resolution grating instruments. To avoid the potential confusion caused by the absorbers intrinsic to the AGN in identifying the IGM O vi systems and in searching for X-ray lines, we do not use the five AGNs (Mrk 279, NGC 5548, Mrk 509, Ark 564, and NGC 7469) toward which warm absorbers have been detected in the X-ray band (e.g., Scott et al. 2004, 2005; Kaastra et al. 2002; Yaqoob et al. 2003; Matsumoto et al. 2004). Chandra observations of the five AGNs could contribute an additional less than 10% in total to the spectral counts of our final spectra (see below), so excluding them in the data analysis does not significantly affect our results. To avoid the confusion caused by possible absorber misidentification or noise features, we only used the O vi systems that were detected at greater than 3σ significance levels. We further excluded the proximate absorbers that are within ∼2,000 km s-1 to the AGN redshifts; Tripp et al. (2008) have shown that the majority of the proximate O vi systems are found within this velocity interval at low redshifts. The proximate systems can be high-velocity AGN ejecta and/or can be located near the central engine of the AGN. Consequently, these systems should be excluded when searching for WHIM absorbers.

Table 1 summarizes our O vi absorber sample, including the redshifts of the O vi systems (zOvi), λ1031.93 equivalent widths, dispersion velocities (b), and column densities (NOvi) along the six sight lines used in this work. For those absorbers at very close redshifts along a single sight line that cannot be distinguished with the Chandra resolution, only the systems with higher EWs were used. Galaxy surveys at low-z show that most of these O vi absorbers are usually found within several hundred kpc of intervening galaxies along the sight lines (Tripp & Savage 2000; Savage et al. 2002; Shull et al. 2003; Sembach et al. 2004; Tumlinson et al. 2005; Prochaska et al. 2006; Stocke et al. 2006; Tripp et al. 2006; Cooksey et al. 2008). The investigation of the connection between galaxy environment and X-ray absorbers is underway and will be presented in a separate paper. In this work, we focus on searching for corresponding X-ray absorption of the O vi-bearing gas.

Table 1. Properties of the Identified O vi Absorbers

Source Name $z_{\rm{AGN}}$ zOvi EW (mÅ) b (km s−1) $\log [N_{{\rm O}\mathsc{vi}}(\rm cm^{-2})]$ Classb
H1821+643 0.2970 0.0244 27 ± 8 21+9−6 13.44 ± 0.1 Sm
    0.1214 97 ± 14 76+13−11 13.97 ± 0.06 Cm, St
    0.2133 39 ± 9 28 ± 5 13.54 ± 0.06 Cm
    0.2250a 190 ± 10 45 ± 2 14.27 ± 0.02 Cm, St
    0.2264a 32 ± 4 16 ± 2 13.51 ± 0.04 Cm
    0.2453 51 ± 7 26 ± 2 13.71 ± 0.03 Sm
    0.2666 45 ± 7 25 ± 3 13.63 ± 0.04 Sm
3C 273 0.1583 0.0033 31 ± 7 51+11−9 13.44 ± 0.07 Cm
    0.0902 16 ± 3 22 ± 6 13.18 ± 0.06 Cm
    0.1200 24 ± 3 8 ± 3 13.37 ± 0.04 Sm
PG 1116+215 0.1765 0.0593 63 ± 9 36+11−8 13.52 ± 0.08 Sm
    0.1385 83 ± 16 36 ± 7 13.97 ± 0.06 Cm, St
    0.1655 111 ± 9 32+21−13 14.08 ± 0.04 Cm, St
PKS 2155-304 0.1165 0.0540a 32 ± 5 14 ± 6 13.63 ± 0.12 Cm
    0.0572a 44 ± 11 24 ± 7 13.57 ± 0.09 Cm
Ton S180 0.0620 0.0456 62 ± 14 20 13.74 ± 0.06 Sm
PG 1211+143 0.0809 0.0511 187 53+10−9 14.21 ± 0.08 Cm, St
    0.0645a 144 ± 32 54 ± 11 14.16 ± 0.08 Cm, St
    0.0649a 55 ± 10 21 13.83 ± 0.11 Cm
Averagedc     45.4 28 13.62 All
Averagedd     49.0 32 13.65 Cm
Averagede     131.7 59 14.11 St

Notes. Errors are listed in 1σ range. The values of the first three sources (H1821+643, 3C 273, and PG1116+215) are adopted from Tripp et al. (2008), and the remaining values are adopted from Danforth & Shull (2008) except for those at the redshift 0.0511 along the PG 1211+143 sight line, which are adopted from Tumlinson et al. (2005). aToward an individual sight line, these absorbers cannot be distinguished with Chandra resolution, and the system with higher equivalent width of O vi is used. bClassification of the O vi absorbers; "Cm," "Sm," and "St" denote complex, simple, and strong O vi systems, respectively. c,d,eValues are weight-averaged with respect to the spectral counts at the O vii Kα wavelength (21.602 Å) in the shifted spectra for all, complex, and strong O vi absorbers, respectively (Section 4).

Download table as:  ASCIITypeset image

We follow the criteria used in Tripp et al. (2008) to classify the O vi absorbers along the six AGN sight lines into simple and complex systems. In simple systems, the velocity centroids of H i and O vi are well aligned (within 2σ of their velocity uncertainties). In complex absorbers, these velocity centroids are well separated (greater than 2σ apart) and/or there are multiple low- and high-ionization stages (e.g., C ii-C iv, Si ii-Si iii, or N v, in addition to O vi) that indicate the presence of multiple ionization phases. We also classify the O vi absorbers with equivalent width (EW) >80 mÅ as strong systems (e.g., Fang & Bryan 2001). Tripp et al. (2008) classified the O vi absorbers along three of the six sight lines. Using these criteria, we classified all the O vi absorbers reported by Danforth et al. (2006) and Danforth & Shull (2008) toward the other three sight lines. Figure 1 presents examples of each classified system.

Figure 1.

Figure 1. From left to right, examples of the simple, complex, and strong O vi absorbers and the corresponding Lyα systems, along the 3C 273, H1821+643, and PG 1211+143 sight lines (Tripp et al. 2008; Danforth & Shull 2008), respectively. Absorbers that are not associated with the labeled systems are marked with "X."

Standard image High-resolution image

As of 2008 May 6, there were ∼100 archived Chandra grating observations of the six selected targets. In this work, we did not use nonstandard observations with uncertain calibrations. We also excluded several short observations (with an individual exposure ≲10 ks) of PKS 2155-304 that contribute less than 5% to the total spectral counts of the source. Table 2 summarizes the number of observations and the total exposure used in this work.

Table 2. Chandra Grating Observations of the Selected Targets

Src. Name $z_{\rm{AGN}}$ No. of Obs. Exp. (ks)
H1821+643 0.297 5 600
3C 273 0.158 16 360
PG 1116+215 0.176 1 89
PKS 2155-304 0.117 36 760
Ton S180 0.062 1 80
PG 1211+143 0.081 3 141
Sub Total: 62 2030

Download table as:  ASCIITypeset image

For ACIS grating observations, we followed the procedures described in Yao & Wang (2005, 2007) to calibrate the data, extract spectra, and calculate corresponding response functions (RSPs). The same energy grid was applied throughout these procedures. We utilized the medium energy grating (MEG) spectra of the High-Energy Transmission Grating (HETG; Canizares et al. 2005) observations and only used the first grating order spectra of all observations.

For each HRC grating observation, we followed the steps presented in Wang et al. (2005) to obtain the first-order RSP (FRSP) and order-overlapped RSP (ORSP) and spectra. We then extracted the first-order spectra by subtracting the difference of the best-fit-model predicted channel counts between the ORSP and the FRSP from the order-overlapped spectra.

Spectra from positive and negative grating arms of each observation were then coadded, and multiple observations toward a single sight line were further stacked to enhance the counting statistics. Visual inspection revealed no significant and consistent Kα lines of O vii, O viii, Ne ix, N vii, or C vi at the corresponding redshift $z_{{\rm{O}{\mathsc{vi}}}}$ along individual sight lines. Figure 2 shows parts of the final spectra of the six AGNs around the expected IGM O vii Kα lines.

Figure 2.

Figure 2. Chandra spectra of the six sources around the O vii Kα line and the best-fit continuua. The vertical dotted lines mark the Galactic O vii, O i, and O ii Kα lines at 21.602, 23.508, and 23.348 Å, respectively. The thick dashed lines mark the expected positions of the intergalactic O vii Kα lines at the corresponding $z_{{\rm{O}{\mathsc{vi}}}}$ (Table 1). The bin size is 12.5 mÅ.

Standard image High-resolution image

Spectra along the PKS 2155–304 sight line deserve to be mentioned particularly because of their excellent spectral quality (Figure 2). Toward this sight line, Shull et al. (2003) detected two O vi absorption systems at redshifts z = 0.054 and 0.057, associated with a cluster of seven H i absorbers and a group of galaxies (Shull et al. 1998). Figure 3 shows the detected O vi and other nondetected metal absorption lines along with the cluster of H i lines. On the Chandra ACIS-grating spectrum, the claimed z = 0.0553 O viii line at λobs = 20.02 Å is clearly visible with an EW of 5.9(3.1, 8.6) mÅ (90% confidence range), which is consistent with that reported by Fang et al. (2002, 2007). There is no other X-ray line at this redshift. The discrepancy between the HRC-grating, the XMM-Newton RGS (with an accumulated exposure time of 980 ks), and the ACIS-grating spectra still exists; the former two spectra yield upper limits of EW < 3.2 mÅ and <3.7 mÅ, respectively (see also Section 1 and references therein). The line in the coadded HRC+ACIS spectrum is severely diluted with EW < 4.2 mÅ.

Figure 3.

Figure 3. Metal absorption lines associated with a cluster of H i Lyα lines detected in the PKS 2155-304 spectrum, plotted vs. velocity with v = 0 km s-1 at zabs = 0.0542. Both lines of the O vi doublet are clearly detected at zabs = 0.0542 and are marked. Only the O vi λ1031.93 transition is detected at zabs = 0.0572, but the weaker O vi λ1037.62 line is blended with a strong Galactic Fe ii transition.

Standard image High-resolution image

3. AN EFFECTIVE STACKING AND SEARCHING STRATEGY

To facilitate an effective search for X-ray IGM absorption lines, we first blueshifted the X-ray spectral data by the observed O vi absorber redshifts $z_{{\rm{O}{\mathsc{vi}}}}$ and reconstructed the corresponding RSPs accordingly. The data are originally in forms of detected counts distributed in a wavelength grid that was chosen to be the same for all the X-ray spectra. The rest-frame wavelength λ0 is related to the observed wavelength λobs by

Equation (1)

For each absorber, we first use this relationship to recalculate the grid boundaries of the corresponding X-ray spectrum and then cast the spectrum and the RSP onto a rest-frame spectral grid, which is chosen to be the original one. The recalculated spectrum, both shifted and compressed, generally no longer matches this chosen rest-frame grid. Therefore, some of the recalculated spectral bins fall completely within individual rest-frame grid intervals, whereas others enclose a grid boundary. In this latter case, we split the counts in such a bin into two parts and assign them to the two grid intervals adjacent to the enclosed boundary; the weights used in the splitting are proportional to the wavelength overlaps of the bin with the two intervals. In the same way, we recalculate the wavelength boundaries and the distribution probability in the corresponding RSP file. We repeat this procedure for all the absorbers to obtain their respective rest-frame spectral and RSP files; for individual sight lines with n intervening O vi systems (Table 1), we obtained n shifted spectra and RSPs.

To enhance the counting statistics, we coadded these shifted spectra to form a single stacked spectrum and RSP. We constructed three different stacked spectra based on the various O vi absorber classifications (Section 2). Figure 4 illustrates the stacked spectrum that contains all of the shifted spectra corresponding to the 16 O vi systems used in this work (Table 1) regardless of their classification. Figure 5 presents the stacked spectrum for only the complex O vi systems, and Figure 6 is the stacked spectrum for only the strong O vi systems (Table 1). To avoid any potential bias in our results caused by the large contribution of the PKS 2155-304 sight line in the stacked spectra (Figure 2) and by its peculiarity in detecting only the O viii line (Section 1), we also obtained two similar stacked spectra and RSPs like those presented in Figures 4 and 5, but without the PKS 2155-304 contribution. For the same reasons, we also report two sets of results in the following sections when applicable, with and without the contribution of PKS 2155-304. It is interesting to point out that, because of the poor spectral resolution of the X-ray instruments, absorbers at very close redshifts/velocities along a single sight line that cannot be distinguished with the Chandra LETG resolution (∼750 km s-1) are automatically stacked/merged as one velocity component.

Figure 4.

Figure 4. Four ranges of the final stacked spectrum that includes shifted spectra corresponding to all 16 O vi systems (Table 1). The thick lines mark the best-fit continuum, the red lines and text mark the blue-shifted Galactic absorptions (mainly contributed from PKS 2155–304 and 3C 273 sight lines), and blue vertical lines mark the positions of the expected intergalactic neon, oxygen, nitrogen, and carbon absorption lines at the corresponding rest-frame wavelengths. The bin-size is 12.5 mÅ.

Standard image High-resolution image
Figure 5.

Figure 5. Same as Figure 4, except that stacked spectrum contains only spectra corresponding to the complex O vi systems (Table 1).

Standard image High-resolution image
Figure 6.

Figure 6. Same as Figure 4, except that stacked spectrum contains only spectra corresponding to the strong O vi systems (Table 1).

Standard image High-resolution image

Because each spectrum has been shifted by the corresponding $z_{{\rm{O}{\mathsc{vi}}}}$ before being stacked, the X-ray IGM Kα absorption lines of Ne ix, O viii, O vii, N vii, and C vi, if they are associated with the O vi-bearing gas, are expected to be at rest-frame wavelengths of 13.448, 18.967, 21.602, 24.781, and 33.736 Å, respectively, in the stacked spectra.

4. ANALYSIS AND RESULTS

We now search for and measure the X-ray absorption lines in the stacked spectra. There are no extragalactic O vii, O viii, Ne ix, N vii, or C vi absorption lines apparent in any of the stacked spectra (Figures 46). Adding Gaussian profiles at the corresponding rest-frame wavelengths, we obtain upper limits to EWs of these lines. To obtain column densities of these ions, we first calculate the mean dispersion velocity ($\overline{b}$) and column density $\overline{N_{{\rm{O}{\mathsc{vi}}}}}$ of different classified O vi systems (Table 1) by weighting them with respect to the spectral counts at the wavelengths of the expected O vii Kα lines (Figure 2). Replacing the Gaussian profiles with the absorption-line model absline4 (Yao & Wang 2005) and fixing the b of different ions to $\overline{b}$, we then estimate the column density upper limits of these ions as reported in Table 3.

Table 3. The 95% Confidence Upper Limits of Kα Absorption Lines

line λ (Å) Ne ix 13.448 O viii 18.967 O vii 21.602 N vii 24.781 C vi 33.736
EW (mÅ)a 1.0(1.5) 2.2(2.6) 1.1(2.6) 1.2(2.0) 3.6(8.1)
log[N(cm−2)]a 15.04(15.29) 15.51(15.57) 14.62(15.12) 14.85(15.06) 15.06(15.71)
EW (mÅ)b 1.5(3.0) 2.3(3.0) 1.4(2.9) 1.3(3.0) 4.4(8.2)
log[N(cm−2)]b 15.28(15.71) 15.53(15.62) 14.75(15.15) 14.90(15.21) 15.15 (15.76)
EW (mÅ)c 3.8 3.1 3.95 4.3 8.0
log[N(cm−2)]c 15.75 15.63 15.19 15.33 15.36

Notes. Superscripts a, b, and c denote the upper limits obtained from stacked X-ray spectra corresponding to all, complex, and strong O vi absorbers, respectively. Values in parenthesis indicate limits constrained from the spectra without contribution from the PKS 2155-304 sight line.

Download table as:  ASCIITypeset image

With these upper limits to the ionic column densities and calculated $\overline{N_{{\rm{O}{\mathsc{vi}}}}}$, we can probe thermal properties of the O vi-bearing gas. For a gas in collisional ionization equilibrium (CIE), the column density ratio between O vii and O vi provides a diagnostic of gas temperature (Sutherland & Dopita 1993; Yao & Wang 2005). Assuming a single gas temperature, we jointly analyze the nondetected O vii with the $\overline{N_{{\rm{O}{\mathsc{vi}}}}}$ (Tables 1 and 3), and obtain a temperature upper limit of log [T(K)] < 5.7 at the 95% confidence level from the three spectra corresponding to all, complex, and strong O vi absorbers. Including the nondetected lines of Ne ix, N vii, and C vi in the joint analysis (assuming relative solar abundances of Ne/O, N/O and C/O; Anders & Grevesse 1989) does not further constrain the temperature. From the spectrum without the contribution of PKS 2155–304, we obtain the same upper limit. The reason why these spectra yield nearly the same temperature upper limit is the tightly constrained ratio of $N_{{\rm{O}{\mathsc{vii}}}}/N_{{\rm{O}{\mathsc{vi}}}}$ from the various spectra from ≲10 to ≲35 (see Tables 1 and 3), compared to 8.9, 26.2, and 63.1 for log [T(K)] = 5.6, 5.7, and 5.8, respectively, for gas in CIE (see Sutherland & Dopita 1993 and Figure 6 in Tripp et al. 2001).

5. DISCUSSION

We have presented a search for X-ray absorption lines produced in the WHIM by using the identified far-UV O vi absorbers as tracers for stacking the archived Chandra observations. The three final stacked spectra, corresponding to all, complex, and strong O vi absorbers, have signal-to-noise ratios of ∼32, 28, and 10, respectively, per 12.5 mÅ spectral bin around the O vii Kα wavelength. There are no detectable X-ray absorption lines at the expected wavelengths in these spectra. We have obtained upper limits to EWs of the Kα lines of O vii, O viii, Ne ix, N vii, and C vi and their column densities. Combining these nondetected lines with the average O vi column density, we have also derived an upper limit to the temperature of the O vi-bearing gas.

The X-ray measurements of the O vi systems with different physical properties in principle can reveal the ionization mechanism of the intervening gas. While most observed O vi absorbers could arise in interfaces between the low- and high-ionization phases (e.g., Danforth & Shull 2008; Tripp et al. 2008), compared to the velocity well aligned simple O vi absorbers, the dynamically complex O vi systems provide more direct evidence for the multiphase nature of the absorbers (e.g., Tripp et al. 2001, 2008; Shull et al. 2003; Sembach et al. 2004; Savage et al. 2005). Cosmological simulations of large-structure formation also predict that the O vi systems with larger equivalent widths are more likely to be collisionally ionized (Cen et al. 2001; Fang & Bryan 2001). If the O vi absorbers utilized in this work are indeed multiphase and trace the shock-heated IGM, the hot O vii-bearing gas is expected to be either surrounding the cool absorbers (traced by the absorption lines of Lyα, C ii, Si ii, etc.) or itself containing the observed O vi or a combination of both. The column densities of O vii (and/or O viii) are expected to be at least an order of magnitude more than that of O vi. But the current X-ray observations only constrain the $N_{{\rm O}{\mathsc{vii}}}/N_{{\rm O}{\mathsc{vi}}}\lesssim 10$ on average (Tables 1 and 3), which does not strongly favor the collisional ionization scenario of the O vi absorbers for shock temperatures T>105.7 K (Furlanetto et al. 2005). However, these data still cannot rule out the photoionization scenario, in which $N_{{\rm O}{\mathsc{vii}}}/N_{{\rm O}{\mathsc{vi}}}$ is expected to be ≲3 (see Figure 6 in Tripp et al. 2001 and Furlanetto et al. 2005).

There are several possibilities or their combinations that may cause the nondetections of the X-ray lines:

  • 1.  
    Some of the O vi systems may be photoionized or partly photoionized. The purely photoionized systems may not contribute to the baryon inventory (Section 1) since they could have already been counted through Lyα absorbers. Since the density of the IGM is very low, photoionization could play an important role in producing the highly ionized oxygen species (especially O vi5) in the WHIM besides the gravitational shocks (e.g., Chen et al. 2003; Tripp et al. 2008). In this case, the $N_{{\rm O}{\mathsc{vii}}}/N_{{\rm O}{\mathsc{vi}}}$ is expected to be lower than that in a solely collisionally ionized gas. Simulations show that, while the EWs of O vi and O vii are fairly well correlated in the putative WHIM in a broad overdensity range (δ∼10–100), for the O vi-absorbers with EW(O vi)>34 mÅ, there are ≲20% of the associated O vii-absorbers with EW(O vii)>2 mÅ (see Figures 13–15 in Chen et al. 2003 and Table 3).
  • 2.  
    Thermal properties of the WHIM may vary among different systems and sight lines. The O vii and/or O viii absorptions are expected to be strong in some systems and weak in others, but the "strong" X-ray absorption signals that originate from the shock-heated gas could have been diluted in the stacked spectra.
  • 3.  
    The characteristic temperature of the hot gas may be too high (e.g., T>3 × 106 K, as observed in the intragroup and intracluster medium) to produce observable O vii or O viii. In this case, the interface between hot and cool IGM could still be present to be responsible for the observed O vi absorbers, but the average temperature upper limit (log [T(K)] < 5.7) derived from the ratio of $N_{{\rm O}{\mathsc{vii}}}/N_{{\rm O}{\mathsc{vi}}}$ (Section 4) only applies to the vicinity of the O vi-bearing gas. This may be the case for the z = 0.0553 absorber along the PKS 2155-304 sight line where the O viii line was detected but the O vii was not (Section 1 and references therein).
  • 4.  
    Lastly, the O vi absorbers might only trace the warm part of the WHIM at temperatures T ≲ 3 × 105 K, and the hot (T ∼ 106 K) IGM is not necessarily always colocated with the warm gas. If this is the case, because of the small ionization fraction of O vi at high temperatures, the O vi absorption lines produced in the hot gas might not have been detected along most of the sight lines. The possible large dispersion velocity of the hot gas could also make the O vi absorption feature shallower and harder to detect. The high sensitivity of the Cosmic Origins Spectrograph, which is scheduled to be installed on HST in spring 2009, will provide a great opportunity to search for such weak O vi absorbers (at z>0.12) and to examine this possibility.

The results obtained in this work also have important implications regarding future observations of the WHIM with the current and the future X-ray observatories. With the still limited spectral resolution [even in the proposed International X-ray Observatory (IXO), EE ∼ 3000, which is still ∼10–15 times less than those of current far-UV instruments], it will be challenging to conduct a blind search for the X-ray absorption line "forest" produced in the WHIM. A more effective strategy is still to use the identified O vi and/or Lyα absorbers as references, as implemented in this work (also see Chen et al. 2003; Danforth & Shull 2005, 2008). However, if the six sight lines (Table 1) used in this work fairly sample the WHIM filaments over the whole sky, the gas only contains ≲1015cm-2 of O vii in column density and the O vii Kα EW ≲2.5 mÅ in absorption on average (Table 3). To detect a weak line with EW∼2 mÅ in an AGN spectrum with a flux of 4 × 10−12 erg s-1 cm−2 keV−1 around the line centroid, our 1000-run simulations for an exposure of 50 Ms with the Chandra ACIS-LETG only yield 375 detections at EWEW ≳ 3σ significance level. In comparison, for the IXO with EE ∼ 3000 and effective area of ∼3000 cm2, our simulations indicate that a 100 ks observation will have a 90% probability of detecting such a weak line at ≳3.5σ significance level.

We are grateful to Renyue Cen for useful discussion on nondetection of X-ray absorption lines. This work was made possible by Chandra archival research grant AR7-8014. Additional support for this research was provided by NASA grant NNX08AC14G, provided to the University of Colorado to support data analysis and scientific discoveries related to the Cosmic Origins Spectrograph on the Hubble Space Telescope, and partly through the Smithsonian Astrophysical Observatory contract SV3-73016 to MIT for support of the Chandra X-Ray Center under contract NAS 08-03060. T.M.T. and L.S. also acknowledge support for this work from NASA grant NNX08AJ44G.

Footnotes

  • In modeling a single absorption line, this model is similar to the curve-of-growth analysis. But it can be used to jointly analyze multiple absorption lines at the same time. For a detailed description of the model, see Yao & Wang (2005).

  • Higher ionization stages are harder to produce by photoionization, but we note that Hellsten et al. (1998) have argued that in the IGM, even O vii and O viii are predominantly photoionized by the X-ray background.

Please wait… references are loading.
10.1088/0004-637X/697/2/1784