The VLA-COSMOS 3 GHz Large Project: Average radio spectral energy distribution of active galactic nuclei

As the SKA is expected to be operational in the next decade, investigations of the radio sky in the range of 100 MHz to 10 GHz have become important for simulations of the SKA observations. In determining physical properties of galaxies from radio data, the radio 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 cancel out in large samples. We constructed the average radio SED of radio-excess active galactic nuclei (RxAGN), defined as those that exhibit a 3 $sigma$ radio luminosity excess with respect to the value expected only from contribution from star formation, out to z~4. We combined VLA observations of the COSMOS field at 1.4 GHz and 3 GHz with GMRT observations at 325 MHz and 610 MHz. To account for nondetections in the GMRT maps, we employed the survival analysis technique. We selected a sample of RxAGN out to z~4. We find that a sample of RxAGN can be described by a spectral index of $alpha_1=0.28pm0.03$ below the break frequency $ u_b=(4.1pm0.2)$ GHz and $alpha_2=1.16pm0.04$ above, while a simple power-law model yields a single spectral index of $alpha=0.64pm0.07$. By binning in 1.4 GHz radio luminosity and redshift, we find that the power-law spectral index, as well as broken power-law spectral indices, may increase for larger source sizes, while the power-law spectral index and lower-frequency (<4 GHz) broken power-law spectral index are additionally positively correlated with redshift.