Based on a sample of over 1, 800 radio AGN at redshifts out to z ∼ 5, which have typical stellar masses within ∼ 3 × (10 10 − 10 11 ) M ⊙ , and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L 1.4GHz ∼ 10 22 − 10 27 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.

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 /N 3 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 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 highredshift 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 provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) -stellar mass (M * ) relation, called the main sequence of star-forming galaxies (MS), for starforming and all galaxies out to z ∼ 5. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for ∼ 200,000 mass-selected galaxies at z > 0.3 in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log(M * /M )<10) and a flattening at high stellar mass that becomes more prominent at low redshift (z < 1.5). We find that the SFR density peaks at 1.5 < z < 2 and at each epoch there is a characteristic stellar mass (M * = 1−4×10 10 M ) that contributes the most to the overall SFR density. This characteristic mass increases with redshift, at least to z ∼ 2.5. We find no significant evidence for variations in the MS relation for galaxies in different environments traced by the galaxy number density at 0.3 < z < 3, nor for galaxies in X-ray groups at z ∼ 0.75. We confirm that massive bulgedominated galaxies have lower SFRs than disk-dominated galaxies at a fixed stellar mass at z < 1.2. As a consequence, the increase in bulge-dominated galaxies in the local star-forming population leads to a flattening of the MS at high stellar masses. This indicates that "mass-quenching" is linked with changes in the morphological composition of galaxies at a fixed stellar mass.

We present the results of a search for galaxy clusters and groups in the ∼ 2 square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of 3 × 10 −16 ergs cm −2 s −1 in 0.5-2 keV range, and identify 247 X-ray groups with M 200c = 8 × 10 12 − 3 × 10 14 M at a redshift range of 0.08 z < 1.53, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At z < 0.5, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both lowand large-offset BGGs at fixed halo mass is ∼ 0.3dex.