Search

We present the first results of the JWST Emission Line Survey (JELS). Utilising the first NIRCam narrow-band imaging at 4.7$\mu$m, over 63 arcmin$^{2}$ in the PRIMER/COSMOS field, we identified 609 emission line galaxy candidates. From these, we robustly selected 35 H$\alpha$ star-forming galaxies at $z \sim 6.1$, with H$\alpha$ star-formation rates ($\rm{SFR_{H\alpha}}$) $\sim0.9-15\ \rm{M_{\odot} \ yr^{-1}}$. Combining our unique H$\alpha$ sample with the exquisite panchromatic data in the field, we explored their physical properties and star-formation histories, and compared these to a broad-band selected sample at $z\sim 6$ which offered vital new insights into the nature of high-redshift galaxies. UV-continuum slopes ($\beta$) were considerably redder for our H$\alpha$ sample ($\langle\beta\rangle\sim-1.92$) compared to the broad-band sample ($\langle\beta\rangle\sim-2.35$). This was not due to dust attenuation as our H$\alpha$ sample was relatively dust-poor (median $A_V=0.23$); instead, we argued the reddened slopes could be due to nebular continuum. We compared $\rm{SFR_{H\alpha}}$ and the UV-continuum-derived $\rm{SFR_{UV}}$ to SED-fitted measurements averaged over canonical timescales of 10 and 100 Myr ($\rm{SFR_{10}}$ and $\rm{SFR_{100}}$). We found an increase in recent SFR for our sample of H$\alpha$ emitters, particularly at lower stellar masses ($<10^9 \ \rm{M_{\odot}}$). We also found $\rm{SFR_{H\alpha}}$ strongly traced SFR averaged over 10 Myr timescales, whereas the UV-continuum over-predicts SFR on 100 Myr timescales at low stellar masses. These results point to our H$\alpha$ sample undergoing `bursty' star formation. Our F356W $z \sim 6$ sample showed a larger scatter in $\rm{SFR_{10}/SFR_{100}}$ across all stellar masses, which highlighted how narrow-band photometric selections of H$\alpha$ emitters are key to quantifying the burstiness of star-formation activity., Comment: Submitted to MRNAS; 27 pages, 19 figures. Accompanying data will be publicly released upon acceptance of this manuscript

We present the JWST Emission Line Survey (JELS), a Cycle 1 JWST imaging programme exploiting the wavelength coverage and sensitivity of NIRCam to extend narrow-band rest-optical emission line selection into the epoch of reionization (EoR) for the first time, and to enable unique studies of the resolved ionised gas morphology in individual galaxies across cosmic history. The primary JELS observations comprise $\sim4.7\mu$m narrow-band imaging over $\sim63$ arcmin$^{2}$ designed to enable selection of H$\alpha$ emitters at $z\sim6.1$, as well as the selection of a host of novel emission-line samples including [OIII] at $z\sim8.3$ and Pa $\alpha/\beta$ at $z\sim1.5/2.8$. For the prime F466N and F470N narrow-band observations, the emission-line sensitivities achieved are up to $\sim2\times$ more sensitive than current slitless spectroscopy surveys (5$\sigma$ limits of 1.1-1.6$\times10^{-18}\text{erg s}^{-1}\text{cm}^{-2}$), corresponding to unobscured H$\alpha$ star-formation rates (SFRs) of 1-1.6 $\text{M}_{\odot}\,\text{yr}^{-1}$ at $z\sim6.1$ and extending emission-line selections in the EoR to fainter populations. Simultaneously, JELS also obtained F200W broadband and F212N narrow-band imaging (H$\alpha$ at $z\sim2.23$) that probes SFRs $\gtrsim5\times$ fainter than previous ground-based narrow-band studies ($\sim0.2 \text{M}_{\odot}\text{yr}^{-1}$), offering an unprecedented resolved view of star formation at cosmic noon. In this paper we describe the detailed JELS survey design, key data processing steps specific to the survey observations, and demonstrate the exceptional data quality and imaging sensitivity achieved. We then summarise the key scientific goals of JELS and present some early science results, including examples of spectroscopically confirmed H$\alpha$ and [OIII] emitters discovered by JELS that illustrate the novel parameter space probed., Comment: Submitted to MRNAS; 13 pages, 8 figures. Accompanying data will be publicly released upon acceptance of this and companion paper (Pirie et al., subm.)

We present an analysis of 152,355 radio sources identified in the second data release of the LOFAR Two Metre Sky Survey (LoTSS-DR2) with Sloan Digital Sky Survey (SDSS) spectroscopic redshifts in the range 0.00 < z < 0.57. Using Monte Carlo simulations we determine the reliability of each source exhibiting an excess in radio luminosity relative to that predicted from their Ha emission, and, for a subset of 124,023 sources we combine this measurement with a full BPT analysis. Using these two independent diagnostics we determine the reliability of each source hosting a supermassive black hole of high or low Eddington-scaled accretion rate, and combine the measurements to determine the reliability of sources belonging to each of four physical classes of objects: star forming galaxies (SFGs), radio-quiet active galactic nuclei (RQAGN), and high- or low-excitation radio galaxies (HERGs or emission-line LERGs). The result is a catalogue which enables user-defined samples of radio sources with a reliability threshold suited to their science goal e.g. prioritising purity or completeness. Here we select high-confidence samples of radio sources (>90% reliability) to report: 38,588 radio-excess AGN in the LoTSS DR2 sample (362 HERGs, and 12,648 emission-line LERGs), together with 38,729 SFGs, and 18,726 RQAGN. We validate these results through comparison to literature using independent emission-line measurements, and to widely-adopted WISE photometric selection techniques. While our use of SDSS spectroscopy limits our current analysis to ~4 percent of the LoTSS-DR2 catalogue, our method is directly applicable to data from the forthcoming WEAVE-LOFAR survey which will obtain over a million spectra of 144 MHz selected sources., Comment: Accepted for publication in MNRAS

Studies show that both radio jets from the active galactic nuclei (AGN) and the star formation (SF) activity in quasar host galaxies contribute to the quasar radio emission; yet their relative contributions across the population remain unclear. Here, we present an improved parametric model that allows us to statistically separate the SF and AGN components in observed quasar radio flux density distributions, and investigate how their relative contributions evolve with AGN bolometric luminosity ($L_\mathrm{bol}$) and redshift ($z$) using a fully Bayesian method. Based on the newest data from LOFAR Two-Metre Sky Survey Data Release 2, our model gives robust fitting results out to $z\sim4$, showing a quasar host galaxy SFR evolution that increases with bolometric luminosity and with redshift out to $z\sim4$. This differs from the global cosmic SFR density, perhaps due to the importance of galaxy mergers. The prevalence of radio AGN emissions increases with quasar luminosity, but has little dependence on redshift. Furthermore, our new methodology and large sample size allow us to subdivide our dataset to investigate the role of other parameters. Specifically, in this paper, we explore quasar colour and demonstrate that the radio excess in red quasars is due to an enhancement in AGN-related emission, since the host galaxy SF contribution to the total radio emission is independent of quasar colour. We also find evidence that this radio enhancement occurs mostly in quasars with weak or intermediate radio power., Comment: Accepted for publication in MNRAS. 19 pages, 13 figures, 1 table

The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3\sigma$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths., Comment: 32 pages, 11 figures, 5 tables, 1 Appendix

We present an analysis of the quenching of star formation in massive galaxies ($M_* > 10^{9.5} M_\odot$) within the first 0.5 - 3 Gyr of the Universe's history utilizing JWST-CEERS data. We utilize a combination of advanced statistical methods to accurately constrain the intrinsic dependence of quenching in a multi-dimensional and inter-correlated parameter space. Specifically, we apply Random Forest (RF) classification, area statistics, and a partial correlation analysis to the JWST-CEERS data. First, we identify the key testable predictions from two state-of-the-art cosmological simulations (IllustrisTNG & EAGLE). Both simulations predict that quenching should be regulated by supermassive black hole mass in the early Universe. Furthermore, both simulations identify the stellar potential ($\phi_*$) as the optimal proxy for black hole mass in photometric data. In photometric observations, where we have no direct constraints on black hole masses, we find that the stellar potential is the most predictive parameter of massive galaxy quenching at all epochs from $z = 0 - 8$, exactly as predicted by simulations for this sample. The stellar potential outperforms stellar mass, galaxy size, galaxy density, and S\'ersic index as a predictor of quiescence at all epochs probed in JWST-CEERS. Collectively, these results strongly imply a stable quenching mechanism operating throughout cosmic history, which is closely connected to the central gravitational potential in galaxies. This connection is explained in cosmological models via massive black holes forming and growing in deep potential wells, and subsequently quenching galaxies through a mix of ejective and preventative active galactic nucleus (AGN) feedback., Comment: Accepted to ApJ. 32 pages; 10 figures

Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $\mu$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies within the survey. We discuss systematics in the data and an improved methodology for generating random catalogues, compared to that used for LoTSS-DR1, before presenting the angular clustering for $\sim$900,000 sources $\geq$$1.5$ mJy and a peak signal-to-noise $\geq$$7.5$ across $\sim$$80\%$ of the observed area. Using the clustering we infer the bias assuming two evolutionary models. When fitting {angular scales of $0.5 \leq\theta<5\,\deg$, using a linear bias model, we find LoTSS-DR2 sources are biased tracers of the underlying matter, with a bias of $b_{C}= 2.14^{+0.22}_{-0.20}$ (assuming constant bias) and $b_{E}(z=0)= 1.79^{+0.15}_{-0.14}$ (for an evolving model, inversely proportional to the growth factor), corresponding to $b_E= 2.81^{+0.24}_{-0.22}$ at the median redshift of our sample, assuming the LoTSS Deep Fields redshift distribution is representative of our data. This reduces to $b_{C}= 2.02^{+0.17}_{-0.16}$ and $b_{E}(z=0)= 1.67^{+0.12}_{-0.12}$ when allowing preferential redshift distributions from the Deep Fields to model our data. Whilst the clustering amplitude is slightly lower than LoTSS-DR1 ($\geq$2 mJy), our study benefits from larger samples and improved redshift estimates., Comment: Accepted for publication in MNRAS. 29 pages, 24 figures

We present the low-frequency radio spectra of 9 high-redshift quasars at $5.6 \leq z \leq 6.6$ using the Giant Metre Radio Telescope band-3, -4, and -5 observations ($\sim$300-1200 MHz), archival Low Frequency Array (LOFAR; 144 MHz), and Very Large Array (VLA; 1.4 and 3 GHz) data. Five of the quasars in our sample have been discovered recently, representing some of the highest redshift radio bright quasars known at low-frequencies. We model their radio spectra to study their radio emission mechanism and age of the radio jets by constraining the spectral turnover caused by synchrotron self-absorption (SSA) or free-free absorption (FFA). Besides J0309+2717, a blazar at $z=6.1$, our quasars show no sign of a spectral flattening between 144 MHz and a few GHz, indicating there is no strong SSA or FFA absorption in the observed frequency range. However, we find a wide range of spectral indices between $-1.6$ and $0.05$, including the discovery of 3 potential ultra-steep spectrum quasars. Using further archival VLBA data, we confirm that the radio SED of the blazar J0309+2717 likely turns over at a rest-frame frequency of 0.6-2.3 GHz (90-330 MHz observed frame), with a high-frequency break indicative of radiative ageing of the electron population in the radio lobes. Ultra-low frequency data below 50 MHz are necessary to constrain the absorption mechanism for J0309+2717 and the turnover frequencies for the other high-$z$ quasars in our sample. A relation between linear radio jet size and turnover frequency has been established at low redshifts. If this relation were to hold at high redshifts, the limits on the turnover frequency of our sample suggest the radio jet sizes must be more extended than the typical sizes observed in other radio-bright quasars at similar redshift. To confirm this deep radio follow-up observations with high spatial resolution are required., Comment: 15 pages, 7 figures, accepted for publication in A&A, minor correction Tab. 1

We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities ($\sim10^{17-18}$ erg s$^{-1}$ Hz$^{-1}$) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant ($>$ 50 pc) radio emitting M dwarfs., Comment: 10 pages, 5 figures, accepted for publication in A&A

The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many extended radio sources found in LOFAR images as a result of its excellent sensitivity to extended structure. In this paper we present source associations and identifications for sources in the second data release based on optical and near-infrared data, using a combination of a likelihood-ratio cross-match method developed for our first data release, our citizen science project Radio Galaxy Zoo: LOFAR, and new approaches to algorithmic optical identification, together with extensive visual inspection by astronomers. We also present spectroscopic or photometric redshifts for a large fraction of the optical identifications. In total 4,116,934 radio sources lie in the area with good optical data, of which 85% have an optical or infrared identification and 58% have a good redshift estimate. We demonstrate the quality of the dataset by comparing it with earlier optically identified radio surveys. This is by far the largest ever optically identified radio catalogue, and will permit robust statistical studies of star-forming and radio-loud active galaxies., Comment: 29 pages. Accepted by A&A; data products available at https://lofar-surveys.org/dr2_release.html

We present a detailed study of the cosmic star formation history over $90$ per cent of cosmic time ($0\lesssim z\lesssim4$), using deep, radio continuum observations that probe star formation activity independent of dust. The Low Frequency Array Two Metre Sky Survey has imaged three well-studied extragalactic fields, Elais-N1, Bo\"otes and the Lockman Hole, reaching $\sim20\,\mu\rm{Jy/beam}$ rms sensitivity at $150\,\rm{MHz}$. The availability of high-quality ancillary data from ultraviolet to far-infrared wavelengths has enabled accurate photometric redshifts and the robust separation of radio-bright AGN from their star-forming counterparts. We capitalise on this unique combination of deep, wide fields and robustly-selected star-forming galaxies to construct radio luminosity functions and derive the cosmic star formation rate density. We carefully constrain and correct for scatter in the $L_{150\,\rm{MHz}}-\rm{SFR}$ relation, which we find to be $\sim0.3\,\rm{dex}$. Our derived star formation rate density lies between previous measurements at all redshifts studied. We derive higher star formation rate densities between $z\sim0$ and $z\sim3$ than are typically inferred from short wavelength emission; at earlier times, this discrepancy is reduced. Our measurements are generally in good agreement with far-infrared and radio-based studies, with small offsets resulting from differing star formation rate calibrations., Comment: 22 pages, 14 figures; accepted for publication in MNRAS

Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested., Comment: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

We present a sample of 373 peaked-spectrum (PS) sources with spectral peaks around 150MHz, selected using a subset of two LOFAR all-sky surveys, the LOFAR Two Meter Sky Survey and the LOFAR LBA Sky Survey. These surveys are the most sensitive low-frequency widefield surveys to date, allowing us to select low-luminosity PS sources. Our sample increases the number of known PS sources in our survey area by a factor 50. The 5GHz luminosity distribution of our PS sample shows we sample the lowest luminosity PS sources to-date by nearly an order of magnitude. Since high-frequency PS sources and compact steep-spectrum sources are hypothesised to be the precursors to large radio galaxies, we investigate whether this is also the case for our sample of low-frequency PS sources. Using optical line emission criteria, we find that our PS sources are predominately high-excitation radio galaxies instead of low-excitation radio galaxies, corresponding to a quickly evolving population. We compute the radio source counts of our PS sample, and find they are scaled down by a factor of $\sim$40 compared to a general sample of radio-loud active galactic nuclei (AGN). This implies that the lifetimes of PS sources are 40 times shorter than large scale radio galaxies, if their luminosity functions are identical. To investigate this, we compute the first radio luminosity function for a homogeneously-selected PS sample. We find that for 144MHz luminosities $\gtrsim 10^{25}$W Hz$^{-1}$, the PS luminosity function has the same shape as an unresolved radio-loud AGN population but shifted down by a factor of $\sim$10. We interpret this as strong evidence that these high-luminosity PS sources evolve into large-scale radio-loud AGN. For local, low-luminosity PS sources, there is a surplus of PS sources, which we hypothesise to be the addition of frustrated PS sources that do not evolve into large-scale AGN., Comment: 21 pages, accepted in Astronomy & Astrophysics

High redshift quasars ($z>5$) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources at $z\ge5$ are extremely rare and therefore we have started a campaign to search for new high-$z$ quasars by combining an optical dropout selection driven by the $g$, $r$, and $z$ bands from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys with low-frequency radio observations from the LOFAR Two-metre Sky Survey (LoTSS). Currently, LoTSS covers a large fraction of the northern sky (5720 deg$^2$) to such a depth (median noise level of 83 $\mu$Jy beam$^{-1}$) that about 30% of the general quasar population is detected $-$ which is a factor of 5-10 more than previous large sky radio surveys such as NVSS and FIRST, respectively. In this paper, we present the discovery of 20 new quasars (and the independent confirmation of 4) between $4.9\leq z\leq 6.6$. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion of $R=f_{5\text{GHz}}/f_{4400A} > 10$, with the full sample spanning $R\sim$6-1000, thereby more than doubling the sample of known radio-loud quasars at $z \ge 5$. Our radio detection requirement strongly decreases the contamination of stellar sources and allows one to select these quasars in a broad redshift range. Despite selecting our quasar candidates using fewer and less conservative colour restrictions, both the optical and near-infrared colours, Ly$\alpha$ emission line properties, and dust reddening, $E(B-V)$, measurements of our quasar sample do not deviate from the known radio-quiet quasar population, suggesting similar optical quasar properties of the radio-loud and radio-quiet quasar population at high-$z$. Our campaign demonstrates the potential for discovering new high-$z$ quasar populations through next generation radio continuum surveys., Comment: 19 pages, 11 figures, Tab A1 corrected, accepted for publication in A&A

Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate., Comment: Accepted for publication in MNRAS. 19 pages, 15 figures, 2 tables

In this paper we present photometric redshifts for 2.7 million galaxies in the XMM-LSS and COSMOS fields, both with rich optical and near-infrared data from VISTA and HyperSuprimeCam. Both template fitting (using galaxy and Active Galactic Nuclei templates within LePhare) and machine learning (using GPz) methods are run on the aperture photometry of sources selected in the Ks-band. The resulting predictions are then combined using a Hierarchical Bayesian model, to produce consensus photometric redshift point estimates and probability distribution functions that outperform each method individually. Our point estimates have a root mean square error of ~0.08-0.09, and an outlier fraction of ~3-4 percent when compared to spectroscopic redshifts. We also compare our results to the COSMOS2020 photometric redshifts, which contains fewer sources, but had access to a larger number of bands and greater wavelength coverage, finding that comparable photo-z quality can be achieved (for bright and intermediate luminosity sources where a direct comparison can be made). Our resulting redshifts represent the most accurate set of photometric redshifts (for a catalogue this large) for these deep multi-square degree multi-wavelength fields to date., Comment: Published in MNRAS, 15 figures. Photo-z catalogue described in the paper available online here: https://entangled.physics.ox.ac.uk/index.php/s/T1NQAMardpoOQMQ

Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies., Comment: 26 pages, 16 figures, accepted for publication in MNRAS. LF data available at: https://github.com/rohitk-10/AGN_LF_Kondapally22 . Author list error corrected in updated version

We present a novel method to identify candidate high redshift quasars (HzQs; ($z\gtrsim5.5$), which are unique probes of supermassive black hole growth in the early Universe, from large area optical/infrared photometric surveys. Using Gaussian Mixture Models to construct likelihoods and incorporate informed priors based on population statistics, our method uses a Bayesian framework to assign posterior probabilities that differentiate between HzQs and contaminating sources. We additionally include deep radio data to obtain informed priors. Using existing HzQ data in the literature, we set a posterior threshold that accepts ${\sim}90\%$ of known HzQs while rejecting $>99\%$ of contaminants such as dwarf stars or lower redshift galaxies. Running the probability selection on test samples of simulated HzQs and contaminants, we find that the efficacy of the probability method is higher than traditional colour cuts, decreasing the fraction of accepted contaminants by 86% while retaining a similar fraction of HzQs. As a test, we apply our method to the Pan-STARRS Data Release 1 (PS1) source catalogue within the HETDEX Spring field area on the sky, covering 400 sq. deg. and coinciding with deep radio data from the LOFAR Two-metre Sky Survey Data Release 1 (LoTSS DR1). From an initial sample of ${\sim}5\times10^5$ sources in PS1, our selection shortlists 251 candidate HzQs, which are further reduced to 63 after visual inspection. Shallow spectroscopic follow-up of 13 high probability HzQs resulted in the confirmation of a previously undiscovered quasar at $z=5.66$ with photometric colours $i-z = 1.4$, lying outside the typically probed regions when selecting HzQs based on colours. This discovery demonstrates the efficacy of our probabilistic HzQ selection method in selecting more complete HzQ samples, which holds promise when employed on large existing and upcoming photometric data sets., Comment: 15 pages, 8 figures. Accepted for publication in Astronomy & Astrophysics

Radio-loud active galaxies have two accretion modes [radiatively inefficient (RI) and radiatively efficient (RE)], with distinct optical and infrared signatures, and two jet dynamical behaviours, which in arcsec- to arcmin-resolution radio surveys manifest primarily as centre- or edge-brightened structures [Fanaroff-Riley (FR) class I and II]. The nature of the relationship between accretion mode and radio morphology (FR class) has been the subject of long debate. We present a comprehensive investigation of this relationship for a sample of 286 well-resolved radio galaxies in the LOFAR Two-metre Sky Survey Deep Fields (LoTSS-Deep) first data release, for which robust morphological and accretion mode classifications have been made. We find that two-thirds of luminous FRII radio galaxies are RI, and identify no significant differences in the visual appearance or source dynamic range (peak/mean surface brightness) of the RI and RE FRIIs, demonstrating that both RI and RE systems can produce FRII structures. We also find a significant population of low-luminosity FRIIs (predominantly RI), supporting our earlier conclusion that FRII radio structures can be produced at all radio luminosities. We demonstrate that in the luminosity range where both morphologies are present, the probability of producing FRI or FRII radio morphology is directly linked to stellar mass, while across all morphologies and luminosities, RE accretion occurs in systems with high specific star formation rate, presumably because this traces fuel availability. In summary, the relationship between accretion mode and radio morphology is very indirect, with host-galaxy environment controlling these two key parameters in different ways., Comment: 22 pages, 17 figures, 7 tables. Accepted for publication in MNRAS. Update with minor changes to match published version

Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $\mu$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 \sigma$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z<3$. We compare the radio loudness of the high-$z$ quasar sample to a lower-$z$ quasar sample at $z\sim2$ and find that the two radio loudness distributions are consistent with no evolution, although the low number of high-z quasars means that we cannot rule out weak evolution. Furthermore, we make a first order empirical estimate of the $z=6$ quasar radio luminosity function, which is used to derive the expected number of high-$z$ sources that will be detected in the completed LoTSS survey. This work highlights the fact that new deep radio observations can be a valuable tool in selecting high-$z$ quasar candidates for follow-up spectroscopic observations by decreasing contamination of stellar dwarfs and reducing possible selection biases introduced by strict colour cuts., Comment: 11 pages, 10 figures, accepted for publication in A&A

We examine the distribution of radio emission from ~42,000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-Metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star-formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the 'radio-quiet' and 'radio-loud' quasar regimes, without need for any explicit bimodality. The best-fit model parameters indicate that the star-formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z~2. For a model where star-formation rate scales as $SFR \propto L_{bol}^\alpha (1+z)^\beta$, we find $\alpha = 0.47 \pm 0.01$ and $\beta = 1.61 \pm 0.05$, in agreement with far-infrared studies. Quasars contribute ~0.15 per cent of the cosmic star-formation rate density at z=0.5, rising to 0.4 per cent by z=2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets., Comment: Accepted for publication in MNRAS. 21 pages

We investigate the radio properties of a sample of 53 sources selected at 850 $\mu$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their similarly bright submillimetre flux densities. We characterise their spectral shapes in terms of multi-band radio spectral indices. Finding strong spectral flattening at low frequencies in ~20% of sources, we investigate the differences between sources with extremely flat low-frequency spectra and those with `normal' radio spectral indices. As there are no other statistically significant differences between the two subgroups of our sample as split by the radio spectral index, we suggest that any differences are undetectable in galaxy-averaged properties that we can observe with our unresolved images, and likely relate to galaxy properties that we cannot resolve, on scales $\lesssim$ 1 kpc. We attribute the observed spectral flattening in the radio to free-free absorption, proposing that those sources with significant low-frequency spectral flattening have a clumpy distribution of star-forming gas. We estimate an average spatial extent of absorbing material of at most several hundred parsecs to produce the levels of absorption observed in the radio spectra. This estimate is consistent with the highest-resolution observations of submillimetre galaxies in the literature, which find examples of non-uniform dust distributions on scales of ~100 pc, with evidence for clumps and knots in the interstellar medium. Additionally, we find two bright (> 6 mJy) submm sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4., Comment: 15 pages, 10 figures, accepted by A&A

Lyman alpha emitters (LAEs) in the Epoch of Reionization (EoR) offer valuable probes of early galaxy evolution and the process of reionization; however, the exact evolution of their abundance and the nature of their emission remain open questions. We combine samples of 229 and 349 LAE candidates at $z=5.7$ and $z=6.6,$ respectively, from the SILVERRUSH narrowband survey with deep Low Frequency Array (LOFAR) radio continuum observations in the ELAIS-N1 field to search for radio galaxies in the EoR and study the low-frequency radio properties of $z\gtrsim5.7$ LAE emitters. Our LOFAR observations reach an unprecedented noise level of $\sim20\,\mu$Jy beam$^{-1}$ at 150MHz, and we detect five candidate LAEs at $>5\sigma$ significance. Based on detailed spectral energy distribution modelling of independent multi-wavelength observations, we conclude that these sources are likely [OII] emitters at $z=1.47$, yielding no reliable $z\gtrsim5.7$ radio galaxy candidates. We examine the 111 $z=5.7$ and $z=6.6$ LAE candidates from our panchromatic photometry catalogue that are undetected by LOFAR, finding contamination rates of 81-92% for the $z=5.7$ and $z=6.6$ subset of the LAE candidate samples. This subset is biased towards brighter magnitudes and redder near-infrared colours. The contamination rates of the full sample will therefore likely be lower than the reported values. Contamination is lowered significantly through constraints on the near-infrared colours, highlighting the need for infrared observations to robustly identify bright LAEs in narrowband surveys. Finally, the stacking of radio continuum observations for the robust LAE samples yields 2$\sigma$ upper limits on radio luminosity of 8.2$\times$10$^{23}$ and 8.7$\times$10$^{23}$ W Hz$^{-1}$ at $z=5.7$ and $6.6$, respectively, corresponding to limits on their median star-formation rates of $<$53 and $<$56 M$_{\odot}$ yr$^{-1}$., Comment: 17 pages, 11 figures, accepted for publication in A&A. The abstract is abridged. This paper is part of the 1st data release of the LoTSS Deep Fields

The Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey split across multiple tiers over the northern sky. The first LoTSS Deep Fields data release consists of deep radio continuum imaging at 150 MHz of the Bo\"{o}tes, European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1), and Lockman Hole fields, down to rms sensitivities of $\sim$32, 20, and 22 $\mu$Jy beam$^{-1}$, respectively. In this paper we present consistent photometric redshift (photo-$z$) estimates for the optical source catalogues in all three fields - totalling over 7 million sources ($\sim5$ million after limiting to regions with the best photometric coverage). Our photo-$z$ estimation uses a hybrid methodology that combines template fitting and machine learning and is optimised to produce the best possible performance for the radio continuum selected sources and the wider optical source population. Comparing our results with spectroscopic redshift samples, we find a robust scatter ranging from 1.6 to 2% for galaxies and 6.4 to 7% for identified optical, infrared, or X-ray selected active galactic nuclei (AGN). Our estimated outlier fractions ($\left | z_{\text{phot}} - z_{\text{spec}} \right | / (1+z_{\text{spec}}) > 0.15$) for the corresponding subsets range from 1.5 to 1.8% and 18 to 22%, respectively. Replicating trends seen in analyses of previous wide-area radio surveys, we find no strong trend in photo-$z$ quality as a function of radio luminosity for a fixed redshift. We exploit the broad wavelength coverage available within each field to produce galaxy stellar mass estimates for all optical sources at $z < 1.5$. Stellar mass functions derived for each field are used to validate our mass estimates, with the resulting estimates in good agreement between each field and with published results from the literature., Comment: 23 pages, 14 figures. Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields - Electronic data catalogues will be made available on journal publication

In this paper, we investigate the relationship between 150MHz luminosity and star formation rate (the SFR-L150 relation) using 150MHz measurements for a near-infrared selected sample of 118,517 $z<1$ galaxies. New radio survey data offer compelling advantages for studying star formation in galaxies, with huge increases in sensitivity, survey speed and resolution over previous generation surveys, and remaining impervious to extinction. The LOFAR Surveys Key Science Project is transforming our understanding of the low-frequency radio sky, with the 150MHz data over the ELAIS-N1 field reaching an RMS sensitivity of 20uJy/beam over 10 deg$^2$ at 6" resolution. All of the galaxies studied have SFR and stellar mass estimates derived from energy balance SED fitting, using redshifts and aperture-matched forced photometry from the LOFAR Two-metre Sky Survey (LoTSS) deep fields data release. The impact of active galactic nuclei is minimised by leveraging the deep ancillary data alongside outlier-resistant median-likelihood methods. We find a linear and non-evolving SFR-L150 relation, apparently consistent with expectations based on calorimetric arguments, down to the lowest SFRs. However, we also recover compelling evidence for stellar mass dependence in line with previous work on this topic, in the sense that higher mass galaxies have a larger 150MHz luminosity at a given SFR, suggesting that the overall agreement with calorimetric arguments may be a coincidence. We conclude that in the absence of AGN, 150MHz observations can be used to measure accurate galaxy SFRs out to $z=1$ at least, but it is necessary to account for stellar mass in order to obtain 150MHz-derived SFRs accurate to <0.5 dex. Our best-fit relation is $\log_{10} (L_\mathrm{150 MHz} / W\,Hz^{-1}) = (0.90\pm 0.01) \log_{10}(\psi/M_\odot\,\mathrm{yr}^{-1}) + (0.33 \pm 0.04) \log_{10} (M/10^{10}M_\odot) + 22.22 \pm 0.02$. (Abridged), Comment: 17 pages, 17 figures, accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields

We present the source associations, cross-identifications, and multi-wavelength properties of the faint radio source population detected in the deep tier of the LOFAR Two Metre Sky Survey (LoTSS): the LoTSS Deep Fields. The first LoTSS Deep Fields data release consists of deep radio imaging at 150~MHz of the ELAIS-N1, Lockman Hole, and Bo\"{o}tes fields, down to RMS sensitives of around 20, 22, and 32$~\mu$Jy\,beam$^{-1}$, respectively. These fields are some of the best studied extra-galactic fields in the northern sky, with existing deep, wide-area panchromatic photometry from X-ray to infrared wavelengths, covering a total of $\approx$~26~\mbox{deg$^{2}$}. We first generated improved multi-wavelength catalogues in ELAIS-N1 and Lockman Hole; combined with the existing catalogue for Bo\"{o}tes, we present forced, matched aperture photometry for over 7.2 million sources across the three fields. We identified multi-wavelength counterparts to the radio detected sources, using a combination of the Likelihood Ratio method and visual classification, which greatly enhances the scientific potential of radio surveys and allows for the characterisation of the photometric redshifts and the physical properties of the host galaxies. The final radio-optical cross-matched catalogue consists of 81\,951 radio-detected sources, with counterparts identified and multi-wavelength properties presented for 79\,820 ($>$97\%) sources. We also examine the properties of the host galaxies, and through stacking analysis find that the radio population with no identified counterpart is likely dominated by AGN at $z\sim3-4$. This dataset contains one of the largest samples of radio-selected star-forming galaxies and active galactic nuclei (AGN) at these depths, making it ideal for studying the history of star-formation, and the evolution of galaxies and AGN across cosmic time., Comment: 34 pages, 14 figures and 14 tables. Accepted for publication in A&A special issue. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication

We present a new, ambitious survey performed with the Chandra X-ray Observatory of the 9.3 deg$^2$ Bo\"otes field of the NOAO Deep Wide-Field Survey. The wide field probes a statistically representative volume of the Universe at high redshift. The Chandra Deep Wide-Field Survey exploits the excellent sensitivity and angular resolution of Chandra over a wide area, combining 281 observations spanning 15 years, for a total exposure time of 3.4 Ms, and detects 6891 X-ray point sources down to limiting fluxes of $4.7\times10^{-16}$, $1.5\times10^{-16}$, and $9\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$, in the $0.5-7$ keV, $0.5-2$ keV, and $2-7$ keV bands, respectively. The robustness and reliability of the detection strategy is validated through extensive, state-of-the-art simulations of the whole field. Accurate number counts, in good agreement with previous X-ray surveys, are derived thanks to the uniquely large number of point sources detected, which resolve $65.0 \pm 12.8\%$ of the cosmic X-ray background between $0.5-2$ keV and $81.0 \pm 11.5\%$ between $2-7$ keV. Exploiting the wealth of multi-wavelength data available on the field, we assign redshifts to $\sim 94\%$ of the X-ray sources, estimate their obscuration and derive absorption-corrected luminosities. We provide an electronic catalog containing all the relevant quantities needed for future investigations., Comment: Accepted for publication in ApJS. The catalogs associated to this paper can be accessed at http://www.dartmouth.edu/~hickox/cdwfs.php

Various studies have laid claim to finding an alignment of the polarization vectors or radio jets of active galactic nuclei (AGN) over large distances, but these results have proven controversial and so far, there is no clear explanation for this observed alignment. To investigate this case further, we tested the hypothesis that the position angles of radio galaxies are randomly oriented in the sky by using data from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS). A sample of 7,555 double-lobed radio galaxies was extracted from the list of 318,520 radio sources in the first data release of LoTSS at 150 MHz. We performed statistical tests for uniformity of the two-dimensional (2D) orientations for the complete 7,555 source sample. We also tested the orientation uniformity in three dimensions (3D) for the 4,212 source sub-sample with photometric or spectroscopic redshifts. Our sample shows a significant deviation from uniformity (p-value < $10^{-5}$) in the 2D analysis at angular scales of about four degrees, mainly caused by sources with the largest flux densities. No significant alignment was found in the 3D analysis. Although the 3D analysis has access to fewer sources and suffers from uncertainties in the photometric redshift, the lack of alignment in 3D points towards the cause of the observed effect being unknown systematics or biases that predominantly affect the brightest sources, although this has yet to be demonstrated irrefutably and should be the subject of subsequent studies., Comment: 10 pages, 13 figures. Accepted for publication in Astronomy & Astrophysics

Radio galaxies are known to go through cycles of activity, where phases of apparent quiescence can be followed by repeated activity of the central supermassive black hole. A better understanding of this cycle is crucial for ascertaining the energetic impact that the jets have on the host galaxy, but little is known about it. We used deep LOFAR images at 150 MHz of the Lockman Hole extragalactic field to select a sample of 158 radio sources with sizes $> 60^{\prime\prime}$ in different phases of their jet life cycle. Using a variety of criteria (e.g. core prominence combined with low-surface brightness of the extended emission and steep spectrum of the central region) we selected a subsample of candidate restarted radio galaxies representing between 13% and 15% of the 158 sources of the main sample. We compare their properties to the rest of the sample, which consists of remnant candidates and active radio galaxies. Optical identifications and characterisations of the host galaxies indicate similar properties for candidate restarted, remnant, and active radio galaxies, suggesting that they all come from the same parent population. The fraction of restarted radio galaxies is slightly higher with respect to remnants, suggesting that the restarted phase can often follow after a relatively short remnant phase (the duration of the remnant phase being a few times 10$^{7}$ years). This confirms that the remnant and restarted phases are integral parts of the life cycle of massive elliptical galaxies. A preliminary investigation does not suggest a strong dependence of this cycle on the environment surrounding any given galaxy., Comment: 16 pages, 9 figures, 4 tables, A&A accepted, companion paper to Shabala et al., also in today arXiv posting

We present VLT/SINFONI observations of 35 quasars at 2.1 < z < 3.2, the majority of which were selected from the Clusters Around Radio-Loud AGN (CARLA) survey. CARLA quasars have large CIV-based black hole (BH) masses (M(BH) > 10^9 Msun) and powerful radio emission (P(500MHz) > 27.5 W/Hz). We estimate Ha-based M(BH), finding a scatter of 0.35 dex compared to CIV. We evaluate several recipes for correcting CIV-based masses, which reduce the scatter to 0.24 dex. The radio power of the radio-loud quasars is at most weakly correlated with the interconnected quantities Ha-width, L(5100A) and M(BH), suggesting that it is governed by different physical processes. However, we do find a strong inverse correlation between CIV blueshift and radio power linked to higher Eddington ratios and L(5100A). Under standard assumptions, the BH growth time is longer than the cosmic age for many CARLA quasars, suggesting that they must have experienced more efficient growth in the past. If these BHs were growing from seeds since the epoch of reionization, it is possible that they grew at the Eddington limit like the quasars at z ~ 6-7, and then continued to grow at the reduced rates observed until z ~ 2-3. Finally, we study the relation between M(BH) and environment, finding a weak positive correlation between M(BH) and galaxy density measured by CARLA., Comment: MNRAS, accepted (28 pages, 22 figures)

Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift z~0:05) and improve the use of the rest-frame 150-MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction. Methods. We cross-match the 60-um selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150-MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers. Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFR_Halpha based on the Halpha line luminosity, SFR_60 based on the rest-frame 60-um luminosity and SFR_IR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (Lsun) = 1.35(0.06) x log SFR_Halpha (Msun/yr) + 3.20(0.06), log L150 (Lsun) = 1.31(0.05) x log SFR_60 (Msun/yr) + 3.14(0.06), and log L150 (Lsun) = 1.37(0.05) x log SFR_IR (Msun/yr) + 3.09(0.05), which show excellent agreement with each other., Comment: 13 pages, 13 figures, accepted for publication in A&A

We present spectral energy distributions (SEDs) of 41 active galactic nuclei, derived from multiwavelength photometry and archival spectroscopy. All of the SEDs span at least 0.09 to 30 micron, but in some instances wavelength coverage extends into the X-ray, far-infrared and radio. For some AGNs we have fitted the measured far-infrared photometry with greybody models, while radio flux density measurements have been approximated by power-laws or polynomials. We have been able to fill some of the gaps in the spectral coverage using interpolation or extrapolation of simple models. In addition to the 41 individual AGN SEDs, we have produced 72 Seyfert SEDs by mixing SEDs of the central regions of Seyferts with galaxy SEDs. Relative to the literature, our templates have broader wavelength coverage and/or higher spectral resolution. We have tested the utility of our SEDs by using them to generate photometric redshifts for 0 < z < 6.12 AGNs in the Bootes field (selected with X-ray, IR and optical criteria) and, relative to SEDs from the literature, they produce comparable or better photometric redshifts with reduced flux density residuals., Comment: Accepted for publication in MNRAS. 19 pages, 14 figures, 4 tables. The AGN spectral energy distributions, photometry and images are available via https://archive.stsci.edu/hlsp/agnsedatlas

The relative positions of the high and low surface brightness regions of radio-loud active galaxies in the 3CR sample were found by Fanaroff and Riley to be correlated with their luminosity. We revisit this canonical relationship with a sample of 5805 extended radio-loud AGN from the LOFAR Two-Metre Sky Survey (LoTSS), compiling the most complete dataset of radio-galaxy morphological information obtained to date. We demonstrate that, for this sample, radio luminosity does *not* reliably predict whether a source is edge-brightened (FRII) or centre-brightened (FRI). We highlight a large population of low-luminosity FRIIs, extending three orders of magnitude below the traditional FR break, and demonstrate that their host galaxies are on average systematically fainter than those of high-luminosity FRIIs and of FRIs matched in luminosity. This result supports the jet power/environment paradigm for the FR break: low-power jets may remain undisrupted and form hotspots in lower mass hosts. We also find substantial populations that appear physically distinct from the traditional FR classes, including candidate restarting sources and ``hybrids''. We identify 459 bent-tailed sources, which we find to have a significantly higher SDSS cluster association fraction (at $z<0.4$) than the general radio-galaxy population, similar to the results of previous work. The complexity of the LoTSS faint, extended radio sources demonstrates the need for caution in the automated classification and interpretation of extended sources in modern radio surveys, but also reveals the wealth of morphological information such surveys will provide and its value for advancing our physical understanding of radio-loud AGN., Comment: 21 pages, 4 tables, 16 Figures. Accepted for publication on Monthly Notices of the Royal Astronomical Society (MNRAS). Our LoTSS morphological catalogue containing classifications for 5805 extended radio-loud AGN is available from www.lofar-surveys.org/releases.html

We present spectra and near-infrared images of a sample of faint radio sources initially selected as promising high-redshift radio galaxy (HzRG) candidates. We have determined redshifts for a total of 13 radio galaxies with redshifts ranging from $0.52\le z \le 5.72$. Our sample probes radio luminosities that are almost an order of magnitude fainter than previous large samples at the highest redshifts. We use near-infrared photometry for a subsample of these galaxies to calculate stellar masses using simple stellar population models, and find stellar masses to be in the range $10^{10.8} - 10^{11.7} M_\odot$. We then compare our faint radio galaxies with brighter radio galaxies at $z\ge2$ from the literature. We find that fainter radio galaxies have lower Ly$\alpha$ luminosities and narrower full width at half maxima compared to the bright ones, implying photoionisation by weaker AGN. We also rule out the presence of strong shocks in faint HzRGs. The stellar masses determined for faint HzRGs are lower than what has been observed for brighter ones. We find that faint HzRG population in the redshift range $2-4$ forms a bridge between star-forming and narrow-line AGN, whereas the ones at $z>4$ are dominated by star-formation, and may be building up their stellar mass through cold accretion of gas. Finally, we show that the overall redshift evolution of radio sizes at $z>2$ is fully compatible with increased inverse Compton scattering losses at high redshifts., Comment: 26 pages (including 2 appendices), 17 figures, accepted for publication in MNRAS

Giant radio galaxies (GRGs) are a subclass of radio galaxies which have grown to megaparsec scales. GRGs are much rarer than normal sized radio galaxies (< 0.7 Mpc) and the reason for their gigantic sizes is still debated. Here, we report the biggest sample of GRGs identified to date. These objects were found in the LOFAR Two-metre Sky Survey (LoTSS) first data release images, which cover a 424 square degrees region. Of the 239 GRGs found, 225 are new discoveries. The GRGs in our sample have sizes ranging from 0.7 to 3.5 Mpc and have redshifts (z) between 0.1 and 2.3. Seven GRGs have sizes above 2 Mpc and one has a size of ~ 3.5 Mpc. The sample contains 40 GRGs hosted by spectroscopically confirmed quasars. Here, we present the search techniques employed and the resulting catalogue of the newly discovered large sample of GRGs along with their radio properties. We, here also show for the first time that the spectral index of GRGs is similar to that of normal sized radio galaxies, indicating that most of the GRG population is not dead or is not like remnant type radio galaxy. We find 20/239 GRGs in our sample are located at the centres of clusters and present our analysis on their cluster environment and radio morphology., Comment: Accepted for publication in Astronomy & Astrophysics

We report optical observations of TGSS J1054+5832, a candidate high-redshift ($z=4.8\pm2$) steep-spectrum radio galaxy, in $r$ and $i$ bands using the faint object spectrograph and camera mounted on 3.6-m Devasthal Optical Telescope (DOT). The source previously detected at 150 MHz from Giant Meterwave Radio Telescope and at 1420 MHz from Very Large Array has a known counterpart in near-infrared bands with $K$-band magnitude of AB 22. The source is detected in $i$-band with AB$24.3\pm0.2$ magnitude in the DOT images presented here. The source remains undetected in the $r$-band image at a 2.5$\sigma$ depth of AB 24.4 mag over an $1.2''\times1.2''$ aperture. An upper limit to $i-K$ color is estimated to be $\sim$2.3, suggesting youthfulness of the galaxy with active star formation. These observations highlight the importance and potential of the 3.6-m DOT for detections of faint galaxies., Comment: 7 pages, 2 figures, Accepted in Jr. of Astrophysics & Astronomy (JApA)

We present a study of the low-frequency radio properties of broad absorption line quasars (BALQSOs) from the LOFAR Two-metre Sky-Survey Data Release 1 (LDR1). The value-added LDR1 catalogue contains Pan-STARRS counterparts, which we match with the Sloan Digital Sky Survey (SDSS) DR7 and DR12 quasar catalogues. We find that BALQSOs are twice as likely to be detected at 144$\,$MHz than their non-BAL counterparts, and BALQSOs with low-ionisation species present in their spectra are three times more likely to be detected than those with only high-ionisation species. The BALQSO fraction at 144$\,$MHz is constant with increasing radio luminosity, which is inconsistent with previous results at 1.4$\,$GHz, indicating that observations at the different frequencies may be tracing different sources of radio emission. We cross-match radio sources between the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey and LDR1, which provides a bridge via the LDR1 Pan-STARRS counterparts to identify BALQSOs in SDSS. Consequently we expand the sample of BALQSOs detected in FIRST by a factor of three. The LDR1-detected BALQSOs in our sample are almost exclusively radio-quiet (\logr $\,<2$), with radio sizes at 144$\,$MHz typically less than $200\,$kpc; these radio sizes tend to be larger than those at 1.4$\,$GHz, suggesting more extended radio emission at low frequencies. We find that although the radio detection fraction increases with increasing balnicity index (BI), there is no correlation between BI and either low-frequency radio power or radio-loudness. This suggests that both radio emission and BI may be linked to the same underlying process, but are spatially distinct phenomena., Comment: 16 figures, 3 tables and 19 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation., Comment: 21 pages, 11 colour figures. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication of the journal

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections., Comment: 7 figures, 2 tables and 10 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release., Comment: 8 figures, 5 tables and 23 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication in the journal

An understanding of the relationship between radio-loud active galaxies and their large-scale environments is essential for realistic modelling of radio-galaxy evolution and environmental impact, for understanding AGN triggering and life cycles, and for calibrating galaxy feedback in cosmological models. We use the LOFAR Two-Metre Sky Survey (LoTSS) Data Release 1 catalogues to investigate this relationship. We cross-matched 8,745 radio-loud AGN with 0.08

This paper presents a study of the local radio source population, by cross-comparing the data from the first data release (DR1) of the LOFAR Two-Metre Sky Survey (LoTSS) with the Sloan Digital Sky Survey (SDSS) DR7 main galaxy spectroscopic sample. The LoTSS DR1 provides deep data (median rms noise of 71 $\mathrm{\mu}$Jy at 150 MHz) over 424 square degrees of sky, which is sufficient to detect 10615 (32 per cent) of the SDSS galaxies over this sky area. An improved method to separate active galactic nuclei (AGN) accurately from sources with radio emission powered by star formation (SF) is developed and applied, leading to a sample of 2121 local ($z < 0.3$) radio AGN. The local 150 MHz luminosity function is derived for radio AGN and SF galaxies separately, and the good agreement with previous studies at 1.4 GHz suggests that the separation method presented is robust. The prevalence of radio AGN activity is confirmed to show a strong dependence on both stellar and black hole masses, remarkably reaching a fraction of 100 per cent of the most massive galaxies ($> 10^{11} \mathrm{M_{\odot}}$) displaying radio-AGN activity with $L_{\rm 150 MHz} \geq 10^{21}$W Hz$^{-1}$; thus, the most massive galaxies are always switched on at some level. The results allow the full Eddington-scaled accretion rate distribution (a proxy for the duty cycle) to be probed for massive galaxies. More than 50 per cent of the energy is released during the $\le 2$ per cent of the time spent at the highest accretion rates, $L_{\mathrm{mech}}/L_{\mathrm{Edd}} > 10^{-2.5}$. Stellar mass is shown to be a more important driver of radio-AGN activity than black hole mass, suggesting a possible connection between the fuelling gas and the surrounding halo. This result is in line with models in which these radio AGN are essential for maintaining the quenched state of galaxies at the centres of hot gas haloes., Comment: 14 pages, 8 figures, 3 tables. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

We report the discovery of the most distant radio galaxy to date, TGSS1530 at a redshift of $z=5.72$ close to the presumed end of the Epoch of Reionisation. The radio galaxy was selected from the TGSS ADR1 survey at 150 MHz for having to an ultra-steep spectral index, $\alpha^{\textrm{150 MHz}}_{\textrm{1.4 GHz}} = -1.4$ and a compact morphology obtained using VLA imaging at 1.4 GHz. No optical or infrared counterparts for it were found in publicly available sky surveys. Follow-up optical spectroscopy at the radio position using GMOS on Gemini North revealed the presence of a single emission line. We identify this line as Lyman alpha at $z=5.72$, because of its asymmetric line profile, the absence of other optical/UV lines in the spectrum and a high equivalent width. With a Ly$\alpha$ luminosity of $5.7 \times 10^{42}$ erg s$^{-1}$ and a FWHM of $370$ km s$^{-1}$, TGSS1530 is comparable to `non-radio' Lyman alpha emitters (LAEs) at a similar redshift. However, with a radio luminosity of $\log L_{\textrm{150 MHz}} = 29.1$ W Hz$^{-1}$ and a deconvolved physical size $3.5$ kpc, its radio properties are similar to other known radio galaxies at $z>4$. Subsequent $J$ and $K$ band imaging using LUCI on the Large Binocular Telescope resulted in non-detection of the host galaxy down to $3\sigma$ limits of $J>24.4$ and $K>22.4$ (Vega). The $K$ band limit is consistent with $z>5$ from the $K-z$ relation for radio galaxies, suggesting stellar mass limits using simple stellar population models of $M_{\textrm{stars}}< 10^{10.5}$ $M_\odot$. Its high redshift coupled with relatively small radio and Ly$\alpha$ sizes suggest that TGSS1530 may be a radio galaxy in an early phase of its evolution., Comment: 11 pages, 7 figures. Latest published version 10 October 2018

Ultra-steep spectrum (USS) radio sources are good tracers of powerful radio galaxies at $z > 2$. Identification of even a single bright radio galaxy at $z > 6$ can be used to detect redshifted 21cm absorption due to neutral hydrogen in the intervening IGM. Here we describe a new sample of high-redshift radio galaxy (HzRG) candidates constructed from the TGSS ADR1 survey at 150 MHz. We employ USS selection ($\alpha \le -1.3$) in $\sim10000$ square degrees, in combination with strict size selection and non-detections in all-sky optical and infrared surveys. We apply flux density cuts that probe a unique parameter space in flux density ($50 < S_{\textrm{150}} < 200$ mJy) to build a sample of 32 HzRG candidates. Follow-up Karl G. Jansky Very Large Array (VLA) observations at 1.4 GHz with an average beam size of $1.3$ arcseconds ($"$) revealed $\sim 48\%$ of sources to have a single radio component. P-band (370 MHz) imaging of 17 of these sources revealed a flattening radio SED for ten sources at low frequencies, which is expected from compact HzRGs. Two of our sources lie in fields where deeper multi-wavelength photometry and ancillary radio data are available and for one of these we find a best-fit photo-z of $4.8 \pm 2.0$. The other source has $z_{\textrm{phot}}=1.4 \pm 0.1$ and a small angular size ($3.7"$), which could be associated with an obscured star forming galaxy or with a "dead" elliptical. One USS radio source not part of the HzRG sample but observed with the VLA nonetheless is revealed to be a candidate giant radio galaxy with a host galaxy photo-z of $1.8\pm0.5$, indicating a size of 875 kpc., Comment: 17 pages, 13 figures, accepted for publication in MNRAS

This paper presents a study of the redshift evolution of radio-loud active galactic nuclei (AGN) as a function of the properties of their galaxy hosts in the Bo\"otes field. To achieve this we match low-frequency radio sources from deep $150$-MHz LOFAR observations to an $I$-band-selected catalogue of galaxies, for which we have derived photometric redshifts, stellar masses and rest-frame colours. We present spectral energy distribution (SED) fitting to determine the mid-infrared AGN contribution for the radio sources and use this information to classify them as High- versus Low-Excitation Radio Galaxies (HERGs and LERGs) or Star-Forming galaxies. Based on these classifications we construct luminosity functions for the separate redshift ranges going out to $z = 2$. From the matched radio-optical catalogues, we select a sub-sample of $624$ high power ($P_{150\mathrm{\,MHz}}>10^{25}$ W Hz$^{-1}$) radio sources between $0.5 \leq z < 2$. For this sample, we study the fraction of galaxies hosting HERGs and LERGs as a function of stellar mass and host galaxy colour. The fraction of HERGs increases with redshift, as does the fraction of sources in galaxies with lower stellar masses. We find that the fraction of galaxies that host LERGs is a strong function of stellar mass as it is in the local Universe. This, combined with the strong negative evolution of the LERG luminosity functions over this redshift range, is consistent with LERGs being fuelled by hot gas in quiescent galaxies., Comment: 25 pages, 21 figures, accepted for publication in MNRAS

In these proceedings we highlight the primary scientific goals and design of the WEAVE-LOFAR survey, which will use the new WEAVE spectrograph on the 4.2m William Herschel Telescope to provide the primary source of spectroscopic information for the LOFAR Surveys Key Science Project. Beginning in 2018, WEAVE-LOFAR will generate more than 10$^6$ R=5000 365-960 nm spectra of low-frequency selected radio sources, across three tiers designed to efficiently sample the redshift-luminosity plane, and produce a data set of enormous legacy value. The radio frequency selection, combined with the high multiplex and throughput of the WEAVE spectrograph, make obtaining redshifts in this way very efficient, and we expect that the redshift success rate will approach 100 per cent at $z < 1$. This unprecedented spectroscopic sample - which will be complemented by an integral field component - will be transformational in key areas, including studying the star formation history of the Universe, the role of accretion and AGN-driven feedback, properties of the epoch of reionisation, cosmology, cluster haloes and relics, as well as the nature of radio galaxies and protoclusters. Each topic will be addressed in unprecedented detail, and with the most reliable source classifications and redshift information in existence., Comment: 10 pages, 4 figures, proceedings of the SF2A conference, Lyon, 2016