- Split View
-
Views
-
Cite
Cite
Thiago S. Gonçalves, Antara Basu-Zych, Roderik A. Overzier, Laura Pérez, D. Christopher Martin, Molecular gas properties of UV-bright star-forming galaxies at low redshift, Monthly Notices of the Royal Astronomical Society, Volume 442, Issue 2, 1 August 2014, Pages 1429–1439, https://doi.org/10.1093/mnras/stu852
- Share Icon Share
Abstract
Lyman break analogues (LBAs) are a population of star-forming galaxies at low redshift (z ∼ 0.2) selected in the ultraviolet (UV). These objects present higher star formation rates and lower dust extinction than other galaxies with similar masses and luminosities in the local universe. In this work, we present results from a survey with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) array to detect CO(1–0) emission in LBAs, in order to analyse the properties of the molecular gas in these galaxies. Our results show that LBAs follow the same Schmidt–Kennicutt law as local galaxies. On the other hand, they have higher gas fractions (up to 66 per cent) and faster gas depletion time-scales (below 1 Gyr). These characteristics render these objects more akin to high-redshift star-forming galaxies. We conclude that LBAs are a great nearby laboratory for studying the cold interstellar medium in low-metallicity, UV-bright compact star-forming galaxies.
1 INTRODUCTION
More recently, Bigiel et al. (2008) analysed data from nearby SFGs to resolve the S-K relation down to scales of 750 pc. These authors found that the relation in the nearby universe can be divided into two: one for atomic gas alone, which saturates at approximately 10 M⊙ pc−2; and a linear relation (N = 1) for molecular gas, mostly above that saturation level. The linearity has an interesting consequence: since the S-K law relates star formation rates (SFR) and the gas reservoir, an N = 1 index means a constant gas depletion time-scale for molecular gas in observed SFGs, |$t_{\rm depletion}=\Sigma _{\text{H}_2}/\Sigma _{\rm SFR}$|. This value was measured to be approximately 2 × 109 Gyr, with 1σ scatter of around 0.24 dex, or 75 per cent (Bigiel et al. 2008).
In an attempt to offer a theoretical scenario for the formation of stars in giant molecular clouds, Krumholz, McKee & Tumlinson (2009) proposed a simple model listing three main factors that determined the gas depletion time-scales: gas fractions, internal star formation feedback and turbulence. At first the fraction of gas available for star formation is dependent on self-shielding from the interstellar radiation field – it is this factor that sets the threshold of conversion from atomic to molecular hydrogen. Once molecular clouds are formed, feedback determines their properties, since the internal pressure is higher at this stage than the average interstellar medium (ISM) pressure. Ongoing star formation in these objects is regulated by turbulence to a universal rate of 1 per cent of the mass per free-fall time. These can be combined into an analytic formalism capable of predicting the linear correlation between gas surface densities and constant depletion time-scales.
Understanding the detailed process through which molecular gas translates into stars becomes even more complicated for extreme objects, as has been shown for ULIRGs (Lbol > 1012 L⊙; Sanders & Mirabel 1996) both at low and high redshifts. The observed gas density in these galaxies is well above the limit of ∼100 M⊙ pc−2 mentioned in the work of Krumholz et al. (2009), which means the assumption of the properties in giant molecular clouds being determined by internal dynamics is no longer valid, and the end result is likely a distinct scenario for conversion of gasbrk into stars.
In an effort to understand the conversion of gas into stars in these environments, Bouché et al. (2007) and Tacconi et al. (2008) have measured the surface density of molecular gas in submillimetre-selected high-redshift ULIRGs, and have determined that the exponent in the S-K law seems to be higher for these objects, i.e. the star formation surface density is higher than expected for a given value of gas surface density. The same is true for intermediate-redshift ULIRGs: Combes et al. (2011) have measured CO luminosities of 30 galaxies between 0.2 < z < 0.6 and have found star formation efficiencies three times higher than those found in the local universe.
Nevertheless, these studies are still biased towards more massive galaxies, simply because smaller objects are faint and difficult to observe. These massive objects, in turn, tend to be dynamically colder, with more significant rotational support and a more metal-rich ISM than less massive ones (Maiolino et al. 2008; Gonçalves et al. 2010). Moreover, they probably do not represent the ‘typical’ SFG at high-z, since the mass function of galaxies is steep at these epochs (Reddy & Steidel 2009).
Until recently, only two Lyman break galaxies (LBGs) had been observed, both of which are strongly lensed and therefore highly magnified (Baker et al. 2004; Coppin et al. 2007; Riechers et al. 2010). Later, Magdis et al. (2012a) have detected one unlensed LBG; this object, however, is very massive, with log M*/M⊙ ∼ 1011. In order to understand how the star formation process operates on all scales, on the other hand, it is important to study SFGs at all masses. We conclude that our understanding of the conversion of gas into stars in such conditions could benefit greatly from the detailed study of the cold molecular gas in low-redshift galaxies with properties similar to those commonly found around z ∼ 2, allowing us to target multiple objects with less extreme CO luminosities which are likely more representative of the typical SFGs in these epochs.
Lyman break analogues (LBAs) offer an elegant solution to this problem. These objects, present at z ∼ 0.2, show a number of properties that render them more akin to high-redshift UV-selected SFGs than normal spirals in the local universe (see Section 2.1). Therefore, we are able to study the cold ISM in these galaxies in order to infer how the conversion of the molecular gas reservoir into stars happens under such conditions.
In this work, we present results from a molecular gas survey of five LBAs, observed with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) interferometer. We target the lowest CO transition, CO(1–0), to most accurately trace the cold molecular gas in these galaxies. This paper is divided as follows: in Section 2, we present details on sample selection, observations and data reduction; in Section 3, we present the results from our survey, comparing properties such as gas masses and gas fractions with those found at high redshift; and in Section 4 we discuss our results, with caveats (including the implications of different αCO factors at low metallicities) and future perspectives on our work. We summarize our results in Section 5.
2 OBSERVATIONS AND DATA REDUCTION
2.1 The sample of LBAs
Heckman et al. (2005) and Hoopes et al. (2007) have defined a sample of ultraviolet-luminous galaxies (UVLGs) as those objects with extreme far-ultraviolet (FUV) luminosities, i.e. LFUV > 2 × 1010 L⊙ at redshifts z ∼ 0.2 from a combination of the SDSS (Abazajian et al. 2009) and GALEX catalogues (Martin et al. 2005). In order to exclude large spirals and select those galaxies with similar properties to high-redshift SFGs, these authors have further refined their criteria to include only galaxies with high values of surface brightness, IFUV > 109 L⊙ kpc−2.
Subsequent work has shown that these supercompact UVLGs do indeed resemble high-redshift UV-bright SFGs, especially the well-studied sample of LBGs (see Shapley 2011, and references therein). Hoopes et al. (2007) demonstrated that supercompact UVLGs shared many properties with LBGs: generally speaking, they are more luminous in UV bands, more metal poor, present higher SFR and less dust extinction than typical galaxies with similar masses at low redshift, being instead more similar to their counterparts at z ∼ 2. This justifies our denomination of the sample as LBAs.
Overzier et al. (2010) showed that the morphologies of LBAs are also distinct. By using quantitative measurements of their morphological structure (such as Gini and M20; Lotz et al. 2006), the authors have shown that these galaxies are irregular, and smaller than galaxies with similar masses at low redshift (with typical half-light radii R50 ∼ 1 kpc). Simulated observations of the sample at z ≳ 2 showed that LBAs are very similar to LBGs, with many of the asymmetries and giant star-forming clumps regularly seen at high redshift (e.g. Elmegreen et al. 2009).
Similarly, IFU studies of the ionized gas in LBAs have shown similarities in the kinematic structures between both samples – with ratios between rotational velocities and dispersion in the gas of vcirc/σ ≲ 1, much less than typical values above 10 in local spirals – although there are indications that observational limitations at high redshift, such as the loss in spatial resolution, makes objects appear more symmetric than in reality (Gonçalves et al. 2010).
More recently, a study on the infrared properties of LBAs has shown that, while these objects have bolometric luminosities similar to local luminous infrared galaxies (LIRGs), with typical values above LIR > 1011 L⊙, their dust obscuration, as measured by the infrared excess LIR/LFUV, are at least an order of magnitude weaker (Overzier et al. 2011). Again, these values are very similar to what is found in LBGs.
Finally, recent work by Basu-Zych et al. (2013b) using the Chandra telescope has shown that the same analogy applies to the X-ray properties of LBAs. These galaxies show elevated X-ray luminosities, with LX ∼ 4 times the typical value of galaxies with similar SFR in the local universe. The authors argue that this is a result of the low metallicity of the objects driving the formation of luminous high-mass X-ray binaries. Once more, this is in good agreement with findings for high-redshift SFGs (e.g. Basu-Zych et al. 2013a).
2.2 CARMA sample selection
Our sample was drawn from the same parent sample as that first presented in Overzier et al. (2009). These objects have been extensively studied with multiple instruments and in a wide variety of wavelengths (e.g. Basu-Zych et al. 2007; Gonçalves et al. 2010; Overzier et al. 2010; Heckman et al. 2011; Alexandroff et al. 2012), and as such will provide useful data in future work comparing other physical properties, such as gas kinematics, with the molecular gas observations presented here.
In this work we have assumed αCO = 4.6, to ensure consistency with mass values determined at low and high redshift. This is also consistent with most high-redshift studies of ‘main-sequence’ galaxies (Daddi et al. 2010a; Tacconi et al. 2010; Magdis et al. 2012a). We emphasize that, given their SFR and stellar masses, galaxies in our sample would be classified as starbursts, with typical specific star formation rates, sSFR ≡ SFR/M* ∼1 Gyr−1, well above other ‘main-sequence’ objects at z ∼ 0.2 (Elbaz et al. 2011). Nevertheless, they are remarkably distinct from typical dusty starbursts at the same epoch, instead following the same relation between infrared excess and UV slope as other low-redshift SFGs and high-redshift LBGs (Meurer, Heckman & Calzetti 1999; Overzier et al. 2011). Therefore, we consider our galaxies to be more akin to early ‘main-sequence’ objects, justifying our choice of αCO. We discuss the implications and possible caveats associated with this selection in Section 4.3.
Following equation (4), we have inferred gas masses for all galaxies in our sample from observed SFR and optical half-light radii, and converted those values into expected CO fluxes. Table 1 summarizes the physical properties of our observed sample, and Fig. 1 shows Hubble Space Telescope (HST) images of all six objects. The stellar masses were taken from the SDSS/DR7 MPA-JHU value-added catalogue.1 These masses were calculated by fitting a large grid of spectral synthesis models from Bruzual & Charlot (2003) to the SDSS u′, g′, r′, i′, z′ photometry, using a Chabrier (2003) initial mass function. Gas-phase metallicities were calculated using SDSS spectra from the same catalogue, following the methodology in Tremonti et al. (2004). Half-light radii are derived from HST imaging in optical bands (Overzier et al. 2010). SFR are measured from combined Hα and MIPS 24 μm data; they typically present an uncertainty up to 0.3 dex (Overzier et al. 2009). FIR luminosities are calculated by fitting the model library of Siebenmorgen & Krügel (2007) to Spitzer IRAC MIPS and IRS data points (Overzier et al. 2011). Finally, dynamical time-scales are calculated by taking the ratio between half-light radii and measured circular velocities (Gonçalves et al. 2010), in order to ensure consistency with similar studies at high redshift (Daddi et al. 2010a).
ID . | z . | log LFUV . | log IFUV . | log LFIR . | r1/2 . | SFR . | log M* . | 12 + log(O/H) . | τdyn . |
---|---|---|---|---|---|---|---|---|---|
. | . | (L⊙) . | (L⊙ kpc−2) . | (L⊙) . | (kpc) . | (M⊙ yr−1) . | (M⊙) . | . | (Myr) . |
001054 | 0.243 | 10.42 | 9.01 | ... | 4.24 | 26.9 | – | 8.7 | 27 |
015028 | 0.147 | 10.62 | 9.38 | 11.4 | 1.83 | 50.7 | 10.3 | 8.5 | 17 |
080844 | 0.091 | 10.45 | 10.25 | 10.5 | 0.88 | 16.1 | 9.8 | 8.8 | 7 |
092159 | 0.235 | 10.82 | 9.65 | 11.6 | 1.80 | 55.1 | 10.8 | 8.8 | 13 |
210358 | 0.137 | 10.52 | 9.58 | 11.4 | 2.70 | 108.3 | 10.9 | 8.8 | 15 |
231812 | 0.252 | 10.85 | 9.33 | 11.2 | 2.54 | 63.1 | 10.0 | 8.4 | 30 |
ID . | z . | log LFUV . | log IFUV . | log LFIR . | r1/2 . | SFR . | log M* . | 12 + log(O/H) . | τdyn . |
---|---|---|---|---|---|---|---|---|---|
. | . | (L⊙) . | (L⊙ kpc−2) . | (L⊙) . | (kpc) . | (M⊙ yr−1) . | (M⊙) . | . | (Myr) . |
001054 | 0.243 | 10.42 | 9.01 | ... | 4.24 | 26.9 | – | 8.7 | 27 |
015028 | 0.147 | 10.62 | 9.38 | 11.4 | 1.83 | 50.7 | 10.3 | 8.5 | 17 |
080844 | 0.091 | 10.45 | 10.25 | 10.5 | 0.88 | 16.1 | 9.8 | 8.8 | 7 |
092159 | 0.235 | 10.82 | 9.65 | 11.6 | 1.80 | 55.1 | 10.8 | 8.8 | 13 |
210358 | 0.137 | 10.52 | 9.58 | 11.4 | 2.70 | 108.3 | 10.9 | 8.8 | 15 |
231812 | 0.252 | 10.85 | 9.33 | 11.2 | 2.54 | 63.1 | 10.0 | 8.4 | 30 |
ID . | z . | log LFUV . | log IFUV . | log LFIR . | r1/2 . | SFR . | log M* . | 12 + log(O/H) . | τdyn . |
---|---|---|---|---|---|---|---|---|---|
. | . | (L⊙) . | (L⊙ kpc−2) . | (L⊙) . | (kpc) . | (M⊙ yr−1) . | (M⊙) . | . | (Myr) . |
001054 | 0.243 | 10.42 | 9.01 | ... | 4.24 | 26.9 | – | 8.7 | 27 |
015028 | 0.147 | 10.62 | 9.38 | 11.4 | 1.83 | 50.7 | 10.3 | 8.5 | 17 |
080844 | 0.091 | 10.45 | 10.25 | 10.5 | 0.88 | 16.1 | 9.8 | 8.8 | 7 |
092159 | 0.235 | 10.82 | 9.65 | 11.6 | 1.80 | 55.1 | 10.8 | 8.8 | 13 |
210358 | 0.137 | 10.52 | 9.58 | 11.4 | 2.70 | 108.3 | 10.9 | 8.8 | 15 |
231812 | 0.252 | 10.85 | 9.33 | 11.2 | 2.54 | 63.1 | 10.0 | 8.4 | 30 |
ID . | z . | log LFUV . | log IFUV . | log LFIR . | r1/2 . | SFR . | log M* . | 12 + log(O/H) . | τdyn . |
---|---|---|---|---|---|---|---|---|---|
. | . | (L⊙) . | (L⊙ kpc−2) . | (L⊙) . | (kpc) . | (M⊙ yr−1) . | (M⊙) . | . | (Myr) . |
001054 | 0.243 | 10.42 | 9.01 | ... | 4.24 | 26.9 | – | 8.7 | 27 |
015028 | 0.147 | 10.62 | 9.38 | 11.4 | 1.83 | 50.7 | 10.3 | 8.5 | 17 |
080844 | 0.091 | 10.45 | 10.25 | 10.5 | 0.88 | 16.1 | 9.8 | 8.8 | 7 |
092159 | 0.235 | 10.82 | 9.65 | 11.6 | 1.80 | 55.1 | 10.8 | 8.8 | 13 |
210358 | 0.137 | 10.52 | 9.58 | 11.4 | 2.70 | 108.3 | 10.9 | 8.8 | 15 |
231812 | 0.252 | 10.85 | 9.33 | 11.2 | 2.54 | 63.1 | 10.0 | 8.4 | 30 |
2.3 CO observations and data analysis
We conducted our observations between 2010 April and 2011 June. The summary of observations, including observing date and total integration time, can be found in Table 2. We have selected the CO(1–0) transition (νrest = 115.271 GHz) as the line of choice for these observations. Although higher transitions would provide us with higher signal-to-noise ratios (|$\text{S/N}$|) for the same integration times, targeting the lowest CO transition allows us to better trace the cold molecular mass, without the need to rely on any assumptions about gas excitations levels. To complicate matters, excitation measurements for high-redshift colour-selected galaxies are scarce (e.g. Dannerbauer et al. 2009; Aravena et al. 2010; Riechers et al. 2010). At zeroth order, one could infer that CO luminosities are constant throughout all transitions, assuming the gas is thermalized. However, higher transitions require higher gas densities and may be tracing distinct regions of the ISM (Harris et al. 2010); furthermore, the gas may be subthermally excited (e.g. Weiss et al. 2007; Riechers et al. 2009). In addition, observations in the 3 mm band are simpler, since weather requirements are not as strict as observations in shorter wavelengths.
ID . | RA . | Dec. . | νobs . | Integration . |
---|---|---|---|---|
. | . | . | (GHz) . | time (h) . |
001054 | 00:10:54.85 | 00:14:51.35 | 92.7 | 7.2 |
015028 | 01:50:28.41 | 13:08:58.40 | 100.5 | 20.7 |
080844 | 08:08:44.27 | 39:48:52.36 | 105.6 | 3.4 |
092159 | 09:21:59.39 | 45:09:12.38 | 93.3 | 5.4 |
210358 | 21:03:58.75 | −07:28:02.45 | 101.4 | 15.3 |
231812 | 23:18:13.00 | −00:41:26.10 | 92.1 | 16.1 |
ID . | RA . | Dec. . | νobs . | Integration . |
---|---|---|---|---|
. | . | . | (GHz) . | time (h) . |
001054 | 00:10:54.85 | 00:14:51.35 | 92.7 | 7.2 |
015028 | 01:50:28.41 | 13:08:58.40 | 100.5 | 20.7 |
080844 | 08:08:44.27 | 39:48:52.36 | 105.6 | 3.4 |
092159 | 09:21:59.39 | 45:09:12.38 | 93.3 | 5.4 |
210358 | 21:03:58.75 | −07:28:02.45 | 101.4 | 15.3 |
231812 | 23:18:13.00 | −00:41:26.10 | 92.1 | 16.1 |
ID . | RA . | Dec. . | νobs . | Integration . |
---|---|---|---|---|
. | . | . | (GHz) . | time (h) . |
001054 | 00:10:54.85 | 00:14:51.35 | 92.7 | 7.2 |
015028 | 01:50:28.41 | 13:08:58.40 | 100.5 | 20.7 |
080844 | 08:08:44.27 | 39:48:52.36 | 105.6 | 3.4 |
092159 | 09:21:59.39 | 45:09:12.38 | 93.3 | 5.4 |
210358 | 21:03:58.75 | −07:28:02.45 | 101.4 | 15.3 |
231812 | 23:18:13.00 | −00:41:26.10 | 92.1 | 16.1 |
ID . | RA . | Dec. . | νobs . | Integration . |
---|---|---|---|---|
. | . | . | (GHz) . | time (h) . |
001054 | 00:10:54.85 | 00:14:51.35 | 92.7 | 7.2 |
015028 | 01:50:28.41 | 13:08:58.40 | 100.5 | 20.7 |
080844 | 08:08:44.27 | 39:48:52.36 | 105.6 | 3.4 |
092159 | 09:21:59.39 | 45:09:12.38 | 93.3 | 5.4 |
210358 | 21:03:58.75 | −07:28:02.45 | 101.4 | 15.3 |
231812 | 23:18:13.00 | −00:41:26.10 | 92.1 | 16.1 |
The array was used at D-configuration, providing a synthesized beam size of approximately 5 arcsec at the observed frequencies. This is slightly larger than the size of the observed galaxies, and we do not expect to resolve any of them at this stage of the survey. Proposed integration times were calculated aiming for a 5σ detection at the emission peak assuming the aforementioned expected CO fluxes and velocity dispersion as found in Gonçalves et al. (2010), smoothed over 10 km s−1 channels. Actual times, however, varied according to instrument availability and weather.
We used the miriad package2 to perform the CO(1–0) data analysis, choosing baseline solutions for each individual observed track. Data were mainly flagged according to variations in system temperature and resulting flux for a given gain calibrator, although specifics vary greatly on an individual basis. When data were flagged only for a fraction of a track, care was taken always to include a gain calibrator observation up to 10 min before or after each astronomical data point. In the case of multiple observed tracks for a single source, those were later combined through use of the uvcat task in miriad.
Visibility files were later inverted from UV plane into real spatial/velocity flux data cubes. Velocity resolution varies according to S/N in the observations, but in most cases spectra shown here are defined in steps of 10 km s−1, with 30 km s−1 smoothing. ‘Dirty’ maps were cleaned through use of the ‘clean’ procedure, and care was taken never to use a cutoff below the observed peak flux in Jy/beam of the galaxy.
We used the casa package to determine the total flux in each velocity channel. The integrating area was equivalent to an ellipse with same size as the beam. In Fig. 2 we show the CO(1–0) line profile of galaxies in our sample as a function of velocity with respect to the measured systemic redshift in each case.
3 RESULTS
Fig. 2 shows the spectra of all observed galaxies, in mJy. Spectral channels are binned at 10 km s−1, and smoothed over 30 km s−1 (three channels) each. To determine the total CO flux in the galaxy, we fit a Gaussian to the spectrum of each object. Flux is then given simply by the area below the curve. Total flux for galaxies with multiple components (such as 015028 and 231812) is given by the sum of the integrated areas of each Gaussian. All galaxies that show multiple components show clear signs of merger activity, as can be seen in Fig. 1. High-resolution interferometry should be able to decompose different spectral features into distinct spatial components. In the case of 210358, we extrapolate this integral to frequencies not covered by the observations, since bandwidth was not enough to cover the whole curve – this might explain why the measured CO luminosity is higher than estimated (Fig. 3). We also show the spectrum for 001054, but recent work (Overzier et al., in preparation) has shown that this object is instead an AGN. We show the CO(1–0) results for this galaxy in Table 3 for completeness, but exclude it from subsequent analysis.
ID . | SCOΔν . | |$L^{\prime }_{\rm CO}$| . | M(H2) . | μgas . |
---|---|---|---|---|
. | (Jy km s−1) . | (109 K km s−1 pc2) . | (109 M⊙) . | . |
001054 | 2.19 ± 0.41 | 6.5 | 29.7 | 0.27 |
015028 | 2.37 ± 0.28 | 2.5 | 11.4 | 0.36 |
080844 | 1.74 ± 0.56 | 0.7 | 3.1 | 0.33 |
092159 | 1.04 ± 0.31 | 2.8 | 13.1 | 0.17 |
210358 | 11.5 ± 1.5 | 10.3 | 47.7 | 0.37 |
231812 | 1.34 ± 0.23 | 4.2 | 1.95 | 0.66 |
ID . | SCOΔν . | |$L^{\prime }_{\rm CO}$| . | M(H2) . | μgas . |
---|---|---|---|---|
. | (Jy km s−1) . | (109 K km s−1 pc2) . | (109 M⊙) . | . |
001054 | 2.19 ± 0.41 | 6.5 | 29.7 | 0.27 |
015028 | 2.37 ± 0.28 | 2.5 | 11.4 | 0.36 |
080844 | 1.74 ± 0.56 | 0.7 | 3.1 | 0.33 |
092159 | 1.04 ± 0.31 | 2.8 | 13.1 | 0.17 |
210358 | 11.5 ± 1.5 | 10.3 | 47.7 | 0.37 |
231812 | 1.34 ± 0.23 | 4.2 | 1.95 | 0.66 |
ID . | SCOΔν . | |$L^{\prime }_{\rm CO}$| . | M(H2) . | μgas . |
---|---|---|---|---|
. | (Jy km s−1) . | (109 K km s−1 pc2) . | (109 M⊙) . | . |
001054 | 2.19 ± 0.41 | 6.5 | 29.7 | 0.27 |
015028 | 2.37 ± 0.28 | 2.5 | 11.4 | 0.36 |
080844 | 1.74 ± 0.56 | 0.7 | 3.1 | 0.33 |
092159 | 1.04 ± 0.31 | 2.8 | 13.1 | 0.17 |
210358 | 11.5 ± 1.5 | 10.3 | 47.7 | 0.37 |
231812 | 1.34 ± 0.23 | 4.2 | 1.95 | 0.66 |
ID . | SCOΔν . | |$L^{\prime }_{\rm CO}$| . | M(H2) . | μgas . |
---|---|---|---|---|
. | (Jy km s−1) . | (109 K km s−1 pc2) . | (109 M⊙) . | . |
001054 | 2.19 ± 0.41 | 6.5 | 29.7 | 0.27 |
015028 | 2.37 ± 0.28 | 2.5 | 11.4 | 0.36 |
080844 | 1.74 ± 0.56 | 0.7 | 3.1 | 0.33 |
092159 | 1.04 ± 0.31 | 2.8 | 13.1 | 0.17 |
210358 | 11.5 ± 1.5 | 10.3 | 47.7 | 0.37 |
231812 | 1.34 ± 0.23 | 4.2 | 1.95 | 0.66 |
Three of the five remaining objects are classified as Dominant Central Objects (DCOs) by Overzier et al. (2009), which means the presence of a very compact source in their centre: 080844, 092159, and 210358. It has been speculated that these galaxies could host AGN, since they are classified as composite starburst-AGN in typical emission line ratios used as diagnostics for nuclear activity. Jia et al. (2011) have studied the X-ray properties of DCOs and concluded that, while they do show signs of obscured AGN activity, the black hole masses estimated are low (between 105 and 106 M⊙, compared to total masses up to 109 M⊙ of individual DCOs) and the bolometric luminosities are dominated by the star-forming activity. Likewise, Alexandroff et al. (2012) have observed four DCOs with the European very long baseline interferometry Network – including all three objects in this work – in search of point sources in the galaxy centre. Only one of them (092159) was detected, at 1.7 GHz. The black hole mass is consistent with estimates by Jia et al. (2011), with an upper limit of 107 M⊙. This leads us to believe that DCOs do host the seeds of black holes commonly found in galactic bulges, but the UV light is still predominantly due to star-forming activity, generating strong superwinds with velocities up to 1500 km s−1 (Heckman et al. 2011). In this work, we treat all five objects indistinctly, but high-resolution interferometry should clarify whether the gas distribution is also concentrated in the centre of these three objects.
Fluxes are converted to CO luminosities according to the relation given by equation (5). Errors in flux measurements are the sum of the standard deviation on each velocity channel over 2σ (where σ is the velocity dispersion measured from the Gaussian fitting). In Table 3 we present the total CO(1–0) flux measured for each galaxy, along with inferred CO luminosities, molecular gas mass and gas fractions.
3.1 Gas masses and fractions
In Fig. 4 we show the LBA gas fractions as a function of stellar mass, including also additional LBAs that have not yet been observed with gas fraction values inferred from inverting the S-K relation [again following equation (4), shown as red circles]. The first thing to notice is that measured values are in agreement with expectations, given the uncertainties. For comparison, we also include gas fractions estimated in the same manner for LBGs from Erb et al. (2006, black dots), as well as estimates at redshifts between z = 1 and 2.5 from the observed evolution in gas depletion times and sSFR of SFGs (Tacconi et al. 2013, and references therein). There is a clear trend in all cases of more massive galaxies showing smaller gas fractions (as found in Tacconi et al. 2013), possibly indicating these galaxies have already consumed a larger portion of their gas reservoir and converted it into stellar mass. Gas fractions are also remarkably higher than typically found in the local universe, where most SFGs have less than 10 per cent of their baryonic mass in the form of molecular hydrogen (Saintonge et al. 2011, indicated by the yellow shaded area in Fig. 4). Another possibility is that star formation feedback plays an important role in regulating the infall of gas from the intergalactic medium (Tacconi et al. 2013). This trend – especially given the high gas fractions found – is well reproduced by our LBA sample, again reinforcing the validity of the analogy between both galaxy populations. We estimate lower gas fractions for less massive galaxies than Erb et al. (2006); this may be due to a lower limit in detectable gas fractions in their sample.
3.2 The Schmidt–Kennicutt relation at low and high redshift
Whether the S-K relation holds true for high-redshift SFGs – where higher SFR densities have been shown to exist (Shapley 2011) – is still subject to some debate, as we have noted in Section 1. With reliable measurement of CO(1–0) luminosities for some of the densest ‘normal’ galaxies (i.e. not ULIRGs) ever observed, we expect that the LBA sample can shed some light on the issue.
In Fig. 5, we show a compilation of the S-K relation for SFGs. The low-density end of the plot shows low-redshift spiral galaxies from the sample of Kennicutt (1998). The remainder of the points show low- and high-redshift infrared bright samples (Kennicutt 1998; Bouché et al. 2007) and ‘main-sequence’ SFGs at z > 1 from Tacconi et al. (2010), Daddi et al. (2010a) and Tacconi et al. (2013). To ensure consistency, we have estimated gas surface densities for all galaxies on the lower relation [excluding submillimeter galaxies (SMGs) and ULIRGs] with a uniform αCO = 4.6 value, using the total gas and size measurements mentioned in the original studies.
As we can see, the relation is apparently bimodal, as indicated by both the solid and the dashed lines. The solid line is a fit to the ‘main-sequence’ galaxies at all redshifts, while the dashed line has the same slope, only displaced by 0.9 dex to fit through the ULIRGs and SMGs. Physically speaking, that means the latter are almost 10 times as efficient in converting cold gas into stars, although the clear separation between the relations could be a consequence of using two distinct αCO values (Carilli & Walter 2013, and references therein).
Our sample, however, presents short dynamical time-scales (of order a few tens of Myr), and appear to lie below the relation inferred by Daddi et al. (2010a). This might be a result of how dynamical time-scales are defined; we discuss in Section 4.2 whether the choice of the rotational periods of galaxies as dynamical time-scales is valid for all samples.
4 DISCUSSION
One must be careful to take into account possible caveats in order not to overinterpret the data. In this section, we list possible caveats that might affect how we view our results: the distinction between atomic and molecular gas, the adopted definition for dynamical time-scale and the αCO conversion factor.
4.1 Atomic and molecular hydrogen
First of all, one needs to be careful regarding total gas versus the molecular gas alone. At high redshifts, in particular, where the direct detection of atomic hydrogen is nearly impossible, it seems that total and molecular gas are used indistinguishably.
In fact, at the observed scales, it is true we expect the molecular gas to dominate the baryonic mass in the galaxy. As Bigiel et al. (2008) pointed out, the surface density of H i saturates at about 10 M⊙ pc−2, at least in the local universe. In our LBA sample the surface density of (molecular) gas is about a hundred times higher; therefore we would expect the atomic hydrogen to dominate the bulk of the mass only at distances greater than 10 times its optical radii, assuming constant surface density at this saturation level over the whole area.
4.2 The dynamical time-scale
Likewise, we need to ensure that we understand the meaning of dynamical time in these galaxies. In principle, the dynamical time represents the amount of time a single star takes to orbit around the galaxy. This time-scale can be linked to the duration of processes that can spread throughout a given radius; for instance, a simple model for star formation in a galaxy through the collapse of cold molecular gas requires the cooling time to be shorter than the dynamical time-scale within a given radius.
In the case of extreme SFGs, however, the definition of the dynamical time-scale is not so clear. In Fig. 6, for example, all ‘main-sequence’ galaxies at high redshift (and our own LBA sample, for consistency) have dynamical time-scales defined as the rotation time-scale at the half-light radius, while for the local sample of spirals and ULIRGs (Kennicutt 1998) this quantity is defined at the outer radii of galaxies.
More importantly, though, we question the validity of characterizing dynamical time-scales in all SFGs by use of circular velocities. Velocity dispersions in SMGs, in particular, are remarkably high, so perhaps using tdyn ∼ r/σ might be a better description (for a detailed discussion of the resolved gas kinematics in SMGs, see Menéndez-Delmestre et al. 2013). In the case of LBAs, Gonçalves et al. (2010) have shown that not all galaxies can be accurately described as rotating discs. Likewise, Förster Schreiber et al. (2009) have shown that about 1/3 of their SINS sample are dispersion dominated. Since most of our galaxies present v/σ values smaller than 1, a different definition of tdyn would place our data on the same relation as Daddi et al. (2010a). Furthermore, it has been found that there is a correlation between stellar masses and the amount of rotational dynamical support (Gonçalves et al. 2010; Newman et al. 2013); this result, combined with the steepness of the stellar mass function at high redshift (Reddy & Steidel 2009), would indicate that the percentage of dispersion-dominated galaxies in the high-redshift population could be even higher.
According to Krumholz, Dekel & McKee (2012), the time-scale for high-redshift galaxies is related to the free-fall time for gas, which is dominated not by internal dynamics of giant molecular clouds, but instead by the bulk motion of the ISM. While the formation of instabilities in a rotating disc can describe well the advent of clumps in massive galaxies, we argue that turbulent motion makes for a more complex scenario in smaller irregular objects, in particular given strong indications that turbulent pressure might be responsible for self-regulating star formation in these cases (Le Tiran et al. 2011; Lehnert et al. 2013).
4.3 The CO→H2 conversion factor
One of the key ingredients in determining the total molecular gas mass in galaxies is converting CO luminosities into molecular hydrogen masses. Although a typical value of αCO = 4.6 M⊙ (K km s−1 pc2)−1 appears appropriate for normal SFGs in the local universe, the picture is less clear regarding low-metallicity objects and compact SFGs at high redshift.
On one hand, Leroy et al. (2011) have found that, in the Local Group, low-metallicity galaxies such as the Small Magellanic Cloud show αCO values as high as 70 – at this level, CO is no longer self-shielded from ionizing radiation in the galaxy, while hydrogen is still present in molecular form; therefore CO emission is much fainter for the same amount of gas (Wolfire, Hollenbach & McKee 2010). On the other hand, dusty starburst galaxies typically show αCO values closer to 0.8, probably due to turbulence in the ISM (Downes & Solomon 1998; Papadopoulos et al. 2012). For a thorough review on the subject, we refer the reader to Bolatto et al. (2013).
This presents a particular conundrum for the study of molecular gas in LBAs; while in the local universe low-metallicity galaxies are typically dwarf irregulars and compact starbursts tend to be (dusty) LIRGs or ULIRGs, many LBAs can be classified both as low-metallicity and compact, dense SFGs. On one hand, as discussed in Section 2.2, many LBAs – and 4 out of the 5 objects presented here, see Fig. 3 – are LIRGs, with IR luminosities above LFIR > 1011 L⊙, and all of them would be considered as ‘starburst’ galaxies in the local universe. Narayanan et al. (2012) have shown that increased gas temperatures and high velocity dispersions, both of which are clearly present in LBAs, could have the effect of decreasing the αCO factor to values more similar to infrared-bright merger-dominated starbursts. On the other hand, LBAs present lower metallicities and extinction values than galaxies within the same stellar mass range, and that would drive αCO values up. We emphasize that this combination of physical conditions of the ISM should be prevalent in SFGs in earlier epochs.
This determination is particularly important because the bimodality observed in the S-K relation can be attributed, at least in part, to the use of different αCO values for dusty and UV bright samples. In the case of high-redshift galaxies, direct measurements of the conversion factor are nearly impossible, and one has to rely on empirical results for low-redshift galaxies, which present dramatically different ISM properties and might misrepresent the problem.
To illustrate this statement, we compared estimated values for αCO in all objects in our sample [estimating gas masses from equation (4) and CO luminosities from equation (6); Fig. 7]. We converted between metallicities following the empirical relation found in Kewley & Ellison (2008). We then compared these estimates with values found in Genzel et al. (2012) and Magdis et al. (2012b). At high metallicities there is good agreement between our sample and the high-redshift one; however, at low metallicities, our inferred αCO are smaller than local galaxies with similar metallicities almost by a factor of 5. For galaxies observed thus far with CARMA (red points), most of which are more metal rich, we conclude that the Galactic value for |$\alpha _{\text{CO}}$| is appropriate. The estimates for lower metallicity galaxies still carry considerable uncertainties, given that we have not directly measured CO luminosities in these objects yet.
Another way to look at the problem is pictured in Fig. 8, where we show values of |$L^{\prime }_{\rm CO}$| versus LFIR, as in Fig. 3, but this time we infer CO luminosities assuming an αCO conversion factor following equation (11). We notice that the inferred CO luminosities are much lower than expected for a given FIR luminosity, because at lower luminosities (i.e. lower metallicities) the same gas masses would produce much fainter CO emission. The conclusion is that one of three assumptions must be wrong: either (1) the |$L^{\prime }_{\rm CO}-L_{\rm FIR}$| is not valid for LBAs; (2) the S-K relation is not valid for LBAs; or (3) the |$\alpha _{\rm CO}-Z$| relation is not valid for LBAs. Atacama Large Millimeter/submillimeter Array (ALMA) should be able to provide a definite answer to this question, since measurement of CO(1–0) with the finished array is achievable after less than an hour of integration even for the least massive, most metal-poor objects in our sample. The study of low-mass dust-deficient galaxies is of particular relevance, since these objects are believed to dominate the star formation density of the universe at high redshifts (Reddy & Steidel 2009).
5 SUMMARY
We have initiated a survey measuring CO(1–0) emission in LBAs, in an attempt to estimate gas fractions in these galaxies and obtain further insight towards the relation between star formation and the cold gas reservoir in UV-bright SFGs. We discuss expectations using well-established relations for the local universe (such as the S-K relation) and observed surface densities and gas fractions of high-redshift SFGs.
Using the CARMA interferometer, we have detected five objects at high S/N. We find CO luminosities that agree well with expectations from LFIR (Fig. 3). These galaxies show strong emission, indicative of high gas masses and gas fractions, between 20 and 60 per cent (Fig. 4). They also show shorter gas consumption time-scales than typically found in SFGs in the local universe (between 0.1 and 1 Gyr). These values agree with our expectations, further supporting our use of LBAs as low-redshift proxies for exploring the connections between star formation, molecular gas and triggering mechanism in sources that are more similar to typical high-redshift galaxies than any other local samples studied thus far.
Furthermore, we have shown that LBAs follow the local S-K relation, albeit at much higher surface densities than typical spirals at z ∼ 0 (Fig. 5). This is in accordance with other high-redshift populations (Daddi et al. 2010b), and the result distinguishes our sample from infrared-bright objects such as ULIRGs and SMGs. The authors in that work have argued that the bimodality can be removed if one takes into account the dynamical times of each object. This is still uncertain in the case of our galaxies; there is an apparent shift to smaller SFR surface densities even when taking into account dynamical times, but this may be the result of the use of rotational periods as the dynamical time-scales of galaxies.
We note the inherent difficulties in determining gas masses from CO luminosities, in particular given the uncertainties in |$\alpha _{\text{CO}}$| (especially at low metallicities). Given that low- and intermediate-mass galaxies will be more metal poor than similar objects in the local universe, this could strongly affect our understanding of molecular gas at earlier epochs. Our LBA sample, on the other hand, offers an excellent opportunity to study molecular gas in low metallicity, compact SFGs and ALMA observations will allow for accurate estimates of αCO in such objects.
We thank the anonymous referee for suggestions that helped improve this paper. We would also like to thank Andrew Baker and Tim Heckman for useful comments. TSG gratefully acknowledges CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) for financial support. Support for CARMA construction was derived from the states of California, Illinois, and Maryland, the James S. McDonnell Foundation, the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the University of Chicago, the Associates of the California Institute of Technology, and the National Science Foundation. Ongoing CARMA development and operations are supported by the National Science Foundation under a cooperative agreement, and by the CARMA partner universities.
Jansky Fellow.
The package is available at http://bima.astro.umd.edu/miriad/
LFIR is defined as the luminosity integrated between 40 and 120 μm, and is typically a factor of 1.7 smaller than total infrared luminosity LIR(3–1000 μm).