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.
We present very long baseline interferometry (VLBI) observations of 179 radio sources in the COSMOS field with extremely high sensitivity using the Green Bank Telescope (GBT) together with the Very ...Long Baseline Array (VLBA) (VLBA+GBT) at 1.4 GHz, to explore the faint radio population in the flux density regime of tens of μJy. Here, the identification of active galactic nuclei (AGN) is based on the VLBI detection of the source, meaning that it is independent of X-ray or infrared properties. The milli-arcsecond resolution provided by the VLBI technique implies that the detected sources must be compact and have large brightness temperatures, and therefore they are most likely AGN (when the host galaxy is located at z ≥ 0.1). On the other hand, this technique only allows us to positively identify when a radio-active AGN is present, in other words, we cannot affirm that there is no AGN when the source is not detected. For this reason, the number of identified AGN using VLBI should be always treated as a lower limit. We present a catalogue containing the 35 radio sources detected with the VLBA+GBT, ten of which were not previously detected using only the VLBA. We have constructed the radio source counts at 1.4 GHz using the samples of the VLBA and VLBA+GBT detected sources of the COSMOS field to determine a lower limit for the AGN contribution to the faint radio source population. We found an AGN contribution of >40−75% at flux density levels between 150 μJy and 1 mJy. This flux density range is characterised by the upturn of the Euclidean-normalised radio source counts, which implies a contribution of a new population. This result supports the idea that the sub-mJy radio population is composed of a significant fraction of radio-emitting AGN, rather than solely by star-forming galaxies, in agreement with previous studies.
We explore the redshift evolution of the specific star formation rate (SSFR) for galaxies of different stellar mass by drawing on a deep 3.6 Delta *mm selected sample of >105 galaxies in the 2 deg2 ...COSMOS field. The average star formation rate (SFR) for subsets of these galaxies is estimated with stacked 1.4 GHz radio continuum emission. We separately consider the total sample and a subset of galaxies that shows evidence for substantive recent star formation in the rest-frame optical spectral energy distributions. At redshifts 0.2 < z < 3 both populations show a strong and mass-independent decrease in their SSFR toward the present epoch. It is best described by a power law (1 + z) n , where n ~ 4.3 for all galaxies and n ~ 3.5 for star-forming (SF) sources. The decrease appears to have started at z>2, at least for high-mass (M * 4 X 1010 M ) systems where our conclusions are most robust. Our data show that there is a tight correlation with power-law dependence, SSFR M * Delta *b, between SSFR and stellar mass at all epochs. The relation tends to flatten below M * 1010 M if quiescent galaxies are included; if they are excluded from the analysis a shallow index Delta *bSFG --0.4 fits the correlation. On average, higher mass objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper threshold in SSFR that an average galaxy cannot exceed, possibly due to gravitationally limited molecular gas accretion. It is suggested by a flattening of the SSFR-M * relation (also for SF sources), but affects massive (>1010 M ) galaxies only at the highest redshifts. Since z = 1.5 there thus is no direct evidence that galaxies of higher mass experience a more rapid waning of their SSFR than lower mass SF systems. In this sense, the data rule out any strong 'downsizing' in the SSFR. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of an SF galaxy: we find that since z ~ 3 the majority of all new stars were always formed in galaxies of M * = 1010.6?0.4 M . In this sense, too, there is no 'downsizing.' Finally, our analysis constitutes the most extensive SFR density determination with a single technique out to z = 3. Recent Herschel results are consistent with our results, but rely on far smaller samples.
We present a sub-kiloparsec localization of the sites of supermassive black hole (SMBH) growth in three active galactic nuclei (AGNs) at z ∼ 3 in relation to the regions of intense star formation in ...their hosts. These AGNs are selected from Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA) observations in the Hubble Ultra-Deep Field and COSMOS, with the centimetric radio emission tracing both star formation and AGN, and the sub/millimeter emission by dust tracing nearly pure star formation. We require radio emission to be more luminous than the level associated with the sub/millimeter star formation to ensure that the radio emission is AGN-dominated, thereby allowing localization of the AGN and star formation independently. In all three galaxies, the AGNs are located within the compact regions of gas-rich, heavily obscured, intense nuclear star formation, with Re = 0.4-1.1 kpc and average star formation rates of 100-1200 M yr−1. If the current episode of star formation continues at such a rate over the stellar mass doubling time of their hosts, 0.2 Gyr, the newly formed stellar mass will be of the order of 1011 M within the central kiloparsec region, concurrently and cospatially with significant growth of the SMBH. This is consistent with a picture of in situ galactic bulge and SMBH formation. This work demonstrates the unique complementarity of VLA and ALMA observations to unambiguously pinpoint the locations of AGNs and star formation down to 30 mas, corresponding to 230 pc at z = 3.
ABSTRACT
Ly α λ1216 (Ly α) emission extending over $\gtrsim \, \rm 10\, kilo\,parsec\, (kpc)$ around dusty, massive starbursts at z ≳ 3 might represent a short-lived phase in the evolution of ...present-day, massive quiescent galaxies. To obtain empirical constraints on this emerging scenario, we present Ly α, C ivλ1550 (C iv), and He ii λ1640 (He ii) observations taken with the Multi-Unit Spectroscopic Explorer towards J1000+0234: a galaxy pair at z = 4.5 composed of a low-mass starburst (J1000+0234−South) neighbouring a massive Submillimeter Galaxy (SMG; J1000+0234−North) that harbours a rotationally supported gas disc. Based on the spatial distribution and relative strength of Ly α, C iv, and He ii, we find that star formation in J1000+0234−South and an active galactic nucleus in J1000+0234−North are dominant factors in driving the observed 40 kiloparsec-scale Ly α blob (LAB). We use the non-resonant He ii line to infer kinematic information of the LAB. We find marginal evidence for two spatially and spectrally separated He ii regions, which suggests that the two-peaked Ly α profile is mainly a result of two overlapping and likely interacting H i clouds. We also report the serendipitous identification of three Ly α emitters spanning over a redshift bin Δz ≤ 0.007 (i.e. $\lesssim 380\, \rm km\, s^{-1}$) located at $\lesssim 140\, \rm kpc$ from J1000+0234. A galaxy overdensity analysis confirms that J1000+0234 lies near the centre of a Megaparsec-scale galaxy overdensity at z = 4.5 that might evolve into a galaxy cluster at z = 0. The properties of J1000+0234 and its large-scale environment strengthen the link between SMGs within LABs, tracing overdense regions, as the progenitors of local massive ellipticals in galaxy clusters.
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.
Future radio surveys will generate catalogs of tens of millions of radio sources, for which redshift estimates will be essential to achieve many of the science goals. However, spectroscopic data will ...be available for only a small fraction of these sources, and in most cases even the optical and infrared photometry will be of limited quality. Furthermore, radio sources tend to be at higher redshift than most optical sources (most radio surveys have a median redshift greater than 1) and so a significant fraction of radio sources hosts differ from those for which most photometric redshift templates are designed. We therefore need to develop new techniques for estimating the redshifts of radio sources. As a starting point in this process, we evaluate a number of machine-learning techniques for estimating redshift, together with a conventional template-fitting technique. We pay special attention to how the performance is affected by the incompleteness of the training sample and by sparseness of the parameter space or by limited availability of ancillary multiwavelength data. As expected, we find that the quality of the photometric-redshift degrades as the quality of the photometry decreases, but that even with the limited quality of photometry available for all-sky-surveys, useful redshift information is available for the majority of sources, particularly at low redshift. We find that a template-fitting technique performs best in the presence of high-quality and almost complete multi-band photometry, especially if radio sources that are also X-ray emitting are treated separately, using specific templates and priors. When we reduced the quality of photometry to match that available for the EMU all-sky radio survey, the quality of the template-fitting degraded and became comparable to some of the machine-learning methods. Machine learning techniques currently perform better at low redshift than at high redshift, because of incompleteness of the currently available training data at high redshifts.
Context. Quantifying the fraction of active galactic nuclei (AGN) in the faint radio population and understanding their relation with star-forming activity are fundamental to studies of galaxy ...evolution. Very long baseline interferometry (VLBI) observations are able to identify AGN above relatively low redshifts (z> 0.1) since they provide milli-arcsecond resolution. Aims. We have created an AGN catalogue from 2865 known radio sources observed in the Cosmic Evolution Survey (COSMOS) field, which has exceptional multi-wavelength coverage. With this catalogue we intend to study the faint radio sky with statistically relevant numbers and to analyse the AGN – host galaxy co-evolution, making use of the large amount of ancillary data available in the field. Methods. Wide-field VLBI observations were made of all known radio sources in the COSMOS field at 1.4 GHz to measure the AGN fraction, in particular in the faint radio population. We describe in detail the observations, data calibration, source detection and flux density measurements, parts of which we have developed for this survey. The combination of number of sources, sensitivity, and area covered with this project are unprecedented. Results. We have detected 468 radio sources, expected to be AGN, with the Very Long Baseline Array (VLBA). This is, to date, the largest sample assembled of VLBI detected sources in the sub-mJy regime. The input sample was taken from previous observations with the Very Large Array (VLA). We present the catalogue with additional optical, infrared and X-ray information. Conclusions. We find a detection fraction of 20 ± 1%, considering only those sources from the input catalogue which were in principle detectable with the VLBA (2361). As a function of the VLA flux density, the detection fraction is higher for higher flux densities, since at high flux densities a source could be detected even if the VLBI core accounts for a small percentage of the total flux density. As a function of redshift, we see no evolution of the detection fraction over the redshift range 0.5 <z< 3. In addition, we find that faint radio sources typically have a greater fraction of their radio luminosity in a compact core – ~70% of the sub-mJy sources detected with the VLBA have more than half of their total radio luminosity in a VLBI-scale component, whereas this is true for only ~30% of the sources that are brighter than 10 mJy. This suggests that fainter radio sources differ intrinsically from brighter ones. Across our entire sample, we find the predominant morphological classification of the host galaxies of the VLBA detected sources to be early type (57%), although this varies with redshift and at z> 1.5 we find that spiral galaxies become the most prevalent (48%). The number of detections is high enough to study the faint radio population with statistically significant numbers. We demonstrate that wide-field VLBI observations, together with new calibration methods such as multi-source self-calibration and mosaicing, result in information which is difficult or impossible to obtain otherwise.
Abstract
Dark matter haloes in which galaxies reside are likely to have a significant impact on their evolution. We investigate the link between dark matter haloes and their constituent galaxies by ...measuring the angular two-point correlation function of radio sources, using recently released 3 GHz imaging over ∼2 deg2 of the Cosmological Evolution Survey (COSMOS) field. We split the radio source population into star-forming galaxies (SFGs) and active galactic nuclei (AGN), and further separate the AGN into radiatively efficient and inefficient accreters. Restricting our analysis to z < 1, we find SFGs have a bias, $b = 1.5 ^{+0.1}_{-0.2}$, at a median redshift of z = 0.62. On the other hand, AGN are significantly more strongly clustered with b = 2.1 ± 0.2 at a median redshift of 0.7. This supports the idea that AGN are hosted by more massive haloes than SFGs. We also find low accretion rate AGN are more clustered (b = 2.9 ± 0.3) than high accretion rate AGN ($b = 1.8^{+0.4}_{-0.5}$) at the same redshift (z ∼ 0.7), suggesting that low accretion rate AGN reside in higher mass haloes. This supports previous evidence that the relatively hot gas that inhabits the most massive haloes is unable to be easily accreted by the central AGN, causing them to be inefficient. We also find evidence that low accretion rate AGN appear to reside in halo masses of Mh ∼ 3–4 × 1013 h−1 M⊙ at all redshifts. On the other hand, the efficient accreters reside in haloes of Mh ∼ 1–2 × 1013 h−1 M⊙ at low redshift but can reside in relatively lower mass haloes at higher redshifts. This could be due to the increased prevalence of cold gas in lower mass haloes at z ≥ 1 compared to z < 1.
Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio active galactic nuclei (AGNs) to be used to ...their fullest potential. We study the composition of the Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of ∼350 arcmin2. This ultra-deep survey consists of a single pointing in the well-studied COSMOS field at both 3 and 10 GHz and reaches rms sensitivities of 0.53 and 0.41 Jy beam−1, respectively. We find multiwavelength counterparts for 97% of radio sources, based on a combination of near-UV/optical to sub-millimeter data, and through a stacking analysis at optical/near-IR wavelengths we further show that the sources lacking such counterparts are likely to be high-redshift in nature (typical z ∼ 4−5). Utilizing the multiwavelength data over COSMOS, we identify AGNs through a variety of diagnostics and find these to make up 23.2 1.3% of our sample, with the remainder constituting uncontaminated star-forming galaxies. However, more than half of the AGNs exhibit radio emission consistent with originating from star formation, with only 8.8 0.8% of radio sources showing a clear excess in radio luminosity. At flux densities of ∼30 Jy at 3 GHz, the fraction of star formation-powered sources reaches ∼90%, and this fraction is consistent with unity at even lower flux densities. Overall, our findings imply that ultra-deep radio surveys such as COSMOS-XS constitute a highly effective means of obtaining clean samples of star formation-powered radio sources.