sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Sargent, M.

Normalized to: Sargent, M.

133 article(s) in total. 728 co-authors, from 1 to 66 common article(s). Median position in authors list is 9,0.

[1]  [pdf] - 2052172
The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$
Comments: 32 pages, 13 figures, Accepted for publication in ApJ
Submitted: 2020-02-20
Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.
[2]  [pdf] - 2065422
ALMA reveals the molecular gas properties of 5 star-forming galaxies across the main sequence at 3 < z < 3.5
Comments: 17 pages, 8 figures; accepted for publication in The Astrophysical Journal
Submitted: 2020-02-10
We present the detection of CO(5-4) with S/N> 7 - 13 and a lower CO transition with S/N > 3 (CO(4-3) for 4 galaxies, and CO(3-2) for one) with ALMA in band 3 and 4 in five main sequence star-forming galaxies with stellar masses 3-6x10^10 M/M_sun at 3 < z < 3.5. We find a good correlation between the total far-infrared luminosity LFIR and the luminosity of the CO(5-4) transition L'CO(5-4), where L'CO(5-4) increases with SFR, indicating that CO(5-4) is a good tracer of the obscured SFR in these galaxies. The two galaxies that lie closer to the star-forming main sequence have CO SLED slopes that are comparable to other star-forming populations, such as local SMGs and BzK star-forming galaxies; the three objects with higher specific star formation rates (sSFR) have far steeper CO SLEDs, which possibly indicates a more concentrated episode of star formation. By exploiting the CO SLED slopes to extrapolate the luminosity of the CO(1-0) transition, and using a classical conversion factor for main sequence galaxies of alpha_CO = 3.8 M_sun(K km s^-1 pc^-2)^-1, we find that these galaxies are very gas rich, with molecular gas fractions between 60 and 80%, and quite long depletion times, between 0.2 and 1 Gyr. Finally, we obtain dynamical masses that are comparable with the sum of stellar and gas mass (at least for four out of five galaxies), allowing us to put a first constraint on the alpha_CO parameter for main sequence galaxies at an unprecedented redshift.
[3]  [pdf] - 2061726
Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$
Comments: Comments: 6 pages, 4 figures, accepted for publication in Astronomy and Astrophysics Letters on February 4th 2020
Submitted: 2020-02-07
We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.
[4]  [pdf] - 2034616
The ALMaQUEST Survey: II. What drives central starbursts at z~0?
Comments: Accepted for publication in MNRAS
Submitted: 2020-01-14
Starburst galaxies have elevated star formation rates (SFRs) for their stellar mass. In Ellison et al. (2018) we used integral field unit (IFU) maps of star formation rate surface density (Sigma_SFR) and stellar mass surface density (Sigma_*) to show that starburst galaxies in the local universe are driven by SFRs that are preferentially boosted in their central regions. Here, we present molecular gas maps obtained with the Atacama Large Millimeter Array (ALMA) observatory for 12 central starburst galaxies at z~0 drawn from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The ALMA and MaNGA data are well matched in spatial resolution, such that the ALMA maps of molecular gas surface density (Sigma_H2) can be directly compared with MaNGA maps at kpc-scale resolution. The combination of Sigma_H2, Sigma_* and Sigma_SFR at the same resolution allow us to investigate whether central starbursts are driven primarily by enhancements in star formation efficiency (SFE) or by increased gas fractions. By computing offsets from the resolved Kennicutt-Schmidt relation (Sigma_H2 vs. Sigma_SFR) and the molecular gas main sequence (Sigma_* vs. Sigma_H2), we conclude that the primary driver of the central starburst is an elevated SFE. We also show that the enhancement in Sigma_SFR is accompanied by a dilution in O/H, consistent with a triggering that is induced by metal poor gas inflow. These observational signatures are found in both undisturbed (9/12 galaxies in our sample) and recently merged galaxies, indicating that both interactions and secular mechanisms contribute to central starbursts.
[5]  [pdf] - 2057862
A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48
Comments: 14 pages, 15 figures, accepted for publication in A&A
Submitted: 2020-01-10
[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.
[6]  [pdf] - 2038358
The ALMaQUEST Survey: III. Scatter in the resolved star forming main sequence is primarily due to variations in star formation efficiency
Comments: Accepted for publication in MNRAS
Submitted: 2019-11-26
Using a sample of 11,478 spaxels in 34 galaxies with molecular gas, star formation and stellar maps taken from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the parameters that correlate with variations in star formation rates on kpc scales. We use a combination of correlation statistics and an artificial neural network to quantify the parameters that drive both the absolute star formation rate surface density (Sigma_SFR), as well as its scatter around the resolved star forming main sequence (Delta Sigma_SFR). We find that Sigma_SFR is primarily regulated by molecular gas surface density (Sigma_H2) with a secondary dependence on stellar mass surface density (Sigma_*), as expected from an `extended Kennicutt-Schmidt relation'. However, Delta Sigma_SFR is driven primarily by changes in star formation efficiency (SFE), with variations in gas fraction playing a secondary role. Taken together, our results demonstrate that whilst the absolute rate of star formation is primarily set by the amount of molecular gas, the variation of star formation rate above and below the resolved star forming main sequence (on kpc scales) is primarily due to changes in SFE.
[7]  [pdf] - 1987505
Automated Mining of the ALMA Archive in the COSMOS Field (A3COSMOS): I. Robust ALMA Continuum Photometry Catalogs and Stellar Mass and Star Formation Properties for ~700 Galaxies at z=0.5-6
Comments: 56 pages, 36 figures, 5 tables, published on ApJS. Catalogs are available at journal website in FITS format. Current and future catalogs and continuum images will be released at COSMOS IRSA server and A3COSMOS website
Submitted: 2019-10-28
The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.
[8]  [pdf] - 2026167
A$^3$COSMOS: The Dust Attenuation of Star-Forming Galaxies at $z=2.5-4.0$ from the COSMOS-ALMA Archive
Comments: 12 pages, 5 figures, MNRAS submitted
Submitted: 2019-10-28
We present an analysis of the dust attenuation of star forming galaxies at $z=2.5-4.0$ through the relationship between the UV spectral slope ($\beta$), stellar mass ($M_{\ast}$) and the infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$) based on far-infrared continuum observations from the Atacama Large Millimeter/sub-millimeter Array (ALMA). Our study exploits the full ALMA archive over the COSMOS field processed by the A$^3$COSMOS team, which includes an unprecedented sample of $\sim1500$ galaxies at $z\sim3$ as primary or secondary targets in ALMA band 6 or 7 observations with a median continuum sensitivity of 126 $\rm{\mu Jy/beam}$ (1$\sigma$). The detection rate is highly mass dependent, decreasing drastically below $\log (M_{\ast}/M_{\odot})=10.5$. The detected galaxies show that the IRX-$\beta$ relationship of massive ($\log M_{\ast}/M_{\odot} > 10$) main sequence galaxies at $z=2.5-4.0$ is consistent with that of local galaxies, while starbursts are generally offset by $\sim0.5\,{\rm dex}$ to larger IRX values. At the low mass end, we derive upper limits on the infrared luminosities through stacking of the ALMA data. The combined IRX-$M_{\ast}$ relation at $\rm{log\,(M_{\ast}/M_{\odot})>9}$ exhibits a significantly steeper slope than reported in previous studies at similar redshifts, implying little dust obscuration at $\log M_{\ast}/M_{\odot}<10$. However, our results are consistent with early measurements at $z\sim5.5$, indicating a potential redshift evolution between $z\sim2$ and $z\sim6$. Deeper observations targeting low mass galaxies will be required to confirm this finding.
[9]  [pdf] - 2026166
Automated Mining of the ALMA Archive in the COSMOS Field (A3COSMOS): II. Cold Molecular Gas Evolution out to Redshift 6
Comments: 41 pages, 23 figures, 3 tables, resubmitted to ApJ after addressing referee's comments. Meanwhile we provide a Python package to calculate our gas evolution function ( )
Submitted: 2019-10-28
We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.
[10]  [pdf] - 1974205
Morphological evolution and galactic sizes in the L-Galaxies SA model
Submitted: 2018-10-11, last modified: 2019-10-07
In this work we update the L-Galaxies semi-analytic model (SAM) to better follow the physical processes responsible for the growth of bulges via disc instabilities (leading to pseudo-bulges) and mergers (leading to classical bulges). We address the former by considering the contribution of both stellar and gaseous discs in the stability of the galaxy, and we update the latter by including dissipation of energy in gas-rich mergers. Furthermore, we introduce angular momentum losses during cooling and find that an accurate match to the observed correlation between stellar disc scale length and mass at z ~ 0.0 requires that the gas loses 20% of its initial specific angular momentum to the corresponding dark matter halo during the formation of the cold gas disc. We reproduce the observed trends between the stellar mass and specific angular momentum for both disc- and bulge-dominated galaxies, with the former rotating faster than the latter of the same mass. We conclude that a two-component instability recipe provides a morphologically diverse galaxy sample which matches the observed fractional breakdown of galaxies into different morphological types. This recipe also enables us to obtain an excellent fit to the morphology-mass relation and stellar mass function of different galactic types. Finally, we find that energy dissipation during mergers reduces the merger remnant sizes and allows us to match the observed mass-size relation for bulge-dominated systems.
[11]  [pdf] - 1962379
JINGLE V: Dust properties of nearby galaxies derived from hierarchical Bayesian SED fitting
Comments: MNRAS in press, 31 pages
Submitted: 2019-09-11
We study the dust properties of 192 nearby galaxies from the JINGLE survey using photometric data in the 22-850micron range. We derive the total dust mass, temperature T and emissivity index beta of the galaxies through the fitting of their spectral energy distribution (SED) using a single modified black-body model (SMBB). We apply a hierarchical Bayesian approach that reduces the known degeneracy between T and beta. Applying the hierarchical approach, the strength of the T-beta anti-correlation is reduced from a Pearson correlation coefficient R=-0.79 to R=-0.52. For the JINGLE galaxies we measure dust temperatures in the range 17-30 K and dust emissivity indices beta in the range 0.6-2.2. We compare the SMBB model with the broken emissivity modified black-body (BMBB) and the two modified black-bodies (TMBB) models. The results derived with the SMBB and TMBB are in good agreement, thus applying the SMBB, which comes with fewer free parameters, does not penalize the measurement of the cold dust properties in the JINGLE sample. We investigate the relation between T and beta and other global galaxy properties in the JINGLE and Herschel Reference Survey (HRS) sample. We find that beta correlates with the stellar mass surface density (R=0.62) and anti-correlates with the HI mass fraction (M(HI)/M*, R=-0.65), whereas the dust temperature correlates strongly with the SFR normalized by the dust mass (R=0.73). These relations can be used to estimate T and beta in galaxies with insufficient photometric data available to measure them directly through SED fitting.
[12]  [pdf] - 1962324
Sunyaev-Zel'dovich detection of the galaxy cluster Cl J1449+0856 at z = 1.99: the pressure profile in uv space
Comments: 11 pages, 6 figures (+appendices); accepted for publication in Astronomy & Astrophysics
Submitted: 2019-07-25
We present Atacama Large Millimetre Array and Atacama Compact Array observations of the Sunyaev-Zel'dovich effect in the z = 2 galaxy cluster Cl J1449+0856, an X-ray-detected progenitor of typical massive clusters in the present day Universe. While in a cleaned but otherwise untouched 92 GHz map of this cluster, little to no negative signal is visible, careful subtraction of known sub-millimetre emitters in the uv plane reveals a decrement at 5$\sigma$ significance. The total signal is -190$\pm$36 $\mu$Jy, with a peak offset by 5"-9" ($\sim$50 kpc) from both the X-ray centroid and the still-forming brightest cluster galaxy. A comparison of the recovered uv-amplitude profile of the decrement with different pressure models allows us to derive total mass constraints consistent with the $\sim$6$\times$10$^{13}$ M$_{\odot}$ estimated from X-ray data. Moreover, we find no strong evidence for a deviation of the pressure profile with respect to local galaxy clusters, although a slight tension at small-to-intermediate spatial scales suggests a flattened central profile, opposite to what seen in a cool core and possibly an AGN-related effect. This analysis of the lowest mass single SZ detection so far illustrates the importance of interferometers when observing the SZ effect in high-redshift clusters, the cores of which cannot be considered quiescent, such that careful subtraction of galaxy emission is necessary.
[13]  [pdf] - 1900837
Early- and late-stage mergers among main sequence and starburst galaxies at 0.2<z<2
Comments: Published version
Submitted: 2018-09-03, last modified: 2019-06-16
We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2$\leq z\leq$2.0. The novelty of our work lies in the use of a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to $\geq$70% above the MS, implying that -- already at $z{\gtrsim}1$ -- starburst (SB) events ($\Delta_{\rm MS}\geq$0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS-offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of $\sim$5-10% at all masses ${\gtrsim} 10^{10}M_{\odot}$. The MS major merger fraction roughly doubles between $z=0.2$ and $z=2$, with morphological mergers driving the overall increase at $z{\gtrsim}1$. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability timescales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with $\Delta_{\rm MS}$ follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different SFR enhancements. This points to a starburst-merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.
[14]  [pdf] - 1894355
The Main Sequence at $z \sim 1.3$ contains a sizable fraction of galaxies with compact star formation sizes: a new population of early post-starbursts?
Comments: 7 pages, 5 figures. Accepted for publication in ApJ Letters
Submitted: 2019-05-08
ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$.
[15]  [pdf] - 1875477
JINGLE, a JCMT legacy survey of dust and gas for galaxy evolution studies: II. SCUBA-2 850 {\mu}m data reduction and dust flux density catalogues
Comments: Accepted by MNRAS; data available at
Submitted: 2019-04-23
We present the SCUBA-2 850 ${\mu}m$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 ${\mu}m$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $\beta$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample.
[16]  [pdf] - 1890483
The IRAM/GISMO two-millimeter survey in the COSMOS field
Comments: 27 pages, 14 figures, accepted for publication in ApJ
Submitted: 2019-04-22
We present deep continuum observations at a wavelength of 2mm centered on the COSMOS field using the Goddard IRAM Superconducting Millimeter Observer (GISMO) at the IRAM 30m-telescope. These data constitute the widest deep 2mm survey to-date, reaching a uniform $\sigma\sim0.23$mJy beam$^{-1}$ sensitivity over $\sim250$arcmin$^2$ at $\sim\,24''$ resolution. We detect four sources at high significance (S/N$\,\geq\,$4.4) with an expected number of false detection of 0.09 sources, and five sources at $4.4>\,$S/N$\,\geq\,3.7$ with an expected number of false detection of 1.65 sources. Combined with deep GISMO observations in GOODS-N, we constrain the 2mm number counts over one decade in flux density. These measurements agree with most galaxy evolution models tested here, except those with large population of dusty star-forming galaxies at $z>7$. Five GISMO sources have counterparts in (sub-)millimeter catalogs available in COSMOS. Their redshifts suggest that all but one lie above $z\sim3$. These four high-redshift ($z>3$) galaxies have $\tilde{z}=3.9$, SFRs $\sim\,$400 - 1200M$_{\odot}\,$yr$^{-1}$ and $M_{\rm dust}\sim10^{9.5}\,$M$_{\odot}$. They provide a relatively complete selection ($\sim66\%$) of the most luminous ($L_{\rm IR}>10^{12.6}\,$L$_{\odot}$) and highest redshift ($z>3$) galaxies detected within our survey area by AzTEC at 1.1mm. We thus conclude that 2mm surveys favor the selection of massive, vigorously star-forming, high-redshift galaxies. This is corroborated by GISMO-C4, a source with a low false detection probability ($\sim\,6.2\%$), for which the absence of a (sub-)millimeter counterpart supports a high redshift origin ($z\gtrsim3$).
[17]  [pdf] - 1886566
Radio continuum size evolution of star-forming galaxies over 0.35 < z < 2.25
Comments: 21 pages, 19 figures, accepted for publication in A&A
Submitted: 2019-03-28
We present the first systematic study of the radio continuum size evolution of star-forming galaxies (SFGs) over the redshift range $0.35<z<2.25$. We use the VLA COSMOS 3GHz map (noise $\rm rms=2.3\,\mu Jy \,beam^{-1}$, $\theta_{\rm beam}=0.75\,\rm arcsec$) to construct a mass-complete sample of 3184 radio-selected SFGs that reside on and above the main-sequence (MS) of SFGs. We find no clear dependence between the radio size and stellar mass, $M_{\star}$, of SFGs with $10.5\lesssim\log(M_\star/\rm M_\odot)\lesssim11.5$. Our analysis suggests that MS galaxies are preferentially extended, while SFGs above the MS are always compact. The median effective radius of SFGs on (above) the MS of $R_{\rm eff}=1.5\pm0.2$ ($1.0\pm0.2$) kpc remains nearly constant with cosmic time; a parametrization of the form $R_{\rm eff}\propto(1+z)^\alpha$ yields a shallow slope of only $\alpha=-0.26\pm0.08\,(0.12\pm0.14)$ for SFGs on (above) the MS. The size of the stellar component of galaxies is larger than the extent of the radio continuum emission by a factor $\sim$2 (1.3) at $z=0.5\,(2)$, indicating star formation is enhanced at small radii. The galactic-averaged star formation rate surface density $(\Sigma_{\rm SFR})$ scales with the distance to the MS, except for a fraction of MS galaxies ($\lesssim10\%$) that harbor starburst-like $\Sigma_{\rm SFR}$. These "hidden" starbursts might have experienced a compaction phase due to disk instability and/or merger-driven burst of star formation, which may or may not significantly offset a galaxy from the MS. We thus propose to jointly use $\Sigma_{\rm SFR}$ and distance to the MS to better identify the galaxy population undergoing a starbursting phase.
[18]  [pdf] - 1966676
The ALMA Spectroscopic Survey in the HUDF: CO luminosity functions and the molecular gas content of galaxies through cosmic history
Comments: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee
Submitted: 2019-03-21
We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).
[19]  [pdf] - 1853644
The evolution of the cosmic molecular gas density
Comments: 7 pages, 2 figures, Science White paper submitted to Astro2020 Decadal Survey
Submitted: 2019-03-20
One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.
[20]  [pdf] - 1848626
Simultaneous Measurements of Star Formation and Supermassive Black Hole Growth in Galaxies
Comments: Science White paper submitted to Astro2020 Decadal Survey
Submitted: 2019-03-12
Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.
[21]  [pdf] - 1916806
A closer look at the deep radio sky: Multi-component radio sources at 3-GHz VLA-COSMOS
Comments: A&A accepted; 32 pages, 14 figures, 5 tables, 3 Appendices
Submitted: 2019-01-29
In this data paper we present and characterise the multi-component radio sources identified in the VLA-COSMOS Large Project at 3 GHz (0.75 arcsec resolution, 2.3 {\mu}Jy/beam rms), i.e. the radio sources which are composed of two or more radio blobs.The classification of objects into multi-components was done by visual inspection of 351 of the brightest and most extended blobs from a sample of 10,899 blobs identified by the automatic code blobcat. For that purpose we used multi-wavelength information of the field, such as the 1.4-GHz VLA-COSMOS data and the UltraVISTA stacked mosaic available for COSMOS. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with AGN powered radio emission and 9 star-forming galaxies. We report 8 new detections that were not observed by the VLA-COSMOS Large Project at 1.4 GHz, due to the slightly larger area coverage at 3 GHz. The increased spatial resolution of 0.75 arcsec has allowed us to resolve (and isolate) multiple emission peaks of 28 extended radio sources not identified in the 1.4-GHz VLA-COSMOS map. We report the multi-frequency flux densities (324 MHz, 325 MHz, 1.4 GHz & 3 GHz), star-formation-rates, and stellar masses of these objects. Multi-component objects at 3-GHz VLA-COSMOS inhabit mainly massive galaxies (>10^10.5 Msun). The majority of the multi-component AGN lie below the main-sequence of star-forming galaxies (SFGs), in the green valley and the quiescent region. We provide detailed description of the objects: amongst the AGN there are 2 head-tail, 10 core-lobe, 9 wide-angle-tail (WAT), 8 double-double or Z-/X-shaped, 3 bent-tail radio sources, and 26 symmetric sources, while amongst the SFGs we find the only star-forming ring seen in radio emission in COSMOS. We report a large number (32/58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (0.08 < z < 1.53). [abridged]
[22]  [pdf] - 1834252
Suppressed CO emission and high G/D ratios in z=2 galaxies with sub-solar gas-phase metallicity
Comments: 15 pages, 5 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2019-01-18
We study a population of significantly sub-solar enrichment galaxies at z=1.99, to investigate how molecular gas, dust and star-formation relate in low-metallicity galaxies at the peak epoch of star-formation. We target our sample with several deep ALMA and VLA datasets, and find no individual detections of CO[4-3], CO[1-0] or dust, in stark contrast to the $>$60% detection rate expected for solar-enrichment galaxies with these MS Halpha SFRs. We find that both low and high density molecular gas (traced by CO[1-0] and CO[4-3] respectively) are affected by the low enrichment, showing sample average (stacked) luminosity deficits $>$0.5-0.7 dex below expectations. This is particularly pertinent for the use of high-J CO emission as a proxy of instantaneous star-formation rate. Our individual galaxy data and stacked constraints point to a strong inverse dependence $\propto Z^{\gamma}$ of gas-to-dust ratios (G/D) and CO-to-H$_{2}$ conversion factors (aco) on metallicity at z$\sim$2, with $\gamma_{\rm G/D}<$-2.2 and $\gamma_{\alpha_{\rm CO}}<$-0.8, respectively. We quantify the importance of comparing G/D and aco vs. metallicity trends from the literature on a common, suitably normalised metallicity scale. When accounting for systematic offsets between different metallicity scales, our z$\sim$2 constraints on these scaling relations are consistent with the corresponding relations for local galaxies. However, among those local relations, we favour those with a steep/double power-law dependence of G/D on metallicity. Finally, we discuss the implications of these findings for (a) gas mass measurements for sub-M$^{*}$ galaxies, and (b) efforts to identify the characteristic galaxy mass scale contributing most to the comoving molecular gas density at z=2.
[23]  [pdf] - 1830437
COLDz: Shape of the CO Luminosity Function at High Redshift and the Cold Gas History of the Universe
Comments: 16 pages, 6 figures, 6 tables, to appear in ApJ (accepted January 04, 2019)
Submitted: 2018-08-13, last modified: 2019-01-07
We report the first detailed measurement of the shape of the CO luminosity function at high redshift, based on $>$320 hr of the NSF's Karl G. Jansky Very Large Array (VLA) observations over an area of $\sim$60 arcmin$^2$ taken as part of the CO Luminosity Density at High Redshift (COLDz) survey. COLDz "blindly" selects galaxies based on their cold gas content through CO($J$=1$\to$0) emission at $z$$\sim$2-3 and CO($J$=2$\to$1) at $z$$\sim$5-7 down to a CO luminosity limit of log($L'_{\rm CO}$/K km s$^{-1}$ pc$^2$)$\simeq$9.5. We find that the characteristic luminosity and bright end of the CO luminosity function are substantially higher than predicted by semi-analytical models, but consistent with empirical estimates based on the infrared luminosity function at $z$$\sim$2. We also present the currently most reliable measurement of the cosmic density of cold gas in galaxies at early epochs, i.e., the cold gas history of the universe, as determined over a large cosmic volume of $\sim$375,000 Mpc$^3$. Our measurements are in agreement with an increase of the cold gas density from $z$$\sim$0 to $z$$\sim$2-3, followed by a possible decline towards $z$$\sim$5-7. These findings are consistent with recent surveys based on higher-$J$ CO line measurements, upon which COLDz improves in terms of statistical uncertainties by probing $\sim$50-100 times larger areas and in the reliability of total gas mass estimates by probing the low-$J$ CO lines accessible to the VLA. Our results thus appear to suggest that the cosmic star-formation rate density follows an increased cold molecular gas content in galaxies towards its peak about 10 billion years ago, and that its decline towards the earliest epochs is likely related to a lower overall amount of cold molecular gas (as traced by CO) bound in galaxies towards the first billion years after the Big Bang.
[24]  [pdf] - 1783975
Revealing environmental dependence of molecular gas content in a distant X-ray cluster at z=2.51
Comments: accepted for publication in ApJ Letters
Submitted: 2018-10-24
We present a census of the molecular gas properties of galaxies in the most distant known X-ray cluster, CLJ1001, at z=2.51, using deep observations of CO(1-0) with JVLA. In total 14 cluster members with $M_{*} > 10^{10.5} M_{\odot}$ are detected, including all the massive star-forming members within the virial radius, providing the largest galaxy sample in a single cluster at $z > 2$ with CO(1-0) measurements. We find a large variety in the gas content of these cluster galaxies, which is correlated with their relative positions (or accretion states), with those closer to the cluster core being increasingly gas-poor. Moreover, despite their low gas content, the galaxies in the cluster center exhibit an elevated star formation efficiency (SFE=SFR/$M_{\rm gas}$) compared to field galaxies, suggesting that the suppression on the SFR is likely delayed compared to that on the gas content. Their gas depletion time is around $t_{\rm dep} \sim 400$ Myrs, comparable to the cluster dynamical time. This implies that they will likely consume all their gas within a single orbit around the cluster center, and form a passive cluster core by $z\sim2$. This result is one of the first direct pieces of evidence for the influence of environment on the gas reservoirs and SFE of $z > 2$ cluster galaxies, thereby providing new insights into the rapid formation and quenching of the most massive galaxies in the early universe.
[25]  [pdf] - 1795845
The Effect of Galaxy Interactions on Molecular Gas Properties
Comments: 26 pages, 11 figures, accepted for publication in ApJ
Submitted: 2018-10-23
Galaxy interactions are often accompanied by an enhanced star formation rate (SFR). Since molecular gas is essential for star formation, it is vital to establish whether, and by how much, galaxy interactions affect the molecular gas properties. We investigate the effect of interactions on global molecular gas properties by studying a sample of 58 galaxies in pairs and 154 control galaxies. Molecular gas properties are determined from observations with the JCMT, PMO, CSO telescopes, and supplemented with data from the xCOLD GASS and JINGLE surveys at $^{12}$CO(1-0) and $^{12}$CO(2-1). The SFR, gas mass ($M_\mathrm{H_{2}}$), and gas fraction ($f_{gas}$) are all enhanced in galaxies in pairs by $\sim$ 2.5 times compared to the controls matched in redshift, mass, and effective radius, while the enhancement of star formation efficiency (SFE $\equiv$ SFR/$M_{H_{2}}$) is less than a factor of 2. We also find that the enhancements in SFR, $M_{H_{2}}$ and $f_{gas}$ increase with decreasing pair separation and are larger in systems with smaller stellar mass ratio. Conversely, the SFE is only enhanced in close pairs (separation $<$ 20 kpc) and equal-mass systems; therefore most galaxies in pairs lie in the same parameter space on the SFR-$M_{H_{2}}$ plane as controls. This is the first time that the dependence of molecular gas properties on merger configurations is probed statistically with a relatively large sample and with a carefully-selected control sample for individual galaxies. We conclude that galaxy interactions do modify the molecular gas properties, although the strength of the effect is merger configuration dependent.
[26]  [pdf] - 1767845
Science with an ngVLA. Cold gas in High-z Galaxies: CO as redshift beacon
Comments: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-16
The goal of this science case is to address the use of a ngVLA as a CO redshift machine for dust-obscured high-redshift galaxies which lack of clear counterparts at other wavelengths. Thanks to its unprecedentedly large simultaneous bandwidth and sensitivity, the ngVLA will be able to detect low--J CO transitions at virtually any $z>1$. In particular, at $z>4.76$ two CO transitions will be covered in a single frequency setting, thus ensuring unambiguous line identification. The ngVLA capabilities fill in a redshift range where other approaches (e.g., photometric redshifts, search for optical/radio counterparts, etc) typically fail due to the combination of intrinsically faint emission and increasing luminosity distance. This will allow us to explore the formation of massive galaxies in the early cosmic times.
[27]  [pdf] - 1769133
Science with an ngVLA. Cold gas in High-z Galaxies: The dense ISM
Comments: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: text overlap with arXiv:1510.06411
Submitted: 2018-10-16
The goal of this science case is to study physical conditions of the interstellar medium (ISM) in distant galaxies. In particular, its densest component is associated with the inner cores of clouds -- this is where star formation takes place. Carbon monoxide is usually used to trace molecular gas emission; however, its transitions are practically opaque, thus preventing astronomers from piercing through the clouds, into the deepest layers that are most intimately connected with the formation of stars. Other dense gas tracers are required, although they are typically too faint and/or at too low frequencies to be effectively observed in high redshift galaxies. The ngVLA will offer for the first time the sensitivity at radio frequencies that is needed to target [CI]$_{1-0}$ (at $z>5$), as well as the ground transitions of dense gas tracers of the ISM such as HCN, HNC, HCO+ (at various redshifts $z>1$), beyond the tip of the iceberg of the hyper-luminous sources that could be studied up to now. These new tools will critically contribute to our understanding of the intimate interplay between gas clouds and star formation in different environments and cosmic epochs.
[28]  [pdf] - 1769132
Science with an ngVLA. Cold gas in High-z Galaxies: The molecular gas budget
Comments: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1510.06411
Submitted: 2018-10-16
The goal of this science case is to accurately pin down the molecular gas content of high redshift galaxies. By targeting the CO ground transition, we circumvent uncertainties related to CO excitation. The ngVLA can observe the CO(1-0) line at virtually any $z>1.5$, thus exposing the evolution of gaseous reservoirs from the earliest epochs down to the peak of the cosmic history of star formation. The order-of-magnitude improvement in the number of CO detections with respect to state-of-the-art observational campaigns will provide a unique insight on the evolution of galaxies through cosmic time.
[29]  [pdf] - 1767869
Science with an ngVLA: New Parameter Space for Deep Field Radio Continuum Surveys
Comments: 6 pages, 2 figures. To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-16
Deep continuum surveys at radio wavelengths are able to cover large areas, yield high angular resolution, and do not suffer from dust extinction, thus providing a robust way to measure the star formation history of the universe. However, with the current sensitivities of existing radio telescopes, it remains challenging to detect galaxies that dominate the cosmic star formation history even with extremely long integrations. With the ngVLA, a new portion of parameter space will be opened up for radio continuum surveys: deep (~200 nJy/bm), large-area (~1 deg^2), sub-arcsecond surveys at high frequencies (~8 GHz), where the observed radio emission from high-z galaxies should be dominated by free-free emission, providing a robust measurement for the star formation history of the universe. By being able to image the star formation activity, unbiased by dust, for a large, homogeneous sample of galaxies with a wide range of luminosities into the epoch of reionization, such surveys with the ngVLA will be highly complementary to those conducted by JWST, which will only be sensitive to un-obscured star formation.
[30]  [pdf] - 1769125
Science with an ngVLA: Radio Continuum Emission from Galaxies: An Accounting of Energetic Processes
Comments: 9 pages, 3 figures; To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1510.06431
Submitted: 2018-10-15
Radio continuum observations have proven to be a workhorse in our understanding of the star formation process (i.e., stellar birth and death) from galaxies both in the nearby universe and out to the highest redshifts. In this article we focus on how the ngVLA will transform our understanding of star formation by enabling one to map and decompose the radio continuum emission from large, heterogeneous samples of nearby galaxies on $\gtrsim 10$\,pc scales to conduct a proper accounting of the energetic processes powering it. At the discussed sensitivity and angular resolution, the ngVLA will simultaneously be able to create maps of current star formation activity at $\sim$100\,pc scales, as well as detect and characterize (e.g., size, spectral shape, density, etc.) discrete H{\sc ii} regions and supernova remnants on 10\,pc scales in galaxies out to the distance of the Virgo cluster. Their properties can then be used to see how they relate to the local and global ISM and star formation conditions. Such investigations are essential for understanding the astrophysics of high-$z$ measurements of galaxies, allowing for proper modeling of galaxy formation and evolution.
[31]  [pdf] - 1791097
Concurrent starbursts in molecular gas disks within a pair of colliding galaxies at z = 1.52
Comments: 21 pages, 16 figures, accepted to ApJ
Submitted: 2018-10-03
We report on the discovery of a merger-driven starburst at z = 1.52, PACS-787, based on high signal-to-noise ALMA observations. CO(5-4) and continuum emission (850um) at a spatial resolution of 0.3" reveal two compact (r_1/2 ~ 1 kpc) and interacting molecular gas disks at a separation of 8.6 kpc thus indicative of an early stage in a merger. With a SFR of 991 Msun/yr, this starburst event should occur closer to final coalescence, as usually seen in hydrodynamical simulations. From the CO size, inclination, and velocity profile for both disks, the dynamical mass is calculated through a novel method that incorporates a calibration using simulations of galaxy mergers. Based on the dynamical mass, we measure (1) the molecular gas mass, independent from the CO luminosity, (2) the ratio of the total gas mass and the CO(1 - 0) luminosity (alpha_CO = M_gas/L'_CO(1-0)), and (3) the gas-to-dust ratio, with the latter two being lower than typically assumed. We find that the high star formation, triggered in both galaxies, is caused by a set of optimal conditions: a high gas mass/fraction, a short depletion time (t_depl=85 and 67 Myrs) to convert gas into stars, and the interaction of likely counter-rotating molecular disks that may accelerate the loss of angular momentum. The state of interaction is further established by the detection of diffuse CO and continuum emission, tidal debris that bridges the two nuclei and is associated with stellar emission seen by HST/WFC3. This observation demonstrates the power of ALMA to study the dynamics of galaxy mergers at high redshift.
[32]  [pdf] - 1783945
The molecular gas content and fuel efficiency of starbursts at z ~ 1.6 with ALMA
Comments: 16 pages, 8 figures
Submitted: 2018-10-03
We present an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed.
[33]  [pdf] - 1755149
SMBH accretion properties of radio-selected AGN out to z~4
Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in MNRAS
Submitted: 2018-09-17, last modified: 2018-09-24
Exploring how radio-emitting active galactic nuclei (AGN) behave and evolve with time is critical for understanding how AGN feedback impacts galaxy evolution. In this work, we investigate the relationship between 1.4 GHz radio continuum AGN luminosity ($L^{\rm AGN}_{\rm 1.4}$), specific black hole accretion rate (s-BHAR, defined as the accretion luminosity relative to the galaxy stellar mass) and redshift, for a luminosity-complete sample of radio-selected AGN in the VLA COSMOS 3 GHz Large Project. The sample was originally selected from radio-continuum observations at 3 GHz, and includes about 1800 radio AGN identified via ($>2\sigma$) radio-excess relative to the infrared-radio correlation of star-forming galaxies. We further select a subsample of over 1200 radio AGN that is complete in $L^{\rm AGN}_{\rm 1.4}$ over different redshift ranges, out to z~4, and use X-ray stacking to calculate the average s-BHAR in each $L^{\rm AGN}_{\rm 1.4}$-$z$ bin. We find that the average s-BHAR is independent of $L^{\rm AGN}_{\rm 1.4}$, at all redshifts. However, we see a strong increase of s-BHAR with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. This trend resembles the strong increase in the fraction of star-forming host galaxies (based on the $(NUV-r)$ / $(r-J)$ colours) with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. A possible explanation for this similarity might imply a link between average AGN radiative power and availability of cold gas supply within the host galaxy. This study corroborates the idea that radio-selected AGN become more radiatively efficient towards earlier epochs, independently of their radio power.
[34]  [pdf] - 1756083
JINGLE, a JCMT legacy survey of dust and gas for galaxy evolution studies: I. Survey overview and first results
Comments: MNRAS in press, 25 pages
Submitted: 2018-09-19
JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.
[35]  [pdf] - 1751828
Probing star formation and ISM properties using galaxy disk inclination II: Testing typical FUV attenuation corrections out to z$\sim$0.7
Comments: Accepted for publication in A&A (11/05/2018). Tables of the computed SFRs are only available at the CDS via anonymous ftp to ( or via
Submitted: 2018-05-17, last modified: 2018-09-12
We evaluate dust-corrected far ultraviolet (FUV) star formation rates (SFRs) for samples of star-forming galaxies at $z\sim0$ and $z\sim0.7$ and find significant differences between values obtained through corrections based on UV colour, from a hybrid mid-infrared (MIR) plus FUV relation, and from a radiative transfer based attenuation correction method. The performances of the attenuation correction methods are assessed by their ability to remove the dependency of the corrected SFR on inclination, as well as returning, on average, the expected population mean SFR. We find that combining MIR (rest-frame $\sim$13$\mu$m) and FUV luminosities gives the most inclination independent SFRs and reduces the intrinsic SFR scatter out of the methods tested. However, applying the radiative transfer based method of Tuffs et al. gives corrections to the FUV SFR that are inclination independent and in agreement with the expected SFRs at both $z\sim0$ and $z\sim0.7$. SFR corrections based on the UV-slope perform worse than the other two methods tested. For our local sample, the UV-slope method works on average but does not remove inclination biases. At $z\sim$0.7 we find that the UV-slope correction used locally flattens the inclination dependence compared to the raw FUV measurements but was not sufficient to correct for the large attenuation observed at $z\sim$0.7.
[36]  [pdf] - 1749245
Probing star formation and ISM properties using galaxy disk inclination I: Evolution in disk opacity since $z\sim0.7$
Comments: Accepted for publication in Section 4. Extragalactic astronomy of Astronomy and Astrophysics. 22/12/2017
Submitted: 2018-01-10, last modified: 2018-09-12
Disk galaxies at intermediate redshift ($z\sim0.7$) have been found in previous work to display more optically thick behaviour than their local counterparts in the rest-frame B-band surface brightness, suggesting an evolution in dust properties over the past $\sim$6 Gyr. We compare the measured luminosities of face-on and edge-on star-forming galaxies at different wavelengths (Ultraviolet (UV), mid-infrared (MIR), far-infrared (FIR), and radio) for two well-matched samples of disk-dominated galaxies: a local Sloan Digital Sky Survey (SDSS)-selected sample at $z\sim0.07$ and a sample of disks at $z\sim0.7$ drawn from Cosmic Evolution Survey (COSMOS). We have derived correction factors to account for the inclination dependence of the parameters used for sample selection. We find that typical galaxies are transparent at MIR wavelengths at both redshifts and that the FIR and radio emission is also transparent as expected. However, reduced sensitivity at these wavelengths limits our analysis; we cannot rule out opacity in the FIR or radio. Ultra-violet attenuation has increased between $z\sim0$ and $z\sim0.7$, with the $z\sim0.7$ sample being a factor of $\sim$3.4 more attenuated. The larger UV attenuation at $z\sim0.7$ can be explained by more clumpy dust around nascent star-forming regions. There is good agreement between the fitted evolution of the normalisation of the SFR$_{\text{UV}}$ versus 1-cos(i) trend (interpreted as the clumpiness fraction) and the molecular gas fraction/dust fraction evolution of galaxies found out to $z<1$.
[37]  [pdf] - 1748019
The [C II] emission as a molecular gas mass tracer in galaxies at low and high redshift
Comments: Accepted for publication in MNRAS. 28 pages, including 12 figures, 5 tables, 3 appendices
Submitted: 2018-08-30
We present ALMA Band 9 observations of the [C II]158um emission for a sample of 10 main-sequence galaxies at redshift z ~ 2, with typical stellar masses (log M*/Msun ~ 10.0 - 10.9) and star formation rates (~ 35 - 115 Msun/yr). Given the strong and well understood evolution of the interstellar medium from the present to z = 2, we investigate the behaviour of the [C II] emission and empirically identify its primary driver. We detect [C II] from six galaxies (four secure, two tentative) and estimate ensemble averages including non detections. The [C II]-to-infrared luminosity ratio (L[C II]/LIR) of our sample is similar to that of local main-sequence galaxies (~ 2 x 10^-3), and ~ 10 times higher than that of starbursts. The [C II] emission has an average spatial extent of 4 - 7 kpc, consistent with the optical size. Complementing our sample with literature data, we find that the [C II] luminosity correlates with galaxies' molecular gas mass, with a mean absolute deviation of 0.2 dex and without evident systematics: the [C II]-to-H2 conversion factor (alpha_[C II] ~ 30 Msun/Lsun) is largely independent of galaxies' depletion time, metallicity, and redshift. [C II] seems therefore a convenient tracer to estimate galaxies' molecular gas content regardless of their starburst or main-sequence nature, and extending to metal-poor galaxies at low and high redshifts. The dearth of [C II] emission reported for z > 6 - 7 galaxies might suggest either a high star formation efficiency or a small fraction of UV light from star formation reprocessed by dust.
[38]  [pdf] - 1743860
The CO Luminosity Density at High-z (COLDz) Survey: A Sensitive, Large Area Blind Search for Low-J CO Emission from Cold Gas in the Early Universe with the Karl G. Jansky Very Large Array
Comments: 53 pages, 22 figures, 8 tables, ApJ in press. Full resolution images and publicly available data may be found on the survey webpage at
Submitted: 2018-08-13
We describe the CO Luminosity Density at High-z (COLDz) survey, the first spectral line deep field targeting CO(1-0) emission from galaxies at $z=1.95-2.85$ and CO(2-1) at $z=4.91-6.70$. The main goal of COLDz is to constrain the cosmic density of molecular gas at the peak epoch of cosmic star formation. By targeting both a wide ($\sim$51 arcmin$^2$) and a deep area ($\sim$9 arcmin$^2$), the survey is designed to robustly constrain the bright end and the characteristic luminosity of the CO(1-0) luminosity function. An extensive analysis of the reliability of our line candidates, and new techniques provide detailed completeness and statistical corrections as necessary to determine the best constraints to date on the CO luminosity function. Our blind search for CO(1-0) uniformly selects starbursts and massive Main Sequence galaxies based on their cold molecular gas masses. Our search also detects CO(2-1) line emission from optically dark, dusty star-forming galaxies at $z>5$. We find a range of spatial sizes for the CO-traced gas reservoirs up to $\sim40$ kpc, suggesting that spatially extended cold molecular gas reservoirs may be common in massive, gas-rich galaxies at $z\sim2$. Through CO line stacking, we constrain the gas mass fraction in previously known typical star-forming galaxies at $z=2$-3. The stacked CO detection suggests lower molecular gas mass fractions than expected for massive Main Sequence galaxies by a factor of $\sim3-6$. We find total CO line brightness at $\sim34\,$GHz of $0.45\pm0.2\,\mu$K, which constrains future line intensity mapping and CMB experiments.
[39]  [pdf] - 1795680
GOODS-ALMA: 1.1 mm galaxy survey - I. Source catalogue and optically dark galaxies
Comments: 26 pages, 23 figures, accepted to A&A
Submitted: 2018-02-28, last modified: 2018-08-03
We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes
[40]  [pdf] - 1743794
"Super-deblended" Dust Emission in Galaxies: II. Far-IR to (sub)millimeter photometry and high redshift galaxy candidates in the full COSMOS field
Comments: 49 pages (+appendix), 48 figures, 3 tables, accepted for publication in ApJ
Submitted: 2018-07-12
We present a "super-deblended" far-infrared to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al. 2018, with key adaptations. We obtain point spread function (PSF) fitting photometry at fixed prior positions including 88,008 galaxies detected in either VLA 1.4~GHz, 3~GHz and/or MIPS 24~$\mu$m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed ``active' removal of non relevant priors at FIR/(sub)mm bands using spectral energy distribution (SED) fitting and redshift information. In order to cope with the shallower COSMOS data we subtract from the maps the flux of faint non-fitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer, Herschel, SCUBA2, AzTEC, MAMBO and NSF's Karl G. Jansky Very Large Array at 3~GHz and 1.4~GHz) displays well behaved quasi-Gaussian uncertainties, calibrated from Monte Carlo simulations and tailored to observables (crowding, residual maps). Comparison to ALMA photometry for hundreds of sources provide a remarkable validation of the technique. We detect 11,220 galaxies over the 100--1200~$\mu$m range, extending to $z_{\rm phot}\sim7$. We conservatively selected a sample of 85 $z>4$ high redshift candidates, significantly detected in the FIR/(sub)mm, often with secure radio and/or Spitzer/IRAC counterparts. This provides a chance to investigate the first generation of vigorous starburst galaxies (SFRs$\sim1000\mathrm{M}_\odot$~yr$^{-1}$). The photometric and value added catalogs are publicly released.
[41]  [pdf] - 1783817
The molecular gas content in obscured AGN at z > 1
Comments: 20 pages, 7 figures, accepted for publication in A&A
Submitted: 2018-07-09
The standard AGN-galaxy co-evolutionary scenario predicts a phase of deeply buried supermassive black hole growth coexisting with a starburst (SB) before feedback phenomena deplete the cold molecular gas reservoir of the galaxy and an optically luminous QSO is revealed ('SB-QSO evolutionary sequence'). The aim of this work is to measure the cold gas reservoir of three highly obscured QSOs to test if their gas fraction is similar to that of sub-millimeter galaxies (SMGs), as expected by some models, and place these measurements in the context of the SB-QSO framework. We target CO(1-0) transition in BzK4892, a Compton Thick (CT) QSO at z=2.6, CO(1-0) in BzK8608 and CO(2-1) in CDF153, two highly obscured QSOs at z=2.5 and z=1.5, respectively. For all these targets, we place 3$\sigma$ upper limits on the CO, with $L'_{CO} < (1.5\div 2.8)\times 10^{10}$ K km/s pc$^2$. We also compare the molecular gas conditions of our targets with those of other systems at z>1, considering normal star forming galaxies and SMGs, unobscured and obscured AGN from the literature. For the AGN samples, we provide an updated and (almost) complete collection of targets with CO follow-up. BzK4892 displays a high star formation efficiency (SFE$=L_{IR}/L'_{CO}>410$ L$_{\odot}$/(K km s$^{-1}$ pc$^2$)) and a gas fraction $f_{gas}<0.1$. Less stringent constraints are derived for the other two targets ($f_{gas}<0.5$ and SFE$>10$). From the comparison with literature data, we found that a) obscured AGN at z>1 are associated with higher SFE and lower $f_{gas}$ with respect to star forming galaxies; b) mildly and highly obscured active galaxies have comparable gas fractions; c) the SFE of CT and obscured AGN are similar to those of unobscured AGN. Within the SB-QSO framework, these findings could be consistent with a scenario where feedback can impact the host galaxy already from the early phases of the SB-QSO sequence.
[42]  [pdf] - 1729062
Deciphering the activity and quiescence of high-redshift cluster environments: ALMA observations of ClJ1449+0856 at z=2
Comments: 11 pages, 6 figures, ApJ in press
Submitted: 2018-06-18
We present ALMA observations of the 870$\mu$m continuum and CO(4-3) line emission in the core of the galaxy cluster ClJ1449+0856 at z=2, a NIR-selected, X-ray detected system in the mass range of typical progenitors of today's massive clusters. The 870$\mu$m map reveals six F$_{870\mu m}$ > 0.5 mJy sources spread over an area of 0.07 arcmin$^2$, giving an overdensity of a factor ~10 (6) with respect to blank field counts down to F$_{870\mu m}$ > 1 (0.5) mJy. On the other hand, deep CO(4-3) follow-up confirms membership of three of these sources, but suggests that the remaining three, including the brightest 870$\mu$m sources in the field (F$_{870\mu m}\gtrsim$2 mJy), are likely interlopers. The measurement of 870$\mu$m continuum and CO(4-3) line fluxes at the positions of previously-known cluster members provides a deep probe of dusty star formation occurring in the core of this high-redshift structure, adding up to a total SFR~700$\pm$100 M$_{\odot}$/yr and yielding an integrated star formation rate density of ~10$^4$ M$_{\odot}$/yr/Mpc$^3$, five orders of magnitude larger than in the field at the same epoch, due to the concentration of star-forming galaxies in the small volume of the dense cluster core. The combination of these observations with previously available HST imaging highlights the presence in this same volume of a population of galaxies with already suppressed star formation. This diverse composition of galaxy populations in ClJ1449+0856 is especially highlighted at the very cluster center, where a complex assembly of quiescent and star-forming sources is likely forming the future Brightest Cluster Galaxy.
[43]  [pdf] - 1698236
Merger driven star-formation activity in Cl J1449+0856 at z=1.99 as seen by ALMA and JVLA
Comments: 29 pages, 15 figures, 7 tables + Appendix. Accepted for publication in MNRAS
Submitted: 2018-05-24
We use ALMA and JVLA observations of the galaxy cluster Cl J1449+0856 at z=1.99, in order to study how dust-obscured star-formation, ISM content and AGN activity are linked to environment and galaxy interactions during the crucial phase of high-z cluster assembly. We present detections of multiple transitions of $^{12}$CO, as well as dust continuum emission detections from 11 galaxies in the core of Cl J1449+0856. We measure the gas excitation properties, star-formation rates, gas consumption timescales and gas-to-stellar mass ratios for the galaxies. We find evidence for a large fraction of galaxies with highly-excited molecular gas, contributing $>$50% to the total SFR in the cluster core. We compare these results with expectations for field galaxies, and conclude that environmental influences have strongly enhanced the fraction of excited galaxies in this cluster. We find a dearth of molecular gas in the galaxies' gas reservoirs, implying a high star-formation efficiency (SFE) in the cluster core, and find short gas depletion timescales $\tau$<0.1-0.4 Gyrs for all galaxies. Interestingly, we do not see evidence for increased specific star-formation rates (sSFRs) in the cluster galaxies, despite their high SFEs and gas excitations. We find evidence for a large number of mergers in the cluster core, contributing a large fraction of the core's total star-formation compared with expectations in the field. We conclude that the environmental impact on the galaxy excitations is linked to the high rate of galaxy mergers, interactions and active galactic nuclei in the cluster core.
[44]  [pdf] - 1694048
Revolutionizing Our Understanding of AGN Feedback and its Importance to Galaxy Evolution in the Era of the Next Generation Very Large Array
Comments: 26 pages, 9 figures, accepted for publication in ApJ
Submitted: 2018-03-06
Energetic feedback by Active Galactic Nuclei (AGNs) plays an important evolutionary role in the regulation of star formation (SF) on galactic scales. However, the effects of this feedback as a function of redshift and galaxy properties such as mass, environment and cold gas content remain poorly understood. The broad frequency coverage (1 to 116 GHz), high sensitivity (up to ten times higher than the Karl G. Jansky Very Large Array), and superb angular resolution (maximum baselines of at least a few hundred km) of the proposed next generation Very Large Array (ngVLA) are uniquely poised to revolutionize our understanding of AGNs and their role in galaxy evolution. Here, we provide an overview of the science related to AGN feedback that will be possible in the ngVLA era and present new continuum ngVLA imaging simulations of resolved radio jets spanning a wide range of intrinsic extents. We also consider key computational challenges and discuss exciting opportunities for multi-wavelength synergy with other next-generation instruments, such as the Square Kilometer Array and the James Webb Space Telescope. The unique combination of high-resolution, large collecting area, and wide frequency range will enable significant advancements in our understanding of the effects of jet-driven feedback on sub-galactic scales, particularly for sources with extents of a few pc to a few kpc such as young and/or lower-power radio AGNs, AGNs hosted by low-mass galaxies, radio jets that are interacting strongly with the interstellar medium of the host galaxy, and AGNs at high redshift.
[45]  [pdf] - 1712855
Molecular gas in AzTEC/C159: a star-forming disk galaxy 1.3Gyr after the Big Bang
Comments: 10 pages, 4 figures, Accepted for publication in A&A
Submitted: 2017-10-27, last modified: 2018-02-12
We studied the molecular gas properties of AzTEC/C159, a star-forming disk galaxy at $z=4.567$. We secured $^{12}$CO molecular line detections for the $J=2\to1$ and $J=5\to4$ transitions using the Karl G. Jansky VLA and the NOEMA interferometer. The broad (FWHM$\sim750\,{\rm km\,s}^{-1}$) and tentative double-peaked profiles of both $^{12}$CO lines are consistent with an extended molecular gas reservoir, which is distributed in a rotating disk as previously revealed from [CII] 158$\mu$m line observations. Based on the $^{12}$CO(2$\to$1) emission line we derived $L'_{\rm{CO}}=(3.4\pm0.6)\times10^{10}{\rm \,K\,km\,s}^{-1}{\rm \,pc}^{2}$, that yields a molecular gas mass of $M_{\rm H_2 }(\alpha_{\rm CO}/4.3)=(1.5\pm0.3)\times 10^{11}{\rm M}_\odot$ and unveils a gas-rich system with $\mu_{\rm gas}(\alpha_{\rm CO}/4.3)\equiv M_{\rm H_2}/M_\star=3.3\pm0.7$. The extreme star formation efficiency (SFE) of AzTEC/C159, parametrized by the ratio $L_{\rm{IR}}/L'_{\rm{CO}}=(216\pm80)\, {\rm L}_{\odot}{\rm \,(K\,km\,s}^{-1}{\rm \,pc}^{2})^{-1}$, is comparable to merger-driven starbursts such as local ultra-luminous infrared galaxies (ULIRGs) and SMGs. Likewise, the $^{12}$CO(5$\to$4)/CO(2$\to$1) line brightness temperature ratio of $r_{52}= 0.55\pm 0.15$ is consistent with high excitation conditions, similar to that observed in SMGs. We constrained the value for the $L'_{\text{CO}}-{\rm H}_2$ mass conversion factor in AzTEC/C159, i.e. $\alpha_{\text{CO}}=3.9^{+2.7}_{-1.3}{\rm \,M}_{\odot}{\rm \,K}^{-1}{\rm \,km}^{-1}{\rm \,s\,pc}^{-2}$, that is consistent with a self-gravitating molecular gas distribution as observed in local star-forming disk galaxies. Cold gas streams from cosmological filaments might be fueling a gravitationally unstable gas-rich disk in AzTEC/C159, which breaks into giant clumps forming stars as efficiently as in merger-driven systems and generate high gas excitation.
[46]  [pdf] - 1641504
Galaxy pairs in the SDSS - XIII. The connection between enhanced star formation and molecular gas properties in galaxy mergers
Comments: 15 pages, 9 figures
Submitted: 2018-02-06
We investigate the connection between star formation and molecular gas properties in galaxy mergers at low redshift (z$\leq$0.06). The study we present is based on IRAM 30-m CO(1-0) observations of 11 galaxies with a close companion selected from the Sloan Digital Sky Survey (SDSS). The pairs have mass ratios $\leq$4, projected separations r$_{\mathrm{p}} \leq$30 kpc and velocity separations $\Delta$V$\leq$300 km s$^{-1}$, and have been selected to exhibit enhanced specific star formation rates (sSFR). We calculate molecular gas (H$_{2}$) masses, assigning to each galaxy a physically motivated conversion factor $\alpha_{\mathrm{CO}}$, and we derive molecular gas fractions and depletion times. We compare these quantities with those of isolated galaxies from the extended CO Legacy Data base for the GALEX Arecibo SDSS Survey sample (xCOLDGASS, Saintonge et al. 2017) with gas quantities computed in an identical way. Ours is the first study which directly compares the gas properties of galaxy pairs and those of a control sample of normal galaxies with rigorous control procedures and for which SFR and H$_{2}$ masses have been estimated using the same method. We find that the galaxy pairs have shorter depletion times and an average molecular gas fraction enhancement of 0.4 dex compared to the mass matched control sample drawn from xCOLDGASS. However, the gas masses (and fractions) in galaxy pairs and their depletion times are consistent with those of non-mergers whose SFRs are similarly elevated. We conclude that both external interactions and internal processes may lead to molecular gas enhancement and decreased depletion times.
[47]  [pdf] - 1629243
"Super-deblended" Dust Emission in Galaxies: I. The GOODS-North Catalog and the Cosmic Star Formation Rate Density out to Redshift 6
Comments: Published in The Astrophysical Journal, Volume 853, Number 2. The entire machine-readable catalog is available from ApJ at . We also provide the same catalog in FITS format at
Submitted: 2017-03-15, last modified: 2018-02-05
We present a new technique to measure multi-wavelength "Super-deblended" photometry from highly confused images, which we apply to Herschel and ground-based far-infrared (FIR) and (sub-)millimeter (mm) data in the northern field of the Great Observatories Origins Deep Survey (GOODS). There are two key novelties. First, starting with a large database of deep Spitzer 24{\mu}m and VLA 20cm detections that are used to define prior positions for fitting the FIR/submm data, we perform an active selection of useful priors independently at each frequency band, moving from less to more confused bands. Exploiting knowledge of redshift and all available photometry, we identify hopelessly faint priors that we remove from the fitting pool. This approach significantly reduces blending degeneracies and allows reliable photometry to be obtained for galaxies in FIR+mm bands. Second, we obtain well-behaved, nearly Gaussian flux density uncertainties, individually tailored to all fitted priors in each band. This is done by exploiting extensive simulations that allow us to calibrate the conversion of formal fitting uncertainties to realistic uncertainties depending on quantities directly measurable. We achieve deeper detection limits with high fidelity measurements and uncertainties at FIR+mm bands. As an illustration of the utility of these measurements, we identify 70 galaxies with z>3 and reliable FIR+mm detections. We present new constraints on the cosmic star formation rate density at 3<z<6, finding a significant contribution from z>3 dusty galaxies that are missed by optical-to-near-infrared color selection. Photometric measurements for 3306 priors, including over 1000 FIR+mm detections are released publicly with our catalog.
[48]  [pdf] - 1674796
Molecular outflow and feedback in the obscured Quasar XID2028 revealed by ALMA
Comments: 15 pages, 11 figures (including Appendix). Version matches A&A published version (after proof)
Submitted: 2017-12-12, last modified: 2018-01-31
We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.
[49]  [pdf] - 1659667
Red, redder, reddest: SCUBA-2 imaging of colour-selected \textit{Herschel} sources
Comments: 22 pages, 14 figures and 3 tables
Submitted: 2018-01-22
High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).
[50]  [pdf] - 1617633
The unexpectedly large dust and gas content of quiescent galaxies at z>1.4
Comments: 16 pages, 4 main figures main, 7 supplementary figures, 2 supplementary tables; accepted for publication in Nature Astronomy
Submitted: 2017-03-06, last modified: 2018-01-15
Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time.
[51]  [pdf] - 1626352
The infrared-radio correlation of spheroid- and disc-dominated star-forming galaxies to z $\sim$ 1.5 in the COSMOS field
Comments: Resubmitted to MNRAS after addressing reviewer feedback
Submitted: 2017-10-20
Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.
[52]  [pdf] - 1587378
Radio selection of the most distant galaxy clusters
Comments: ApJ Letters in press
Submitted: 2017-08-23
We show that the most distant X-ray detected cluster known to date, ClJ1001 at z=2.506, hosts a strong overdensity of radio sources. Six of them are individually detected (within 10") in deep 0.75" resolution VLA 3GHz imaging, with S(3GHz)>8uJy. Of the six, AGN likely affects the radio emission in two galaxies while star formation is the dominant source powering the remaining four. We searched for cluster candidates over the full COSMOS 2-square degree field using radio-detected 3GHz sources and looking for peaks in Sigma5 density maps. ClJ1001 is the strongest overdensity by far with >10sigma, with a simple z_phot>1.5 preselection. A cruder photometric rejection of z<1 radio foregrounds leaves ClJ1001 as the second strongest overdensity, while even using all radio sources ClJ1001 remains among the four strongest projected overdensities. We conclude that there are great prospects for future, deep and wide-area radio surveys to discover large samples of the first generation of forming galaxy clusters. In these remarkable structures widespread star formation and AGN activity of massive galaxy cluster members, residing within the inner cluster core, will ultimately lead to radio continuum as one of the most effective means for their identification, with detection rates expected in the ballpark of 0.1-1 per square degree at z>2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models.
[53]  [pdf] - 1598002
The implications of the surprising existence of a large, massive CO disk in a distant protocluster
Comments: Accepted by A&A. 15 pages, 8 figures and 4 tables
Submitted: 2017-01-18, last modified: 2017-07-09
It is not yet known if the properties of molecular gas in distant protocluster galaxies are significantly affected by their environment as galaxies are in local clusters. Through a deep, 64 hours of effective on-source integration with the Australian Telescope Compact Array (ATCA), we discovered a massive, M_mol=2.0+-0.2x10^11 M_sun, extended, ~40 kpc, CO(1-0)-emitting disk in the protocluster surrounding the radio galaxy, MRC1138-262. The galaxy, at z_CO=2.1478, is a clumpy, massive disk galaxy, M_star~5x10^11 M_sun, which lies 250 kpc in projection from MRC1138-262 and is a known H-alpha emitter, HAE229. HAE229 has a molecular gas fraction of ~30%. The CO emission has a kinematic gradient along its major axis, centered on the highest surface brightness rest-frame optical emission, consistent with HAE229 being a rotating disk. Surprisingly, a significant fraction of the CO emission lies outside of the UV/optical emission. In spite of this, HAE229 follows the same relation between star-formation rate and molecular gas mass as normal field galaxies. HAE229 is the first CO(1-0) detection of an ordinary, star-forming galaxy in a protocluster. We compare a sample of cluster members at z>0.4 that are detected in low-order CO transitions with a similar sample of sources drawn from the field. We confirm findings that the CO-luminosity and FWHM are correlated in starbursts and show that this relation is valid for normal high-z galaxies as well as those in overdensities. We do not find a clear dichotomy in the integrated Schmidt-Kennicutt relation for protocluster and field galaxies. Not finding any environmental dependence in the "star-formation efficiency" or the molecular gas content, especially for such an extended CO disk, suggests that environmentally-specific processes such as ram pressure stripping are not operating efficiently in (proto)clusters. (abridged)
[54]  [pdf] - 1581802
ALMA constraints on star-forming gas in a prototypical z=1.5 clumpy galaxy: the dearth of CO(5-4) emission from UV-bright clumps
Comments: Accepted for publication in MNRAS. Minor revisions to match published version
Submitted: 2017-03-07, last modified: 2017-06-08
We present deep ALMA CO(5-4) observations of a main sequence, clumpy galaxy at z=1.5 in the HUDF. Thanks to the ~0.5" resolution of the ALMA data, we can link stellar population properties to the CO(5-4) emission on scales of a few kpc. We detect strong CO(5-4) emission from the nuclear region of the galaxy, consistent with the observed $L_{\rm IR}$-$L^{\prime}_{\rm CO(5-4)}$ correlation and indicating on-going nuclear star formation. The CO(5-4) gas component appears more concentrated than other star formation tracers or the dust distribution in this galaxy. We discuss possible implications of this difference in terms of star formation efficiency and mass build-up at the galaxy centre. Conversely, we do not detect any CO(5-4) emission from the UV-bright clumps. This might imply that clumps have a high star formation efficiency (although they do not display unusually high specific star formation rates) and are not entirely gas dominated, with gas fractions no larger than that of their host galaxy (~50%). Stellar feedback and disk instability torques funnelling gas towards the galaxy centre could contribute to the relatively low gas content. Alternatively, clumps could fall in a more standard star formation efficiency regime if their actual star-formation rates are lower than generally assumed. We find that clump star-formation rates derived with several different, plausible methods can vary by up to an order of magnitude. The lowest estimates would be compatible with a CO(5-4) non-detection even for main-sequence like values of star formation efficiency and gas content.
[55]  [pdf] - 1583087
The Dust Attenuation of Star-forming Galaxies at $z\sim3$ and Beyond: New Insights from ALMA Observations
Comments: 8 pages, 4 figures, submitted to MNRAS; v2 - Small changes to text
Submitted: 2017-05-03, last modified: 2017-05-25
We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.
[56]  [pdf] - 1583552
Dust and Gas in Star Forming Galaxies at z~3 - Extending Galaxy Uniformity to 11.5 Billion Years
Comments: Accepted for publication at A&A
Submitted: 2017-05-17
We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.
[57]  [pdf] - 1582197
VLA-COSMOS 3GHz Large Project: The infrared-radio correlation of star-forming galaxies and AGN to $z\lesssim6$
Comments: 17 pages, 22 figures, 3 tables. Accepted for publication in A&A
Submitted: 2017-03-28, last modified: 2017-05-04
We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.
[58]  [pdf] - 1582198
The VLA-COSMOS 3 GHz Large Project: The cosmic star formation history since z~5
Comments: 18 pages, 8 figures, 2 tables, accepted for publication in A&A
Submitted: 2017-03-28, last modified: 2017-03-30
We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.
[59]  [pdf] - 1582190
The VLA-COSMOS 3 GHz Large Project: Continuum data and source catalog release
Comments: 19 pages, 18 figures, 5 tables, accepted for publication in A&A
Submitted: 2017-03-28
We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.
[60]  [pdf] - 1552495
ALMA observations of cold molecular gas in AGN hosts at z~1.5 - Evidence of AGN feedback?
Comments: 12 pages, 5 figures, 3 tables, accepted for publication in MNRAS
Submitted: 2017-03-23
Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 1<z<3. This cosmic epoch is hence best suited for investigating the effects of radiative feedback from AGN. Observational efforts are underway to quantify the impact of AGN feedback, if any, on their host galaxies. Here we present a study of the molecular gas content of AGN hosts at z~1.5 using CO[2-1] line emission observed with ALMA for a sample of 10 AGNs. We compare this with a sample of galaxies without an AGN matched in redshift, stellar mass, and star formation rate. We detect CO in 3 AGNs with $\mathrm{L_{CO} \sim 6.3-25.1\times 10^{9} L_{\odot}}$ which translates to a molecular hydrogen gas mass of $\mathrm{2.5-10\times 10^{10} M_{\odot}}$ assuming conventional conversion factor of $\mathrm{\alpha_{CO}}\sim3.6$. Our results indicate a >99% probability of lower depletion time scales and lower molecular gas fractions in AGN hosts with respect to the non-AGN comparison sample. We discuss the implications of these observations on the impact that AGN feedback may have on star formation efficiency of z>1 galaxies.
[61]  [pdf] - 1548688
Further evidence for a quasar-driven jet impacting its neighbour galaxy: The saga of HE0450-2958 continues
Comments: 11 pages including 5 figures, published in MNRAS
Submitted: 2017-01-11, last modified: 2017-03-17
HE0450-2958, an interacting quasar-starburst galaxy pair at $z$ = 0.285, is one of the best known examples of strong star formation activity in the presence of a quasar-driven jet. We present new multi-band JVLA-imaging covering 1 to 6 GHz and reaching an angular resolution of up to $0{}_{{}^.}^{"}6$ (a 6-fold improvement over existing radio data). We confirm the previous detection of a spatially extended radio component around the quasar indicating that there is on-going star formation activity in the quasar host galaxy. For the first time, we directly detect a jet-like bipolar outflow from the quasar aligned with its companion star-forming galaxy (SFG) and several blobs of ionized gas in its vicinity identified in previous studies. Within the companion SFG we find evidence for a flattening of the synchrotron spectral index towards the point of intersection with the jet axis, further suggesting that the outflow may actually be impacting its interstellar medium (ISM). We discuss two possible mechanisms that could have triggered the starburst in the companion SFG: a wet-dry merger with the quasar and jet-induced star formation. While triggering through interaction-driven gas dynamics cannot be excluded with current data, our new observations make HE0450-2958 a strong candidate for jet-induced star formation, and one of the rare links between local systems (like Minkowski's Object or Centaurus A) and the high-z regime where radio-optical alignments suggest that this phenomenon could be more common.
[62]  [pdf] - 1542861
In and out star formation in z~1.5 quiescent galaxies from rest-frame UV spectroscopy and the far-infrared
Comments: 12 pages, 10 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2016-12-20
We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.
[63]  [pdf] - 1532017
Gas Fraction and Depletion Time of Massive Star Forming Galaxies at z~3.2: No Change in Global Star Formation Process out to z>3
Comments: 24 pages plus 16 pages references and 5 figures, ApJ, accepted
Submitted: 2016-10-12
The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$\mu$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm \mu_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1\sigma$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm \mu_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.
[64]  [pdf] - 1528296
The evolution of the star formation rate function and cosmic star formation rate density of galaxies at $z \sim 1-4$
Comments: 20 pages, 7 figures, 10 tables, accepted for publication in MNRAS
Submitted: 2016-10-11
We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 1-4 $ galaxies, using cosmological Smoothed Particle Hydrodynamic (SPH) simulations and a compilation of UV, IR and H$\alpha$ observations. These tracers represent different populations of galaxies with the IR light being a probe of objects with high star formation rates and dust contents, while UV and H$\alpha$ observations provide a census of low star formation galaxies where mild obscuration occurs. We compare the above SFRFs with the results of SPH simulations run with the code {\small{P-GADGET3(XXL)}}. We focus on the role of feedback from Active Galactic Nuclei (AGN) and supernovae in form of galactic winds. The AGN feedback prescription that we use decreases the simulated CSFRD at $z < 3$ but is not sufficient to reproduce the observed evolution at higher redshifts. We explore different wind models and find that the key factor for reproducing the evolution of the observed SFRF and CSFRD at $z \sim1-4$ is the presence of a feedback prescription that is prominent at high redshifts ($z \ge 4$) and becomes less efficient with time. We show that variable galactic winds which are efficient at decreasing the SFRs of low mass objects are quite successful in reproducing the observables.
[65]  [pdf] - 1531094
ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas
Comments: Accepted for publication in the Astrophysical Journal
Submitted: 2016-07-22, last modified: 2016-09-12
In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).
[66]  [pdf] - 1531095
ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies
Comments: Accepted for publication in the Astrophysical Journal
Submitted: 2016-07-22, last modified: 2016-09-12
We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.
[67]  [pdf] - 1475316
Discovery of a galaxy cluster with a violently starbursting core at z=2.506
Comments: matched to the published version by ApJ
Submitted: 2016-04-25, last modified: 2016-08-30
We report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$\sigma$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $H\alpha$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation.
[68]  [pdf] - 1531092
ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description
Comments: Submitted to the Astrophysical Journal, incorporating the first round of referee comments
Submitted: 2016-07-22
We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.
[69]  [pdf] - 1486852
A giant Ly$\alpha$ nebula in the core of an X-ray cluster at $z=1.99$: implications for early energy injection
Comments: 15 pages, 11 figures, 1 appendix. Accepted for publication in ApJ
Submitted: 2016-05-10, last modified: 2016-07-04
We present the discovery of a giant $\gtrsim$100~kpc Ly$\alpha$ nebula detected in the core of the X-ray emitting cluster CL~J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$\alpha$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Ly$\alpha$ to X-ray luminosity ratio ($L_{\mathrm{Ly\alpha}}/L_{\mathrm{X}} \approx0.3$, $\gtrsim10-1000\times$ higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Ly$\alpha$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$\alpha$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ M$_{\odot}$ yr$^{-1}$. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. In our baseline calculation AGN provide up to $85$% of the injected energy and 2/3 of the mass, while the rest is supplied by supernovae-driven winds.
[70]  [pdf] - 1486886
HELP: XID+, The Probabilistic De-blender for Herschel SPIRE maps
Comments: 14 pages, 10 Figures
Submitted: 2016-06-18
We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates. In this paper, we discuss the details of XID+ and demonstrate the basic capabilities and performance by running it on simulated SPIRE maps resembling the COSMOS field, and comparing to the current prior-based source extraction tool DESPHOT. We show that not only does XID+ perform better on metrics such as flux accuracy and flux uncertainty accuracy, we illustrate how obtaining the posterior probability distribution can help overcome some of the issues inherent with maximum likelihood based source extraction routines. We run XID+ on the COSMOS SPIRE maps from HerMES, using a 24 $\mathrm{\mu m}$ catalogue as a prior and show the marginalised SPIRE colour-colour plot and marginalised contribution to the cosmic infrared background at the SPIRE wavelengths. XID+ is a core tool arising from the Herschel Extragalactic Legacy Project (HELP) and we discuss how additional work within HELP providing prior information on fluxes can and will be utilised. The software is available at We also provide the data product for COSMOS. We believe this is the first time that the full posterior probability of galaxy photometry has been provided as a data product.
[71]  [pdf] - 1475375
HELP: star formation as function of galaxy environmentwith Herschel
Comments: 14 pages, 7 figures. MNRAS accepted
Submitted: 2016-06-16
The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$.
[72]  [pdf] - 1411549
The multiplicity of 250-$\mu$m Herschel sources in the COSMOS field
Comments: 12 pages, 7 figures. MNRAS accepted
Submitted: 2016-04-29
We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$\mu$m catalogues, in conjunction with 250-$\mu$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $\mu$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$\mu$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$\mu$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$\mu$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$\mu$m source flux. At the highest 250-$\mu$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$\mu$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent.
[73]  [pdf] - 1374241
The Chandra COSMOS Legacy survey: overview and point source catalog
Comments: 20 pages, 16 figures, accepted for publication on ApJ on December 24, 2015
Submitted: 2016-01-05
The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.
[74]  [pdf] - 1298118
A higher efficiency of converting gas to stars push galaxies at z ~ 1.6 well above the star-forming main sequence
Comments: Published in ApJ Letters; 6 pages, 4 figures
Submitted: 2015-05-19, last modified: 2015-10-22
Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the molecular gas content of seven high-redshift (z ~ 1.6) starburst galaxies with the Atacama Large (sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the sample of Herschel far-infrared detected galaxies having star formation rates (~300-800 Msolar/yr) elevated (>4x) above the star-forming `main sequence', and included in the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect CO emission in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more compelling, we firmly establish with a clean and systematic selection that starbursts, identified as main-sequence outliers, at high redshift generally have a lower ratio of CO to total infrared luminosity as compared to typical main-sequence star-forming galaxies, although with a smaller offset than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. Along with a heightened star formation efficiency, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the interstellar medium.
[75]  [pdf] - 1298204
Next Generation Very Large Array Memo No. 8 Science Working Group 3: Galaxy Assembly through Cosmic Time
Comments: 52 pages, 13 figures, NRAO Next Generation Very Large Array Memos Series:
Submitted: 2015-10-21
The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.
[76]  [pdf] - 1308351
Quantifying the UV continuum slopes of galaxies to z~10 using deep Hubble and Spitzer/IRAC observations
Comments: 10 pages, 10 figures, accepted in MNRAS
Submitted: 2015-10-06
Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes.
[77]  [pdf] - 1449859
(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - II. The spatial extent of the radio-emitting regions
Comments: 18 pages (incl. 5 appendices), 7 figures, and 5 tables; accepted for publication in A&A; abstract abridged for arXiv
Submitted: 2015-09-23
Radio emission at cm wavelengths from highly star-forming galaxies, such as SMGs, is dominated by synchrotron radiation arising from supernova activity. Using deep, high-resolution ($1\sigma=2.3$ $\mu$Jy beam$^{-1}$; $0.75^{"}$) cm radio-continuum observations taken by the VLA-COSMOS 3 GHz Large Project, we studied the radio-emitting sizes of a flux-limited sample of SMGs in the COSMOS field. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them ($\sim46\pm11\%$) with S/N ratios in the range of ${\rm S/N=4.2-37.4}$. Using 2D elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of $0.54^{"}\pm 0.11^{"}$ for our 18 SMGs detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of $4.2\pm0.9$ kpc. No clear correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, there is a hint of larger radio sizes at $z\sim2.5-5$ compared to lower redshifts. The sizes we derived are consistent with previous SMG sizes measured at 1.4 GHz and in mid-$J$ CO emission, but significantly larger than those seen in the (sub)mm continuum emission. One possible scenario is that SMGs have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-$J$ CO line emission and radio continuum emission, and ii) a warmer, compact starburst region giving rise to the high-$J$ line emission of CO, which could dominate the dust continuum size measurements. Because of the rapid cooling of CR electrons in dense starburst galaxies ($\sim10^4-10^5$ yr), the more extended synchrotron radio-emitting size being a result of CR diffusion seems unlikely. Instead, if SMGs are driven by galaxy mergers the radio synchrotron emission might arise from more extended magnetised ISM around the starburst region.
[78]  [pdf] - 1273162
ALMA and Herschel reveal that X-ray selected AGN and main-sequence galaxies have different star formation rate distributions
Comments: 5 pages, 3 figures, 1 table; Accepted for publication in MNRAS Letters after moderate changes
Submitted: 2015-06-17, last modified: 2015-07-29
Using deep Herschel and ALMA observations, we investigate the star formation rate (SFR) distributions of X-ray selected AGN host galaxies at 0.5<z<1.5 and 1.5<z<4, comparing them to that of normal, star-forming (i.e., "main-sequence", or MS) galaxies. We find 34--55 per cent of AGNs in our sample have SFRs at least a factor of two below that of the average MS galaxy, compared to ~15 per cent of all MS galaxies, suggesting significantly different SFR distributions. Indeed, when both are modelled as log-normal distributions, the mass and redshift-normalised SFR distributions of X-ray AGNs are roughly twice as broad, and peak ~0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the normalised SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and thus with previous results derived from stacked Herschel data. This apparent contradiction is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily represent a true characterisation of the typical SFR of X-ray AGNs.
[79]  [pdf] - 1269386
Satellite content and quenching of star formation in galaxy groups at z~1.8
Comments: 12 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics
Submitted: 2015-07-14
We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.
[80]  [pdf] - 1182859
A direct constraint on the gas content of a massive, passively evolving elliptical galaxy at z = 1.43
Comments: 6 pages, 3 figures & 1 table. Accepted for publication in ApJL
Submitted: 2015-01-30, last modified: 2015-05-29
Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.
[81]  [pdf] - 1245616
GOODS-HERSCHEL: star formation, dust attenuation and the FIR-radio correlation on the Main Sequence of star-forming galaxies up to z~4
Comments: 19 pages, 12 figure, accepted for publication in ApJ
Submitted: 2014-07-18, last modified: 2015-05-18
We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.
[82]  [pdf] - 1055817
Evidence for feedback in action from the molecular gas content in the z~1.6 outflowing QSO XID2028
Comments: 8 pages, 4 figures. Astronomy and Astrophysics, in press
Submitted: 2015-03-05
Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host
[83]  [pdf] - 979066
The influence of wavelength, flux, and lensing selection effects on the redshift distribution of dusty, star-forming galaxies
Comments: 4 pages, 1 table, 4 figures, accepted by A&A
Submitted: 2015-02-16
We interpret the large variety of redshift distributions of galaxies found by far-infrared and (sub-)millimeter deep surveys depending on their depth and wavelength using the B\'ethermin et al. (2012) phenomenological model of galaxy evolution. This model reproduces without any new parameter tuning the observed redshift distributions from 100 $\mu$m to 1.4 mm, and especially the increase of the median redshift with survey wavelength. This median redshift varies also significantly with the depth of the surveys, and deeper surveys do necessarily not probe higher redshifts. Paradoxically, at fixed wavelength and flux limit, the lensed sources are not always at higher redshift. We found that the higher redshift of 1.4 mm-selected south pole telescope (SPT) sources compared to other SMG surveys is not only caused by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a population of lensed main-sequence galaxies and a minor contribution ($\sim$10\%) of unlensed extreme starbursts.
[84]  [pdf] - 1043010
CO excitation of normal star forming galaxies out to z=1.5 as regulated by the properties of their interstellar medium
Comments: Accepted for publication in A&A
Submitted: 2014-09-29, last modified: 2015-01-18
We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.
[85]  [pdf] - 982096
Physical properties of z>4 submillimeter galaxies in the COSMOS field
Comments: 16 pages, 8 figures, 6 tables; to appear in A&A
Submitted: 2014-12-11, last modified: 2015-01-15
We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.
[86]  [pdf] - 942314
The Herschel view of the dominant mode of galaxy growth from z=4 to the present day
Comments: 30 pages, 26 figures, accepted for publication in A&A
Submitted: 2014-09-18, last modified: 2015-01-07
We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.
[87]  [pdf] - 915527
The Astrophysics of Star Formation Across Cosmic Time at $\gtrsim$10 GHz with the Square Kilometre Array
Comments: 15 pages, 3 figures, to appear as part of 'Continuum Science' in Proceedings "Advancing Astrophysics with the Square Kilometre Array" PoS(AASKA14)085; "Kilometer" replaced by "Kilometre" in updated version
Submitted: 2014-12-17, last modified: 2014-12-30
In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz.
[88]  [pdf] - 910664
Astronomy below the Survey Threshold
Comments: 13 pages, 5 figures, to appear as part of 'Continuum Science' in Proceedings 'Advancing Astrophysics with the SKA (AASKA14)' [PoS(AASKA14)172]
Submitted: 2014-12-18
Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.
[89]  [pdf] - 910674
Enabling the next generation of cm-wavelength studies of high-redshift molecular gas with the SKA
Comments: 10 pages, to be published in Proceedings of Science: "Advancing Astrophysics with the Square Kilometre Array", PoS(AASKA14)161
Submitted: 2014-12-18
The Square Kilometre Array will be a revolutionary instrument for the study of gas in the distant Universe. SKA1 will have sufficient sensitivity to detect and image atomic 21 cm HI in individual galaxies at significant cosmological distances, complementing ongoing ALMA imaging of redshifted high-J CO line emission and far-infrared interstellar medium lines such as [CII] 157.7 um. At frequencies below ~50 GHz, observations of redshifted emission from low-J transitions of CO, HCN, HCO+, HNC, H2O and CS provide insight into the kinematics and mass budget of the cold, dense star-forming gas in galaxies. In advance of ALMA band 1 deployment (35 to 52 GHz), the most sensitive facility for high-redshift studies of molecular gas operating below 50~GHz is the Karl G. Jansky Very Large Array (VLA). Here, we present an overview of the role that the SKA could play in molecular emission line studies during SKA1 and SKA2, with an emphasis on studies of the dense gas tracers directly probing regions of active star-formation.
[90]  [pdf] - 1223457
COLDz: Karl G. Jansky Very Large Array discovery of a gas-rich galaxy in COSMOS
Submitted: 2014-12-03
The broad spectral bandwidth at mm and cm-wavelengths provided by the recent upgrades to the Karl G. Jansky Very Large Array (VLA) has made it possible to conduct unbiased searches for molecular CO line emission at redshifts, z > 1.31. We present the discovery of a gas-rich, star-forming galaxy at z = 2.48, through the detection of CO(1-0) line emission in the COLDz survey, through a sensitive, Ka-band (31 to 39 GHz) VLA survey of a 6.5 square arcminute region of the COSMOS field. We argue that the broad line (FWHM ~570 +/- 80 km/s) is most likely to be CO(1-0) at z=2.48, as the integrated emission is spatially coincident with an infrared-detected galaxy with a photometric redshift estimate of z = 3.2 +/- 0.4. The CO(1-0) line luminosity is L'_CO = (2.2 +/- 0.3) x 10^{10} K km/s pc^2, suggesting a cold molecular gas mass of M_gas ~ (2 - 8)x10^{10}M_solar depending on the assumed value of the molecular gas mass to CO luminosity ratio alpha_CO. The estimated infrared luminosity from the (rest-frame) far-infrared spectral energy distribution (SED) is L_IR = 2.5x10^{12} L_solar and the star-formation rate is ~250 M_solar/yr, with the SED shape indicating substantial dust obscuration of the stellar light. The infrared to CO line luminosity ratio is ~114+/-19 L_solar/(K km/s pc^2), similar to galaxies with similar SFRs selected at UV/optical to radio wavelengths. This discovery confirms the potential for molecular emission line surveys as a route to study populations of gas-rich galaxies in the future.
[91]  [pdf] - 920390
Evolution of the dust emission of massive galaxies up to z=4 and constraints on their dominant mode of star formation
Comments: 17 pages, 14 figures, 2 tables, accepted by A&A
Submitted: 2014-09-19, last modified: 2014-11-19
We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.
[92]  [pdf] - 1165348
Regularity underlying complexity: a redshift-independent description of the continuous variation of galaxy-scale molecular gas properties in the mass-star formation rate plane
Comments: 37 pages, 18 figures. Published in ApJ, 793, 19
Submitted: 2013-03-18, last modified: 2014-09-11
Star-forming galaxies (SFGs) display a continuous distribution of specific star formation rates (sSFR) which can be approximated by the superposition of two log-normal distributions. The 1st of these encompasses the main sequence (MS) of SFGs, the 2nd one a rarer population of starbursts (SB). We show that the sSFR-distribution of SBs can be regarded as the result of a physical process (plausibly merging) taking the mathematical form of a log-normal boosting kernel and enhancing star formation activity. We explore the utility of splitting the star-forming population into MS and SB galaxies - an approach we term "2-Star Formation Mode" (2-SFM) framework - for understanding their molecular gas properties. Variations of star formation efficiency (SFE) and gas fractions among SFGs take a simple, redshift-independent form, once these quantities are normalized to the value of an average MS galaxy. The change in SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected for mergers. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane that separates more clearly into loci for SBs and normal galaxies than observed in SFR vs. M* space. SBs with the largest deviations (>10-fold) from the MS, like many local ULIRGs, are not average SBs, but even rarer events with progenitors having larger gas fractions than typical MS galaxies. We statistically infer that the gas fractions of typical SBs decrease by a factor of 2 to 3 with respect to their direct MS progenitors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 conversion factor in the SFR-M* plane and provide evidence that the higher sSFR of distant galaxies is a direct consequence of larger gas fraction in these systems.
[93]  [pdf] - 873817
Dust and gas in luminous proto-cluster galaxies at z=4.05: the case for different cosmic dust evolution in normal and starburst galaxies
Comments: 17 pages, 9 figures. Updated to A&A accepted version
Submitted: 2014-03-31, last modified: 2014-08-15
We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.
[94]  [pdf] - 861922
Inter-comparison of Radio-Loudness Criteria for Type 1 AGNs in the XMM-COSMOS Survey
Comments: 12 pages, 7 figures, MNRAS submitted
Submitted: 2014-08-05
Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24\mu m}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500\AA})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(\nu L_{\nu}(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24\mu m}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.
[95]  [pdf] - 843100
CO(1-0) line imaging of massive star-forming disc galaxies at z=1.5-2.2
Comments: Accepted 26 April 2014; Title modified to match published version
Submitted: 2014-04-29, last modified: 2014-06-26
We present detections of the CO(J= 1-0) emission line in a sample of four massive star-forming galaxies at z~1.5-2.2 obtained with the Karl G. Jansky Very Large Array (VLA). Combining these observations with previous CO(2-1) and CO(3-2) detections of these galaxies, we study the excitation properties of the molecular gas in our sample sources. We find an average line brightness temperature ratios of R_{21}=0.70+\-0.16 and R_{31}=0.50+\-0.29, based on measurements for three and two galaxies, respectively. These results provide additional support to previous indications of sub-thermal gas excitation for the CO(3-2) line with a typically assumed line ratio R_{31}~0.5. For one of our targets, BzK-21000, we present spatially resolved CO line maps. At the resolution of 0.18 arcsec (1.5 kpc), most of the emission is resolved out except for some clumpy structure. From this, we attempt to identify molecular gas clumps in the data cube, finding 4 possible candidates. We estimate that <40 % of the molecular gas is confined to giant clumps (~1.5 kpc in size), and thus most of the gas could be distributed in small fainter clouds or in fairly diffuse extended regions of lower brightness temperatures than our sensitivity limit.
[96]  [pdf] - 1209389
Modeling the connection between ultraviolet and infrared galaxy populations across cosmic times
Comments: Accepted for publication in MNRAS; 14 pages; 7 figures
Submitted: 2014-05-06
Using a phenomenological approach, we self-consistently model the redshift evolution of the ultraviolet (UV) and infrared (IR) luminosity functions across cosmic time, as well as a range of observed IR properties of UV-selected galaxy population. This model is an extension of the 2SFM (2 star-formation modes) formalism, which is based on the observed "main-sequence" of star-forming galaxies, i.e. a strong correlation between their stellar mass and their star formation rate (SFR), and a secondary population of starbursts with an excess of star formation. The balance between the UV light from young, massive stars and the dust-reprocessed IR emission is modeled following the empirical relation between the attenuation (IRX for IR excess hereafter) and the stellar mass, assuming a scatter of 0.4\,dex around this relation. We obtain a good overall agreement with the measurements of the IR luminosity function up to z~3 and the UV luminosity functions up to z~6, and show that a scatter on the IRX-M relation is mandatory to reproduce these observables. We also naturally reproduce the observed, flat relation between the mean IRX and the UV luminosity at L$_{\rm UV}>$10$^{9.5}$ L$_\odot$. Finally, we perform predictions of the UV properties and detectability of IR-selected samples and the vice versa, and discuss the results in the context of the UV-rest-frame and sub-millimeter surveys of the next decade.
[97]  [pdf] - 1208158
Gas reservoir of a hyper-luminous QSO at z=2.6
Comments: Submitted to A&A
Submitted: 2014-03-04
Understanding the relationship between the formation and evolution of galaxies and their central super massive black holes (SMBH) is one of the main topics in extragalactic astrophysics. Links and feedback may reciprocally affect both black hole and galaxy growth. Observations of the CO line at redshifts of 2-4 are crucial to investigate the gas mass, star formation activity and accretion onto SMBHs, as well as the effect of AGN feedback. Potential correlations between AGN and host galaxy properties can be highlighted by observing extreme objects. Despite their luminosity, hyper-luminous QSOs at z=2-4 are still little studied at mm wavelengths. We targeted CO(3-2) in ULAS J1539+0557, an hyper-luminos QSO (Lbol> 10^48 erg/s) at z=2.658, selected through its unusual red colors in the UKIDSS Large Area Survey (ULAS). We find a molecular gas mass of 4.1+-0.8 10^10 Msun, and a gas fraction of 0.4-0.1, depending mostly on the assumed source inclination. We also find a robust lower limit to the star-formation rate (SFR=250-1600 Msun/yr) and star-formation efficiency (SFE=25-350 Lsun/(K km s-1 pc2) by comparing the observed optical-near-infrared spectral energy distribution with AGN and galaxy templates. The black hole gas consumption timescale, M(H_2)/dM(accretion)/dt, is ~160 Myr, similar or higher than the gas consumption timescale. The gas content and the star formation efficiency are similar to those of other high-luminosity, highly obscured QSOs, and at the lower end of the star-formation efficiency of unobscured QSOs, in line with predictions from AGN-galaxy co-evolutionary scenarios. Further measurements of the (sub)-mm continuum in this and similar sources are mandatory to obtain a robust observational picture of the AGN evolutionary sequence.
[98]  [pdf] - 1202426
A Molecular Line Scan in the Hubble Deep Field North: Constraints on the CO Luminosity Function and the Cosmic H2 Density
Comments: Accepted for publication in ApJ
Submitted: 2013-12-22
We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $\rho_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.01<z<1.89 and z>2. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $\rho_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $\rho_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $\rho_{\rm H2}$(z) in the redshift bins <z>~1.5 and <z>~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $\rho_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory.
[99]  [pdf] - 1202425
A molecular line scan in the Hubble Deep Field North
Comments: Accepted for publication in ApJ
Submitted: 2013-12-22
We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 1<z<1.9 and all z>2. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $\mu$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1<z<3, unless high dust temperatures are invoked. The present study represents a first, fundamental step towards an unbiased census of molecular gas in `normal' galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.
[100]  [pdf] - 1179400
A deep search for molecular gas in two massive Lyman break galaxies at z=3 and 4: vanishing CO-emission due to low metallicity?
Comments: 9 pages, 3 figures, accepted for publication in Astrophysical Journal Letters
Submitted: 2013-09-21
We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA.
[101]  [pdf] - 1164763
The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4
Comments: 32 pages, 20 figures, 9 tables. Published in MNRAS. Replaced Fig. 8 (small scaling bug in the previous version)
Submitted: 2013-02-21, last modified: 2013-07-25
We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.
[102]  [pdf] - 1165958
The redshift-evolution of the distribution of star formation among dark matter halos as seen in the infrared
Comments: 25 pages, 18 figures, 3 tables, revised version accepted by A&A
Submitted: 2013-04-14, last modified: 2013-07-09
[Abridged] Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of star-forming galaxies, the so-called star-forming main sequence. An empirical modeling approach (2-SFM) which distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT and SPT. We use this model to infer the redshift distribution these observables, as well as the level of correlation between CIB-anisotropies at different wavelengths. We also predict that more than 90% of cosmic star formation activity occurs in halos with masses between 10^11.5 and 10^13.5 Msun. Taking into account subsequent mass growth of halos, this implies that the majority of stars were initially (at z>3) formed in the progenitors of clusters, then in groups at 0.5<z<3 and finally in Milky-Way-like halos at z<0.5. At all redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~10^12 Msun, in which the instantaneous star formation efficiency is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material available at
[103]  [pdf] - 1159488
On the effect of the Cosmic Microwave Background in high-redshift (sub-)millimeter observations
Comments: 13 pages, 11 figures. Accepted for publication in the Astrophysical Journal
Submitted: 2013-02-04
Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.
[104]  [pdf] - 1159042
Empirical predictions for (sub-)millimeter line and continuum deep fields
Comments: 25 pages, 14 figures, 7 tables. Accepted for publication in the Astrophysical Journal. Predictions in Tables 6 and 7 are available at:
Submitted: 2013-01-14
[abridged] Modern (sub-)millimeter/radio interferometers will enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will provide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission. In this paper, we present empirical predictions for such (sub-)millimeter line and continuum deep fields. We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble Ultra Deep Field. Using a physically-motivated spectral energy distribution model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z=5, and obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminosity. We derive statistical constraints on the dust emission in the infrared and (sub-)millimeter which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA and JVLA bands. Using empirical relations between the observed CO/[CII] line luminosities and the infrared luminosity, we infer the flux of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our sample. We then predict the fluxes of higher CO transition lines CO(2-1) to CO(7-6) bracketing two extreme gas excitation scenarios. We use our predictions to discuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct benchmark for future deep field campaigns in the (sub-)millimeter regime.
[105]  [pdf] - 1151963
GOODS-Herschel: radio-excess signature of hidden AGN activity in distant star-forming galaxies
Comments: 28 pages, 13 figures, (abstract abridged). Accepted for publication in A&A
Submitted: 2012-10-09
We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.
[106]  [pdf] - 1151829
The Evolving Interstellar Medium of Star Forming Galaxies Since z=2 as Probed by Their Infrared Spectral Energy Distributions
Comments: To appear in ApJ, in the November 10, 2012, v759 -2 issue
Submitted: 2012-10-03, last modified: 2012-10-04
Using data from the mid-infrared to millimeter wavelengths for individual galaxies and for stacked ensembles at 0.5<z<2, we derive robust estimates of dust masses (Mdust) for main sequence (MS) galaxies, which obey a tight correlation between star formation rate (SFR) and stellar mass (M*), and for star-bursting galaxies that fall outside that relation. Exploiting the correlation of gas to dust mass with metallicity (Mgas/Mdust -Z), we use our measurements to constrain the gas content, CO-to-H2 conversion factors (a_co) and star formation efficiencies (SFE) of these distant galaxies. Using large statistical samples, we confirm that a_co and SFE are an order of magnitude higher and lower, respectively, in MS galaxies at high redshifts compared to the values of local galaxies with equivalently high infrared luminosities. For galaxies within the MS, we show that the variations of specific star formation rates (sSFR=SFR/M*) are driven by varying gas fractions. For relatively massive galaxies like those in our samples, we show that the hardness of the radiation field, <U>, which is proportional to the dust mass weighted luminosity (LIR/Mdust), and the primary parameter defining the shape of the SED, is equivalent to SFE/Z. For MS galaxies we measure this quantity, <U>, showing that it does not depend significantly on either the stellar mass or the sSFR. This is explained as a simple consequence of the existing correlations between SFR-M*, M*-Z and Mgas-SFR. Instead, we show that <U> (or LIR/Mdust) does evolve, with MS galaxies having harder radiation fields and thus warmer temperatures as redshift increases from z=0 to 2, a trend which can also be understood based on the redshift evolution of the M*-Z and SFR-M* relations. These results motivate the construction of a universal set of SED templates for MS galaxies which vary as a function of redshift with only one parameter, <U>.
[107]  [pdf] - 1151203
The Molecular Gas Content of z = 3 Lyman Break Galaxies; Evidence of a non Evolving Gas Fraction in Main Sequence Galaxies at z > 2
Comments: Accepted for publication in ApJL
Submitted: 2012-09-07, last modified: 2012-09-10
We present observations of the CO[3-2] emission towards two massive and infrared luminous Lyman Break Galaxies at z = 3.21 and z = 2.92, using the IRAM Plateau de Bure Interferometer, placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall properties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star formation rates ~1.6 and ~2.2 Gyr^(-1), despite their large infrared luminosities L_IR ~2-3 x 10^12 Lsun derived from Herschel. With one plausible CO detection (spurious detection probability of 10^(-3)) and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M*) with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z~2 is followed by a flattening at higher redshifts, providing supporting evidence for the existence of a plateau in the evolution of the specific star formation rate at z > 2.5.
[108]  [pdf] - 1151201
Spectral Energy Distributions of Type 1 AGN in the COSMOS Survey I - The XMM-COSMOS Sample
Comments: 22 pages, 22 figures, ApJ accepted, scheduled to be published October 20th, 2012, v758
Submitted: 2012-09-07, last modified: 2012-09-10
The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.
[109]  [pdf] - 1151021
A unified empirical model for infrared galaxy counts based on observed physical evolution of distant galaxies
Comments: 7 pages, 4 figures, 1 table, ApJL in press
Submitted: 2012-08-31, last modified: 2012-09-04
We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 {\mu}Jy. Material of the model (predictions, SED library, mock catalogs...) is available online at
[110]  [pdf] - 1123740
Millimeter imaging of submillimeter galaxies in the COSMOS field: Redshift distribution
Comments: 35 pages, 18 figures, 10 tables; accepted for publication in A&A
Submitted: 2012-05-29, last modified: 2012-06-14
We present new IRAM PdBI 1.3mm continuum observations at ~1.5" resolution of 28 SMGs previously discovered with the 870um bolometer LABOCA at APEX within the central 0.7deg2 of the COSMOS field. 19 out of the 28 LABOCA sources were detected with the PdBI at a >~3sigma level of ~1.4mJy/b. A combined analysis of this new sample with existing interferometrically identified SMGs in the COSMOS field yields the following results: 1) >~15%, and possibly up to ~40% of single-dish detected SMGs consist of multiple sources, 2) statistical identifications of multi-wavelength counterparts to the single-dish SMGs yield that only ~50% of these single-dish SMGs have real radio or IR counterparts, 3) ~18% of interferometric SMGs have only radio or even no multi-wavelength counterpart at all, and 4) ~50-70% of z>~3 SMGs have no radio counterparts down to an rms of 7-12uJy at 1.4GHz. Using the exact interferometric positions to identify proper multi-wavelength counterparts allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and the individual 1.1mm and 870um selected samples have a higher mean and broader width than the redshift distributions derived in previous studies. Our sample supports the previous tentative trend that on average brighter and/or mm-selected SMGs are located at higher redshifts. There is a tentative offset between the mean redshift for the 1.1mm (<z>=3.1+/-0.4) and 870um (<z>=2.6+/-0.4) selected samples, with the 1.1mm sources lying on average at higher redshifts. Based on our nearly complete sample of AzTEC 1.1mm SMGs within a uniform 0.15deg2 area we infer a higher surface density of z>~4 SMGs than predicted by current cosmological models. In summary, our findings imply that (sub-)millimeter interferometric identifications are crucial to build statistically complete and unbiased samples of SMGs.
[111]  [pdf] - 1123978
Dissecting the stellar mass-SFR correlation in z=1 star-forming disk galaxies
Comments: 6 pages, 4 figures, 1 table. ApJ Letters in press
Submitted: 2012-06-08
Using a mass-limited sample of 24um-detected, star-forming galaxies at 0.5<z<1.3, we study the mass-star formation rate (SFR) correlation and its tightness. The correlation is well defined (sigma=0.28dex) for disk galaxies (n_sersic<1.5), while more bulge-dominated objects often have lower specific SFRs. For disk galaxies, a much tighter correlation (sigma=0.19dex) is obtained if the rest-frame H-band luminosity is used instead of stellar mass derived from multicolor photometry. The specific SFR (sSFR) correlates strongly with rest-frame optical colors (hence luminosity-weighted stellar age) and also with clumpiness (which likely reflects the molecular gas fraction). This implies that most of the observed scatter is real, despite its low level, and not dominated by random measurement errors. After correcting for these differential effects a remarkably small dispersion remains (sigma=0.14dex), suggesting that measurement errors in mass or SFR are ~0.10dex, excluding systematic uncertainties. Measurement errors in stellar masses, the thickening of the correlation due to real sSFR variations, and varying completeness with stellar mass, can spuriously bias the derived slope to lower values due to the finite range over which observables (mass and SFR) are available. When accounting for these effects, the intrinsic slope for the main sequence for disk galaxies gets closer to unity.
[112]  [pdf] - 1117933
The hidden "AGN main sequence": Evidence for a universal black hole accretion to star formation rate ratio since z~2 producing a M_BH-M* relation
Comments: 6 Pages, 3 Figures, ApJL accepted after minor changes
Submitted: 2012-04-12, last modified: 2012-06-06
Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average SMBH growth rate follows remarkably similar trends with M* and redshift as the average star-formation rates (SFRs) of their host galaxies (i.e., (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to M* (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to M* ratio of ~10^{-3} as also inferred from today's M_BH-M_bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism. As such, we suggest that the same secular processes that drive the bulk of star formation are also responsible for the majority of SMBH growth. From this, we speculate that it is the availability of gas reservoirs that regulate both cosmological SMBH growth and star formation.
[113]  [pdf] - 1116331
The contribution of starbursts and normal galaxies to infrared luminosity functions at z < 2
Comments: 7 pages, 4 figures & 1 table. Accepted for publication in ApJL. Minor typos corrected in v2 following receipt of proofs
Submitted: 2012-02-01, last modified: 2012-02-14
We present a parameter-less approach to predict the shape of the infrared (IR) luminosity function (LF) at redshifts z < 2. It requires no tuning and relies on only three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate (sSFR) of main-sequence galaxies, and (3) the double-Gaussian decomposition of the sSFR-distribution at fixed stellar mass into a contribution (assumed redshift- and mass-invariant) from main-sequence and starburst activity. This self-consistent and simple framework provides a powerful tool for predicting cosmological observables: observed IR LFs are successfully matched at all z < 2, suggesting a constant or only weakly redshift-dependent contribution (8-14%) of starbursts to the star formation rate density. We separate the contributions of main-sequence and starburst activity to the global IR LF at all redshifts. The luminosity threshold above which the starburst component dominates the IR LF rises from log(LIR/Lsun) = 11.4 to 12.8 over 0 < z < 2, reflecting our assumed (1+z)^2.8-evolution of sSFR in main-sequence galaxies.
[114]  [pdf] - 1092537
Environmental effects in the interaction and merging of galaxies in zCOSMOS
Comments: 38 pages, 17 figures, submitted to ApJ
Submitted: 2011-12-20
(Abridged) We analyze the environments and galactic properties (morphologies and star-formation histories) of a sample of 153 close kinematic pairs in the redshift range 0.2 < z < 1 identified in the zCOSMOS-bright 10k spectroscopic sample of galaxies. Correcting for projection effects, the fraction of close kinematic pairs is three times higher in the top density quartile than in the lowest one. This translates to a three times higher merger rate because the merger timescales are shown, from mock catalogues based on the Millennium simulation, to be largely independent of environment once the same corrections for projection is applied. We then examine the morphologies and stellar populations of galaxies in the pairs, comparing them to control samples that are carefully matched in environment so as to remove as much as possible the well-known effects of environment on the properties of the parent population of galaxies. Once the environment is properly taken into account in this way, we find that the early-late morphology mix is the same as for the parent population, but that the fraction of irregular galaxies is boosted by 50-75%, with a disproportionate increase in the number of irregular-irregular pairs (factor of 4-8 times), due to the disturbance of disk galaxies. Future dry-mergers, involving elliptical galaxies comprise less than 5% of all close kinematic pairs. In the closest pairs, there is a boost in the specific star-formation rates of star-forming galaxies of a factor of 2-4, and there is also evidence for an increased incidence of post star-burst galaxies. Although significant for the galaxies involved, the "excess" star-formation associated with pairs represents only about 5% of the integrated star-formation activity in the parent sample. Although most pair galaxies are in dense environments, the effects of interaction appear to be largest in the lower density environments.
[115]  [pdf] - 1077452
GOODS-Herschel: The far-infrared view of star formation in AGN host galaxies since z~3
Comments: 21 pages; 14 Figures; MNRAS accepted; some figures have been degraded, full resolution version at:
Submitted: 2011-06-21, last modified: 2011-09-07
We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.
[116]  [pdf] - 1083856
GOODS-Herschel~: Gas-to-dust mass ratios and CO-to-H_2 conversion factors in normal and starbursting galaxies at high-z
Comments: Accepted for publication in ApJ
Submitted: 2011-09-06
We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results
[117]  [pdf] - 1083710
Dissecting Photometric redshift for Active Galactic Nuclei using XMM- and Chandra-COSMOS samples
Comments: ApJ, accepted for publication. Data also available at
Submitted: 2011-08-30
With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.
[118]  [pdf] - 1052243
The redshift and nature of AzTEC/COSMOS 1: A starburst galaxy at z=4.6
Comments: 7 pages, 4 figures, 1 table, accepted for publication in ApJL; one reference corrected
Submitted: 2011-02-21, last modified: 2011-08-15
Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.
[119]  [pdf] - 955613
The star formation history of mass-selected galaxies in the COSMOS field
Comments: 37 pages, 14 figures, 7 tables; accepted for publication in the Astrophysical Journal; High resolution versions of all figures available at
Submitted: 2010-11-29, last modified: 2011-01-18
We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.2<z<3 both populations show a strong and M_*-independent decrease in their SSFR towards z=0.2, best described by a power- law (1+z)^n, where n~4.3 for all galaxies and n~3.5 for SF sources. The decrease appears to have started at z>2, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples.
[120]  [pdf] - 1026760
The VLA-COSMOS Survey. IV. Deep Data and Joint Catalog
Comments: 34 pages, 20 figures (3 at reduced resolution). Accepted for publication in ApJS.
Submitted: 2010-05-10
In the context of the VLA-COSMOS Deep project additional VLA A array observations at 1.4 GHz were obtained for the central degree of the COSMOS field and combined with the existing data from the VLA-COSMOS Large project. A newly constructed Deep mosaic with a resolution of 2.5" was used to search for sources down to 4 sigma with 1 sigma ~ 12 microJy/beam in the central 50'x50'. This new catalog is combined with the catalog from the Large project (obtained at 1.5"x1.4" resolution) to construct a new Joint catalog. All sources listed in the new Joint catalog have peak flux densities of >5 sigma at 1.5" and/or 2.5" resolution to account for the fact that a significant fraction of sources at these low flux levels are expected to be slighty resolved at 1.5" resolution. All properties listed in the Joint catalog such as peak flux density, integrated flux density and source size are determined in the 2.5" resolution Deep image. In addition, the Joint catalog contains 43 newly identified multi-component sources.
[121]  [pdf] - 1025865
No Evolution in the IR-Radio Relation for IR-Luminous Galaxies at z<2 in the COSMOS Field
Comments: 9 pages, 4 figures. Accepted for publication in ApJL.
Submitted: 2010-03-22
Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.
[122]  [pdf] - 1025776
The Opacity of Galactic Disks at z~0.7
Comments: 7 pages, 3 figures. Accepted for publication in ApJL.
Submitted: 2010-03-17
We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.
[123]  [pdf] - 1019069
The VLA-COSMOS Perspective on the IR-Radio Relation. I. New Constraints on Selection Biases and the Non-Evolution of the IR/Radio Properties of Star Forming and AGN Galaxies at Intermediate and High Redshift
Comments: 52 pages, 23 figures (11 at reduced resolution). Accepted for publication in ApJS
Submitted: 2010-01-11
VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.
[124]  [pdf] - 1018865
The X-ray to optical-UV luminosity ratio of X-ray selected Type 1 AGN in XMM-COSMOS
Comments: 14 pages, 15 figures, to appear in A&A
Submitted: 2009-12-21
We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.
[125]  [pdf] - 1018249
The Build-Up of the Hubble Sequence in the COSMOS Field
Comments: 5 pages, 3 figures, submitted to ApJL
Submitted: 2009-11-06
We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.
[126]  [pdf] - 901680
The 10k zCOSMOS: morphological transformation of galaxies in the group environment since z~1
Comments: 21 pages, 13 figures, submitted to ApJ
Submitted: 2009-09-11
We study the evolution of galaxies inside and outside of the group environment since z=1 using a large well defined set of groups and galaxies from the zCOSMOS-bright redshift survey in the COSMOS field. The fraction of galaxies with early-type morphologies increases monotonically with M_B luminosity and stellar mass and with cosmic epoch. It is higher in the groups than elsewhere, especially at later epochs. The emerging environmental effect is superposed on a strong global mass-driven evolution, and at z~0.5 and log(M*/Msol)~10.2, the "effect" of group environment is equivalent to (only) about 0.2 dex in stellar mass or 2 Gyr in time. The stellar mass function of galaxies in groups is enriched in massive galaxies. We directly determine the transformation rates from late to early morphologies, and for transformations involving colour and star formation indicators. The transformation rates are systematically about twice as high in the groups as outside, or up to 3-4 times higher correcting for infall and the appearance of new groups. The rates reach values, for masses around the crossing mass 10^10.5 Msol, as high as (0.3-0.7)/Gyr in the groups, implying transformation timescales of 1.4-3 Gyr, compared with less than 0.2/Gyr, i.e. timescales >5 Gyr, outside of groups. All three transformation rates decrease at higher stellar masses, and must decrease also at the lower masses below 10^10 Msol which we cannot well probe. The rates involving colour and star formation are consistently higher than those for morphology, by a factor of about 50%. Our conclusion is that the transformations which drive the evolution of the overall galaxy population since z~1 must occur at a rate 2-4 times higher in groups than outside of them.
[127]  [pdf] - 20075
The Evolution of the Star Formation of zCOSMOS and SDSS galaxies at z<0.7 as a Function of Mass and Structural Parameters
Comments: Proceedings of the Subaru International Conference: "Panoramic Views of Galaxy Formation and Evolution", Hayama, Japan
Submitted: 2009-01-07, last modified: 2009-09-10
We present in these proceedings some preliminary results we have obtained studying the evolution of the specific star formation rate as a function of surface mass density and Sersic indices at z<0.7. These results are based on the consistent comparison of the properties of ~ 650 massive zCOSMOS galaxies in a mass-complete sample at 0.5<z<0.7 with a mass-complete sample of ~ 21500 SDSS local galaxies.
[128]  [pdf] - 1001486
Spitzer 70 and 160-micron Observations of the COSMOS Field
Comments: Accepted AJ, 15 Aug. 2009. Data available at until released by IRSA
Submitted: 2009-02-18, last modified: 2009-08-25
We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band.
[129]  [pdf] - 315454
The Dependence of Star Formation Activity on Stellar Mass Surface Density and Sersic Index in zCOSMOS Galaxies at 0.5<z<0.9 Compared with SDSS Galaxies at 0.04<z<0.08
Comments: Published 2009 in ApJ, 694, 1099
Submitted: 2009-01-05, last modified: 2009-06-29
One of the key unanswered questions in the study of galaxy evolution is what physical processes inside galaxies drive the changes in the SFRs in individual galaxies that, taken together, produce the large decline in the global star-formation rate density (SFRD) to redshifts since z~2. Many studies of the SFR at intermediate redshifts have been made as a function of the integrated stellar mass of galaxies but these did not use information on the internal structural properties of the galaxies. In this paper we present a comparative study of the dependence of SFRs on the average surface mass densities (SigmaM) of galaxies of different morphological types up to z~1 using the zCOSMOS and SDSS surveys. The main findings about the evolution of these relatively massive galaxies are: 1) There is evidence that, for both SDSS ans zCOSMOS galaxies, the mean specific SFR within a given population (either disk-dominated or bulge-dominated) is independent of SigmaM; 2) The observed SSFR - SigmaM step-function relation is due, at all investigated redshifts, to the changing mix of disk-dominated and bulge-dominated galaxies as surface density increases and the strong difference in the average SSFR between disks and bulges. We also find a modest differential evolution in the size-mass relations of disk and spheroid galaxies; 3) The shape of the median SSFR - SigmaM relation is similar, but with median SSFR values that are about 5-6 times higher in zCOSMOS galaxies than for SDSS, across the whole range of SigmaM, and in both spheroid and disk galaxies. This increase matches that of the global SFRD of the Universe as a whole, emphasizing that galaxies of all types are contributing, proportionally, to the global increase in SFRD in the Universe back to these redshifts (abridged).
[130]  [pdf] - 900434
Chasing highly obscured QSOs in the COSMOS field
Comments: Paper submitted to ApJ. This is the revised version after all referee's comments have been addressed
Submitted: 2008-10-03
(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.
[131]  [pdf] - 88875
The redshift evolution of early-type galaxies in COSMOS: Do massive early-type galaxies form by dry mergers?
Comments: Replaced with accepted version (ApJ COSMOS special issue)
Submitted: 2007-01-25, last modified: 2007-02-07
ABRIDGED: We study the evolution since z~1 of the rest-frame B luminosity function of the early-type galaxies (ETGs) in ~0.7 deg^2 in the COSMOS field. In order to identify ALL progenitors of local ETGs we construct the sample of high-z galaxies using two complementary criteria: (i) A morphological selection based on the Zurich Estimator of Structural Types, and (ii) A photometric selection based on the galaxy properties in the (U-V)-M_V color-magnitude diagram. We furthermore constrain both samples so as to ensure that the selected progenitors of ETGs are compatible with evolving into systems which obey the mu_B-r_{hl} Kormendy relation. Assuming the luminosity evolution derived from studies of the fundamental plane for high-z ETGs, our analysis shows no evidence for a decrease in the number density of the most massive ETGs out to z~ 0.7: Both the morphologically- and the photometrically-selected sub-samples show no evolution in the number density of bright (~L>2.5L*) ETGs. Allowing for different star formation histories, and cosmic variance, we estimate a maximum decrease in the number density of massive galaxies at that redshift of ~30%. We observe, however, in both the photometrical and morphological samples, a deficit of up to ~2-3 of fainter ETGs over the same cosmic period. Our results argue against a significant contribution of recent dissipationless ``dry'' mergers to the formation of the most massive ETGs. We suggest that the mass growth in low luminosity ETGs can be explained with a conversion from z~0.7 to z=0 of blue, irregular and disk galaxies into low- and intermediate-mass ``red'' ETGs, possibly also through gas rich mergers.
[132]  [pdf] - 87030
COSMOS morphological classification with ZEST (the Zurich Estimator of Structural Types) and the evolution since z=1 of the Luminosity Function of early-, disk-, and irregular galaxies
Comments: Accepted for publication in the ApJ COSMOS special issue. A version with high resolution figures is available at
Submitted: 2006-11-20
(ABRIDGED) Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) Analysis. We use a sample of ~56,000 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the Luminosity Functions of COSMOS galaxies of early, disk and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure-luminosity evolution (of about 0.95 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.
[133]  [pdf] - 1938974
The evolution of the number density of large disk galaxies in COSMOS
Comments: Accepted for publication in the ApJ COSMOS special issue. A version with figures in higher resolution is available at
Submitted: 2006-09-01
We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.