MOLECULAR GAS CONTENTS AND SCALING RELATIONS FOR MASSIVE PASSIVE GALAXIESAT INTERMEDIATE REDSHIFTS FROM THE LEGA-C SURVEY

Molecular Gas Contents and Scaling Relations for Massive Passive Galaxies
At Intermediate Redshifts From the Lega-C Survey

[    [    [    [    [    [    [    [    [    [    [    [    [    [    [    [
Abstract

A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to , suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the local universe. Here we present deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2–1) emission in 8 massive () galaxies at selected to lie a factor of 3–10 below the star-forming sequence at this redshift, drawn from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. We significantly detect half the sample, finding molecular gas fractions 0.1. We show that the molecular and stellar rotational axes are broadly consistent, arguing that the molecular gas was not accreted after the galaxies became quiescent. We find that scaling relations extrapolated from the star-forming population over-predict both the gas fraction and gas depletion time for passive objects, suggesting the existence of either a break or large increase in scatter in these relations at low specific star formation rate. Finally, we show that the gas fractions of the passive galaxies we have observed at intermediate redshifts are naturally consistent with evolution into local massive early-type galaxies by continued low-level star formation, with no need for further gas accretion or dynamical stabilization of the gas reservoirs in the intervening 6 billion years.

galaxies: evolution — galaxies: ISM — galaxies: high-redshift
Corresponding author: Justin Spilkerspilkerj@gmail.com

0000-0003-3256-5615]Justin Spilker \move@AU\move@AF\@affiliationDepartment of Astronomy, University of Texas at Austin, 2515 Speedway, Stop C1400, Austin, TX 78712, USA

0000-0001-5063-8254]Rachel Bezanson \move@AU\move@AF\@affiliationDepartment of Physics and Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA

0000-0001-6371-6274]Ivana Barišić \move@AU\move@AF\@affiliationMax-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany

0000-0002-5564-9873]Eric Bell \move@AU\move@AF\@affiliationDepartment of Astronomy, University of Michigan, 1085 South University Avenue, Ann Arbor, MI 48109-1107, USA

0000-0003-3021-8564]Claudia del P. Lagos \move@AU\move@AF\@affiliationInternational Centre for Radio Astronomy Research, M468, University of Western Australia, 35 Stirling Hwy, Crawley, WA 6009, Australia \move@AU\move@AF\@affiliationARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), 44 Rosehill Street Redfern, NSW 2016, Australia

0000-0003-0695-4414]Michael Maseda \move@AU\move@AF\@affiliationLeiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands

0000-0002-9330-9108]Adam Muzzin \move@AU\move@AF\@affiliationDepartment of Physics and Astronomy, York University, 4700 Keele St., Toronto, Ontario, M3J 1P3, Canada

0000-0003-4196-0617]Camilla Pacifici \move@AU\move@AF\@affiliationAstrophysics Science Division, Goddard Space Flight Center, Code 665, Greenbelt, MD 20771, USA \move@AU\move@AF\@affiliationSpace Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

0000-0001-8823-4845]David Sobral \move@AU\move@AF\@affiliationPhysics Department, Lancaster University, Lancaster LA1 4YB, UK

0000-0001-5937-4590]Caroline Straatman \move@AU\move@AF\@affiliationSterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium

0000-0002-5027-0135]Arjen van der Wel \move@AU\move@AF\@affiliationSterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium \move@AU\move@AF\@affiliationMax-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany

0000-0002-8282-9888]Pieter van Dokkum \move@AU\move@AF\@affiliationAstronomy Department, Yale University, New Haven, CT 06511, USA

0000-0001-6065-7483]Benjamin Weiner \move@AU\move@AF\@affiliationSteward Observatory, University of Arizona, 933 North Cherry Ave., Tucson, AZ 85721, USA

0000-0001-7160-3632]Katherine Whitaker \move@AU\move@AF\@affiliationDepartment of Physics, University of Connecticut, 2152 Hillside Road, Unit 3046, Storrs, CT 06269, USA

0000-0003-2919-7495]Christina C. Williams \@AF@joinNSF Fellow \move@AU\move@AF\@affiliationSteward Observatory, University of Arizona, 933 North Cherry Ave., Tucson, AZ 85721, USA

0000-0002-9665-0440]Po-Feng Wu \move@AU\move@AF\@affiliationMax-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany

1 Introduction

The processes by which galaxies grow and evolve are intimately linked to the accretion and conversion of gas into stars. In particular, because stars form from molecular gas (e.g., schruba11), the heating, cooling, and transport of gas from outside and within galaxies play a large role in determining how efficiently a galaxy can form stars, and the overall mass of stars that can be formed. By and large, galaxies form stars across cosmic time in equilibrium with the supply of fresh gas from accretion and mergers; the interplay between gas accretion, outflows, star formation, and mergers naturally regulates the growth of galaxies (e.g., dave11; dave12; lilly13; peng14).

Surveys of increasingly large numbers of galaxies have shown that the majority of galaxies exhibit a relatively tight and nearly-linear relationship between the current star formation rate (SFR) and the mass of stars already formed (; e.g., noeske07; franx08; whitaker12; whitaker14; speagle14; schreiber16). The intrinsic scatter in the relationship is  dex, and objects that lie near it are generally considered to be ‘normal’ galaxies. The normalization of the ‘star-forming sequence’ increases rapidly with redshift, implying much more rapid galaxy growth in the early universe as compared to today, or an overall increase in the SFR density function (sobral14). At the most massive end, however, an ever-larger fraction of galaxies exhibit markedly depressed SFRs (or specific SFR, sSFR  SFR). This transition occurs near the break in the stellar mass function, , and does not appear to evolve significantly with redshift (e.g., peng10b, though see also gavazzi15; tomczak16). Star formation in these massive galaxies appears to be efficiently shut off (‘quenched’) as they transition to the red sequence, but the physical mechanisms responsible for this quenching are still unclear.

After a decade of extensive observational investment, the relationship between gas supply and the growth of star-forming galaxies has become more clear, both at low redshift (e.g., saintonge11; saintonge17; bothwell14) and in the distant universe (e.g., tacconi13; magnelli14; papovich16; scoville16). For galaxies with SFRs near and above the star-forming sequence, scaling relations have been derived relating the properties of the molecular interstellar medium (ISM) with other galaxy properties (e.g., genzel15; scoville17; tacconi18), with remarkably good agreement between various tracers of the molecular gas. In general, the molecular gas fraction increases rapidly with redshift, is elevated for galaxies well above the star-forming sequence, and shows a weak decline with increasing . The gas depletion time, , a measure of how long a galaxy could continue to form stars at its current rate before exhausting its gas supply, shows a much weaker evolution with redshift, is shorter for galaxies above the star-forming sequence, and is either constant or weakly increases towards high . Typical values for Milky Way-mass star-forming galaxies at are few  and  Gyr, and at , and  Gyr (tacconi18).

The recent advent of large samples of ‘normal’ star-forming galaxies with molecular gas measurements has been particularly useful at high redshifts, providing extremely valuable reference samples for other galaxies that may be less ‘normal.’ Of particular relevance for the buildup of the quiescent galaxy population, molecular gas observations have now also targeted smaller samples of galaxies thought to be actively quenching star formation for various reasons (e.g., geach13; spilker16b; popping17) or that show spectral signatures of quenching in the past 1 Gyr (suess17). These studies indicate that, while the suppression of star formation does not require the complete removal or depletion of the molecular gas, the quenching processes do appear to lower at fixed mass compared to the reference samples in most cases.

Because of the relationship between and SFR, however, much less is known about the molecular ISM in galaxies well below the star-forming sequence; the expected molecular masses require very sensitive observations and substantial integration times even with the supreme sensitivity of ALMA. Observations of local massive and passive early-type galaxies, for example, reveal gas fractions an imposing 1–2 orders of magnitude lower than their star-forming counterparts (young11; davis13; davis16). Because these local galaxies show few signs of recent star formation for the past many Gyr, however, it is not clear to what extent inferences about the suppression of star formation in these galaxies also apply in the distant universe, closer to the epoch at which galaxies first became quiescent.

The main focus of our work here is to determine whether the scaling relations derived from observations of star-forming galaxies and cosmological simulations can accurately predict or be extrapolated down to passive galaxies at intermediate redshift. We present ALMA observations of a sample of 8 massive () galaxies at from the LEGA-C survey selected to lie a factor of 3–10 below the star-forming sequence at this redshift. We observed the CO(2–1) transition, a tracer of the molecular ISM. In Section 2, we describe the parent LEGA-C sample, our ALMA observations, and our measurements of . Section 3 provides a broad overview of our basic results, and we compare the stellar and molecular dynamics of our detected galaxies. In Section 4, we compare the gas fractions and depletion times derived for our sample with observationally-based scaling relations and with the EAGLE cosmological simulation. We discuss the implications our observations have for the understanding of galaxy quenching in Section 5, and conclude in Section 6. Throughout, we assume a flat CDM cosmology with and  Mpc (planck15).

2 Data and Analysis

2.1 LEGA-C and Selection of Quiescent Galaxies

We selected galaxies for CO(2–1) observations from the LEGA-C survey of galaxies. The survey is described in detail by vanderwel16. Briefly, galaxies were selected based on -band magnitude from the UltraVISTA catalog described in muzzin13a. The LEGA-C survey consists of 3200 galaxies observed with 20 hour integrations using the VIMOS spectrograph on the Very Large Telescope. The LEGA-C spectra yield high signal-to-noise detections of the stellar absorption features and continuum, allowing determinations of the age of the stellar populations, metallicities, and stellar velocity dispersions. Because the LEGA-C survey targets the COSMOS extragalactic survey field, full panchromatic spectral energy distribution (SED) information is available, along with morphological information from Hubble Space Telescope observations. The galaxies studied in this work were selected from the first data release (DR1) catalog for which VIMOS observations were completed by early 2016, consisting of 644 objects in the primary sample with good spectra. This early sample represents a subset of the full survey area and is not biased with respect to the full LEGA-C sample in terms of, e.g., K-band magnitude. Note that, while the galaxies were selected from the initial DR1 catalog, the figures and table in this work use values from the updated DR2 catalog (1989 galaxies; Straatman et al., in prep.). Stellar masses of the LEGA-C sample have been measured by fitting the photometric SED using FAST (kriek09) assuming a chabrier03 initial mass function.

2.1.1 SFR Estimation

The primary SFR estimates for the LEGA-C sample come from modeling of the spectral energy distribution of our target galaxies, from the ultraviolet to the mid-infrared. SFRs include the unobscured and obscured components based on UV and IR (24 ) fluxes. The low resolution of the Spitzer/MIPS 24  imaging requires the use of a deblending procedure to assign the measured 24  flux to K-band detected galaxies (muzzin13a). In our final ALMA sample, described further below, two objects (IDs 130284 and 132776) are near bright 24  sources, where the deblending is potentially unreliable. For all objects in the ALMA sample, the inclusion of the 24  fluxes increases the inferred SFR by a factor of 2.5, on average. Unsurprisingly, none of these objects with low SFR are detected in Herschel Space Observatory imaging of the COSMOS field.

Several authors have noted that the IR luminosity (or observed-frame 24  luminosity, generally the only available tracer of dust emission at low SFR and high redshift) can overestimate the obscured SFR of quiescent galaxies, in some cases quite severely (e.g., salim09; hayward14; utomo14; man16). This can be due to several effects, including the fact that the IR luminosity is a long-lived tracer compared to the instantaneous SFR, and 24  emission can be boosted by dust heating unrelated to star formation, including circumstellar dust heated by intermediate-age AGB stars, extended cirrus dust heated by old stellar populations, or weak nuclear activity. The latter three result in a higher /SFR ratio, or equivalently only a fraction of the observed should be considered in the calculation of the SFR.

Of these options, weak nuclear activity is unlikely, with the 24  emission expected from stacking of X-ray images falling three orders of magnitude below the observed emission for intermediate-redshift quiescent galaxies (fumagalli14). Dust heating from AGB stars and cirrus dust are potentially more relevant. fumagalli14 used a stacking analysis of MIPS/24  images of quiescent galaxies at , and found that the SFR determined from these images could be overestimated by an order of magnitude. However, there is reason to believe the situation is not so dire for our own selection. fumagalli14 focused on galaxies much more quiescent than ours, 20–40 below the star-forming sequence, compared to our own selection of 3–10 below. The objects in our sample do not reach such low sSFR even if none of the IR emission is related to star formation and only the SFR based on the rest-UV is considered. The SFR overestimation becomes less severe at higher sSFR because true obscured star formation rapidly outshines the lower-level IR emission unrelated to SF. For the sSFRs typical of our sample, the work by fumagalli14 indicates that the SFRs should be overestimated by less than a factor of .

2.1.2 ALMA Selection

We selected galaxies from the LEGA-C sample with spectroscopic redshift , for which the CO(2–1) line is observable with the ALMA Band 4 receivers. We selected massive galaxies, with and SFR 3–10 below the star-forming sequence as defined by whitaker12. Given the intrinsic scatter in this sequence ( dex), this selects objects with SFRs lower than 90% of all galaxies at this redshift. We further required that SFR  2.5 /yr, intended to restrict the sample to galaxies that could plausibly be detected in CO in a reasonable integration time with ALMA. This resulted in a sample of 65 galaxies, from which we selected 8 for observation with ALMA, listed in Table 2.3. We aimed to include galaxies with a variety of morphologies, both obvious spiral galaxies and bulge-dominated elliptical galaxies. Two of the eight galaxies (IDs 110509 and 169076) are identified as radio-loud AGN (barisic17); both lack optical emission lines and are classified as low-excitation radio galaxies. The sensitivity of the best-available radio imaging of the COSMOS field (smolcic17) is not sufficient to detect radio continuum emission related to star formation for our sample (approximate 3 limiting SFR18 /yr). The VIMOS spectra of our 8 galaxies are shown in Figure 2.1.2.

For two objects in the ALMA sample, the LEGA-C spectra clearly detect the [OII] 3727 doublet, which can also be used as an SFR indicator, although one highly susceptible to dust extinction and low-level nuclear activity. Converting this line flux to an SFR yields values less than the UVIR SFRs, as expected, but within 2 the UV-only SFR. The remaining objects show little evidence for nebular emission lines also used as SFR indicators (e.g., H), with any emission generally overwhelmed by the strong absorption features.

In Figure 2.1.2 we show the selection of the ALMA sample in the SFR– plane with respect to the full LEGA-C sample and all other LEGA-C objects that meet our selection criteria. We also show the size-mass relation for star-forming and quiescent galaxies, with the best-fit relations determined by vanderwel14. Because we chose objects with a range of morphologies, it is unsurprising that the ALMA sample includes objects that both follow and do not follow the expected size-mass relation for quiescent galaxies.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref The selection of the ALMA-observed massive, passive sample with respect to the full LEGA-C sample. In both panels, the LEGA-C sample galaxies are color-coded blue (red) if they are classified as star-forming (quiescent) in rest-frame space (see Figure 2.1.2). Red circles are all LEGA-C galaxies that meet our selection criteria, while the ALMA-observed objects are shown with large red diamonds; CO-undetected objects are also marked with a black ‘’. Top: The ALMA sample was primarily selected based on stellar mass and SFR; see text for details. The blue line shows the star-forming sequence at from whitaker12, and the blue shaded region encompasses SFRs a factor of 3 above and below the relation. The red shaded region shows our selection box of massive and passive galaxies. Bottom: Size-mass relation for the LEGA-C and ALMA samples. Blue and red lines and regions show the size-mass relations for star-forming and quiescent galaxies at from vanderwel14.

Figure 2.1.2 shows two other relevant diagrams that give the ALMA sample additional context. The top panel shows the ALMA sample in color-color space, with the division between star-forming and quiescent objects from muzzin13a also shown. The boundary between these two populations in colors is not a definitive one, with varying definitions found in the literature as well as an expectation that galaxies must transition from one region to the other over time. Although this diagram was not used in the selection of the ALMA sample, it is reassuring that 6 of the 8 targets lie within (given the uncertainties) the region of the color-color diagram generally occupied by quiescent galaxies. Of the two remaining objects, one (ID 74512) lies in the region generally occupied by star-forming galaxies; the other (ID 138718) lies in the region generally occupied by dusty galaxies.

The lower panel of Figure 2.1.2 shows the distribution of D4000 index against the equivalent width of H for the LEGA-C and ALMA-observed samples (wu18; an emission line template has been subtracted from the spectra). This diagram is useful as a proxy for the stellar age, older galaxies showing lower H equivalent width and higher D4000. These quantities are also correlated with , indicating older stellar populations in massive galaxies (e.g., kauffmann03; maltby16). Based on this diagram, we infer typical stellar ages for the ALMA-observed galaxies of 1–3 Gyr, depending on the assumed star formation history; a similar result is obtained through fitting to the full LEGA-C spectra (Chauke et al., submitted). In other words, the objects we observed with ALMA have been passive for quite some time. They show little evidence for significant recent star formation, unlike galaxies selected with ‘post-starburst’ criteria that emphasize the recent or impending quenching of star formation (e.g., tremonti07; sell14; suess17).

Finally, we examined the local galaxy environments of our target objects using the estimates of darvish16. These estimates use photometric redshifts because the spectroscopic completeness in the COSMOS field is low. darvish16 note that at , a comparison of spectroscopic and photometric redshifts shows that the redshift uncertainties are generally small, with a dispersion , sufficiently small that the use of photometric redshifts does not wash out line-of-sight galaxy structures. For the LEGA-C sample as a whole, the median overdensity is , while for the ALMA-observed sources, the median is 0.16, with large uncertainties on both values. The ALMA targets are not in systematically overdense environments compared to the full LEGA-C sample or the field as a whole.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Symbols are plotted as in Figure 2.1.2. Top: Rest-frame color-color diagram for the LEGA-C galaxies and the ALMA sample, with the division between star-forming and quiescent galaxies of muzzin13a shown with the black line. The gray shaded band around this line represents both the differences in divisions found in the literature and an expectation that galaxies must transition from one region to the other over a period of time. Although this diagram was not used in the selection of the ALMA sample, six of the eight targets lie within (given the uncertainties) the quiescent region of the diagram at upper left. Bottom: The D4000 index against the H equivalent width, a proxy for the age of the stellar populations, with older and more massive systems located towards the lower right. We infer typical stellar ages of 1–3 Gyr for the ALMA sample; these objects are not recently quenched. For two objects (IDs 110509 and 169076), the LEGA-C spectra do not extend sufficiently blueward to measure D4000, shown arbitrarily at D4000 with arrows in each direction.

2.2 ALMA Observations

Basic details of our ALMA observations and target galaxies are summarized in Table 2.3. ALMA observations of our sample of 8 galaxies were carried out in project 2016.1.00790.S (PI: Spilker) in separate observing sessions from 17 January to 12 March 2017 using the Band 4 (2 mm) receivers (asayama14). All observations were conducted with the array in a compact configuration, with 40 antennas separated by maximum baselines ranging from 272–330 m. Observing sessions varied between 57–80 minutes in duration, with 37–49 minutes spent on-source per target. The precipitable water vapor levels varied between 2.7–4.5 mm, resulting in typical system temperatures of 70–95 K. Quasars J1058+0133 or J0854+2006 served as bandpass calibrators, while the quasar J0948+0022 was observed for complex gain calibration for all sources. The absolute flux scale was determined using observations of Ganymede or one of J0854+2006 or J1058+0133, both of which are monitored regularly by ALMA.

The correlator was configured to observe the CO(2–1) line at the known redshift of each target with one baseband with 7.812 MHz channelization () after correlator pre-averaging by a factor of 8. Three further basebands with 1.875 GHz usable bandwidth each were placed to higher frequencies for continuum observations. All data were reduced using the standard ALMA pipeline, with manual inspection of the quality of the reduction.

For each object, we produce continuum images and CO(2–1) spectral cubes with various frequency channelization. All data were imaged with a natural weighting of the visibilities, which maximizes sensitivity to faint emission in exchange for slightly lower spatial resolution. The effective spatial resolution of the data imaged in this way ranges from 1.9–2.4 (13–17 kpc). No 2 mm continuum emission is detected in any target. We extract integrated spectra for all targets by fitting point source models to the visibility data directly, averaging 6 or 12 channels to create spectra with velocity resolution 100–200 . When averaging over the full line profiles, we see evidence that most of the detected sources are marginally spatially resolved (e.g., by comparing the peak pixel values with spatially-integrated fluxes), but this effect is negligible in narrower channels.

Finally, because all of our targets are located in the COSMOS extragalactic deep field, we checked the ALMA archive to determine if any other observations of our targets were publicly available or fell within the footprint of other projects. No other observations were found.

2.3 Molecular Gas Masses

\H@refstepcounter table \hyper@makecurrenttable\hb@xt@ Table 0. \Hy@raisedlink\hyper@@anchor\@currentHrefALMA-Observed LEGA-C Sample Target Properties

LEGA-C ID Right Ascension Declination z SFR S /yr Jy/beam Jy  74512 0.7330 10.96 6.3 132 0.16 0.04 9.82 0.13 110509 0.6671 11.00 6.5 160 0.24 0.04 9.92 0.07 130284 0.6017 10.96 6.8 151 0.36 0.04 10.00 0.06 132776 0.7500 10.98 7.9 163 0.33 0.07 10.16 0.11 138718 0.6558 11.13 5.6 188 0.21 9.84 169076 0.6772 11.49 5.1 256 0.23 9.91 210210 0.6544 11.38 3.6 212 0.21 9.84 211409 0.7140 11.13 6.6 188 0.13 9.72 Stacked Non-detections 0.6754 11.31 5.3 106 0.093 9.51

ALMA rms sensitivity in 100  channels, naturally-weighted images Note—LEGA-C ID numbers are the same as in the UltraVISTA catalog of muzzin13a. Stellar masses are determined by fitting to multiwavelength photometry using FAST. SFRs are based on a weighted sum of UV and IR (24 ) fluxes. Integrated CO(2–1) line fluxes are converted to molecular gas masses under the assumptions described in Section 2.3. Upper limits for non-detections are 3, and molecular gas masses can be rescaled under different assumptions as .

\twocolumngrid

Spectra and integrated line images of each target are shown in Figure 2.2. We clearly detect four of the eight targets in CO(2–1) emission. Two of the four detections show double-peaked line profiles, usually indicative of rotating disks; the other two sources are centrally-peaked. We discuss kinematics in more detail in Section 3.2. Integrated line fluxes are determined by fitting either one or two Gaussian profiles, also shown in Figure 2.2. For the undetected sources, we estimate upper limits on the CO emission by determining the noise in channels 800  wide, which would fully encompass all of the line emission in all detected sources. The upper limits on the integrated line flux are proportional to , where is the velocity interval over which the spectrum is integrated, so using very wide channels results in conservative upper limits. The integrated line fluxes are given in Table 2.3.

We convert the observed CO luminosities and upper limits to estimates of the molecular gas masses under standard assumptions about the CO excitation and the CO–H conversion factor . The effects of the unknown CO line excitation are minimal, as the CO(2–1) transition we have observed is very close to the ground state CO(1–0) line typically used for molecular gas estimation. Observations of the Milky Way, nearby quiescent and star-forming galaxies, and high-redshift star-forming galaxies ubiquitously show that a line ratio (in brightness temperature units) encompasses the plausible expected range (fixsen99; combes07; dannerbauer09; young11; spilker14; saintonge17). In this work, we assume .

The subsequent conversion between CO luminosity and molecular gas mass is also uncertain (for a recent review, see bolatto13). The conversion factor is known to vary with the gas metallicity (e.g., leroy11), which affects the formation and destruction of CO molecules. The metallicities (either gas-phase or stellar) of our targets have not been measured, but are expected to be solar or near-solar based on the mass-metallicity relation. For example, gallazzi14 predict for quiescent galaxies in our mass range, with a scatter of . We therefore expect only minor variations in due to metallicity effects. The CO-H conversion factor also depends on the gas conditions and kinematics, which affect the optical depth of the CO transitions through radiative trapping. This is chiefly relevant for mergers and other high-SFR systems, in which increased gas turbulence and/or bulk motions lower the effective CO optical depth and also (e.g., narayanan12; spilker15). These effects are also expected to be minor for our targeted objects, which have low SFRs, no signs of interaction, and evidence for disk-like rotation in many cases (see Section 3.2).

In this work, we adopt a ‘Milky Way-like’ value,  (K   pc) (e.g, solomon87; bolatto13; sandstrom13). This value agrees with standard dust-based methods we describe further as part of our stacking analysis below. While this choice is justified for the reasons already mentioned, it does still carry significant systematic uncertainty, likely of order 50%. While this may result in adjustments to the absolute values of the gas masses we derive, the relative values are robust, and it does not affect the overall trends we find. Of relevance to our subsequent discussion in Section 4, the true gas masses of our sample are unlikely to be significantly larger than the values we infer. The molecular gas masses we derive can easily be rescaled using different assumptions, as . When comparing to other samples and galaxies observed in CO by other authors, we also adjust their derived gas masses to match our adopted value of . As we are interested in normal star-forming and quiescent galaxies, this adjustment is minor, no larger than % for the comparison samples.

2.4 ALMA Stacking Analysis

With CO non-detections constituting half the observed sample, it is worth considering whether these sources are detected on average through a stacking analysis. This would imply that somewhat deeper integrations would have been necessary to detect the objects individually. Similarly, although we did not detect (nor expect to detect) dust continuum emission in any individual target, we expect that it should be detectable in a stacked continuum image, if the assumptions made about in the previous section are correct.

2.4.1 2 mm Continuum Stack

Because all targets were observed in similar array configurations with similar synthesized beam sizes, for similar durations, and at similar sky frequencies, we opt to simply average together inverted images of the targets rather than perform this analysis in the Fourier plane. The similarity in observed frequency is of particular importance for the continuum stacked image, minimizing the effects of the steep thermal dust spectral index. We create a stacked 2 mm continuum image simply by averaging together the signal-free continuum maps of each source. In the stacked image, a faint source with is detected at the center of the field. By subdividing the stack into sources either detected or undetected in CO(2–1), it is apparent that the signal seen in the full stack is due entirely to the sources individually detected in CO. The CO-detected sources are detected in the 2 mm continuum stack with , while no continuum emission is detected in the stack of CO non-detections, with a 3 upper limit .

We can use the detection of 2 mm continuum emission in the stack of CO-detected sources to cross-check our assumptions about the CO-H conversion factor described in Section 2.3. The average molecular gas mass of the detected objects is . Assuming a typical dust emissivity,  cm g (e.g., dunne00), where 235 GHz is the rest-frame continuum frequency, a dust temperature  K, and a gas-to-dust ratio typical of solar metallicity systems, (e.g., sandstrom13), we find an average molecular gas mass for the CO-detected sources of . This value is in excellent agreement with our CO-based estimates. The same assumptions applied to the 2 mm upper limit from the CO-undetected sources yields a 3 upper limit of   , unsurprisingly in agreement with the CO-based upper limits for individual sources. We note that the dust-based estimate of is also uncertain by about a factor of two due to uncertainties in the dust emissivity and mass-weighted temperature (e.g., draine07). This independent check indicates that our assumptions about are reasonable.

2.4.2 CO Spectral Stack

We create a stacked CO image cube in similar fashion to the simpler stacked continuum maps. In this case, however, we rely on the fact that all targets (detected and undetected) have precisely-known redshifts from the LEGA-C spectra. The redshifts are accurate to , much less than the typical line widths seen in CO or the stellar absorption features in the LEGA-C spectra. We also do not see significant velocity offsets between the CO emission and stellar absorption features in the CO-detected sources. We stack the spectra of individual sources using image cubes with velocity resolution ranging from 50–800 ; the choice of channel width does not affect our conclusions. We have also verified that this procedure accurately recovers the sample-average line flux by stacking only the sources individually detected in CO.

No CO(2–1) emission is detected in the stack of individually-undetected sources, regardless of the velocity resolution used in the stack. As with the individually-undetected sources, we place an upper limit on the CO emission in this stack using a single 800  wide channel. This results in a 3 upper limit on the CO luminosity  K   pc. Under the same assumptions about the CO excitation and as before, this is equivalent to . Given the average stellar mass of the CO-undetected sources, , the resulting upper limit on the gas fraction of these passive sources is .

3 Results

3.1 Basic Gas Properties

The derived properties of our sample of passive galaxies are given in Table 2.3. Figure 3.1 plots a number of basic correlations from our data, comparing galaxy molecular masses, gas fractions, and depletion times with their stellar masses, SFRs, and specific SFRs. We draw comparison samples from the literature. We include the xCOLDGASS sample (saintonge17), a large sample of galaxies in the local universe that contains both star-forming and passive objects over several orders of magnitude in stellar mass. We also include higher-redshift star-forming galaxies at compiled from the PHIBSS survey (tacconi13), consisting of galaxies, supplemented by the observations of papovich16, which extend to lower mass (). Finally, we include the recent ALMA observations by suess17 of two galaxies selected to have post-starburst-like optical spectra, , and little ongoing star formation. With similar stellar masses to our own sample, these galaxies represent objects that quenched recently, as opposed to the several-Gyr-old stellar populations in our sample. Gas masses from all samples have been renormalized to , but we preserve the original authors’ assumptions about the CO excitation, applicable to the samples.

A number of well-known strong correlations are immediately apparent in Figure 3.1, including the spatially-integrated Schmidt-Kennicutt star formation relation between SFR and (schmidt59; kennicutt98b) and the declining depletion times observed for galaxies with high sSFR (e.g., saintonge17). Because all of these quantities evolve with redshift, it is generally more instructive to interpret the variations seen in Figure 3.1 after removing or otherwise accounting for the redshift evolution. We perform this exercise in Sections 4.1 and 4.2.

At the most basic level, we observe that massive, passive galaxies contain of molecular gas, and may contain significantly less, depending on the true masses of the undetected sources. This is nevertheless over an order of magnitude higher (also ) than observed for early-type and other massive, low-sSFR galaxies in the local universe (e.g., combes07; young11; davis16). Unlike the local samples, however, we see no evidence for increased gas depletion times that might indicate that the molecular material has been stabilized against collapse; instead, the 1–2 Gyr depletion times are typical of, or even shorter than, galaxies near the star-forming sequence.

In terms of basic molecular gas properties, the galaxies we have observed are not so dissimilar from ‘normal’ star-forming galaxies in the local universe. While they were selected to be massive galaxies, the absolute values of , , and in the ‘passive’ galaxies are very similar to star-forming galaxies today. In some sense, the galaxies we have observed can be considered ahead of their time – they formed the bulk of their stellar mass at early times, and reached an evolutionary stage similar to that currently experienced by galaxies in the local universe  Gyr early.

It is rather curious that the detected sources are very clearly detected, at the level, while the others are not detected even in a stacked spectrum. Why were some objects in our sample detected, and others not? With only 8 target galaxies in total, we cannot draw strong conclusions on this point. The four non-detections are both the highest- and lowest sSFR galaxies in the sample, although the absolute SFRs are not significantly different from the full sample. The non-detections tend to lie towards higher stellar surface density than the detections (Figure 2.1.2), but are not obviously disparate in morphology or other structural parameters (e.g., Sérsic index). They also have the lowest H equivalent widths in the sample and the reddest colors (Figure 2.1.2), but do not show large differences in best-fit stellar age, presence or lack of radio AGN, or presence or lack of emission lines in the LEGA-C spectra. Finally, the undetected sources are in regions more overdense than the detected sources by  dex. There is no statistical difference in the relative overdensity between detected and undetected sources. We stress that a much larger sample is required to begin to understand which of these properties, if any, are good or useful predictors of CO line luminosity.

We see an interesting contrast between our massive, passive galaxies, and the post-starburst galaxies at similar redshift and observed by suess17. While SFRs for the post-starburst sample are measured from the extinction-susceptible [OII] doublet, and thus may be underestimated, the absolute quantity of molecular gas in these two objects still serves as a useful reference. While one post-starburst has a similar molecular mass to the objects we have studied, the other has higher than any of our objects, and more than an order of magnitude higher than the stack of CO-undetected sources in our sample. If the SFR of this object is confirmed to be as low as inferred, it would indicate a large diversity in the gas masses, fractions, and depletion times among massive galaxies at intermediate redshifts.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Summary of basic results derived from our data, with comparison samples detailed in the text. From top to bottom, each row plots the molecular gas mass , gas fraction , and gas depletion time ; from left to right each column shows the stellar mass , star formation rate SFR, and specific SFR. The comparison samples are drawn from tacconi13; papovich16; saintonge17; suess17. All upper limits are 3. For clarity of presentation, we do not show upper limits for the xCOLDGASS sample in the bottom () row.

3.2 Molecular and Stellar Dynamics

The dynamical state of the molecular gas in passive galaxies can be a powerful probe of its origin and the processes that led to the cessation of star formation. In early-type galaxies, for example, approximately one in four contain detectable amounts of molecular gas, much more common among the fast rotator galaxies than slow rotator population (young11). Among these objects, about 40% show signs that the molecular gas discs are significantly misaligned (more than 30) with respect to the stellar rotation (alatalo13; davis13). This has been interpreted as evidence that the molecular gas was gained from external sources after the galaxies quenched star formation, for example from the accretion of gas-rich satellites or cold streams from the intergalactic medium. Misaligned rotation axes can also arise during major mergers if gas makes up a significant fraction of the total mass (lagos18). The dynamics of the molecular gas can thus provide insight into the origin of this material and its relationship to the physical processes that quenched star formation.

From the integrated spectra alone, two of our CO-detected objects show double-horned line profiles characteristic of rotating gas discs, while the other two show centrally-peaked line profiles. This is similar to the results from local early-type galaxies, for which % show evidence of rotation from the CO line profiles alone (young11). The remaining objects may also contain rotating gas discs, depending on the spatial distribution of the molecular gas and the inclination with respect to the line of sight.

To investigate the velocity fields of the CO-detected objects in more detail, we re-imaged the ALMA data, splitting the total CO emission from each object into two velocity channels covering the red- and blue-shifted emission. The results are shown in Figure 3.2 (top row). While narrower velocity bins can also be imaged, a single blue and red channel yield the highest signal-to-noise for presentation purposes; we have verified that the velocity gradients apparent in Figure 3.2 are consistent with imaging the data in narrower velocity channels. We fit the red and blue channels of each source in both the visibility and image domains with point source models, finding that the centroids of each velocity component can be determined to 0.3, on average. Given the uncertainties, significant velocity gradients are apparent in three of the four detected sources (IDs 74512, 110509, and 130284). The remaining source (ID 132776) may also show a velocity gradient in CO, but the spatial separation between blueshifted and redshifted velocities is not significant at the signal-to-noise of the current data. Our finding that a high fraction of passive galaxies show measurable rotation mirrors the results for higher-SFR massive galaxies at equivalent and higher redshifts (e.g., forsterschreiber09; tacconi13; wuyts16).

A quantitative comparison between the molecular gas and stellar dynamics is challenging given the data in hand. The LEGA-C spectra provide a high signal-to-noise measurement of the stellar rotation curve (bezanson18), but all VIMOS slits for the ALMA sample were 1 wide and oriented north-south, and so are effectively randomly aligned with respect to the galaxy major axes. Additionally, the typical 1 seeing during the VLT observations leads the rotation curve velocities to be correlated on scales of a kpc. On the other hand, the ALMA CO observations have both lower spatial resolution (2–2.5) and lower signal-to-noise compared to the LEGA-C spectra, but contain full two-dimensional spatial information and independent velocity channelization.

Given the differences between the available data and the modest signal-to-noise of the current CO maps, we choose not to make a detailed comparison between the CO and stellar dynamics. Instead, we address a simpler question: are the LEGA-C stellar rotation curves consistent with the observed CO velocity gradients, given the misalignment between the projected galaxy rotation axes and the north-south oriented VIMOS slits? In Figure 3.2 (bottom row), we show the velocity gradients seen in the CO data, where the offset between red and blue halves of the emission is measured along the galaxy rotation axis derived from the ALMA data. We also show the LEGA-C stellar rotation curves, which, as mentioned, are sensitive only to the north-south projected component of the velocity field. Finally, we show the portion of the ALMA velocity gradient also projected onto the north-south orientation of the LEGA-C slits.

In three of four cases, we find excellent agreement between the LEGA-C rotation curve and the north-south projection of the CO velocity gradient. The final object, ID 74512, is the most compact galaxy in the sample and does not obviously show signs of rotation in the LEGA-C spectrum. This may be at least in part due to the compact size of the galaxy in comparison to the VIMOS slit and the low signal-to-noise of the stellar rotation curve. It is also possible that the stellar component of this galaxy truly has little rotation, which would imply that the angular momentum axes of the stars and molecular gas are misaligned. Given the large uncertainties involved, however, this scenario is neither supported nor unsupported by the data in hand.

In summary, the rotational axes of the stellar and molecular components of the CO-detected objects are consistent in at least three of four cases. This argues against an external origin for the molecular gas in passive galaxies (either recently-accreted cold streams or gas-rich mergers), because these processes should commonly result in misaligned stellar and molecular rotation. Instead, the molecular gas in our sample is probably left over from the formation epoch of the bulk of the galaxies’ stellar mass or replenished directly by stellar mass loss, resulting in stellar and gas discs with aligned rotational axes.

4 Discussion

4.1 Gas Scaling Relations

The past decade has seen a large investment of single-dish and interferometer time devoted to understanding how the molecular gas properties of galaxies vary with other galaxy properties from (e.g., bothwell14; saintonge17) out to (e.g., tacconi13; genzel15; scoville17; tacconi18). Because of the well-known correlation between and SFR (Figure 3.1), the moderate– to high-redshift samples have focused nearly exclusively on fairly massive star-forming galaxies on or above the star-forming sequence, the objects with the largest gas masses and most easily detectable. With our sample of passive galaxies at intermediate redshifts, we are in a position to determine whether these scaling relations extend to lower sSFR, previously unexplored parameter space.

The above studies provide prescriptions for the variation of molecular gas fraction and depletion time , parameterized in terms of overall redshift evolution, sSFR (generally with respect to the sSFR expected from the star-forming sequence at a given epoch), , and/or galaxy size . Based on samples of hundreds of galaxies from , consensus has emerged that the gas fraction evolves steeply with redshift, , with shallower dependencies on star-forming sequence offset, , and stellar mass, . The evolution of with redshift appears to be slightly less rapid than the normalization of the star-forming sequence itself. This implies that to first order, the higher SFRs observed in typical galaxies at high redshift are simply due to larger gas masses, with perhaps a somewhat higher efficiency of star formation (equivalently, lower ) also needed. The depletion time exhibits an overall smaller dynamic range, and varies less steeply with redshift, , star-forming sequence offset, , and very shallow dependence on stellar mass, . Because  Gyr for star-forming galaxies, with shallow redshift evolution, galaxies must also have had high gas accretion rates in order to reconcile with the evolution of , reaching /yr at (scoville17).

It is not clear whether these scaling relations should or do extend to passive galaxies significantly below the star-forming sequence. On one hand, the scaling relations are very successful over a very wide parameter space, with residual scatter in and of just  dex (tacconi18). On the other hand, the scaling relations have been derived using only star-forming galaxies, and thus may not account for the physical mechanisms that quench galaxies, or the diversity of these mechanisms that may induce increased scatter in the relations.

In Figure 4.1, we compare the gas fractions and depletion times for our observed sample with two empirical scaling relations from the literature.333Both scoville17 and tacconi18 provide scaling relations using the prescription for the star-forming sequence of speagle14, which we also adopt for Figure 4.1. The difference between the speagle14 and whitaker12 formulations explains why two objects in our sample are 10 below the star-forming sequence, lower than our nominal selection threshold. The relations derived by scoville17 are based on ALMA observations of long-wavelength dust continuum emission and assumptions about the dust emissivity. The work by tacconi18, on the other hand, is a meta-analysis of studies that represents the largest overall sample size and dynamic range to date, and also accounts for different normalizations in the major techniques used to measure gas masses. tacconi18 adopt a metallicity-dependent CO-H conversion factor based on the observed mass-metallicity relation, while we simply use . We note that using their prescription for our sample gives a smaller value for , resulting in even lower gas masses, fractions, and shorter depletion times; in other words, the differences seen in Figure 4.1 would be even further accentuated.

The darkly shaded regions in Figure 4.1 show the approximate lowest sSFR objects at contained in the observational samples; at lower sSFR (lightly shaded), the scaling relations are extrapolated. For both and , Figure 4.1 shows that the extrapolation of the tacconi18 scaling relation (which also incorporates the data from scoville17) is a better match to our data, while the scoville17 relations predict significantly higher and longer compared to our observations. The tacconi18 scaling relations perform reasonably well for the CO-detected sources closest to the star-forming sequence. However, both scaling relations perform poorly for the CO-undetected sources in our sample, especially so for the stacked non-detections. Individual non-detections have upper limits on lower than the extrapolations by , while the upper limit on in the stack is a factor of 5 lower than expected. A similar discrepancy applies to , with individual sources showing depletion times shorter by and the stack too short by a factor of 4.

Considering either or , the CO-undetected sources contain significantly less molecular gas than expected. This result holds even if we have overestimated the SFRs in our sample by , because the scaling of with sSFR is very steep, and altering the SFR moves galaxies largely parallel to the scaling relation. At the sSFRs of our sample, we expect the SFRs to be overestimated by at most about a factor of two (fumagalli14, and Section 2.1). It is interesting to note that the four undetected sources are the four furthest below the star-forming sequence, but with a small sample size it is unclear that this is significant (Section 3.1).

It is apparent that the literature scaling relations cannot be extrapolated and lose predictive power at the low sSFRs probed by our sample. Instead, there appears to be either a break in the scaling relations towards lower at sSFRs 4–5 below the star-forming sequence, or a sharp increase in the scatter of the scaling relations, with our small sample size explaining the lack of detections closer to the expected relations. Equivalently, a curvature term in the dependency could be required in the scaling relations in order to bring them into agreement with our observations. Only a larger sample of galaxies well below the star-forming sequence would reveal whether this is necessary.

It is not clear why half of the observed sample deviates so significantly from the extrapolated scaling relations. If this trend is real, the answer likely involves the physics of galaxy quenching, including feedback processes and galactic winds, gas starvation or lowered accretion rates, environmental effects, and/or morphological transformation. One or more of these processes likely acted on the observed galaxies in order to disrupt their equilibrium growth, leaving them with significantly less molecular gas than otherwise expected. We return to this question in Section 5.1.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref The molecular gas fraction (left) and depletion time (right) versus the specific SFR, compared with observationally-derived literature scaling relations from scoville17 (green densely dashed line and region) and tacconi18 (navy dashed line and region), evaluated for . Both literature scaling relations are parameterized in terms of the offset of the sSFR with respect to the expectation from the star-forming sequence, and so the left-hand y-axis is also labeled in these terms. The right-hand y-axis is labeled with the corresponding values of the sSFR, evaluated at . The shaded regions encompass the quoted approximate 1 residual scatter around the relations. The darker shaded bands correspond to the approximate sSFR lower bounds for galaxies at that were used in the derivation of the scaling relations; below these regions the scaling relations are extrapolations (lightly shaded). In the right-hand panel, a grey arrow indicates the shift if the SFRs for our galaxies have been overestimated by a factor of 2; the shift is largely parallel to the scaling relations. In general, the CO-detected objects in our sample are in reasonable agreement with the scaling relations, while the undetected objects, and especially the stacked non-detections, are significantly offset from the extrapolated relations. This could indicate either a strong break in the scaling relations at 3–5 lower sSFR than the star-forming sequence, or a large increase in the scatter of and .

4.2 Comparison to Cosmological Simulations

Recent large cosmological simulations have been able to produce populations of galaxies with properties broadly consistent with many observations (e.g., vogelsberger14; schaye15). Typical target benchmarks include matching the galaxy stellar mass function at and higher redshifts, the star formation history of the universe, galaxy colors, and other observables. The generally good agreement between the simulations and these basic observed properties then allows us to explore more detailed aspects of the simulations, which can help probe whether the sub-grid physical processes modeled by the simulations are accurately calibrated. Passive galaxies in particular are useful in this comparison, as they are the end products of powerful galactic feedback and quenching processes that encapsulate the most significant differences between various simulations. Here we compare the molecular gas properties of our observed passive galaxies to predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulation (schaye15).

Molecular gas masses for galaxies in the Ref-L100N1504 EAGLE run were calculated by lagos15; lagos16 using two different prescriptions for the molecular fraction that depend on, e.g., the metallicity, density, pressure, and SFR of each gas particle (gnedin11; krumholz13). In comparing the simulated galaxies’ SFRs and gas masses with observations at , lagos15 find good agreement for ‘normal’ star-forming galaxies. On the other hand, EAGLE underpredicts the number of gas-rich early-type galaxies compared to observations; early-type galaxies from the ATLAS survey are over-represented at the gas-rich end compared to the simulation, with 17% of observed galaxies lying 2 from the median SFR– relation in EAGLE. Although no morphological selection of early types was performed in the simulation, clouding the comparison, the difference could indicate that the star formation ‘law’ employed by EAGLE is too simplistic.

In Figure 4.2 we make a similar comparison to the simulations using our ALMA observations. This figure shows the SFR– star formation relation using a snapshot from EAGLE; we obtain a virtually identical result using a snapshot. For this comparison, we simply average the results from the two prescriptions used by lagos15 to calculate for the simulated galaxies; using one or the other alone also has no impact on our conclusions. The blue line and shaded region in Figure 4.2 show the median and 90% interval of SFR as a function of for all galaxies in EAGLE with , approximately the stellar mass cutoff of the LEGA-C survey.

Making a clean selection of simulated galaxies similar to our own selection criteria for ALMA observations is difficult. While we primarily selected objects from the LEGA-C survey to have sSFRs 3–10 below the star-forming sequence at , it is known that EAGLE does not accurately track the evolution of the star-forming sequence to high redshifts. In particular, the simulated galaxies underpredict the normalization of the star-forming sequence by 0.2–0.4 dex at (furlong15). This issue is not unique to EAGLE – the same is seen in other cosmological simulations (e.g., Illustris; sparre15). It is also worth noting that even the observations are not self-consistent on this point, with the integral of the cosmic SFR density over-predicting the observed stellar mass density by  dex (e.g., madau14, though see also driver18), while the simulations obviously must be self-consistent. Given the uncertainty in the normalizations of these relations in the literature, instead of selecting galaxies with the same range of sSFRs as the observations, we instead simply select the objects with sSFRs lower than 95% of the total population of simulated galaxies. Although the quantitative selection is different, this preserves the spirit of our selection, aiming for objects much more passive than the overall population. We also restrict the comparison to galaxies with , below which gas masses are unreliable due to the baryonic particle mass resolution of EAGLE. We preserve our other two main selection criteria, and SFR2.5 /yr, as in the observed sample, resulting in 91 simulated massive, passive galaxies. While these are absolute, rather than relative, selection cuts, they have little influence on our conclusions because, as we will show, massive and passive galaxies in EAGLE do not show a systematically different star formation relation compared to the rest of the simulated galaxies. The objects mimicking our selection are shown with small orange squares in Figure 4.2.

We find that our observed sample tends to lie near, but slightly offset from the bulk of simulated galaxies. The CO-detected sources, in particular, seem to have either too much molecular gas given their SFRs, or too little star formation given the amount of gas they contain. In other words, their depletion times are slightly longer than seen in the simulation. We draw the same conclusion from comparing only to those objects that mimic our selection criteria, which are not systematically offset from the rest of the galaxies in the simulation. The differences would be further exaggerated if the SFRs in the simulation or observations have been poorly estimated – either too low in EAGLE due to the underprediction of the star-forming sequence, or too high in our targets due to non-SF dust heating at 24 . If the SFRs of the observed and comparison simulated galaxies are generally correct, however, it appears that the differences are mostly in , and not in SFR. The simulated galaxies have SFRs comparable to the observed sample (partly by construction, due to our SFR selection threshold), but a median , with a dispersion of 0.2 dex and a tail that extends to lower molecular masses. Indeed, even the most gas-rich EAGLE galaxy matching our selection criteria has lower than any of our CO-detected objects.

Amusingly, unlike the observationally-derived scaling relations of Section 4.1, in this comparison it is the CO-undetected LEGA-C objects that are in better agreement with the scaling relation. While we have only upper limits on the depletion time of these objects, they would be consistent with the EAGLE galaxies even if they contained more than an order of magnitude less molecular gas than our present upper limits.

The better agreement for the non-detections appears to be due to the fact that EAGLE predicts a weak decline of towards high masses, in contrast to the increasing depletion times for massive galaxies predicted by the tacconi18 scaling relations. In EAGLE, the decline in is largely due to two linked effects. First, the probability distribution function of ISM gas densities is shifted to higher mean density in more massive galaxies. Second, EAGLE assumes a super-linear Schmidt-Kennicutt star formation relation. Combined with the higher gas densities in more massive galaxies, this leads to more star formation per unit molecular gas (and lower depletion times) at high masses. This effect would be muted if instead a linear relationship between SFR and were used in the simulations. EAGLE therefore does not and cannot reproduce the diversity of values that are unveiled by our observations of passive galaxies, suggesting that adopting a universal Schmidt-Kennicutt relation in the simulation is not sufficient to capture the complexity these observations reveal.

In Figure 4.2 (right), we examine the depletion times explicitly, comparing the predicted correlation between and from tacconi18 with that seen in the EAGLE simulation, where the navy shaded region corresponds to the prediction for galaxies a factor of 3 below (upper edge) and above (lower edge) the star-forming sequence at these redshifts. The two scaling relations predict opposite behavior for high-mass galaxies. The differences are especially stark for massive and passive galaxies (again shown with orange squares), which follow the same envelope as star-forming galaxies in EAGLE but are predicted to have depletion times longer by a factor of by the tacconi18 scaling relation.

Our observations, in summary, provide some support for both the observationally– and theoretically-derived scaling relations. While the gas-rich, CO-detected objects in our sample agree reasonably well with the extrapolations from tacconi18, our undetected objects have depletion times that are much more consistent with the EAGLE simulation at these redshifts (lagos15; lagos16). This highlights the importance of building a larger sample of passive galaxies with molecular mass estimates, which would indicate whether one or both of these scaling relations should be revised – how these scaling relations should be revised is unclear without additional information from a larger sample.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Left: The star formation relation from a snapshot of the EAGLE simulation, compared with our results. The individual objects and stacked non-detection from our observations are shown with red diamonds. The blue line and shaded region show the median and 90% range of all galaxies in EAGLE with . The small orange squares show galaxies in EAGLE that mimic our selection criteria; see text for details. The dotted grey line shows a constant depletion time  Gyr. Right: The gas depletion times as a function of stellar mass from the same EAGLE snapshot. The blue line and shaded region again shows the median and 90% range of depletion times for all EAGLE galaxies, with those that mimic our selection criteria as orange squares. The dashed navy line and region, on the other hand, show the predicted scaling relation from tacconi18 for galaxies on and a factor of 3 above (lower edge of region) and below (upper edge) the star-forming sequence at these redshifts. Note also that this Figure does not contradict our results in Figure 4.1, even though the upper limits from the CO non-detections overlap with the shaded region, because these galaxies have significantly lower sSFR than merely 3 below the star-forming sequence.

5 Implications for Galaxy Quenching

5.1 Active and Dynamical Quenching

In contrast to the nearly-linear relationship between and SFR observed for lower-mass galaxies from to , high-mass galaxies () tend to show a sub-linear relationship (e.g., whitaker14; schreiber15). This is generally interpreted as an effect of the ‘mass quenching’ of high-mass galaxies that results in a nearly-constant break in the stellar mass function at (peng10b). The various physical mechanisms that cause this mass quenching predict different observable properties of the ISM in passive galaxies, and can potentially be disentangled by observations of large samples of massive, quenched objects such as those presented here.

At , the local environment of a galaxy is correlated with whether or not it is quenched, with galaxies in dense regions much more likely to be quiescent than those in less dense environments (e.g., sobral11; darvish16). Numerous physical mechanisms can explain this environmental dependence, including ram pressure stripping, starvation, and tidal interactions. The effects of ‘maintenance-mode’ AGN feedback are also most prominent in denser environments (correa18). We found no difference in local overdensity between the CO-detected and undetected galaxies, although the large uncertainties in individual estimates of the local overdensities and the small sample size of our study makes it difficult to determine whether environmental effects are responsible for the variations in and we observe in our sample.

Dynamical processes can suppress the SFR in massive galaxies below the expected rates (e.g., martig09; genzel14b), even in the presence of large molecular masses, generally referred to as morphological or gravitational quenching. In this scenario the formation of a stellar bulge stabilizes an embedded gas disk against collapse, and the stabilization can proceed from the inner regions of a galaxy outwards as the stellar surface density increases over time. These dynamical processes can accommodate relatively high gas fractions, up to , but predict long depletion times reaching 10 Gyr or more; the molecular material is still present within the galaxies, but is prevented from forming stars.

Active, internal processes, on the other hand, such as feedback from AGN or star formation act to suppress star formation by heating and/or removing the molecular material (see fabian12 for a review). The energetics and driving mechanisms of these feedback processes are being actively debated, but generally predict low cold gas fractions and short depletion times.

Our observations reveal low gas fractions, , and short depletion times,  Gyr. The CO-undetected objects, in particular, show lower gas fractions and depletion times comparable to or shorter than expected for star-forming galaxies of similar mass, a conclusion that is robust even if the SFRs have been overestimated by . This demonstrates that intermediate-redshift passive galaxies are genuinely gas-poor, and the low SFRs are not the result of increasing . This suggests that dynamical processes are not acting to suppress the star formation in this sample. A similar conclusion was reached by sargent15 at slightly higher redshift, who failed to detect a single very massive galaxy at with , and concluded that gravitational stabilization of gas discs is not common among the most massive high-redshift quiescent galaxies.

Although there is some evidence that the morphological quenching scenario is consistent with the early-type population (martig13), we do not find evidence of it at higher redshift. There are several key differences between our sample and local massive early-types that are likely responsible for this. First, our (detected) sample galaxies contain 1–2 orders of magnitude more molecular gas than the population, and the upper limits on the undetected galaxies are still consistent with order of magnitude higher . While morphological quenching can accommodate gas fractions of , this may be difficult for galaxies to achieve in practice without resuming star formation. Second, our galaxies have stellar ages several Gyr younger than the early types (wu18; Chauke et al., submitted), and so are much closer to their epoch of quenching. martig09 show that morphological quenching can act when a galaxy has already assembled a stellar spheroid and nearly depleted its gas reservoir, but is slowly re-accreting gas. It is not clear that the same should hold if the galaxy never fully depletes its gas supply, but merely continues to form stars from its existing gas reservoir.

5.2 Connection to Local Quiescent Galaxies

Finally, we consider the possible connection between the passive galaxies we have observed at intermediate redshifts and their counterparts in the local universe. Given the molecular masses we have observed, it is immediately clear that if these objects are indeed progenitors of local massive quiescent galaxies, they must consume, expel, or heat 90–99% of the molecular material over the next 6 Gyr in order to match the range of seen in nearby ellipticals. This implies a much steeper drop in gas fraction compared to star-forming galaxies over the same redshift range, which decline by less than an order of magnitude over the same interval (e.g., tacconi18).

Part of the reason for the slower decline in for the star-forming population is that galaxies continue to accrete substantial amounts of gas over time from minor mergers and streams along the cosmic web (e.g., keres05). Because is observed to be much less than a Hubble time, galaxies’ gas reservoirs must be replenished if they are to continue forming stars at the observed high rates. The implied gas accretion rates are /yr at , declining by a factor of from to , and a further factor of to (scoville17). Galaxies on and near the star-forming sequence then grow largely in equilibrium with the supply of gas made available to them (e.g., dave11; lilly13; peng14).

Massive quiescent galaxies, in contrast, are generally consistent with passive evolution over many Gyr, with growth mostly influenced by gas-poor minor mergers but relatively little further star formation. Of course, the present-day quiescent population is comprised of galaxies that quenched star formation at a variety of redshifts, with a wide variety of evolutionary histories before and after the epoch of quenching. Here, we ask whether the galaxies we have observed at are consistent with the population of present-day massive quiescent galaxies, assuming no external gas accretion further influences their evolution after the epoch at which we observe them.

To address this question, we construct a closed-box toy model, in which a galaxy at with continues to deplete its reservoir of molecular gas in a similar fashion to the galaxies we have observed. Namely, we assume a range of initial gas fractions of at , which decline with time as low-level star formation continues. The evolving (declining) SFRs are determined by assuming a range of depletion times similar to our observations,  Gyr, that remain constant over time. We assume that a long-term average of 30% of the mass of stars formed is returned to the ISM through stellar mass loss (leitner11; this is also the same assumption made by scoville17 to infer high gas accretion rates for star-forming galaxies). Our conclusions are not sensitive to this assumption because the gas masses, and hence SFRs and additional stellar mass formed, are low. The result of this calculation is shown as the densely-dashed red line and region in Figure 5.2, and we compare to nearby quiescent galaxies with selected from the ATLAS and MASSIVE surveys (young11; davis16).

We find that a population of intermediate-redshift galaxies with properties similar to our observed sample can very naturally reproduce the range of seen in nearby passive galaxies, with no need for further gas accretion or dramatic increase in . Alternatively, if additional gas is accreted, it likely remains ionized or neutral, instead of molecular. The same conclusion does not hold for more ‘normal’ star-forming galaxies (blue and white points, and blue dashed line and region in Figure 5.2), which show gas fractions elevated well above what would be expected in the case of negligible gas accretion. If indeed galaxies similar to those we have observed are progenitors of the local quiescent population, the molecular gas fraction must fall more rapidly with redshift than it does for star-forming galaxies in order to reconcile the observations at and .

It is also worth noting that the low-level star formation in our toy model does not contradict conclusions about the ages of the stellar populations of quiescent galaxies at and lower redshift. wu18 compared the stellar age distributions for quiescent galaxies from LEGA-C at intermediate redshift with the ages of similarly-massive quiescent galaxies at from the Sloan Digital Sky Survey. They found a typical stellar age difference of only  Gyr, while the Universe aged  Gyr – in other words, the low-redshift quiescent galaxies are younger than would be expected in the case of pure passive evolution. wu18 surmised that either galaxies with younger stellar ages must join the quiescent population in the intervening time, or that low-level star formation in the existing quiescent galaxies must occur to lower the typical light-weighted stellar age. A true quantitative examination of this point is plagued by degeneracies, but at least qualitatively, the molecular gas fractions and depletion times observed in our sample should improve, rather than worsen, the differences in the stellar ages of quiescent objects at and 0.1.

Our finding that passive galaxies at and can be very naturally linked through internal processes continuing to deplete the gas reservoirs provides another independent line of evidence that preventing gas accretion and/or cooling to molecular form is key to maintaining low SFRs over many Gyr. Several proposed physical mechanisms can accomplish this, including preventing accretion of cold streams and/or satellite galaxies through the hot circumgalactic medium in massive halos (e.g., birnboim03; johansson09), and radio- or ‘maintenance-mode’ feedback from AGN to prevent cooling (e.g., croton06). Indeed, barisic17 concluded that radio-mode feedback is likely at work in the most massive galaxies in the LEGA-C sample, acting to maintain the observed low SFRs over several Gyr.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Evolution of with redshift for galaxies on and below the star-forming sequence. Our observations are shown with red diamonds, blue circles show galaxies near the star-forming sequence at (bauermeister13; tacconi13), white circles are star-forming galaxies with from xCOLDGASS (saintonge17), and small red symbols show massive () quiescent galaxies from the ATLAS and MASSIVE surveys (young11; davis16). For clarity of presentation, all detections and non-detections from these last two surveys are plotted at and , respectively; in reality all these galaxies are within 110 Mpc (). The densely-dashed red line and region sketch a plausible evolutionary path for the passive galaxies we observed from to , assuming initial starting and constant  Gyr, and no subsequent accretion of gas. The navy dashed line and region show the evolution for galaxies with and within a factor of 3 of the star-forming sequence from tacconi18, where the wide range of these parameters is intended to smooth over the trajectory of individual galaxies around the star-forming sequence. The continuity of this relation and the star-forming sequence requires substantial gas accretion over cosmic time. This Figure demonstrates that the galaxies we have observed at can naturally reproduce the observed gas fractions of the local massive, quiescent population simply by maintaining the same (relatively short) we measure, with no need for further gas accretion or dynamical stabilization.

6 Conclusions and Outlook

In this work we have presented ALMA observations of the CO(2–1) line in a sample of 8 massive galaxies at selected to lie significantly below the star-forming sequence of galaxies at this redshift, drawn from the LEGA-C spectroscopic survey. Our observations provide a first foray into understanding the molecular ISM properties of passive galaxies outside the local universe. With modest integration by ALMA, we significantly detect CO emission from four objects at , but do not detect the other four even after stacking, yielding upper limits on the molecular mass . At least three of the four detected galaxies show velocity gradients in the CO emission typical of rotation, and this rotation appears consistent with the stellar rotation seen in the LEGA-C spectra.

The main focus of this work has been to understand the molecular ISM properties of these passive galaxies in the context of star-forming and quiescent galaxies at . In particular, given the large investment in understanding how gas properties such as and vary with other galaxy properties for the star-forming population, we examined whether galaxies further below the star-forming sequence also follow these scaling relations. We find evidence that passive galaxies at intermediate redshifts have both lower gas fractions and shorter depletion times than expected based on extrapolations from observations of the star-forming population, particularly for the objects undetected in CO emission. This suggests either a break in the scaling relations at sSFRs below the star-forming sequence, or a large increase in the scatter of individual galaxy gas properties. On the other hand, a comparison to the EAGLE cosmological simulation indicates that the CO-undetected objects are in reasonable agreement with comparable simulated galaxies, while the CO-detected objects in our sample are somewhat more gas-rich than those seen in the simulation. In both cases, the differences between observations and the predicted scaling relations are likely linked to the physical processes responsible for galaxy quenching, which are not incorporated into the purely-empirical observational scaling relations and are uncertain in the simulation.

The generally short depletion times we observe are more consistent with active feedback processes preventing star formation, rather than dynamical or gravitational stabilization against collapse of the existing gas reservoirs. We also show that a population of galaxies like those we have observed can very naturally reproduce the observed gas fractions of early-type galaxies while maintaining the depletion times we observe, with no need for additional gas accretion or increased depletion times predicted by dynamical stabilization. We see no evidence of environmental effects causing the variations in and we observe, though we stress that the small sample size of our study makes this inconclusive.

Developing a fuller understanding of the range of gas masses, fractions and depletion times in high-redshift quiescent galaxies will require observations of a significantly larger sample of objects than we have conducted here. This will allow us to determine whether passive galaxies truly show a break towards low compared to expectations or whether this is merely an artifact of limited sample size. While the dynamical information provided by CO spectroscopic observations is valuable, it is also clear that a larger sample of passive galaxies will require alternate methods of determining gas masses, such as observing long-wavelength dust continuum emission. Detecting gas fractions at the level of the stacked non-detections in CO with ALMA, for example, would require over an order of magnitude more observing time than we have invested here; reaching the same level in 870  dust continuum would require about the same amount of on-source integration as our present observations. Clearly the detection of molecular gas in passive galaxies at high-redshift remains a significant challenge, even with the sensitivity of ALMA. Nevertheless, these observations present a new avenue to probe the physics of galaxy quenching and the emergence of the dichotomy between star-forming and quiescent galaxies.


JS thanks the McDonald Observatory at the University of Texas at Austin for support through a Harlan J. Smith Fellowship. CL is funded by a Discovery Early Career Researcher Award (DE150100618) of the Australian Research Council. CCW acknowledges support from the National Science Foundation Astronomy and Astrophysics Fellowship grant AST-1701546. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.00790.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 194-A.2005 (The LEGA-C Public Spectroscopic Survey).

Facility: ALMA

\bibliography@latex

legac_alma.bbl

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
187321
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description