arxiv: v2 [astro-ph] 2 Jan 2008

Transcription

1 DRAFT VERSION OCTOBER 25, 2018 Preprint typeset using LATEX style emulateapj v. 4/12/04 ON THE HALO OCCUPATION OF DARK BARYONS JEREMY L. TINKER& HSIAO-WEN CHEN Kavli Institute for Cosmological Pysics & Department of Astronomy and Astropysics, University of Cicago Draft version October 25, 2018 arxiv: v2 [astro-p] 2 Jan 2008 ABSTRACT We introduce a new tecnique tat adopts te alo occupation framework for understanding te origin of QSO absorption-line systems. Our initial study focuses specifically on Mg II absorbers. We construct a model of te gaseous content in wic te absorption equivalent widt W r is determined by te te amount of cold gas, in te form of discrete clouds, along a sigtline troug a alo. Te two quantities tat we specify per alo in te model are (1) te mean absorption strengt per unit surface mass density A W (M), and (2) te mean covering factor κ g (M) of te gaseous clouds. Tese parameters determine te conditional probability distribution of W r as a function of alo mass, P(W r M). Two empirical measurements are applied to constrain te model: (i) te absorber frequency distribution function and (ii) te W r -dependent clustering amplitude. We find tat te data demand a rapid transition in te gas content of alos at M, below wic alos contain predominantly cold gas and beyond wic gas becomes predominantly ot. In order to reproduce te observed overall strong clustering of te absorbers and te anti-correlation between W r and alo mass M, rougly 5% of gas in alos up to M is required to be cold. Te gas covering factor is near unity over a wide range of alo mass, supporting tat Mg II systems probe an unbiased sample of typical galaxies. We discuss te implications of our study in te contexts of mass assembly of distant galaxies and te origin of QSO absorption line systems. Subject eadings: Cosmology: teory dark matter alos quasars: absorption lines galaxies: evolution 1. INTRODUCTION Te forest of Lyα absorption line systems observed in te spectra of background QSOs is a sensitive probe of intervening baryonic matter tat is oterwise invisible (see Rauc 1998 for a compreensive review). Te gaseous clouds uncovered by tese absorption features span a wide range of ionization state and neutral ydrogen column density over N(H I) = cm 2 and are understood to originate in a range of environments in te dark matter distribution, from moderate overdensities to fully collapsed structures (e.g. Davé et al. 1999). Strong Lyα absorbers wit N(H I) cm 2 are tougt to be te ig-redsift analogue to present-day galaxies (Wolfe et al. 2005). Studies of tese ig-density neutral gas clouds sow tat tey all contain low ions suc as Si +, Fe +, and Mg + (e.g. Lu & Wolfe 1994; Procaska et al. 2003). Conversely, studies indicate tat Mg II absorbers of rest-frame absorption equivalent widt W r (2796) 0.3 Å arise in gaseous clouds of N(H I) = cm 2 (Curcill et al. 2000; Rao et al. 2006). Over te redsift range z = , te Mg II λλ2796,2803 doublet transitions are sifted to te optical spectral range λ obs = ,000 Å and terefore provide a convenient optical probe of dark baryons in extended gaseous alos of typical galaxies along random lines of sigt. Extensive studies of te pysical properties of Mg II absorbers ave been carried out by various autors. Early studies tat compared abundance ratios between different ions associated wit Mg II absorbers sowed tat tese absorbers arise in poto-ionized gas of temperature T 10 4 K (e.g. Bergeron & Stasinska 1986; Hamann 1997). In addition, ig-resolution spectra of strong absorbers also sow tat tese systems are multi-component, wit W r rougly proportional to te number of components in te system (Petitjean & Bergeron 1990; Curcill & Vogt 2001; Procter et al. 2006). Finally, comparisons of galaxies and Mg II absorbers along common lines of sigt also indicate tat Mg II absorbers at z = 0.65 are associated wit luminous galaxies (Bergeron 1986; Steidel et al. 1994) and tat field galaxies possess extended gaseous alos to projected distances of r 50 1 kpc (Lanzetta & Bowen 1990; Steidel 1993). Extended alos of cool gas around galaxies were predicted by Spitzer (1956) and discussed in more detail in Bacall & Spitzer (1969). Teoretical models tat explain te origin of extended gaseous alos around galaxies include (i) te stripping of gas from te accretion of small, gasric satellites into a larger system (Wang 1993); (ii) cold gas gravitationally bound to small-scale substructure witin te dark matter alo (e.g., Sternberg et al. 2002); and (iii) a two-pase medium, cold clouds condense out of a otter alo troug termal instability (Mo & Miralda-Escude 1996; Maller & Bullock 2004; Celouce et al. 2007). Tese models imply tat W r reflects te potential well of te dark matter alo more massive systems sustain bigger gaseous alos wit a iger velocity dispersion. However, recent studies ave yielded conflicting results tat callenge tis classical picture. Wile some autors present observations tat suggest strong Mg II absorbers originating primarily in galactic superwinds from low-mass galaxies (Procter et al. 2006; Boucé et al. 2006; Murpy et al. 2007), oters sow tat Mg II absorbing galaxies resemble typical field galaxies (Zibetti et al. 2005; Kacprzak et al. 2007; Nestor et al. 2007) and te covering fraction of Mg + ions is rougly 50% (Tripp & Bowen 2005). A particularly puzzling result is te anti-correlation between W r (2796) (ereafter W r ) and clustering amplitude by Boucé et al. (2006) (ereafter B06). Tese autors measured te cross-correlation function of Mg II absorbers and luminous red galaxies (LRGs) in data release tree of te Sloan Digital Sky Survey (SDSS DR3; Abazajian 2005). From te relative bias wit respect to LRGs, B06 inferred te mean alo mass of absorbers as a

2 2 function of W r, finding tat strong absorbers of W r = Å on average arise in less massive dark matter alos of logm = 11.11, wile absorbers of W r = Å on average arise in dark matter alos of logm = Te autors consider tis anti-correlation in strong favor of a superwind origin for W r > 1 Å absorbers; te large equivalent widt in low-mass systems is attributed to te outflow velocity of te systems undergoing starburst episodes. Te cross-correlation function of absorbers and galaxies on large scales ( 1 1 Mpc) provides a quantitative caracterization of te origin of tese absorbers. Te clustering of absorbers is a consequence of te alos in wic tey are found. Low-mass alos form in a wide range of environments wile ig-mass alos can only form in te most dense regions of te dark matter distribution, tus ig-mass alos are igly clustered. Te clustering of alos depends only on simple gravitational pysics, and te bias of dark matter alos as a function of mass is a well-establised relation in te standard CDM cosmology. Tus it is straigtforward to infer a mean alo mass from large-scale clustering, as presented in B06. However, wen interpreting te clustering of gaseous clouds based on te presence of a specific ion, it is necessary to consider te state of te gas in te dark matter alos. For example, if te observed Mg + ions originate in potoionized gas of temperature T 10 4 K, ten te observed anti-correlation between W r and M may be produced wen te alo gas in massive dark matter alos becomes too ot for abundant Mg + to survive. To investigate ow te content of warm/cold gas depends on alo mass, we introduce a new tecnique tat adopts te alo occupation framework for studying te origin of Mg II absorbers. Tis same tecnique can also be applied for studying te nature of QSO absorptionline systems. Our alo occupation approac is similar to te Halo Occupation Distribution (HOD) metod tat was developed to establis an empirical mapping between galaxies and dark matter alos (e.g., Seljak 2000; Scoccimarro et al. 2001; Berlind & Weinberg 2002; see Zeng & Weinberg 2007 and references terein). Tis mapping is purely statistical in nature; it interprets galaxy bias based on te probability tat a alo of mass M contains N galaxies of a given class, P(N M). Te probability, P(N M), is constrained by measurements of te space density of te class of objects and teir two-point auto-correlation function. HOD analysis of galaxy clustering data can provide insigt into galaxy formation and evolution (van den Bosc et al. 2003b; Cooray 2006; Zeng et al. 2007; Wite et al. 2007; Tinker et al. 2007) and to constrain cosmological parameters (van den Bosc et al. 2003a, 2007; Tinker et al. 2005; Tinker 2007; Zeng & Weinberg 2007). In tis paper, we extend tis framework to caracterize te cold gas detected troug QSO absorption lines and focus specifically on Mg II absorbers. To study te dark matter alo population of Mg II absorbers, we adopt a conditional probability distribution of equivalent widts P(W r M) to describe te alo occupation of cold gas. Our tecnique is fundamentally te same as te conditional luminosity function for interpreting galaxy bias (Yang et al. 2003; van den Bosc et al. 2007). In our approac, we coose to parameterize P(W r M) te distribution of W r witin a alo of given mass according to expectations of a uniform gaseous alo, in wic te equivalent widt is proportional to te density-weigted pat-lengt along a sigt-line troug te alo. Te parameters of P(W r M) are left free and constrained by empirical data. Te available observations are (1) te frequency distribution function d 2 N/(dW r dz), te number density of absorbers per unit absorption equivalent widt interval per unit redsift pat, and (2) te large-scale bias relative to LRGs. Te alo occupation analysis allows us to go beyond simply estimating te mean mass of dark matter alos tat ost absorbers and constrain te distribution of tese alos as a function of W r. We will consider two scenarios for P(W r M). In te first scenario, te fraction of cold gas relative to dark matter mass is a smootly varying function of alo mass. Tis scenerio is motivated by te canonical model for alo gas in wic accreting gas is always sock-eated to te virial temperature wen it crosses te virial radius (e.g., Wite & Frenk 1991). Te gas ten cools from te inside out, creating a cooling radius tat varies wit te properties of te alo. We will demonstrate tat tis model cannot simultaneously fit te frequency and bias data; models of tis type tat reproduce te frequency distribution predict a positive correlation between W r and M because more massive alos are expected to contain more absorbing gas. Taking into account te fact tat te Mg II transitions probe primarily warm, poto-ionized gas, we include in te second model a transitional mass scale, at wic sock eated gas progressively takes over iger mass alos and te absorption efficiency is significantly reduced. In te transition region between entirely cold-mode and ot-mode alos, eating occurs as an inside-out process. Namely, a igly ionized core expands until it envelopes te entire gas alo (e.g., Dekel & Birnboim 2006). We allow te mass scale and widt of tis transition to be free parameters. In addition, we allow some fraction of cold gas in te ot alos tat is also constrained by empirical data. We demonstrate tat te second scenario reproduces bot te observed d 2 N/(dW r dz) and bias data well. We obtain a best-fit transitional mass scale of M and a cold gas fraction (as represented by te Mg + ions) of<10% in massive alos. Given te agreement between observations and our alo occupation model, we argue tat it is not necessary to invoke additional, complicated star formation feedback for explaining te inverse correlation between W r and M. We demonstrate tat te alo occupation approac allows us to gain pysical insigts for (1) understanding te origin of intervening absorbers observed in te spectra of background QSOs, and (2) constraining te distribution of baryons in dark matter alos. We discuss tese results in te context of gas accretion in recent teoretical developments (i.e., Birnboim & Dekel 2003; Dekel & Birnboim 2006; Kereš et al. 2005; Birnboim et al. 2007). Tis paper is organized as te following. In 2, we outline te teoretical framework for constructing te conditional probability distribution of equivalent widts P(W r M). In 3, we derive te observables tat caracterize te statistical properties of Mg II absorbers, based on known dark matter alo statistics and P(W r M). In 4, we compare te model predictions wit empirical measurements for Mg II absorbers. In 5, we apply te best-fit model to predict additional statistical quantities on te galaxy absorber correlation tat can be tested wit future observations. Finally, we discuss in 6 te interpretations of our analysis. We adopt a flat ΛCDM cosmology wit Ω M = 0.25 and σ 8 = 0.8 and a dimensionless Hubble constant = H 0 /(100 km s 1 Mpc 1 ) trougout te paper. 2. THEORETICAL FRAMEWORK

3 3 In tis section, we describe te approac to populate dark matter alos wit cold baryons tat are represented by te presence of Mg + ions. An implicit assumption of our alo occupation approac is tat all absorbers originate in alos, rater tan in metal enriced intergalactic medium (IGM) outside of individual alos. Tis assumption is justified by te strong correlation between absorbers and galaxies; te LRG Mg II cross-correlation function exibits a strong signal at separations less tan alo radii (Boucé et al. 2004). In addition, for eac absorber a galaxy is nearly always found at a separation significantly less tan te typical alo radius expected given te galaxy luminosity. Finally, te multi-component nature of Mg II systems demonstrates tat te absorbing medium is igly clumped, in contrast to te smooter gas in te IGM. For eac dark matter alo, we will specify two quantities: (1) te mean absorption strengt per unit surface mass density A W (M), and (2) te mean covering factorκ g (M) of te gaseous clouds. Tese parameters determine te conditional probability distribution of W r as a function of alo mass, P(W r M) wic, wen combined wit known statistics of dark matter alos, allows us to derive teoretical predictions of various statistical properties of Mg II absorbers for direct comparisons wit observations Constructing te Gaseous Halo A dark matter alo is defined as a collapsed and virialized object wit a mean interior density equal to 200 times te background density. Tus te radius of a alo is [ R = R 200 = 3M 4π (200 ρ m ) ] 1/3. (1) Te coice of 200 ρ m is motivated by estimates of te virial radius from te sperical collapse model, wic predict a value of 180. Wile tese estimates are idealized, tey also matc te alos identified in N-body simulations, from wic te global properties of alos are calibrated. According to Equation (1), a M alo as a co-moving radius of kpc, wile R 200 for a M alo is Mpc. Wit tis definition, all alos of mass M ave te same comoving radius at any redsift. We adopt te alo mass profile of Navarro et al. (1997) and te concentration-mass relation of Bullock et al. (2001) (wit updated parameters from Wecsler et al. 2006) to calculate te alo mass as a function of radius. We begin wit te empirical result tat te equivalent widt W r of an absorption system is proportional to te number of absorption components, wic we interpret as gas clouds or clumps encountered along te sigtline troug te alo 1. We adopt an isotermal profile wit a core radius a to approximate te mean gas density tat follows were f g is te gas fraction, ρ g (r) = f g G 0 ( r 2 + a 2 ) 1, (2) G 0 = M(< R g )/4π R g a tan 1 (R g /a ), (3) and M(< R g ) is te mass of dark matter witin te effective radius R g of te gaseous alo. Beyond R g, no absorption is 1 We note tat an absorption line of widt W r implies a maximum velocity difference between components of 100 km s 1 per Angstrom. Te maximum circular velocity of a M NFW alo is 150 km s 1. found. We note tat an isotermal distribution of discrete gas clumps is expected to ave a mean absorption strengt tat declines wit radius. At a certain impact parameter te mean number of clouds per line of sigt will rapidly decline and become negligible (see, e.g., Figure 1 in Celouce et al. 2007). We approximate tis effect by creating an effective gas radius, R g, beyond wic no gas clumps are present. Considering te total absorption equivalent widt W r as a sum over all clumps encountered along a sigtline leads to [ R 2σ 2 g s2 cl W r (s M) =W 0 ρ g ( ] s M 2 + l 2 )dl (4) cl 0 ( ) W0 σ cl f g 2G0 = tan 1 R 2 g s 2 M cl s 2 + a 2 s 2 + a 2, (5) were s is te impact parameter of te sigtline wit respect to te center of te alo,σ cl and M cl are te cross section and mean mass of individual gas clumps, and W 0 is te absorption per clump. It is clear tat te parameters tat control te relation between impact parameter and W r, suc as te gas fraction f g, and te mean absorption strengt W 0 and crosssection σ cl of individual clumps, are degenerate. We express teir product as a single free parameter A W W 0σ cl f g M cl, (6) wic represents te mean absorption equivalent widt per unit surface mass density of te cold gas. In te limit tat all te components of te absorption line are optically tin, A W would be a measure of te mean particle density of Mg + ions. Some fraction of te lines will be optically tick, owever, so we are unable to constrain te density from te values of A W. Te relation between te predicted absorption equivalent widt W r at impact parameter s of a given alo is described by 2G 0 W r (s M) = A W (M) tan 1 R 2 g s 2 s 2 + a 2 s 2 + a 2. (7) Equation (7) sows tat if te fraction of cold gas and its properties do not cange across all alo masses, ten A W is expected to be a constant and te total mean equivalent widt is expected to increase wit alo mass due to te larger gas mass in more massive alos. We coose a value of a = 0.2R g, noting tat our results are insensitive to te exact coice for tis parameter. We set R g = 80 comoving 1 kpc (50 1 kpc pysical at z = 0.6) for alos of mass M. Te gas radius scales as R g = 80(M/10 12 ) 1/3 1 kpc, te same scaling as te alo radius in equation (1). Wit tis normalization, te gas radius is rougly one tird of te alo radius for all alo masses. Te coice of 50 1 kpc for our fiducial radius is motivated by observations of te radial extent to wic Mg II absorption of W r 0.3 Å is detected around L galaxies (e.g., Steidel 1993). We will demonstrate in 5 tat te value of R g makes distinct predictions for te distribution of impact parameters, and tat our coice for R g matces available observations well Te Conditional Probability Distribution of W r A Smoot Gaseous Halo Model To construct te conditional probability distribution of W r, P(W r M), we consider te probability distribution of a random

4 4 line of sigt intercepting a alo at impact parameter s and te mean absorption expected at te impact parameter s. For a random line of sigt, te probability tat a sigtline passes troug a alo of mass M wit impact parameter s is { 2s/R 2 P(s M) = g if s R g (8) 0 if s>r g. Te distribution of W r is related to te distribution of impact parameters according to P(W r M)dW r =κ g (M)P(s M)ds, (9) were κ g is te total probability of detecting an absorber in a alo of mass M, wic may be less tan unity. Tis integrated probability tat a sigt line troug a alo yields an absorber is governed by two pysical quantities: te mean covering fraction of cold gas per alo and te fraction of alos tat ost cold gas clouds. Tese two quantities are degenerate, and we parameterize teir product wit te quantityκ g. From equations (8) and (9), te probability distribution function of equivalent widts is P(W r M) =κ g (M) 2s(W r M) ds R 2, (10) g dw r were s(w r M) is te inversion of equation (7), wic is performed numerically. Te derivative of s wit respect to W r is ds = (2A W G 0 ) 1 s tan 1 x 1 + x2 (a 2 1 dw r (s 2 + a R2 g )(s x) )3/2 s 2 + a 2 (R 2 g s2 )(a 2 +, s2 ) (11) were x (R 2 g s2 )/(a 2 + s2 ). Instead of parameterizing te mass dependence ofκ g using power-laws or polynomials, we adopt a non-parametric approac and specify te values of κ g at four different masses and spline interpolate between tem. Tus, te form of κ g is fully non-parametric. Te four mass values, M 1 toug M 4 are logm i = 10.0, 11.33, 12.66, and We denote te four parameters of te κ g function as κ 1, κ 2, κ 3, and κ 4. We perform te interpolation in logm and logκ g and extrapolate to iger and lower masses. In tis smoot gaseous alo model, wic we will refer to as te classical model, we adopt a double power law to account for te mass dependence of te mean absorption per unit surface mass density A W, i.e., { A0 M12 α if M 12 1 A W (M) = A 0 M β (12) 12 if M 12 > 1. As a result, tis classical model contains seven free parameters to be constrained by te data. Tree govern te amount of absorption as a function of alo mass, A 0,α,β, and four caracterize te mass dependence of gas covering fractionκ 1 -κ 4. In summary, te gas density profile determines te relative absorption between different impact parameters in a alo te density profile sets te sape of P(W r M). Te relative absorption between alos of different masses is determined by A W ; at s = 0, te line of sigt directly troug te center of a alo, te absorption scales as W r (0 M) M 1/3 12 (A W/10.9) were A W is in units of ( Å cm 2 /gm). Te incidence of te absorption is set by κ g. Some sigtlines tat pass troug a alo will encounter no absorbing gas, eiter due to te fact tat te absorbing gas is igly clumpy and te covering fraction is low or due to te fact tat te alo does not contain cold gas Incorporating a Cold-Mode/Hot-Mode Transition As we summarized in 1, te Mg II absorbers are known to originate in primarily cool, poto-ionized gas. Recent numerical simulations sow tat stable socks develop as gas accretes onto a dark matter alo, and above a critical mass tresold all gas is sock eated to te virial temperature or above. At low mass scales, owever, gas accreted onto alos is not sock eated, but remains cold as it sinks to te center of te alo (e.g., Kereš et al. 2005). Dekel & Birnboim (2006) demonstrate tat tis is due to te gas cooling time being too sort wit respect to te compression time in low mass alos, tus socks are not stable at any radius witin te alo. It is terefore expected tat alos above a certain mass scale may not contribute significantly to te observed Mg II absorbers due to sock eating of teir gas. In tis second scenario, we include a transitional mass scale, above wic te gas becomes predominantly sock-eated and unable to contribute to Mg II absorption. We denote te sock radius R s. Te absorption efficiency at r < R s is reduced to a fraction f cold, wic represents te fraction of gas tat is cold witin te sock radius relative to gas outside te sock. For impact parameters smaller tan te sock radius, equation (7) is modified as te following, 2G 0 W r (s M) =W r (s M) (Rs=0) (1 f cold )A W tan 1 s 2 + a 2 R 2 s s2 s 2 + a 2, (13) were te first term on te left and side is te result from equation (7), in wic R s = 0. Equation (11) is calculated for s<r s in te same manner. Figure 1 sows P(W r ) (normalized to κ g = 1) for alos of M (left panel) and M (rigt panel). Te value of A W is set suc tat a sigtline at an impact parameter of R g /2 troug a M alo yields an equivalent widt of 1.5 Å for comparison purpose. In eac panel, we present tree different probability distribution functions (PDFs). Te dotted curve sows equation (10), wic represents te W r distribution for a continuous gas density profile wit no sock. As te impact parameter approaces te edge of te alo, W r becomes small. It is more likely to go troug at alo at large impact parameters, but te projected density drops steeply as s approaced R g, tus P(W r ) = 0 at W r = 0. Te mode of te distribution occurs at s 0.9R g, and at larger W r te curve follows Wr 3. Te distribution reaces a maximum value of W r 3.6 Å at s = 0. Te dased curve sows te distribution of equivalent widts for a alo wit R s = 0.3R g, and f cold = 5%. Te primary effect of sock eated gas in te inner spere is to reduce te maximum W r. In te case of a M alo, te largest equivalent widt is reduced by a factor of two. For s < R s, P(W r ) as a very different form. As s falls just inside te sock radius, W decreases rapidly. At s R s /2, W r (s) reaces a minimum value as te small fraction of cold gas witin R s accumulates enoug to contribute to te total absorption. Te minimum value of W r depends on f cold, but around te minimum, te derivative dw r /ds becomes small, producing a sarp feature in te total P(W r ) curve at W r 1 Å. Tis feature is even more pronounced in te PDF for te iger mass alo (rigt panel), also wit R s = 0.3R g. Te sarp feature is due to te idealized nature of te calculation. In practice, te propagation of te socks is unlikely to be sperically symmetric, and a alo of Poisson-distributed

5 5 FIG. 1. Te distribution of equivalent widts P(W r M) for alos of M (left panel) and M (rigt panel). In bot panels, te dotted curve is P(W r M) for a model in wic te cold gas follows an isotermal density profile trougout te alo. Te dased curves represent at model in wic te cold gas witin R = 0.3R g is reduced by a factor of 1 f cold = Te sarp feature in te PDF is due slow cange in W r as s approaces 0. Te solid curve represents a model were tis feature is smooted out, as expressed in equation (14). clouds will also smear out tis feature. To facilitate te calculations, owever, we approximate te distribution of W r (s < R s ) as a Gaussian centered on µ W = W r (s = R s /2) wit widtσ W = [W r (R s )/2 W r (0)]/2. Te amplitude of te Gaussian function is set by te fraction of te cross section witin te sock radius, wic for R s = 0.3R g is 9%. Te Gaussian function is ten added to te PDF for s > R s. For tis socked alo model, te probability distribution is P(W r M) =κ g (1 f σ ) s(w r M) ds R 2 + g dw r [ κ g f σ (Wr µ W ) 2 ] exp, W r W max,(14) 2πσ 2 W 2σ 2 W were f σ = R 2 s /R2 g is te probability tat a sigt line as an impact parameter witin R s, and W max = max[w r (R s ),W r (0)] from equation (4). We note tat P(W r ) = 0 for W r > W max. Te results are sown in te solid lines in bot panels of Figure 1. Instead of a sarp peak in P(W r ), te smoot model creates a soulder in te PDF, smooting out te discontinuity in te dased curve. Te parameters of te Guassian are cosen suc tat tis soulder tracks te peak in te PDF as R s canges. In te one-dimensional simulations of Birnboim & Dekel (2003) and Dekel & Birnboim (2007), sock eating first become stable at te inner regions of te alo, ten propagates outward as te alo evolves and becomes more massive. Tus at a given epoc, te widt of te transition region in teir simulations is rougly a factor of 2 5 in alo mass. In more realistic tree-dimensional simulations, but in wic te sock fronts are not fully resolved, te fraction of cold gas decreases linearly wit logm over 1 dex in mass (Kereš et al. 2005; Birnboim et al. 2007). We parameterize te mass dependence of R s as R s = ˆR 0 s +γ s log R 10 M 12, (15) g and restrict R s /R g to be no smaller ten 0 and no larger tan 1. We will refer to te transition scale as te mass at wic R s /R g = 0.5. At low masses te sock radius is zero, wile in te transition region te fraction of te gas radius tat is witin te socked core increases linearly wit logm until R s = R g. For alos wit R s = R g, te equivalent widt is calculated as 2G 0 W r (s M) = f cold A W tan 1 R 2 g s 2 s 2 + a 2 s 2 + a 2. (16) Nominally, tis transition model as eigt free parameters: A W, κ 1 troug κ 4, ˆR 0 s, γ s, and f cold. In tis scenario, we leave A W to be invariant wit alo mass and attribute relevant mass-dependent gas fraction to R s and f cold, in order to simplify te model analysis. Removing freedom from A W in favor of R s is beneficial in tat it igligts te differences between te two scenarios and allows us to isolate features from eac model. As we will sow in 4, te classical model and te transition model ave distinct predictions for various statistical quantities of te Mg II absorbers. Tese distinctions can be attributed to te different spatial distributions of te cold gas in alos of different mass, and serve as strong constraints of te models. 3. DERIVATION OF THE OBSERVED STATISTICAL PROPERTIES OF MG II ABSORBERS Bot te frequency distribution of Mg II absorbers as a function of W r and teir clustering amplitude wit respect to LRG are known from previous work. Here we summarize tese known statistical properties from observations and derive teoretical formulae based on te two scenarios described in Te Frequency Distribution Function Te frequency distribution function is defined as te number of absorbers per absorption equivalent widt interval

6 6 per unit redsift pat lengt. Here we combine te results from Procter et al. (2006) for W r > 1Å wit te results from Steidel & Sargent (1992) for 0.3Å W r 1Å. Te Procter et al. (2006) data is nearly complete (> 95%) over te redsift range 0.4 z 2.2 for absorbers of W r > 1 Å. Te Steidel & Sargent (1992) data is complete to W r = 0.3 Å over te redsift range 0.2 z 2.2. For empirical measurements, we take te raw number counts from Procter et al. (2006) and divide te number counts by te total co-moving pat lengt of all te lines of sigt, l tot. Te total co-moving pat lengt is calculated by integrating te distance-redsift relation, weigted by te number of lines of sigt g(z) at eac redsift interval dz as follows, l tot = dzg(z) [ Ω m (1 + z) 3 ] 1/2. +Ω Λ (17) Equation (17) yields te total co-moving radial distance covered by te lines of sigt to 45,000 quasars available in SDSS DR3. Te total number of absorbers is N abs = 4835, wit l tot = Mpc. Tere is no systematic or cosmic variance error estimate for te frequency data, so we apply Poisson fluctuations for estimating random errors, including a 5% systematic uncertainty for eac of te 38 bins tat ave at least one absorber. Tis additional systematic error only alters te total uncertainty for te lowest few W r bins from te Procter et al. sample (W r = Å) and prevents tose data from driving te total χ 2 wen evaluating a model. For W r 1Å, we use te results from Steidel & Sargent (1992). Tis is a muc smaller sample, but is valuable for constraining te sape of d 2 N/(dWdl). Te Steidel & Sargent measurements are scaled to tose of te Procter et al. (2006) data for 1Å< W r 1.3Å to account for te difference in te total redsift survey pats of te two samples. To obtain teoretical predictions, we calculate te conditional probability of a alo containing an absorber of a given strengt P(W r M) as specified in 2 for a given set of parameters. Te integral of te product of te cross-section weigted alo mass function and te probability of detecting an absorber of W r in a alo of a given mass leads to te frequency of te absorbers versus W r, i.e. f (W r ) d2 N dw r dl = dm dn dm σ g(m)p(w r M), (18) were σ g (M) = πr 2 g is te cross section of te gas alo and dn/dm is te alo mass function. Te frequency distribution in equation (18) as units of Å 1 ( 1 Mpc) 1 in co-moving coordinates. We note tat two scaling factors need to be accounted for before a direct comparison between empirical measurements and model predictions can be done. First, te redsift distributions of te frequency data and te bias data are different. As discussed in 3.2 below, te bias data are limited by te redsift range for LRGs at 0.35<z<0.8. We terefore must re-normalize te frequency function to te effective redsift of te B06 sample at z = 0.6. Te amplitude of te frequency function rises between z = 0.5 and z = 1, and te ratio of dn/dz at z = 0.6 to te full sample is Second, te amplitude of f (W r ) based on a SDSS DR5 Mg II absorber sample is found to increase by a factor of 1.1 (J. X. Procaska, private communication). We ave terefore included tis correction in our empirical data for comparisons wit models Te Mg II LRG Cross-Correlation Function Te cross-correlation function between absorbers and galaxies along common lines of sigt determines te clustering strengt (and terefore te relative bias) of te absorbers wit respect to te galaxies. Given te mean mass scale of te alos tat ost te galaxies, one can derive te mean mass scale of te alos osting te absorbers. Boucé et al. (2006) presented te projected cross-correlation function of Mg II absorbers and LRGs on scales of r p = Mpc. Te Mg II -LRG sample covers a redsift range of 0.35<z<0.8. Tese autors calculate te relative bias of absorbers to LRGs by taking te ratio of absorber-lrg cross-correlation function to te LRG autocorrelation function, ˆb W b W b G = ξ W G ξ GG. (19) Equation (19) is only exact at large scales because it assumes linear bias, but te scale dependence of alo bias is largely independent of alo mass (Tinker et al. 2005) and will cancel in te ratio. We take te absolute bias of LRGs wit respect to dark matter to be b LRG = 1.85 at z = 0.6 (Wite et al. 2007, wo adopt te same cosmology as wat is assumed ere), and estimate te absolute bias of te absorbers following Equation (19). To compare te models to te data, we calculate te absolute bias of absorbers and divide te calculation by b LRG. Te bias data are in four bins of W r : [0.3,1.15] Å, [1.15,2.0] Å, [2.0,2.85] Å, and [2.85,4.0] Å. Te igest-w r bin as little influence on te parameter constraints due to te small number of objects in tat bin. Wen comparing models to bias data, we take te number-weigted average bias over te given range in W r. Te B06 sample is incomplete for W r < 1 Å. Te true frequency distribution monotonically increases wit smaller W r, but te [0.3,1.15] Å bin is dominated by 1 Å absorbers. Wen calculating ˆb W in tis bin, we introduce te same incompleteness in te model by multiplying P(W r M) by 21.4(1 W r ) for W r < 1 Å. Wit tis factor, a model tat accurately re-produces te true d 2 N/(dW r dl) data will produce a frequency distribution tat matces te sape of te B06 sample. In te context of te alo model, te bias of absorbers wit respect to te dark matter is equal to te mean alo bias, weigted by te probability of finding an absorber W r at mass M; b W = 1 f (W r ) dm dn dm σ g(m)b (M)P(W r M), (20) were b is te alo bias. For te alo mass function we use te fitting function of Warren et al. (2006), wile for te alo bias we use te fitting function of Tinker et al. (2005). All alo properties are calculated at te effective redsift of te B06 absorber sample. Te redsift range of te data is large enoug suc tat integrating equations (18) and (20) over redsift will yield sligtly different results, but for te purpose of tis study it is sufficient to coose an effective redsift for all 2 We note tat te slope of te full Procter et al. (2006) sample is d log N/dW = 0.72, a value consistent wit te results from Nestor et al. (2005), wo calculate d log N/dW = 0.75 over 0.4 < z < 0.8 for a smaller sample culled from te SDSS early data release.

7 7 calculations. In our subsequent analyses, te redsift evolution of f (W r ) will be used as an additional constraint on te models. 4. COMPARISONS BETWEEN OBSERVATIONS AND BEST-FIT MODELS Te scenarios discussed in 2 ave distinct predictions for te empirical data discussed in 3. We constrain te models using a Monte Carlo Markov Cain tecnique tat determines te best-fit parameters based on comparisons of te data wit random realizations from various model predictions. For eac realization in te cain, te totalχ 2 is te sum of teχ 2 values from bot te frequency and bias data. We will refer to tese quantities as χ 2 f and χ 2 b, respectively. Because of te larger number of frequency data and smaller errors on tose data, χ 2 f dominates te total χ 2. However, te χ 2 f values sould be taken wit caution, because we ave no true estimate of te cosmic variance Te Classical Gaseous Halo Model As described in 2.2, te classical model contains seven free parameters, including tree tat caracterize te absorbing gas column, A W, for alos of different mass scales and four tat caracterize gas covering fraction κ 1 -κ 4 of individual alos. Te mass dependence in A W allows te predicted W r versus impact parameter s relation (Equation 7) to vary smootly wit alo mass. Te best-fit parameters are listed in Table 1. Figure 2 compares observational data wit predictions based on te classical model (dotted curves). We present in panel (a) te frequency data and in panel (b) te residuals normalized by te corresponding error of eac bin. Te frequency data are adequately fit by tis model, yielding χ 2 f = 43.3 for 41 data points and seven free parameters. Te classical model matces te abundance of strong absorbers, W r > 2 Å, but at lower equivalent widts te residuals become large. Specifically, tis model overproduces te number of absorbers by nearly 60% at W 0.3 Å due to a large contribution from low-mass alos. We present te comparison for te alo bias of te absorbers in panel (c) and te inferred mean alo mass in panel (d). Te bias data (points wit errorbars) are derived from Mg II LRG cross-correlation measurements. Te relative bias ˆb W derived from observations decreases from 0.8 at W r 0.8 Å to nearly 0.4 at W r 2.4 Å. Te measured bias for W r 2.8 Å increases sligtly to 0.5, but te number of objects in tis bin is small and it as little influence on our fitting results. Using a simple mean-mass approximation, B06 infer a decrease in mean alo mass of nearly 1.5 dex over tis range in W r. Te best-fit classical model is unable to reproduce te inverse correlation between ˆb W and W r. As expected, te bias monotonically increases wit W r, placing te weakest absorbers in M = M alos wile locating te strongest absorbers in te nearly same alos as LRGs. Figures 3a-3b present te parameter constraints for tis model. In bot panels, te inner and outer saded regions sow te allowed parameter range, constrained respectively by χ 2 < 1 and < 4 from te best-fit model. Figure 3a presents te constraints on te mean gas covering fraction κ g as a function of alo mass. Te function peaks at logm 12; nearly all alos of logm are expected to produce an absorber if te impact parameter is less tan R g to a background quasar. At lower and iger masses, κ g rapidly falls off, but te constraints become weak at te extrema of te mass range. Given te constraints, tere is little contribution to Mg II absorption from alos more massive tan M and less massive tan M. Figure 3b sows constraints for te absorption column A W (M). For low-mass alos, te constraints on te slope and normalization are strong; M alos produce 1Å absorbers. Te best-fit model as a single power-law, wit α β Te negative values of α and β increase te absorption efficiency in lower mass alos, but not sufficiently to produce an anti-correlation between W r and M. A slope more negative tan 1/3 in A W (M) versus M is required to make te mean W r constant as a function of M, but tese negative values would violate te constraints on te frequency distribution by overproducing ig-w r absorbers. At M > M, te constraints on A W are poor because κ g falls off rapidly at tese masses. Tere is a degeneracy between κ 4 and β; if κ g falls off for M > M, ig-w r systems can only be produced ifβ is large. However, tis steepens te slope ofˆb W and tese models fare even worse in comparison to te bias data. Te models wic produce te lowestχ 2 b ave negative values of botαandβ. Figure 3c sows te relative contribution to absorber counts from alos of different mass based on te best-fit model. We refer to tis mass distribution as te occupation function. Te curves sow tree values of W r : 0.5 Å, 1 Å, and 2 Å, and te relative areas under te curves represent te relative number of absorbers. Te y-axis, in arbitrary units, is te probability of detecting an absorber of W r in a alo of mass M, i.e., P(W r M)P(M), were P(M) is te probability tat a given line of sigt intersects a alo of mass M based on a crosssection weigted alo mass function, σ g (M)dn/dM. Tese curves ave te Gaussian-type sape of te κ g function, wit te sarp edges near te peak resulting from te mass range over wic κ g is truncated at unity. Te mass dispersion is 1 dex in logm for all tree curves. Te modes of eac curve move to iger mass, producing te increasing trend of W r wit M sown in Figure 2d. In summary, te classical model wit a smoot gaseous alo cannot fit bot te frequency data and te bias data simultaneously. Under tis scenario, more massive alos are expected to contain more absorbing gas, resulting in a positive correlation between W r and alo mass M. Forcing te mean absorption A W per alo to decline steeply wit increasing alo mass would improve te model fit to te bias data, but at te same time it would over-produce te abundance of strong (W r > 2 Å) absorbers Te Significance of Sock-Heated Gas in Massive Halos Te best-fit parameters of te transition model are also listed in Table 1. Tis transition model as eigt free parameters: A W, κ 1 troug κ 4, ˆR 0 s, γ s, and f cold. As discussed in 2.3, we ave left A W to be invariant wit alo mass and attributed relevant mass-dependent gas fraction to R s and f cold. Comparisons between observations and predictions based on te transition model are also presented in Figure 2 (solid curves). Te addition of te cold-ot transition produces a best-fit model tat precisely fits te frequency data over te entire W r range, including tose points at W < 1Å, witχ 2 f = Te model also reproduces te anti-correlation between W r and ˆb W for W r < 2.5 Å, sown in Figure 2c. Altoug te slope

8 8 FIG. 2. Panel (a) Te frequency distribution of Mg II absorbers as a function of equivalent widt W r. Te points wit error bars represent a combination of te data sets of Procter et al. (2006) (W r > 1Å) and Steidel & Sargent (1992) (W r < 1Å). Te transition between te two data sets is seen easily in te relative error bars between low and ig W r. Te solid and dotted curves represent te best-fit models listed in Table 1. Te solid curve is te model wit a cold-ot transition in alo gas, wile te dotted curve as a smootly varying cold gas fraction as a function of alo mass. Te χ 2 values are broken into te χ 2 for te frequency and bias data, respectively. Panel (b): Te residuals of te model f (W r) curves, normalized by te errors on te data. Panel (c): Te bias of absorbers relative to LRGs,ˆb W, as a function of W r. Te points wit errors are te measurements of Boucé et al. (2006). Te solid and dotted curves are te two models from panel (a). Panel (d): Te mean logaritmic alo mass as a function of W r. Points wit error bars are te values inferred from te bias data by B06 (not used in te fit). Solid and dotted curves are te two best-fit models. TABLE 1 PARAMETERS OF THE BEST-FIT MODELS Model A W ˆR 0 s γ s logκ 1 logκ 2 logκ 3 logκ 4 α β f cold χ 2 f χ 2 b Classical Transition NOTE. Te units of A W are [ Å cm 2 /gm]. of te trend is not as strong as tat seen in te B06 measurements, χ 2 b = 3.5 for four data points. Te relative bias in te model decreases fromˆb W = 0.7 at W r = 0.4Å to ˆb W = 0.55 for W r = 2Å. Tis declining trend in te absorber bias is due to a decrease in te mean alo mass wit increasing W r (Figure 2d). Te mean mass falls by nearly a decade from W r = 0.4Å

9 9 FIG. 3. Panel (a): Constraints on te covering fraction κ g as a function of alo mass from te classical model. Inner and outer contours represent te range in models tat produce χ 2 < 1 and 4 wit respect to te best-fit model. Panel (b): Constraints on te mean absorption strengt per unit gas mass density A W (M). Contours are for te range of models described in panel (a). Panel (c): Te occupation function of absorbers for tree different equivalent widt values: 0.5 Å, 1 Å, and 2 Å. Te y-axis, in arbitrary units, is te probability tat a sigtline passes troug a alo of mass M (P(M)) multiplied by te probability tat a alo of mass M yields a W r absorber (P(W r M)). to W r = 4Å. 3 Te total χ 2 for tis model is 23.3 for 37 degrees of freedom. In a purely statistical sense tis suggests tat te model as too many free parameters. However, we note once again tat te errors on f (W r ) are estimates only, and we wis to include as muc freedom in te model to prevent te frequency data from driving te constraints on te alo occupation. In tis sense, χ 2 b is te most relevant quantity, noting tat eigt free parameters is barely sufficient to reproduce te bias measurements. Figures 4a and 4b present te parameter constraints for te transition model. Te covering fraction is narrower and 3 We note tat te masses derived by B06 depend on bot te model used and te cosmology adopted, and tey are not utilized in te fit. Even a perfect fit to te bias data in panel (c) will not necessarily yield te B06 masses in panel (d). sarper tan te classical model: κ g = 1 for alos of logm and cuts off sarply at lower masses. Over nearly a decade in alo mass, all models wit χ 2 < 4 yieldκ g = 1. However, te transition from cold mode (R s /R g = 0) to ot mode (R s /R g = 1) begins at M M and ends at M M. Tis indicates tat massive alos of log M 12 contribute to te observed Mg II statistics troug te presence of cold flow in te sock-eated alos. Te constraints on te cold fraction witin te sock are strong, wit te best-fit model yielding f cold = 0.06 wit an allowed range of ±0.03 for models wit χ 2 < 4. Bot te transitional mass scale and cold gas fraction in massive alos agree well wit te simulations of Kereš et al. (2005) and Birnboim et al. (2007). In te post-sock regime log M > 12.5, te constraints on te covering factor are poor; models wit a low κ 4 yield better fits to te frequency data wile ig covering factors at

10 10 FIG. 4. Panel (a): Constraints on te covering fraction κ g as a function of alo mass from te transition model. Inner and outer contours represent te range in models tat produce χ 2 < 1 and 4 wit respect to te best-fit model. Panel (b): Constraints on te sock radius R s (relative to te gas radius). Contours are for te range of models described in panel (a). Te orizontal bands to te rigt of te R s contours represent te constraints on f cold. Panel (c): Te occupation function of absorbers from te best-fit model. Te abrupt transition at M = M is due to R s becoming larger tan R g M produce iger bias factors and lowerχ 2 b. Figure 4c sows te occupation functions for te best-fit model. For alos of M = M, te absorber counts rises asκ g increases wit alo mass as sown in panel (a). For W r = 0.5 Å, te occupation function peaks are M, ten abruptly drops wen R s = R g and all alos enter te ot pase. For W r = 1 and 2 Å, te peak in te occupation functions occur at M and M, respectively. Tis trend is due to te inside-out nature of te eating; iger-w r systems are produced at smaller impact parameters and te sock radius envelopes tat radius at lower alo masses. Above M, te occupation as a tail out to M. In te best-fit model,κ 4 = 0.67, implying tat a significant fraction of cluster-sized alos contain some cold gas wit a ig covering fraction. However, we note again tat te constraint on κ 4 is poor due to te low frequency of suc massive alos. For W r 1 Å absorbers, alos in te post-sock regime do not initially ave enoug gas to produce absorbers of tis strengt, and so te distribution of alo masses occupied by tese absorbers becomes bimodal; only at M > M is tere enoug cold gas to produce absorbers of tis strengt. For W r = 2 Å, for example, alos between and M alos do not produce strong absorbers. Tis bimodal distribution is wat causes ˆb W to increase for W r > 3 Å, wile logm remains constant in Figure 2; bias increases non-linearly wit mass, weigting ig-mass alos more inˆb W tan in logm. In summary, te transition model can simultaneously reproduce te observed frequency data and bias data of Mg II absorbers. Te anti-correlation between W r and M in te transition model is produced troug two mecanisms. First, as te sock propagates outward from te center of a massive

11 11 FIG. 5. Te distribution of equivalent widt as a function of impact parameter. In bot panels, te points represent te mean W r as a function of log s from a combination of four sets of data, wile te errors sow te dispersion about te mean. Panel (a) sows te predictions for te transition model. Te solid curve plots te mean W r wile te saded region represents te predicted dispersion. Te dotted curve is te upper-99% bound on te distribution of W r at eac log s. Panel (b) presents te predictions for te classical model. Te small open squares represent eac individual system in te observational sample. Te x-axis is in pysical units. alo, ig-w r systems are no longer produced, but low-w r systems are produced more frequently because sigt lines at low impact parameters (s R s ) now produce weaker absorbers. Tis is demonstrated in Figure 4c by te location of te cutoff in te occupation function as W r increases te location of te cutoff moves to lower masses. Second, te small amount of cold gas in te post-sock regime (masses at wic R s = R g ) yields mostly low-w r absorbers in ig mass alos. An instantaneous transition between te cold-mode and ot-mode wit a non-zero cold fraction could produce a weak W r M anti-correlation. A broad transition region (a small value of γ s ) wit f cold = 0 could also produce a weak anticorrelation. In tis model te two effects work togeter to fit te measurements. 5. ADDITIONAL OBSERVABLES PREDICTED BY THE MODELS We ave demonstrated tat te transition model can reproduce bot te observed frequency distribution function and te Mg II LRG relative bias data. Given te known W r (s M) constrained by te data in 4, and known mass distribution function of te dark matter alos, we can also predict (1) te correlation between W r and galaxy impact parameter s, and (2) te impact parameter distribution of Mg II absorbing galaxies P(s). Here we derive tese properties on te basis of te gaseous alo model described in 2 and te best-fit parameters constrained by te frequency and bias data in Table Te Absorber W r vs. Galaxy Impact Parameter Correlation To estimate te correlation between W r and impact parameter, as well as te distribution of impact parameters, we adopt te center of eac alo as te location of te corresponding absorbing galaxy. At a given impact parameter, te distribution of equivalent widts is ten determined by FIG. 6. Te distribution of impact parameters for all absorbers wit W r > 0.3 Å. Te saded istogram represents te data presented in Figure 5. Te solid and dotted curves sow te predictions from te transition and classical models, respectively. Te x-axis is in pysical units. P(W r s) = = dm dn dm σ g(m)p(w r M,s)P(s M) (21) dm dn dm σ g(m)δ(w r (s M) W r )κ g (M)P(s M)