We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored.

The nature of dark matter sets the timeline for the formation of first collapsed halos and thus affects the sources of reionization. Here, we consider two different models of dark matter: cold dark matter (CDM) and thermal warm dark matter (WDM), and study how they impact the epoch of reionization (EoR) and its 21-cm observables. Using a suite of simulations, we find that in the WDM scenarios, the structure formation on small scales gets suppressed resulting in a smaller number of low mass dark matter halos compared to the CDM scenario. Assuming that the efficiency of sources in producing ionizing photons remain the same, this leads to a lower number of total ionizing photons produced at any given cosmic time and thus in a delay in the reionization process. We also find visual differences in the neutral hydrogen (HI) topology and in 21-cm maps in case of the WDM compared to the CDM. However, differences in the 21-cm power spectra, at the same neutral fraction, are found to be small. Thus, we focus on the non-Gaussianity in the EoR 21-cm signal, quantified through its bispectrum. We find that the 21-cm bispectra (driven by the HI topology) are significantly different in WDM models when compared with CDM, even for same mass averaged neutral fractions. This establishes that the 21-cm bispectrum is a unique and promising way to differentiate between different dark matter models, and can be used to constrain the nature of the dark matter in the future EoR observations.

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $\lambda 1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

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

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

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

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

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

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

High-redshift radio galaxies (HzRGs) are some of the rarest objects in the Universe. They are often found to be the most massive galaxies observed at any epoch and are known to harbour active supermassive black holes that give rise to powerful relativistic jets. Finding such galaxies at high redshifts can shed light on the processes that shaped the most massive galaxies very early in the Universe. We have started a new campaign to identify and follow-up promising radio sources selected at 150 MHz in a bid to identify the most distant radio galaxies and study their properties, both intrinsic and environmental. Here we describe the progress of our campaign so far, highlighting in particular the discovery of the most distant radio galaxy known till date, at z = 5.72.

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

We report on the measurement of the $^{7}$Be($n, p$)$^{7}$Li cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux experimental area (EAR2) of the n_TOF facility at CERN. This reaction plays a key role in the lithium yield of the Big Bang Nucleosynthesis (BBN) for standard cosmology. The only two previous time-of-flight measurements performed on this reaction did not cover the energy window of interest for BBN, and showed a large discrepancy between each other. The measurement was performed with a Si-telescope, and a high-purity sample produced by implantation of a $^{7}$Be ion beam at the ISOLDE facility at CERN. While a significantly higher cross section is found at low-energy, relative to current evaluations, in the region of BBN interest the present results are consistent with the values inferred from the time-reversal $^{7}$Li($p, n$)$^{7}$Be reaction, thus yielding only a relatively minor improvement on the so-called Cosmological Lithium Problem (CLiP). The relevance of these results on the near-threshold neutron production in the p+$^{7}$Li reaction is also discussed.

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

We present a model to predict the luminosity function for radio galaxies and their linear size distribution at any redshift. The model takes a black hole mass function and Eddington ratio distribution as input and tracks the evolution of radio sources, taking into account synchrotron, adiabatic and inverse Compton energy losses. We first test the model at z = 2 where plenty of radio data is available and show that the radio luminosity function (RLF) is consistent with observations. We are able to reproduce the break in luminosity function that separates locally the FRI and FRII radio sources. Our prediction for linear size distribution at z = 2 matches the observed distribution too. We then use our model to predict a RLF and linear size distribution at z = 6, as this is the epoch when radio galaxies can be used as probes of reionisation. We demonstrate that higher inverse Compton losses lead to shorter source lifetimes and smaller sizes at high redshifts. The predicted sizes are consistent with the generally observed trend with redshift. We evolve the z = 2 RLF based on observed quasar space densities at high redshifts, and show that our RLF prediction at z = 6 is consistent. Finally, we predict the detection of 0.63, 0.092 and 0.0025 z>=6 sources per sq. degree at flux density limits of 0.1, 0.5 and 3.5 mJy. We assess the trade-off between coverage area and depth and show that LOFAR surveys with flux density limits of 0.1 and 0.5 mJy would are the most efficient at detecting a large number of z>=6 radio sources.

Low radio frequency surveys are important for testing unified models of radio-loud quasars and radio galaxies. Intrinsically similar sources that are randomly oriented on the sky will have different projected linear sizes. Measuring the projected linear sizes of these sources provides an indication of their orientation. Steep-spectrum isotropic radio emission allows for orientation-free sample selection at low radio frequencies. We use a new radio survey of the Bo\"otes field at 150 MHz made with the Low Frequency Array (LOFAR) to select a sample of radio sources. We identify 44 radio galaxies and 16 quasars with powers $P>10^{25.5}$ W Hz$^{-1}$ at 150 MHz using cross-matched multi-wavelength information from the AGN and Galaxy Evolution Survey (AGES), which provides spectroscopic redshifts. We find that LOFAR-detected radio sources with steep spectra have projected linear sizes that are on average 4.4$\pm$1.4 larger than those with flat spectra. The projected linear sizes of radio galaxies are on average 3.1$\pm$1.0 larger than those of quasars (2.0$\pm$0.3 after correcting for redshift evolution). Combining these results with three previous surveys, we find that the projected linear sizes of radio galaxies and quasars depend on redshift but not on power. The projected linear size ratio does not correlate with either parameter. The LOFAR data is consistent within the uncertainties with theoretical predictions of the correlation between the quasar fraction and linear size ratio, based on an orientation-based unification scheme.

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