We examine the behaviour of the infrared-radio correlation (IRRC) over the range 0 <z ≲ 6 using new, highly sensitive 3 GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared ...data from the Herschel Space Observatory in the 2 deg2 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.4 GHz radio luminosity ratio (qTIR), decreases with increasing redshift: qTIR(z) = (2.88 ± 0.03)(1 + z)− 0.19 ± 0.01. This is consistent with several previous results from the literature. Moderate-to-high radiative luminosity AGN do not follow the same qTIR(z) trend as star-forming galaxies, 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 qTIR(z) trend of the star-forming galaxy population. We demonstrate that the choice of the average radio spectral index directly affects the normalisation, as well as the derived trend with redshift of the IRRC. An increasing fractional contribution to the observed 3 GHz 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 trend of decreasing qTIR with increasing redshift. A more detailed understanding of the radio spectrum is therefore required for robust K corrections in the radio and to fully understand the cosmic evolution of the IRRC. Lastly, we present a redshift-dependent relation between rest-frame 1.4 GHz radio luminosity and star formation rate taking the derived redshift trend into account.
We study the composition of the faint radio population selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3 GHz Large Project, which is a radio continuum survey ...performed at 10 cm wavelength. The survey covers a 2.6 square degree area with a mean rms of ~ 2.3 μJy/beam, cataloging 10 830 sources above 5σ, and enclosing the full 2 square degree COSMOS field. By combining these radio data with optical, near-infrared (UltraVISTA), and mid-infrared (Spitzer/IRAC) data, as well as X-ray data (Chandra), we find counterparts to radio sources for ~93% of the total radio sample reaching out to z ≲ 6; these sources are found in the unmasked areas of the COSMOS field, i.e., those not affected by saturated or bright sources in the optical to near-infrared (NIR) bands. We further classify the sources as star-forming galaxies or AGN based on various criteria, such as X-ray luminosity; observed mid-infrared color; UV–far-infrared spectral energy distribution; rest-frame, near-UV optical color that is corrected for dust extinction; and radio excess relative to that expected from the star formation rate of the hosts. We separate the AGN into subsamples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, i.e., candidates for high-redshift analogs to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these subpopulations down to radio flux levels of ~11 μJy at 3 GHz (or ~20 μJy at 1.4 GHz assuming a spectral index of –0.7). We find that the dominant fraction at 1.4 GHz flux densities above ~200 μJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100 μJy the fraction of star-forming galaxies increases to ~ 60%, followed by the moderate-to-high radiative luminosity AGN (~ 20%) and low-to-moderate radiative luminosity AGN (~ 20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the Square Kilometer Array (SKA). Our estimates suggest that at the faint flux limits to be reached by the (Wide, Deep, and UltraDeep) SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.
We study the incidence of nuclear obscuration on a complete sample of 1310 active galactic nuclei (AGN) selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in ...the redshift range 0.3 < z < 3.5. We classify the AGN as obscured or unobscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70 per cent of the objects, and the remaining 30 per cent can be further subdivided into two distinct classes: at low luminosities X-ray unobscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high-luminosities broad-line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and unobscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.
We investigate the star formation properties of ∼800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (∼0.1–4) and about four orders of ...magnitude in star formation rate (SFR). It includes both star-forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet (RQ) and radio-loud objects. We compare the SFR derived from the far-infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFG population and the occurrence of AGNs in sources with different level of star formation. Finally, we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.
We use deep Herschel observations taken with both PACS and SPIRE imaging cameras to estimate the dust mass of a sample of galaxies extracted from the GOODS-S, GOODS-N and the COSMOS fields. We divide ...the redshift–stellar mass (Mstar)–star formation rate (SFR) parameter space into small bins and investigate average properties over this grid. In the first part of the work we investigate the scaling relations between dust mass, stellar mass and SFR out to z = 2.5. No clear evolution of the dust mass with redshift is observed at a given SFR and stellar mass. We find a tight correlation between the SFR and the dust mass, which, under reasonableassumptions, is likely a consequence of the Schmidt-Kennicutt (S-K) relation. The previously observed correlation between the stellar content and the dust content flattens or sometimes disappears when considering galaxies with the same SFR. Our finding suggests that most of the correlation between dust mass and stellar mass obtained by previous studies is likely a consequence of the correlation between the dust mass and the SFR combined with the main sequence, i.e., the tight relation observed between the stellar mass and the SFR and followed by the majority of star-forming galaxies. We then investigate the gas content as inferred from dust mass measurements. We convert the dust mass into gas mass by assuming that the dust-to-gas ratio scales linearly with the gas metallicity (as supported by many observations). For normal star-forming galaxies (on the main sequence) the inferred relation between the SFR and the gas mass (integrated S-K relation) broadly agrees with the results of previous studies based on CO measurements, despite the completely different approaches. We observe that all galaxies in the sample follow, within uncertainties, the same S-K relation. However, when investigated in redshift intervals, the S-K relation shows a moderate, but significant redshift evolution. The bulk of the galaxy population at z ~ 2 converts gas into stars with an efficiency (star formation efficiency, SFE = SFR/Mgas, equal to the inverse of the depletion time) about 5 times higher than at z ~ 0. However, it is not clear what fraction of such variation of the SFE is due to an intrinsic redshift evolution and what fraction is simply a consequence of high-z galaxies having, on average, higher SFR, combined with the super-linear slope of the S-K relation (while other studies find a linear slope). We confirm that the gas fraction (fgas = Mgas/(Mgas + Mstar)) decreases with stellar mass and increases with the SFR. We observe no evolution with redshift once Mstarand SFR are fixed. We explain these trends by introducing a universal relation between gas fraction, stellar mass and SFR that does not evolve with redshift, at least out to z ~ 2.5. Galaxies move across this relation as their gas content evolves across the cosmic epochs. We use the 3D fundamental fgas–Mstar–SFR relation, along with the evolution of the main sequence with redshift, to estimate the evolution of the gas fraction in the average population of galaxies as a function of redshift and as a function of stellar mass: we find that Mstar ≳ 1011 M⊙ galaxies show the strongest evolution at z ≳ 1.3 and a flatter trend at lower redshift, while fgas decreases more regularly over the entire redshift range probed in Mstar ≲ 1011 M⊙ galaxies, in agreement with a downsizing scenario.
We use Herschel 70 to 160 μm images to study the size of the far-infrared emitting region in about 400 local galaxies and quasar (QSO) hosts. The sample includes normal “main-sequence” star-forming ...galaxies, as well as infrared luminous galaxies and Palomar-Green QSOs, with different levels and structures of star formation. Assuming Gaussian spatial distribution of the far-infrared (FIR) emission, the excellent stability of the Herschel point spread function (PSF) enables us to measure sizes well below the PSF width, by subtracting widths in quadrature. We derive scalings of FIR size and surface brightness of local galaxies with FIR luminosity, with distance from the star-forming main-sequence, and with FIR color. Luminosities LFIR~ 1011 L⊙ can be reached with a variety of structures spanning 2 dex in size. Ultraluminous LFIR≳ 1012 L⊙ galaxies far above the main-sequence inevitably have small Re,70~ 0.5 kpc FIR emitting regions with large surface brightness, and can be close to optically thick in the FIR on average over these regions. Compared to these local relations, first ALMA sizes for the dust emission regions in high redshift galaxies, measured at somewhat longer rest wavelengths, suggest larger sizes at the same IR luminosity. We report a remarkably tight relation with 0.15 dex scatter between FIR surface brightness and the ratio of Cii 158 μm emission and FIR emission – the so-called Cii-deficit is more tightly linked to surface brightness than to FIR luminosity or FIR color. Among 33 z ≤ 0.1 PG QSOs with typical LFIR/LBol,AGN ≈ 0.1, 19 have a measured 70 μm half light radius, with median Re,70 = 1.1 kpc. This is consistent with the FIR size for galaxies with similar LFIR but lacking a QSO, in accordance with a scenario where the rest FIR emission of these types of QSOs is, in most cases, due to host star formation.
Combining the deepest Herschel extragalactic surveys (PEP, GOODS-H, HerMES), and Monte Carlo mock catalogs, we explore the robustness of dust mass estimates based on modeling of broadband spectral ...energy distributions (SEDs) with two popular approaches: Draine & Li (2007, ApJ, 657, 810; DL07) and a modified blackbody (MBB). We analyze the cause, drivers, and trends of uncertainties and systematics in thorough detail. As long as the observed SED extends to at least 160−200 μm in the rest frame, Mdust can be recovered with a >3σ significance and without the occurrence of systematics. An average offset of a factor ~1.5 exists between DL07- and MBB-based dust masses, based on consistent dust properties. The performance of DL07 modeling turns out to be more robust than that of MBB since relative errors on Mdust are more mildly dependent on the maximum covered rest-frame wavelength and are less scattered. At the depth of the deepest Herschel surveys (in the GOODS-S field), it is possible to retrieve dust masses with a signal-to-noise ratio, S/N ≥ 3 for galaxies on the main sequence of star formation (MS) down to M∗ ~ 1010 M⊙ up to z ~ 1. At higher redshift (z ≤ 2), the same result is only achieved for objects at the tip of the MS or for those objects lying above the tip owing to sensitivity and wavelength coverage limitations. Molecular gas masses, obtained by converting Mdust through the metallicity-dependent gas-to-dust ratio δGDR, are consistent with those based on the scaling of depletion time, τdep, and on CO sub-mm spectroscopy. Focusing on CO-detected galaxies at z> 1, the δGDR dependence on metallicity is consistent with the local relation, provided that a sufficient SED coverage is available. Once we established that Herschel-only and sub-mm-only estimates of dust masses can be affected by large uncertainties and possibly systematics in some cases, we combined far-IR Herschel data and sub-mm ALMA expected fluxes to study the advantages of a full SED coverage. The uncertainty on Mdust reduces to <30% for more than 85% of Herschel galaxies, thus potentially facilitating a fast statistical study of Mdust,gas for large samples, at least up to z ~ 2.
Using spectroscopic observations taken for the Visible Multi-Object Spectrograph (VIMOS) Ultra-Deep Survey (VUDS) we report here on the discovery of PCl J1001+0220, a massive proto-cluster of ...galaxies located at zspec ~ 4.57 in the COSMOS field. With nine spectroscopic members, the proto-cluster was initially detected as a ~12σ spectroscopic overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe (2 < z ≲ 6) performed by VUDS. It was further mapped using a new technique developed which statistically combines spectroscopic and photometric redshifts, the latter derived from a recent compilation of incredibly deep multi-band imaging performed on the COSMOS field. Through various methods, the descendant mass of PCl J1001+0220 is estimated to be log (Mh/M⊙)z=0 $\log(\mathcal{M}_{h}/\mathcal{M}_{\odot})_{z=0}$log(Mh/M⊙)z=0 ~ 14.5–15 with a large amount of mass apparently already in place at z ~ 4.57. An exhaustive comparison was made between the properties of various spectroscopic and photometric member samples and matched samples of galaxies inhabiting less dense environments at the same redshifts. Tentative evidence is found for a fractional excess of older galaxies more massive in their stellar content amongst the member samples relative to the coeval field, an observation which suggests the pervasive early onset of vigorous star formation for proto-cluster galaxies. No evidence is found for the differences in the star formation rates (SFRs) of member and coeval field galaxies either through estimating by means of the rest-frame ultraviolet or through separately stacking extremely deep Very Large Array 3 GHz imaging for both samples. Additionally, no evidence for pervasive strong active galactic nuclei (AGN) activity is observed in either environment. Analysis of Hubble Space Telescope images of both sets of galaxies as well as their immediate surroundings provides weak evidence for an elevated incidence of galaxy–galaxy interaction within the bounds of the proto-cluster. The stacked and individual spectral properties of the two samples are compared, with a definite suppression of Lyα seen in the average member galaxy relative to the coeval field ( fesc, Lyα = 1.8−1.7+0.3 $f_{esc, \, \textrm{Ly}\alpha} = 1.8^{+0.3}_{-1.7}$fesc, Lyα=1.8−1.7+0.3% and 4.0−0.8+1.0 $4.0^{+1.0}_{-0.8}$4.0−0.8+1.0%, respectively). This observation along with other lines of evidence leads us to infer the possible presence of a large, cool, diffuse medium within the proto-cluster environment evocative of a nascent intracluster medium forming in the early universe.
We compare various star formation rate (SFR) indicators for star-forming galaxies at 1.4 < z < 2.5 in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those ...from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR–stellar mass (M
*) relation, compared to that of the dominant main-sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR–M
* relation is ∼0.8–0.9, in agreement with that derived from the dust-corrected UV luminosity. Exploiting deeper 24 μm Spitzer data, we have characterized a subsample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the B band. The combination of Herschel with Spitzer data has allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the Hα luminosity of a sample of star-forming galaxies at 1.4 < z < 1.7. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.
Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the
Herschel
Astrophysical Large Area Survey with
S
500
μ
m
...≥ 80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range 2.08 <
z
< 4.05 with a median value of
z
= 2.9 ± 0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ∼140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700 ± 300 km s
−1
and a median of 800 km s
−1
. We analyze the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous (
L
FIR
> 10
13
L
⊙
) galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-
z
galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest
Herschel
galaxies.