Search

The role of radio mode AGN feedback on galaxy evolution is still under debate. In this study, we utilize a combination of radio continuum observations and optical integral field spectroscopic (IFS) data to explore the impact of radio AGN on the evolution of their host galaxies at both global and sub-galactic scales. We construct a comprehensive radio-IFS sample comprising 5578 galaxies with redshift z < 0.15 by cross-matching the LOFAR Two-Metre Sky Survey (LoTSS) with the MaNGA survey. We revisit the tight linear radio continuum - star formation relation and quantify its intrinsic scatter, then use the relation to classify 616 radio-excess AGN with excessive radio luminosities over that expected from their star formation rate. Massive radio AGN host galaxies are predominantly quiescent systems, but the quenching level shows no correlation with the jet luminosity. The mass assembly histories derived from the stellar population synthesis model fitting agree with the cosmological simulations incorporating radio-mode AGN feedback models. We observe that radio AGN hosts grow faster than a control sample of galaxies matched in stellar mass, and the quenching age (~ 5 Gyr) is at larger lookback times than the typical radio jet age (< 1 Gyr). By stacking the spectra in different radial bins and comparing results for radio AGN hosts and their controls, we find emission line excess features in the nuclear region of radio AGN hosts. This indicates that radio AGN are ionizing the surrounding interstellar medium in the vicinity of the nucleus. We also find that ongoing star formation in the outer regions of the galaxy is weaker if a radio jet is detected. Our findings support the scenario that the observed present-day radio AGN activity is not responsible for the past quenching of their hosts, but may help the host galaxies maintain quiescence through ionizing and heating the surrounding gas., Comment: 16 pages, 15 figures, submitted to A&A

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

We present the first deep polarimetric study of Galactic synchrotron emission at low radio frequencies. Our study is based on 21 observations of the European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1) field using the Low-Frequency Array (LOFAR) at frequencies from 114.9 to 177.4 MHz. These data are a part of the LOFAR Two-metre Sky Survey Deep Fields Data Release 1. We used very low-resolution ($4.3'$) Stokes QU data cubes of this release. We applied rotation measure (RM) synthesis to decompose the distribution of polarised structures in Faraday depth, and cross-correlation RM synthesis to align different observations in Faraday depth. We stacked images of about 150 hours of the ELAIS-N1 observations to produce the deepest Faraday cube at low radio frequencies to date, tailored to studies of Galactic synchrotron emission and the intervening magneto-ionic interstellar medium. This Faraday cube covers $\sim36~{\rm deg^{2}}$ of the sky and has a noise of $27~{\rm \mu Jy~PSF^{-1}~RMSF^{-1}}$ in polarised intensity. This is an improvement in noise by a factor of approximately the square root of the number of stacked data cubes ($\sim\sqrt{20}$), as expected, compared to the one in a single data cube based on five-to-eight-hour observations. We detect a faint component of diffuse polarised emission in the stacked cube, which was not detected previously. Additionally, we verify the reliability of the ionospheric Faraday rotation corrections estimated from the satellite-based total electron content measurements to be of $~\sim0.05~{\rm rad~m^{-2}}$. We also demonstrate that diffuse polarised emission itself can be used to account for the relative ionospheric Faraday rotation corrections with respect to a reference observation., Comment: 15 pages, 15 figures, accepted for publication in A&A

Feedback from both supermassive black holes and massive stars plays a fundamental role in the evolution of galaxies and the inter-galactic medium. In this paper we use available data to estimate the total amount of kinetic energy and momentum created per co-moving volume element over the history of the universe from three sources: massive stars and supernovae, radiation pressure and winds driven by supermassive black holes, and radio jets driven by supermassive black holes. Kinetic energy and momentum injection from jets peaks at z ~ 1, while the other two sources peak at z ~ 2. Massive stars are the dominant global source of momentum injection. For supermassive black holes, we find that the amount of kinetic energy from jets is about an order-of-magnitude larger than that from winds. We also find that amount of kinetic energy created by massive stars is about 2.5 epsilon times that carried by jets (where epsilon is the fraction of injected energy not lost to radiative cooling). We discuss the implications of these results for the evolution of galaxies and the IGM. Because the ratio of black hole mass to galaxy mass is a steeply increasing function of mass, we show that the relative importance of black hole feedback to stellar feedback likewise increases with mass. We show that there is a trend in the present-day universe which, in the simplest picture, is consistent with galaxies that have been dominated by black hole feedback being generally quenched, while galaxies that have been dominated by stellar feedback are star-forming. We also note that the amount of kinetic energy carried by jets and winds appears sufficient to explain the properties of hot gas in massive halos (> 10^13 solar masses).

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

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

Deep narrow-band surveys have revealed a large population of faint Lyman-alpha (Lya) emitters (LAEs) in the distant Universe, but relatively little is known about the most luminous sources ($L_{Ly\alpha}>10^{42.7}$ erg/s; $L_{Ly\alpha}>L^*_{Ly\alpha}$). Here we present the spectroscopic follow-up of 21 luminous LAEs at z~2-3 found with panoramic narrow-band surveys over five independent extragalactic fields (~4x10$^6$ Mpc$^{3}$ surveyed at z~2.2 and z~3.1). We use WHT/ISIS, Keck/DEIMOS and VLT/X-SHOOTER to study these sources using high ionisation UV lines. Luminous LAEs at z~2-3 have blue UV slopes ($\beta=-2.0^{+0.3}_{-0.1}$), high Lya escape fractions ($50^{+20}_{-15}$%) and span five orders of magnitude in UV luminosity ($M_{UV}\approx-19$ to -24). Many (70%) show at least one high ionisation rest-frame UV line such as CIV, NV, CIII], HeII or OIII], typically blue-shifted by ~100-200 km/s relative to Lya. Their Lya profiles reveal a wide variety of shapes, including significant blue-shifted components and widths from 200 to 4000 km/s. Overall, 60+-11% appear to be AGN dominated, and at $L_{Ly\alpha}>10^{43.3}$ erg/s and/or $M_{UV}<-21.5$ virtually all LAEs are AGN with high ionisation parameters (log U=0.6+-0.5) and with metallicities of ~0.5-1 Zsun. Those lacking signatures of AGN (40+-11%) have lower ionisation parameters ($\log U=-3.0^{+1.6}_{-0.9}$ and $\log\xi_{\rm ion}=25.4\pm0.2$) and are apparently metal-poor sources likely powered by young, dust-poor "maximal" starbursts. Our results show that luminous LAEs at z~2-3 are a diverse population and that 2xL$^*_{Ly\alpha}$ and 2xM$_{UV}^*$ mark a sharp transition in the nature of LAEs, from star formation dominated to AGN dominated., Comment: Accepted for publication in MNRAS. 1D spectra available with the final published version on MNRAS (Supplementary data)

Radio emission is a key indicator of star-formation activity in galaxies, but the radio luminosity-star formation relation has to date been studied almost exclusively at frequencies of 1.4 GHz or above. At lower radio frequencies the effects of thermal radio emission are greatly reduced, and so we would expect the radio emission observed to be completely dominated by synchrotron radiation from supernova-generated cosmic rays. As part of the LOFAR Surveys Key Science project, the Herschel-ATLAS NGP field has been surveyed with LOFAR at an effective frequency of 150 MHz. We select a sample from the MPA-JHU catalogue of SDSS galaxies in this area: the combination of Herschel, optical and mid-infrared data enable us to derive star-formation rates (SFRs) for our sources using spectral energy distribution fitting, allowing a detailed study of the low-frequency radio luminosity--star-formation relation in the nearby Universe. For those objects selected as star-forming galaxies (SFGs) using optical emission line diagnostics, we find a tight relationship between the 150 MHz radio luminosity ($L_{150}$) and SFR. Interestingly, we find that a single power-law relationship between $L_{150}$ and SFR is not a good description of all SFGs: a broken power law model provides a better fit. This may indicate an additional mechanism for the generation of radio-emitting cosmic rays. Also, at given SFR, the radio luminosity depends on the stellar mass of the galaxy. Objects which were not classified as SFGs have higher 150-MHz radio luminosity than would be expected given their SFR, implying an important role for low-level active galactic nucleus activity., Comment: 21 pages, 12 figures, accepted for publication in MNRAS

Feedback by radio-loud active galactic nuclei (AGN) in galaxy groups is not fully understood. Open questions include the duty cycle of the AGN, the spatial extent of the radio lobes, the effect they have on the intragroup medium, and the fate of the cosmic rays. We present the discovery of a 650 kpc-radio galaxy embedded in steep diffuse emission at $z = 0.18793 \pm 5 \times 10^{-5}$ located at the center of the galaxy group MaxBCG J199.31832+51.72503 using an observation from the LOFAR Two-meter Sky Survey (LoTSS) at the central frequency of 144 MHz. Subsequently, we performed a GMRT observation at the central frequency of 607 MHz to study the spectral properties of the source. The observations reveal a radio galaxy with a total radio power $P_{\rm tot, 1.4} \sim 2.1 \times 10^{24}$ W Hz$^{-1}$, exhibiting two asymmetrical jets and lobes. The derived spectral index map shows a steepening toward the inner regions and a steep-spectrum core region. We model the integrated radio spectrum, providing two possible interpretations: the radio source is evolved but still active or it is just at the end of its active phase. Finally, in the same field of view we have discovered Mpc-sized emission surrounding a close pair of AGN located at a redshift $z = 0.0587 \pm 2 \times 10^{-4}$ (SDSS J131544.56+521213.2 and SDSS J131543.99+521055.7) which could be a radio remnant source., Comment: 14 pages, 9 figures

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., Comment: 20 pages, 17 figures. Accepted for publication in MNRAS

We present new results from near-infrared spectroscopy with Keck/MOSFIRE of [OIII]-selected galaxies at $z\sim3.2$. With our $H$ and $K$-band spectra, we investigate the interstellar medium (ISM) conditions, such as ionization states and gas metallicities. [OIII] emitters at $z\sim3.2$ show a typical gas metallicity of $\mathrm{12+log(O/H) = 8.07\pm0.07}$ at $\mathrm{log(M_*/M_\odot) \sim 9.0-9.2}$ and $\mathrm{12+log(O/H) = 8.31\pm0.04}$ at $\mathrm{log(M_*/M_\odot) \sim 9.7-10.2}$ when using the empirical calibration method. We compare the [OIII] emitters at $z\sim3.2$ with UV-selected galaxies and Ly$\alpha$ emitters at the same epoch and find that the [OIII]-based selection does not appear to show any systematic bias in the selection of star-forming galaxies. Moreover, comparing with star-forming galaxies at $z\sim2$ from literature, our samples show similar ionization parameters and gas metallicities as those obtained by the previous studies using the same calibration method. We find no strong redshift evolution in the ISM conditions between $z\sim3.2$ and $z\sim2$. Considering that the star formation rates at a fixed stellar mass also do not significantly change between the two epochs, our results support the idea that the stellar mass is the primary quantity to describe the evolutionary stages of individual galaxies at $z>2$., Comment: 20 pages, 11 figures, Accepted for publication in ApJ

We investigate the clustering properties of $\sim 7000$ H$\beta$+[OIII] and [OII] narrowband-selected emitters at $z \sim 0.8 - 4.7$ from the High-$z$ Emission Line Survey. We find clustering lengths, $r_0$, of $1.5 - 4.0h^{-1}$ Mpc and minimum dark matter halo masses of $10^{10.7 - 12.1}\rm{M}_\odot$ for our $z = 0.8 - 3.2$ H$\beta$+[OIII] emitters and $r_0 \sim 2.0 - 8.3h^{-1}$ Mpc and halo masses of $10^{11.5 - 12.6}\rm{M}_\odot$ for our $z = 1.5 - 4.7$ [OII] emitters. We find $r_0$ to strongly increase both with increasing line luminosity and redshift. By taking into account the evolution of the characteristic line luminosity, $L^\star(z)$, and using our model predictions of halo mass given $r_0$, we find a strong, redshift-independent increasing trend between $L/L^\star(z)$ and minimum halo mass. The faintest H$\beta$+[OIII] emitters are found to reside in $10^{9.5}\rm{M}_\odot$ halos and the brightest emitters in $10^{13.0}\rm{M}_\odot$ halos. For [OII] emitters, the faintest emitters are found in $10^{10.5} \rm{M}_\odot$ halos and the brightest emitters in $10^{12.6}\rm{M}_\odot$ halos. A redshift-independent stellar mass dependency is also observed where the halo mass increases from $10^{11}\rm{M}_\odot$ to $10^{12.5} \rm{M}_\odot$ for stellar masses of $10^{8.5}\rm{M}_\odot$ to $10^{11.5}\rm{M}_\odot$, respectively. We investigate the interdependencies of these trends by repeating our analysis in a $L_\textrm{line} - \rm{M}_\textrm{star}$ grid space for our most populated samples (H$\beta$+[OIII] $z = 0.84$ and [OII] $z = 1.47$) and find that the line luminosity dependency is stronger than the stellar mass dependency on halo mass. For $L > L^\star$ emitters at all epochs, we find a relatively flat trend with halo masses of $10^{12.5 - 13}\rm{M}_\odot$ which may be due to quenching mechanisms in massive halos which is consistent with a transitional halo mass predicted by models., Comment: 19 pages, 11 figures, 4 tables, submitted to MNRAS

We use new near-infrared spectroscopic observations to investigate the nature and evolution of the most luminous H\alpha (Ha) emitters at z~0.8-2.23, which evolve strongly in number density over this period, and compare them to more typical Ha emitters. We study 59 luminous Ha emitters with $L_{H\alpha}>L_{H\alpha}^*$, roughly equally split per redshift slice at z~0.8, 1.47 and 2.23 from the HiZELS and CF-HiZELS surveys. We find that, overall, 30$\pm$8% are AGN (80$\pm$30% of these AGN are broad-line AGN, BL-AGN), and we find little to no evolution in the AGN fraction with redshift, within the errors. However, the AGN fraction increases strongly with Ha luminosity and correlates best with $L_{H\alpha}/L_{H\alpha}^*(z)$. While $L_{H\alpha}80%), the most luminous Ha emitters ($L_{H\alpha}>10L_{H\alpha}^*(z)$) at any cosmic time are essentially all BL-AGN. Using our AGN-decontaminated sample of luminous star-forming galaxies, and integrating down to a fixed Ha luminosity, we find a factor of ~1300x evolution in the star formation rate density from z=0 to z=2.23. This is much stronger than the evolution from typical Ha star-forming galaxies and in line with the evolution seen for constant luminosity cuts used to select "Ultra-Luminous" Infrared Galaxies and/or sub-millimetre galaxies. By taking into account the evolution in the typical Ha luminosity, we show that the most strongly star-forming Ha-selected galaxies at any epoch ($L_{H\alpha}>L^*_{H\alpha}(z)$) contribute the same fractional amount of ~15% to the total star-formation rate density, at least up to z=2.23., Comment: 16 pages, 10 figures. Accepted for publication in MNRAS

We investigate the evolution of the H$\beta$+[OIII] and [OII] luminosity functions from $z \sim 0.8$ to $\sim5$ in four redshift slices per emission line using data from the High-{\it z} Emission Line Survey (HiZELS). This is the first time that the H$\beta$+[OIII] and [OII] luminosity functions have been studied at these redshifts in a self-consistent analysis. This is also the largest sample of [OII] and H$\beta$+[OIII] emitters (3475 and 3298 emitters, respectively) in this redshift range, with large co-moving volumes $\sim 1 \times 10^6$ Mpc$^{-3}$ in two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic variance. The emitters were selected by a combination of photometric redshift and color-color selections, as well as spectroscopic follow-up, including recent spectroscopic observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on Subaru. We find a strong increase in $L_\star$ and a decrease in $\phi_\star$ for both H$\beta$+[OIII] and [OII] emitters. We derive the [OII] star-formation history of the Universe since $z\sim5$ and find that the cosmic SFRD rises from $z \sim 5$ to $\sim 3$ and then drops towards $z \sim 0$. We also find that our star-formation history is able to reproduce the evolution of the stellar mass density up to $z\sim 5$ based only on a single tracer of star-formation. When comparing the H$\beta$+[OIII] SFRDs to the [OII] and H$\alpha$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$\beta$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs., Comment: 24 pages, 14 figures, 7 Tables, MNRAS accepted

We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 $\mu$m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H$\alpha$ emitters at z=0.8, 1343 [OIII]+H$\beta$ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5$\sigma$ (spectroscopically confirmed) over-density of H$\alpha$ emitters at z=0.81. We derive H$\alpha$, [OIII]+H$\beta$ and [OII] luminosity functions at z=0.8,1.4,2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $\phi^*$ due to sample (cosmic) variance on surveys probing ~10$^4$ Mpc$^{3}$ and ~10$^5$ Mpc$^{3}$ are typically very high: ~300% and ~40-60%, respectively., Comment: 22 pages, 15 figures, accepted for publication in MNRAS. Full catalogue of line emitters available in FITS format with final MNRAS published paper

We have observed a sample of typical z=1 star forming galaxies, selected from the HiZELS survey, with the new KMOS near-infrared, multi-IFU instrument on the VLT, in order to obtain their dynamics and metallicity gradients. The majority of our galaxies have a metallicity gradient consistent with being flat or negative (i.e. higher metallicity cores than outskirts). Intriguingly, we find a trend between metallicity gradient and specific star formation rate (sSFR), such that galaxies with a high sSFR tend to have relatively metal-poor centres, a result which is strengthened when combined with datasets from the literature. This result appears to explain the discrepancies reported between different high redshift studies and varying claims for evolution. From a galaxy evolution perspective, the trend we see would mean that a galaxy's sSFR is governed by the amount of metal poor gas that can be funnelled into its core, triggered either by merging or through efficient accretion. In fact merging may play a significant role as it is the starburst galaxies at all epochs, which have the more positive metallicity gradients. Our results may help to explain the origin of the fundamental metallicity relation, in which galaxies at a fixed mass are observed to have lower metallicities at higher star formation rates, especially if the metallicity is measured in an aperture encompassing only the central regions of the galaxy. Finally, we note that this study demonstrates the power of KMOS as an efficient instrument for large scale resolved galaxy surveys., Comment: Accepted for publication in MNRAS. Version 2 corrected typo in equation 2

We explore a large uniformly selected sample of H${\alpha}$ selected star-forming galaxies (SFGs) at z=0.40,0.84,1.47,2.23 to unveil the evolution of the star formation rate (SFR) function and the stellar mass function. We find strong evolution in the SFR function, with the typical SFR of SFGs declining exponentially with time in the last 11Gyrs as SFR$^*$(T[Gyr])=10$^{4.23/T+0.37}$ M$_{\odot}$yr$^{-1}$, but with no evolution in the faint-end slope, ${\alpha}\approx$-1.6. The stellar mass function of SFGs, however, reveals little evolution: ${\alpha} \approx$-1.4, M$^*$$\approx$10$^{11.2\pm0.2}$ M$_{\odot}$ and just a slight increase of $\approx$2.3x in ${\Phi}^*$ from z=2.23 to z=0.4. The stellar mass density within SFGs has been roughly constant since z=2.23 at $\approx$10$^{7.65\pm0.08}$ M$_{\odot}$Mpc$^{-3}$, comprising $\approx$100% of the stellar mass density in all galaxies at z=2.23, and declining to $\approx$20% by z=0.4, driven by the rise of the passive population. We find that SFGs with M$\approx$10$^{10.0\pm0.2}$ M$_{\odot}$ contribute most to the SFR density (${\rho}_{\rm SFR}$) per dlog$_{10}$M, and that there is no significant evolution in the fractional contribution from SFGs of different masses to ${\rho}_{\rm SFR}$ or ${\rho}_{\rm SFR}$(dlog$_{10}$M)$^{-1}$ since z=2.23. Instead, we show that the decline of SFR$^*$ and of ${\rho}_{\rm SFR}$ are primarily driven by an exponential decline in SFRs at all masses. Our results have important implications not only on how SFGs need to be quenched across cosmic time, but also on the driver(s) of the exponential decline in SFR$^*$ from $\approx$66 M$_{\odot}$yr$^{-1}$ to $\approx$5 M$_{\odot}$yr$^{-1}$ since z=2.23., Comment: MNRAS in press, 14 pages, 9 figures

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru and the VLT; a combined effort to select large, robust samples of H-alpha (Ha) emitters at z=0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs) in a uniform manner over ~2 deg^2 in the COSMOS and UDS fields. The deep Ha surveys reach ~3M_sun/yr out to z=2.2 for the first time, while the wide area and the coverage over two independent fields allow to greatly overcome cosmic variance. A total of 1742, 637, 515 and 807 Ha emitters are homogeneously selected at z=0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Ha luminosity function and its evolution. The faint-end slope is found to be -1.60+-0.08 over z=0-2.23, showing no evolution. The characteristic luminosity of SF galaxies, L*, evolves significantly as log[L*(z)]=0.45z+log[L*(z=0)]. This is the first time Ha has been used to trace SF activity with a single homogeneous survey at z=0.4-2.23. Overall, the evolution seen in Ha is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards the bulk of the evolution in the last 11 Gyrs being driven by a strong luminosity increase from z~0 to z~2.2. Our uniform analysis reveals an Ha star formation history of the Universe which can be simply parameterised by log(SFRD)=-2.1/(1+z) for z<2.23. Both the shape and normalisation of the Ha star formation history are consistent with the measurements of the stellar mass density growth, confirming that our Ha analysis traces the bulk of the formation of stars in the Universe up to z~2.2. Star formation activity since z~2.2 is responsible for ~95% of the total stellar mass density observed locally today (with half being assembled between z=1.2-2.2, and the rest since z<1.2)., Comment: 20 pages, 11 Figures, accepted for publication in MNRAS. Sample of z=2.23 emitters increased to 807 after completing the COSMOS survey

(Abridged) The clustering properties of a well-defined sample of 734 H-alpha emitters at z=0.84 obtained as part of the Hi-z Emission Line Survey (HiZELS) are investigated. The spatial correlation function is very well-described by (r/r_0)^-1.8, with r_0=2.7+-0.3Mpc/h. The correlation length r_0 increases strongly with H-alpha luminosity, L_H-alpha, from r_0~2Mpc/h for the most quiescent galaxies (star-formation rates of ~4M_sun/yr), up to r_0>5Mpc/h for the brightest galaxies in H-alpha. The correlation length also increases with increasing rest-frame K-band luminosity (M_K), but the r_0-L_H-alpha correlation maintains its full statistical significance at fixed M_K. At z=0.84, star-forming galaxies classified as irregulars or mergers are much more clustered than discs and non-mergers, but once the samples are matched in L_H-alpha and M_K, the differences vanish, implying that the clustering is independent of morphological type at z~1. The typical H-alpha emitters found at z=0.84 reside in dark-matter haloes of ~10^12M_sun, but those with the highest SFRs reside in more massive haloes of ~10^13M_sun. Comparing the results with those of H-alpha surveys at different redshifts, it is seen that although the break of the H-alpha luminosity function, L*, evolves by a factor of ~30 from z=0.24 to z=2.23, galaxies with the same L_H-alpha/L*(z) are found in dark matter haloes of similar masses, independently of cosmic time. This not only confirms that star-formation is more efficient at higher redshift, but also suggests a fundamental connection between the strong decrease of L* since z=2.23 and the quenching of star-formation in galaxies residing within dark-matter haloes significantly more massive than 10^12M_sun at any given epoch., Comment: Accepted for publication in MNRAS, 13 pages, 10 figures, revised version following the reviewer's report

[Abridged] We compare H-alpha, ultraviolet (UV) and infrared (IR) indicators of star formation rate (SFR) for a sample of z = 0.84 galaxies from the High-z Emission Line Survey (HiZELS). Using multi-wavelength diagnostics, we estimate that 5 - 11 per cent of H-alpha emitters at this redshift are active galactic nuclei. We detect 35 per cent of the H-alpha emitters individually at 24 microns, and stack the star-forming emitters on deep 24-micron images in order to calculate the typical SFRs of our galaxies. These are compared to the observed H-alpha line fluxes in order to estimate the extinction at z = 0.84, and we find a significant increase in dust extinction for galaxies with higher SFRs. We demonstrate that the local relationship between SFR and extinction is suitable for our sample, and attribute the overall increase in dust extinction for z = 0.84 galaxies to an increase in the average SFR, rather than to a change in dust properties. We calculate the UV extinction, and find a similar dependence on SFR, but no evidence for a 2175 angstrom UV bump in the dust attenuation law. We calculate the conversion between the dust attenuation of nebular and stellar radiation, gamma, and show that gamma = 0.50 +/- 0.14. The extinction / SFR relationship is shown to be applicable to galaxies with a range of morphologies and bulge-to-disk ratios, to both merging and non-merging galaxies, and to galaxies within high- and low-density environments, implying that it is a fundamental property of star-forming regions. In order to allow future studies to correct for a SFR-dependent amount of dust extinction, we present an equation to predict the extinction of a galaxy, based solely upon its observed H-alpha luminosity, and use this to recalculate the H-alpha luminosity function and star formation rate density at z = 0.84., Comment: Accepted for publication in MNRAS. 15 pages, 8 figures, 4 tables