MAGiiCAT V. ORIENTATION OF OUTFLOWS AND ACCRETION DETERMINE THE KINEMATICS AND COLUMN DENSITIES OF THE CIRCUMGALACTIC MEDIUM

, , , , and

Published 2015 October 9 © 2015. The American Astronomical Society. All rights reserved.
, , Citation Nikole M. Nielsen et al 2015 ApJ 812 83 DOI 10.1088/0004-637X/812/1/83

0004-637X/812/1/83

ABSTRACT

We investigate the dependence of gas kinematics and column densities in the ${\rm{Mg}}\;{\rm{II}}$-absorbing circumgalactic medium on galaxy color, azimuthal angle, and inclination to trace baryon cycle processes. Our sample of 30 foreground isolated galaxies at $0.3\lt {z}_{{\rm{gal}}}\lt 1.0,$ imaged with the Hubble Space Telescope, is probed by background quasars within a projected distance of $20\lt D\lt 110\;{\rm{kpc}}.$ From the high-resolution (${\rm{\Delta }}v\simeq 6.6$ km s−1) quasar spectra, we quantify the extent of the absorber velocity structure with pixel-velocity two-point correlation functions. Absorbers with the largest velocity dispersions are associated with blue, face-on ($i\lt 57^\circ $) galaxies probed along the projected minor axis (${\rm{\Phi }}\geqslant 45^\circ $), while those with the smallest velocity dispersions belong to red, face-on galaxies along the minor axis. The velocity structure is similar for edge-on ($i\geqslant 57^\circ $) galaxies regardless of galaxy color or azimuthal angle, for red galaxies with azimuthal angle, and for blue and red galaxies probed along the projected major axis (${\rm{\Phi }}\lt 45^\circ $). The cloud column densities for face-on galaxies and red galaxies are smaller than for edge-on galaxies and blue galaxies, respectively. These results are consistent with biconical outflows along the minor axis for star-forming galaxies and accreting and/or rotating gas, which is most easily observed in edge-on galaxies probed along the major axis. Gas entrained in outflows may be fragmented with large velocity dispersions, while gas accreting onto or rotating around galaxies may be more coherent due to large path lengths and smaller velocity dispersions. Quiescent galaxies may exhibit little-to-no outflows along the minor axis, while accretion/rotation may exist along the major axis.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The circumgalactic medium (CGM) is generally defined as the bound gaseous halo surrounding galaxies that extends out to a few hundred kiloparsecs (e.g., Kacprzak et al. 2008, 2013; Steidel et al. 2010; Tumlinson et al. 2011, 2013; Rudie et al. 2012; Nielsen et al. 2013a, 2013b; Werk et al. 2013). This region has increasingly been found to host some of the most important mechanisms involved in galaxy evolution through the baryon cycle, including gas accretion via the intergalactic medium and/or recycled accretion via the galactic fountain mechanism, galactic-scale outflowing winds, and merging satellite galaxies. The diffuse, multiphase nature of the CGM lends itself to study by way of absorption lines found in bright, background objects such as quasars or galaxies, or even in the spectrum of the host galaxy itself (i.e., the "down-the-barrel" approach).

Most work examining the low-ionization, cool (T ∼ 104 K) component of the CGM has been focused on ${\rm{Mg}}\;{\rm{II}}\lambda \lambda 2796,2803$ doublet absorption in background quasar spectra, as it is easily observed in the optical at redshifts $0.1\lt z\lt 2.5$ (e.g., Bergeron & Boissé 1991; Steidel et al. 1994, 1997; Guillemin & Bergeron 1997; Churchill et al. 2005; Barton & Cooke 2009; Kacprzak et al. 2011; Lan et al. 2014), and for a range of ${\rm{H}}\;{\rm{I}}$ column densities ($16\leqslant \mathrm{log}N({\rm{H}}\;{\rm{I}})\leqslant 22,$ e.g., Bergeron & Stasińska 1986; Steidel & Sargent 1992; Churchill et al. 2000; Rao & Turnshek 2000; Rigby et al. 2002). The ${\rm{Mg}}\;{\rm{II}}$ ion is well-known to be a tracer of the steps involved in the baryon cycle including accretion, rotating material merging onto the galaxy, and outflows. ${\rm{Mg}}\;{\rm{II}}$ absorbing gas may also be associated with merging satellites, which are in the process of being tidally stripped and/or have ongoing star formation driven outflows.

Accreting gas in the form of filaments from the cosmic web and/or recycled accretion from past outflows has been found to lie near the plane of the galaxy disk both in observations (e.g., Steidel et al. 2002; Kacprzak et al. 2010; Martin et al. 2012; Rubin et al. 2012; Bouché et al. 2013) and in simulations (e.g., Stewart et al. 2011; Ford et al. 2014). Simulations have also shown that this material forms an extended (out to $\sim 0.3{R}_{{\rm{vir}}}$), warped disk that co-rotates with the galaxy when viewed in edge-on orientations (e.g., Stewart et al. 2011; Danovich et al. 2012, 2015). Direct observations of infalling material have been few due to the small covering fraction of the accreting material (at least $6\%;$ e.g., Martin et al. 2012; Rubin et al. 2012) and because outflows dominate the absorption profile, though the spectra used in these works have low resolution with ${\rm{\Delta }}v\gt 150\;\mathrm{km}\;{{\rm{s}}}^{-1},$ which may contribute to the low detection rate. Nonetheless, accretion and rotation signatures may include velocities that are bound to the host galaxy, but are greater than or comparable to projected rotational velocities of the disk for edge-on galaxies probed along the projected major axis. The velocities of these absorbers likely lie to one side of the galaxy's systemic velocity (Steidel et al. 2002; Kacprzak et al. 2010; Stewart et al. 2011; Bouché et al. 2013), especially when viewing galaxies in nearly edge-on inclinations.

Based on low-resolution (${\rm{\Delta }}v\gt 150$ km s−1) spectra, outflows from galactic-scale winds due to star formation and/or supernovae feedback are often invoked to explain the presence of ${\rm{Mg}}\;{\rm{II}}$ absorption (e.g., Rubin et al. 2010, 2014; Bouché et al. 2012; Martin et al. 2012; Bordoloi et al. 2014a, 2014b; Kacprzak et al. 2014). A down-the-barrel approach was used in Weiner et al. (2009) and Rubin et al. (2010), and these authors show that in stacked galaxy spectra, outflows are commonly observed in face-on orientations due to the observed blueshift of material being pushed out of the host galaxy in the direction of the observer. Using both a down-the-barrel sightline and a transverse sightline through a single galaxy, Kacprzak et al. (2014) find that the sightlines are kinematically coupled and associated with outflowing material despite the large projected distance (58 kpc) between the two. In this case, the transverse quasar sightline is aligned with the galaxy projected minor axis. Finally, by modeling galactic winds to reproduce the absorption found in seven transverse sightlines, most of which are aligned with the galaxy projected minor axis, Bouché et al. (2012) were able to estimate properties of the winds such as wind speeds. All of this work is consistent with the picture of biconical, polar outflows whose signatures include broad, complex absorption profiles spanning hundreds of km s−1 that cannot be explained by rotation (Veilleux et al. 2005). Such a biconical outflow model is also applicable to our own Milky Way Galaxy (Fox et al. 2015).

Merging satellite galaxies may also be a source of ${\rm{Mg}}\;{\rm{II}}$ absorption in the CGM. For $z\gt 0.3$ the low luminosities of satellite galaxies makes it difficult to directly observe such satellites. In fact, Martin et al. (2012) found a redshifted absorber with respect to the targeted galaxy and associated it with a satellite that was clearly detected in their galaxy images, but this was the only case out of their sample of over 200 galaxies, indicating that satellites as the source of absorption is likely a rare scenario ($\lt 1\%$ probability). Examining the azimuthal angle dependence of satellites around host galaxies, Yang et al. (2006) found that satellites within D = 700 kpc tend to be isotropically distributed around blue galaxies, but are preferentially located along the major axis of red galaxies. This result has been found with a variety of data sets, which are summarized in Yang et al. (2006), as well as in simulations (most recently by Dong et al. 2014).

Detailed absorber kinematics for a large absorber–galaxy sample, which are required to provide insight into the motions of the gas involved in the baryon cycle and constrain simulations, have thus far remained elusive due to the need to stack low-resolution spectra. Even without needing to stack the spectra, much of the velocity structure is washed out by the low resolution. Most of the previous work examining ${\rm{Mg}}\;{\rm{II}}$ absorbers associated with baryon cycle processes have reported that the ${\rm{Mg}}\;{\rm{II}}$ equivalent width depends strongly on inclination and/or azimuthal angle (e.g., Bordoloi et al. 2011, 2014a, 2014b; Kacprzak et al. 2011, 2012; Bouché et al. 2012); however, equivalent width does not provide a detailed picture of the gas kinematics. This is especially true since equivalent width depends on the velocity spread of absorption and the column densities, which in turn depend on the line of sight geometry, metallicity, and ionization conditions of the absorbers. We therefore make use of a subsample of galaxies in the ${\rm{Mg}}\;{\rm{II}}$ Absorber–Galaxy Catalog (${\rm{MAG}}{\rm{II}}{\rm{CAT}}$; Churchill et al. 2013; Nielsen et al. 2013a, 2013b) for which we have high-resolution (R ∼ 45,000, ${\rm{\Delta }}v\simeq 6.6$ km s−1) background quasar spectra, in addition to the colors and orientations of the galaxies themselves. With these line of sight data, we can study the detailed kinematics of ${\rm{Mg}}\;{\rm{II}}$ absorbers and examine if the enhanced detection rates and larger equivalent widths for particular orientations are due to velocity spreads, column densities, or some combination of both.

In a previous work, Kacprzak et al. (2012) examined the equivalent widths of ${\rm{Mg}}\;{\rm{II}}$ absorbers with galaxy orientation for a subset of absorber–galaxy pairs from ${\rm{MAG}}{\rm{II}}{\rm{CAT}}$ (Nielsen et al. 2013a, 2013b) by modeling the galaxies and measuring the azimuthal angles at which background quasars probe the CGM. They found that ${\rm{Mg}}\;{\rm{II}}$ absorption prefers to be located along the major and minor axes of blue galaxies, while there is no such preference for absorbers located around red galaxies (a result confirmed with a larger, statistical sample by Lan et al. 2014), or for sightlines in which no absorption is detected (nonabsorbers) regardless of galaxy color. Kacprzak et al. (2012) also found a tendency for ${\rm{Mg}}\;{\rm{II}}$ equivalent widths to be larger along the minor axis than the major axis, possibly an indication of probing more enriched outflowing material along the minor axis.

In this paper, we expand upon our previous work and examine the kinematics of absorbing gas as a function of galaxy color and inclination as well as the azimuthal angle at which the galaxy is probed by using ${\rm{Mg}}\;{\rm{II}}$ absorber pixel-velocity two-point correlation functions (TPCFs, described in detail by Nielsen et al. 2015) and cloud column densities. This paper is organized as follows. Section 2 describes the sample, our methods for analyzing the spectra and data, and how we calculate the TPCFs. Section 3 details how the absorber kinematics (TPCFs) and cloud column densities differ for galaxies of various colors, azimuthal angles, and inclinations. In Section 4 we place our kinematics results in the context of the baryon cycle and discuss the implications. Finally, we summarize and conclude our findings in Section 5. We adopt a ΛCDM cosmology (${H}_{0}=70$ km s−1 Mpc−1, ${{\rm{\Omega }}}_{M}=0.3,$ ${{\rm{\Omega }}}_{{\rm{\Lambda }}}=0.7$) throughout this paper.

2. SAMPLE AND DATA ANALYSIS

2.1. Data

We use a sample of 30 spectroscopically confirmed ($0.3\lt {z}_{{\rm{gal}}}\lt 1.0$) ${\rm{Mg}}\;{\rm{II}}$ absorption-selected galaxies from the ${\rm{Mg}}\;{\rm{II}}$ Absorber–Galaxy Catalog (${\rm{MAG}}{\rm{II}}{\rm{CAT}}$; Nielsen et al. 2013a, 2013b). The galaxies are isolated to the limits of the data, where isolation is defined as having no spectroscopically identified neighbor within a projected distance of 100 kpc or a line of sight velocity of 500 km s−1. Each galaxy has been imaged with the WFPC2/Hubble Space Telescope (HST) in the F702W band and we have a rest-frame BK color for each galaxy that was determined as described in Nielsen et al. (2013b).

All galaxies were modeled using GIM2D (Simard et al. 2002) to obtain quantified morphological parameters, inclinations, and position angles. Full details of the method used and the morphological properties of most galaxies in our sample are presented in Kacprzak et al. (2011). We define inclinations of i = 0° as face-on and i = 90° as edge-on. We convert position angles to an "azimuthal angle," which describes where a background quasar sightline is located with respect to the projected major axis of the galaxy. An azimuthal angle of Φ = 0° is defined as having the quasar line of sight along the projected galaxy major axis and Φ = 90° as having the sightline along the projected minor axis. Figure 1 illustrates our galaxy orientation definitions with azimuthal angle in panel (a) and inclination in panel (b).

Figure 1.

Figure 1. Diagram demonstrating galaxy azimuthal angle and inclination subsamples. Panel (a) presents an on-the-sky view of the azimuthal angle around an inclined galaxy (gray ellipse). An azimuthal angle of Φ = 0° is defined as the projected galaxy major axis, while Φ = 90° is the projected galaxy minor axis. Galaxies that are probed by quasar sightlines at $0^\circ \leqslant {\rm{\Phi }}\lt 45^\circ $ are included in the "major axis" subsample (green shaded region), while those probed at $45^\circ \leqslant {\rm{\Phi }}\leqslant 90^\circ $ are included in the "minor axis" subsample (pink shaded region). Colors in this panel correspond to the colors in Figures 2, 3, and 6. Panel (b) presents a side view of three galaxies with different inclinations whose CGM is probed by a background quasar. We define "face-on" galaxies as those with $0^\circ \leqslant i\lt 57^\circ $ (blue ellipse and arrow) while "edge-on" galaxies are those with $57^\circ \leqslant i\leqslant 90^\circ $ (orange ellipse and arrow, including the gray ellipse). An inclination of i = 57° is the value for which we cut our sample and corresponds to the minimum inclination for an "edge-on" galaxy (gray ellipse). Colors and arrow line widths in this panel correspond to the colors in Figures 2, 4, and 5.

Standard image High-resolution image

The CGM of each galaxy is probed by a nearby (projected on the sky distance of 20 < D < 110 kpc) background quasar for which we have a high-resolution HIRES/Keck or UVES/VLT spectrum. We refer the reader to Churchill (1997), Churchill & Vogt (2001), Churchill et al. (2003), Evans (2011), and Kacprzak et al. (2011) for the spectra, full details of the reduction and analysis, and the detection and Voigt profile (VP) fitting of ${\rm{Mg}}\;{\rm{II}}\lambda \lambda 2796,2803$ absorption. We measure ${\rm{Mg}}\;{\rm{II}}$ equivalent widths, and from the VP fitting we obtain VP component (cloud) column densities, velocities, and Doppler b parameters. Velocity zero points are defined as the median velocity of the apparent optical depth distribution of ${\rm{Mg}}\;{\rm{II}}$ absorption (Churchill 1997).

For our TPCF analysis, we use only the velocities of pixels that reside in regions where absorption is formally detected. These regions are defined as "kinematic subsystems" in Churchill & Vogt (2001) and their velocity bounds are determined by searching the spectra to either side of the subsystem centroids for the point at which the significance in the per pixel equivalent width falls below 1σ. We enforce a sensitivity cut to these regions of ${W}_{r}(2796)\geqslant 0.04\;\mathring{\rm A} $ to account for differences in the quality of our spectra (Nielsen et al. 2015). We are ∼95% complete to this detection sensitivity within ±800 km s−1 for all absorbers in our sample. This sensitivity cut is used to ensure that we can detect absorption uniformly throughout our sample. Thus, the results in this paper apply to ${\rm{Mg}}\;{\rm{II}}$ absorption with ${W}_{r}(2796)\geqslant 0.04\;\mathring{\rm A} ,$ and therefore samples only those regions of the CGM for which the temperature, metallicity, ionization, and line of sight conditions give rise to gas detected with ${\rm{Mg}}\;{\rm{II}}$ for our sample.

We slice the sample into several subsamples based on galaxy rest-frame BK color, azimuthal angle, Φ, and inclination, i. Galaxy color and azimuthal angle cuts were determined by calculating the median values for the sample we present here (similar to the method used in Nielsen et al. 2013a, 2013b). Our inclination cut was defined by the average inclination of galaxies in the universe, i = 57fdg3 (for a derivation, see the appendix of Law et al. 2009). We tabulate the characteristics of each subsample in Table 1 including the subsample names, number of galaxies in each subsample, and the median value(s) by which the sample was cut. The BK color cut we use to separate our sample into "blue" and "red" galaxies is strictly to indicate whether the galaxy is more likely to be star-forming (blue) or passive (red), rather than indicating morphological types such as early-type or late-type galaxies. Figure 1 illustrates our subsample cuts for azimuthal angle (panel (a)) and inclination (panel (b)).

Table 1.  TPCF Velocity Measurements

Galaxy Sample # Gals Cut Cut ${\rm{\Delta }}v(50)$ a ${\rm{\Delta }}v(90)$ a
Figure 2
Face-on 17 $i\lt 57^\circ $ ${69}_{-9}^{+8}$ ${200}_{-21}^{+15}$
Edge-on 13 $i\geqslant 57^\circ $ ${64}_{-6}^{+4}$ ${149}_{-14}^{+10}$
Major Axis 15 ${\rm{\Phi }}\lt 45^\circ $ ${55}_{-13}^{+7}$ ${146}_{-32}^{+13}$
Minor Axis 15 ${\rm{\Phi }}\geqslant 45^\circ $ ${73}_{-7}^{+5}$ ${192}_{-21}^{+15}$
Face-on—Major Axis 10 $i\lt 57^\circ $ ${\rm{\Phi }}\lt 45^\circ $ ${45}_{-6}^{+6}$ ${124}_{-22}^{+20}$
Face-on—Minor Axis 7 $i\lt 57^\circ $ ${\rm{\Phi }}\geqslant 45^\circ $ ${84}_{-13}^{+8}$ ${220}_{-25}^{+13}$
Edge-on—Major Axis 5 $i\geqslant 57^\circ $ ${\rm{\Phi }}\lt 45^\circ $ ${73}_{-20}^{+10}$ ${171}_{-45}^{+23}$
Edge-on—Minor Axis 8 $i\geqslant 57^\circ $ ${\rm{\Phi }}\geqslant 45^\circ $ ${60}_{-5}^{+4}$ ${141}_{-11}^{+8}$
Figure 3
Blue—Major Axis 5 $B-K\lt 1.4$ ${\rm{\Phi }}\lt 45^\circ $ ${61}_{-11}^{+9}$ ${162}_{-23}^{+20}$
Blue—Minor Axis 10 $B-K\lt 1.4$ ${\rm{\Phi }}\geqslant 45^\circ $ ${81}_{-9}^{+8}$ ${209}_{-23}^{+16}$
Red—Major Axis 10 $B-K\geqslant 1.4$ ${\rm{\Phi }}\lt 45^\circ $ ${51}_{-13}^{+10}$ ${135}_{-36}^{+20}$
Red—Minor Axis 5 $B-K\geqslant 1.4$ ${\rm{\Phi }}\geqslant 45^\circ $ ${52}_{-4}^{+1}$ ${123}_{-9}^{+2}$
Figure 4
Blue—Face-on 8 $B-K\lt 1.4$ $i\lt 57^\circ $ ${88}_{-10}^{+8}$ ${227}_{-19}^{+12}$
Blue—Edge-on 7 $B-K\lt 1.4$ $i\geqslant 57^\circ $ ${61}_{-8}^{+4}$ ${143}_{-19}^{+9}$
Red—Face-on 9 $B-K\geqslant 1.4$ $i\lt 57^\circ $ ${40}_{-4}^{+3}$ ${97}_{-9}^{+6}$
Red—Edge-on 6 $B-K\geqslant 1.4$ $i\geqslant 57^\circ $ ${62}_{-14}^{+7}$ ${147}_{-33}^{+15}$

Note.

akm s−1.

Download table as:  ASCIITypeset image

We caution that, in our sample, a weak trend exists between BK and halo mass, $\mathrm{log}({M}_{{\rm{h}}}/{M}_{\odot }),$ (Kendall τ rank correlation test, 2.1σ) where redder galaxies tend to be more massive. We find that, with the exception of eight galaxies, blue galaxies tend to be low mass while red galaxies tend to be high mass (see Nielsen et al. 2015, for details). Therefore, the color dependencies in our results are more accurately color–mass dependencies.

To rule out the possibility that any differences in our results are due to biased distributions in azimuthal angle and inclination, we ran a one-dimensional Kolmogorov–Smirnov (KS) test for both orientation measures. We find that the azimuthal angles and inclinations of galaxies in our sample are consistent with unbiased samples at the 0.6σ and 2.3σ levels, respectively. Additionally, rank-correlation tests between Φ or i and galaxy properties such as rest-frame BK color show no correlations. Therefore, any differences we see in our results are likely not due to underlying sample biases.

2.2. Pixel-velocity Two-point Correlation Functions

Throughout this paper we examine absorber pixel-velocity two-point correlation functions (TPCFs), which are a measure of the internal absorber velocity dispersion. We remind the reader that the absorber–galaxy sample presented here is an absorption-selected sample with an equivalent width detection threshold of ${W}_{r}(2796)\geqslant 0.4\;\mathring{\rm A} $ (see Section 2.1).

To construct the TPCF, we first define a subsample of galaxies, e.g., blue galaxies, and examine the associated absorbers. We obtain the velocities of all pixels that reside only in spectral regions where ${\rm{Mg}}\;{\rm{II}}$ λ2796 absorption is formally detected using the detection methods of Churchill & Vogt (2001), who refer to these regions as "kinematic subsystems." After pooling all of the absorbing pixels from every line of sight in the galaxy subsample together, we then calculate the velocity differences of each possible pair of pixels and take their absolute value to get ${\rm{\Delta }}{v}_{{\rm{pixel}}}.$ We bin up these pixel-velocity separations and normalize each bin by the total number of pixel pairs in the subsample for comparison between different galaxy subsamples. Thus, the TPCF is a probability distribution function. We use a bin size of 10 km s−1, which corresponds to roughly one resolution element of both the HIRES/Keck and UVES/VLT spectrographs (three pixels per resolution element, with a resolution of ∼6.6 km s−1).

To obtain the uncertainties on the TPCF, we conduct a bootstrap analysis for 100 realizations, where the maximum number of realizations allowed for a sample size of five (our smallest sample) with replacement is 126. More realizations than this begin repeating permutations too often. We calculate 1σ standard deviations from the mean of the bootstrap realizations in each bin.

We quantitatively characterize the TPCFs by measuring the velocity separations within which 50% and 90% of the area under the TPCF distribution is located, ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$, respectively. Uncertainties on these values are obtained from the bootstrap analysis and represent 1σ deviations from the mean of the bootstrap realizations. The values for ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ are tabulated in Table 1 for each subsample.

3. RESULTS

3.1. TPCFs: Galaxy Inclinations and Azimuthal Angles

In Figure 2 we present TPCFs for galaxy subsamples cut by galaxy azimuthal angle, Φ, and/or inclination, i, for all galaxy colors, BK. Solid lines represent the TPCF for each subsample while shaded regions indicate the 1σ uncertainties on the TPCFs from the bootstrap analysis described in Section 2.2. Colors for azimuthal angle and inclination subsamples correspond to those in Figures 1(a) and (b), respectively. In the top row we compare face-on and edge-on galaxy subsamples at all azimuthal angles (panel (a)), along the projected minor axis (panel (b)), and along the projected major axis (panel (c)). In the bottom row, we compare galaxies probed along the major and minor axes for galaxies at all inclinations (panel (d)), edge-on galaxy subsamples (panel (e)), and face-on galaxy subsamples (panel (f)). In each panel we present the results of a chi-squared test comparing the binned data (TPCFs) for each plotted subsample pair, including the significance, reduced chi-squared value, ${\chi }_{\nu }^{2},$ and degrees of freedom, ν.

Figure 2.

Figure 2. Pixel-velocity two-point correlation functions examining how the spread in the pixel velocities of absorbers differs when probing galaxies of different inclinations and azimuthal angles (for all galaxy colors, BK). Panels (a)–(c) compare TPCFs of face-on ($i\lt 57^\circ $) and edge-on ($i\geqslant 57^\circ $) galaxies for all Φ, ${\rm{\Phi }}\geqslant 45^\circ $, and ${\rm{\Phi }}\lt 45^\circ $, respectively. Panels (d)–(f) examine how the TPCFs along the major axis (${\rm{\Phi }}\lt 45^\circ $) and minor axis (${\rm{\Phi }}\geqslant 45^\circ $) differ for all i, $i\geqslant 57^\circ $, and $i\lt 57^\circ $, respectively. In each panel the TPCF is represented as a solid line while shaded regions are the 1σ bootstrap uncertainties. We list the significance of a chi-squared test comparing the TPCFs in each panel in addition to the reduced chi-squared value, ${\chi }_{\nu }^{2},$ and the degrees of freedom, ν. We see a broader velocity dispersion for face-on galaxies than for edge-on galaxies when all azimuthal angles are considered (panel (a)). We also see a larger velocity dispersion for galaxies probed along the minor axis than along the major axis when all inclinations are considered (panel (d)). This latter result becomes highly significant for face-on galaxies where the minor axis dispersion is greater than that for the major axis (panel (f)). Conversely, there is no difference in the velocity dispersion with azimuthal angle for edge-on orientations (panel (e)). Thus, the velocity dispersion is the greatest for face-on galaxies probed along the projected minor axis.

Standard image High-resolution image

In Figure 2(a) we examine bivariate trends in the absorber TPCFs for galaxies with "face-on" ($i\lt 57^\circ $) and "edge-on" ($i\geqslant 57^\circ $) inclinations for all colors and probed at all azimuthal angles. Absorbers in face-on galaxies have larger velocity dispersions than those in edge-on galaxies. A chi-squared test yields the result that the null hypothesis that the TPCFs are drawn from the same population can be ruled out at the 8.7σ level. While the chi-squared test result is significant, we find that the ${\rm{\Delta }}v(50)$ measurements are consistent within uncertainties for face-on and edge-on galaxies but the value of ${\rm{\Delta }}v(90)$ is larger for face-on galaxies. Therefore, the difference between these two distributions is due to the large velocity dispersion tail in face-on galaxies.

We also examine absorption associated with galaxies probed along the "major axis" (${\rm{\Phi }}\lt 45^\circ $) and "minor axis" (${\rm{\Phi }}\geqslant 45^\circ $) in Figure 2(d) for all inclinations and colors and find that absorbers located in galaxies probed along the minor axis have larger velocity dispersions than those along the major axis (4.3σ). Both ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ are larger for the minor axis sample than the major axis. For this TPCF pair and the rest in this paper, the measurements of both ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ reflect the chi-squared test results, i.e., where we find an insignificant chi-squared value when comparing galaxy subsamples, we also find values of ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ that are consistent within uncertainties between subsamples.

While we find significant differences in the absorber TPCFs with inclination or azimuthal angle, both orientation measures are important for describing the precise location of an absorber around a galaxy. Therefore, in Figure 2 we present a multivariate analysis comparing the absorber TPCFs for face-on and edge-on galaxies probed along the minor axis (panel (b)) and the major axis (panel (c)). In both panels we find that the velocity structure for absorbers probed in face-on galaxies is significantly different from those in edge-on galaxies. Along the minor (major) axis, the velocity dispersion of absorbers is greater (smaller) for face-on galaxies than edge-on galaxies, with a significance of 14σ (3.1σ) in panel (b) (panel (c)).

We also compare absorber TPCFs for galaxies probed along the major and minor axes for edge-on galaxies (panel (e)) and face-on galaxies (panel (f)). We find that absorbers located in edge-on galaxies have similar velocity dispersions regardless of the azimuthal angle at which the galaxy is probed (panel (e), 0.7σ). On the other hand, we find a highly significant difference (13σ) for face-on galaxies probed along the major and minor axes in Figure 2(f) where the absorber velocity dispersions are much greater along the minor axis than the major axis. In fact, both ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ for the face-on, minor axis subsample are roughly twice as large as the face-on, major axis subsample.

3.2. TPCFs: Galaxy Colors and Orientations

We also study the differences in absorber TPCFs for galaxies of various orientations and galaxy color. In Figure 3, we present subsamples cut by azimuthal angle, Φ, and galaxy rest-frame color, BK, for all galaxy inclinations. We find no significant difference (0.2σ) in absorber velocity dispersions along the major axis in "blue" ($B-K\lt 1.4$) and "red" ($B-K\geqslant 1.4$) galaxies (panel (a)) or for absorbers hosted by red galaxies along the major and minor axes (panel (d)). The only TPCF that is significantly different from the rest of the subsamples is the blue, minor axis subsample. This blue, minor axis subsample has a significantly larger absorber velocity dispersion than either the red, minor axis subsample (panel (b), 16σ) or the blue, major axis subsample (panel (c), 3.6σ). While ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ are nearly identical for the red, major axis subsample, the red, minor axis subsample, and the blue, major axis subsample, the measurements for the blue, minor axis subsample are larger than the rest.

Figure 3.

Figure 3. Pixel-velocity two-point correlation functions for subsamples cut by galaxy rest-frame BK color and azimuthal angle. Lines, shading, and the results of a chi-squared test comparing subsamples are plotted as in Figure 2. For all inclinations, blue and red galaxy TPCFs are compared along the major axis (${\rm{\Phi }}\lt 45^\circ $) and the minor axis (${\rm{\Phi }}\geqslant 45^\circ $) in panels (a) and (b), respectively. Panels (c) and (d) present the TPCFs of blue and red galaxies, respectively, at different azimuthal angles. The velocity dispersions are all statistically consistent (panels (a) and (d)) with the exception of large dispersions for the blue, minor axis subsample.

Standard image High-resolution image

In Figure 4, we examine absorber TPCFs of subsamples sliced by inclination, i, and BK color. Absorbers located around edge-on galaxies have similar velocity dispersions (for all Φ) regardless of galaxy color (0σ, Figure 4(b)). The values of ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ for edge-on blue and red subsamples are nearly identical and consistent within uncertainties. On the other hand, absorbers around face-on galaxies at all azimuthal angles in Figure 4(a) are dramatically different for blue and red galaxies. The absorber velocity dispersions are much greater (29σ) for face-on, blue galaxies than for face-on, red galaxies, which are more highly peaked at low velocity separations. The ${\rm{\Delta }}v(50)$ and ${\rm{\Delta }}v(90)$ measurements are roughly twice as large for the face-on, blue subsample than the face-on, red subsample. The velocity dispersions of absorbers are hosted by both blue and red galaxies in Figures 4(c) and (d) depend on the observed inclination, but differ in their trends; blue (red) galaxies have larger velocity dispersions for face-on (edge-on) inclinations than for edge-on (face-on) inclinations. This is due to the significant differences in face-on inclinations (Figure 4(a)) between blue and red galaxies (larger velocity dispersions for blue galaxies than red galaxies) and the nearly identical TPCFs for edge-on blue and red galaxies in panel (b).

Figure 4.

Figure 4. Pixel-velocity two-point correlation functions for subsamples cut by galaxy rest-frame BK color and inclination. Lines, shading, and the results of a chi-squared test comparing subsamples are plotted as in Figure 2. For all azimuthal angles, we compare the TPCFs of blue and red face-on ($i\lt 57^\circ $) galaxies in panel (a) and edge-on ($i\geqslant 57^\circ $) galaxies in panel (b). Panels (c) and (d) compare the TPCFs of blue and red galaxies, respectively, at different inclinations. We find dramatic differences between the absorber velocity dispersions for blue and red galaxies for face-on orientations (panel (a)) where blue galaxies have much larger velocity dispersions than red galaxies, but we find no differences between blue and red for edge-on orientations (panel (b)).

Standard image High-resolution image

3.3. Cloud Column Densities and Velocities

To obtain further understanding of the CGM material being probed for each orientation presented in the TPCFs, we plot cloud (VP component) column densities and velocities in Figure 5. Clouds plotted in panel (a) are associated with galaxies probed along the minor axis. In panel (b), clouds are associated with galaxies probed along the major axis. Scatter plots show the cloud column densities and velocities from VP fitting, where vertical dashed lines on the scatter plots at ± 100 km s−1 are to help guide the eye between panels. Points are colored by galaxy color, BK, while point types represent galaxy inclination, i. Histograms to the left and right in each panel present the distribution of cloud column densities for galaxies sliced by (left) galaxy color and (right) inclination, where the indicated significance is the result of a KS test between subsamples for the unbinned data. Finally, the histograms on top and bottom show the distribution of pixel velocities (rather than cloud velocities) used to calculate the TPCFs for subsamples cut by galaxy color (bottom) and inclination (top). The result of an F-test3 between the unbinned pixel velocities of each subsample pair is indicated next to the top and bottom histograms. In each case, the pixel velocity F-test results show similar results as the TPCF chi-squared test, i.e., a significant chi-squared test corresponds to a significant F-test, and vice versa.

Figure 5.

Figure 5. Cloud (VP component) column densities and velocities, and pixel velocities for galaxy subsamples sliced by inclination and galaxy color for galaxies probed along (a) the minor axis (${\rm{\Phi }}\geqslant 45^\circ $) and (b) the major axis (${\rm{\Phi }}\lt 45^\circ $). Scatter plots show the modeled cloud velocities and column densities. Blue points represent blue galaxies, red points are red galaxies, circles are face-on galaxies, and diamonds are edge-on galaxies. Vertical dashed lines at ±100 km s−1 are provided to help guide the reader's eye between panels. Histograms on the left and right of each panel show the distribution of cloud column densities for subsamples sliced by (left) galaxy color, BK, and (right) inclination, i. The quoted significance near the column density histograms is the result of a KS test between plotted subsamples. Top and bottom histograms of each panel present the distribution of pixel (not cloud) velocities for subsamples cut by (bottom) galaxy color and (top) inclination. The significances quoted for the top and bottom histograms are for an F-test comparing the variance in the distribution of pixel velocities for each pair of plotted subsamples. All histograms are normalized by the total number of data points in each subsample. In general, we find that clouds with larger velocities relative to the optical depth median velocity have smaller column densities, while those with larger column densities have smaller velocities relative to the optical depth median velocity. We find that the cloud column density distributions are statistically consistent for galaxies probed along the major axis (panel (b)), regardless of whether we compare subsamples sliced by galaxy color or inclination. However, we find statistically larger column densities for blue and face-on galaxies probed along the minor axis (panel (a)).

Standard image High-resolution image

In general, we find that higher velocity material has lower column density values, while higher column density material is only found at a lower velocity. This is consistent with previous works (e.g., Churchill et al. 2003). The highest velocity ($| v| \gt 100\;\mathrm{km}\;{{\rm{s}}}^{-1})$, low column density clouds are mostly associated with blue, face-on galaxies probed along the minor axis. There is also a population of $| v| \sim 100$ km s−1 clouds that are associated with edge-on galaxies probed along the minor axis, though we find a few for major axis galaxies as well. All of these clouds at high velocity are found in absorbers that have higher column density material located near the velocity zero point.

For minor axis galaxies in Figure 5(a), we find that blue galaxies (all inclinations) have significantly larger column densities than red galaxies (3.2σ) and edge-on galaxies (all colors) have significantly larger column densities than face-on galaxies (3.4σ). For major axis galaxies in Figure 5(b), we find only a suggestion of larger column densities for edge-on galaxies than face-on galaxies (2.9σ, all colors), though clouds at $\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\gt 14$ are mainly hosted by edge-on galaxies. Conversely, we find no difference in the column density distribution between blue and red galaxies (1.9σ, all inclinations) along the major axis.

In Figure 6 we compare the cloud column densities and velocities, and the pixel velocities for (a) face-on galaxies and (b) edge-on galaxies, slicing the samples by galaxy color and azimuthal angle. Similar to Figure 5, point colors and the left and bottom histograms represent blue and red galaxy subsamples, while point types and top and right histograms represent major axis and minor axis subsamples. The indicated significances correspond to the same statistical tests as in Figure 5. Again, the pixel velocity F-test results in Figure 6 follow the TPCF chi-squared test results, with the exception of the test comparing edge-on galaxies probed along the major versus minor axes (panel (b)). In this case the F-test result is significant, where edge-on, major axis galaxies have a slightly larger dispersion in their pixel velocities than for edge-on, minor axis galaxies. Due to the large bootstrap uncertainties on the major axis TPCF, the chi-squared test on the corresponding TPCFs is not significant (see Figure 2(e)).

Figure 6.

Figure 6. Cloud column densities and velocities, and pixel velocities for (a) face-on galaxies ($i\lt 57^\circ $) and (b) edge-on galaxies ($i\geqslant 57^\circ $) split by galaxy color and azimuthal angle. Blue galaxies are represented as blue points, red galaxies as red points, galaxies probed along the major axis as triangles, and galaxies probed along the minor axis as stars. Vertical dashed lines are plotted to help guide the reader's eye between panels. Histograms in each panel represent the distributions of cloud column densities split by (left) galaxy color, BK, (right) azimuthal angle, Φ, and pixel velocities (rather than cloud velocities) split by (bottom) galaxy color, and (top) azimuthal angle. The quoted significances near the histograms report the results of a KS test between plotted subsamples for the cloud column densities (left and right) and the results of an F-test comparing the variance in the pixel velocity distributions between plotted subsamples for the top and bottom histograms. We find that the cloud column densities for absorbers around face-on galaxies rarely exceed $\mathrm{log}N({\rm{Mg}}\;{\rm{II}})=14$ (panel (a)), in contrast to edge-on galaxies (panel (b)). For both face-on and edge-on galaxies, we see no difference in the column densities with regard to the azimuthal angle (right histograms). Lastly, we find that the cloud column densities in absorbers around blue galaxies with edge-on inclinations are greater than those for absorbers around blue, face-on galaxies (4.7σ, not plotted) or red, edge-on galaxies (4.0σ, panel (b)).

Standard image High-resolution image

We find that absorbers around face-on galaxies consist of clouds with column densities that rarely exceed $\mathrm{log}N({\rm{Mg}}\;{\rm{II}})=14.$ This is in contrast to edge-on galaxies, which are more likely to have these higher column density clouds. For face-on galaxies, the column densities of blue and red galaxies are statistically consistent (0.3σ). Comparing blue galaxies with face-on and edge-on inclinations (for all azimuthal angles), we find that the column densities of the edge-on subsample are larger than those in the face-on subsample (4.7σ, not plotted). We also find that the column densities of the blue, edge-on subsample are larger than the red, edge-on sample (4.0σ). There are no such differences in the column densities between galaxies probed along the major and minor axes for face-on (0.7σ), edge-on (0.9σ), blue (1.4σ, not plotted) or red (1.5σ, not plotted) galaxies.

4. DISCUSSION

The results in the previous section demonstrate that absorber velocity dispersions and cloud column densities of the gas traced by ${\rm{Mg}}\;{\rm{II}}$ absorption detected with ${W}_{r}(2796)\geqslant 0.04\;\mathring{\rm A} $ depend strongly on where the gas is located around isolated galaxies. In particular, we find that the absorbers with the largest velocity spreads are associated with blue galaxies that have "face-on" inclinations and are probed along the galaxy projected minor axis. For "edge-on" galaxies, the velocity structure is similar regardless of galaxy color or azimuthal angle. Thus, the absorber velocity spreads (TPCFs) depend mostly on if the galaxy is "star-forming" (using color as a proxy for the star formation rate), or if the galaxy is "face-on." The cloud column densities depend mostly on galaxy inclination or color, where face-on or red galaxies host absorbers with smaller cloud column densities than those hosted by edge-on or blue galaxies.

Mechanisms in the baryon cycle that could give rise to these differences in velocity structure and column densities of the ${\rm{Mg}}\;{\rm{II}}$-absorbing gas include (1) merging satellite galaxies whose gas is either being ejected due to star formation or being stripped by tidal forces, (2) IGM accretion, recycled accretion, and/or rotating material merging onto the galaxy disk, and (3) outflowing gas in star-forming galaxies. In this section, we discuss which mechanisms may most likely contribute to the differences in the observed velocity dispersions and cloud column densities as a function of host galaxy color (as a proxy for star formation) and orientation.

4.1. Merging Satellite Galaxies

Regardless of our results with orientation, previous works have repeatedly shown that the probability of probing a satellite galaxy with quasar absorption lines is low. Martin et al. (2012) found only one satellite galaxy contributing to ${\rm{Mg}}\;{\rm{II}}$ absorption out of a full sample of over 200 galaxies, corresponding to a probability of finding absorption due to a satellite galaxy of <1%. Tumlinson et al. (2013) used simulations to estimate the cross-section of ${\rm{H}}\;{\rm{I}}$ gas bound to satellites and assumed that the internal velocity dispersion of gas in the satellites cannot exceed the maximum circular velocity of the satellite. They found that the mean number of satellites per sightline is well below the number of absorption components per sightline at a given velocity for their sample of ${\rm{H}}\;{\rm{I}}$ absorbers. For ${\rm{Mg}}\;{\rm{II}}$ absorption, Gauthier et al. (2010) examined covering fractions for absorbers around luminous red galaxies (which are comparable to our reddest galaxies) and compared them to the estimated maximum cross-section of satellites in host halos. They found that this maximum cross-section is much lower than the absorber covering fractions at all impact parameters, and therefore satellites are not sufficient for explaining the presence of the majority of absorption around galaxies.

If we assume that the cross-section of satellites is great enough to explain the presence of the majority of ${\rm{Mg}}\;{\rm{II}}$ absorption, we can investigate the preferential location of satellite galaxies around a larger, central galaxy to determine if satellites could drive our results. Many previous works have examined satellite distributions and the consensus appears to be that for red galaxies, satellites within the virial radius are preferentially aligned with the central galaxy's major axis, while there is no preferred alignment for satellites around blue galaxies (with the exception of galaxies in the Local Group; Yang et al. 2006, and references therein). This has been confirmed in recent simulations by Dong et al. (2014). Our color cut defining blue and red galaxies is slightly bluer than that used by Yang et al. (2006; our BK = 1.399 cut corresponds to gr = 0.62, compared to their cut of gr = 0.83). Our bluer color cut may result in more blue galaxies in our "red" sample than in the Yang et al. (2006) red sample. Since the blue galaxies in Yang et al. have been found to have an isotropic distribution of satellites, we may have a somewhat more isotropic distribution of satellites around red galaxies for our sample.

We cannot rule out that we may be observing some amount of ${\rm{Mg}}\;{\rm{II}}$ absorption due to satellites. However, even though satellites align with the major axis of red galaxies, we find no difference between the TPCFs or column densities of our red galaxy samples along the major and minor axes. The isotropic distribution of satellites around blue galaxies is not consistent with the highly significant differences between the TPCFs of our blue galaxy samples along the major and minor axes. If satellites are the dominant source of ${\rm{Mg}}\;{\rm{II}}$ absorption, we would expect the TPCFs to show differences along the major and minor axes in red galaxies (which could be mitigated by our bluer color cut), but show no differences in blue galaxies (which we do find).

We do not have velocity information on the host galaxy–satellite galaxy samples published in Yang et al. (2006). However, the velocity distribution of the satellites around galaxies is unlikely to be the cause of the different absorber kinematic distributions for several reasons. We are measuring the internal velocity dispersions of the absorbers themselves, not the material with respect to the galaxy systemic velocity. We assume our sightlines are unlikely to pass through multiple satellites with a single sightline since the probability of passing through a single galaxy is low (<1%, from above), and probing multiple satellites at once is even lower. Therefore the internal gas velocity distributions of the satellite galaxies would have to show differences for blue central galaxies (to match the TPCFs), but not for red central galaxies. However, if ${\rm{Mg}}\;{\rm{II}}$ absorption is due to tidal stripping of or outflows from the satellites, we would expect the satellites around red galaxies (which tend to be located in more massive halos in our sample) to be more disturbed and therefore have larger absorber velocity dispersions than if they were located around blue galaxies, which tend to be less massive in our sample. This is the reverse of our findings.

The galaxies in our sample were drawn from ${\rm{MAG}}{\rm{II}}{\rm{CAT}}$ (Nielsen et al. 2013b), which contains only galaxies that are isolated to the limits of the data, where no spectroscopically identified neighbor was found within $D\lt 100\;{\rm{kpc}}$ or a line of sight velocity of 500 km s−1. If satellites are present in our sample, they are undetected within both velocity space and impact parameter. Additionally, all of the galaxies were imaged with HST and satellite galaxies were not detected in any of the images to the limits of our data. Therefore, if we are probing satellite galaxies, they are below our detection limits (i.e., below roughly $0.1{L}_{B}^{*},$ corresponding to a circular velocity of ∼80 km s−1 for these satellites; Churchill et al. 2013; Nielsen et al. 2013b).

4.2. Accretion and Rotation

Previous work examining the low-ionization CGM has found accreting and/or rotating material for edge-on galaxies probed along the projected major axis in the form of absorption that is shifted to one side of the galaxy systemic velocity (e.g., Steidel et al. 2002; Kacprzak et al. 2010; Rubin et al. 2012; Bouché et al. 2013). This is also commonly observed in simulations (e.g., Stewart et al. 2011). Since we do not know whether our absorbers are blueshifted or redshifted with respect to the galaxy as can be determined through the down-the-barrel approach used by Rubin et al. (2012), we cannot say for certain whether our absorbers are accreting or rotating. However, we have examined the properties of this material using absorber velocity dispersions and cloud column densities for orientations in which accreting and rotating material is expected to be found.

Edge-on galaxies probed along the major axis are most likely to exhibit accreting or rotating material and are the least likely orientation to be contaminated with outflowing material. The line of sight velocity for rotation is greatest in edge-on galaxies since the rotation is in the plane of the galaxy. In our sample, this orientation is dominated by absorbers that contain clouds with larger column densities ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\gt 13$), most of which have velocities $| v| \lt 100$ km s−1. This is most easily seen as the triangle points in Figure 5(b). We might expect to only see blue galaxies with absorption in this orientation due to the red galaxies having their star formation quenched from a lack of a gas reservoir, but the subsample is populated by both blue and red galaxies.

We could be seeing recycled accretion in galaxies of all colors from past outflows (e.g., Oppenheimer et al. 2010). However, accretion onto galaxies provides fuel for star formation, which we do not infer for our red galaxies. It is possible that this material is rotating around the galaxy and not accreting onto the galaxy, especially for the red galaxies, though we cannot tell if this is the case for the data presented here as our velocities are not shifted with respect to the galaxy systemic velocity. Some mechanism may be suppressing star formation in the red galaxies, but this mechanism is not present in blue galaxies. It is also possible that some of the red galaxies we have observed are in fact dusty, star-forming galaxies that are currently accreting material, especially since our BK color cut does include Sbc SED types in with the "red" galaxy subsample (for details concerning our BK color calculations, see Nielsen et al. 2013b). This is less likely because we would have observed signatures of outflows in both blue and red galaxies, unless dusty, star-forming galaxies are rare in our sample. We do not yet have the information needed to calculate the star formation rates to determine if we have dusty, star-forming galaxies in our sample.

It is interesting that we find no significant difference in the TPCFs of blue and red galaxies probed along the major axis (for all i) or in edge-on inclinations (for all Φ) in Figures 3(a) and 4(b), respectively. This suggests that the velocity structure of absorbers along the major axis or for edge-on galaxies does not depend on the star formation rate of the host galaxy, and likely depends only on mass. This fits in with the accretion or rotating gas picture, which, unlike outflows, should not depend on the star formation rate unless the star formation rate is so great as to prevent accretion and/or remove the gas reservoir. Also interesting is that, for all BK, the absorber velocity dispersions depend on inclination for the major axis sample in Figure 2(c). Along the major axis, the absorber velocity dispersion is greater for edge-on galaxies than for face-on galaxies, and this may be due to rotating gas whose line of sight velocity is maximized in edge-on inclinations, while the vertical velocity dispersions in the disks of galaxies are small. The column densities may also be greater for edge-on galaxies than for face-on galaxies along the major axis (Figure 5(b)). This could indicate that the accreting or rotating material for the edge-on galaxies is more coherent, i.e., the path lengths or amount of material are larger than for face-on galaxies.

4.3. Outflows

Absorption line studies have frequently found evidence for outflows as blueshifted absorption in face-on galaxies (e.g., Rubin et al. 2014) or as enhanced equivalent widths along the projected minor axis for edge-on galaxies (e.g., Bordoloi et al. 2011; Kacprzak et al. 2012). Therefore we also can associate the kinematics results for our blue, face-on, and minor axis subsamples with outflows. Again, we cannot say for certain if our absorbers are actually entrained in outflowing material so we instead focus on possible absorption properties of the absorbers for orientations in which outflows are most likely to occur.

For blue, edge-on galaxies probed along the minor axis (blue triangles in Figure 5(a)), we found a group of lower column density ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\lt 13$), higher velocity ($| v| \sim 100$ km s−1) material. This is similar to Fox et al. (2015), who used a background quasar whose line of sight passes through the biconical Fermi Bubbles at the Milky Way Galactic center, where the Milky Way in this case is similar to our edge-on, minor axis subsample. They reported a complex absorption profile with two high-velocity metal absorption components at ${v}_{{\rm{LSR}}}\sim \pm 250$ km s−1 with $\mathrm{log}N({\rm{Si}}\;{\rm{II}})\sim 13.2.$ They associate these components with cool gas that has been entrained in the near and far sides of an outflow from the Galactic center. In our sample, comparable components have similar column densities, but a lower velocity with respect to the bulk of the absorption, which may be due to the fact that we are probing outflows at much larger impact parameters than Fox et al. (2015); D > 20 kpc compared to $D=2.3\;{\rm{kpc}}.$ This suggests that the low column density clouds that we find in our edge-on, minor axis galaxies may also be associated with fragmented cool gas entrained in outflows.

In our subsample consisting of blue, face-on galaxies probed along the minor axis (blue circles in Figure 5(a)), the absorbers are dominated by low column density material ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\lt 14,$ where most of the clouds are $\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\leqslant 13$), with high velocity components at $| v| \gt 100$ km s−1. This material is also similar to the absorber found by Fox et al. (2015). We suggest that the velocities are larger for our face-on galaxies than for our edge-on galaxies along the minor axis because we are observing down into the outflow (i.e., close to the down-the-barrel approach) as opposed to across the outflow. We expect the line of sight velocity dispersions to be larger for the face-on sample, as the biconical outflows are pointed toward (or away from) the observer, and this is observed in our face-on, blue galaxy subsample TPCFs.

The larger column density ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\gt 13$), lower velocity ($| v| \lesssim 50$ km s−1) clouds for blue, edge-on galaxies probed along the minor axis (seen in Figure 5(a) as blue triangles), suggest that we are observing more material, a larger path length, a larger metallicity, or ionization conditions that are more conducive to ${\rm{Mg}}\;{\rm{II}}$ absorption than for the face-on, minor axis galaxies (blue circles) described in the previous paragraph. Since the material probed in both of these orientations is most likely outflows, we do not expect the metallicity or the ionization conditions to differ for these two subsamples, and therefore should not be the cause of the enhanced column densities for the edge-on sample compared to the face-on sample. However, Veilleux et al. (2005) listed large variations in the metallicity and ionization levels between clouds as a possible signature of outflows in quasar absorption lines. If this were the case here, we should also see the large column density clouds ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\gt 13$) for face-on galaxies, but they are nonexistent in our sample. Differences in the column densities with inclination therefore hint at probing larger amounts of gas or larger path lengths for the lowest velocity clouds associated with outflows in edge-on galaxies (compared to face-on galaxies), an effect that may be due to the geometry of the gas flow.

In contrast to the blue galaxies, we find very small absorber velocity dispersions for red galaxies in face-on, minor axis orientations (red circles in Figure 5(a)). We suggest that this is due to a lack of current outflows. This is supported by the fact that redder galaxies typically have lower star formation rates and do not typically drive outflows. We caution that this subsample is small, which may indicate that outflows in red galaxies for this orientation are rare rather than nonexistent in general.

5. SUMMARY AND CONCLUSIONS

We examined the dependence of gas kinematics and cloud column densities on galaxy rest-frame BK color, azimuthal angle, and inclination for a subset of 30 isolated ${\rm{MAG}}{\rm{II}}{\rm{CAT}}$ (Nielsen et al. 2013b) galaxies. Each galaxy was imaged with HST and was modeled with GIM2D to obtain galaxy orientations. In each case, we have a high-resolution optical spectrum of an associated background quasar within $20\lt D\lt 110\;{\rm{kpc}}$ for detailed kinematic analysis. To characterize the absorber gas kinematics, we examined only those absorbers detected with ${W}_{r}(2796)\geqslant 0.04\;\mathring{\rm A} $ and calculated pixel-velocity TPCFs for galaxy color, azimuthal angle, and inclination subsamples.

Our findings include the following:

  • 1.  
    Absorption TPCFs with the longest high velocity separation tails (the largest velocity dispersions) are associated with blue galaxies, face-on galaxies, and galaxies probed along the minor axis. Conversely, the narrowest TPCF is associated with red, face-on galaxies. The velocity structure in the absorbers hosted by edge-on galaxies is similar regardless of galaxy color or azimuthal angle, and we find a lack of differences for red galaxies with azimuthal angle as well as for blue and red galaxies probed along the projected major axis.
  • 2.  
    Examining cloud column densities and velocities, we find, in general, that the highest velocity clouds with respect to the optical depth weighted mean of absorption (${z}_{{\rm{abs}}}$) have small column densities, while the largest column density clouds are located at small velocities. When slicing our sample by galaxy color and inclination, we find that the column densities are larger for edge-on galaxies and blue galaxies than for face-on galaxies and red galaxies, respectively. The column densities show no dependence on azimuthal angle, however.
  • 3.  
    We find large absorber velocity dispersions and large column density clouds at low velocity for edge-on galaxies probed along the major axis, an orientation that is associated with gas accreting onto and/or rotating around galaxies. The velocity structure of the absorbers for this orientation does not depend on galaxy color and, by proxy, star formation rate. The large absorber velocity dispersions for this orientation may be due to rotating gas whose line of sight velocity is maximized due to the edge-on inclination, especially compared to the smaller velocity dispersions for face-on galaxies along the major axis. The large column density clouds may indicate that accreting and/or rotating material is fairly coherent, where the path lengths or amount of gas probed is larger than for other subsamples.
  • 4.  
    We associate the largest absorber velocity dispersions and smaller column density clouds at higher velocity for blue, face-on galaxies probed along the minor axis and for blue, edge-on galaxies along the minor axis with bipolar outflows, which are commonly observed in these orientations. The behavior of this material and the fact that it is located in regions that are associated with outflows may be an indication that the material probed by ${\rm{Mg}}\;{\rm{II}}$ absorption is fragmented material entrained in outflows. These absorbers are similar to the absorption associated with the Fermi Bubbles in the Milky Way Galaxy (Fox et al. 2015). Larger column density clouds at low velocity are present for the blue, edge-on galaxies probed along the minor axis, and since comparable components are lacking for the face-on sample, we suggest that the path lengths of these structures are larger, possibly a geometry effect. For the smallest column density clouds ($\mathrm{log}N({\rm{Mg}}\;{\rm{II}})\lesssim 13$), the cloud velocities (and absorber velocity dispersions) are greater for the face-on galaxies than edge-on galaxies. We attribute this to looking down into the outflow for face-on galaxies, and across the outflow for edge-on galaxies. Conversely, the very small absorber velocity dispersions for red galaxies at these orientations suggest that outflows are not active, a result of star formation having possibly been quenched.
  • 5.  
    Merging satellite galaxies may contribute to the observed ${\rm{Mg}}\;{\rm{II}}$ absorption but are not expected to dominate the absorption due to their cross-section being too low. If they did dominate, then we would expect differences in the TPCFs to correspond to the preferred locations of the satellites around galaxies (isotropic distribution for blue galaxies and along the major axis for red galaxies), but this is not the case. We also do not expect the velocities of the satellites to have much of an effect on our TPCFs because it is unlikely that our sightlines are hitting multiple satellites at once, as the probability of detecting a single satellite is $\lesssim 1\%$ (Martin et al. 2012). Since red galaxies tend to be more massive than blue, we would expect larger internal gas velocity dispersions for satellites hosted by red central galaxies due to tidal stripping or outflows from the satellites themselves, but we do not find this to be the case.

This work shows that the velocity structure and cloud column densities of ${\rm{Mg}}\;{\rm{II}}$ absorbers depend on a combination of galaxy orientations and colors (and likely other galaxy properties not included in this work), not just one or two properties at a time as has been commonly examined in the past. For example, when observing edge-on galaxies, absorption could be probing accretion along the major axis, or outflows along the minor axis. In this case, the outflows will dominate the absorption signature and the information for accretion will likely be lost, which may contribute to the low covering fractions found for accreting material.

We also find that while the equivalent width of an absorption profile provides a useful diagnostic for generally determining what type of material is being studied, considering both the velocity structure (in terms of the velocity spread or dispersion) and cloud column densities of absorbers with galaxy orientation is important for probing the details of gas flows. This is evident in the fact that absorbers probing outflows in face-on galaxies may have large equivalent widths because of the large velocity spread in the absorbers, whereas outflows in edge-on, minor axis galaxies may have large equivalent widths mainly due to large column density clouds. Modeling the ionization conditions and metallicities in the clouds would also be beneficial, as large cloud–cloud variations in outflows have been suggested by Veilleux et al. (2005), but this is beyond the scope of the present work.

Future work might include examining mock quasar absorption line spectra in simulations for a variety of galaxy colors (or star formation rate and mass) and orientations to better understand the origin of the absorber velocity structure and cloud column densities. By applying our analysis to simulations, we not only obtain a more physical explanation for what we observe, but we also help constrain the CGM in simulations. Observationally, a larger sample of galaxies with HST images and high-resolution spectra would allow us to slice the sample by galaxy color, azimuthal angle, and inclination to better explore the multivariate dependence of absorber properties on the host galaxy and its baryon cycle. Our current sample is only large enough to slice by two galaxy properties in the TPCFs, though we have presented all three properties in Figures 5 and 6.

It would also be useful to obtain estimates of the absorber metallicities to further understand whether we are probing cold-mode accretion, recycled accretion, satellite material, and/or outflows. Metallicities are important to help distinguish between accretion and outflows, where accreting material is expected to have lower metallicities than outflowing material (see Lehner et al. 2013, and references therein). Comparable metallicities across all galaxy orientations might indicate that ${\rm{Mg}}\;{\rm{II}}$ absorption is primarily associated with outflows and recycled accretion, not cold-mode accretion. Lastly, determining the star formation rates of the galaxies in our sample is important for identifying any possible dusty, star forming galaxies which may be contaminating the results of our red galaxy samples. The galaxy star formation rates provide a more physical indicator of ongoing star formation than galaxy color.

We thank C. Steidel for providing a reduced HIRES/Keck quasar spectrum. This material is based upon work supported by the National Science Foundation under grant No. 1210200 (NSF East Asia and Pacific Summer Institutes). N.M.N. was also partially supported through an NMSGC Graduate Fellowship and a Graduate Research Enhancement Grant (GREG) sponsored by the Office of the Vice President for Research at New Mexico State University. G.G.K. acknowledges the support of the Australian Research Council through the award of a Future Fellowship (FT140100933). M.T.M. thanks the Australian Research Council for Discovery Project grant DP130100568, which supported this work.

Footnotes

  • The F-test is designed to measure whether two sample distributions have similar or different dispersions; it is sensitive to the tails of the distributions.

Please wait… references are loading.
10.1088/0004-637X/812/1/83