arXiv:1309.6513v1 [astro-ph.GA] 25 Sep 2013 mCO mm 3 fnab aais h ao esn o sn Oas CO using for reasons components major gas The molecular and . atomic nearby the of study to ods ie(.. beckwe l 09 aemitre al. et Bauermeister 2009; cosmic through al. ISM et 2010). Obreschkow phase the 2010) ISM (e.g., of al. evolution of time et the theories Ostriker 2009; galax- constrain proposed al. nearby and test et of to Krumholz variety (e.g., used molecu- wide balance 2012). be a of al. over can distribution et gas ies Lada spatial atomic (e.g., and the lar rate of formation studies and global Resolved is fraction gas metallicity, the dense to the and hence predicting related morphology, to critical is mass, therefore gas e.g., pre- molecular conversion properties, is to the yet (ISM) atomic how 2007), medium of Understanding interstellar Ostriker atomic. & the dominantly McKee galaxies (e.g., most in formation star for oyWong Tony AM UVYTWR NRRDBIH ERYGLXE (STIN GALAXIES NEARBY INFRARED-BRIGHT TOWARD SURVEY CARMA edleg Germany Heidelberg, So 2006, Johannesburg Park, Africa Auckland Johannesburg, of sity USA Str Albert-Ueberle Germany Heidelberg, Heidelberg, Universit¨at 69120 der tronomie Canada 2E1, T6G AB USA 94720, CA ley, USA 22903, USA 20742, MD Park, USA 61801, rpittpstuigL using typeset Letters Preprint ApJ by Accepted vrteps e eae,te2 mH cm 21 the decades, few past the Over material raw the provides galaxies in gas Molecular 10 9 8 7 6 5 4 3 2 1 ot fiaSAPsdcoa Fellow Postdoctoral Univer- SKA Africa 524, South Box PO 140, Lab C1 Physics, of Department 87801, NM Socorro, Observatory, Astronomy As- Radio f¨ur National Zentrum Astrophysik, theoretische Edmonton, f¨ur Alberta, Institut of University Physics, of Berke- Department California, of University Astronomy, of Department V Charlottesville, Observatory, Astronomy Radio College National Maryland, of University Astronomy, of Department I Urbana, Illinois, of University Astronomy, of Department a-lnkIsiu ¨rAtooi,Kngth 7 6911 17, Konigstuhl f¨ur Astronomie, Max-Planck-Institut bevdH observed oundniymaue nrgosehbtn Oeiso.Sc depe H a balances Such which Krumholz, emission. & meta CO McKee a exhibiting for of evidence regions model reveals equilibrium in metallicity, the and measured mass density stellar column in and range sub-kpc at wide (STING) a Galaxies Nearby IR-Bright Toward Survey oursle lmig lhuhvle ls otemdlpredictio model the to close values although H observe clumping, unresolved to H fmlclrgsi ihl euae yteselrdisk. stellar the by dens regulated surface headings: tightly Subject stellar is the gas with molecular strongly of correlates CO from inferred eivsiaetecreainbtenC n H and CO between correlation the investigate We FAOI N OEUA A UFC ESTE NGLX PROP GALAXY ON DENSITIES SURFACE GAS MOLECULAR AND ATOMIC OF J I opnn htddntcnrbt oH to contribute not did that component 1 = ai .Fisher B. David 1 u Xue Rui , − 1. I iehv entepiaymeth- primary the been have line 0 I oundniismc agrta rdce,a ih eepce w expected be might as predicted, than larger much densities column A INTRODUCTION T oundniyi fe mle hnpeitdb h oe,a effect an model, the by predicted than smaller often is density column E tl mltajv 5/2/11 v. emulateapj style X 1 let .Bolatto D. Alberto , aais S S:aos—IM molecules ISM: — atoms ISM: — ISM galaxies: 2 J , re Ott urgen ¨ I ieadthe and line 7 uu Rahman Nurur , 2 cetdb p Letters ApJ by Accepted dmK Leroy K. Adam , 2, . ABSTRACT uth 2 A L 7 hedn.W lofidta h H the that find also We shielding. I msini 8nab aaisfo h CARMA the from galaxies nearby 18 in emission ad eodrl,rdainfil)gvr h H pressure the interstellar govern field) that radiation (1993) secondarily, Elmegreen (and, by model f rniin oethat Note transition. a et stogtt esgicn.Tepiaychal- H primary opti- The mapping the significant. systematically where be for to cases lenge thought cm in 21 is except from depth mass con- direct, cal gas In fairly atomic is of galaxies. inference method emission in the content preferred CO, gas to the trast molecular 2013 remains the al. CO survey et review), to Bolatto recent see a (“X-factor”; for factor scaling stant hs ept h netit ncnetn Ointensity CO H converting in to uncertainty the despite Thus, s xiainrqieet,hg bnac eaieto relative mod- abundance brightness, high intrinsic H requirements, its are excitation tracer est gas molecular a eeti hoeia oeigo h H the of modeling theoretical in terest exhibits but range. fraction dynamic gas molecular less the much as to referred erly a led niae htH that indicated already had sensitivity. and resolution angular high both achieving is ela osrc nltceulbimmdl suitable as models 2010), equilibrium al. et analytic Glover construct time-dependent and 2006; as a physical al. well in et the scales (Pelupessy small model manner on to processes attempts chemical including sition, oteselrds,tecmaio frda H radial of comparison confined the largely disk, CO stellar with the galaxies, to in differently very eo ta.(08 eeldsrn ailted nthe den- surface in mass trends hydrogen radial atomic sity, to strong molecular revealed of and (2008) ratio BR06), Blitz al. hereafter & et (2006, Wong Leroy Rosolowsky by & samples Blitz significant (2002), across profiles CO ftemdln a rsuebsdo yrsai equi- hydrostatic on based estimate pressure simple gas librium, a midplane with the correlation of tight a exhibit which 8,9 mol 2 nrcn er hr a enarsrec fin- of resurgence a been has there years recent In hl al bevtos(.. acn on 1986) Young & Tacconi (e.g., observations early While n cesblt ogon-ae ai telescopes. radio ground-based to accessibility and , tatN Vogel N. Stuart , 2 3 ≡ e Blitz Leo , oundniy yial efre sn con- a using performed typically density, column Σ H P 2 mp / t,sgetn httelclsupply local the that suggesting ity, Σ 2 gas ∝ p cls u ape spanning sample, Our scales. kpc omto n iscain The dissociation. and formation 4 = Σ rkRosolowsky Erik , r lose.W onot do We seen. also are n R lct eedneo h H the of dependence llicity gas R mol mol 2 √ ainWalter Fabian , R ρ dnei rdce by predicted is ndence ) I.TEDEPENDENCE THE III. G). / ≡ ∗ mol (1+ hsspotda earlier an supported This . Σ H srltdt h quantity the to related is r hr diffuse a there ere R 2 2 I / mol n Oaedistributed are CO and Σ oundensity column HI ,wihi oeprop- more is which ), eattribute we , I ERTIES 5 nnab galaxies nearby in rn Bigiel Frank , 10 I oH to I I 2 oH to I tran- and 6 , 2 2 Wong et al. for cloud scales and larger. The latter approach has performed using CASA12 (McMullin et al. 2007), taking been outlined in a series of papers by Krumholz et al. advantage of the multi-scale clean feature of CASA to (2008), Krumholz et al. (2009, hereafter KMT09), and better recover extended structures by searching for model McKee & Krumholz (2010, hereafter MK10). Roughly components on multiple scales (in our case, point sources speaking, an equilibrium H2 abundance arises from a as well as Gaussians twice and five times the width of competition between H2 formation on grains and pho- the synthesized beam). For four galaxies (NGC 3198, todissociation by UV radiation. Thus, for a given cloud 4254, 4536, and 5713) with CO J =2 1 single-dish the molecular gas fraction depends on the dust opti- maps from the HERACLES project (Leroy− et al. 2009, cal depth of the cloud (since dust excludes dissociat- Leroy et al. 2013b), a detailed comparison in X14 sug- ing radiation and encourages molecule formation) and gests that the CARMA maps recover &75% of the flux the strength of the external radiation field. A depen- within the field of view. 13 dence of fmol on metallicity enters via the correlation We searched the NRAO Data Archive for public H I of metallicity with the dust abundance. Already there 21 cm data of STING galaxies obtained in B, C, or D have been attempts to apply the empirical Rmol–Pmp re- configurations of the VLA with adequate velocity res- −1 lation of BR06 and the theoretical predictions for fmol olution (.20 kms ) to separate line from continuum. of KMT09 and MK10 to cosmological contexts (e.g., Archival data were found for 19 out of the 23 galaxies, Obreschkow et al. 2009; Fu et al. 2010; Lagos et al. 2011; so the present study is restricted to this sub-sample. In Kuhlen et al. 2012), underscoring the importance of test- addition, NGC 3593 shows unusually low deprojected H I ing their applicability across a broader range of galaxies. column densities, which we attribute to its high inclina- The CARMA11 Survey Toward IR-Bright Nearby tion and significant optical depth in the 21 cm line, so we Galaxies (STING) is a major effort to measure the molec- do not include it in the present analysis. To ensure the ular gas properties and their relation to star forma- most uniform possible data processing, we reprocessed tion in a sample of 20 nearby star-forming galaxies. the visibility data using CASA. We employed multi-scale Unlike the precursor∼ BIMA Survey of Nearby Galaxies clean for well-resolved galaxies (with map resolution (SONG; Helfer et al. 2003), the STING sample is specif- < 1 kpc), and flagged regions with strong continuum ically chosen to cover a wide range in stellar mass, and emission (Tb > 50 K) as potentially susceptible to opti- thus covers a wide range in associated galaxy properties cal depth effects. Further details on the pipeline as well (star formation rate, color, luminosity, metallicity, etc.). as data sources can be found in X14. The FWHM beam Supplementing STING CO data with ancillary multi- parameters of the resulting H I maps (major and minor wavelength data from the ultraviolet (UV) to infrared axis and position angle) are given in Table 1, and set the (IR), Rahman et al. (2011) studied the spatially resolved limiting resolution for our analysis. star formation law and its sensitivity to different analysis A smoothing and masking process was used to gen- methods in NGC 4254, and have recently extended this erate integrated line intensity images from the CO and study to the full STING sample (Rahman et al. 2012). H I data cubes. The cube was first smoothed along the The molecular gas data from the STING survey also per- spatial and velocity axes to increase the signal-to-noise mit a study of atomic-to-molecular gas transition, which ratio. A mask was then generated from the smoothed is the subject of this Letter. We compare the STING cube by starting at emission peaks defined by at least CO J =1 0 dataset with archival H I data from the two consecutive >4σ channels (where σ is the RMS noise NRAO Very− Large Array (VLA) to examine the atomic in the smoothed cube) and expanding outward to a 2σ and molecular gas relation in 18 galaxies of the STING edge. This mask was then applied to the original cube, sample on sub-kpc to kpc scales. Our principal observa- blanking regions outside the mask, before summing the tional result is a clear dependence of H I column density velocity channels. Finally, the integrated intensity im- on metallicity, and an equally striking dependence of H2 ages were deprojected to face-on using the orientation column density on stellar surface density. A more de- parameters given in Table 1 after being regridded to a tailed presentation of these data and comparison with common pixel scale and smoothed to a common resolu- theoretical models will be presented in a future work tion (the latter chosen to achieve a circular beam when (Xue et al. 2014, in preparation, hereafter X14). deprojected). For reference, Table 1 also gives, as a mea- 2. sure of the physical resolution, the length scale rphy cor- OBSERVATIONS AND DATA REDUCTION responding to the deprojected H I beam using the me- The STING sample is composed of northern (δ > dian -independent distance (e.g., from the Tully- 20◦), moderately inclined (i < 75◦) galaxies within 45 14 − Fisher method) reported by NED. Mpc selected from the IRAS Revised Bright Galaxy Sur- Face-on mass surface densities were calculated by us- vey (Sanders et al. 2003) to sample a wide range of stellar ing the optically thin assumption for H I and a CO-to-H2 −2 −1 −1 masses (and consequently, ISM metallicity). CO obser- conversion factor of αCO = 4.4 M⊙ pc (K kms ) vations were conducted in the C and D configurations of for CO (e.g., Leroy et al. 2008). Reported surface densi- CARMA to yield a robustly weighted synthesized beam ′′ ′′ ties include a factor of 1.36 to account for helium. The of 3 –4 across a field of view extending to roughly half typical 1σ uncertainty in the deprojected surface density, the optical radius. Calibration was performed using the MIRIAD package, but imaging and deconvolution were 12 Common Astronomy Software Applications, http://casa.nrao.edu 11 The Combined Array for Research in Millimeter-wave Astron- 13 http://archive.nrao.edu omy (CARMA) is operated by the Universities of California (Berke- 14 The NASA/IPAC Extragalactic Database (NED) is operated ley), Chicago, Illinois, and Maryland, and the California Institute by the Jet Propulsion Laboratory, California Institute of Tech- of Technology, under a cooperative agreement with the University nology, under contract with the National Aeronautics and Space Radio Observatory program of the National Science Foundation. Administration. Atomic and Molecular Gas Surface Density 3

TABLE 1 Sample Properties

× Name Dist. i P.A. MB 12+log(O/H) Metal Ref. bmaj bmin (P.A.) rphy σΣHI σΣH2 ◦ ◦ ′′ ′′ ◦ −2 −2 (Mpc) ( ) ( ) (mag) (dex) ( × ) ( ) (kpc) (M⊙ pc ) (M⊙ pc ) NGC 337 21.3 52 121 −19.99 ± 0.24 8.60 ± 0.34 1 28.5 × 20.5 (195) 4.74 0.09 0.61 NGC 772 30.2 37 315 −21.57 ± 0.49 8.87 ± 0.13 6 48.1 × 43.0 ( 88) 8.50 0.04 0.67 NGC 1156 7.0 52 272 −17.73 ± 0.34 8.18 ± 0.10 3 7.0 × 6.2 (336) 0.38 0.40 2.97 NGC 1569 2.5 63 112 −17.65 ± 0.96 8.13 ± 0.12 3 6.6 × 5.6 ( 17) 0.18 0.25 1.01 NGC 1637 9.8 39 213 −18.63 ± 0.56 8.80 ± 0.34 4 23.3 × 19.2 (169) 1.34 0.08 0.94 NGC 2782 16.1 30 75 −18.79 ± 0.97 8.59 ± 0.10 3 7.5 × 6.2 ( 87) 0.61 0.57 4.48 NGC 2976 3.6 65 335 −17.20 ± 0.50 8.67 ± 0.31 1 8.9 × 7.7 ( 46) 0.37 0.12 0.98 NGC 3147 40.9 32 147 −21.72 ± 0.37 9.02 ± 0.36 2 21.3 × 18.0 (136) 4.40 0.10 1.12 NGC 3198 14.0 72 215 −19.90 ± 0.28 8.62 ± 0.28 1 7.7 × 6.5 ( 66) 1.53 0.10 0.77 NGC 3486 15.6 36 263 −19.99 ± 0.59 8.75 ± 0.32 2 59.2 × 49.8 ( 79) 4.69 0.07 0.37 NGC 3593 5.5 67 90 −16.92 ± 0.08 8.29 ± 0.26 2 15.4 × 14.6 (346) 1.06 0.03 0.37 NGC 4151 6.6 21 22 −17.71 ± 1.40 8.41 ± 0.28 2 9.6 × 8.8 ( 66) 0.33 0.59 3.22 NGC 4254 15.6 31 69 −20.67 ± 0.23 8.79 ± 0.34 1 10.6 × 9.6 ( 95) 0.89 0.40 1.54 NGC 4536 14.7 68 301 −19.74 ± 0.37 8.61 ± 0.40 1 10.1 × 9.4 ( 81) 1.86 0.14 0.72 NGC 4605 5.5 67 300 −17.85 ± 0.41 8.43 ± 0.28 2 5.5 × 4.9 ( 72) 0.36 0.30 1.24 NGC 4654 16.1 62 125 −20.02 ± 0.28 8.83 ± 0.27 5 16.9 × 16.4 ( 55) 2.82 0.06 0.45 NGC 5371 27.7 48 11 −20.93 ± 0.55 8.90 ± 0.34 2 85.1 × 45.6 ( 81) 16.68 0.04 0.11 NGC 5713 21.4 33 203 −19.95 ± 0.53 8.64 ± 0.40 1 16.9 × 16.2 (328) 2.08 0.14 0.77 NGC 6951 23.3 46 138 −21.55 ± 0.27 8.99 ± 0.36 2 11.1 × 6.5 ( 48) 1.83 0.37 2.53

Note. — Metallicity references: 1: Moustakas et al. (2010); 2: Moustakas et al. (2010) (L-Z); 3: Engelbracht et al. (2008); 4: van Zee et al. (1998); 5: Pilyugin et al. (2004), Dors & Copetti (2006); 6: Anderson et al. (2010)

−1 based on a typical integration window of 20 km s for this is rather insensitive to the actual value of ΣH2 (and −1 CO and 40 kms for H I, is given in Table 1 and is ap- thus the adopted CO-to-H2 conversion factor), as indi- propriate for maps smoothed to the H I resolution. The cated by the nearly vertical shape of the model curves deprojected intensity maps were sampled onto a hexago- drawn in Figure 1(a). For the sample as a whole we find nal grid with a spacing between grid centers equal to the a reduced dispersion in ΣHI (from 0.23 to 0.17 dex) af- (deprojected) FWHM resolution of the image, following ter applying the MK10 prediction, as well as a shift of the methodology of Leroy et al. (2013b). some of the lower metallicity outliers into the main con- centration of points. We therefore see evidence for the 3. RESULTS metallicity dependence predicted by the model. Never- theless, the distribution is strongly skewed to the left of Figure 1(a) displays a point-by-point comparison of the dashed line, where observed and predicted column Σ and Σ at the limiting resolution of the H I data. HI H2 densities are equal. Thus, the model frequently overpre- Only regions where both H I and CO are detected at dicts Σ by factors of up to 5, especially for lower values the >3σ level are shown. Open triangle symbols indi- HI of Σ . The trends and scatter we observe remain un- cate H I measurements which are potentially contami- H2 changed when we repeat the analysis for 9 galaxies at nated by absorption in front of bright continuum emis- a fixed deprojected resolution of 2 kpc, suggesting they sion, and a diagonal line labeled R =1 divides H I and mol are not strongly influenced by varying r across our H -dominated regimes. Points are color coded by metal- phy 2 sample. licity, determined as follows. For six objects lacking spec- The tendency for the MK10 model to overpredict Σ troscopic gas-phase metallicities (NGC 3147, 3486, 4151, HI may be linked to a reduction in measured peak gas col- 4605, 5371, and 6951), we derive a fiducial metallicity for umn densities (relative to what would be measured on the galaxy by applying the luminosity-metallicity rela- GMC-like scales) as a result of spatial averaging over tion of Moustakas et al. (2010). For four galaxies (NGC kpc-scale beams. Because f is a steeply increasing 1637, 3198, 4254, and 4654) for which a metallicity gra- mol function of Σ in the model, the suppression of peaks dient has been determined from multiple spectroscopic gas in Σ by beam averaging causes the model to underes- measurements, we determine a local metallicity appro- gas timate the molecular fraction and thus overestimate Σ . priate to each grid sample. For the remaining galaxies, HI Higher values of Σ may be subject to less beam dilu- we assign a single metallicity value (given in Table 1) H2 tion and may better reflect the conditions under which across the galaxy, based on published H II region spec- the MK10 model is directly applicable (i.e., a resolved troscopy. (Since the gradients, when measured, are typi- GMC complex and its atomic envelope). This may ac- cally around 0.5 dex per R , and our observations are 25 count for the tendency for Σ /Σ to approach confined to R− 0.5 R , inclusion or omission of a gra- HI HI,MK10 25 unity for larger Σ . Further discussion of this effect dient has little≪ effect on our results.) We find that points H2 can be found in KMT09. for a given galaxy tend to be distributed vertically, con- To investigate trends in Σ independent of Σ , sistent with a nearly constant value of Σ when CO is H2 HI HI Figures 1(c) and (d) plot values of Σ against stel- detectable. However, horizontal offsets between galax- H2 lar mass surface density (Σ∗) and scaled galactocen- ies are clearly apparent, with lower metallicity galaxies tric radius, with points still color-coded by metallicity. shifted to the right. The stellar mass surface density is inferred from Spitzer In Figure 1(b), Σ has been normalized by the H I HI 3.6 µm imaging (S4G project; Sheth et al. 2010) using surface density predicted by MK10 (their Equations 13, Σ∗ = 317.7 I3.6 (Eskew et al. 2012), where I3.6 is the 86, 91, and 93) for the observed ΣH2 and metallicity; × 4 Wong et al.

Fig. 1.— (a) Comparison of H I and H2 surface densities sampled at the resolution of the H I data. Points are color coded by gas-phase metallicity, assumed uniform except for four galaxies (NGC 1637, 3198, 4254, and 4654) where a measured gradient has been applied. A diagonal line corresponds to the locus of equal H I and H2 surface densities. The correlation of characteristic H I column density with metallicity is apparent as a color gradient in the plotted points. The MK10 predictions for three different metallicities (12+log(O/H)= 9.1, 8.65, and 8.2) are shown as colored lines. (b) The same data and color coding but with the abscissa normalized by the predicted value from MK10 as appropriate to each point. Panels (c) and (d) compare ΣH2 with stellar surface density and normalized galactocentric radius respectively. background-subtracted 3.6 µm intensity in units of MJy factor and uniform mass-to-light ratio in the IR.) Aver- −1 sr . A strong correlation between ΣH2 and Σ∗ is appar- age values of ΣHI and ΣH2 (filled triangles), calculated in ent, whereas no relationship with R/R25 is discernible in uniformly spaced logarithmic bins of Σ∗, underscore the the sample as a whole. much greater sensitivity of ΣH2 to Σ∗. The same trend

Figure 2 provides another view of the ΣH2 –Σ∗ cor- is seen at a fixed resolution of 2 kpc (right panel), for the relation by showing the same data as Figure 1(a), but 9 galaxies observed at sufficiently high resolution, except with points color coded by Σ∗. At the full H I resolution that several of the low-metallicity systems drop out of

(left panel), it is clear that ΣH2 correlates well with Σ∗, the sample. We stress that the CO measurements are except for a few low-metallicity, H I-dominated galax- taken over areas much larger than an individual GMC, ies where ΣH2 is low despite high Σ∗. (Such exceptions so the observed trend probably originates from changes may reveal inaccuracies in our assumption a uniform X- in the covering fraction of GMCs rather than the actual Atomic and Molecular Gas Surface Density 5

Fig. 2.— (a) The same data points are shown as in Figure 1(a), but the points are color coded by stellar surface density as inferred from

Spitzer 3.6 µm surface brightness. Solid triangles denote average values of ΣHI and ΣH2 , calculated in deciles of Σ∗. There is a strong gradient in stellar surface density along the ΣH2 axis. (b) ΣH2 vs. ΣHI, again color coded by Σ∗, determined at a fixed resolution of 2 kpc for galaxies with rphy < 2 kpc. surface density of GMCs. (Leroy et al. 2013a), and be inconsistent with the bulk of the H I being confined to the GMC envelopes modeled 4. DISCUSSION AND SUMMARY by MK10.

We have compared the CO and H I surface brightnesses The correlation we observe between ΣH2 and Σ∗ (set- in individual apertures across 18 galactic disks from the ting aside low-metallicity systems for which the CO-to- CARMA STING sample spanning a range of metallicity H2 conversion factor or 3.6 µm mass-to-light ratio may and stellar mass. For simplicity we have plotted only in- change) is hardly unexpected: a correspondence between dependent quantities, avoiding intrinsic correlations be- CO and stellar radial profiles has been noted in many pre- tween axes and assumptions about which combination vious studies (e.g., Regan et al. 2001; Leroy et al. 2008, of observables is most physically relevant for star for- 2009; Schruba et al. 2011; Bigiel & Blitz 2012). Still, a mation or ISM physics. We also exclude non-detections comparison of Figures 1(a) and 2(a) vividly illustrates (upper limits) from our analysis, so we are sensitive only that the H I and H2 column densities of galactic disks to trends in CO-detected regions, which may be quite appear to be governed by distinct processes. It is possi- distinct from trends related to the detectability of CO ble that the two trends we observe, lower ΣHI with higher (e.g., Saintonge et al. 2011). metallicity (which often associates with higher Σ∗) and

Our principal result is a clear dependence of the the higher ΣH2 with higher Σ∗, reinforce each other to pro- characteristic H I column density on metallicity, as pre- duce a tight correlation of Rmol with Σ∗ and estimates of dicted by the models of KMT09 and MK10. A similar hydrostatic pressure. However, since Σ∗ is dominated by result had been obtained for a sample by the old stellar population, the processes that maintain

Fumagalli et al. (2010), but with less of a clear metal- its tight correlation with ΣH2 remain unclear. licity dependence. We find that a major success of the Studies probing GMC size scales find that the mass MK10 model is its ability to predict the upper envelope surface densities of GMCs are confined to a fairly narrow of the observed ΣHI values, although the full distribu- range (Bolatto et al. 2008; Roman-Duval et al. 2010), tion of points appears skewed below the predicted values implying that kiloparsec-scale variations in ΣH2 are by unresolved clumping within the beam. In particular, largely reflecting the CO areal covering factor. We there- we do not see evidence for a diffuse H I component that fore hypothesize that GMCs are more abundant in re- does not contribute to H2 shielding and would therefore gions of galaxies with higher Σ∗. A simple explanation cause the observed ΣHI to exceed the model prediction. is that star formation in GMCs is responsible for build- On the other hand, the fact that the predicted value of ing up the stellar disk, resulting in a stellar disk that

ΣHI is sometimes reached in low ΣH2 regions suggests naturally mimics the H2 distribution. However, obser- that H I and CO need not suffer equally from beam di- vational evidence for GMC lifetimes of only a few tens lution. This would imply a difference in how H I and of Myr (Kawamura et al. 2009; Murray 2011) raises the CO are spatially distributed, as found by recent stud- question of why the GMCs should be regenerated at ies of the SMC (Bolatto et al. 2011) and other galaxies the same locations in the disk. In present-day, highly 6 Wong et al. evolved galaxies with gas fractions 0.1, it seems more proved the presentation of results. T. W. and R. X. ∼ plausible to argue that the H2 distribution is governed acknowledge support from NSF University Radio Facil- by the stellar distribution. For instance, in the thermal ities Program and the NASA Astrophysics Data Analy- equilibrium model of Ostriker et al. (2010), higher stel- sis Program. Support for CARMA construction was de- lar density leads to higher pressure in the diffuse gas (as rived from the Gordon and Betty Moore Foundation, the required by hydrostatic equilibrium), increasing the cool- Eileen and Kenneth Norris Foundation, the Caltech As- ing rate and leading to the formation of bound clouds. sociates, the states of California, Illinois, and Maryland, Alternatively, H2 consumption could be regulated by Σ∗- and the NSF. Funding for ongoing CARMA development dependent gravitational instabilities (an idea explored and operations are supported by NSF and CARMA part- further by Zheng et al. 2013), or a significant fraction ner universities. The National Radio Astronomy Obser- of the total H2 supply may be recycled from the old stel- vatory is a facility of the National Science Foundation op- lar disk (an idea explored further by Leitner & Kravtsov erated under cooperative agreement by Associated Uni- 2011). versities, Inc. Facilities: VLA, CARMA, Spitzer We thank the referee for suggestions which have im-

REFERENCES Anderson, J. P., Covarrubias, R. A., James, P. A., Hamuy, M., & Leroy, A. K., Walter, F., Sandstrom, K., Schruba, A., Habergham, S. M. 2010, MNRAS, 407, 2660 Munoz-Mateos, J.-C., Bigiel, F., Bolatto, A., Brinks, E., de Bauermeister, A., Blitz, L., & Ma, C.-P. 2010, ApJ, 717, 323 Blok, W. J. G., Meidt, S., Rix, H.-W., Rosolowsky, E., Bigiel, F. & Blitz, L. 2012, ApJ, 756, 183 Schinnerer, E., Schuster, K.-F., & Usero, A. 2013b, AJ, 146, 19 Blitz, L. & Rosolowsky, E. 2006, ApJ, 650, 933 McKee, C. F. & Krumholz, M. R. 2010, ApJ, 709, 308 Bolatto, A. D., Leroy, A. K., Jameson, K., Ostriker, E., Gordon, McKee, C. F. & Ostriker, E. C. 2007, ARA&A, 45, 565 K., Lawton, B., Stanimirovi´c, S., Israel, F. P., Madden, S. C., McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, Hony, S., Sandstrom, K. M., Bot, C., Rubio, M., Winkler, K. 2007, in Astronomical Society of the Pacific Conference P. F., Roman-Duval, J., van Loon, J. T., Oliveira, J. M., & Series, Vol. 376, Astronomical Data Analysis Software and Indebetouw, R. 2011, ApJ, 741, 12 Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell, 127 Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., & Blitz, Moustakas, J., Kennicutt, Jr., R. C., Tremonti, C. A., Dale, L. 2008, ApJ, 686, 948 D. A., Smith, J.-D. T., & Calzetti, D. 2010, ApJS, 190, 233 Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, Murray, N. 2011, ApJ, 729, 133 207 Obreschkow, D., Croton, D., De Lucia, G., Khochfar, S., & Dors, Jr., O. L. & Copetti, M. V. F. 2006, A&A, 452, 473 Rawlings, S. 2009, ApJ, 698, 1467 Elmegreen, B. G. 1993, ApJ, 411, 170 Ostriker, E. C., McKee, C. F., & Leroy, A. K. 2010, ApJ, 721, 975 Engelbracht, C. W., Rieke, G. H., Gordon, K. D., Smith, Pelupessy, F. I., Papadopoulos, P. P., & van der Werf, P. 2006, J.-D. T., Werner, M. W., Moustakas, J., Willmer, C. N. A., & ApJ, 645, 1024 Vanzi, L. 2008, ApJ, 678, 804 Pilyugin, L. S., V´ılchez, J. M., & Contini, T. 2004, A&A, 425, 849 Eskew, M., Zaritsky, D., & Meidt, S. 2012, AJ, 143, 139 Rahman, N., Bolatto, A. D., Wong, T., Leroy, A. K., Walter, F., Fu, J., Guo, Q., Kauffmann, G., & Krumholz, M. R. 2010, Rosolowsky, E., West, A. A., Bigiel, F., Ott, J., Xue, R., MNRAS, 409, 515 Herrera-Camus, R., Jameson, K., Blitz, L., & Vogel, S. N. 2011, Fumagalli, M., Krumholz, M. R., & Hunt, L. K. 2010, ApJ, 722, ApJ, 730, 72 919 Rahman, N., Bolatto, A. D., Xue, R., Wong, T., Leroy, A. K., Glover, S. C. O., Federrath, C., Mac Low, M.-M., & Klessen, Walter, F., Bigiel, F., Rosolowsky, E., Fisher, D. B., Vogel, R. S. 2010, MNRAS, 404, 2 S. N., Blitz, L., West, A. A., & Ott, J. 2012, ApJ, 745, 183 Helfer, T. T., Thornley, M. D., Regan, M. W., Wong, T., Sheth, Regan, M. W., Thornley, M. D., Helfer, T. T., Sheth, K., Wong, K., Vogel, S. N., Blitz, L., & Bock, D. C.-J. 2003, ApJS, 145, T., Vogel, S. N., Blitz, L., & Bock, D. C.-J. 2001, ApJ, 561, 218 259 Roman-Duval, J., Jackson, J. M., Heyer, M., Rathborne, J., & Kawamura, A., Mizuno, Y., Minamidani, T., Filipovi´c, M. D., Simon, R. 2010, ApJ, 723, 492 Staveley-Smith, L., Kim, S., Mizuno, N., Onishi, T., Mizuno, Saintonge, A., Kauffmann, G., Kramer, C., Tacconi, L. J., A., & Fukui, Y. 2009, ApJS, 184, 1 Buchbender, C., Catinella, B., Fabello, S., Graci´a-Carpio, J., Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2008, ApJ, 689, Wang, J., Cortese, L., Fu, J., Genzel, R., Giovanelli, R., Guo, 865 Q., Haynes, M. P., Heckman, T. M., Krumholz, M. R., —. 2009, ApJ, 693, 216 Lemonias, J., Li, C., Moran, S., Rodriguez-Fernandez, N., Kuhlen, M., Krumholz, M. R., Madau, P., Smith, B. D., & Wise, Schiminovich, D., Schuster, K., & Sievers, A. 2011, MNRAS, J. 2012, ApJ, 749, 36 415, 32 Lada, C. J., Forbrich, J., Lombardi, M., & Alves, J. F. 2012, ApJ, Sanders, D. B., Mazzarella, J. M., Kim, D.-C., Surace, J. A., & 745, 190 Soifer, B. T. 2003, AJ, 126, 1607 Lagos, C. D. P., Baugh, C. M., Lacey, C. G., Benson, A. J., Kim, Schruba, A., Leroy, A. K., Walter, F., Bigiel, F., Brinks, E., de H.-S., & Power, C. 2011, MNRAS, 418, 1649 Blok, W. J. G., Dumas, G., Kramer, C., Rosolowsky, E., Leitner, S. N. & Kravtsov, A. V. 2011, ApJ, 734, 48 Sandstrom, K., Schuster, K., Usero, A., Weiss, A., & Leroy, A. K., Bolatto, A., Bot, C., Engelbracht, C. W., Gordon, Wiesemeyer, H. 2011, AJ, 142, 37 K., Israel, F. P., Rubio, M., Sandstrom, K., & Stanimirovi´c, S. Sheth, K., Regan, M., Hinz, J. L., et al. 2010, PASP, 122, 1397 2009, ApJ, 702, 352 Tacconi, L. J. & Young, J. S. 1986, ApJ, 308, 600 Leroy, A. K., Lee, C., Schruba, A., Bolatto, A., Hughes, A., Pety, van Zee, L., Salzer, J. J., Haynes, M. P., O’Donoghue, A. A., & J., Sandstrom, K., Schinnerer, E., & Walter, F. 2013a, ApJ, Balonek, T. J. 1998, AJ, 116, 2805 769, L12 Wong, T. & Blitz, L. 2002, ApJ, 569, 157 Leroy, A. K., Walter, F., Brinks, E., Bigiel, F., de Blok, W. J. G., Zheng, Z., Meurer, G. R., Heckman, T. M., Thilker, D. A., & Madore, B., & Thornley, M. D. 2008, AJ, 136, 2782 Zwaan, M. A. 2013, MNRAS, 434, 3389