Brought to you by:

The Impact of the Group Environment on the O vi Circumgalactic Medium

, , , , , and

Published 2017 July 18 © 2017. The American Astronomical Society. All rights reserved.
, , Citation Stephanie K. Pointon et al 2017 ApJ 844 23 DOI 10.3847/1538-4357/aa7743

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/844/1/23

Abstract

We present a study comparing O viλλ1031, 1037 doublet absorption found toward group galaxy environments to that of isolated galaxies. The O vi absorption in the circumgalactic medium (CGM) of isolated galaxies has been studied previously by the "Multiphase Galaxy Halos" survey, where the kinematics and absorption properties of the CGM have been investigated. We extend these studies to group environments. We define a galaxy group as having two or more galaxies with a line-of-sight velocity difference of no more than 1000 km s−1 and located within 350 kpc (projected) of a background quasar sightline. We identified a total of six galaxy groups associated with O vi absorption Wr > 0.06 Å that have a median redshift of $\langle {z}_{\mathrm{gal}}\rangle =0.1669$ and a median impact parameter of $\langle D\rangle =134.1\,\mathrm{kpc}$. An additional 12 non-absorbing groups were identified with a median redshift of $\langle {z}_{\mathrm{gal}}\rangle =0.2690$ and a median impact parameter of $\langle D\rangle =274.0\,\mathrm{kpc}$. We find the average equivalent width to be smaller for group galaxies than for isolated galaxies (3σ). However, the covering fractions are consistent with both samples. We used the pixel-velocity two-point correlation function method and find that the velocity spread of O vi in the CGM of group galaxies is significantly narrower than that of isolated galaxies (10σ). We suggest that the warm/hot CGM does not exist as a superposition of halos; instead, the virial temperature of the halo is hot enough for O vi to be further ionized. The remaining O vi likely exists at the interface between hot diffuse gas and cooler regions of the CGM.

Export citation and abstract BibTeX RIS

1. Introduction

The environment in which galaxies exist plays a significant part in the way they will evolve. While unperturbed isolated galaxies have a disk-like structure, interactions between galaxies can result in other interesting phenomena such as tidal streams, shells, and increased star formation rates (e.g., Barnes & Hernquist 1992; Veilleux et al. 2005; Poggianti et al. 2016). In the most extreme cases interactions remove the galaxy gas reservoir, leading to the quenching of star formation, and the original disk-like structures of the merging galaxies deform into elliptical galaxies (e.g., Gunn & Gott 1972; Cowie & Songaila 1977; Larson et al. 1980; Nulsen 1982; Moore et al. 1996; Cen & Ostriker 1999; Oppenheimer & Davé 2008; Lilly et al. 2013). These phenomena are observed in images taken of interacting galaxies with, for example, bursts of star formation activity (e.g., Keel et al. 1985; Barnes 2004) or streams of HI gas connecting galaxies that can be seen (e.g., Bridge et al. 2010). However, we know very little about the impact of a merger or interaction on the diffuse gas around galaxies.

Of particular interest is the circumgalactic medium (CGM), which is a vast halo of diffuse gas surrounding galaxies out to radii ≈200 kpc (e.g., Kacprzak et al. 2008; Chen et al. 2010; Steidel et al. 2010; Kacprzak et al. 2011; Tumlinson et al. 2011; Rudie et al. 2012; Burchett et al. 2013; Nielsen et al. 2013b, 2013a; Werk et al. 2013). The mass of this halo is comparable to the gas mass of the galaxy itself (i.e., the ISM; Thom et al. 2011; Tumlinson et al. 2011; Werk et al. 2013) and hence plays an important role in the evolution of galaxies. Models of the CGM indicate that gas can flow in and out of this reservoir, which in turn controls the star formation rate of the galaxy and the metallicity of stars formed from this gas (Oppenheimer & Davé 2008; Lilly et al. 2013; Kacprzak et al. 2016). While a more concrete model of how this gas drives isolated galaxy evolution is being built, we are still only just beginning to study the effects that an interaction or merger can have on the CGM.

Since the visible components of galaxies are clearly affected by galaxy interactions, it is reasonable to expect that such effects would also take place in the CGM. Indeed, due to the large radius of the CGM, two galaxies could have overlapping gaseous halos even if they do not yet show any signs that an interaction is taking place (Tully et al. 2009; Bordoloi et al. 2011; Stocke et al. 2014). Thus, it is possible to use the CGM to understand the very first processes that take place in an interaction. Similarly, as interactions often result in the ISM being stripped from the galaxy leading to gas depletion (e.g., Gunn & Gott 1972; Larson et al. 1980; Fujita & Nagashima 1999), the CGM could also be susceptible to gas removal through tidal interactions (Chen & Mulchaey 2009). This, in turn, would affect the ability of galaxies to form stars later in life as the accretion of gas from the CGM is what sustains star formation.

Previous studies of the CGM in group environments have used the Mg iiλλ2796, 2803 absorption doublet. Bordoloi et al. (2011) stacked background galaxy spectra and found that Mg ii absorption at a given equivalent width was located at higher impact parameters in group environments compared to isolated environments. This could be explained by constructing a simple model where the CGM of the group galaxies were superimposed. Thus the extended CGM is due to more gas along the line of sight, and hence, the Mg ii component of the CGM is minimally affected by the interactions in the group. Recent work by N. M. Nielsen et al. (2017, in preparation) has also compared Mg ii absorption between group and isolated environments. They found that the absorption associated with group galaxies had larger equivalent widths than for isolated environments. The authors also found an increased fraction of Mg ii components with higher velocities in group environments, which indicated that there could be some interaction between the CGM halos of the group galaxies.

Another ideal absorption doublet to trace the CGM is C iv, which traces gas with temperatures T ≈ 104–5 K. Burchett et al. (2016) compared C iv absorption in isolated and group/cluster environments and found that for a massive and dense group environment with seven group member galaxies (Mh > 1012.7 M and Mr ≤ −19), no C iv was detected despite associated H i absorbers. They find that the cause of C iv depletion in group environments is not clear due to the limited sample size. However, if larger studies of C iv absorption follow similar trends, the authors suggest that the CGM of the group galaxies may be experiencing ram pressure or tidal stripping, and that C iv traces gas with higher temperatures than HI absorbers.

The warmer, more diffuse phase of the CGM gas is traced by O vi due to its higher ionization state and temperature T ≈ 105.5 K. O vi can be both collisionally ionized and photoionized, which makes interpretation more difficult because it traces multiphase structures (Mo & Miralda-Escude 1996; Maller & Bullock 2004; Dekel & Birnboim 2006). Previous studies by Stocke et al. (2014, 2017) compared O vi absorption in group environments to that found in isolated environments. They found that O vi absorption profiles had to be modeled using fewer, broader components in group environments, indicating a warmer environment. However, the total absorption profile was narrower, from which the authors reasoned that O vi could not be distributed over the "circum-group" medium, especially since O vi is unstable due to rapid cooling. This suggests that a diffuse halo of warm CGM gas would be difficult to maintain, and that O vi exists at the interface of the diffuse hot (T > 106 K) and cooler photoionized regions embedded in the"circum-group" medium. The existence of interfaces between hot and cold gas within the CGM has also been found by Churchill et al. (2012) and Stern et al. (2016) for isolated environments.

Recently, Oppenheimer et al. (2016) used the EAGLE simulations to compare the column density of various oxygen ionization states as a function of halo mass and found that the presence and strength of O vi absorption are strongly related to the mass of the halo. Halos of sub-L* galaxies contain little O vi absorption because the virial temperature is too low to ionize oxygen to such a high state. As the halo mass increases, the presence and strength of O vi absorption increases until the virial temperature approaches the temperature where the ionization fraction of O vi peaks, which occurs in the CGM surrounding L* galaxies. Finally, as the halo mass increases and virial temperature beyond this point, the presence of O vi absorption decreases as oxygen is ionized to higher states. They find that oxygen is more likely to exist as O vii and O viii for more massive halos, in particular halos significantly larger than L*.

To investigate the simulation predictions from Oppenheimer et al. (2016) and compare to Mg ii and CIV absorption results, we use two O vi samples, one which contains only isolated galaxies and the other which contains galaxies in group environments. The isolated sample, which represents L* galaxies, has been presented in the "Multiphase Galaxy Halos" survey (Kacprzak et al. 2015; Muzahid et al. 2015, 2016; Nielsen et al. 2017). Here we present the group sample, representing halos larger than that of L* environments, which is the combination of data from the "Multiphase Galaxy Halos" survey and the literature (Ellingson et al. 1991; Chen et al. 2001; Chen & Mulchaey 2009; Prochaska et al. 2011; Werk et al. 2012, 2013, 2014; Johnson et al. 2013; Savage et al. 2014). A combination of HST images and HST/COS spectra allows us to investigate the equivalent widths, covering fractions, and kinematics of group environments and compare them to existing surveys of O vi absorption toward isolated galaxies.

This paper is organized as follows: In Section 2 we describe our sample of galaxies, quasars, and the selection criterion used to classify group and isolated environments. We also describe how we obtain a physical description of the CGM in such environments using quasar spectra. We present the results of our analysis with absorption properties and kinematics in Section 3. In Section 4 we discuss the differences between O vi absorption for group and isolated environments, and in Section 5 we summarize the work of this paper. We use standard a ΛCDM cosmology (Ho = 70 km s−1 Mpc−1, ΩM = 0.3, ΩΛ = 0.7).

2. Sample and Data

We identified group environments in a total of 10 fields. All but three were identified in quasar spectra from the Cosmic Origins Spectrograph (COS) on HST. The remaining three were from the FUSE telescope, and the spectra were provided by B. Wakker (2016, private communication). The quasar locations, redshifts, and the observation PIDs are shown in Table 1. The quasar spectra that sample foreground galaxy groups with O vi absorption are shown in columns two and four of Figure 1 as a black line.

Figure 1.

Figure 1. Sky locations and HST/COS or FUSE spectra absorption profiles for each group environment with O vi absorption. The absorbers J040748−121136 at zabs = 0.091917 and J111908+211918 at zabs = 0.059188 are from FUSE spectra. The relative positions of the galaxies in each group (red circles) are shown with respect to the locations of the quasar (black cross). The top profile shows the O vi λ1031 absorption line, while the bottom profile shows the O vi λ1037 line, where the data is shown in black, the error spectrum in green, and the O vi model in red. The center of each model component is indicated by a red tick. The blue line shows the modeled fit to the data where other transitions were modeled in addition to O vi due to the presence of blends. The shaded region contains the absorbing pixels within the velocity regions defined in Section 2.2.

Standard image High-resolution image

Table 1.  Quasar Observations

J-Name R.A. (J2000) Decl. (J2000) zqso Instrument Referencesa COS/FUSE PID
J004706+031955 00: 47: 05.93 +03: 19: 54.90 0.6233 COS (3) 12275
J012529−000556 01: 25: 28.84 −00: 05: 55.93 1.0748 COS (1) 13398
J022815−405714 02: 28: 15.17 −40: 57: 14.29 0.4934 COS (2) 11541
J035129−142909 03: 51: 28.54 −14: 29: 08.71 0.6163 COS (5) 13398
J040748−121137 04: 07: 48.43 −12: 11: 36.66 0.5726 COS, FUSEb (4) 11541, B087
J045609−215909 04: 56: 08.92 −21: 59: 09.40 0.5335 COS (7) 13398
J092838+602521 09: 28: 37.98 +60: 25: 21.02 0.2959 COS (3) 11598
J111909+211918 11: 19: 08.68 +21: 19: 18.01 0.1765 FUSEb (3) P101
J113328+032719 11: 33: 27.78 +03: 27: 19.17 0.5245 COS (6) 11598
J113910−135044 11: 39: 10.70 −13: 50: 43.64 0.5565 COS (1) 12275
J130112+590206 13: 01: 12.93 +59: 02: 06.75 0.4778 COS (8) 11541
J131956+272808 13: 19: 56.23 +27: 28: 08.22 1.0147 COS (8) 11667
J135704+191907 13: 57: 04.43 +19: 19: 07.37 0.7200 COS (8) 13398
J170441+604430 17: 04: 41.37 +60: 44: 30.50 0.3719 COS (8) 12276
J182157+642037 18: 21: 57.31 +64: 20: 36.37 0.2970 COS (8) 12038

Notes.

aReference for the position and redshift of the quasars are: (1) Beasley et al. (2002); (2) Beuermann et al. (1999); (3) Zickgraf et al. (2003); (4) Fey et al. (2004); (5) Li & Jin (1996); (6) Wu et al. (2012); (7) Healey et al. (2007); (8) Chen et al. (2001). bThese quasar spectra were obtained from B. Wakker (2016, private communication).

Download table as:  ASCIITypeset image

2.1. Galaxy Properties

We use a sample of galaxies located in group environments drawn from the "Multiphase Galaxy Halos" survey (Kacprzak et al. 2015; Muzahid et al. 2015, 2016; Nielsen et al. 2017) or obtained from the literature. A group is defined as having the nearest of two or more galaxies within 205 to 350 kpc of the quasar sightline, and the galaxies must have a line-of-sight velocity separation less than 1000 km s−1. With these criteria, we compiled a sample of six group environments with detected O vi absorption above a rest-frame equivalent width of Wr ≥ 0.06 Å in background quasar spectra, and an additional 12 group environments with no detected absorption. Two of these group environments detected with O vi absorption were identified by Savage et al. (2014) from FUSE data. The median redshift and impact parameter of the galaxies in the group environment sample with O vi absorption are $\langle {z}_{\mathrm{gal}}\rangle =0.1669$ and $\langle D\rangle =134.1$ kpc, respectively. In comparison, the median values for the equivalent width and impact parameter for group environments with no O vi absorption are $\langle {z}_{\mathrm{gal}}\rangle =0.2690$ and $\langle D\rangle =274.0$ kpc. We show the group galaxy redshift, position, and impact parameter in Table 2 as well as the same values for the non-absorbing groups. The values for the right ascension and declination offsets and the angular separation of the galaxies from the quasar sightline, as well at the impact parameter, were recalculated using methods from Nielsen et al. (2013b). In Figure 1 we show the locations of the group galaxy members associated with O vi absorption (red points) with relation to the quasar sightline (central black cross) in the first and third columns. Similarly, we show in Figure 2 the locations of galaxies that were not associated with O vi absorption (red points) relative to the location of the quasar sightline (central black cross).

Figure 2.

Figure 2. Sky locations for each group environment with no associated O vi absorption ${W}_{r}\lt 0.06\,\mathring{\rm A} $. The relative positions of the galaxies in each group (red circles) are shown with respect to the locations of the quasar (black cross).

Standard image High-resolution image

Table 2.  Group Galaxies

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
Quasar J-name zgal Referencesa Δα Δδ θ Db zabs Wrc log N (O vi)c vc v+c
      ('') ('') ('') (kpc)   (Å)   (km s−1) (km s−1)
Properties of Group Galaxies Associated with O vi Absorption (Wr > 0.06 Å)
J012328−000555 0.3787 5 −8.8 −12.3 15.07 78.3 0.378964 0.32 ± 0.04 14.60 ± 0.07 −158 143
  0.3792 3 27.7 −36.5 45.80 238.1          
  0.2065 4 −9.1 −8.4 10.87 36.8          
J022815−405714 0.2078 4 −24.9 −25.9 32.04 108.9 0.207058 0.199 ± 0.007 14.38 ± 0.01 −117 129
  0.2077 4 −69.9 63.2 82.32 279.7          
  0.0923 2 13.6 −39.9 42.08 72.3          
  0.0908 2 −61.9 13.6 62.01 104.9          
J040748−121136 0.0914 2 −78.9 −10.6 77.89 132.6 0.091924 0.07 ± 0.02 13.82 ± 0.06 −75 74
  0.0917 2 −123.5 −127.3 175.44 299.6
  0.0908 2 62.5 −252.3 259.64 439.4
J040748−121136 0.1670 3 −1.1 34.8 34.81 99.4 0.166984 0.245 ± 0.004d 14.70 ± 0.01d −170d 156d
  0.1669 3 41.3 −1.8 40.36 115.2          
J040748−121136 0.3614 2 −38.5 27.9 46.87 236.4 0.361484 0.070 ± 0.003 13.82 ± 0.01 −73 74
  0.3608 2 61.0 −19.3 62.71 315.9          
J111908+211918 0.0600 1 −48.1 −104.9 117.01 135.6 0.059184 0.06 ± 0.02 13.91 ± 0.13 −115 52
  0.0594 1 −83.1 −174.6 190.99 219.3          
Properties of Group Galaxies Associated with O vi Non-absorption (Wr < 0.06 Å)e
  0.3142 6 −4.0 63.0 63.13 290          
J004705+031954 0.3130 6 56.0 39.0 68.13 312 0.3135 <0.046 <13.61
  0.3133 6 61.0 71.0 93.54 429          
J004705+031954 0.3810 6 33.0 41.0 52.60 274 0.3798 <0.054 <13.68
  0.3786 6 −96.0 8.0 96.17 499          
  0.1992 4 84.0 24.6 68.04 220  
J022815−405714 0.1997 4 −110.6 −32.4 89.60 296 0.1998 <0.011 <12.94
  0.2006 4 133.7 78.4 127.88 419          
J022815−405714 0.2678 4 16.9 −13.0 18.21 70 0.2677 <0.016 <13.11
  0.2690 4 8.5 −36.7 37.25 150          
  0.2680 4 36.2 −29.2 39.98 160          
  0.2654 4 93.4 −8.5 71.03 284          
  0.2664 4 105.7 9.4 80.36 323          
  0.2706 4 108.4 15.4 83.28 339          
  0.2662 4 116.0 16.9 89.24 345          
  0.2683 4 61.1 −72.9 86.25 349          
  0.3236 3 13.0 −23.5 26.72 127          
J035128−142908 0.3244 3 −29.9 18.5 34.33 164 0.3251 <0.032 <13.41
  0.3273 3 −59.0 19.5 60.31 288          
  0.3506 2 −26.0 −40.3 47.66 235          
J040748−121136 0.3521 2 37.0 55.7 66.40 329 0.3515 <0.004 <12.51
  0.3517 2 80.5 −141.3 161.78 802          
  0.4838 3 −1.1 −18.0 18.03 108          
J045608−215909 0.4836 3 −20.0 −16.6 24.88 284 0.4837 <0.032 <13.73
  0.4837 3 40.1 −36.8 52.29 315          
J085334+434902 0.0915 4 −1.8 40.8 40.81 53 0.0909 <0.054 <13.68
  0.0903 7 14.1 −34.0 35.45 79          
J092554+400414 0.2475 8 −8.0 −20.8 21.64 81 0.2471 <0.085 <13.91
  0.2467 8 −7.2 −24.1 24.69 92          
  0.1537 8 −3.5 −14.7 14.82 38          
J092837+602521 0.1542 8 30.2 −12.1 19.19 48 0.1540 <0.111 <14.06
  0.1540 8 67.2 −12.3 35.38 89          
J113327+032719 0.2367 8 4.5 −1.7 4.79 18 0.2365 <0.067 <13.79
  0.2364 8 −4.1 −9.6 10.39 35          
  0.3598 3 −52.4 13.3 52.57 265          
J113910−135043 0.3604 3 −28.0 51.7 58.41 296 0.3599 <0.027 <13.36
  0.3595 3 −45.9 38.5 58.95 298          

Notes.

aGalaxy identification reference: (1) Prochaska et al. (2011); (2) Johnson et al. (2013); (3) Chen et al. (2001); (4) Chen & Mulchaey (2009); (5) Muzahid et al. (2015); (6) Ellingson et al. (1991); (7) Lanzetta et al. (1995); (8) Werk et al. (2012). bWe include known galaxies with D > 350 kpc for completeness. cUnless otherwise noted, values were calculated from the O vi λ1031 line. dValues were calculated from the O vi λ1037 line. eWe quote the 3σ limits for Wr and log N (O vi).

Download table as:  ASCIITypeset images: 1 2

We compare this group environment sample to the isolated galaxy sample presented in Kacprzak et al. (2015), who identified 53 isolated galaxies associated with 29 absorbers and 24 non-absorbers. However, the Kacprzak et al. (2015) sample only extends out to 200 kpc. We supplement this sample with an additional 14 isolated galaxies from the "Multiphase Galaxy Halos" survey to extend the impact parameter range out to 350 kpc. These additional galaxies are presented in Table 3. The galaxies in this sample have a median redshift and impact parameter of $\langle {z}_{\mathrm{gal}}\rangle =0.2608$ and $\langle D\rangle =259\,\mathrm{kpc}$, respectively. We find that the galaxies in this extended impact parameter range are all identified as O vi non-absorbers.

Table 3.  Properties of Isolated Galaxies with 200 kpc ≤ D ≤ 350 kpc

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11)
Quasar J-name zgal R.A. Decl. Referencesa Δα Δδ θ D Wr logN(O vi)
    (J2000) (J2000)   ('') ('') ('') (kpc) (Å) (kpc)
J012528−000555 0.4299 01: 25: 31.548 −00: 06: 04.22 1 40.6 −8.3 41.46 233.6 <0.061 <13.74
J022815−405714 0.1248 02: 28: 10.244 −40: 55: 48.24 2 −73.8 86.1 102.53 208.6 <0.019 <13.18
J022815−405714 0.2804 02: 28: 11.316 −40: 56: 36.94 2 −57.7 37.4 57.41 248.1 <0.011 <12.95
J040748−121136 0.2032 04: 07: 42.680 −12: 11: 31.87 3 −86.3 4.8 84.45 282.0 <0.005 <12.61
J040748−121136 0.2978 04: 07: 50.585 −12: 12: 24.23 3 32.3 −47.6 57.10 257.9 <0.005 <12.61
J040748−121136 0.4100 04: 07: 45.851 −12: 12: 13.82 3 −38.7 −37.2 53.03 291.5 <0.006 <12.69
J100402+285535 0.2143 10: 04: 03.600 +28: 54: 33.00 1 14.8 −62.4 63.72 221.6 <0.020 <13.22
J113910−135043 0.4007 11: 39: 12.076 −13: 51: 37.45 1 20.6 −53.8 57.41 310.6 <0.059 <13.73
J130112+590206 0.2412 13: 01: 09.806 +59: 03: 16.35 1 −46.9 69.6 73.67 279.9 <0.017 <13.15
J131956+272808 0.6717 13: 19: 59.582 +27: 28: 06.86 1 50.3 −1.4 44.63 316.2 <0.044 <13.59
J135704+191907 0.5293 13: 57: 03.470 +19: 18: 19.59 1 −14.4 −47.8 49.67 313.8 <0.098 <13.99
J170441+604430 0.1877 17: 04: 51.297 +60: 44: 38.66 1 148.8 8.2 73.19 230.9 <0.031 <13.42
J170441+604430 0.2260 17: 04: 31.779 +60: 44: 42.73 1 −144.0 12.2 71.41 259.9 <0.022 <13.26
J182157+642037 0.2110 18: 22: 02.759 +64: 21: 38.92 4 81.7 62.6 71.85 248.4 <0.005 <12.63

Note.

aGalaxy identification reference: (1) Chen et al. (2001); (2) Chen & Mulchaey (2009); (3) Johnson et al. (2013); (4) Stocke et al. (2013).

Download table as:  ASCIITypeset image

2.2. Modeling

We performed a Voigt profile analysis using the software VPFIT6 (Carswell & Webb 2014) to fit the absorption profiles of the O vi absorbers. As the COS instrument on the HST has a non-Gaussian line spread function (LSF), leading to significant wings that are prominent at short wavelengths, we use the non-Gaussian LSF from Kriss (2011). The LSF for each absorption profile was calculated by interpolating the LSF from Kriss (2011) to the central wavelength of the profile. This LSF was then convolved with the model profile in the fitting process. The velocity resolution of the FUSE observations is ∼15–20 km s−1 (FWHM) and we assume a Gaussian LSF.

Where possible, both lines of the O vi λλ1031, 1037 doublet were used to fit the absorption profile to better constrain the number of components required and determine any blends that might be present. The number of components was based on the minimum number required to produce a reasonably minimized chi-squared value. We present the results of the fitting in the second and fourth columns of Figure 1. The black line represents the data, while the red line is the model of the O vi absorption profile made using VPFIT. Red ticks indicate the central position of each component used to fit the O vi absorption profile. Three absorbers were found to have blends present. In these cases, we used the unblended transition to fit the profile. The blended transition was then modeled by adding generic components. The total fit in these two instances is shown as a blue line.

From the fitting process, we extracted the column density, Doppler b parameter, and redshift of each component of the absorption profiles. We then calculated the equivalent width by integrating the absorption features between the points where the model deviated from the continuum, which we defined as the edges of the absorption profile. These were converted into velocity boundaries using zabs, which is the median velocity of the apparent optical depth distribution of the absorption and defines the velocity zero point (Churchill 1997). The regions contained by these velocity boundaries are shown as the shaded regions in Figure 1. We list the equivalent width, column density, zabs, and the absorption velocity range in Table 2.

We also calculated 3σ upper limits on the equivalent widths and column densities of the non-absorbing galaxies and group environments. To compute the column density limits, we assumed a single cloud with a Doppler b parameter of 30 km s−1, which we found to be typical of O vi absorption in our sample. The equivalent width and column density limits for group environments and isolated galaxy are shown in Tables 2 and 3, respectively.

3. Results

Here we present an analysis of the equivalent width, covering fractions and kinematics of both group and isolated galaxies to determine whether environment plays a role in shaping the properties of the CGM.

3.1. Equivalent Width and Covering Fraction

Kacprzak et al. (2015) found that the equivalent width decreases as the distance of the absorber from the galaxy increases for the isolated sample. Similarly, Tumlinson et al. (2011) and Johnson et al. (2015) found that the column density decreases as the ratio of the impact parameter and the virial radius increases. We investigate this trend in group galaxies by showing the equivalent width as a function of the impact parameter in Figure 3. In group environments, the impact parameter is taken to be the projected distance from the quasar to the nearest galaxy. The nearest group members are shown in orange, with circles indicating measurements of the equivalent width and crosses representing limits. In our sample, we detect absorbers with equivalent widths greater than 0.06 Å, which we choose to be the cut-off between absorbers and non-absorbers, represented as a black dashed line. The nearest group galaxies also appear to exhibit decreasing equivalent widths with increasing impact parameter, similar to the trend found in isolated environments. We also investigate the effect of measuring the impact parameter to the geometric center and the most massive (luminous) galaxy and find that the variations did not change the relationship significantly. Figure 3 suggests that the equivalent width of group galaxies is lower than what is expected from isolated galaxies.

Figure 3.

Figure 3. Equivalent widths of each O vi λ1031 profile as a function of impact parameter to the nearest member galaxy. The nearest group galaxies are shown in orange and isolated galaxies in purple. Circles indicate where measurements of the equivalent width were possible, while crosses show where only a limit could be obtained. The black dashed line indicates the cut-off between absorbers and non-absorbers (0.06 Å).

Standard image High-resolution image

We compare the equivalent width distributions of the group and isolated galaxies in Figure 4. The peak of the equivalent width distribution of group environments appears to be shifted lower than the peak for isolated galaxies. The average equivalent width for group environments is $\langle {W}_{r}\rangle =0.16\pm 0.02\,\mathring{\rm A} $, while the average equivalent width for isolated galaxies is $\langle {W}_{r}\rangle =0.31\pm 0.03\,\mathring{\rm A} $. This is a 3σ difference in the mean equivalent widths. We also test the removal of the outlier from the isolated galaxy sample at the top of the plot in Figure 3 due to its unusually large equivalent width (Muzahid et al. 2015). The average equivalent width for isolated galaxies then is $\langle {W}_{r}\rangle =0.28\pm 0.02\,\mathring{\rm A} $. The difference between the means of the equivalent widths of the group and isolated environments is then 2.4σ.

Figure 4.

Figure 4. Equivalent width distribution of group absorbers (orange histogram) compared to the isolated absorbers (purple histogram). The group absorbers may be offset from the main peak of the isolated absorbers.

Standard image High-resolution image
Figure 5.

Figure 5. Pixel-velocity two-point correlation function (TPCF) comparing the velocity spread of the O vi λ1031 absorption profiles for the group and isolated environments. The group galaxy TPCF is shown in orange, while the isolated galaxy TPCF is shown in purple. The result of the chi-squared test is 10σ, indicating that the O vi velocity spread of group absorption profiles is significantly narrower than the isolated TPCF. The shaded region shows the 1σ errors calculated from bootstrapping.

Standard image High-resolution image

The covering fraction, which is the ratio of absorbers to the total number of systems, is an indication of a lack of gas or ideal ionizing conditions in the CGM. We find that the covering fraction for group galaxies is ${C}_{f}={0.33}_{-0.11}^{+0.06}$, while the isolated galaxies have ${C}_{f}={0.43}_{-0.06}^{+0.03}$, where the 1σ errors are derived from 10,000 bootstrap realizations. These are consistent within the errors of the measurements (0.8σ), although the smaller covering fraction for the group galaxies suggests that the CGM in these environments appears to have a lack of gas or insufficient ionizing conditions compared to isolated environments. Lower equivalent widths and covering fractions for group galaxies could indicate that these environments are being depleted of O vi and a larger sample would be invaluable to investigate this effect.

3.2. Pixel-velocity Two-point Correlation Function

To compare the kinematics of the O vi absorbers for group and isolated environments, we use the pixel-velocity two-point correlation function (TPCF) method from Nielsen et al. (2015, 2016, 2017). To construct the isolated galaxy TPCF, Nielsen et al. (2017) compiled the pixel velocities within the absorption velocity bounds (i.e., shaded regions in Figure 1) for every isolated absorber in the sample. With this list, they calculated the absolute value of the velocity separations between every pixel pair, binned these velocity separations into 20 km s−1 bins, and normalized the distribution by the total number of pixel pairs in the sample. This produced a probability distribution function, which describes the probability that any two pixels will have a given velocity separation. They used a bootstrap analysis with 100 realizations to obtain 1σ uncertainties on the TPCF. We use the same method to calculate the group TPCF.

The group galaxy TPCF (orange) and the Nielsen et al. (2017) isolated galaxy TPCF (purple) are shown in Figure 5, where the TPCF is the thick solid line and the shaded regions represent the one sigma errors calculated from the bootstrap analysis. Any two pixels for a typical group absorption environment are more likely to be closer together in velocity space than two randomly selected pixels from an isolated galaxy. A chi-squared test run on the binned TPCFs, including the uncertainties, found that the null hypothesis that the two samples were drawn from the same population can be ruled out at the 10σ level. Thus, the O vi absorption velocity dispersion is smaller in group environments than in isolated environments.

To characterize the differences in the TPCFs quantitatively, we calculate Δv(90) and Δv(50) values, which are defined as the velocity separations within which the integrated TPCF curve contains 90% and 50% of the data, respectively. Thus, Δv(90) and Δv(50) are also measures of the velocity spreads of the distributions. We find that Δv(90) is ${153}_{-18}^{+21}$ km s−1 for group galaxy absorption, compared to 236 ± 15 km s−1 for isolated galaxy absorption (2.3σ). Similarly, the group absorption Δv(50) is ${64}_{-7}^{+9}$ km s−1, while the isolated galaxy absorption is 100 ± 6 km s−1 (2.4σ). Even though the uncertainties on these values do not overlap, they are still only 2.4σ different. Thus, combined with the chi-squared test of the TPCF profiles, which is more sensitive to velocity differences in the absorption profiles (Nielsen et al. 2015, 2016, 2017), we find that group environments are more likely to have narrower absorption line profiles than isolated environments.

4. Discussion

The UV spectra from HST/COS and FUSE have allowed us to compare the O vi CGM in the environment surrounding groups of galaxies to the CGM of isolated galaxies. We have investigated the effects of group environments on the equivalent widths, the covering fractions, and the velocity spread of the O vi absorption profiles.

Differences between the CGM in group and isolated environments have been explored using Mg ii. Bordoloi et al. (2011) found that for a given equivalent width, the absorption associated with group galaxies exists out to larger impact parameters than for isolated galaxies. This is supported by Chen et al. (2010), who did not find a strong inverse correlation between equivalent width and impact parameter. Bordoloi et al. (2011) suggested that this is a result of the CGM of group members being superimposed with little interaction between the halos. N. M. Nielsen (2017, in preparation) found that the velocity spread of absorption profiles for group environments is larger than the velocity spread for isolated galaxies. They also found that a higher fraction of high-velocity components (relative to zabs) exist in group halos. This could suggest some interaction occurring between the halos of the group galaxies. We study O vi absorption to investigate the influence of group environments on the warm/hot CGM.

Using the EAGLE simulations, Oppenheimer et al. (2016) found interesting results with mass that may explain the difference between O vi absorption from the CGM surrounding group and isolated environments. By examining simulated O vi column densities, the authors found that halo mass, and therefore virial temperature, governs the presence and strength of O vi, where sub-L* halos do not have sufficient virial temperatures to create significant amounts of O vi. The most favorable halo masses, and hence virial temperatures, for the creation of O vi are from L* galaxies, which are represented by our isolated sample. As the halo mass increases above that of a typical L* galaxy to group and cluster environments, the amount of O vi absorption decreases due to the temperatures in the halo preferentially ionizing O vi to higher states, such as O vii, which are difficult to observe. In the case of studying group galaxies, the halo mass is certainly larger that this limit, hence we should start to see this effect.

The effects of group environments on absorption properties such as covering fractions and equivalent widths were investigated to determine if the trends suggested by Oppenheimer et al. (2016) were visible in our sample. Taking these results into account, we expected that the equivalent widths and covering fractions for our group sample would be lower than for isolated sample if the virial temperatures were higher in the group sample. The average equivalent width of the group sample is lower than the isolated galaxies (3σ). However, the covering fractions of the absorption associated with group environments, while lower than what were found for isolated galaxies, are consistent with the values from Kacprzak et al. (2015; 0.8σ). This suggests that the mass trends with the presence of O vi could be present but are limited by the small number of group systems in this study. The unchanged covering fraction for group samples may be due to interactions between the group members stripping gas from the galaxies into the CGM. Further investigations should focus on improving the sample size of group environments to determine if groups environments have significantly smaller covering fractions.

The effects of group environments on C iv was investigated by Burchett et al. (2016), who found that groups with more galaxy members have a deficiency of C iv absorption. Similarly, we find that the larger groups of galaxies tend to have little or no O vi absorption. In our study, the largest group galaxy environment (with eight members), J022815−405714 at zabs = 0.2677, had an upper limit of Wr = 0.016 Å and was classed as a O vi non-absorber. Additionally, the second largest group (with five members), J040748−121136 at zabs = 0.0919, was a weak O vi absorber with Wr = 0.08 Å. Thus, it is possible that the gas traced by O vi in dense group environments could also be undergoing ram pressure or tidal stripping, leading to the gas being removed from the CGM. Future studies could focus on building a complete sample of groups with a range of member galaxies to investigate this trend.

While the effect of group environments on the CGM was less evident in the covering fractions, we further support the differences in the average equivalent width of the two samples with the kinematics, which were found to be significantly different from isolated environments. The TPCFs showed that, on average, group environment absorption is narrower than isolated galaxy absorption with a difference between the two populations of 10σ. To quantify this, Δv(90) for the group environments $({153}_{-18}^{+21}$ km s−1) is roughly two-thirds that of the isolated sample (236 ± 15 km s−1). Stocke et al. (2014) similarly found that the total absorption profiles of O vi absorbers were narrower for group environments. However, the individual clouds within the profiles were represented by broader components, which indicates higher virial temperatures. This narrowing of the absorption profiles, along with potentially smaller average equivalent widths and covering fractions, strongly suggests that the warm CGM does not exist as a superposition of halos as found in Mg ii by Bordoloi et al. (2011). This is further supported by Stocke et al. (2014), who found that O vi absorbers were not associated with the nearest galaxy in the group environment and favored a mixed CGM scenario. If the group environment CGM properties were due to a superposition of halos, we would expect to see similar or larger velocity widths in group environments compared to the isolated galaxies, which we did not observe.

Similarly to Stocke et al. (2014), we suggest that the regions of warm O vi CGM gas would exist at the interface between the higher virial temperature, diffuse medium of the combined group CGM, and the cooler gas associated with individual group galaxies. We suggest that the narrow absorption profiles of O vi in group environments is due to the ionization of lower column density gas in the higher velocity wings of the absorbers. Thus, the O vi gas regions are ionized from the outside of the region's mass, effectively forming an "onion skin" effect, with higher ionization states on the outside and decreasing inwards. This would create a shell of higher ionized oxygen surrounding the O vi, leading to the observed smaller velocity spread. Similar ideas have been found by Shull et al. (1998, 2003), Tripp et al. (2001), Aracil et al. (2006), Churchill et al. (2012), Kacprzak et al. (2013), Stocke et al. (2014), Stocke et al. (2017), and Stern et al. (2016). Thus, narrower profiles could be the result of a higher virial temperature, causing regions of O vi to be ionized from the outside inward.

The destruction of O vi and its existence at interface regions means that O vi is an interesting probe of the processes which occur in group environments. More information of the diffuse medium surrounding the interface could be found from higher ionization states such as O vii (Shull et al. 2003). However, at this redshift, O vii is in the x-ray regime and not detectable with current telescopes. Since O vi accounts for such a small fraction of the total baryons and oxygen budget (Peeples et al. 2014; Oppenheimer et al. 2016) investigating higher ionization states such as O vii, O viii, Ne viii, and Mg x could further the understanding of the total baryon contents of the CGM.

5. Summary and Conclusions

We compare the CGM for a total of six group galaxy environments to the CGM of 29 isolated galaxies from Kacprzak et al. (2015) associated with O vi absorption. We also have another 12 groups and 38 isolated galaxies with only an upper limit on O vi absorption. Group environments were classified as having at least two galaxies with a line-of-sight velocity separation of less than 1000 km s−1, and the nearest member located within 20–350 kpc of the quasar sightline.

We modeled the absorption profiles for the group sample using VPFIT to obtain the column densities, Doppler (b) parameters, and redshifts, from which we calculated the equivalent width and velocity spreads. We directly compared the equivalent widths and covering fractions of the two samples, and used the pixel-velocity TPCF method from Nielsen et al. (2015, 2016) to compare the velocity spreads of the group and isolated galaxy samples. Our findings are as follows.

  • 1.  
    The equivalent width of the absorption in group environments follows the previously established trend of decreasing as the impact parameter increases. Larger samples of group galaxies could determine if an anti-correlation exists similar to that found for isolated galaxies.
  • 2.  
    The average equivalent width of absorption associated with group environments is 0.16 ± 0.02 Å compared to the larger value for isolated galaxies, which is 0.31 ± 0.03 Å. The difference between the two populations is 3σ, which is significant.
  • 3.  
    The covering fraction is ${C}_{f}={0.33}_{-0.17}^{+0.06}$ for group galaxies and ${C}_{f}={0.43}_{-0.06}^{+0.03}$ for isolated galaxies. The covering fraction of absorption associated with group environments is smaller than that for isolated galaxies, but they are consistent within the errors of the populations with a difference of only 0.8σ.
  • 4.  
    The velocity spread of the absorption profiles of group environments was found to be significantly narrower than the velocity spread of the absorption profiles of isolated galaxies using the TPCF analysis, and a Chi-squared test found that this result is significant at the 10σ level. A comparison of the velocity spread between the group sample (Δv(90) = 154 ± 19 km s−1, Δv(50) = 66 ± 8 km s−1) and the isolated sample (Δv(90) = 236 ± 15 km s−1, Δv(50) = 100 ± 6 km s−1) found that the difference in these values between the two samples was at least 2.3σ.

The narrow velocity spreads in group environments indicate that O vi exists at boundaries between the more highly ionized and diffuse CGM and the cooler gas, which is in agreement with Shull et al. (1998, 2003) and Stocke et al. (2014, 2017). We find that is it unlikely that the warm and hot CGM exists as a superposition of the galaxy group member halos, unlike the scenario suggested by Bordoloi et al. (2011) for the cool CGM. This idea is consistent with Oppenheimer et al. (2016), who found that O vi is strongly dependent on mass as more massive halos have larger virial temperatures, indicating that the CGM in group environments is warm enough to ionize oxygen to O vii or higher. The lower average equivalent width and slightly smaller covering fraction for group environments also suggests that warm gas is being ionized further in group environments.

The O vi absorption toward group environments is markedly different from that detected toward isolated galaxies, although we would benefit from a larger sample size to better understand some of the differences. A study of the higher ionization states of the CGM in group environments would help us to understand the processes which are occurring in group environments; in particular, an investigation of the narrowing of the O vi absorption profiles and the non-superposition of warm CGM halos surrounding group galaxies. Such a study is not currently feasible, as these features would be found in the x-ray regime. However, investigating other high ionization states, located in the UV, would also help to broaden our understanding of the CGM in group environments.

Support for this research was provided by NASA through grants HST GO-13398 from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. S.K.P. acknowledges support through the Australian Government Research Training Program Scholarship. G.G.K. acknowledges the support of the Australian Research Council through the award of a Future Fellowship (FT140100933). S.M. acknowledges support from European Research Council (ERC), Grant Agreement 278594-GasAroundGalaxies. C.W.C. and J.C. acknowledge the support from award 1517831 from the National Science Foundation.

Footnotes

Please wait… references are loading.
10.3847/1538-4357/aa7743