2025MNRAS.539.1856H

The wide-area component of the LOFAR Two-Metre Sky Survey (LoTSS) is currently the largest radio survey ever carried out, and a large fraction of the 4.5 million radio sources it contains have been optically identified with galaxies or quasars with spectroscopic or photometric redshifts. Identification of radio-luminous active galactic nucleus (AGN) from this LoTSS source catalogue is not only important from the point of view of understanding the accretion history of the universe, but also enables a wide range of other science. However, at present the vast majority of the optical identifications lack spectroscopic information or well-sampled spectral energy distributions. We show that colour and absolute magnitude information from the Wide-Field Infrared Survey Explorer (WISE) allows for the robust and efficient selection of radio AGN candidates, generating a radio AGN candidate sample of around 600 000 objects with flux density $> 1.1$ mJy, spanning 144-MHz luminosities between $10^{21}$ and $10^{29}$ W Hz$^{-1}$. We use the catalogue to constrain the total sky density of radio-luminous AGN and the evolution of their luminosity function between $z=0$ and $z\approx 1$, and show that the typical mass of their host galaxies, around $10^{11} {\rm M}_\odot$, is essentially independent of radio luminosity above around $L_{144} \approx 10^{24}$ W Hz$^{-1}$. Combining with Very Large Array Sky Survey (VLASS) data, we show that the core prominences, radio spectral indices and variability of extended sources from the sample are qualitatively consistent with the expectations from unified models. A catalogue of the radio AGN candidates is released with this paper.