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.
Based on a sample of over 1800 radio AGN at redshifts out to z ~ 5, which have typical stellar masses within ~3 × (1010 − 1011)M⊙, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz ...radio luminosity functions for radio AGN (L1.4 GHz~ 1022 − 1027 W Hz-1) out to z ~ 5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Φ∗ ∝ (1 + z)(2.00 ± 0.18) − (0.60 ± 0.14)z, and L∗ ∝ (1 + z)(2.88 ± 0.82) − (0.84 ± 0.34)z, respectively, with a turnover in number and luminosity densities of the population at z ≈ 1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z ~ 5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.
The pre-merging system of galaxy clusters Abell 3391-Abell 3395 located at a mean redshift of 0.053 has been observed at 1 GHz in an ASKAP/EMU Early Science observation as well as in X-rays with ...eROSITA. The projected separation of the X-ray peaks of the two clusters is ~50′ or ~3.1 Mpc. Here we present an inventory of interesting radio sources in this field around this cluster merger. While the eROSITA observations provide clear indications of a bridge of thermal gas between the clusters, neither ASKAP nor MWA observations show any diffuse radio emission coinciding with the X-ray bridge. We derive an upper limit on the radio emissivity in the bridge region of 〈
J
〉
1 GHz
< 1.2 × 10
−44
W Hz
−1
m
−3
. A non-detection of diffuse radio emission in the X-ray bridge between these two clusters has implications for particle-acceleration mechanisms in cosmological large-scale structure. We also report extended or otherwise noteworthy radio sources in the 30 deg
2
field around Abell 3391-Abell 3395. We identified 20 Giant Radio Galaxies, plus 7 candidates, with linear projected sizes greater than 1 Mpc. The sky density of field radio galaxies with largest linear sizes of >0.7 Mpc is ≈1.7 deg
−2
, three times higher than previously reported. We find no evidence for a cosmological evolution of the population of Giant Radio Galaxies. Moreover, we find seven candidates for cluster radio relics and radio halos.
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.
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.
We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs), predominantly from the AS2UDS survey at z 2, with bright rest-frame optical counterparts ( ) to ...investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio ( ) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii = 0.6). This ratio does not change with specific star formation rate over the factor of 30 spanned by our targets, sampling the locus of "normal" main-sequence galaxies up to the starburst regime, . Unlike typical spiral galaxies in the local universe, our results imply that massive SMGs are experiencing centrally enhanced star formation. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z = 1.5, which is consistent with these systems being the descendants of z 2 SMGs.
We have studied a sample of 1604 moderate-to-high radiative luminosity active galactic nuclei (HLAGN) selected at 3 GHz within the VLA-COSMOS 3 GHz Large Project. These were classified by combining ...multiple AGN diagnostics: X-ray data, mid-infrared data and broadband spectral energy distribution fitting. We decomposed the total radio 1.4 GHz luminosity (L1.4 GHz, TOT) into the emission originating from star formation and AGN activity by measuring the excess in L1.4 GHz, TOT relative to the infrared-radio correlation of star-forming galaxies. To quantify the excess, for each source we calculated the AGN fraction (fAGN) defined as the fractional contribution of AGN activity to L1.4 GHz, TOT. The majority of the HLAGN, (68.0 ± 1.5)%, are dominated by star-forming processes (fAGN ≤ 0.5), while (32.0 ± 1.5)% are dominated by AGN-related radio emission (0.5 < fAGN ≤ 1). We used the AGN-related 1.4 GHz emission to derive the 1.4 GHz AGN luminosity functions of HLAGN. By assuming pure density and pure luminosity evolution models we constrained their cosmic evolution out to z ∼ 6, finding Φ*(z)∝(1 + z)(2.64 ± 0.10)+(−0.61 ± 0.04)z and L*(z)∝(1 + z)(3.97 ± 0.15)+(−0.92 ± 0.06)z. These evolutionary laws show that the number and luminosity density of HLAGN increased from higher redshifts (z ∼ 6) up to a maximum in the redshift range 1 < z < 2.5, followed by a decline toward local values. By scaling the 1.4 GHz AGN luminosity to kinetic luminosity using the standard conversion, we estimate the kinetic luminosity density as a function of redshift. We compare our result to the semi-analytic models of radio mode feedback, and find that this feedback could have played an important role in the context of AGN-host co-evolution in HLAGN which shows evidence of AGN-related radio emission (fAGN > 0).
To better constrain the physical mechanisms driving star formation, 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 3 GHz map (noise rms = 2.3 μJy beam−1, θbeam = 0.75 arcsec) to construct a mass-complete sample of 3184 radio-selected SFGs that reside on and above the main sequence (MS) of SFGs. We constrain the overall extent of star formation activity in galaxies by applying a 2D Gaussian model to their radio continuum emission. Extensive Monte Carlo simulations are used to validate the robustness of our measurements and characterize the selection function. We find no clear dependence between the radio size and stellar mass, M⋆, of SFGs with 10.5 ≲ log(M⋆/M⊙) ≲ 11.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 Reff = 1.5 ± 0.2 (1.0 ± 0.2) kpc remains nearly constant with cosmic time; a parametrization of the form Reff ∝ (1 + z)α yields a shallow slope of only α = −0.26 ± 0.08 (0.12 ± 0.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 ∼2 (1.3) at z = 0.5 (2), indicating star formation is enhanced at small radii. The galactic-averaged star formation rate surface density (ΣSFR) scales with the distance to the MS, except for a fraction of MS galaxies (≲10%) that harbor starburst-like ΣSFR. These “hidden” starbursts might have experienced a compaction phase due to disk instability and/or a merger-driven burst of star formation, which may or may not significantly offset a galaxy from the MS. We thus propose to use ΣSFR and distance to the MS in conjunction to better identify the galaxy population undergoing a starbursting phase.
We construct the average radio spectral energy distribution (SED) of highly star-forming galaxies (HSFGs) up to z ∼ 4. Infrared and radio luminosities are bound by a tight correlation that is defined ...by the so-called q parameter. This infrared–radio correlation provides the basis for the use of radio luminosity as a star-formation tracer. Recent stacking and survival analysis studies find q to be decreasing with increasing redshift. It was pointed out that a possible cause of the redshift trend could be the computation of rest-frame radio luminosity via a single power-law assumption of the star-forming galaxies’ (SFGs) SED. To test this, we constrained the shape of the radio SED of a sample of HSFGs. To achieve a broad rest-frame frequency range, we combined previously published Very Large Array observations of the COSMOS field at 1.4 GHz and 3 GHz with unpublished Giant Meterwave Radio Telescope (GMRT) observations at 325 MHz and 610 MHz by employing survival analysis to account for non-detections in the GMRT maps. We selected a sample of HSFGs in a broad redshift range (z ∈ 0.3, 4, SFR ≥ 100 M⊙ yr−1) and constructed the average radio SED. By fitting a broken power-law, we find that the spectral index changes from α1 = 0.42 ± 0.06 below a rest-frame frequency of 4.3 GHz to α2 = 0.94 ± 0.06 above 4.3 GHz. Our results are in line with previous low-redshift studies of HSFGs ( SFR > 10 M⊙ yr−1) that show the SED of HSFGs to differ from the SED found for normal SFGs ( SFR < 10 M⊙ yr−1). The difference is mainly in a steeper spectrum around 10 GHz, which could indicate a smaller fraction of thermal free–free emission. Finally, we also discuss the impact of applying this broken power-law SED in place of a simple power-law in K-corrections of HSFGs and a typical radio SED for normal SFGs drawn from the literature. We find that the shape of the radio SED is unlikely to be the root cause of the q − z trend in SFGs.
Context.
As the Square Kilometer Array (SKA) is expected to be operational in the next decade, investigations of the radio sky in the range of 100 MHz–10 GHz have become important for simulating SKA ...observations. In determining physical properties of galaxies from radio data, the radio spectral energy distribution (SED) is often assumed to be described by a simple power law, usually with a spectral index of 0.7 for all sources. Even though radio SEDs have been shown to exhibit deviations from this assumption, both in differing spectral indices and complex spectral shapes, it is often presumed that their individual differences can be canceled out in large samples.
Aims.
Since the average spectral index around 1 GHz (observed-frame) is important for determining physical properties of large samples of galaxies, we aim to test whether individual differences in the spectra of radio-identified active galactic nuclei align with the simple assumption of
α
= 0.7 and test the evolution of the parameters of the synchrotron aging model with redshift and radio luminosity.
Methods.
We use a sample of 744 radio-excess active galactic nuclei (RxAGN), defined as those that exhibit more than a 3
σ
radio luminosity excess with respect to the value expected only from the contribution from star formation, out to
z
∼ 4. We constructed their average radio SED by combining Very Large Array (VLA) observations of the COSMOS field at 1.4 GHz and 3 GHz with Giant Meterwave Radio Telescope (GMRT) observations at 325 MHz and 610 MHz. To account for nondetections in the GMRT maps, we employed the survival analysis technique. We binned the RxAGN sample into luminosity- and redshift-complete subsamples. In each bin, we constrained the shape of the average radio SED by fitting a broken power-law model.
Results.
We find that the RxAGN sample can be described by a spectral index of
α
1
= 0.28 ± 0.03 below the break frequency
ν
b
= (4.1 ± 0.2) GHz and
α
2
= 1.16 ± 0.04 above it, while a simple power-law model, capturing fewer spectral features, yields a single spectral index of 0.64 ± 0.07. By binning in 1.4 GHz of radio luminosity and redshift, we find that the power-law spectral index is positively correlated with redshift and that the broken power-law spectral index above 4 GHz is positively correlated with both the redshift and source size. By selecting sources with sizes less than 1 kpc, we find a subsample of flat-spectrum sources, which can be described by a spectral index of
α
= 0.41 ± 0.07 and a broken power-law spectral index of
α
1
= 0.1 ± 0.1 (
α
2
= 0.55 ± 0.09) below (above) a break frequency of
ν
b
= (2.7 ± 0.5) GHz.
Conclusions.
We have constrained the radio SED for a sample of RxAGN in the COSMOS field using available VLA and GMRT data, corresponding to the rest-frame frequency range from ∼0.3 GHz to ∼10 GHz. We describe our derived average radio SED of RxAGN using power-law and broken power-law models, yielding a radio SED that steepens above ∼4 GHz.