Showing 1–50 of 66 results for author: Duncan, K J

Flux dependence of redshift distribution and clustering of LOFAR radio sources

Authors: Nitesh Bhardwaj, Dominik J. Schwarz, Catherine L. Hale, Kenneth J. Duncan, Stefano Camera, Caroline S. Heneka, Szymon J. Nakoneczny, Huub J. A. Röttgering, Thilo M. Siewert, Prabhakar Tiwari, Jinglan Zheng, George Miley, Cyril Tasse

Abstract: In this work we study the flux density dependence of the redshift distribution of low-frequency radio sources observed in the LOFAR Two-metre Sky Survey (LoTSS) deep fields and apply it to estimate the clustering length of the large-scale structure of the Universe, examining flux density limited samples (1 mJy, 2 mJy, 4 mJy and 8 mJy) of LoTSS wide field radio sources. We utilise and combine the p… ▽ More In this work we study the flux density dependence of the redshift distribution of low-frequency radio sources observed in the LOFAR Two-metre Sky Survey (LoTSS) deep fields and apply it to estimate the clustering length of the large-scale structure of the Universe, examining flux density limited samples (1 mJy, 2 mJy, 4 mJy and 8 mJy) of LoTSS wide field radio sources. We utilise and combine the posterior probability distributions of photometric redshift determinations for LoTSS deep field observations from three different fields (Boötes, Lockman hole and ELAIS-N1, together about $26$ square degrees of sky), which are available for between $91\%$ to $96\%$ of all sources above the studied flux density thresholds and observed in the area covered by multi-frequency data. We estimate uncertainties by a bootstrap method. We apply the inferred redshift distribution on the LoTSS wide area radio sources from the HETDEX field (LoTSS-DR1; about $424$ square degrees) and make use of the Limber approximation and a power-law model of three dimensional clustering to measure the clustering length, $r_0$, for various models of the evolution of clustering. We find that the redshift distributions from all three LoTSS deep fields agree within expected uncertainties. We show that the radio source population probed by LoTSS at flux densities above $1$ mJy has a median redshift of at least $0.9$. At $2$ mJy, we measure the clustering length of LoTSS radio sources to be $r_0 = (10.1\pm 2.6) \ h^{-1}$Mpc in the context of the comoving clustering model. Our findings are in agreement with measurements at higher flux density thresholds at the same frequency and with measurements at higher frequencies in the context of the comoving clustering model. △ Less

Submitted 25 June, 2024; originally announced June 2024.

The LOFAR Two-metre Sky Survey: The nature of the faint source population and SFR-radio luminosity relation using Prospector

Authors: Soumyadeep Das, Daniel J. B. Smith, Paul Haskell, Martin J. Hardcastle, Philip N. Best, Kenneth J. Duncan, Marina I. Arnaudova, Shravya Shenoy, Rohit Kondapally, Rachel K. Cochrane, Alyssa B. Drake, Gülay Gürkan, Katarzyna Małek, Leah K. Morabito, Isabella Prandoni

Abstract: Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospecto… ▽ More Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospector, to model up to 26 bands of UV$-$far-infrared aperture-matched photometry for $\sim$31,000 sources in the ELAIS-N1 field from the LOFAR Two-Metre Sky Survey (LoTSS) Deep fields first data release. One of a new generation of SED-fitting codes, Prospector accounts for potential contributions from radiative active galactic nuclei (AGN) when estimating galaxy properties, including star formation rates (SFRs) derived using non-parametric star formation histories. Combining this information with radio luminosities, we classify 92 per cent of the radio sources as a star-forming galaxy, high/low-excitation radio galaxy, or radio-quiet AGN and study the population demographics as a function of 150 MHz flux density, luminosity, SFR, stellar mass, redshift and apparent $r$-band magnitude. Finally, we use Prospector SED fits to investigate the SFR$-$150 MHz luminosity relation for a sample of $\sim$$133,000~3.6~μ$m-selected $z<1$ sources, finding that the stellar mass dependence is significantly weaker than previously reported, and may disappear altogether at $\log_{10} (\mathrm{SFR}/M_\odot~\mathrm{yr}^{-1}) > 0.5$. This approach makes it significantly easier to classify radio sources from LoTSS and elsewhere, and may have important implications for future studies of star-forming galaxies at radio wavelengths. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 20 pages, 13 figures, accepted for publication in MNRAS

Constraining the giant radio galaxy population with machine learning and Bayesian inference

Authors: Rafaël I. J. Mostert, Martijn S. S. L. Oei, B. Barkus, Lara Alegre, Martin J. Hardcastle, Kenneth J. Duncan, Huub J. A. Röttgering, Reinout J. van Weeren, Maya Horton

Abstract: Large-scale sky surveys at low frequencies, like the LOFAR Two-metre Sky Survey (LoTSS), allow for the detection and characterisation of unprecedented numbers of giant radio galaxies (GRGs, or 'giants'). In this work, by automating the creation of radio--optical catalogues, we aim to significantly expand the census of known giants. We then combine this sample with a forward model to constrain GRG… ▽ More Large-scale sky surveys at low frequencies, like the LOFAR Two-metre Sky Survey (LoTSS), allow for the detection and characterisation of unprecedented numbers of giant radio galaxies (GRGs, or 'giants'). In this work, by automating the creation of radio--optical catalogues, we aim to significantly expand the census of known giants. We then combine this sample with a forward model to constrain GRG properties of cosmological interest. In particular, we automate radio source component association through machine learning and optical host identification for resolved radio sources. We create a radio--optical catalogue for the full LoTSS Data Release 2 (DR2) and select all possible giants. We combine our candidates with an existing catalogue of LoTSS DR2 crowd-sourced GRG candidates and visually confirm or reject them. To infer intrinsic GRG properties from GRG observations, we develop further a population-based forward model that takes into account selection effects and constrain its parameters using Bayesian inference. We confirm 5,647 previously unknown giants from the crowd-sourced catalogue and 2,597 previously unknown giants from the ML-driven catalogue. Our confirmations and discoveries bring the total number of known giants to at least 11,585. We predict a comoving GRG number density $n_\mathrm{GRG} = 13 \pm 10\ (100\ \mathrm{Mpc})^{-3}$, close to a recent estimate of the number density of luminous non-giant radio galaxies. We derive a current-day GRG lobe volume-filling fraction $V_\mathrm{GRG-CW}(z = 0) = 1.4 \pm 1.1 \cdot 10^{-5}$ in clusters and filaments of the Cosmic Web. Our analysis suggests that giants are more common than previously thought. Moreover, tentative results imply that it is possible that magnetic fields once contained in giants pervade a significant ($\gtrsim 10\%$) fraction of today's Cosmic Web. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 25 pages, 15 figures, 3 tables, accepted for publication in Astronomy & Astrophysics

A novel Bayesian approach for decomposing the radio emission of quasars: I. Modelling the radio excess in red quasars

Authors: B. -H. Yue, P. N. Best, K. J. Duncan, G. Calistro-Rivera, L. K. Morabito, J. W. Petley, I. Prandoni, H. J. A. Röttgering, D. J. B. Smith

Abstract: 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 dist… ▽ More 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. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS. 19 pages, 13 figures, 1 table

Widespread AGN feedback in a forming brightest cluster galaxy at $z=4.1$ unveiled by JWST

Authors: Aayush Saxena, Roderik A. Overzier, Montserrat Villar-Martín, Tim Heckman, Namrata Roy, Kenneth J. Duncan, Huub Röttgering, George Miley, Catarina Aydar, Philip Best, Sarah E. I. Bosman, Alex J. Cameron, Krisztina Éva Gabányi, Andrew Humphrey, Sandy Morais, Masafusa Onoue, Laura Pentericci, Victoria Reynaldi, Bram Venemans

Abstract: We present rest-frame optical spectroscopy using JWST/NIRSpec IFU for the radio galaxy TN J1338-1942 at z=4.1, one of the most luminous galaxies in the early Universe with powerful extended radio jets. Previous observations showed evidence for strong, large-scale outflows on the basis of its large (~150 kpc) halo detected in Ly-alpha, and high velocity [O II] emission features detected in ground-b… ▽ More We present rest-frame optical spectroscopy using JWST/NIRSpec IFU for the radio galaxy TN J1338-1942 at z=4.1, one of the most luminous galaxies in the early Universe with powerful extended radio jets. Previous observations showed evidence for strong, large-scale outflows on the basis of its large (~150 kpc) halo detected in Ly-alpha, and high velocity [O II] emission features detected in ground-based IFU data. Our NIRSpec/IFU observations spatially resolve the emission line properties across the host galaxy in great detail. We find at least five concentrations of line emission, coinciding with discrete continuum features previously detected in imaging from HST and JWST, over an extent of ~2'' (~15 kpc). The spectral diagnostics enabled by NIRSpec unambiguously trace the activity of the obscured AGN plus interaction between the interstellar medium and the radio jet as the dominant mechanisms for the ionization state and kinematics of the gas in the system. A secondary region of very high ionization lies at roughly 5 kpc distance from the nucleus, and within the context of an expanding cocoon enveloping the radio lobe, this may be explained by strong shock-ionization of the entrained gas. However, it could also signal the presence of a second obscured AGN, which may also offer an explanation for an intriguing outflow feature seen perpendicular to the radio axis. The presence of a dual SMBH system in this galaxy would support that large galaxies in the early Universe quickly accumulated their mass through the merging of smaller units (each with their own SMBH), at the centers of large overdensities. The inferred black hole mass to stellar mass ratio of 0.01-0.1 for TNJ1338 points to a more rapid assembly of black holes compared to the stellar mass of galaxies at high redshifts, consistent with other recent observations. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 17 pages, 11 figures, submitted to MNRAS, comments welcome!

TREASUREHUNT: Transients and Variability Discovered with HST in the JWST North Ecliptic Pole Time Domain Field

Authors: Rosalia O'Brien, Rolf A. Jansen, Norman A. Grogin, Seth H. Cohen, Brent M. Smith, Ross M. Silver, W. P. Maksym III, Rogier A. Windhorst, Timothy Carleton, Anton M. Koekemoer, Nimish P. Hathi, Christopher N. A. Willmer, Brenda L. Frye, M. Alpaslan, M. L. N. Ashby, T. A. Ashcraft, S. Bonoli, W. Brisken, N. Cappelluti, F. Civano, C. J. Conselice, V. S. Dhillon, S. P. Driver, K. J. Duncan, R. Dupke , et al. (34 additional authors not shown)

Abstract: 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… ▽ More 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σ$ 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. △ Less

Submitted 2 May, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 32 pages, 11 figures, 5 tables, 1 Appendix

Galaxy quenching at the high redshift frontier: A fundamental test of cosmological models in the early universe with JWST-CEERS

Authors: Asa F. L. Bluck, Christopher J. Conselice, Katherine Ormerod, Joanna M. Piotrowska, Nathan Adams, Duncan Austin, Joseph Caruana, K. J. Duncan, Leonardo Ferreira, Paul Goubert, Thomas Harvey, James Trussler, Roberto Maiolino

Abstract: 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 For… ▽ More 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 ($φ_*$) 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érsic 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. △ Less

Submitted 4 November, 2023; originally announced November 2023.

Comments: Accepted to ApJ. 32 pages; 10 figures

Cosmology from LOFAR Two-metre Sky Survey Data Release 2: Angular Clustering of Radio Sources

Authors: C. L. Hale, D. J. Schwarz, P. N. Best, S. J. Nakoneczny, D. Alonso, D. Bacon, L. Böhme, N. Bhardwaj, M. Bilicki, S. Camera, C. S. Heneka, M. Pashapour-Ahmadabadi, P. Tiwari, J. Zheng, K. J. Duncan, M. J. Jarvis, R. Kondapally, M. Magliocchetti, H. J. A. Rottgering, T. W. Shimwell

Abstract: Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$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 wit… ▽ More Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$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θ<5\,°$, 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. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 29 pages, 24 figures

No strong radio absorption detected in the low-frequency spectra of radio-loud quasars at z > 5.6

Authors: A. J. Gloudemans, A. Saxena, H. Intema, J. R. Callingham, K. J. Duncan, H. J. A. Rottgering, S. Belladitta, M. J. Hardcastle, Y. Harikane, C. Spingola

Abstract: 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 brigh… ▽ More 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. △ Less

Submitted 19 September, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 15 pages, 7 figures, accepted for publication in A&A, minor correction Tab. 1

Journal ref: A&A 678, A128 (2023)

Plausible association of distant late M dwarfs with low-frequency radio emission

Authors: A. J. Gloudemans, J. R. Callingham, K. J. Duncan, A. Saxena, Y. Harikane, G. J. Hill, G. R. Zeimann, H. J. A. Rottgering, M. J. Hardcastle, J. S. Pineda, T. W. Shimwell, D. J. B. Smith, J. D. Wagenveld

Abstract: 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… ▽ More 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. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 10 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 678, A161 (2023)

The LOFAR Two-Metre Sky Survey (LoTSS): VI. Optical identifications for the second data release

Authors: M. J. Hardcastle, M. A. Horton, W. L. Williams, K. J. Duncan, L. Alegre, B. Barkus, J. H. Croston, H. Dickinson, E. Osinga, H. J. A. Röttgering, J. Sabater, T. W. Shimwell, D. J. B. Smith, P. N. Best, A. Botteon, M. Brüggen, A. Drabent, F. de Gasperin, G. Gürkan, M. Hajduk, C. L. Hale, M. Hoeft, M. Jamrozy, M. Kunert-Bajraszewska, R. Kondapally , et al. (27 additional authors not shown)

Abstract: 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 ext… ▽ More 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. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: 29 pages. Accepted by A&A; data products available at https://lofar-surveys.org/dr2_release.html

Journal ref: A&A 678, A151 (2023)

The LOFAR Two-metre Sky Survey: the radio view of the cosmic star formation history

Authors: R. K. Cochrane, R. Kondapally, P. N. Best, J. Sabater, K. J. Duncan, D. J. B. Smith, M. J. Hardcastle, H. J. A. Röttgering, I. Prandoni, P. Haskell, G. Gürkan, G. K. Miley

Abstract: 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ötes and the Lockman Hole, reaching $\sim20\,μ\rm{Jy/beam}$ rms sensi… ▽ More 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ötes and the Lockman Hole, reaching $\sim20\,μ\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. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 22 pages, 14 figures; accepted for publication in MNRAS

The LOFAR Two-metre Sky Survey Deep Fields Data Release 1: V. Survey description, source classifications and host galaxy properties

Authors: P. N. Best, R. Kondapally, W. L. Williams, R. K. Cochrane, K. J. Duncan, C. L. Hale, P. Haskell, K. Malek, I. McCheyne, D. J. B. Smith, L. Wang, A. Botteon, M. Bonato, M. Bondi, G. Calistro Rivera, F. Gao, G. Gurkan, M. J. Hardcastle, M. J. Jarvis, B. Mingo, H. Miraghaei, L. K. Morabito, D. Nisbet, I. Prandoni, H. J. A. Rottgering , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

Finding AGN remnant candidates based on radio morphology with machine learning

Authors: Rafael I. J. Mostert, Raffaella Morganti, Marisa Brienza, Kenneth J. Duncan, Martijn S. S. L. Oei, Huub J. A. Rottgering, Lara Alegre, Martin J. Hardcastle, Nika Jurlin

Abstract: Remnant radio galaxies represent the dying phase of radio-loud active galactic nuclei (AGN). Large samples of remnant radio galaxies are important for quantifying the radio galaxy life cycle. The remnants of radio-loud AGN can be identified in radio sky surveys based on their spectral index, or, complementary, through visual inspection based on their radio morphology. However, this is extremely ti… ▽ More Remnant radio galaxies represent the dying phase of radio-loud active galactic nuclei (AGN). Large samples of remnant radio galaxies are important for quantifying the radio galaxy life cycle. The remnants of radio-loud AGN can be identified in radio sky surveys based on their spectral index, or, complementary, through visual inspection based on their radio morphology. However, this is extremely time-consuming when applied to the new large and sensitive radio surveys. Here we aim to reduce the amount of visual inspection required to find AGN remnants based on their morphology, through supervised machine learning trained on an existing sample of remnant candidates. For a dataset of 4107 radio sources, with angular sizes larger than 60 arcsec, from the LOw Frequency ARray (LOFAR) Two-Metre Sky Survey second data release (LoTSS-DR2), we started with 151 radio sources that were visually classified as 'AGN remnant candidate'. We derived a wide range of morphological features for all radio sources from their corresponding Stokes-I images: from simple source catalogue-derived properties, to clustered Haralick-features, and self-organising map (SOM) derived morphological features. We trained a random forest classifier to separate the 'AGN remnant candidates' from the not yet inspected sources. The SOM-derived features and the total to peak flux ratio of a source are shown to be most salient to the classifier. We estimate that $31\pm5\%$ of sources with positive predictions from our classifier will be labelled 'AGN remnant candidates' upon visual inspection, while we estimate the upper bound of the $95\%$ confidence interval for 'AGN remnant candidates' in the negative predictions at $8\%$. Visual inspection of just the positive predictions reduces the number of radio sources requiring visual inspection by $73\%$. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 23 pages; accepted for publication in A&A

Journal ref: A&A 674, A208 (2023)

JWST's PEARLS: TN J1338-1942 -- I. Extreme jet triggered star-formation in a $z=4.11$ luminous radio galaxy

Authors: Kenneth J. Duncan, Rogier A. Windhorst, Anton M. Koekemoer, Huub J. A. Röttgering, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Scott Tompkins, Taylor A. Hutchison, Christopher J. Conselice, Simon P. Driver, Haojing Yan, Nathan J. Adams, Cheng Cheng, Dan Coe, Jose M. Diego, Hervé Dole, Brenda Frye, Hansung B. Gim, Norman A. Grogin, Benne W. Holwerda, Jeremy Lim, Madeline A. Marshall, Mario Nonino, Nor Pirzkal , et al. (3 additional authors not shown)

Abstract: We present the first JWST observations of the $z=4.11$ luminous radio galaxy TN J1338-1942, obtained as part of the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') project. Our NIRCam observations, designed to probe the key rest-frame optical continuum and emission line features at this redshift, enable resolved spectral energy distribution modelling that incorporat… ▽ More We present the first JWST observations of the $z=4.11$ luminous radio galaxy TN J1338-1942, obtained as part of the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') project. Our NIRCam observations, designed to probe the key rest-frame optical continuum and emission line features at this redshift, enable resolved spectral energy distribution modelling that incorporates both a range of stellar population assumptions and radiative shock models. With an estimated stellar mass of $\log_{10}(M/\text{M}_{\odot}) \sim 10.9$, TN J1338--1942 is confirmed to be one of the most massive galaxies known at this epoch. Our observations also reveal extremely high equivalent-width nebular emission coincident with the luminous AGN jets that is best fit by radiative shocks surrounded by extensive recent star-formation. We estimate the total star-formation rate (SFR) could be as high as $\sim1600\,\text{M}_{\odot}\,\text{yr}^{-1}$, with the SFR that we attribute to the jet induced burst conservatively $\gtrsim500\,\text{M}_{\odot}\,\text{yr}^{-1}$. The mass-weighted age of the star-formation, $t_{\text{mass}} <4$ Myr, is consistent with the likely age of the jets responsible for the triggered activity and significantly younger than that measured in the core of the host galaxy. The extreme scale of the potential jet-triggered star-formation activity indicates the potential importance of positive AGN feedback in the earliest stages of massive galaxy formation, with our observations also illustrating the extraordinary prospects for detailed studies of high-redshift galaxies with JWST. △ Less

Submitted 25 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 17 pages, 15 figures, Accepted for publication in MNRAS after minor revisions

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

Extragalactic Peaked-Spectrum Radio Sources at Low-Frequencies are Young Radio Galaxies

Authors: M. M. Slob, J. R. Callingham, H. J. A. Röttgering, W. L. Williams, K. J. Duncan, F. de Gasperin, M. J. Hardcastle, G. K. Miley

Abstract: 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 a… ▽ More 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. △ Less

Submitted 29 October, 2022; originally announced October 2022.

Comments: 21 pages, accepted in Astronomy & Astrophysics

Journal ref: A&A 668, A186 (2022)

JWST's PEARLS: A JWST/NIRCam view of ALMA sources

Authors: Cheng Cheng, Jia-Sheng Huang, Ian Smail, Haojing Yan, Seth H. Cohen, Rolf A. Jansen, Rogier A. Windhorst, Zhiyuan Ma, Anton Koekemoer, Christopher N. A. Willmer, S. P. Willner, Jose M. Diego, Brenda Frye, Christopher J. Conselice, Leonardo Ferreira, Andreea Petric, Min Yun, Hansung B. Gim, Maria del Carmen Polletta, Kenneth J. Duncan, Rachel Honor, Benne W. Holwerda, Huub J. A. Röttgering, Nimish P. Hathi, Patrick S. Kamieneski , et al. (11 additional authors not shown)

Abstract: We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its k… ▽ More We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts' spectral energy distributions from rest-frame ultraviolet to near infrared provide photometric redshifts ($1<z<4.5$) and stellar masses ($M_*>10^{10.5}$ Msol), which are similar to sub-millimeter galaxy (SMG) hosts studied previously. However, our sample is fainter in submm/mm than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star-formation rates as low as 10 Msol yr$^{-1}$, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few hosts have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question. △ Less

Submitted 9 December, 2022; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 19 pages, 11 figures, accepted by ApJL

Discovery of 24 radio-bright quasars at $4.9 \leq z \leq6.6$ using low-frequency radio observations

Authors: A. J. Gloudemans, K. J. Duncan, A. Saxena, Y. Harikane, G. J. Hill, G. R. Zeimann, H. J. A. Rottgering, D. Yang, P. N. Best, E. Banados, A. Drabent, M. J. Hardcastle, J. F. Hennawi, G. Lansbury, M. Magliocchetti, G. K. Miley, R. Nanni, T. W. Shimwell, D. J. B. Smith, B. P. Venemans, J. D. Wagenveld

Abstract: 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 Spectro… ▽ More 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 $μ$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$α$ 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. △ Less

Submitted 21 October, 2022; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 19 pages, 11 figures, Tab A1 corrected, accepted for publication in A&A

Radio source-component association for the LOFAR Two-metre Sky Survey with region-based convolutional neural networks

Authors: Rafaël I. J. Mostert, Kenneth J. Duncan, Lara Alegre, Huub J. A. Röttgering, Wendy L. Williams, Philip N. Best, Martin J. Hardcastle, Raffaella Morganti

Abstract: Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sk… ▽ More Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) at $144$ MHz, is a daunting, time-consuming process, even with extensive manpower. Using a machine learning approach, we aim to automate the radio component association of large ($> 15$ arcsec) radio components. We turned the association problem into a classification problem and trained an adapted Fast region-based convolutional neural network to mimic the expert annotations from the first LoTSS data release. We implemented a rotation data augmentation to reduce overfitting and simplify the component association by removing unresolved radio sources that are likely unrelated to the large and bright radio components that we consider using predictions from an existing gradient boosting classifier. For large ($> 15$ arcsec) and bright ($> 10$ mJy) radio components in the LoTSS first data release, our model provides the same associations for $85.3\%\pm0.6$ of the cases as those derived when astronomers perform the association manually. When the association is done through public crowd-sourced efforts, a result similar to that of our model is attained. Our method is able to efficiently carry out manual radio-component association for huge radio surveys and can serve as a basis for either automated radio morphology classification or automated optical host identification. This opens up an avenue to study the completeness and reliability of samples of radio sources with extended, complex morphologies. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 22 pages; accepted for publication in A&A

Journal ref: A&A 668, A28 (2022)

JWST's PEARLS: Prime Extragalactic Areas for Reionization and Lensing Science: Project Overview and First Results

Authors: Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Scott Tompkins, Christopher J. Conselice, Simon P. Driver, Haojing Yan, Dan Coe, Brenda Frye, Norman Grogin, Anton Koekemoer, Madeline A. Marshall, Rosalia O'Brien, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Christopher N. A. Willmer, Timothy Carleton, Jose M. Diego, William C. Keel, Paolo Porto, Caleb Redshaw, Sydney Scheller, Stephen M. Wilkins , et al. (60 additional authors not shown)

Abstract: We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters;… ▽ More We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $μ$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 μm. PEARLS is designed to be of lasting benefit to the community. △ Less

Submitted 28 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

Identifying active galactic nuclei via brightness temperature with sub-arcsecond International LOFAR Telescope observations

Authors: Leah K. Morabito, F. Sweijen, J. F. Radcliffe, P. N. Best, Rohit Kondapally, Marco Bondi, Matteo Bonato, K. J. Duncan, Isabella Prandoni, T. W. Shimwell, W. L. Williams, R. J. van Weeren, J. E. Conway, G. Calistro Rivera

Abstract: 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 LOF… ▽ More 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. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS. 19 pages, 15 figures, 2 tables

Hybrid photometric redshifts for sources in the COSMOS and XMM-LSS fields

Authors: P. W. Hatfield, M. J. Jarvis, N. Adams, R. A. A. Bowler, B. Häußler, K. J. Duncan

Abstract: 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 p… ▽ More 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. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 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

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 3, pp.3719-3733, 2022

Cosmic evolution of low-excitation radio galaxies in the LOFAR Two-meter Sky Survey Deep Fields

Authors: R. Kondapally, P. N. Best, R. K. Cochrane, J. Sabater, K. J. Duncan, M. J. Hardcastle, P. Haskell, B. Mingo, H. J. A. Röttgering, D. J. B. Smith, W. L. Williams, M. Bonato, G. Calistro Rivera, F. Gao, C. L. Hale, K. Małek, G. K. Miley, I. Prandoni, L. Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 22 April, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 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

All-purpose, all-sky photometric redshifts for the Legacy Imaging Surveys Data Release 8

Authors: Kenneth J. Duncan

Abstract: In this paper we present photometric redshift (photo-$z$) estimates for the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, currently the most sensitive optical survey covering the majority of the extra-galactic sky. Our photo-$z$ methodology is based on a machine-learning approach, using sparse Gaussian processes augmented with Gaussian mixture models (GMMs) that allow regions… ▽ More In this paper we present photometric redshift (photo-$z$) estimates for the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, currently the most sensitive optical survey covering the majority of the extra-galactic sky. Our photo-$z$ methodology is based on a machine-learning approach, using sparse Gaussian processes augmented with Gaussian mixture models (GMMs) that allow regions of parameter space to be identified and trained separately in a purely data-driven way. The same GMMs are also used to calculate cost-sensitive learning weights that mitigate biases in the spectroscopic training sample. By design, this approach aims to produce reliable and unbiased predictions for all parts of the parameter space present in wide area surveys. Compared to previous literature estimates using the same underlying photometry, our photo-$z$s are significantly less biased and more accurate at $z > 1$, with negligible loss in precision or reliability for resolved galaxies at $z < 1$. Our photo-$z$ estimates offer accurate predictions for rare high-value populations within the parent sample, including optically selected quasars at the highest redshifts ($z > 6$), as well as X-ray or radio continuum selected populations across a broad range of flux (densities) and redshift. Deriving photo-$z$ estimates for the full Legacy Imaging Surveys Data Release 8, the catalogues provided in this work offer photo-$z$ estimates predicted to be high quality for $\gtrsim9\times10^{8}$ galaxies over $\sim 19\,400\,\text{deg}^{2}$ and spanning $0 < z \lesssim 7$, offering one of the most extensive samples of redshift estimates ever produced. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 22 pages, 19 figures - Accepted for publication in MNRAS. Catalogues produced will be made available through queryable public databases - users interested in the full catalogues or early access to subsets are also encouraged to contact the author directly

Revealing new high redshift quasar populations through Gaussian mixture model selection

Authors: J. D. Wagenveld, A. Saxena, K. J. Duncan, H. J. A. Röttgering, M. Zhang

Abstract: 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 poste… ▽ More 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. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 15 pages, 8 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A22 (2022)

Accretion mode versus radio morphology in the LOFAR Deep Fields

Authors: B. Mingo, J. H. Croston, P. N. Best, K. J. Duncan, M. J. Hardcastle, R. Kondapally, I. Prandoni, J. Sabater, T. W. Shimwell, W. L. Williams, R. D. Baldi, M. Bonato, M. Bondi, P. Dabhade, G. Gürkan, J. Ineson, M. Magliocchetti, G. Miley, J. C. S. Pierce, H. J. A. Röttgering

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 22 pages, 17 figures, 7 tables. Accepted for publication in MNRAS. Update with minor changes to match published version

Low frequency radio properties of the $z>5$ quasar population

Authors: A. J. Gloudemans, K. J. Duncan, H. J. A. Röttgering, T. W. Shimwell, B. P. Venemans, P. N. Best, M. Brüggen, G. Calistro Rivera, A. Drabent, M. J. Hardcastle, G. K. Miley, D. J. Schwarz, A. Saxena, D. J. B. Smith, W. L. Williams

Abstract: 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… ▽ More 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 $μ$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 σ$ 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. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 11 pages, 10 figures, accepted for publication in A&A

The radio loudness of SDSS quasars from the LOFAR Two-metre Sky Survey: ubiquitous jet activity and constraints on star formation

Authors: C. Macfarlane, P. N. Best, J. Sabater, G. Gurkan, M. J. Jarvis, H. J. A. Rottgering, R. D. Baldi, G. Calistro Rivera, K. J. Duncan, L. K. Morabito, I. Prandoni, E. Retana-Montenegro

Abstract: 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… ▽ More 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}^α(1+z)^β$, we find $α= 0.47 \pm 0.01$ and $β= 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. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS. 21 pages

Low-frequency radio spectra of submillimetre galaxies in the Lockman Hole

Authors: J. Ramasawmy, J. E. Geach, M. J. Hardcastle, P. N. Best, M. Bonato, M. Bondi, G. Calistro Rivera, R. K. Cochrane, J. E. Conway, K. Coppin, K. J. Duncan, J. S. Dunlop, M. Franco, C. García-Vergara, M. J. Jarvis, R. Kondapally, I. McCheyne, I. Prandoni, H. J. A. Röttgering, D. J. B. Smith, C. Tasse, L. Wang

Abstract: We investigate the radio properties of a sample of 53 sources selected at 850 $μ$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… ▽ More We investigate the radio properties of a sample of 53 sources selected at 850 $μ$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. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 15 pages, 10 figures, accepted by A&A

Journal ref: A&A 648, A14 (2021)

A hyperluminous obscured quasar at a redshift of z ~ 4.3

Authors: Andreas Efstathiou, Katarzyna Malek, Denis Burgarella, Peter Hurley, Seb Oliver, Veronique Buat, Raphael Shirley, Steven Duivenvoorden, Vicky Papadopoulou Lesta, Duncan Farrah, Kenneth J. Duncan, Maria del Carmen Campos Varillas

Abstract: In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies a… ▽ More In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$. △ Less

Submitted 8 April, 2022; v1 submitted 14 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. The CYGNUS models used in this work are available at https://arc.euc.ac.cy/cygnus-project-arc/

Journal ref: 2021MNRAS503L..11E

The LOFAR Two-metre Sky Survey Deep Fields -- Data Release 1: IV. Photometric redshifts and stellar masses

Authors: K. J. Duncan, R. Kondapally, M. J. I. Brown, M. Bonato, P. N. Best, H. J. A. Röttgering, M. Bondi, R. A. A. Bowler, R. K. Cochrane, G. Gürkan, M. J. Hardcastle, M. J. Jarvis, M. Kunert-Bajraszewska, S. K. Leslie, K. Małek, L. K. Morabito, S. P. O'Sullivan, I. Prandoni, J. Sabater, T. W. Shimwell, D. J. B. Smith, L. Wang, A. Wołowska

Abstract: 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ötes, European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1), and Lockman Hole fields, down to rms sensitivities of… ▽ More 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ö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 $μ$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. △ Less

Submitted 19 November, 2020; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: 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

Journal ref: A&A 648, A4 (2021)

LOFAR properties of SILVERRUSH Ly$α$ emitter candidates in the ELAIS-N1 field

Authors: A. J. Gloudemans, K. J. Duncan, R. Kondapally, J. Sabater, R. K. Cochrane, H. J. A. Röttgering, P. N. Best, M. Bonato, M. Bondi, K. Malek, I. McCheyne, D. J. B. Smith, I. Prandoni, L. Wang

Abstract: 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 Arr… ▽ More 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\,μ$Jy beam$^{-1}$ at 150MHz, and we detect five candidate LAEs at $>5σ$ 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$σ$ 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}$. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 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 LOFAR Two Metre Sky Survey: Deep Fields Data Release 1 -- III. Host-galaxy identifications and value added catalogues

Authors: R. Kondapally, P. N. Best, M. J. Hardcastle, D. Nisbet, M. Bonato, J. Sabater, K. J. Duncan, I. McCheyne, R. K. Cochrane, R. A. A. Bowler, W. L. Williams, T. W. Shimwell, C. Tasse, J. H. Croston, A. Goyal, M. Jamrozy, M. J. Jarvis, V. H. Mahatma, H. J. A. R\öttgering, D. J. B. Smith, A. Wo\łowska, M. Bondi, M. Brienza, M. J. I. Brown, M. Br\üggen , et al. (18 additional authors not shown)

Abstract: 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ötes fields, down to RMS sensitives of around 20, 22, and 32… ▽ More 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ötes fields, down to RMS sensitives of around 20, 22, and 32$~μ$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ö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. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 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

Journal ref: A&A 648, A3 (2021)

The LOFAR Two-metre Sky Survey Deep fields: The star formation rate - radio luminosity relation at low frequencies

Authors: D. J. B. Smith, P. Haskell, G. Gürkan, P. N. Best, M. J. Hardcastle, R. Kondapally, W. Williams, K. J. Duncan, R. K. Cochrane, I. McCheyne, H. J. A. Röttgering, J. Sabater, T. W. Shimwell, C. Tasse, M. Bonato, M. Bondi, M. J. Jarvis, S. K. Leslie, I. Prandoni, L. Wang

Abstract: 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, a… ▽ More 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}(ψ/M_\odot\,\mathrm{yr}^{-1}) + (0.33 \pm 0.04) \log_{10} (M/10^{10}M_\odot) + 22.22 \pm 0.02$. (Abridged) △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 17 pages, 17 figures, accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields

Journal ref: A&A 648, A6 (2021)

Unveiling the rarest morphologies of the LOFAR Two-metre Sky Survey radio source population with self-organised maps

Authors: Rafaël I. J. Mostert, Kenneth J. Duncan, Huub J. A. Röttgering, Kai L. Polsterer, Philip N. Best, Marisa Brienza, Marcus Brüggen, Martin J. Hardcastle, Nika Jurlin, Beatriz Mingo, Raffaella Morganti, Tim Shimwell, Dan Smith, Wendy L. Williams

Abstract: The Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) is a low-frequency radio continuum survey of the Northern sky at an unparalleled resolution and sensitivity. In order to fully exploit this huge dataset and those produced by the Square Kilometre Array in the next decade, automated methods in machine learning and data-mining will be increasingly essential both for morphological classific… ▽ More The Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) is a low-frequency radio continuum survey of the Northern sky at an unparalleled resolution and sensitivity. In order to fully exploit this huge dataset and those produced by the Square Kilometre Array in the next decade, automated methods in machine learning and data-mining will be increasingly essential both for morphological classifications and for identifying optical counterparts to the radio sources. Using self-organising maps (SOMs), a form of unsupervised machine learning, we created a dimensionality reduction of the radio morphologies for the $\sim$25k extended radio continuum sources in the LoTSS first data release, which is only $\sim$2 percent of the final LoTSS survey. We made use of \textsc{PINK}, a code which extends the SOM algorithm with rotation and flipping invariance, increasing its suitability and effectiveness for training on astronomical sources. After training, the SOMs can be used for a wide range of science exploitation and we present an illustration of their potential by finding an arbitrary number of morphologically rare sources in our training data (424 square degrees) and subsequently in an area of the sky ($\sim$5300 square degrees) outside the training data. Objects found in this way span a wide range of morphological and physical categories: extended jets of radio active galactic nuclei, diffuse cluster haloes and relics, and nearby spiral galaxies. Finally, to enable accessible, interactive, and intuitive data exploration, we showcase the LOFAR-PyBDSF Visualisation Tool, which allows users to explore the LoTSS dataset through the trained SOMs. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 26 pages; accepted for publication in A&A

Journal ref: A&A 645, A89 (2021)

The Chandra Deep Wide-Field Survey: A New Chandra Legacy Survey in the Boötes Field I. X-ray Point Source Catalog, Number Counts and Multi-Wavelength Counterparts

Authors: A. Masini, R. C. Hickox, C. M. Carroll, J. Aird, D. M. Alexander, R. J. Assef, R. Bower, M. Brodwin, M. J. I. Brown, S. Chatterjee, C. -T. J. Chen, A. Dey, M. A. DiPompeo, K. J. Duncan, P. R. M. Eisenhardt, W. R. Forman, A. H. Gonzalez, A. D. Goulding, K. N. Hainline, B. T. Jannuzi, C. Jones, C. S. Kochanek, R. Kraft, K. -S. Lee, E. D. Miller , et al. (8 additional authors not shown)

Abstract: We present a new, ambitious survey performed with the Chandra X-ray Observatory of the 9.3 deg$^2$ Boötes 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… ▽ More We present a new, ambitious survey performed with the Chandra X-ray Observatory of the 9.3 deg$^2$ Boötes 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. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJS. The catalogs associated to this paper can be accessed at http://www.dartmouth.edu/~hickox/cdwfs.php

Alignment in the orientation of LOFAR radio sources

Authors: E. Osinga, G. K. Miley, R. J. van Weeren, T. W. Shimwell, K. J. Duncan, M. J. Hardcastle, A. P. Mechev, H. J. A. Röttgering, C. Tasse, W. L. Williams

Abstract: 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 sk… ▽ More 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. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 10 pages, 13 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 642, A70 (2020)

The MOSDEF Survey: Calibrating the relationship between H$α$ star-formation rate and radio continuum luminosity at $1.4 < z < 2.6$

Authors: Kenneth J. Duncan, Irene Shivaei, Alice E. Shapley, Naveen A. Reddy, Bahram Mobasher, Alison L. Coil, Mariska Kriek, Brian Siana

Abstract: The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$α$ SFR. Our sample con… ▽ More The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$α$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$α$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$α$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 11 pages, 5 figures. Submitted for publication in MNRAS - comments welcome

The life cycle of radio galaxies in the LOFAR Lockman Hole field

Authors: N. Jurlin, R. Morganti, M. Brienza, S. Mandal, N. Maddox, K. J. Duncan, S. S. Shabala, M. J. Hardcastle, I. Prandoni, H. J. A. Röttgering, V. Mahatma, P. N. Best, B. Mingo, J. Sabater, T. W. Shimwell, C. Tasse

Abstract: 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 fi… ▽ More 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. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures, 4 tables, A&A accepted, companion paper to Shabala et al., also in today arXiv posting

Journal ref: A&A 638, A34 (2020)

VLT/SINFONI study of black hole growth in high redshift radio-loud quasars from the CARLA survey

Authors: M. Marinello, R. A. Overzier, H. J. A. Röttgering, J. D. Kurk, C. De Breuck, J. Vernet, D. Wylezalek, D. Stern, K. J. Duncan, N. Hatch, N. Kashikawa, Y. -T. Lin, R. S. Nemmen, A. Saxena

Abstract: 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 co… ▽ More 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. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: MNRAS, accepted (28 pages, 22 figures)

A LOFAR-IRAS cross-match study: the far-infrared radio correlation and the 150-MHz luminosity as a star-formation rate

Authors: L. Wang, F. Gao, K. J. Duncan, W. L. Williams, M. Rowan-Robinson, J. Sabater, T. W. Shimwell, M. Bonato, G. Calistro-Rivera, K. T. Chyzy, D. Farrah, G. Gurkan, M. J. Hardcastle, I. McCheyne, I. Prandoni, S. C. Read, H. J. A. Rottgering, D. J. B. Smith

Abstract: 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… ▽ More 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. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: 13 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 631, A109 (2019)

HELP: A catalogue of 170 million objects, selected at 0.36-4.5 $μ$m, from 1270 deg.$^{2}$ of prime extragalactic fields

Authors: Raphael Shirley, Yannick Roehlly, Peter D Hurley, Veronique Buat, María del Carmen Campos Varillas, Steven Duivenvoorden, Kenneth J Duncan, Andreas Efstathiou, Duncan Farrah, Eduardo González Solares, Katarzyna Małek, Lucia Marchetti, Ian McCheyne, Andreas Papadopoulos, Estelle Pons, Roberto Scipioni, Mattia Vaccari, Seb Oliver

Abstract: We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170… ▽ More We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5σ$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky. △ Less

Submitted 22 October, 2019; v1 submitted 9 September, 2019; originally announced September 2019.

Journal ref: MNRAS 490.1 (2019): 634-656

One- and Two-point Source Statistics from the LOFAR Two-metre Sky Survey First Data Release

Authors: Thilo M. Siewert, Catherine Hale, Nitesh Bhardwaj, Marian Biermann, David J. Bacon, Matt Jarvis, Huub Röttgering, Dominik J. Schwarz, Timothy Shimwell, Philip N. Best, Kenneth J. Duncan, Martin J. Hardcastle, Jose Sabater, Cyril Tasse, Glenn J. White, Wendy L. Williams

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-cell statistics and differential source counts statistic and measure the angular two-point correlation function of the LoTSS radio sources. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies and with simulated catalogues from the SKA design study sky and the Tiered Radio Extragalactic Continuum Simulation. The angular two-point correlation is $<10^{-2}$ at angular scales $> 1$ deg. Restricting the value added source catalogue to low-noise regions and a flux density threshold of 2 mJy provides our most reliable estimate of the angular two-point correlation. For smaller flux density thresholds systematic issues are identified, most likely related to the flux density calibration of the individual pointings. Based on the distribution of photometric redshifts of LoTSS sources and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 6 deg agrees with the measured clustering for a subsample of radio sources with redshift information. We find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level. The angular two-point correlation agrees well with the expectation of the cosmological standard model. (abbreviated) △ Less

Submitted 3 September, 2020; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 30 pages, 35 figures; corrections of affiliations and acknowledgements; version accepted by A&A

Journal ref: A&A 643, A100 (2020)

The Spectral Energy Distributions of Active Galactic Nuclei

Authors: M. J. I. Brown, K. J. Duncan, H. Landt, M. Kirk, C. Ricci, N. Kamraj, M. Salvato, T. Ananna

Abstract: 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 mea… ▽ More 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. △ Less

Submitted 10 August, 2019; originally announced August 2019.

Comments: 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

Revisiting the Fanaroff-Riley dichotomy and radio-galaxy morphology with the LOFAR Two-Metre Sky Survey (LoTSS)

Authors: B. Mingo, J. H. Croston, M. J. Hardcastle, P. N. Best, K. J. Duncan, R. Morganti, H. J. A. Rottgering, J. Sabater, T. W. Shimwell, W. L. Williams, M. Brienza, G. Gurkan, V. H. Mahatma, L. K. Morabito, I. Prandoni, M. Bondi, J. Ineson, S. Mooney

Abstract: 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 o… ▽ More 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. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: 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

The nature of faint radio galaxies at high redshifts

Authors: A. Saxena, H. J. A. Röttgering, K. J. Duncan, G. J. Hill, P. N. Best, B. L. Indahl, M. Marinello, R. A. Overzier, L. Pentericci, I. Prandoni, H. Dannerbauer, R. Barrena

Abstract: 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 redsh… ▽ More 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$α$ 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. △ Less

Submitted 23 September, 2019; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: 26 pages (including 2 appendices), 17 figures, accepted for publication in MNRAS

Giant radio galaxies in the LOFAR Two-metre Sky Survey-I

Authors: P. Dabhade, H. J. A. Rottgering, J. Bagchi, T. W. Shimwell, M. J. Hardcastle, S. Sankhyayan, R. Morganti, M. Jamrozy, A. Shulevski, K. J. Duncan

Abstract: 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 sq… ▽ More 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. △ Less

Submitted 4 December, 2019; v1 submitted 31 March, 2019; originally announced April 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 635, A5 (2020)

The Spectral Energy Distributions of Active Galactic Nuclei

Authors: M. J. I. Brown, K. J. Duncan, H. Landt, M. Kirk, C. Ricci, N. Kamraj

Abstract: We present ongoing work on the spectral energy distributions (SEDs) of active galactic nuclei (AGNs), derived from X-ray, ultraviolet, optical, infrared and radio photometry and spectroscopy. Our work is motivated by new wide-field imaging surveys that will identify vast numbers of AGNs, and by the need to benchmark AGN SED fitting codes. We have constructed 41 SEDs of individual AGNs and 80 addit… ▽ More We present ongoing work on the spectral energy distributions (SEDs) of active galactic nuclei (AGNs), derived from X-ray, ultraviolet, optical, infrared and radio photometry and spectroscopy. Our work is motivated by new wide-field imaging surveys that will identify vast numbers of AGNs, and by the need to benchmark AGN SED fitting codes. We have constructed 41 SEDs of individual AGNs and 80 additional SEDs that mimic Seyfert spectra. All of our SEDs span 0.09 to 30 microns, while some extend into the X-ray and/or radio. We have tested the utility of the SEDs by using them to generate AGN photometric redshifts, and they outperform SEDs from the prior literature, including reduced redshift errors and flux density residuals. △ Less

Submitted 18 February, 2019; originally announced February 2019.

Comments: 5 pages, 3 figures. To be published in the Proceedings of the IAU symposium 341, Osaka, Japan, 2018

Journal ref: Proc. IAU 15 (2019) 21-25