We explore the cosmic evolution of radio luminous active galactic nuclei (AGNs) with low radio powers (L 1.4 GHz 5 X 1025 W Hz-1) out to z = 1.3 using to date the largest sample of ~600 ...low-luminosity radio AGN at intermediate redshift drawn from the VLA-COSMOS survey. We derive the radio-luminosity function for these AGNs, and its evolution with cosmic time assuming two extreme cases: (1) pure luminosity and (2) pure density evolution. The former and latter yield L * (1 + z)0.8 ± 0.1, and * (1 + z)1.1 ± 0.1, respectively, both implying a fairly modest change in properties of low-radio-power AGNs since z = 1.3. We show that this is in stark contrast with the evolution of powerful (L 1.4 GHz > 5 X 1025 W Hz-1) radio AGN over the same cosmic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by Willot et al. We demonstrate that this can be explained through differences in black hole fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful AGNs. Our findings suggest that high- and low-radio-power AGN activities are triggered in different stages during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive galaxies already at z ~ 1, and they may significantly contribute to the heating of their surrounding medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the 'radio mode' in cosmological models.
We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the Atacama Large Millimeter/submillimeter Array (ALMA) public archive in the ...COSMOS deep field (A3COSMOS). Our A3COSMOS data set contains ∼700 galaxies (0.3 z 6) with high-confidence ALMA detections in the (sub)millimeter continuum and multiwavelength spectral energy distributions. Multiple gas mass calibration methods are compared, and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans tail continuum) have been tested. Combining our A3COSMOS sample with ∼1000 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 ( ) each as a function of the stellar mass ( ), offset from the star-forming main sequence ( ) 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 semianalytic modeling. These together provide a coherent picture of cold molecular gas, star formation rate, and stellar mass evolution in galaxies across cosmic time.
We present the VLA-COSMOS 3 GHz Large Project based on 384 h 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 μJy beam-1 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σ, 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/N3 GHz > 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 α = −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.
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 identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on
Spitzer
3.6
μ
m images. At the simplest level, we ...distinguish five environments: centres, bars, spiral arms, interarm regions, and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, which are publicly released but not explicitly used in this paper. We examine trends with environment in the molecular gas content, star formation rate, and depletion time using PHANGS–ALMA CO(2–1) intensity maps and tracers of star formation. The interarm regions and discs without strong spirals clearly dominate in area, whereas molecular gas and star formation are quite evenly distributed among the five basic environments. We reproduce the molecular Kennicutt–Schmidt relation with a slope compatible with unity within the uncertainties and without significant slope differences among environments. In contrast to what has been suggested by early studies, we find that bars are not always deserts devoid of gas and star formation, but instead they show large diversity. Similarly, spiral arms do not account for most of the gas and star formation in disc galaxies, and they do not have shorter depletion times than the interarm regions. Spiral arms accumulate gas and star formation, without systematically boosting the star formation efficiency. Centres harbour remarkably high surface densities and on average shorter depletion times than other environments. Centres of barred galaxies show higher surface densities and wider distributions compared to the outer disc; yet, depletion times are similar to unbarred galaxies, suggesting highly intermittent periods of star formation when bars episodically drive gas inflow, without enhancing the central star formation efficiency permanently. In conclusion, we provide quantitative evidence that stellar structures in galaxies strongly affect the organisation of molecular gas and star formation, but their impact on star formation efficiency is more subtle.
ABSTRACT We present high-resolution (0 16) 870 m Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous ( ) submillimeter galaxies (SMGs) from the ALESS survey of the Extended ...Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these galaxies on ∼1.3 kpc scales. The emission has a median effective radius of Re = 0 24 0 02, corresponding to a typical physical size of 1.8 0.2 kpc. We derive a median Sérsic index of n = 0.9 0.2, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0 12 (∼1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate signal-to-noise ratio interferometric data. We compare our maps to comparable-resolution Hubble Space Telescope -band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive ( 1-2 × 1011 ) early-type galaxies observed locally.
We present subkiloparsec-scale mapping of the 870 m ALMA continuum emission in six luminous (LIR ∼ 5 × 1012 L ) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field ...South. Our high-fidelity 0 07-resolution imaging (∼500 pc) reveals robust evidence for structures with deconvolved sizes of 0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the "ring" and "bar" radii (1.9 0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q < 1). The inferred one-dimensional velocity dispersions ( r 70-160 km s−1) are marginally consistent with the limits implied if the sizes of the largest structures are comparable to the Jeans length. We create maps of the star formation rate density ( SFR) on ∼500 pc scales and show that the SMGs are able to sustain a given (galaxy-averaged) SFR over much larger physical scales than local (ultra)luminous infrared galaxies. However, on 500 pc scales, they do not exceed the Eddington limit set by radiation pressure on dust. If confirmed by kinematics, the potential presence of nonaxisymmetric structures would provide a means for net angular momentum loss and efficient star formation, helping to explain the very high star formation rates measured in SMGs.
We present EMPIRE, an IRAM 30 m large program that mapped λ = 3-4 mm dense gas tracers at ∼1-2 kpc resolution across the whole star-forming disk of nine nearby massive spiral galaxies. We describe ...the EMPIRE observing and reduction strategies and show new whole-galaxy maps of HCN(1−0), HCO+(1−0), HNC(1−0), and CO(1−0). We explore how the HCN-to-CO and IR-to-HCN ratios, observational proxies for the dense gas fraction and dense gas star formation efficiency, depend on host galaxy and local environment. We find that the fraction of dense gas correlates with stellar surface density, gas surface density, molecular-to-atomic gas ratio, and dynamical equilibrium pressure. In EMPIRE, the star formation rate per unit dense gas is anticorrelated with these same environmental parameters. Thus, although dense gas appears abundant in the central regions of many spiral galaxies, this gas appears relatively inefficient at forming stars. These results qualitatively agree with previous work on nearby galaxies and the Milky Way's Central Molecular Zone. To first order, EMPIRE demonstrates that the conditions in a galaxy disk set the gas density distribution and that the dense gas traced by HCN shows an environment-dependent relation to star formation. However, our results also show significant ( 0.2 dex) galaxy-to-galaxy variations. We suggest that gas structure below the scale of our observations and dynamical effects likely also play an important role.
We study the spectral energy distribution (SED) of the radio continuum (RC) emission from the Key Insight in Nearby Galaxies Emitting in Radio (KINGFISHER) sample of nearby galaxies to understand the ...energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4, 4.8, 8.4, and 10.5 GHz combined with archive data allow us, for the first time, to determine the mid-RC (1-10 GHz, MRC) bolometric luminosities and further present calibration relations versus the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo technique leading to measurements for the nonthermal spectral index ( ) and the thermal fraction ( ) with mean values of for the total spectral index) and = (10 9)% at 1.4 GHz. The MRC luminosity changes over ∼3 orders of magnitude in the sample, MRC . The thermal emission is responsible for ∼23% of the MRC on average. We also compare the extinction-corrected diagnostics of the star-formation rate (SFR) with the thermal and nonthermal radio tracers and derive the first star-formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing SFR surface density, indicating the effect of the star-formation feedback on the cosmic-ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with SFR, due to the amplification of the magnetic fields in star-forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star-forming galaxies are detected.