Source author record

David Sobral

David Sobral appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.IM

Catalog footprint

What is connected

67works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

J-HERTz: J-PLUS Heritage Exploration of Radio Targets at z $<$ 5

We introduce J-HERTz (J-PLUS Heritage Exploration of Radio Targets at $z < 5$), a new multi-wavelength catalog that combines optical narrow-band photometry from J-PLUS, infrared observations from WISE, and deep low-frequency radio data from LoTSS for nearly half a million sources across 2,100 deg$^2$ of the northern sky. Key innovations of J-HERTz include Bayesian neural network classifications for 390,000 galaxies, 31,000 quasars, and 20,000 stars, along with significantly improved photometric redshifts for 235,000 galaxies compared to previous J-PLUS DR3 and LoTSS DR2 estimates. We identify 831 candidate Galactic radio stars, which, if confirmed, would constitute a significant addition to the number of radio-emitting stars identified to date. Among radio-loud galaxies with spectroscopic observations, $\gtrsim$20% lack Seyfert or LINER signatures, indicating a substantial population of optically quiescent radio galaxies, in agreement with previous works. Spectral energy distribution fitting of their host galaxies using J-PLUS photospectra reveals systematically low specific star formation rates, consistent with quenched stellar populations. J-HERTz thus provides a powerful dataset to exploit radio-optical synergies, enabling studies that span from the origin of stellar radio emission to the AGN life cycle and the role of jet activity in shaping host galaxy evolution.

preprint2022arXiv

(Re)Solving Reionization with Lyα: How Bright Lyα Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$α$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$α$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Lyα}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$α$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$α$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $ξ_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$α$-anchored formalism.

preprint2022arXiv

Black hole virial masses from single-epoch photometry: the miniJPAS test case

Precise measurements of black hole masses are essential to understanding the coevolution of these sources and their host galaxies. We develop a novel approach for computing black hole virial masses using measurements of continuum luminosities and emission line widths from partially overlapping, narrow-band observations of quasars; we refer to this technique as single-epoch photometry. This novel method relies on forward-modelling quasar observations for estimating emission line widths, which enables unbiased measurements even for lines coarsely resolved by narrow-band data. We assess the performance of this technique using quasars from the Sloan Digital Sky Survey (SDSS) observed by the miniJPAS survey, a proof-of-concept project of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) collaboration covering $\simeq1\,\mathrm{deg}^2$ of the northern sky using the 56 J-PAS narrow-band filters. We find remarkable agreement between black hole masses from single-epoch SDSS spectra and single-epoch miniJPAS photometry, with no systematic difference between these and a scatter ranging from 0.4 to 0.07 dex for masses from $\log(M_\mathrm{BH})\simeq8$ to 9.75, respectively. Reverberation mapping studies show that single-epoch masses present approximately 0.4 dex precision, letting us conclude that our novel technique delivers black hole masses with only mildly lower precision than single-epoch spectroscopy. The J-PAS survey will soon start observing thousands of square degrees without any source preselection other than the photometric depth in the detection band, and thus single-epoch photometry has the potential to provide details on the physical properties of quasar populations that do not satisfy the preselection criteria of previous spectroscopic surveys.

preprint2022arXiv

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics.

preprint2022arXiv

The LEGA-C and SAMI Galaxy Surveys: Quiescent Stellar Populations and the Mass-Size Plane across 6 Gyr

We investigate the change in mean stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift ($0.60\leq z\leq0.76$) using the LEGA-C Survey, to low redshift ($0.014\leq z\leq0.10$) using the SAMI Galaxy Survey. We find that, similarly to their low-redshift counterparts, the stellar metallicity of quiescent galaxies at $0.60\leq z\leq 0.76$ closely correlates with $M_*/R_\mathrm{e}$ (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals, by determining the escape velocity required by metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. On the other hand, we find no correlation between stellar age and $M_*/R_\mathrm{e}^2$ (stellar mass surface density $Σ$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$Σ$ relation in the context of the redshift evolution of the star-forming and quiescent populations in the mass--size plane, and find our results can be explained as a consequence of galaxies forming more compactly at higher redshifts, and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The $z\sim 0$ age--$Σ$ relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and so a range of surface densities.

preprint2022arXiv

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $σ'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $σ'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a Sérsic profile-based $R$ estimate, then a Sérsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex.

preprint2022arXiv

The miniJPAS survey: Identification and characterization of the emission line galaxies down to $z < 0.35$ in the AEGIS field

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) is expected to map thousands of square degrees of the northern sky with 56 narrowband filters in the upcoming years. This will make J-PAS a very competitive and unbiased emission line survey compared to spectroscopic or narrowband surveys with fewer filters. The miniJPAS survey covered 1 deg$^2$, and it used the same photometric system as J-PAS, but the observations were carried out with the pathfinder J-PAS camera. In this work, we identify and characterize the sample of emission line galaxies (ELGs) from miniJPAS with a redshift lower than $0.35$. Using a method based on artificial neural networks, we detect the ELG population and measure the equivalent width and flux of the $Hα$, $Hβ$, [OIII], and [NII] emission lines. We explore the ionization mechanism using the diagrams [OIII]/H$β$ versus [NII]/H$α$ (BPT) and EW(H$α$) versus [NII]/H$α$ (WHAN). We identify 1787 ELGs ($83$%) from the parent sample (2154 galaxies) in the AEGIS field. For the galaxies with reliable EW values that can be placed in the WHAN diagram (2000 galaxies in total), we obtained that $72.8 \pm 0.4$%, $17.7 \pm 0.4$% , and $9.4 \pm 0.2$% are star-forming (SF), active galactic nucleus (Seyfert), and quiescent galaxies, respectively. Based on the flux of $Hα$ we find that the star formation main sequence is described as $\log$ SFR $[M_\mathrm{\odot} \mathrm{yr}^{-1}] = 0.90^{+ 0.02}_{-0.02} \log M_{\star} [M_\mathrm{\odot}] -8.85^{+ 0.19}_{-0.20}$ and has an intrinsic scatter of $0.20^{+ 0.01}_{-0.01}$. The cosmic evolution of the SFR density ($ρ_{\text{SFR}}$) is derived at three redshift bins: $0 < z \leq 0.15$, $0.15 < z \leq 0.25$, and $0.25 < z \leq 0.35$, which agrees with previous results that were based on measurements of the $Hα$ emission line.

preprint2021arXiv

J-PLUS: Support Vector Machine Applied to STAR-GALAXY-QSOClassification

Context. In modern astronomy, machine learning has proved to be efficient and effective to mine the big data from the newesttelescopes. Spectral surveys enable us to characterize millions of objects, while long exposure time observations and wide surveysconstrain their strides from millions to billions. Aims.In this study, we construct a supervised machine learning algorithm, to classify the objects in the Javalambre Photometric LocalUniverse Survey first data release (J-PLUS DR1). Methods.The sample set is featured with 12-waveband photometry, and magnitudes are labeled with spectrum-based catalogs, in-cluding Sloan Digital Sky Survey spectroscopic data, Large Sky Area Multi-Object Fiber Spectroscopic Telescope, and VERONCAT- Veron Catalog of Quasars & AGN. The performance of the classifier is presented with applications of blind test validations basedon RAdial Velocity Extension, Kepler Input Catalog, 2 MASS Redshift Survey, and the UV-bright Quasar Survey. A new algorithmis applied to constrain the extrapolation that could decrease accuracies for many machine learning classifiers. Results.The accuracies of the classifier are 96.5% in blind test and 97.0% in training cross validation. The F1-scores for each classare presented to show the precision of the classifier. We also discuss different methods to constrain the po

preprint2021arXiv

The LEGA-C of nature and nurture in stellar populations of galaxies at z~0.6-1.0: D4000 and H-delta reveal different assembly histories for quiescent galaxies in different environments

Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z~0.6-1.0 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D4000 and H-delta equivalent widths (EWs) distributions over the stellar mass vs environment 2D spaces for the massive galaxy population (M>10^10 M$_{\odot}$) at z~0.6-1.0. D4000 and H-delta EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D4000 and H-delta EW when controlling for stellar mass are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D4000 and H-delta EWs which depend strongly on their stellar mass and are completely independent of the environment at 0.6<z<1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point towards massive galaxies forming first in over-densities or the most compact dark matter haloes.

preprint2021arXiv

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra.

preprint2020arXiv

Dust, gas, and metal content in star-forming galaxies at $z\sim3.3$ revealed with ALMA and Near-IR spectroscopy

We conducted sub-millimeter observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) of star-forming galaxies at $z\sim3.3$, whose gas-phase metallicities have been previously measured. We investigate the dust and gas contents of the galaxies at $z\sim3.3$ and study how galaxies are interacting with their circumgalactic/intergalactic medium at this epoch by probing their gas mass fractions and gas-phase metallicities. Single-band dust continuum emission tracing dust mass and the relation between the gas-phase metallicity and gas-to-dust mass ratio are used to estimate the gas masses. The estimated gas mass fractions and depletion timescales are $f_{\rm gas}=$ 0.20-0.75 and $t_{\rm dep}=$ 0.09-1.55 Gyr, respectively. Although the galaxies appear to tightly distribute around the star-forming main sequence at $z\sim3.3$, both quantities show a wider spread at a fixed stellar mass than expected from the scaling relation, suggesting a large diversity of fundamental gas properties among star-forming galaxies apparently on the main sequence. Comparing gas mass fraction and gas-phase metallicity between the star-forming galaxies at $z\sim3.3$ and at lower redshifts, star-forming galaxies at $z\sim3.3$ appear to be more metal-poor than local galaxies with similar gas mass fractions. Using the gas regulator model to interpret this offset, we find that it can be explained by a higher mass-loading factor, suggesting that the mass-loading factor in outflows increases at earlier cosmic times.

preprint2020arXiv

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $Hδ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $Hδ_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ Å/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ Å$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius.

preprint2020arXiv

Spectroscopic Confirmation of a Coma Cluster Progenitor at z ~ 2.2

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$α$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $μ$m) and $K$ ($\sim$ 1.92-2.40 $μ$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $σ_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $δ_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $δ_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

preprint2020arXiv

Stellar Kinematics and Environment at z~0.8 in the LEGA-C Survey: Massive, Slow-Rotators are Built First in Overdense Environments

In this letter, we investigate the impact of environment on integrated and spatially-resolved stellar kinematics of a sample of massive, quiescent galaxies at intermediate redshift ($0.6<z<1.0$). For this analysis, we combine photometric and spectroscopic parameters from the UltraVISTA and Large Early Galaxy Astrophysics Census (LEGA-C) surveys in the COSMOS field and environmental measurements. We analyze the trends with overdensity (1+$δ$) on the rotational support of quiescent galaxies and find no universal trends at either fixed mass or fixed stellar velocity dispersion. This is consistent with previous studies of the local Universe; rotational support of massive galaxies depends primarily on stellar mass. We highlight two populations of massive galaxies ($\log M_\star/M_\odot\geq11$) that deviate from the average mass relation. First, the most massive galaxies in the most under-dense regions ($(1+δ)\leq1$) exhibit elevated rotational support. Similarly, at the highest masses ($\log M_\star/M_\odot\geq11.25$) the range in rotational support is significant in all but the densest regions. This corresponds to an increasing slow-rotator fraction such that only galaxies in the densest environments ($(1+δ)\geq3.5$) are primarily (90$\pm$10\%) slow-rotators.This effect is not seen at fixed velocity dispersion, suggesting minor merging as the driving mechanism: only in the densest regions have the most massive galaxies experienced significant minor merging, building stellar mass and diminishing rotation without significantly affecting the central stellar velocity dispersion. In the local Universe, most massive galaxies are slow-rotators, regardless of environment, suggesting minor merging occurs at later cosmic times $(z\lesssim0.6)$ in all but the most dense environments.

preprint2020arXiv

The evolution of rest-frame UV properties, Lya EWs and the SFR-Stellar mass relation at z~2-6 for SC4K LAEs

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $σ$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.

preprint2020arXiv

The nature of CR7 revealed with MUSE: a young starburst powering extended Lyman-$α$ emission at z=6.6

CR7 is among the most luminous Lyman-$α$ emitters (LAEs) known at $z = 6.6$ and consists of at least three UV components that are surrounded by Lyman-$α$ (Ly$α$) emission. Previous studies have suggested that it may host an extreme ionising source. Here, we present deep integral field spectroscopy of CR7 with VLT/MUSE. We measure extended emission with a similar halo scale length as typical LAEs at $z\approx5$. CR7's Ly$α$ halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly$α$ emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly$α$ emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly$α$ equivalent width (EW) is $\approx100$ Å. This is a factor four higher than the EW measured in low-redshift analogues with carefully matched Ly$α$ profiles (and thus arguably HI column density), but this EW can plausibly be explained by star formation. Alternative scenarios requiring AGN powering are also disfavoured by the narrower and steeper Ly$α$ spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly$α$ blobs. CR7's Ly$α$ emission, while extremely luminous, resembles the emission in more common LAEs at lower redshifts very well and is likely powered by a young metal poor starburst.

preprint2020arXiv

The X-ray and radio activity of typical and luminous Lya emitters from z~2 to z~6: evidence for a diverse, evolving population

Despite recent progress in understanding Lyman-alpha (Lya) emitters (LAEs), relatively little is known regarding their typical black hole activity across cosmic time. Here, we study the X-ray and radio properties of ~4000 LAEs at 2.2<z<6 from the SC4K survey in the COSMOS field. We detect 254 (6.8+-0.4%) of the LAEs individually in the X-rays (S/N>3) and find an average luminosity of 10^{44.31+-0.01} erg/s and an average black hole accretion rate (BHAR) of 0.72+-0.01 Msun/yr, consistent with moderate to high accreting AGN. We detect 120 sources in deep radio data (radio AGN fraction of 3.2%+-0.3%). The global AGN fraction (8.6%+-0.4%) rises with Lya luminosity and declines with increasing redshift. For X-ray detected LAEs, Lya luminosities correlate with the BHARs, suggesting that Lya luminosity becomes an accretion rate indicator. Most LAEs (93.1+-0.6%) at 2<z<6 have no detectable X-ray emission (BHARs<0.017 Msun/yr). The median star formation rate (SFR) of star-forming LAEs from Lya and radio luminosities is 7.6^{+6.6}_{-2.8} Msun/yr. The black hole to galaxy growth ratio (BHAR/SFR) for LAEs is <0.0022, consistent with typical star forming galaxies and the local BHAR/SFR relation. We conclude that LAEs at 2<z<6 include two different populations: an AGN population, where Lya luminosity traces BHAR, making them bright in Lya, and another with low SFRs which remain undetected in even the deepest X-ray stacks but is detected in the radio stacks.

preprint2019arXiv

Chandra Observations of the Spectacular A3411-12 Merger Event

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

preprint2019arXiv

Resolved Lyman-α properties of a luminous Lyman-break galaxy in a large ionised bubble at z = 6.53

The observed properties of the Lyman-$α$ (Ly$α$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$α$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$α$ equivalent width (EW$_0\approx38$ Å). These data are combined with deep resolved [CII]$_{\rm 158 μm}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$α$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$α$ and [CII] are similarly extended beyond the UV, with effective radius r$_{\rm eff} = 2.1\pm0.2$ kpc for a single exponential model or r$_{\rm eff, Lyα, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$α$ profile is broader and redshifted with respect to the [CII] line (by 213 km s$^{-1}$), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different HI column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly$α$ halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local ($\sim 10$ kpc) neutral gas conditions that determine the observed Ly$α$ properties in VR7 resemble the conditions in post-re-ionisation galaxies.

preprint2019arXiv

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(Hδ) \geq 4$Å) at $z\sim0.8$ drawn from DR2 of the LEGA-C survey to test the hypothesis that these galaxies experienced compact, central starbursts before quenching. These recently quenched galaxies, usually referred to as post-starburst galaxies, span a wide range of colors and we find a clear correlation between color and half-light radius, such that bluer galaxies are smaller. We build simple toy models to explain this correlation: a normal star-forming disk plus a central, compact starburst component. Bursts with exponential decay timescale of $\sim$~100 Myr that produce $\sim10\%$ to more than 100\% of the pre-existing masses can reproduce the observed correlation. More significant bursts also produce bluer and smaller descendants. Our findings imply that when galaxies shut down star formation rapidly, they generally had experienced compact, starburst events and that the large, observed spread in sizes and colors mostly reflects a variety of burst strengths. Recently quenched galaxies should have younger stellar ages in the centers; multi-wavelength data with high spatial resolution are required to reveal the age gradient. Highly dissipative processes should be responsible for this type of formation history. While determining the mechanisms for individual galaxies is challenging, some recently quenched galaxies show signs of gravitational interactions, suggesting that mergers are likely an important mechanism in triggering the rapid shut-down of star-formation activities at $z\sim0.8$.

preprint2019arXiv

Unveiling the most luminous Lyman-α emitters in the epoch of reionisation

Distant luminous Lyman-alpha emitters are excellent targets for detailed observations of galaxies in the epoch of reionisation. Spatially resolved observations of these galaxies allow us to simultaneously probe the emission from young stars, partially ionised gas in the interstellar medium and to constrain the properties of the surrounding hydrogen in the circumgalactic medium. We review recent results from (spectroscopic) follow-up studies of the rest-frame UV, Lyman-alpha and [CII] emission in luminous galaxies observed ~500 Myr after the Big Bang with ALMA, HST/WFC3 and VLT/X-SHOOTER. These galaxies likely reside in early ionised bubbles and are complex systems, consisting of multiple well separated and resolved components where traces of metals are already present.

preprint2019arXiv

VIS3COS: III. environmental effects on the star formation histories of galaxies at z~0.8 seen in [OII], H$δ$, and Dn4000

[ABRIDGED] We present spectroscopic observations of 466 galaxies in and around a superstructure at $z\sim0.84$ targeted by the VIMOS Spectroscopic Survey of a Supercluster in the COSMOS field (VIS$^{3}$COS). We use [OII]$λ$3727, H$δ$, and $D_n4000$ to trace the recent, mid-, and long-term star formation histories and investigate how stellar mass and the local environment impacts those. By studying trends both in individual and composite galaxy spectra, we find that both stellar mass and environment play a role in the observed galactic properties. We find that the median [OII] equivalent width (|EW$_\mathrm{[OII]}|$) decreases from $27\pm2$ Å to $2.0_{-0.4}^{+0.5}$ Å and $D_n4000$ increases from $1.09\pm0.01$ to $1.56\pm0.03$ with increasing stellar mass (from $\sim10^{9.25}$ to $\sim10^{11.35}\ \mathrm{M_\odot}$). Concerning the dependence on the environment, we find that at fixed stellar mass |EW$_\mathrm{[OII]}|$ is tentatively lower in higher density environments. Regarding $D_n4000$, we find that the increase with stellar mass is sharper in denser environments, hinting that such environments may accelerate galaxy evolution. Moreover, we find larger $D_n4000$ values in denser environments at fixed stellar mass, suggesting that galaxies are on average older and/or more metal-rich in such dense environments. This set of tracers depicts a scenario where the most massive galaxies have, on average, the lowest sSFRs and the oldest stellar populations (age $\gtrsim1$ Gyr, showing a mass-downsizing effect). We also hypothesize that the observed increase in star formation (higher EW$_\mathrm{[OII]|}$, higher sSFR) at intermediate densities may lead to quenching since we find the quenched fraction to increase sharply from the filament to cluster-like regions at similar stellar masses.

preprint2016arXiv

Effects of Local Environment and Stellar Mass on Galaxy Quenching out to z~3

We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($ε_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $ε_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics.

preprint2016arXiv

Lyman Alpha Signatures from Direct Collapse Black Holes

`Direct collapse black holes' (DCBHs) provide possible seeds for supermassive black holes that exist at redshifts as high as z~7. We study Lyman Alpha (Lya) radiative transfer through simplified representations of the DCBH-scenario. We find that gravitational heating of the collapsing cloud gives rise to a Lya cooling luminosity of up to ~ 1e38(M_gas/1e6 Msun)^2 erg/s. The Lya production rate can be significantly larger during the final stages of collapse, but collisional deexcitation efficiently suppresses the emerging Lya flux. Photoionization by a central source boosts the Lya luminosity to L~1e43(M_BH/1e6 M_sun) erg/s during specific evolutionary stages of the cloud, where M_BH denotes the mass of the black hole powering this source. We predict that the width and velocity off-set of the Lya spectral line range from a few tens to few thousands km/s, depending sensitively on the evolutionary state of the cloud. We also compare our predictions to observations of CR7 (Sobral et al. 2015), a luminous Lya emitter at z~7, which is potentially associated with a DCBH. If CR7 is powered by a black hole, then its Lya flux alone requires that M_BH> 1e7 M_sun, which exceeds the mass of DCBHs when they first form. The observed width of the Lya spectrum favors the presence of only a low column density of hydrogen, log [N_HI/cm^-2]~19-20. The shape of the Lya spectrum indicates that this gas is outflowing. These requirements imply that if CR7 harbors a DCBH, then the physical conditions that enabled its formation have been mostly erased, which is in agreement with theoretical expectations.

preprint2016arXiv

The CALYMHA survey: Ly$α$ escape fraction and its dependence on galaxy properties at $z=2.23$

We present the first results from our CAlibrating LYMan-$α$ with H$α$ (CALYMHA) pilot survey at the Isaac Newton Telescope. We measure Ly$α$ emission for 488 H$α$ selected galaxies at $z=2.23$ from HiZELS in the COSMOS and UDS fields with a specially designed narrow-band filter ($λ_c$ = 3918 Å, $Δλ$= 52 Å). We find 17 dual H$α$-Ly$α$ emitters ($f_{\rm Lyα} >5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, of which 5 are X-ray AGN). For star-forming galaxies, we find a range of Ly$α$ escape fractions (f$_{\rm esc}$, measured with 3$"$ apertures) from $2$\%$-30$\%. These galaxies have masses from $3\times10^8$ M$_{\odot}$ to 10$^{11}$ M$_{\odot}$ and dust attenuations E$(B-V)=0-0.5$. Using stacking, we measure a median escape fraction of $1.6\pm0.5$\% ($4.0\pm1.0$\% without correcting H$α$ for dust), but show that this depends on galaxy properties. The stacked f$_{\rm esc}$ tends to decrease with increasing SFR and dust attenuation. However, at the highest masses and dust attenuations, we detect individual galaxies with f$_{\rm esc}$ much higher than the typical values from stacking, indicating significant scatter in the values of f$_{\rm esc}$. Relations between f$_{\rm esc}$ and UV slope are bimodal, with high f$_{\rm esc}$ for either the bluest or reddest galaxies. We speculate that this bimodality and large scatter in the values of f$_{\rm esc}$ is due to additional physical mechanisms such as outflows facilitating f$_{\rm esc}$ for dusty/massive systems. Ly$α$ is significantly more extended than H$α$ and the UV. f$_{\rm esc}$ continues to increase up to at least 20 kpc (3$σ$, 40 kpc [2$σ$]) for typical SFGs and thus the aperture is the most important predictor of f$_{\rm esc}$.

preprint2016arXiv

The growth of typical star-forming galaxies and their super massive black holes across cosmic time since z~2

Understanding galaxy formation and evolution requires studying the interplay between the growth of galaxies and the growth of their black holes across cosmic time. Here we explore a sample of Ha-selected star-forming galaxies from the HiZELS survey and use the wealth of multi-wavelength data in the COSMOS field (X-rays, far-infrared and radio) to study the relative growth rates between typical galaxies and their central supermassive black holes, from z=2.23 to z=0. Typical star-forming galaxies at z~1-2 have black hole accretion rates (BHARs) of 0.001-0.01 Msun/yr and star formation rates (SFRs) of ~10-40 Msun/yr, and thus grow their stellar mass much quicker than their black hole mass (~3.3 orders of magnitude faster). However, ~3% of the sample (the sources detected directly in the X-rays) show a significantly quicker growth of the black hole mass (up to 1.5 orders of magnitude quicker growth than the typical sources). BHARs fall from z=2.23 to z=0, with the decline resembling that of star formation rate density or the typical SFR. We find that the average black hole to galaxy growth (BHAR/SFR) is approximately constant for star-forming galaxies in the last 11 Gyrs. The relatively constant BHAR/SFR suggests that these two quantities evolve equivalently through cosmic time and with practically no delay between the two.

preprint2016arXiv

The KMOS Redshift One Spectroscopic Survey (KROSS): Dynamical properties, gas and dark matter fractions of typical z~1 star-forming galaxies

The KMOS Redshift One Spectroscopic Survey (KROSS) is an ESO guaranteed time survey of 795 typical star-forming galaxies in the redshift range z=0.8-1.0 with the KMOS instrument on the VLT. In this paper we present resolved kinematics and star formation rates for 584 z~1 galaxies. This constitutes the largest near-infrared Integral Field Unit survey of galaxies at z~1 to date. We demonstrate the success of our selection criteria with 90% of our targets found to be Halpha emitters, of which 81% are spatially resolved. The fraction of the resolved KROSS sample with dynamics dominated by ordered rotation is found to be 83$\pm$5%. However, when compared with local samples these are turbulent discs with high gas to baryonic mass fractions, ~35%, and the majority are consistent with being marginally unstable (Toomre Q~1). There is no strong correlation between galaxy averaged velocity dispersion and the total star formation rate, suggesting that feedback from star formation is not the origin of the elevated turbulence. We postulate that it is the ubiquity of high (likely molecular) gas fractions and the associated gravitational instabilities that drive the elevated star-formation rates in these typical z~1 galaxies, leading to the ten-fold enhanced star-formation rate density. Finally, by comparing the gas masses obtained from inverting the star-formation law with the dynamical and stellar masses, we infer an average dark matter to total mass fraction within 2.2$r_e$ (9.5kpc) of 65$\pm$12%, in agreement with the results from hydrodynamic simulations of galaxy formation.

preprint2016arXiv

The KMOS Redshift One Spectroscopic Survey (KROSS): The Tully-Fisher Relation at z ~ 1

We present the stellar mass ($M_{*}$), and K-corrected $K$-band absolute magnitude ($M_{K}$) Tully-Fisher relations (TFRs) for sub-samples of the 584 galaxies spatially resolved in H$α$ emission by the KMOS Redshift One Spectroscopic Survey (KROSS). We model the velocity field of each of the KROSS galaxies and extract a rotation velocity, $V_{80}$ at a radius equal to the major axis of an ellipse containing 80% of the total integrated H$α$ flux. The large sample size of KROSS allowed us to select 210 galaxies with well measured rotation speeds. We extract from this sample a further 56 galaxies that are rotationally supported, using the stringent criterion $V_{80}/σ> 3$, where $σ$ is the flux weighted average velocity dispersion. We find the $M_{K}$ and $M_{*}$ TFRs for this sub-sample to be $M_{K} / \rm{mag}= (-7.3 \pm 0.9) \times [(\log(V_{80}/\rm{km\ s^{-1}})-2.25]- 23.4 \pm 0.2$ , and $\log(M_{*} / M_{\odot})= (4.7 \pm 0.4) \times [(\log(V_{80}/\rm{km\ s^{-1}}) - 2.25] + 10.0 \pm 0.3$, respectively. We find an evolution of the $M_{*}$ TFR zero-point of $-0.41 \pm 0.08$ dex over the last $\sim $8 billion years. However, we measure no evolution in the $M_{K}$ TFR zero-point over the same period. We conclude that rotationally supported galaxies of a given dynamical mass had less stellar mass at $z \sim 1$ than the present day, yet emitted the same amounts of $K$-band light. The ability of KROSS to differentiate, using integral field spectroscopy with KMOS, between those galaxies that are rotationally supported and those that are not explains why our findings are at odds with previous studies without the same capabilities.

preprint2016arXiv

The Lyman-alpha luminosity function at z=5.7-6.6 and the steep drop of the faint end: implications for reionization

We present new results from the widest narrow band survey search for Lyman-alpha (Lya) emitters at z=5.7, just after reionization. We survey a total of 7 deg$^2$ spread over the COSMOS, UDS and SA22 fields. We find over 11,000 line emitters, out of which 514 are robust Lya candidates at z=5.7 within a volume of 6.3x10$^6$ Mpc$^3$. Our Lya emitters span a wide range in Lya luminosities, from faint to bright (L$_{\rm Lyα}\sim10^{42.5-44}$ erg s$^{-1}$) and rest-frame equivalent widths (EW$_0$~25-1000 Å) in a single, homogeneous data-set. By combining all our fields we find that the faint end slope of the z=5.7 Lya luminosity function is very steep, with $α=-2.3^{+0.4}_{-0.3}$. We also present an updated z=6.6 Lya luminosity function, based on comparable volumes and obtained with the same methods, which we directly compare with that at z=5.7. We find a significant decline of the number density of faint Lya emitters from z=5.7 to z=6.6 (by $0.5\pm0.1$ dex), but no evolution at the bright end/no evolution in L*. Faint Lya emitters at z=6.6 show much more extended haloes than those at z=5.7, suggesting that neutral Hydrogen plays an important role, increasing the scattering and leading to observations missing faint Lya emission within the epoch of reionization. All together, our results suggest that we are observing patchy reionization which happens first around the brightest Lya emitters, allowing the number densities of those sources to remain unaffected by the increase of neutral Hydrogen fraction from z~5 to z~7.

preprint2016arXiv

The most luminous H$α$ emitters at z~0.8-2.23 from HiZELS: evolution of AGN and star-forming galaxies

We use new near-infrared spectroscopic observations to investigate the nature and evolution of the most luminous Hα(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α}>L_{Hα}^*$, 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α}/L_{Hα}^*(z)$. While $L_{Hα}<L_{\rm Hα}^*(z)$ Ha emitters are largely dominated by star-forming galaxies (>80%), the most luminous Ha emitters ($L_{Hα}>10L_{Hα}^*(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α}>L^*_{Hα}(z)$) contribute the same fractional amount of ~15% to the total star-formation rate density, at least up to z=2.23.

preprint2016arXiv

The nature of H-alpha star-forming galaxies at z~0.4 in and around Cl 0939+4713: the environment matters

Cluster star-forming galaxies are found to have an excess of Far-Infrared emission relative to H-alpha (Ha), when compared to those in the field, which could be caused by intense AGN activity, dust and/or declining star formation histories. Here we present spectroscopic observations of Ha emitters in the Cl 0939+4713 (Abell 851) super-cluster at z=0.41, using AF2+WYFFOS on the WHT. We measure [OII], Hbeta (Hb), [OIII], Ha and [NII] for a sample of 119 Ha emitters in and around the cluster. We find that 17+-5% of the Ha emitters are AGN, irrespective of environment. For star-forming galaxies, we obtain Balmer decrements, metallicities and ionisation parameters with different methods, individually and by stacking. We find a strong mass-metallicity relation at all environments, with no significant dependence on environment. The ionisation parameter declines with increasing stellar mass for low-mass galaxies. Ha emitters residing in intermediate environments show the highest ionisation parameters (along with high [OIII]/Ha and high [OIII]/[OII] line ratios, typically twice as large as in the highest and lowest densities), which decline with increasing environmental density. Dust extinction (A$_{Hα}$) correlates strongly with stellar mass, but also with environmental density. Star-forming galaxies in the densest environments are found to be significantly dustier (A$_{Hα}$~1.5-1.6) than those residing in the lowest density environments (A$_{Hα}$~0.6), deviating significantly from what would be predicted given their stellar masses.

preprint2016arXiv

The Nature of H$β+$[O{\sc iii}] and [O{\sc ii}] emitters to $z \sim 5$ with HiZELS: stellar mass functions and the evolution of EWs

We investigate the properties of $\sim7000$ narrow-band selected galaxies with strong H$β+$[OIII] and [OII] nebular emission lines from the High-$z$ Emission Line Survey (HiZELS) between $z \sim 0.8 - 5.0$. Our sample covers a wide range in stellar mass ($M\sim10^{7.5 - 12.0}$ M$_\odot$), rest-frame equivalent widths (EW$_\mathrm{rest}\sim 10 - 10^5$ Å), and line luminosities ($L\sim10^{40.5 - 43.2}$ erg s$^{-1}$). We measure the H$β+$[OIII] and [OII]-selected stellar mass functions out to $z \sim 3.5$ and find that both $M_\star$ and $ϕ_\star$ increases with cosmic time, which may be due to the [OIII] selection including an increasing fraction of AGN at lower redshifts. The [OII]-selected stellar mass functions show a constant $M_\star\sim10^{11.6}$ M$_\odot$ and a strong, increasing evolution with cosmic time in $ϕ_\star$ in line with H$α$ studies. We also investigate the EW$_\mathrm{rest}$ evolution as a function of redshift with a fixed mass range (10$^{9.5 - 10.0}$ M$_\odot$) and find an increasing trend best represented by $(1+z)^{3.81\pm0.14}$ and $(1+z)^{2.72\pm0.19}$ up to $z\sim 2$ and $\sim 3$ for H$β$+[OIII] and [OII] emitters, respectively. This is the first time that the EW evolution has been directly measured for H$β+$[OIII] and [OII] emitters up to these redshifts. There is evidence for a slower evolution for $z > 2$ in the H$β+$[OIII] EW and a decreasing trend for $z > 3$ in the [OII] EW evolution, which would imply low [OII] EW at the highest redshifts and higher [OIII]/[OII] line ratios. This suggests that the ionization parameter at higher redshift may be significantly higher than the local Universe. Our results set the stage for future near-IR space-based spectroscopic surveys to test our extrapolated predictions and also produce $z > 5$ measurements to constrain the high-$z$ end of the EW$_\mathrm{rest}$ and [OIII]/[OII] evolution.

preprint2016arXiv

The production and escape of Lyman-Continuum radiation from star-forming galaxies at z~2 and their redshift evolution

We study the production rate of ionizing photons of a sample of 588 H$α$ emitters (HAEs) and 160 Lyman-$α$ emitters (LAEs) at $z=2.2$ in the COSMOS field in order to assess the implied emissivity from galaxies, based on their UV luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/$NUV$ data, we find f$_{\rm esc} < 2.8\, (6.4)$% through median (mean) stacking. By combining the H$α$ luminosity density with IGM emissivity measurements from absorption studies, we find a globally averaged $\langle$f$_{\rm esc}\rangle$ of $5.9^{+14.5}_{-4.2}$ % at $z=2.2$ if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global $\langle$f$_{\rm esc}\rangle$ at $z\approx3-5$, also ruling out a high $\langle$f$_{\rm esc}\rangle$ at $z<5$. These low escape fractions allow us to measure $ξ_{ion}$, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical $ξ_{ion} \approx 10^{24.77\pm0.04}$ Hz erg$^{-1}$ for HAEs and $ξ_{ion} \approx 10^{25.14\pm0.09}$ Hz erg$^{-1}$ for LAEs. LAEs and low mass HAEs at $z=2.2$ show similar values of $ξ_{ion}$ as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing $ξ_{ion}$ with increasing EW(H$α$), $ξ_{ion}$ likely increases with redshift. This evolution alone is fully in line with the observed evolution of $ξ_{ion}$ between $z\approx2-5$, indicating a typical value of $ξ_{ion} \approx 10^{25.4}$ Hz erg$^{-1}$ in the reionization era.

preprint2015arXiv

A Comparative Study of Density Field Estimation for Galaxies: New Insights into the Evolution of Galaxies with Environment in COSMOS out to z~3

It is well-known that galaxy environment has a fundamental effect in shaping its properties. We study the environmental effects on galaxy evolution, with an emphasis on the environment defined as the local number density of galaxies. The density field is estimated with different estimators (weighted adaptive kernel smoothing, 10$^{th}$ and 5$^{th}$ nearest neighbors, Voronoi and Delaunay tessellation) for a K$_{s}<$24 sample of $\sim$190,000 galaxies in the COSMOS field at 0.1$<$z$<$3.1. The performance of each estimator is evaluated with extensive simulations. We show that overall, there is a good agreement between the estimated density fields using different methods over $\sim$2 dex in overdensity values. However, our simulations show that adaptive kernel and Voronoi tessellation outperform other methods. Using the Voronoi tessellation method, we assign surface densities to a mass complete sample of quiescent and star-forming galaxies out to z$\sim$3. We show that at a fixed stellar mass, the median color of quiescent galaxies does not depend on their host environment out to z$\sim$3. We find that the number and stellar mass density of massive ($>$10$^{11}$M$_{\odot}$) star-forming galaxies have not significantly changed since z$\sim$3, regardless of their environment. However, for massive quiescent systems at lower redshifts (z$\lesssim$1.3), we find a significant evolution in the number and stellar mass densities in denser environments compared to lower density regions. Our results suggest that the relation between stellar mass and local density is more fundamental than the color-density relation and that environment plays a significant role in quenching star formation activity in galaxies at z$\lesssim$1.

preprint2015arXiv

A large narrow band H$α$ survey at $z\sim0.2$: the bright end of the luminosity function, cosmic variance and clustering across cosmic time

We carried out the largest ($>3.5\times10^5$ Mpc$^3$, 26 deg$^2$) H$α$ narrow band survey to date at $z\sim0.2$ in the SA22, W2 and XMMLSS extragalactic fields. Our survey covers a large enough volume to overcome cosmic variance and to sample bright and rare H$α$ emitters up to an observed luminosity of $\sim10^{42.4}$ erg s$^{-1}$, equivalent to $\sim11 M_\odot$ yr$^{-1}$. Using our sample of $220$ sources brighter than $>10^{41.4}$ erg s$^{-1}$ ($>1 M_\odot$ yr$^{-1}$), we derive H$α$ luminosity functions, which are well described by a Schechter function with $ϕ^* = 10^{-2.85\pm0.03}$ Mpc$^{-3}$ and $L^*_{Hα} = 10^{41.71\pm0.02}$ erg s$^{-1}$ (with a fixed faint end slope $α=-1.35$). We find that surveys probing smaller volumes ($\sim3\times10^4$ Mpc$^3$) are heavily affected by cosmic variance, which can lead to errors of over $100$ per cent in the characteristic density and luminosity of the H$α$ luminosity function. We derive a star formation rate density of $ρ_\mathrm{SFRD} = 0.0094\pm0.0008$ $M_\odot$ yr$^{-1}$, in agreement with the redshift-dependent H$α$ parametrisation from Sobral et al. (2013). The two-point correlation function is described by a single power law $ω(θ) = (0.159\pm0.012) θ^{(-0.75\pm0.05)}$, corresponding to a clustering length of $r_0 = 3.3\pm0.8$ Mpc/h. We find that the most luminous H$α$ emitters at $z\sim0.2$ are more strongly clustered than the relatively fainter ones. The $L^*_{Hα}$ H$α$ emitters at $z\sim0.2$ in our sample reside in $\sim10^{12.5-13.5}$ $M_\odot$ dark matter haloes. This implies that the most star forming galaxies always reside in relatively massive haloes or group-like environments and that the typical host halo mass of star-forming galaxies is independent of redshift if scaled by $L_\mathrm{Hα}/L^*_{Hα}(z)$, as proposed by Sobral et al. (2010).

preprint2015arXiv

CF-HiZELS, a 10 deg$^2$ emission-line survey with spectroscopic follow-up: Hα, [OIII]+Hβ and [OII] luminosity functions at z=0.8, 1.4 and 2.2

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 $μ$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$α$ emitters at z=0.8, 1343 [OIII]+H$β$ 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$σ$ (spectroscopically confirmed) over-density of H$α$ emitters at z=0.81. We derive H$α$, [OIII]+H$β$ 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 $ϕ^*$ 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.

preprint2015arXiv

Evidence for PopIII-like stellar populations in the most luminous Lyman-$α$ emitters at the epoch of re-ionisation: spectroscopic confirmation

Faint Lyman-$α$ (Ly$α$) emitters become increasingly rarer towards the re-ionisation epoch (z~6-7). However, observations from a very large (~5deg$^2$) Ly$α$ survey at z=6.6 (Matthee et al. 2015) show that this is not the case for the most luminous emitters. Here we present follow-up observations of the two most luminous z~6.6 Ly$α$ candidates in the COSMOS field: `MASOSA' and `CR7'. We used X-SHOOTER, SINFONI and FORS2 (VLT), and DEIMOS (Keck), to confirm both candidates beyond any doubt. We find redshifts of z=6.541 and z=6.604 for MASOSA and CR7, respectively. MASOSA has a strong detection in Ly$α$ with a line width of $386\pm30$ km/s (FWHM) and with high EW$_0$ (>200 Å), but it is undetected in the continuum. CR7, with an observed Ly$α$ luminosity of $10^{43.93\pm0.05}$erg/s is the most luminous Ly$α$ emitter ever found at z>6. CR7 reveals a narrow Ly$α$ line with $266\pm15$ km/s FWHM, being detected in the NIR (rest-frame UV, with $β=-2.3\pm0.1$) with an excess in $J$, and also strongly detected in IRAC/Spitzer. We detect a narrow HeII1640$Å$ emission line ($6σ$) which explains the excess seen in the $J$ band photometry (EW$_0$~80 Å). We find no other emission lines from the UV to the NIR in our X-SHOOTER spectra (HeII/OIII]1663A>3 and HeII/CIII]1908A>2.5). We find that CR7 is best explained by a combination of a PopIII-like population which dominates the rest-frame UV and the nebular emission, and a more normal stellar population which dominates the mass. HST/WFC3 observations show that the light is indeed spatially separated between a very blue component, coincident with Ly$α$ and HeII emission, and two red components (~5 kpc away), which dominate the mass. Our findings are consistent with theoretical predictions of a PopIII wave, with PopIII star formation migrating away from the original sites of star formation.

preprint2015arXiv

Evolution of the H$β$+[OIII] and [OII] luminosity functions and the [OII] star-formation history of the Universe up to $z$ ~ 5 from HiZELS

We investigate the evolution of the H$β$+[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$β$+[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$β$+[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 $ϕ_\star$ for both H$β$+[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$β$+[OIII] SFRDs to the [OII] and H$α$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$β$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs.

preprint2015arXiv

Identification of the brightest Lyα emitters at z=6.6: implications for the evolution of the luminosity function in the re-ionisation era

Using wide field narrow-band surveys, we provide a new measurement of the $z=6.6$ Lyman-$α$ Emitter (LAE) luminosity function (LF), which constraints the bright end for the first time. We use a combination of archival narrow-band NB921 data in UDS and new NB921 measurements in SA22 and COSMOS/UltraVISTA, all observed with the Subaru telescope, with a total area of $\sim 5$ deg$^2$. We exclude lower redshift interlopers by using broad-band optical and near-infrared photometry and also exclude three supernovae with data split over multiple epochs. Combining the UDS and COSMOS samples we find no evolution of the bright end of the Ly$α$ LF between $z=5.7$ and $6.6$, which is supported by spectroscopic follow-up, and conclude that sources with \emph{Himiko}-like luminosity are not as rare as previously thought, with number densities of $\sim 1.5\times10^{-5}$ Mpc$^{-3}$. Combined with our wide-field SA22 measurements, our results indicate a non-Schechter-like bright end of the LF at $z=6.6$ and a different evolution of \emph{observed} faint and bright LAEs. This differential evolution is also seen in the spectroscopic follow-up of UV selected galaxies and is now also confirmed for Ly$α$ emitters, and we argue that it may be an effect of re-ionisation. Using a toy-model, we show that such differential evolution of the LF is expected, since brighter sources are able to ionise their surroundings earlier, such that Ly$α$ photons are able to escape. Our targets are excellent candidates for detailed follow-up studies and provide the possibility to give a unique view on the earliest stages in the formation of galaxies and re-ionisation process.

preprint2015arXiv

MC$^2$: Boosted AGN and star-formation activity in CIZA J2242.8+5301, a massive post-merger cluster at z=0.19

Cluster mergers may play a fundamental role in the formation and evolution of cluster galaxies. Stroe et al. (2014) revealed unexpected over-densities of candidate H$α$ emitters near the ~1 Mpc-wide shock fronts of the massive (~2x10$^{15}$M$_{\odot}$) "Sausage" merging cluster, CIZA J2242.8+5301. We used Keck/DEIMOS and WHT/AF2 to confirm 83 H$α$ emitters in and around the merging cluster. We find that cluster star-forming galaxies in the hottest X-ray gas and/or in the cluster sub-cores (away from the shock fronts) show high [SII]6716/[SII]6761 and high [SII]6716/H$α$, implying very low electron densities (<30x lower than all other star-forming galaxies outside the cluster) and significant contribution from supernovae, respectively. All cluster star-forming galaxies near the cluster centre show evidence of significant outflows (blueshifted Na D~200-300km/s), likely driven by supernovae. Strong outflows are also found for the cluster H$α$ AGN. H$α$ star-forming galaxies in the merging cluster follow the z~0 mass-metallicity relation, showing systematically higher metallicity (~0.15-0.2 dex) than H$α$ emitters outside the cluster (projected R>2.5 Mpc). This suggests that the shock front may have triggered remaining metal-rich gas which galaxies were able to retain into forming stars. Our observations show that the merger of impressively massive (~10$^{15}$M$_\odot$) clusters can provide the conditions for significant star-formation and AGN activity, but, as we witness strong feedback by star-forming galaxies and AGN (and given how massive the merging cluster is), such sources will likely quench in a few 100 Myrs.

preprint2015arXiv

Nebular and Stellar Dust Extinction Across the Disk of Emission-Line Galaxies on Small (KPC) Scales

We investigate resolved kpc-scale stellar and nebular dust distribution in eight star-forming galaxies at z~0.4 in the GOODS fields. Constructing the observed Spectral Energy Distributions (SEDs) per pixel, based on seven bands photometric data from HST/ACS and WFC3, we performed pixel-by-pixel SED fits to population synthesis models and estimated small-scale distribution of stellar dust extinction. We use Halpha / Hbeta nebular emission line ratios from Keck/DEIMOS high resolution spectra at each spatial resolution element to measure the amount of attenuation faced by ionized gas at different radii from the center of galaxies. We find a good agreement between the integrated and median of resolved color excess measurements in our galaxies. The ratio of integrated nebular to stellar dust extinction is always greater than unity, but does not show any trend with stellar mass or star formation rate. We find that inclination plays an important role in the variation of the nebular to stellar excess ratio. The stellar color excess profiles are found to have higher values at the center compared to outer parts of the disk. However, for lower mass galaxies, a similar trend is not found for the nebular color excess. We find that the nebular color excess increases with stellar mass surface density. This explains the absence of radial trend in the nebular color excess in lower mass galaxies which lack a large radial variation of stellar mass surface density. Using standard conversions of star formation rate surface density to gas mass surface density, and the relation between dust mass surface density and color excess, we find no significant variation in the dust to gas ratio in regions with high gas mass surface densities, over the scales probed in this study.

preprint2015arXiv

Neutral hydrogen gas, past and future star-formation in galaxies in and around the 'Sausage' merging galaxy cluster

CIZA J2242.8+5301 ($z = 0.188$, nicknamed 'Sausage') is an extremely massive ($M_{200}\sim 2.0 \times 10^{15}M_\odot$ ), merging cluster with shock waves towards its outskirts, which was found to host numerous emission-line galaxies. We performed extremely deep Westerbork Synthesis Radio Telescope HI observations of the 'Sausage' cluster to investigate the effect of the merger and the shocks on the gas reservoirs fuelling present and future star formation (SF) in cluster members. By using spectral stacking, we find that the emission-line galaxies in the 'Sausage' cluster have, on average, as much HI gas as field galaxies (when accounting for the fact cluster galaxies are more massive than the field galaxies), contrary to previous studies. Since the cluster galaxies are more massive than the field spirals, they may have been able to retain their gas during the cluster merger. The large HI reservoirs are expected to be consumed within $\sim0.75-1.0$ Gyr by the vigorous SF and AGN activity and/or driven out by the out-flows we observe. We find that the star-formation rate in a large fraction of H$α$ emission-line cluster galaxies correlates well with the radio broad band emission, tracing supernova remnant emission. This suggests that the cluster galaxies, all located in post-shock regions, may have been undergoing sustained SFR for at least 100 Myr. This fully supports the interpretation proposed by Stroe et al. (2015) and Sobral et al. (2015) that gas-rich cluster galaxies have been triggered to form stars by the passage of the shock.

preprint2015arXiv

Spectroscopic Study of Star-forming Galaxies in Filaments and the Field at $z\sim$0.5: Evidence for Environmental Dependence of Electron Density

We study the physical properties of a spectroscopic sample of 28 star-forming galaxies in a large filamentary structure in the COSMOS field at $z\sim$0.53, with spectroscopic data taken with the Keck/DEIMOS spectrograph, and compare them with a control sample of 30 field galaxies. We spectroscopically confirm the presence of a large galaxy filament ($\sim$ 8 Mpc), along which five confirmed X-ray groups exist. We show that within the uncertainties, the ionization parameter, equivalent width (EW), EW versus specific star-formation rate (sSFR) relation, EW versus stellar mass relation, line-of-sight velocity dispersion, dynamical mass, and stellar-to-dynamical mass ratio are similar for filament and field star-forming galaxies. However, we show that on average, filament star-forming galaxies are more metal-enriched ($\sim$ 0.1$-$0.15 dex), possibly due to the inflow of the already enriched intrafilamentary gas into filament galaxies. Moreover, we show that electron densities are significantly lower (a factor of $\sim$17) in filament star-forming systems compared to those in the field, possibly because of a longer star-formation timescale for filament star-forming galaxies. Our results highlight the potential pre-processing role of galaxy filaments and intermediate-density environments on the evolution of galaxies, which has been highly underestimated.

preprint2014arXiv

A 10 deg$^2$ Lyman-$α$ survey at z=8.8 with spectroscopic follow-up: strong constraints on the LF and implications for other surveys

Candidate galaxies at redshifts of $z \sim 10$ are now being found in extremely deep surveys, probing very small areas. As a consequence, candidates are very faint, making spectroscopic confirmation practically impossible. In order to overcome such limitations, we have undertaken the CF-HiZELS survey, which is a large area, medium depth near infrared narrow-band survey targeted at $z=8.8$ Lyman-$α$ (Ly$α$) emitters (LAEs) and covering 10 deg$^2$ in part of the SSA22 field with the Canada-France-Hawaii Telescope. We surveyed a comoving volume of $4.7\times 10^6$ Mpc$^3$ to a Ly$α$ luminosity limit of $6.3\times10^{43}$ erg s$^{-1}$. We look for Ly$α$ candidates by applying the following criteria: i) clear emission line source, ii) no optical detections ($ugriz$ from CFHTLS), iii) no visible detection in the optical stack ($ugriz > 27$), iv) visually checked reliable NB$_J$ and $J$ detections and v) $J-K \leq 0$. We compute photometric redshifts and remove a significant amount of dusty lower redshift line-emitters at $z \sim 1.4 $ or $2.2$. A total of 13 Ly$α$ candidates were found, of which two are marked as strong candidates, but the majority have very weak constraints on their SEDs. Using follow-up observations with SINFONI/VLT we are able to exclude the most robust candidates as Ly$α$ emitters. We put a strong constraint on the Ly$α$ luminosity function at $z \sim 9$ and make realistic predictions for ongoing and future surveys. Our results show that surveys for the highest redshift LAEs are susceptible of multiple contaminations and that spectroscopic follow-up is absolutely necessary.

preprint2014arXiv

A relationship between specific star formation rate and metallicity gradient within z=1 galaxies from KMOS-HiZELS

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.

preprint2014arXiv

Cosmic Web and Star Formation Activity in Galaxies at z~1

We investigate the role of the delineated cosmic web/filaments on the star formation activity by exploring a sample of 425 narrow-band selected Hα emitters, as well as 2846 color-color selected underlying star-forming galaxies for a large scale structure (LSS) at z=0.84 in the COSMOS field from the HiZELS survey. Using the scale-independent Multi-scale Morphology Filter (MMF) algorithm, we are able to quantitatively describe the density field and disentangle it into its major components: fields, filaments and clusters. We show that the observed median star formation rate (SFR), stellar mass, specific star formation rate (sSFR), the mean SFR-Mass relation and its scatter for both Hα emitters and underlying star-forming galaxies do not strongly depend on different classes of environment, in agreement with previous studies. However, the fraction of Hα emitters varies with environment and is enhanced in filamentary structures at z~1. We propose mild galaxy-galaxy interactions as the possible physical agent for the elevation of the fraction of Hα star-forming galaxies in filaments. Our results show that filaments are the likely physical environments which are often classed as the "intermediate" densities, and that the cosmic web likely plays a major role in galaxy formation and evolution which has so far been poorly investigated.

preprint2014arXiv

Mapping the large scale structure around a z=1.46 galaxy cluster in 3-D using two adjacent narrow-band filters

We present a novel method to estimate accurate redshifts of star-forming galaxies by measuring the flux ratio of the same emission line observed through two adjacent narrow-band filters. We apply this method to our NB912 and new NB921 data taken with Suprime-Cam on the Subaru Telescope of a galaxy cluster, XMMXCS J2215.9-1738, at z=1.46 and its surrounding structures. We obtain redshifts for 170 [OII] emission line galaxies at z~1.46, among which 41 galaxies are spectroscopically confirmed with MOIRCS and FMOS on the Subaru mainly, showing an accuracy of $σ$((z-z_spec})/(1+z_spec))=0.002. This allows us to reveal filamentary structures that penetrate towards the centre of the galaxy cluster and intersect with other structures, consistent with the picture of hierarchical cluster formation. We also find that the projected celestial distribution does not precisely trace the real distribution of galaxies, indicating the importance of the three dimensional view of structures to properly identify and quantify galaxy environments. We investigate the environmental dependence of galaxy properties with local density, confirming that the median colour of galaxies becomes redder in higher density region while the star-formation rate of star-forming galaxies does not depend strongly on local environment in this structure. This implies that the star-forming activity in galaxies is truncated on a relatively short time scale in the cluster centre.

preprint2014arXiv

MC^2: Constraining the Dark Matter Distribution of the Violent Merging Galaxy Cluster CIZA J2242.8+5301 by Piercing through the Milky Way

CIZA J2242.8+5301 is a merging system with a prominent (~2 Mpc long) radio relic, which together with the morphology of the X-ray emission provides strong evidence for a violent collision along the N-S axis. We present our constraints on the dark matter distribution of this unusual system using Subaru and CFHT imaging data. Measuring a high S/N lensing signal from this cluster is potentially a challenging task because of its proximity to the Milky Way plane (|b|~5 deg). We overcome this challenge with careful observation planning and systematics control, which enables us to successfully map the dark matter distribution of the cluster with high fidelity. The resulting mass map shows that the mass distribution is highly elongated along the N-S merger axis inferred from the orientation of the radio relics. Based on our mass reconstruction, we identify two sub-clusters, which coincide with the cluster galaxy distributions. We determine their masses using MCMC analysis by simultaneously fitting two NFW halos without fixing their centroids. The resulting masses of the northern and southern systems are $M_{200}=11.0_{-3.2}^{+3.7}\times10^{14} M_{\sun}$ and $9.8_{-2.5}^{+3.8}\times10^{14} M_{\sun}$, respectively, indicating that we are witnessing a post-collision of two giant systems of nearly equal mass. When the mass and galaxy centroids are compared in detail, we detect ~ 1' (~190 kpc) offsets in both northern and southern sub-clusters. We find that the galaxy luminosity-mass offset for the northern clump is statistically significant at the ~2 sigma level whereas the detection is only marginal for the southern sub-cluster in part because of a relatively large mass centroid error. We conclude that it is yet premature to uniquely attribute the galaxy-mass misalignment to SIDM and discuss caveats.

preprint2014arXiv

MC^2: Galaxy Imaging and Redshift Analysis of the Merging Cluster CIZA J2242.8+5301

X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band) as well as a comprehensive spectroscopic survey of the cluster area (505 redshifts) using Keck DEIMOS. We use this data to perform a comprehensive galaxy/redshift analysis of the system, which is the first step to a proper understanding the geometry and dynamics of the merger, as well as using the merger to constrain self-interacting dark matter. We find that the system is dominated by two subclusters of comparable richness with a projected separation of 6.9'^{+0.7}_{-0.5} (1.3^{+0.13}_{-0.10} Mpc). We find that the north and south subclusters have similar redshifts of z=0.188 with a relative line-of-sight velocity difference of 69+/-190 km/s. We also find that north and south subclusters have velocity dispersions of 1160^{+100}_{-90} km/s and 1080^{+100}_{-70} km/s, respectively. These correspond to masses of 16.1^{+4.6}_{-3.3}x10^14 M_sun and 13.0^{+4.0}_{-2.5}x10^14 M_sun, respectively. While velocity dispersion measurements of merging clusters can be biased we believe the bias in this system to be minor due to the large projected separation and nearly plane-of-sky merger configuration. CIZA J2242.8+5301 is a relatively clean dissociative cluster merger with near 1:1 mass ratio, which makes it an ideal merger for studying merger associated physical phenomena.

preprint2014arXiv

Spectral and morphological signatures of cluster merger shocks: CIZA J2242.8+5301

The CIZA J2242.8+5301 cluster hosts an extraordinary, narrow, Mpc-wide, diffuse patch of radio emission tracing travelling shock waves, called a relic. We perform a multi-wavelength of observations of the cluster. We discover radio spectral steepening and increasing spectral curvature in the shock downstream area, indicative of emission produced by spectrally-aged electrons, following a JP ageing model. The curved integrated spectrum towards $16$ GHz challenges the favoured relic formation model. Further, we find a boosting the number density of H$α$ emitters by a factor of $10$, possibly caused by interactions of the galaxies with the shock wave.

preprint2014arXiv

The highest-frequency detection of a radio relic: 16-GHz AMI observations of the `Sausage' cluster

We observed the cluster CIZA J2242.8+5301 with the Arcminute Microkelvin Imager at $16$ GHz and present the first high radio-frequency detection of diffuse, non-thermal cluster emission. This cluster hosts a variety of bright, extended, steep-spectrum synchrotron-emitting radio sources, associated with the intra-cluster medium, called radio relics. Most notably, the northern, Mpc-wide, narrow relic provides strong evidence for diffusive shock acceleration in clusters. We detect a puzzling, flat-spectrum, diffuse extension of the southern relic, which is not visible in the lower radio-frequency maps. The northern radio relic is unequivocally detected and measures an integrated flux of $1.2\pm0.3$ mJy. While the low-frequency ($<2$ GHz) spectrum of the northern relic is well represented by a power-law, it clearly steepens towards $16$ GHz. This result is inconsistent with diffusive shock acceleration predictions of ageing plasma behind a uniform shock front. The steepening could be caused by an inhomogeneous medium with temperature/density gradients or by lower acceleration efficiencies of high energy electrons. Further modelling is necessary to explain the observed spectrum.

preprint2014arXiv

The rise and fall of star-formation in $\bf z\sim0.2$ merging galaxy clusters

CIZA J2242.8+5301 (`Sausage') and 1RXS J0603.3+4213 (`Toothbrush') are two low-redshift ($z\sim0.2$), massive ($\sim2\times10^{15}M_\odot$), post-core passage merging clusters, which host shock waves traced by diffuse radio emission. To study their star-formation properties, we uniformly survey the `Sausage' and `Toothbrush' clusters in broad and narrow band filters and select a sample of $201$ and $463$ line emitters, down to a rest-frame equivalent width ($13$Å). We robustly separate between H$α$ and higher redshift emitters using a combination of optical multi-band (B, g, V, r, i, z) and spectroscopic data. We build H$α$ luminosity functions for the entire cluster region, near the shock fronts, and away from the shock fronts and find striking differences between the two clusters. In the dynamically younger, $1$ Gyr old `Sausage' cluster we find numerous ($59$) H$α$ emitters above a star-formation rate (SFR) of $0.17$ M$_{\sun}$ yr$^{-1}$ surprisingly located in close proximity to the shock fronts, embedded in very hot intra-cluster medium plasma. The SFR density for the cluster population is at least at the level of typical galaxies at $z\sim2$. Down to the same star-formation rate, the possibly dynamically more evolved `Toothbrush' cluster has only $9$ H$α$ galaxies. The cluster H$α$ galaxies fall on the SFR-stellar mass relation $z\sim0.2$ for the field. However, the `Sausage' cluster has an H$α$ emitter density $>20$ times that of blank fields. If the shock passes through gas-rich cluster galaxies, the compressed gas could collapse into dense clouds and excite star-formation for a few $100$ Myr. This process ultimately leads to a rapid consumption of the molecular gas, accelerating the transformation of gas-rich field spirals into cluster S0s or ellipticals.

preprint2013arXiv

A fundamental metallicity relation for galaxies at z = 0.84 - 1.47 from HiZELS

We obtained Subaru FMOS observations of Halpha emitting galaxies selected from the HiZELS narrow-band survey, to investigate the relationship between stellar mass, metallicity and star-formation rate at z = 0.84 - 1.47, for comparison with the Fundamental Metallicity Relation seen at low redshift. Our findings demonstrate, for the first time with a homogeneously selected sample, that a relationship exists for typical star-forming galaxies at z = 1 - 1.5 and that it is surprisingly similar to that seen locally. Therefore, star-forming galaxies at z = 1 - 1.5 are no less metal abundant than galaxies of similar mass and star formation rate (SFR) at z = 0.1, contrary to claims from some earlier studies. We conclude that the bulk of the metal enrichment for this star-forming galaxy population takes place in the 4 Gyr before z = 1.5. We fit a new mass-metallicity-SFR plane to our data which is consistent with other high redshift studies. However, there is some evidence that the mass-metallicity component of this high redshift plane is flattened, at all SFR, compared with z = 0.1, suggesting that processes such as star-formation driven winds, thought to remove enriched gas from low mass halos, are yet to have as large an impact at this early epoch. The negative slope of the SFR-metallicity relation from this new plane is consistent with the picture that the elevation in the SFR of typical galaxies at z > 1 is fuelled by the inflow of metal-poor gas and not major merging.

preprint2013arXiv

On the evolution and environmental dependence of the star formation rate versus stellar mass relation since z~2

This paper discusses the evolution of the correlation between galaxy star formation rates (SFRs) and stellar mass (M*) over the last ~10 Gyrs, particularly focusing on its environmental dependence. We first present the mid-infrared (MIR) properties of the H-alpha-selected galaxies in a rich cluster Cl0939+4713 at z=0.4. With a MIR stacking analysis, we find that the median SFR of H-alpha emitters is higher in higher-density environment at z=0.4. We also find that star-forming galaxies in high-density environment tend to have higher specific SFR (SSFR), although the trend is much less significant compared to that of SFR. This increase of SSFR in high-density environment is not visible when we consider the SFR derived from H-alpha alone, suggesting that the dust attenuation in galaxies depends on environment; galaxies in high-density environment tend to be dustier (by up to ~0.5 mag). We then discuss the environmental dependence of the SFR-M* relation for star-forming galaxies since z~2, by compiling our comparable, narrow-band selected, large H-alpha emitter samples in both distant cluster environments (from MAHALO-Subaru) and field environments (from HiZELS). We find that the SSFR of H-alpha-selected galaxies (at the fixed mass of log(M*/Msun)=10) rapidly evolves as (1+z)^3, but the SFR-M* relation is independent of the environment since z~2, as far as we rely on the H-alpha-based SFRs (with M*-dependent extinction correction). Even if we consider the possible environmental variation in the dust attenuation, we conclude that the difference in the SFR-M* relation between cluster and field star-forming galaxies is always small (~<0.2dex level) at any time in the history of the Universe since z~2. [abridged]

preprint2013arXiv

The merger rates and sizes of galaxies across the peak epoch of star formation from the HiZELS survey

We use the HiZELS narrow-band H-alpha survey in combination with CANDELS, UKIDSS and WIRDS near-infrared imaging, to investigate the morphologies, merger rates and sizes of a sample of H-alpha emitting galaxies in the redshift range z=0.40 - 2.23, an epoch encompassing the rise to the peak of the star formation rate density. Merger rates are estimated from space- and ground-based imaging using the M20 coefficient. To account for the increase in the specific star-formation rate (sSFR) of the star forming `main-sequence' with redshift, we normalise the star-formation rate of galaxies at each epoch to the typical value derived from the H-alpha luminosity function. Once this trend in sSFR is removed we see no evidence for an increase in the number density of star-forming galaxies or the merger rate with redshift. We thus conclude that neither is the main driver of the enhanced star-formation rate density at z=1-2, with secular processes such as instabilities within efficiently fuelled, gas-rich discs or multiple minor mergers the most likely alternatives. However, we find that 40-50% of starburst galaxies, those with enhanced specific star formation at their epoch, are major mergers and this fraction is redshift independent. Finally, we find the surprising result that the typical size of a star-forming galaxy of a given mass does not evolve across the redshift range considered, suggesting a universal size-mass relation. Taken in combination, these results indicate a star-forming galaxy population that is statistically similar in physical size, merger rate and mass over the 6 Gyr covered in this study, despite the increase in typical sSFR.

preprint2013arXiv

The role of cluster mergers and travelling shocks in shaping the H$α$ luminosity function at $\bf z\sim0.2$: `sausage' and `toothbrush' clusters

The most extreme cluster mergers can lead to massive cluster-wide travelling shock waves. The CIZA J2242.8+5301 ('sausage') and 1RXS J0603.3+4213 (`toothbrush') clusters ($z\sim0.2$) host enormous radio-emitting shocks with simple geometry. We investigate the role of mergers and shocks in shaping the H$α$ luminosity function, using custom-made narrow-band filters matching the cluster redshifts mounted on the INT. We surveyed $\sim0.28$ deg$^2$ for each cluster and found $181$ line emitters in the `sausage' (volume of $3.371\times10^3$ Mpc$^3$ for H$α$ at $z=0.1945$) and $141$ in the `toothbrush' ($4.546\times10^3$ Mpc$^3$ for H$α$ at $z=0.225$), out of which $49$ (`sausage') and $30$ (`toothbrush') are expected to be H$α$. We build luminosity functions for the field-of-view down to an average limiting star formation rate of $0.14$ M$_{\odot}$ yr$^{-1}$, find good agreement with field luminosity functions at $z=0.2$, but significant differences between the shapes of the luminosity functions for the two clusters. We discover extended, tens-of-kpc-wide H$α$ haloes in galaxies neighbouring relics, which were possibly disrupted by the passage of the shock wave. By comparing the `sausage' cluster with blank fields and other clusters, we also uncover an order of magnitude boost (at $9σ$ level) in the normalisation $ϕ^*$ of the luminosity function in the relic areas. Our results suggest that cluster mergers may play an important role in the evolution of cluster galaxies through shock-induced star formation.

preprint2013arXiv

The stellar mass function of star-forming galaxies and the mass-dependent SFR function since z=2.23 from HiZELS

We explore a large uniformly selected sample of H$α$ 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, $α\approx$-1.6. The stellar mass function of SFGs, however, reveals little evolution: $α \approx$-1.4, M$^*$$\approx$10$^{11.2\pm0.2}$ M$_{\odot}$ and just a slight increase of $\approx$2.3x in $Φ^*$ 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 ($ρ_{\rm SFR}$) per dlog$_{10}$M, and that there is no significant evolution in the fractional contribution from SFGs of different masses to $ρ_{\rm SFR}$ or $ρ_{\rm SFR}$(dlog$_{10}$M)$^{-1}$ since z=2.23. Instead, we show that the decline of SFR$^*$ and of $ρ_{\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.

preprint2012arXiv

A large Hα survey at z=2.23, 1.47, 0.84 & 0.40: the 11 Gyr evolution of star-forming galaxies from HiZELS

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).

preprint2012arXiv

Calibrating [OII] star-formation rates at z>1 from dual Hα-[OII] imaging from HiZELS

We investigate the relationship between Hα and [OII](λ3727) emission in faint star-forming galaxies at z=1.47 with dust uncorrected star formation rates (SFRs) down to 1.4 Msun/yr, using data in two narrow-bands from WFCAM/UKIRT and Suprime-Cam/Subaru. A stacking analysis allows us to investigate Hα emission flux from bright [OII] emitters as well as faint ones for which Hα is not individually detected, and to compare them with a large sample of local galaxies. We find that there is a clear, positive correlation between the average Hα and [OII] luminosities for [OII] emitters at z=1.47, with its slope being consistent with the local relation. [OII] emitters at z=1.47 have lower mean observed ratios of Hα/[OII] suggesting a small but systematic offset (at 2.8σ significance) towards lower values of dust attenuation, A(Hα)~0.35, than local galaxies. This confirms that [OII] selection tends to pick up galaxies which are significantly less dusty on average than Hα selected ones, with the difference being higher at z=1.47 than at z=0. The discrepancy of the observed line ratios between [OII] emitters at z=1.47 and the local galaxies may in part be due to the samples having different metallicities. However, we demonstrate that metallicity is unlikely to be the main cause. Therefore, it is important to take into account that the relations for the dust correction which are derived using Hα emitter samples, and frequently used in many studies of high-z galaxies, may overestimate the intrinsic SFRs of [OII]-selected galaxies, and that surveys of [OII] emission galaxies are likely to miss dusty populations.

preprint2012arXiv

The Properties of the Star-Forming Interstellar Medium at z=0.8-2.2 from HiZELS - II: Star-Formation and Clump Scaling Laws in Gas Rich, Turbulent Disks

We present adaptive optics assisted integral field spectroscopy of nine Halpha-selected galaxies at z=0.84--2.23 selected from the HiZELS narrow-band survey. Our observations map the star-formation and kinematics of these representative star-forming galaxies on ~kpc-scales. We demonstrate that within the ISM of these galaxies, the velocity dispersion of the star-forming gas (σ) follows a scaling relation σ\proptoΣ_SFR^(1/n)+constant (where Σ_SFR is the star formation surface density and the constant includes the stellar surface density). Assuming the disks are marginally stable (Toomre Q=1), we show that this follows from the Kennicutt-Schmidt relation (Σ_SFR = AΣ_gas^n), and we use the data to derive best fit parameters of n=1.34+/-0.15 and A=3.4_(-1.6)^(+2.5)x10^(-4)Mo/yr/kpc^2, consistent with the local relation and implying cold molecular gas masses of M_gas=10^(9-10)Mo and molecular gas fractions M_gas/(M_gas+Mstars)=0.3+/-0.1, with a range of 10-75%. These values confirm the high gas fractions for high-redshift star-forming galaxies, independent of CO-H_2 conversion factor. We also identify eleven ~kpc-scale star-forming regions (clumps) within our sample and show that their sizes are comparable to the wavelength of the fastest growing unstable mode. The luminosities and velocity dispersions of these clumps follow the same scaling relations as local HII regions, although their star formation densities are a factor 15+/-5x higher than typically found locally. We discuss how the clump properties are related to the disk, and show that their high masses and luminosities are a consequence of the high disk surface density.

preprint2012arXiv

The Properties of the Star-Forming Interstellar Medium at z=0.84-2.23 from HiZELS - I: Mapping the Internal Dynamics and Metallicity Gradients in High-Redshift Disk Galaxies

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine H-alpha-selected galaxies at z=0.84--2.23 drawn from the HiZELS narrow-band survey. These galaxies have star-formation rates of 1-27Mo/yr and are therefore representative of the typical high-redshift star-forming population. Our ~kpc-scale resolution observations show that approximately half of the sample have dynamics suggesting that the ionised gas is in large, rotating disks. We model their velocity fields to infer the inclination-corrected, asymptotic rotational velocities. We use the absolute B-band magnitudes and stellar masses to investigate the evolution of the B-band and stellar mass Tully-Fisher relationships. By combining our sample with a number of similar measurements from the literature, we show that, at fixed circular velocity, the stellar mass of star-forming galaxies has increased by a factor 2.5 between z=2 and z=0, whilst the rest-frame B-band luminosity has decreased by a factor ~6 over the same period. Together, these demonstrate a change in mass-to-light ratio in the B-band of Delta(M/L_B)/(M/L_B)_(z=0) \sim 3.5 between z=1.5 and z=0, with most of the evolution occurring below z=1. We also use the spatial variation of [NII]/Halpha to show that the metallicity of the ionised gas in these galaxies declines monotonically with galacto-centric radius, with an average Delta(log O/H)/DeltaR=-0.027+/-0.005dex/kpc. This gradient is consistent with predictions for high-redshift disk galaxies from cosmologically based hydrodynamic simulations.

preprint2011arXiv

Star formation at z=1.47 from HiZELS: An Hα+[OII] double-blind study

This paper presents the results from the first wide and deep dual narrow-band survey to select H-alpha (Ha) and [OII] line emitters at z=1.47+-0.02 (using matched narrow-band filters in the H and z' bands), exploiting synergies between the UKIRT and Subaru telescopes. The Ha survey at z=1.47 reaches a flux limit of ~7x10^-17 erg/s/cm^2 and detects ~200 Ha emitters over 0.7deg^2, while the much deeper [OII] survey reaches an effective flux of ~7x10^-18 erg/s/cm^2, detecting ~1400 z=1.47 [OII] emitters in a matched co-moving volume of ~2.5x10^5 Mpc^3. The combined survey results in the identification of 190 simultaneous Ha and [OII] emitters at z=1.47. Ha and [OII] luminosity functions are derived and both are shown to evolve significantly from z~0 in a consistent way. The star formation rate density of the Universe at z=1.47 is evaluated, with the Ha analysis yielding 0.16+-0.05 M_sun/yr/Mpc^3 and the [OII] analysis 0.17+-0.04 M_sun/yr/Mpc^3. The measurements are combined with other studies, providing a self-consistent measurement of the star formation history of the Universe over the last ~11Gyrs. By using a large comparison sample at z~0.1 (from the SDSS), [OII]/Ha line ratios are calibrated as probes of dust-extinction. Ha emitters at z~1.47 show on average 1 mag of extinction at Ha, similar to the SDSS sources at z~0. Although we find that dust extinction correlates with SFR, the relation evolves by about ~0.5 mag from z~1.5 to z~0, with z~0 relations over-predicting the dust extinction corrections at high-z by that amount. Stellar mass is found to be a much more fundamental extinction predictor, with the relation between mass and extinction being valid at both z~0 and z~1.5. Dust extinction corrections as a function of optical colours are also derived, offering simpler mechanisms for estimating extinction in moderately star-forming systems over the last ~9Gyrs [Abridged].

preprint2010arXiv

From fields to a super-cluster: the role of the environment at z=0.84 with HiZELS

At z=0, clusters are primarily populated by red, elliptical and massive galaxies, while blue, spiral and lower-mass galaxies are common in low-density environments. Understanding how and when these differences were established is of absolute importance for our understanding of galaxy formation and evolution, but results at high-z remain contradictory. By taking advantage of the widest and deepest H-alpha narrow-band survey at z=0.84 over the COSMOS and UKIDSS UDS fields, probing a wide range of densities (from poor fields to rich groups and clusters, including a confirmed super-cluster with a striking filamentary structure), we show that the fraction of star-forming galaxies falls continuously from ~40% in fields to approaching 0% in rich groups/clusters. We also find that the median SFR increases with environmental density, at least up to group densities - but only for low and medium mass galaxies, and thus such enhancement is mass-dependent at z~1. The environment also plays a role in setting the faint-end slope (alpha) of the H-alpha luminosity function. Our findings provide a sharper view on galaxy formation and evolution and reconcile previously contradictory results at z~1: stellar mass is the primary predictor of star formation activity, but the environment also plays a major role.

preprint2010arXiv

HiZELS: the High Redshift Emission Line Survey with UKIRT

In these proceedings we report on HiZELS, the High-z Emission Line Survey, our successful panoramic narrow-band Campaign Survey using WFCAM on UKIRT to detect and study emission line galaxies at z~1-9. HiZELS employs the H2(S1) narrow-band filter together with custom-made narrow-band filters in the J and H-bands, with the primary aim of delivering large, identically-selected samples of H-alpha emitting galaxies at redshifts of 0.84, 1.47 and 2.23. Comparisons between the luminosity function, the host galaxy properties, the clustering, and the variation with environment of these H-alpha-selected samples are yielding unique constraints on the nature and evolution of star-forming galaxies, across the peak epoch of star-formation activity in the Universe. We provide a summary of the project status, and detail the main scientific results obtained so far: the measurement of the evolution of the cosmic star-formation rate density out to z > 2 using a single star-formation indicator, determination of the morphologies, environments and dust-content of the star-forming galaxies, and a detailed investigation of the evolution of their clustering properties. We also summarise the on-going work and future goals of the project.

preprint2010arXiv

The clustering and evolution of H-alpha emitters at z~1 from HiZELS

(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.

preprint2009arXiv

Obscured star formation at z = 0.84 with HiZELS: the relationship between star formation rate and H-alpha or ultra-violet dust extinction

[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.

preprint2009arXiv

The HiZELS/UKIRT large area survey for bright Lyman-alpha emitters at z~9

We present the largest area survey to date (1.4 deg2) for Lyman-alpha emitters (LAEs) at z~9, as part of the Hi-z Emission Line Survey (HiZELS). The survey, which primarily targets H-alpha emitters at z < 3, uses the Wide Field CAMera on the United Kingdom Infrared Telescope and a custom narrow-band filter in the J band to reach a Lyman-alpha luminosity limit of ~10^43.8 erg/s over a co-moving volume of 1.12x10^6 Mpc^3 at z = 8.96+-0.06. Two candidates were found out of 1517 line emitters, but those were rejected as LAEs after follow-up observations. This improves the limit on the space density of bright Lyman-alpha emitters by 3 orders of magnitude and is consistent with suppression of the bright end of the Lyman-alpha luminosity function beyond z~6. Combined with upper limits from smaller but deeper surveys, this rules out some of the most extreme models for high-redshift Lyman-alpha emitters. The potential contamination of narrow-band Lyman-alpha surveys at z>7 by Galactic brown dwarf stars is also examined, leading to the conclusion that such contamination may well be significant for searches at 7.7 < z < 8.0, 9.1 < z < 9.5 and 11.7 < z < 12.2.

David Sobral

What is connected

Connect this record

See the researcher in context

Building this map preview

67 published item(s)

J-HERTz: J-PLUS Heritage Exploration of Radio Targets at z $<$ 5

(Re)Solving Reionization with Lyα: How Bright Lyα Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

Black hole virial masses from single-epoch photometry: the miniJPAS test case

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

The LEGA-C and SAMI Galaxy Surveys: Quiescent Stellar Populations and the Mass-Size Plane across 6 Gyr

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

The miniJPAS survey: Identification and characterization of the emission line galaxies down to $z < 0.35$ in the AEGIS field

J-PLUS: Support Vector Machine Applied to STAR-GALAXY-QSOClassification

The LEGA-C of nature and nurture in stellar populations of galaxies at z~0.6-1.0: D4000 and H-delta reveal different assembly histories for quiescent galaxies in different environments

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Dust, gas, and metal content in star-forming galaxies at $z\sim3.3$ revealed with ALMA and Near-IR spectroscopy

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

Spectroscopic Confirmation of a Coma Cluster Progenitor at z ~ 2.2

Stellar Kinematics and Environment at z~0.8 in the LEGA-C Survey: Massive, Slow-Rotators are Built First in Overdense Environments

The evolution of rest-frame UV properties, Lya EWs and the SFR-Stellar mass relation at z~2-6 for SC4K LAEs

The nature of CR7 revealed with MUSE: a young starburst powering extended Lyman-$α$ emission at z=6.6

The X-ray and radio activity of typical and luminous Lya emitters from z~2 to z~6: evidence for a diverse, evolving population

Chandra Observations of the Spectacular A3411-12 Merger Event

Resolved Lyman-α properties of a luminous Lyman-break galaxy in a large ionised bubble at z = 6.53

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

Unveiling the most luminous Lyman-α emitters in the epoch of reionisation

VIS3COS: III. environmental effects on the star formation histories of galaxies at z~0.8 seen in [OII], H$δ$, and Dn4000

Effects of Local Environment and Stellar Mass on Galaxy Quenching out to z~3

Lyman Alpha Signatures from Direct Collapse Black Holes

The CALYMHA survey: Ly$α$ escape fraction and its dependence on galaxy properties at $z=2.23$

The growth of typical star-forming galaxies and their super massive black holes across cosmic time since z~2

The KMOS Redshift One Spectroscopic Survey (KROSS): Dynamical properties, gas and dark matter fractions of typical z~1 star-forming galaxies

The KMOS Redshift One Spectroscopic Survey (KROSS): The Tully-Fisher Relation at z ~ 1

The Lyman-alpha luminosity function at z=5.7-6.6 and the steep drop of the faint end: implications for reionization

The most luminous H$α$ emitters at z~0.8-2.23 from HiZELS: evolution of AGN and star-forming galaxies

The nature of H-alpha star-forming galaxies at z~0.4 in and around Cl 0939+4713: the environment matters

The Nature of H$β+$[O{\sc iii}] and [O{\sc ii}] emitters to $z \sim 5$ with HiZELS: stellar mass functions and the evolution of EWs

The production and escape of Lyman-Continuum radiation from star-forming galaxies at z~2 and their redshift evolution

A Comparative Study of Density Field Estimation for Galaxies: New Insights into the Evolution of Galaxies with Environment in COSMOS out to z~3

A large narrow band H$α$ survey at $z\sim0.2$: the bright end of the luminosity function, cosmic variance and clustering across cosmic time

CF-HiZELS, a 10 deg$^2$ emission-line survey with spectroscopic follow-up: Hα, [OIII]+Hβ and [OII] luminosity functions at z=0.8, 1.4 and 2.2

Evidence for PopIII-like stellar populations in the most luminous Lyman-$α$ emitters at the epoch of re-ionisation: spectroscopic confirmation

Evolution of the H$β$+[OIII] and [OII] luminosity functions and the [OII] star-formation history of the Universe up to $z$ ~ 5 from HiZELS

Identification of the brightest Lyα emitters at z=6.6: implications for the evolution of the luminosity function in the re-ionisation era

MC$^2$: Boosted AGN and star-formation activity in CIZA J2242.8+5301, a massive post-merger cluster at z=0.19

Nebular and Stellar Dust Extinction Across the Disk of Emission-Line Galaxies on Small (KPC) Scales

Neutral hydrogen gas, past and future star-formation in galaxies in and around the 'Sausage' merging galaxy cluster

Spectroscopic Study of Star-forming Galaxies in Filaments and the Field at $z\sim$0.5: Evidence for Environmental Dependence of Electron Density

A 10 deg$^2$ Lyman-$α$ survey at z=8.8 with spectroscopic follow-up: strong constraints on the LF and implications for other surveys

A relationship between specific star formation rate and metallicity gradient within z=1 galaxies from KMOS-HiZELS

Cosmic Web and Star Formation Activity in Galaxies at z~1

Mapping the large scale structure around a z=1.46 galaxy cluster in 3-D using two adjacent narrow-band filters

MC^2: Constraining the Dark Matter Distribution of the Violent Merging Galaxy Cluster CIZA J2242.8+5301 by Piercing through the Milky Way

MC^2: Galaxy Imaging and Redshift Analysis of the Merging Cluster CIZA J2242.8+5301

Spectral and morphological signatures of cluster merger shocks: CIZA J2242.8+5301

The highest-frequency detection of a radio relic: 16-GHz AMI observations of the `Sausage' cluster

The rise and fall of star-formation in $\bf z\sim0.2$ merging galaxy clusters

A fundamental metallicity relation for galaxies at z = 0.84 - 1.47 from HiZELS

On the evolution and environmental dependence of the star formation rate versus stellar mass relation since z~2

The merger rates and sizes of galaxies across the peak epoch of star formation from the HiZELS survey

The role of cluster mergers and travelling shocks in shaping the H$α$ luminosity function at $\bf z\sim0.2$: `sausage' and `toothbrush' clusters

The stellar mass function of star-forming galaxies and the mass-dependent SFR function since z=2.23 from HiZELS

A large Hα survey at z=2.23, 1.47, 0.84 & 0.40: the 11 Gyr evolution of star-forming galaxies from HiZELS

Calibrating [OII] star-formation rates at z>1 from dual Hα-[OII] imaging from HiZELS

The Properties of the Star-Forming Interstellar Medium at z=0.8-2.2 from HiZELS - II: Star-Formation and Clump Scaling Laws in Gas Rich, Turbulent Disks

The Properties of the Star-Forming Interstellar Medium at z=0.84-2.23 from HiZELS - I: Mapping the Internal Dynamics and Metallicity Gradients in High-Redshift Disk Galaxies

Star formation at z=1.47 from HiZELS: An Hα+[OII] double-blind study

From fields to a super-cluster: the role of the environment at z=0.84 with HiZELS

HiZELS: the High Redshift Emission Line Survey with UKIRT

The clustering and evolution of H-alpha emitters at z~1 from HiZELS

Obscured star formation at z = 0.84 with HiZELS: the relationship between star formation rate and H-alpha or ultra-violet dust extinction

The HiZELS/UKIRT large area survey for bright Lyman-alpha emitters at z~9