Merging is potentially the dominate process in galaxy formation, yet there is still debate about its history over cosmic time. To address this we classify major mergers and measure galaxy merger rates up to z $\sim$ 3 in all five CANDELS fields (UDS, EGS, GOODS-S, GOODS-N, COSMOS) using deep learning convolutional neural networks (CNNs) trained with simulated galaxies from the IllustrisTNG cosmological simulation. The deep learning architecture used is objectively selected by a Bayesian Optmization process over the range of possible hyperparameters. We show that our model can achieve 90% accuracy when classifying mergers from the simulation, and has the additional feature of separating mergers before the infall of stellar masses from post mergers. We compare our machine learning classifications on CANDELS galaxies and compare with visual merger classifications from Kartaltepe et al. (2015), and show that they are broadly consistent. We finish by demonstrating that our model is capable of measuring galaxy merger rates, $\mathcal{R}$, that are consistent with results found for CANDELS galaxies using close pairs statistics, with $\mathcal{R}(z) = 0.02 \pm 0.004 \times (1 +z) ^ {2.76 \pm 0.21}$. This is the first general agreement between major mergers measured using pairs and structure at z < 3.

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

We use photometric redshifts derived from new $u$-band through 4.5$\mu$m Spitzer IRAC photometry in the 4.8\,deg$^2$ of the XMM-LSS field to construct surface density maps in the redshift range 0.1-1.5. Our density maps show evidence for large-scale structure in the form of filaments spanning several tens of Mpc. Using these maps, we identify 339 overdensities that our simulated lightcone analysis suggests are likely associated with dark matter haloes with masses, $M_{\rm halo}$, log($M_{\rm halo}/M_{\odot})>$13.7. From this list of overdensities we recover 43 of 70 known X-ray detected and spectroscopically confirmed clusters. The missing X-ray clusters are largely at lower redshifts and lower masses than our target log($M_{\rm halo}/M_{\odot})>$13.7. The bulk of the overdensities are compact, but a quarter show extended morphologies which include likely projection effects, clusters embedded in apparent filaments as well as at least one potential cluster merger (at $z\sim1.28$). The strongest overdensity in our highest redshift slice (at $z\sim1.5$) shows a compact red galaxy core potentially implying a massive evolved cluster.

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

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 a size analysis of a sample of $\sim$ 49,000 galaxies from the CANDELS GOODS North and South fields using redshift-independent relative surface brightness metrics to determine an unbiased measure of the differential size evolution of galaxies at $1 \leq z \leq 7$. We introduce a novel method of removing foreground objects from distant galaxy ($z > 3$) images that makes use of the Lyman-break at 912{\AA}, in what we call `2-D Lyman-Break Imaging'. The images used are in the rest-frame optical at $z < 3$ and progressively bluer bands at $z > 3$. They are therefore subject to K-correction and cosmological dimming effects which are tested and corrected for. We separately consider a mass-selected sample (with masses in the range 10$^9$M$_{\odot}$$\leq$M$_*$$\leq$10$^{10.5}$M$_{\odot}$) and a number density selected sample (using a constant number density of $n = 1\times10^{-4}$Mpc$^{-3}$). Instead of utilising the commonly used, but potentially biased, effective radii for size measurements, we measure the redshift-independent Petrosian radius, defined by the parameter $\eta$, for each galaxy for three values of $\eta$ and use this as a proxy for size. The evolution of the measured radii can be described by a power-law of the form $R_{Petr} = \alpha(1+z)^\beta$kpc where $\beta < 0$. We find that the outer radius increases more rapidly, suggesting that as a galaxy grows mass is added to its outer regions via an inside-out growth. This growth is stronger for the number density selected sample, with a growth rate of nearly three in the outer radii compared to the inner. We test and confirm these results using a series of image simulations.

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

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.

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

The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-cell statistics and differential source counts statistic and measure the angular two-point correlation function of the LoTSS radio sources. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies. Simulated catalogues from the SKA design study sky match with our findings, in contrast to the more recent T-RECS models. For angular separations between 0.1 deg and 1 deg we find self-consistency amongst different estimates of the angular two-point correlation. The angular two-point correlation is $<10^{-2}$ at angular scales $> 1$ deg and up to the largest scales probed. We see indications of unidentified systematic issues at angular scales $>1$ deg, likely due to data calibration issues. Based on the distribution of photometric redshifts of LoTSS sources and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 1 deg agrees with the measured clustering for a subsample of radio sources with redshift information. We find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level and below angular separations of 1 deg the angular two-point correlation agrees with the expectation of the cosmological standard model. (abbreviated)

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

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

Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $\Gamma \approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation.

We present new measurements of the evolution of the galaxy stellar mass functions (GSMF) and UV luminosity functions (UV LF) for galaxies from $z=6-9$ within the Frontier Field cluster MACSJ0416.1-2403 and its parallel field. To obtain these results, we derive the stellar masses of our sample by fitting synthetic stellar population models to their observed spectral energy distribution with the inclusion of nebular emission lines. This is the deepest and farthest in distance mass function measured to date and probes down to a level of M$_{*} = 10^{6.8}M_{\odot}$. The main result of this study is that the low-mass end of our GSMF to these limits and redshifts appears to become steeper from $-1.98_{-0.07}^{+0.07}$ at $z=6$ to $-2.38_{-0.88}^{+0.72}$ at $z=9$, steeper than previously observed mass functions at slightly lower redshifts, and we find no evidence of turnover in the mass range probed. We furthermore demonstrate that the UV LF for these system also appears to show a steepening at the highest redshifts, without any evidence of turnover in the luminosity range probed. Our $M_{\mathrm{UV}}-M_{*}$ relation exhibit shallower slopes than previously observed and are in accordance with a constant mass-to-light ratio. Integrating our GSMF, we find that the stellar mass density increases by a factor of $\sim15_{-6}^{+21}$ from $z=9$ to $z=6$. We estimate the dust-corrected star formation rates (SFRs) to calculate the specific star formation rates ($\mathrm{sSFR}=\mathrm{SFR/M_{*}}$) of our sample, and find that for a fixed stellar mass of $5\times10^{9}M_{\odot}$, sSFR $\propto(1+z)^{2.01\pm0.16}$. Finally, from our new measurements, we estimate the UV luminosity density ($\rho_{\textrm{UV}}$) and find that our results support a smooth decline of $\rho_{\textrm{UV}}$ towards high redshifts.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

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

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.

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

Double-double radio galaxies (DDRGs) represent a short but unique phase in the life-cycle of some of the most powerful radio-loud active galactic nuclei (RLAGN). These galaxies display large-scale remnant radio plasma in the intergalactic medium left behind by a past episode of active galactic nuclei (AGN) activity, and meanwhile, the radio jets have restarted in a new episode. The knowledge of what causes the jets to switch off and restart is crucial to our understanding of galaxy evolution, while it is important to know if DDRGs form a host galaxy dichotomy relative to RLAGN. We utilised the LOFAR Two-Metre Sky Survey DR1, using a visual identification method to compile a sample of morphologically selected candidate DDRGs, showing two pairs of radio lobes. To confirm the restarted nature in each of the candidate sources, we obtained follow-up observations with the VLA at higher resolution to observe the inner lobes or restarted jets, the confirmation of which created a robust sample of 33 DDRGs. We created a comparison sample of 777 RLAGN from the DR1 catalogue, and compared the optical and infrared magnitudes and colours of their host galaxies. We find that there is no statistically significant difference in the brightness of the host galaxies between double-doubles and single-cycle RLAGN. The DDRG and RLAGN samples also have similar distributions in WISE mid-infrared colours, indicating similar ages of stellar populations and dust levels in the hosts of DDRGs. We conclude that DDRGs and 'normal' RLAGN are hosted by galaxies of the same type, and that DDRG activity is simply a normal part of the life cycle of RLAGN. Restarted jets, particularly for the class of low-excitation radio galaxies, rather than being a product of a particular event in the life of a host galaxy, must instead be caused by smaller scale changes, such as in the accretion system surrounding the black hole.

The radio-loud/radio-quiet (RL/RQ) dichotomy in quasars is still an open question. Although it is thought that accretion onto supermassive black holes in the centre the host galaxies of quasars is responsible for some radio continuum emission, there is still a debate as to whether star formation or active galactic nuclei (AGN) activity dominate the radio continuum luminosity. To date, radio emission in quasars has been investigated almost exclusively using high-frequency observations in which the Doppler boosting might have an important effect on the measured radio luminosity, whereas extended structures, best observed at low radio frequencies, are not affected by the Doppler enhancement. We used a sample of quasars selected by their optical spectra in conjunction with sensitive and high-resolution low-frequency radio data provided by the LOw Frequency ARray (LOFAR) as part of the LOFAR Two-Metre Sky Survey (LoTSS) to investigate their radio properties using the radio loudness parameter ($\mathcal{R} = \frac{L_{\mathrm{144-MHz}}}{L_{\mathrm{i\,band}}}$). The examination of the radio continuum emission and RL/RQ dichotomy in quasars exhibits that quasars show a wide continuum of radio properties (i.e. no clear bimodality in the distribution of $\mathcal{R}$). Radio continuum emission at low frequencies in low-luminosity quasars is consistent with being dominated by star formation. We see a significant albeit weak dependency of $\mathcal{R}$ on the source nuclear parameters. For the first time, we are able to resolve radio morphologies of a considerable number of quasars. All these crucial results highlight the impact of the deep and high-resolution low-frequency radio surveys that foreshadow the compelling science cases for the Square Kilometre Array (SKA).

An understanding of the relationship between radio-loud active galaxies and their large-scale environments is essential for realistic modelling of radio-galaxy evolution and environmental impact, for understanding AGN triggering and life cycles, and for calibrating galaxy feedback in cosmological models. We use the LOFAR Two-Metre Sky Survey (LoTSS) Data Release 1 catalogues to investigate this relationship. We cross-matched 8,745 radio-loud AGN with 0.08<z<0.4, selected from LoTSS, with two Sloan Digital Sky Survey cluster catalogues, and find that only 10 percent of LoTSS AGN in this redshift range have an association, so that the majority of low-redshift AGN (including a substantial fraction of the most radio-luminous objects) must inhabit haloes with M < 10^14 M_sun. We find that the probability of a cluster association, and the richness of the associated cluster, is correlated with AGN radio luminosity, and for the cluster population, the number of associated AGN and the radio luminosity of the brightest associated AGN is richness-dependent. We demonstrate that these relations are not driven solely by host-galaxy stellar mass, supporting models in which large-scale environment is influential in driving AGN jet activity. At the lowest radio luminosities we find that the minority of objects with a cluster association are located at larger mean cluster-centre distances than more luminous AGN, which appears to be driven primarily by host-galaxy mass. Finally, we also find that FRI radio galaxies inhabit systematically richer environments than FRIIs. The work presented here demonstrates the potential of LoTSS for AGN environmental studies. In future, the full northern-sky LoTSS catalogue, together with the use of deeper optical/IR imaging data and spectroscopic follow-up with WEAVE-LOFAR, will provide opportunities to extend this type of work to much larger samples and higher redshifts.

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

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections.

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

Historically, the blazar population has been poorly understood at low frequencies because survey sensitivity and angular resolution limitations have made it difficult to identify megahertz counterparts. We used the LOFAR Two-Metre Sky Survey (LoTSS) first data release value-added catalogue (LDR1) to study blazars in the low-frequency regime with unprecedented sensitivity and resolution. We identified radio counterparts to all $98$ known sources from the Third \textit{Fermi}-LAT Point Source Catalogue (3FGL) or Roma-BZCAT Multi-frequency Catalogue of Blazars ($5^{\mathrm{th}}$ edition) that fall within the LDR1 footprint. Only the 3FGL unidentified $\gamma$-ray sources (UGS) could not be firmly associated with an LDR1 source; this was due to source confusion. We examined the redshift and radio luminosity distributions of our sample, finding flat-spectrum radio quasars (FSRQs) to be more distant and more luminous than BL Lacertae objects (BL Lacs) on average. Blazars are known to have flat spectra in the gigahertz regime but we found this to extend down to $144$ MHz, where the radio spectral index, $\alpha$, of our sample is $-0.17 \pm 0.14$. For BL Lacs, $\alpha = -0.13 \pm 0.16$ and for FSRQs, $\alpha = -0.15 \pm 0.17$. We also investigated the radio-to-$\gamma$-ray connection for the $30$ $\gamma$-ray-detected sources in our sample. We find Pearson's correlation coefficient is $0.45$ ($p = 0.069$). This tentative correlation and the flatness of the spectral index suggest that the beamed core emission contributes to the low-frequency flux density. We compare our sample distribution with that of the full LDR1 on colour-colour diagrams, and we use this information to identify possible radio counterparts to two of the four UGS within the LDR1 field. We will refine our results as LoTSS continues.

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.

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.

The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift.

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

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.

The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We perform SED fitting on the UV/NIR to FIR emission of 42 047 galaxies from the pilot HELP field: ELAIS N1. We use the latest release of CIGALE, a galaxy SED fitting code relying on energy balance, to deliver the main physical parameters such as $M_{star}$, SFR, and $L_{dust}$. We implement additional quality criteria to the fits by calculating $\chi^2$ values for the stellar and dust part of the spectra independently. These criteria allow us to identify the best fits and to identify peculiar galaxies. We perform the SED fitting by assuming three different dust attenuation laws separately allowing us to test the impact of the assumed law on estimated physical parameters. We implemented two additional quality value checks for the SED fitting method based on $M_{star}$ estimation and energy budget. This method allows us to identify possible objects with incorrect matching in the catalogue and peculiar galaxies; we found 351 possible candidates of lensed galaxies using two complementary $\chi^2$s criteria (stellar and IR) and $z_{phot}$ calculated for the IR part of the spectrum only. We find that the attenuation law has an important impact on the $M_{star}$ estimate (on average leading to disparities of a factor of two). We derive the relation between $M_{star}$ estimates obtained by different attenuation laws and we find the best recipe for our sample. We also make independent estimates of the total $L_{dust}$ parameter from stellar emission by fitting the galaxies with and without IR data separately.

The radio and far-infrared luminosities of star-forming galaxies are tightly correlated over several orders of magnitude; this is known as the far-infrared radio correlation (FIRC). Previous studies have shown that a host of factors conspire to maintain a tight and linear FIRC, despite many models predicting deviation. This discrepancy between expectations and observations is concerning since a linear FIRC underpins the use of radio luminosity as a star-formation rate indicator. Using LOFAR 150MHz, FIRST 1.4 GHz, and Herschel infrared luminosities derived from the new LOFAR/H-ATLAS catalogue, we investigate possible variation in the monochromatic (250$\mathrm{\mu m}$) FIRC at low and high radio frequencies. We use statistical techniques to probe the FIRC for an optically-selected sample of 4,082 emission-line classified star-forming galaxies as a function of redshift, effective dust temperature, stellar mass, specific star formation rate, and mid-infrared colour (an empirical proxy for specific star formation rate). Although the average FIRC at high radio frequency is consistent with expectations based on a standard power-law radio spectrum, the average correlation at 150MHz is not. We see evidence for redshift evolution of the FIRC at 150MHz, and find that the FIRC varies with stellar mass, dust temperature and specific star formation rate, whether the latter is probed using MAGPHYS fitting, or using mid-infrared colour as a proxy. We can explain the variation, to within 1$\sigma$, seen in the FIRC over mid-infrared colour by a combination of dust temperature, redshift, and stellar mass using a Bayesian partial correlation technique.

Building on the first paper in this series (Duncan et al. 2018), we present a study investigating the performance of Gaussian process photometric redshift (photo-z) estimates for galaxies and active galactic nuclei detected in deep radio continuum surveys. A Gaussian process redshift code is used to produce photo-z estimates targeting specific subsets of both the AGN population - infrared, X-ray and optically selected AGN - and the general galaxy population. The new estimates for the AGN population are found to perform significantly better at z > 1 than the template-based photo-z estimates presented in our previous study. Our new photo-z estimates are then combined with template estimates through hierarchical Bayesian combination to produce a hybrid consensus estimate that outperforms either of the individual methods across all source types. Photo-z estimates for radio sources that are X-ray sources or optical/IR AGN are signficantly improved in comparison to previous template-only estimates, with outlier fractions and robust scatter reduced by up to a factor of ~4. The ability of our method to combine the strengths of the two input photo-z techniques and the large improvements we observe illustrate its potential for enabling future exploitation of deep radio continuum surveys for both the study of galaxy and black hole co-evolution and for cosmological studies.

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.

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

We present a multi-component structural analysis of the internal structure of $1074$ high redshift massive galaxies at $1<z<3$ from the CANDELS HST Survey. In particular we examine galaxies best-fit by two structural components, and thus likely forming discs and bulges. We examine the stellar mass, star formation rates, and colours of both the inner `bulge' and outer `disc' components for these systems using SED information from the resolved ACS+WFC3 HST imaging. We find that the majority of both inner and outer components lie in the star-forming region of UVJ space ($68$ and $90$ per cent respectively). However, the inner portions, or the likely forming bulges, are dominated by dusty star formation. Furthermore, we show that the outer components of these systems have a higher star formation rate than their inner regions, and the ratio of star formation rate between `disc' and `bulge' increases at lower redshifts. Despite the higher star formation rate of the outer component, the stellar mass ratio of inner to outer component remains constant through this epoch. This suggests that there is mass transfer from the outer to inner components for typical two component forming systems, thus building bulges from disks. Finally, using Chandra data we find that the presence of an AGN is more common in both $1$-component spheroid-like galaxies and $2$-component systems ($13\pm3$ and $11\pm2$ per cent) than in $1$-component disc-like galaxies ($3\pm1$ per cent), demonstrating that the formation of a central inner-component likely triggers the formation of central massive black holes in these galaxies.

We present a study of photometric redshift performance for galaxies and active galactic nuclei detected in deep radio continuum surveys. Using two multi-wavelength datasets, over the NOAO Deep Wide Field Survey Bo\"otes and COSMOS fields, we assess photometric redshift (photo-z) performance for a sample of $\sim 4,500$ radio continuum sources with spectroscopic redshifts relative to those of $\sim 63,000$ non radio-detected sources in the same fields. We investigate the performance of three photometric redshift template sets as a function of redshift, radio luminosity and infrared/X-ray properties. We find that no single template library is able to provide the best performance across all subsets of the radio detected population, with variation in the optimum template set both between subsets and between fields. Through a hierarchical Bayesian combination of the photo-z estimates from all three template sets, we are able to produce a consensus photo-z estimate which equals or improves upon the performance of any individual template set.

We use a large sample of $\sim 350,000$ galaxies constructed by combining the UKIDSS UDS, VIDEO/CFHT-LS, UltraVISTA/COSMOS and GAMA survey regions to probe the major merging histories of massive galaxies ($>10^{10}\ \mathrm{M}_\odot$) at $0.005 < z < 3.5$. We use a method adapted from that presented in Lopez-Sanjuan et al. (2014) using the full photometric redshift probability distributions, to measure pair $\textit{fractions}$ of flux-limited, stellar mass selected galaxy samples using close-pair statistics. The pair fraction is found to weakly evolve as $\propto (1+z)^{0.8}$ with no dependence on stellar mass. We subsequently derive major merger $\textit{rates}$ for galaxies at $> 10^{10}\ \mathrm{M}_\odot$ and at a constant number density of $n > 10^{-4}$ Mpc$^{-3}$, and find rates a factor of 2-3 smaller than previous works, although this depends strongly on the assumed merger timescale and likelihood of a close-pair merging. Galaxies undergo approximately 0.5 major mergers at $z < 3.5$, accruing an additional 1-4 $\times 10^{10}\ \mathrm{M}_\odot$ in the process. Major merger accretion rate densities of $\sim 2 \times 10^{-4}$ $\mathrm{M}_\odot$ yr$^{-1}$ Mpc$^{-3}$ are found for number density selected samples, indicating that direct progenitors of local massive ($>10^{11}\mathrm{M}_\odot$) galaxies have experienced a steady supply of stellar mass via major mergers throughout their evolution. While pair fractions are found to agree with those predicted by the Henriques et al. (2014) semi-analytic model, the Illustris hydrodynamical simulation fails to quantitatively reproduce derived merger rates. Furthermore, we find major mergers become a comparable source of stellar mass growth compared to star-formation at $z < 1$, but is 10-100 times smaller than the SFR density at higher redshifts.

We present a study of the low-frequency radio properties of star forming (SF) galaxies and active galactic nuclei (AGN) up to redshift $z=2.5$. The new spectral window probed by the Low Frequency Array (LOFAR) allows us to reconstruct the radio continuum emission from 150 MHz to 1.4 GHz to an unprecedented depth for a radio-selected sample of $1542$ galaxies in $\sim 7~ \rm{deg}^2$ of the LOFAR Bo\"otes field. Using the extensive multi-wavelength dataset available in Bo\"otes and detailed modelling of the FIR to UV spectral energy distribution (SED), we are able to separate the star-formation (N=758) and the AGN (N=784) dominated populations. We study the shape of the radio SEDs and their evolution across cosmic time and find significant differences in the spectral curvature between the SF galaxy and AGN populations. While the radio spectra of SF galaxies exhibit a weak but statistically significant flattening, AGN SEDs show a clear trend to become steeper towards lower frequencies. No evolution of the spectral curvature as a function of redshift is found for SF galaxies or AGN. We investigate the redshift evolution of the infrared-radio correlation (IRC) for SF galaxies and find that the ratio of total infrared to 1.4 GHz radio luminosities decreases with increasing redshift: $ q_{\rm 1.4GHz} = (2.45 \pm 0.04) \times (1+z)^{-0.15 \pm 0.03} $. Similarly, $q_{\rm 150MHz}$ shows a redshift evolution following $ q_{\rm 150GHz} = (1.72 \pm 0.04) \times (1+z)^{-0.22 \pm 0.05}$. Calibration of the 150 MHz radio luminosity as a star formation rate tracer suggests that a single power-law extrapolation from $q_{\rm 1.4GHz}$ is not an accurate approximation at all redshifts.

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.

The evolution of the number density of galaxies in the universe, and thus also the total number of galaxies, is a fundamental question with implications for a host of astrophysical problems including galaxy evolution and cosmology. However there has never been a detailed study of this important measurement, nor a clear path to answer it. To address this we use observed galaxy stellar mass functions up to $z\sim8$ to determine how the number densities of galaxies changes as a function of time and mass limit. We show that the increase in the total number density of galaxies ($\phi_{\rm T}$), more massive than M$_{*} = 10^{6}$ M_0, decreases as $\phi_{\rm T} \sim t^{-1}$, where $t$ is the age of the universe. We further show that this evolution turns-over and rather increases with time at higher mass lower limits of M$_{*}>10^{7}$ M_0. By using the M$_{*}=10^{6}$ M_0 lower limit we further show that the total number of galaxies in the universe up to $z = 8$ is $2.0^{+0.7}_{-0.6} \times 10^{12}$ (two trillion), almost a factor of ten higher than would be seen in an all sky survey at Hubble Ultra-Deep Field depth. We discuss the implications for these results for galaxy evolution, as well as compare our results with the latest models of galaxy formation. These results also reveal that the cosmic background light in the optical and near-infrared likely arise from these unobserved faint galaxies. We also show how these results solve the question of why the sky at night is dark, otherwise known as Olbers' paradox.

We examine a sample of 1495 galaxies in the CANDELS fields to determine the evolution of two component galaxies, including bulges and discs, within massive galaxies at the epoch 1 < z < 3 when the Hubble sequence forms. We fit all of our galaxies' light profiles with a single S\'ersic fit, as well as with a combination of exponential and S\'ersic profiles. The latter is done in order to describe a galaxy with an inner and an outer component, or bulge and disc component. We develop and use three classification methods (visual, F-test and the RFF) to separate our sample into 1-component galaxies (disc/spheroids-like galaxies) and 2-component galaxies (galaxies formed by an 'inner part' or bulge and an 'outer part' or disc). We then compare the results from using these three different ways to classify our galaxies. We find that the fraction of galaxies selected as 2-component galaxies increases on average 50 per cent from the lowest mass bin to the most massive galaxies, and decreases with redshift by a factor of four from z = 1 to z = 3. We find that single S\'ersic 'disc-like' galaxies have the highest relative number densities at all redshifts, and that 2-component galaxies have the greatest increase and become at par with S\'ersic discs by z = 1. We also find that the systems we classify as 2-component galaxies have an increase in the sizes of their outer components, or 'discs', by about a factor of three from z = 3 to z = 1.5, while the inner components or 'bulges' stay roughly the same size. This suggests that these systems are growing from the inside out, whilst the bulges or protobulges are in place early in the history of these galaxies. This is also seen to a lesser degree in the growth of single 'disc-like' galaxies vs. 'spheroid-like' galaxies over the same epoch.

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 \AA\ ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $\alpha=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $\alpha=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant Universe. Here, we model the spectral energy distributions (SEDs) of galaxies at z = 1.5--3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, $L_{\text{TIR}}/L_{\text{UV}}$) and UV slope ($\beta$). We generalize the shape of the dust law with an empirical model, $A_{\lambda,\delta}=E(B-V)\ k_\lambda\ (\lambda/\lambda_V)^\delta$ where $k_\lambda$ is the dust law of Calzetti et al. (2000), and show that there exists a correlation between the color excess ${E(B-V)}$ and tilt $\delta$ with ${\delta=(0.62\pm0.05)\log(E(B-V))}$+ ${(0.26~\pm~0.02)}$. Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star-formation rate, or $\beta$. The change in the dust law with color excess is consistent with a model where attenuation is caused by by scattering, a mixed star-dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher ($z>3$) redshifts.

Due to significant galaxy contamination and impurity in stellar mass selected samples (up to 95% from z=0-3), we examine the star formation history, quenching time-scales, and structural evolution of galaxies using a constant number density selection with data from the UKIDSS Ultra-Deep Survey field. Using this methodology we investigate the evolution of galaxies at a variety of number densities from $z=0-3$. We find that samples chosen at number densities ranging from $3\times10^{-4}$ to 10$^{-5}$ galaxies Mpc$^{-3}$ (corresponding to $z\sim0.5$ stellar masses of M$_{*}= 10^{10.95-11.6}$ M$_{0}$) have a star forming blue fraction of $\sim50$\% at $z\sim2.5$, which evolves to a nearly $100$\% quenched red and dead population by $z\sim 1$. We also see evidence for number density downsizing, such that the galaxies selected at the lowest densities (highest masses) become a homogeneous red population before those at higher number densities. Examining the evolution of the colours for these systems furthermore shows that the formation redshift of galaxies selected at these number densities is $z_{\rm form}>3$. The structural evolution through size and Sersic index fits reveal that while there remains evolution in terms of galaxies becoming larger and more concentrated in stellar mass at lower redshifts, the magnitude of the change is significantly smaller than for a mass selected sample. We also find that changes in size and structure continues at $z < 1$, and is coupled strongly to passivity evolution. We conclude that galaxy structure is driving the quenching of galaxies, such that galaxies become concentrated before they become passive.

With a large sample of bright, low-redshift galaxies with optical$-$near-IR imaging from the GAMA survey we use bulge-disc decompositions to understand the wavelength-dependent behavior of single-S\'ersic structural measurements. We denote the variation in single-S\'ersic index with wavelength as $\mathcal{N}$, likewise for effective radius we use $\mathcal{R}$. We find that most galaxies with a substantial disc, even those with no discernable bulge, display a high value of $\mathcal{N}$. The increase in S\'ersic index to longer wavelengths is therefore intrinsic to discs, apparently resulting from radial variations in stellar population and/or dust reddening. Similarly, low values of $\mathcal{R}$ ($<$ 1) are found to be ubiquitous, implying an element of universality in galaxy colour gradients. We also study how bulge and disc colour distributions vary with galaxy type. We find that, rather than all bulges being red and all discs being blue in absolute terms, both components become redder for galaxies with redder total colours. We even observe that bulges in bluer galaxies are typically bluer than discs in red galaxies, and that bulges and discs are closer in colour for fainter galaxies. Trends in total colour are therefore not solely due to the colour or flux dominance of the bulge or disc.

We present a new analysis of the ionizing emissivity ($\dot{N}_{\rm{ion}}$, s$^{-1}$ Mpc$^{-3}$) for galaxies during the epoch of reionization and their potential for completing and maintaining reionization. We use extensive SED modelling -- incorporating two plausible mechanisms for the escape of Lyman continuum photon -- to explore the range and evolution of ionizing efficiencies consistent with new results on galaxy colours ($\beta$) during this epoch. We estimate $\dot{N}_{\rm{ion}}$ for the latest observations of the luminosity and star-formation rate density at $z<10$, outlining the range of emissivity histories consistent with our new model. Given the growing observational evidence for a UV colour-magnitude relation in high-redshift galaxies, we find that for any plausible evolution in galaxy properties, red (brighter) galaxies are less efficient at producing ionizing photons than their blue (fainter) counterparts. The assumption of a redshift and luminosity evolution in $\beta$ leads to two important conclusions. Firstly, the ionizing efficiency of galaxies naturally increases with redshift. Secondly, for a luminosity dependent ionizing efficiency, we find that galaxies down to a rest-frame magnitude of $M_{\rm{UV}} \approx -15$ alone can potentially produce sufficient numbers of ionizing photons to maintain reionization as early as $z\sim8$ for a clumping factor of $C_{\rm{H {\small II}}} \leq 3$.

We combine photometry from the UDS, and CANDELS UDS and CANDELS GOODS-S surveys to construct the galaxy stellar mass function probing both the low and high mass end accurately in the redshift range 0.3<z<3. The advantages of using a homogeneous concatenation of these datasets include meaningful measures of environment in the UDS, due to its large area (0.88 deg^2), and the high resolution deep imaging in CANDELS (H_160 > 26.0), affording us robust measures of structural parameters. We construct stellar mass functions for the entire sample as parameterised by the Schechter function, and find that there is a decline in the values of phi and of alpha with higher redshifts, and a nearly constant M* up to z~3. We divide the galaxy stellar mass function by colour, structure, and environment and explore the links between environmental over-density, morphology, and the quenching of star formation. We find that a double Schechter function describes galaxies with high Sersic index (n>2.5), similar to galaxies which are red or passive. The low-mass end of the n>2.5 stellar mass function is dominated by blue galaxies, whereas the high-mass end is dominated by red galaxies. This hints that possible links between morphological evolution and star formation quenching are only present in high-mass galaxies. This is turn suggests that there are strong mass dependent quenching mechanisms. In addition, we find that the number density of high mass systems is elevated in dense environments, suggesting that an environmental process is building up massive galaxies quicker in over densities than in lower densities.

We present a study on the stellar mass growth of the progenitors of local massive galaxies with a variety of number density selections with $n\le1\times10^{-4}\,\rm{Mpc^{-3}}$ (corresponding to $M_*=10^{11.24}\rm{M_{\odot}}$ at z=0.3) in the redshift range $0.3<z<3.0$. We select the progenitors of massive galaxies using a constant number density selection, and one which is adjusted to account for major mergers. We find that the progenitors of massive galaxies grow by a factor of four in total stellar mass over this redshift range. On average the stellar mass added via the processes of star formation, major and minor mergers account for $24\pm8\%$, $17\pm15\%$ and $34\pm14\%$, respectively, of the total galaxy stellar mass at $z=0.3$. Therefore $51\pm20\%$ of the total stellar mass in massive galaxies at $z=0.3$ is created externally to their z=3 progenitors. We explore the implication of these results on the cold gas accretion rate and size evolution of the progenitors of most massive galaxies over the same redshift range. We find an average gas accretion rate of $\sim66\pm32\,\rm{M_{\odot}yr^{-1}}$ over the redshift range of $1.5<z<3.0$. We find that the size evolution of a galaxy sample selected this way is on average lower than the findings of other investigations.

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed.

The ability to resolve all processes which drive galaxy formation is one of the most fundamental goals in extragalactic astronomy. While star formation rates and the merger history are now measured with increasingly high certainty, the role of gas accretion from the intergalactic medium in supplying gas for star formation still remains largely unknown. We present in this paper indirect evidence for the accretion of gas into massive galaxies with initial stellar masses M_*>10^{11} M_sol and following the same merger adjusted co-moving number density at lower redshifts during the epoch 1.5 < z < 3, using results from the GOODS NICMOS Survey (GNS). We show that the measured gas mass fractions of these massive galaxies are inconsistent with the observed star formation history for the same galaxy population. We further demonstrate that this additional gas mass cannot be accounted for by cold gas delivered through minor and major mergers. We also consider the effects of gas outflows and gas recycling due to stellar evolution in these calculations. We argue that to sustain star formation at the observed rates there must be additional methods for increasing the cold gas mass, and that the likeliest method for establishing this supply of gas is by accretion from the intergalactic medium. We calculate that the average gas mass accretion rate into these massive galaxies between 1.5 < z < 3.0, is \dot{M} = 96+/-19 M_sol/yr after accounting for outflowing gas. We show that during this epoch, and for these very massive galaxies, 49+/-20% of baryonic mass assembly is a result of gas accretion and unresolved mergers. However, 66+/-20% of all star formation in this epoch is the result of gas accretion. This reveals that for the most massive galaxies at 1.5< z< 3 gas accretion is the dominant method for instigating new stellar mass assembly.