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

Galaxy Mergers in the Epoch of Reionization I: A JWST Study of Pair Fractions, Merger Rates, and Stellar Mass Accretion Rates at $z = 4.5-11.5$

Authors: Qiao Duan, Christopher J. Conselice, Qiong Li, Duncan Austin, Thomas Harvey, Nathan J. Adams, Kenneth J. Duncan, James Trussler, Leonardo Ferreira, Lewi Westcott, Honor Harris, Rogier A. Windhorst, Benne W. Holwerda, Thomas J. Broadhurst, Dan Coe, Seth H. Cohen, Simon P. Driver, Brenda Frye, Norman A. Grogin, Nimish P. Hathi, Rolf A. Jansen, Anton M. Koekemoer, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III , et al. (7 additional authors not shown)

Abstract: We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLA… ▽ More We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLASS, El-Gordo, SMACS-0723, MACS-0416), covering an unmasked area of 189.36 $\mathrm{arcmin}^2$. We develop a new probabilistic pair-counting methodology that integrates full photometric redshift posteriors and corrects for detection incompleteness to quantify close pairs with physical projected separations between 20 and 50 kpc. Our analysis reveals an increase in pair fractions up to $z = 8$, reaching $0.211 \pm 0.065$, followed by a statistically flat evolution to $z = 11.5$. We find that the galaxy merger rate increases from the local Universe up to $z = 6$ and then stabilizes at a value of $\sim 6$ Gyr$^{-1}$ up to $z = 11.5$. We fit both a power-law and a power-law + exponential model to our pair fraction and merger rate redshift evolution, finding that the latter model describes the trends more accurately, particularly at $z = 8.0 - 11.5$. In addition, we measure that the average galaxy increases its stellar mass due to mergers by a factor of $2.77 \pm 0.99$ from redshift $z = 10.5$ to $z = 5.0$. Lastly, we investigate the impact of mergers on galaxy stellar mass growth, revealing that mergers contribute $71 \pm 25\%$ as much to galaxy stellar mass increases as star formation from gas. This indicates that mergers drive about half of galaxy assembly at high redshift. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 30 Pages, 15 Figures, Submitted to MNRAS

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!

JWST Reveals Powerful Feedback from Radio Jets in a Massive Galaxy at z = 4.1

Authors: Namrata Roy, Timothy Heckman, Roderik Overzier, Aayush Saxena, Kenneth Duncan, George Miley, Montserrat Villar Martín, Krisztina Éva Gabányi, Catarina Aydar, Sarah E. I. Bosman, Huub Rottgering, Laura Pentericci, Masafusa Onoue, Victoria Reynaldi

Abstract: We report observations of a powerful ionized gas outflow in a z = 4.1 luminous ($ L_{1.4GHz} \sim 10^{28.3} \ W \ Hz^{-1}$) radio galaxy TNJ1338-1942 hosting an obscured quasar using the Near Infrared Spectrograph (NIRSpec) on board JWST. We spatially resolve a large-scale (~15 kpc) outflow and measure resolved outflow rates. The outflowing gas shows velocities exceeding 900 $ km \ s^{-1}$ and bro… ▽ More We report observations of a powerful ionized gas outflow in a z = 4.1 luminous ($ L_{1.4GHz} \sim 10^{28.3} \ W \ Hz^{-1}$) radio galaxy TNJ1338-1942 hosting an obscured quasar using the Near Infrared Spectrograph (NIRSpec) on board JWST. We spatially resolve a large-scale (~15 kpc) outflow and measure resolved outflow rates. The outflowing gas shows velocities exceeding 900 $ km \ s^{-1}$ and broad line profiles with line widths exceeding 1200 $ km \ s^{-1}$ located at ~10 kpc projected distance from the central nucleus. The outflowing nebula spatially overlaps with the brightest radio lobe, indicating that the powerful radio jets are responsible for the extraordinary kinematics exhibited by the ionized gas. The ionized gas is possibly ionized by the central obscured quasar with a contribution from shocks. The spatially resolved mass outflow rate shows that the region with the broadest line profiles exhibits the strongest outflow rates, with an integrated mass outflow rate of ~500 $ M_{\odot} \ yr^{-1}$. Our hypothesis is that an over-pressured shocked jet fluid expands laterally to create an expanding ellipsoidal "cocoon" that causes the surrounding gas to accelerate outwards. The total kinetic energy injected by the radio jet is about 3 orders of magnitude larger than the total kinetic energy measured in the outflowing ionized gas. This implies that kinetic energy must be transferred inefficiently from the jets to the gas. The bulk of the deposited energy possibly lies in the form of hot (~$ 10^7$ K) X-ray-emitting gas. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: 24 pages, 15 figures, 2 tables. Under review in ApJ

Exploring the radio-loudness of SDSS quasars with spectral stacking

Authors: M. I. Arnaudova, D. J. B. Smith, M. J. Hardcastle, S. Das, A. Drake, K. Duncan, G. Gürkan, M. Magliocchetti, L. K. Morabito, J. W. Petley, S. Shenoy, C. Tasse

Abstract: We use new 144 MHz observations over 5634 deg$^2$ from the LOFAR Two-metre Sky Survey (LoTSS) to compile the largest sample of uniformly-selected, spectroscopically-confirmed quasars from the 14th data release of the Sloan Digital Sky Survey (SDSS-DR14). Using the classical definition of radio-loudness, $R=\log(L_{\rm{1.4GHz}}/L_{i})$, we identify 3,697 radio-loud (RL) and 111,132 radio-quiet (RQ)… ▽ More We use new 144 MHz observations over 5634 deg$^2$ from the LOFAR Two-metre Sky Survey (LoTSS) to compile the largest sample of uniformly-selected, spectroscopically-confirmed quasars from the 14th data release of the Sloan Digital Sky Survey (SDSS-DR14). Using the classical definition of radio-loudness, $R=\log(L_{\rm{1.4GHz}}/L_{i})$, we identify 3,697 radio-loud (RL) and 111,132 radio-quiet (RQ) sources at $0.6<z<3.4$. To study their properties, we develop a new rest-frame spectral stacking algorithm, designed with forthcoming massively-multiplexed spectroscopic surveys in mind, and use it to create high signal-to-noise composite spectra of each class, matched in redshift and absolute $i$-band magnitude. We show that RL quasars have redder continuum and enhanced [OII] emission than their RQ counterparts. These results persist when additionally matching in black hole mass, suggesting that this parameter is not the defining factor in making a QSO radio-loud. We find that these features are not gradually varying as a function of radio-loudness but are maintained even when probing deeper into the RQ population, indicating that a clear-cut division in radio-loudness is not apparent. Upon examining the star formation rates (SFRs) inferred from the [OII] emission line, with the contribution from AGN removed using the [NeV] line, we find that RL quasars have a significant excess of star-formation relative to RQ quasars out to $z=1.9$ at least. Given our findings, we suggest that radio-loud sources either preferably reside in gas-rich systems with rapidly-spinning black holes, or represent an earlier obscured phase of QSO evolution. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

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)

The Faraday Rotation Measure Grid of the LOFAR Two-metre Sky Survey: Data Release 2

Authors: S. P. O'Sullivan, T. W. Shimwell, M. J. Hardcastle, C. Tasse, G. Heald, E. Carretti, M. Brüggen, V. Vacca, C. Sobey, C. L. Van Eck, C. Horellou, R. Beck, M. Bilicki, S. Bourke, A. Botteon, J. H. Croston, A. Drabent, K. Duncan, V. Heesen, S. Ideguchi, M. Kirwan, L. Lawlor, B. Mingo, B. Nikiel-Wroczyński, J. Piotrowska , et al. (2 additional authors not shown)

Abstract: A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polar… ▽ More A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8σ_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 20 pages, 23 figures, 3 tables. Accepted for publication in MNRAS. Go to https://lofar-mksp.org/data/ to download the catalogue and associated data products

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

Optimisation of the WEAVE target assignment algorithm

Authors: Sarah Hughes, Gavin Dalton, Daniel Smith, Kenneth Duncan, David Terrett, Don Carlos Abrams, J. Alfonso Aguerri, Marc Balcells, Georgia Bishop, Piercarlo Bonifacio, Esperansa Carrasco, Shoko Jin, Ian Lewis, Scott Trager, Antonella Vallenari

Abstract: WEAVE is the new wide-field spectroscopic facility for the prime focus of the William Herschel Telescope in La Palma, Spain. Its fibre positioner is essential for the accurate placement of the spectrograph's ~960-fibre multiplex. To maximise the assignment of its optical fibres, WEAVE uses a simulated annealing algorithm called Configure, which allocates the fibres to targets in the field of view.… ▽ More WEAVE is the new wide-field spectroscopic facility for the prime focus of the William Herschel Telescope in La Palma, Spain. Its fibre positioner is essential for the accurate placement of the spectrograph's ~960-fibre multiplex. To maximise the assignment of its optical fibres, WEAVE uses a simulated annealing algorithm called Configure, which allocates the fibres to targets in the field of view. We have conducted an analysis of the algorithm's behaviour using a subset of mid-tier WEAVE-LOFAR fields, and adjusted the priority assignment algorithm to optimise the total fibres assigned per field, and the assignment of fibres to the higher priority science targets. The output distributions have been examined, to investigate the implications for the WEAVE science teams. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 19 pages, 17 figures

Journal ref: Ground-based and Airborne Instrumentation for Astronomy IX 12184 (2022) 121846J

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

Shapes of galaxies hosting radio-loud AGNs with z \leq 1

Authors: X. C. Zheng, H. J. A. Röttgering, A. van der Wel, K. Duncan

Abstract: Links between the properties of radio-loud active galactic nuclei (RLAGNs) and the morphology of their hosts may provide important clues for our understanding of how RLAGNs are triggered. In this work, focusing on passive galaxies, we study the shape of the hosts of RLAGNs selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3GHz Large Project, and compare them wi… ▽ More Links between the properties of radio-loud active galactic nuclei (RLAGNs) and the morphology of their hosts may provide important clues for our understanding of how RLAGNs are triggered. In this work, focusing on passive galaxies, we study the shape of the hosts of RLAGNs selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3GHz Large Project, and compare them with previous results based on the first data release (DR1) of the LOFAR Two-Metre Sky Survey (LoTSS). We find that, at redshifts of between 0.6 and 1, high-luminosity ($L_{1.4 GHz}\gtrsim10^{24}\rm W Hz^{-1}$) RLAGNs have a wider range of optical projected axis ratios than their low-redshift counterparts, which are essentially all found in round galaxies with axis ratios of higher than 0.7. We construct control samples and show that although the hosts of high-redshift RLAGNs with the highest luminosities still have a rounder shape compared with the non-RLAGNs, they on average have a smaller axis ratio (more elongated) than the local RLAGNs with similar stellar masses and radio luminosities. This evolution can be interpreted as a byproduct of radio luminosity evolution, namely that galaxies at fixed stellar mass are more radio luminous at high redshifts: artificially increasing the radio luminosities of local galaxies ($z\leq$0.3) by a factor of 2 to 4 can remove the observed evolution of the axis ratio distribution. If this interpretation is correct then the implication is that the link between AGN radio luminosity and host galaxy shape is similar at $z\simeq1$ to in the present-day Universe. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 figures. Accepted by A&A

Journal ref: A&A 665, A114 (2022)

The detection of cluster magnetic fields via radio source depolarisation

Authors: E. Osinga, R. J. van Weeren, F. Andrade-Santos, L. Rudnick, A. Bonafede, T. Clarke, K. Duncan, S. Giacintucci, Tony Mroczkowski, H. J. A. Röttgering

Abstract: It has been well established that galaxy clusters have magnetic fields. The exact properties and origin of these magnetic fields are still uncertain even though these fields play a key role in many astrophysical processes. Various attempts have been made to derive the magnetic field strength and structure of nearby galaxy clusters using Faraday rotation of extended cluster radio sources. This appr… ▽ More It has been well established that galaxy clusters have magnetic fields. The exact properties and origin of these magnetic fields are still uncertain even though these fields play a key role in many astrophysical processes. Various attempts have been made to derive the magnetic field strength and structure of nearby galaxy clusters using Faraday rotation of extended cluster radio sources. This approach needs to make various assumptions that could be circumvented when using background radio sources. However, because the number of polarised radio sources behind clusters is low, at the moment such a study can only be done statistically. In this paper, we investigate the depolarisation of radio sources inside and behind clusters in a sample of 124 massive clusters at $z<0.35$ observed with the Karl G. Jansky Very Large Array. We detect a clear depolarisation trend with the cluster impact parameter, with sources at smaller projected distances to the cluster centre showing more depolarisation. By combining the radio observations with ancillary X-ray data from Chandra, we compare the observed depolarisation with expectations from cluster magnetic field models using individual cluster density profiles. The best-fitting models have a central magnetic field strength of $5-10\,μ$G with power-law indices between $n=1$ and $n=4$. We find no strong difference in the depolarisation trend between sources embedded in clusters and background sources located at similar projected radii, although the central region of clusters is still poorly probed by background sources. We also examine the depolarisation trend as a function of cluster properties such as the dynamical state, mass, and redshift. Our findings show that the statistical depolarisation of radio sources is a good probe of cluster magnetic field parameters. [abridged] △ Less

Submitted 20 September, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Replaced with Published version (A&A). 32 pages, 34 figures

Journal ref: A&A 665, A71 (2022)

A direct measurement of galaxy major and minor merger rates and stellar mass accretion histories at $z < 3$ using galaxy pairs in the REFINE survey

Authors: Christopher J. Conselice, Carl J. Mundy, Leonardo Ferreira, Kenneth Duncan

Abstract: We measure the role of major and minor mergers in forming the stellar masses of galaxies over $0<z<3$ using a combination of $\sim 3.25$ deg$^{2}$ of the deepest ground based near-infrared imaging taken to date as part of the REFINE survey. We measure the pair fraction and merger fractions for galaxy mergers of different mass ratios, and quantify the merger rate with newly measured time-scales der… ▽ More We measure the role of major and minor mergers in forming the stellar masses of galaxies over $0<z<3$ using a combination of $\sim 3.25$ deg$^{2}$ of the deepest ground based near-infrared imaging taken to date as part of the REFINE survey. We measure the pair fraction and merger fractions for galaxy mergers of different mass ratios, and quantify the merger rate with newly measured time-scales derived from the Illustris simulation as a function of redshift and merger mass ratio. We find that over $0 < z < 3$ major mergers with mass ratios greater than 1:4 occur $0.85^{+0.19}_{-0.20}$ times on average, while minor mergers down to ratios of 1:10 occur on average $1.43^{+0.5}_{-0.3}$ times per galaxy. We also quantify the role of major and minor mergers in galaxy formation, whereby the increase in mass due to major mergers is $93^{+49}_{-31}$% while minor mergers account for an increase of $29^{+17}_{-12}$%; thus major mergers add more stellar mass to galaxies than minor mergers over this epoch. Overall, mergers will more than double the mass of massive galaxies over this epoch. Finally, we compare our results to simulations, finding that minor mergers are over predicted in Illustris and in semi-analytical models, suggesting a mismatch between observations and theory in this fundamental aspect of galaxy assembly. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Comments: AAS Journals submitted, 29 pages

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

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

Authors: Marcin Glowacki, Jordan D. Collier, Amir Kazemi-Moridani, Bradley Frank, Hayley Roberts, Jeremy Darling, Hans-Rainer Klöckner, Nathan Adams, Andrew J. Baker, Matthew Bershady, Tariq Blecher, Sarah-Louise Blyth, Rebecca Bowler, Barbara Catinella, Laurent Chemin, Steven M. Crawford, Catherine Cress, Romeel Davé, Roger Deane, Erwin de Blok, Jacinta Delhaize, Kenneth Duncan, Ed Elson, Sean February, Eric Gawiser , et al. (43 additional authors not shown)

Abstract: In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exc… ▽ More In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 10 pages, 4 figures. Accepted to ApJ Letters

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

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

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

HELP: The Herschel Extragalactic Legacy Project

Authors: R. Shirley, K. Duncan, M. C. Campos Varillas, P. D. Hurley, K. Malek, Y. Roehlly, M. W. L. Smith, H. Aussel, T. Bakx, V. Buat, D. Burgarella, N. Christopher, S. Duivenvoorden, S. Eales, A. Efstathiou, E. A. Gonzalez Solares, M. Griffin, M. Jarvis, B. Lo Faro, L. Marchetti, I. McCheyne, A. Papadopoulos, K. Penner, E. Pons, M. Prescott , et al. (7 additional authors not shown)

Abstract: We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the H… ▽ More We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Galaxy Evolution in all Five CANDELS Fields and IllustrisTNG: Morphological, Structural, and the Major Merger Evolution to $z \sim 3$

Authors: A. Whitney, L. Ferreira, C. J. Conselice, K. Duncan

Abstract: A fundamental feature of galaxies is their structure, yet we are just now understanding the evolution of structural properties in quantitative ways. As such, we explore the quantitative non-parametric structural evolution of 16,778 galaxies up to $z\sim3$ in all five CANDELS fields, the largest collection of high resolution images of distant galaxies to date. Our goal is to investigate how the str… ▽ More A fundamental feature of galaxies is their structure, yet we are just now understanding the evolution of structural properties in quantitative ways. As such, we explore the quantitative non-parametric structural evolution of 16,778 galaxies up to $z\sim3$ in all five CANDELS fields, the largest collection of high resolution images of distant galaxies to date. Our goal is to investigate how the structure, as opposed to size, surface brightness, or mass, changes with time. In particular, we investigate how the concentration and asymmetry of light evolve in the rest-frame optical. To interpret our galaxy structure measurements, we also run and analyse 300 simulation realisations from IllustrisTNG to determine the timescale of mergers for the CAS system. We measure that from $z=0-3$, the median asymmetry merger timescale is $0.56^{+0.23}_{-0.18}$Gyr, and find it does not vary with redshift. Using this data, we find that galaxies become progressively asymmetric at a given mass at higher redshifts and we derive merger rates which scale as $\sim(1+z)^{1.87\pm0.04}$Gyr$^{-1}$, which agrees well with recent machine learning and galaxy pair approaches, removing previous inconsistencies. We also show that far-infrared selected galaxies that are invisible to \textit{HST} have a negligible effect on our measurements. We also find that galaxies are more concentrated at higher redshifts. We interpret this as a sign of how their formation occurs from a smaller initial galaxy that later grows into a larger one through mergers, consistent with the size growth of galaxies from `inside-out', suggesting that the centres are the oldest parts of most galaxies. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: Submitted to ApJ. 26 pages, 19 figures

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 bright end of the infrared luminosity functions and the abundance of hyperluminous infrared galaxies

Authors: L. Wang, F. Gao, P. N. Best, K. Duncan, M. J. Hardcastle, R. Kondapally, K. Malek, I. McCheyne, J. Sabater, T. Shimwell, C. Tasse, M. Bonato, M. Bondi, R. K. Cochrane, D. Farrah, G. Gurkan, P. Haskell, W. J. Pearson, I. Prandoni, H. J. A. Rottgering, D. J. B. Smith, M. Vaccari, W. L. Williams

Abstract: We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match H… ▽ More We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match Herschel sources and LOFAR sources in Bootes (8.63 deg^2), Lockman Hole (10.28 deg^2), and ELAIS-N1 (6.74 deg^2) with rms sensitivities of around 32, 22, and 20 mJy per beam, respectively. We divide the matched samples into unique and multiple categories. For the multiple matches, we de-blend the Herschel fluxes using the LOFAR positions and the 150-MHz flux densities as priors. We perform spectral energy distribution (SED) fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities. The depth of the LOFAR data allows us to identify highly complete (around 92% completeness) samples of bright Herschel sources with a simple selection based on the 250 micron flux density (45, 40, and 35 mJy in Bootes, Lockman Hole, and ELAIS-N1, respectively). Most of the bright Herschel sources fall into the unique category (i.e. a single LOFAR counterpart). For the multiple matches, there is excellent correspondence between the radio emission and the far-IR emission. We find a good agreement in the IR LFs with a previous study out to z around 6 which used de-blended Herschel data. Our sample gives the strongest and cleanest indication to date that the population of HLIRGs has surface densities of around 5 to 18 / deg^2 (with variations due to a combination of the applied flux limit and cosmic variance) and an uncertainty of a factor of 2. In comparison, the GALFORM semi-analytic model significantly under-predicts the abundance of HLIRGs. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 16 pages, 11 figures, accepted for publication in A&A as part of the LOFAR Deep Fields Paper Splash

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)