Abstract:We model the evolution of active galactic nuclei (AGN) by constructing their radio LFs. We used a set of surveys of varying area and depth, namely, the deep COSMOS survey of $1,916$ AGN sources; the wide, shallow 3CRR, 7C, and 6CE surveys, together containing $356$ AGN; and the intermediate XXL-North and South fields consisting of $899$ and $1,484$ sources, respectively. We also used the CENSORS, BRL, Wall $ $ Peacock, and Config surveys, respectively consisting of $150$, $178$, $233$, and $230$ sources. Together, these surveys account for $5,446$ AGN sources and constrained the LFs at high redshift and over a wide range of luminosities (up to $z 3$ and $ (L / W Hz^ )$. We concentrated on parametric methods within the Bayesian framework, which allowed us to perform model selection between a set of different models. By comparing the marginalised likelihoods and both the Akaike information criterion and the Bayesian information criterion, we show that the luminosity-dependent density evolution (LDDE) model fits the data best, with evidence ratios varying from "strong" ($>10$) to "decisive" ($>100$), according to the Jeffreys' interpretation. The median values of the model parameters, defined via $ (1 + z_c)^ p_1 +(1 + z_c)^ p_2 1+z_c 1+z p_1 1+z_c 1+z p_2 with $z_c = z_c^*$ for $L > L_a$ and $z_c^* L L_a for $L < L_a$ equalled $z^* = 2.01$, $L_a = 27.94$, $a = 0.40$, $p_1 = -0.39$, $p_2 = 4.53$, $ = -3.92$, $L^* = 22.26$, $ = 1.39$ and $ The best-fitting model gives insight into the physical picture of AGN evolution, where AGN evolve differently as a function of their radio luminosity. We determined the number density, luminosity density, and kinetic luminosity density as a function of redshift, and we observed a flattening of these functions at higher redshifts, which is not present in simpler models. We explain these trends by our use of the LDDE model. Finally, we divided our sample into subsets according to the stellar mass of the host galaxies in order to investigate a possible bimodality in evolution. We found a difference in LF shape and evolution between these subsets. All together, these findings point to a physical picture where the evolution and density of AGN cannot be explained well by simple models but require more complex models either via AGN sub-populations, where the total AGN sample is divided into sub-samples according to various properties, such as optical properties and stellar mass, or via luminosity-dependent functions.

I'm in AGNi: A new standard for AGN pluralisation

Active galactic nuclei: what’s in a name?

pAGN : the one-stop solution for AGN disc modeling

pAGN: the one-stop solution for AGN disc modeling

A Brief History of AGN

Amoeba: An AGN Model of Optical Emissions Beyond steady-state Accretion discs

Recommendations from the Ad Hoc Committee on SETI Nomenclature

The global properties of all variety of AGN

The XXL survey. LII. The evolution of radio AGN LF determined viaparametric methods from GMRT, ATCA, VLA, and Cambridge interferometer observations

An active galactic nucleus recognition model based on deep neural network

Nomenclature of The Galactic Center Radio Sources

Super-Size Me: The Big Multi-AGN Catalog (The Big MAC), Data Release 1: The Source Catalog

A Walk-Through of AGN Country -- for the somewhat initiated!

AGN populations in the local Universe: Their alignment with the main sequence, stellar population characteristics, accretion efficiency, and the impact of AGN feedback

Active galactic nuclei with high-resolution X-ray spectroscopy

Towards a Complete Census of AGNs in Nearby Galaxies: A Large Population of Optically Unidentified AGNs

Aegis: Demographics of X-Ray and Optically Selected Active Galactic Nuclei

Changing-look Active Galactic Nuclei from the Dark Energy Spectroscopic Instrument. I.Sample from the Early Data

The innermost regions of Active Galactic Nuclei - from radio to X-rays

AGN Selection and Demographics: A New Age with JWST/MIRI