Search

Dwarf galaxies hold the key to crucial frontiers of astrophysics, however, their faintness renders spectroscopy challenging. Here we present the JWST Cycle 2 survey, All the Little Things (ALT, PID 3516), which is designed to seek late-forming Pop III stars and the drivers of reionization at $z\sim6-7$. ALT has acquired the deepest NIRCam grism spectroscopy yet (7-27 hr), at JWST's most sensitive wavelengths (3-4 $\mu$m), covering the powerful lensing cluster Abell 2744. Over the same 30 arcmin$^2$, ALT's ultra-deep F070W+F090W imaging ($\sim$30 mag) enables selection of very faint sources at $z>6$. We demonstrate the success of ALT's novel ``butterfly" mosaic to solve spectral confusion and contamination, and introduce the ``Allegro" method for emission line identification. By collecting spectra for every source in the field of view, ALT has measured precise ($R\sim1600$) redshifts for 1630 sources at $z=0.2-8.5$. This includes one of the largest samples of distant dwarf galaxies: [1015, 475, 50] sources less massive than the SMC, Fornax, and Sculptor with $\log(M_{*}/M_{\odot})<$[8.5, 7.5, 6.5]. We showcase ALT's discovery space with: (i) spatially resolved spectra of lensed clumps in galaxies as faint as $M_{\rm{UV}}\sim-15$; (ii) large-scale clustering -- overdensities at $z$=[2.50, 2.58, 3.97, 4.30, 5.66, 5.77, 6.33] hosting massive galaxies with striking Balmer breaks; (iii) small-scale clustering -- a system of satellites around a Milky Way analog at $z\sim6$; (iv) spectroscopically confirmed multiple images that help constrain the lensing model underlying all science in this legacy field; (v) sensitive star-formation maps based on dust-insensitive tracers such as Pa$\alpha$; (vi) direct spectroscopic discovery of rare sources such as AGN with ionized outflows. These results provide a powerful proof of concept for how grism surveys maximize the potential of strong lensing fields., Comment: Submitted to the Open Journal of Astrophysics. Figs. 6 and 10 illustrate the quality of the spectra and imaging, while Fig. 12 summarizes the yield of the survey. Comments warmly welcomed and greatly appreciated. The ALT DR1 catalog is available at https://zenodo.org/records/13871850

We investigate the `intrinsic' H$\alpha$ EW distributions of $z \sim 0.4 - 2.2$ narrowband-selected H$\alpha$ samples from HiZELS and DAWN using a forward modeling approach. We find an EW - stellar mass anti-correlation with steepening slopes $-0.18\pm0.03$ to $-0.24^{+0.06}_{-0.08}$ at $z \sim 0.4$ and $z\sim 2.2$, respectively. Typical EW increases as $(1+z)^{1.78^{+0.22}_{-0.23}}$ for a $10^{10}$ M$_\odot$ emitter from $15^{+2.4}_{-2.3}$\r{A} ($z \sim 0.4$) to $67.7^{+10.4}_{-10.0}$\r{A} ($z \sim 2.2$) and is steeper with decreasing stellar mass highlighting the high EW nature of low-mass high-$z$ systems. We model this redshift evolving anti-correlation, $W_0(M,z)$, and find it produces H$\alpha$ luminosity and SFR functions strongly consistent with observations validating the model and allowing us to use $W_0(M,z)$ to investigate the relative contribution of H$\alpha$ emitters towards cosmic SF. We find EW$_0 > 200$ \r{A} emitters contribute significantly to cosmic SF activity at $z \sim 1.5 - 2$ making up $\sim 40$% of total SF consistent with sSFR $> 10^{-8.5}$ yr$^{-1}$ ($\sim 45 - 55$%). Overall, this highlights the importance of high EW systems at high-$z$. Our $W_0(M,z)$ model also reproduces the cosmic sSFR evolution found in simulations and observations and show that tension between the two can simply arise from selection effects in observations. Lastly, we forecast Roman and Euclid grism surveys using $W_0(M,z)$ including observational efficiency and limiting resolution effects where we predict $\sim 24000$ and $\sim 30000$ $0.5 < z < 1.9$ H$\alpha$ emitters per deg$^{-2}$, respectively, down to $>5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ including $10^{7.2 - 8}$ M$_\odot$ galaxies at $z > 1$ with EW$_0 >1000$\r{A}. Both Roman and Euclid will enable us to observe with unprecedented detail some of the most bursty/high EW, low-mass star-forming galaxies near cosmic noon., Comment: 25 pages, 12 Figures, and 12 Tables. Submitted to MNRAS. Comments are welcomed. (Abridged Abstract)

Submillimetre galaxies (SMGs) are some of the most extreme star-forming systems in the Universe, whose place in the framework of galaxy evolution is as yet uncertain. It has been hypothesised that SMGs are progenitors of local early-type galaxies, requiring that SMGs generally reside in galaxy cluster progenitors at high redshift. We test this hypothesis and explore SMG environments using a narrowband VLT/HAWK-I+GRAAL study of H$\alpha$ and [OIII] emitters around an unbiased sample of three ALMA-identified and spectroscopically-confirmed SMGs at $z \sim 2.3$ and $z \sim 3.3$, where these SMGs were selected solely on spectroscopic redshift. Comparing with blank-field observations at similar epochs, we find that one of the three SMGs lies in an overdensity of emission-line sources on the $\sim4$ Mpc scale of the HAWK-I field of view, with overdensity parameter $\delta_{g} = 2.6^{+1.4}_{-1.2}$. A second SMG is significantly overdense only on $\lesssim 1.6$ Mpc scales and the final SMG is consistent with residing in a blank field environment. The total masses of the two overdensities are estimated to be $\log(M_{h}/{\rm M}_{\odot}) =$12.1--14.4, leading to present-day masses of $\log(M_{h,z=0}/{\rm M}_{\odot}) =$12.9--15.9. These results imply that SMGs occupy a range of environments, from overdense protoclusters or protogroups to the blank field, suggesting that while some SMGs are strong candidates for the progenitors of massive elliptical galaxies in clusters, this may not be their only possible evolutionary pathway., Comment: 21 pages, 12 figures. Accepted for publication in MNRAS

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$\alpha$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$\alpha$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35M_{\rm UV}>-19.5$)., Comment: 15 pages, 13 figures. Accepted for publication in A&A. Main figs. 1, 2, 10 (SFR surface density), 11 (spectra and positions of galaxies in the bubble). Corrected error in escape fraction estimate from SFR surface density

We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for both quiescent and star-forming galaxies within this redshift range, though we note the presence of weak evolution in the number densities of star-forming galaxies. We compare both VDFs with previous direct and inferred measurements at local and intermediate redshifts, with the caveat that previous measurements of the VDF for star-forming galaxies are poorly constrained at all epochs. We emphasize that this work is the first to directly push to low-stellar velocity dispersion ($\sigma_\star > 100$ km s$^{-1}$) and extend to star-forming galaxies. We are largely consistent with the high-sigma tail measured from BOSS, and we find that the VDF remains constant from the median redshift of LEGA-C, $z\sim0.8$, to the present day., Comment: 11 pages, 3 figures, accepted for publication in ApJ

We use multi-object near-infrared (NIR) spectroscopy with VLT/KMOS to investigate the role of the environment in the evolution of the ionized gas properties of narrow-band selected H$\alpha$ emitters (HAEs) in the Spiderweb protocluster at $z=2.16$. Based on rest-frame optical emission lines, H$\alpha$ and [NII]$\lambda$6584, we confirm the cluster membership of 39 of our targets (i.e. 93% success rate), and measure their star-formation rates (SFR), gas-phase oxygen abundances and effective radius. We parametrize the environment where our targets reside by using local and global density indicators based on previous samples of spectroscopic and narrow-band cluster members. We find that star-forming galaxies embedded in the Spiderweb protocluster display SFRs compatible with those of the main sequence and morphologies comparable to those of late-type galaxies at $z=2.2$ in the field. We also report a mild gas-phase metallicity enhancement ($0.6\pm0.3$ dex) at intermediate stellar masses. Furthermore, we identify two UVJ-selected quiescent galaxies with residual H$\alpha$-based star formation and find signs of extreme dust obscuration in a small sample of SMGs based on their FIR and H$\alpha$ emission. Interestingly, the spatial distribution of these objects differs from the rest of HAEs, avoiding the protocluster core. Finally, we explore the gas fraction-gas metallicity diagram for 7 galaxies with molecular gas masses measured by ATCA using CO(1-0). In the context of the gas-regulator model, our objects are consistent with relatively low mass-loading factors, suggesting lower outflow activity than field samples at the cosmic noon and thus, hinting at the onset of environmental effects in this massive protocluster., Comment: 29 pages, 15 figures. Accepted for publication in MNRAS

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 \sigma'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $\sigma'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a S\'ersic 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 $\sigma'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $\sigma'_{\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\'ersic profile-based $R$ estimate, then a S\'{e}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., Comment: Published in ApJ

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\alpha$, $H\beta$, [OIII], and [NII] emission lines. We explore the ionization mechanism using the diagrams [OIII]/H$\beta$ versus [NII]/H$\alpha$ (BPT) and EW(H$\alpha$) versus [NII]/H$\alpha$ (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\alpha$ 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 ($\rho_{\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\alpha$ emission line., Comment: 22 pages, 19 figures

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., Comment: 19 pages, 7 figures, accepted for publication by ApJ, in press

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

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., Comment: 16 pages, 10 figures, accepted by A&A

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$\alpha$ 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$\alpha$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Ly\alpha}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$\alpha$ 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$\alpha$ 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\%$, $\xi_{\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$\alpha$-anchored formalism., Comment: Accepted for publication in MNRAS. Figure 2 shows the main result -- the comoving emissivity due to bright LAEs. Our fiducial model is based on results presented in our companion paper -- Naidu & Matthee et al. 2022, arXiv: 2110.11961

The ionizing photon escape fraction (LyC $f_{\rm{esc}}$) of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshift. Here we circumvent this uncertainty by inferring $f_{\rm{esc}}$ with resolved (R>4000) LyA profiles from the X-SHOOTER LyA survey at z=2 (XLS-z2). We select leakers ($f_{\rm{esc}}>20$%) and non-leakers ($f_{\rm{esc}}<5$%) from a representative sample of $>0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 {\AA} (rest-frame) we investigate the conditions for LyC $f_{\rm{esc}}$. We find the following differences between leakers vs. non-leakers: (i) strong nebular CIV and HeII emission vs. non-detections, (ii) O32~8.5 vs. ~3, (iii) Ha/Hb indicating no dust vs. E(B-V)~0.3, (iv) MgII emission close to the systemic velocity vs. redshifted, optically thick MgII, (v) LyA $f_{\rm{esc}}$ of ~50% vs. ~10%. The extreme EWs in leakers (O3+Hb~1100 {\AA}) constrain the characteristic timescale of LyC escape to ~3-10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers -- extremely ionizing stellar populations, low column densities, a dust-free, high ionization state ISM -- occur simultaneously in the $f_{\rm{esc}}>20\%$ stack, suggesting they are causally connected, and motivating why indicators like O32 may suffice to constrain $f_{\rm{esc}}$ at z>6 with JWST. The leakers comprise half our sample, have a median LyC $f_{\rm{esc}}$~50%, and an ionising production efficiency $\log({\xi_{\rm{ion}}/\rm{Hz\ erg^{-1}}})$~25.9. These results show LAEs -- the type of galaxies rare at z=2, but that become the norm at higher redshift -- are highly efficient ionizers, with extreme $\xi_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED), Comment: Submitted to MNRAS. Figure 4 and Table 2 show key results from the stacks. Figure 8 presents a unifying scheme for the LyC duty cycle. Our companion paper (Matthee & Naidu et al. 2021) works out implications of these results for the cosmic UV background at $z\approx2-8$. Comments warmly welcomed and greatly appreciated!

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., Comment: ApJ accepted: 16 pages, 11 figures

We present the first comparison of observed stellar continuum spectra of high-redshift galaxies and mock galaxy spectra generated from hydrodynamical simulations. The mock spectra are produced from the IllustrisTNG TNG100 simulation combined with stellar population models and take into account dust attenuation and realistic observational effects (aperture effects and noise). We compare the simulated $D_n4000$ and EW(H$\delta$) of galaxies with $10.5 \leq \log(M_\ast/M_\odot) \leq 11.5$ at $0.6 \leq z \leq 1.0$ to the observed distributions from the LEGA-C survey. TNG100 globally reproduces the observed distributions of spectral indices, implying that the age distribution of galaxies in TNG100 is generally realistic. Yet there are small but significant differences. For old galaxies, TNG100 shows small $D_n4000$ when compared to LEGA-C, while LEGA-C galaxies have larger EW(H$\delta$) at fixed $D_n4000$. There are several possible explanations: 1) LEGA-C galaxies have overall older ages combined with small contributions (a few percent in mass) from younger ($<1$~Gyr) stars, while TNG100 galaxies may not have such young sub-populations; 2) the spectral mismatch could be due to systematic uncertainties in the stellar population models used to convert stellar ages and metallicities to observables. In conclusion, the latest cosmological galaxy formation simulations broadly reproduce the global age distribution of galaxies at $z\sim1$ and, at the same time, the high quality of the latest observed and simulated datasets help constrain stellar population synthesis models as well as the physical models underlying the simulations., Comment: 17 pages, 12 figures, accepted by AJ

We present spatially resolved stellar kinematics for 797 $z=0.6-1$ galaxies selected from the LEGA-C survey and construct axisymmetric Jeans models to quantify their dynamical mass and degree of rotational support. The survey is $K_s$-band selected, irrespective of color or morphological type, and allows for a first assessment of the stellar dynamical structure of the general $L^*$ galaxy population at large lookback time. Using light profiles from Hubble Space Telescope imaging as a tracer, our approach corrects for observational effects (seeing convolution and slit geometry), and uses well-informed priors on inclination, anisotropy and a non-luminous mass component. Tabulated data include total mass estimates in a series of spherical apertures (1, 5, and 10 kpc; 1$\times$ and 2$\times$\re), as well as rotational velocities, velocity dispersions and anisotropy. We show that almost all star-forming galaxies and $\sim$50\% of quiescent galaxies are rotation-dominated, with deprojected $V/\sigma\sim1-2$. Revealing the complexity in galaxy evolution, we find that the most massive star-forming galaxies are among the most rotation-dominated, and the most massive quiescent galaxies among the least rotation-dominated galaxies. These measurements set a new benchmark for studying galaxy evolution, using stellar dynamical structure for galaxies at large lookback time. Together with the additional information on stellar population properties from the LEGA-C spectra, the dynamical mass and $V/\sigma$ measurements presented here create new avenues for studying galaxy evolution at large lookback time., Comment: Accepted for publication in ApJ. Data table will be published with journal article and is now available upon request

We present the third and final data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO/VLT public spectroscopic survey targeting $0.6 < z < 1.0$, Ks-selected galaxies. The data release contains 3528 spectra with measured stellar velocity dispersions and stellar population properties, a 25-fold increase in sample size compared to previous work. This $K_s$-selected sample probes the galaxy population down to $\sim0.3 L^*$, for all colors and morphological types. Along with the spectra we publish a value-added catalog with stellar and ionized gas velocity dispersions, stellar absorption line indices, emission line fluxes and equivalent widths, complemented with structural parameters measured from HST/ACS imaging. With its combination of high precision and large sample size, LEGA-C provides a new benchmark for galaxy evolution studies., Comment: Accepted for publication in ApJS. Data available at http://archive.eso.org/cms/eso-archive-news/Third-and-final-release-of-the-Large-Early-Galaxy-Census-LEGA-C-Spectroscopic-Public-Survey-published.html or https://users.ugent.be/\string~avdrwel/research.html\#legac

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 $\Sigma$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$\Sigma$ 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--$\Sigma$ 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., Comment: 18 pages, 11 figures, accepted to MNRAS

We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-alpha (Lya) emitters (LAEs) from z~2 to z~6 by exploring ~4000 LAEs from the SC4K sample. We find a correlation between Lya luminosity (LLya) and rest-frame UV (M_UV), with best-fit M_UV=-1.6+-0.2 log10(LLya/erg/s)+47+-12 and a shallower relation between LLya and stellar mass (Mstar), with best-fit log10( Mstar/Msun)=0.9+-0.1 log10(LLya/erg/s)-28+-4.0. An increasing LLya cut predominantly lowers the number density of faint M_UV and low Mstar LAEs. We estimate a proxy for the full UV LFs and SMFs of LAEs with simple assumptions of the faint end slope. For the UV LF, we find a brightening of the characteristic UV luminosity (M_UV*) with increasing redshift and a decrease of the characteristic number density (Phi*). For the SMF, we measure a characteristic stellar mass (Mstar*/Msun) increase with increasing redshift, and a Phi* decline. However, if we apply a uniform luminosity cut of log10 (LLya/erg/s) >= 43.0, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (rho_UV) of the full sample of LAEs shows moderate evolution and the stellar mass density (rho_M) decreases, with both being always lower than the total rho_UV and rho_M of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both rho_UV and rho_M of LAEs slowly approach the measurements of continuum-selected galaxies at z>6, which suggests a key role of LAEs in the epoch of reionisation., Comment: Accepted for publication in MNRAS. 20 pages, 13 figures, 4 tables

We explore the connection between the kinematics, structures and stellar populations of massive galaxies at $0.610.5$) galaxies that span a wide range in morphology, star formation activity and environment, and therefore is representative of the massive galaxy population at $z\sim0.8$. We find that quiescent and star-forming galaxies occupy the parameter space of the $g$-band FP differently and thus have different distributions in the dynamical mass-to-light ratio ($M_{\rm dyn}/L_g$), largely owing to differences in the stellar age and recent star formation history, and, to a lesser extent, the effects of dust attenuation. In contrast, we show that both star-forming and quiescent galaxies lie on the same mass FP at $z\sim 0.8$, with a comparable level of intrinsic scatter about the plane. We examine the variation in $M_{\rm dyn}/M_*$ through the thickness of the mass FP, finding no significant residual correlations with stellar population properties, S\'ersic index, or galaxy overdensity. Our results suggest that, at fixed size and velocity dispersion, the variations in $M_{\rm dyn}/L_g$ of massive galaxies reflect an approximately equal contribution of variations in $M_*/L_g$, and variations in the dark matter fraction or initial mass function., Comment: 32 pages, 19 figures (including appendices); accepted for publication in ApJ

The growth of galaxy clusters is energetic and may trigger and/or quench star formation and black hole activity. The ENISALA project is a collection of multiwavelength observations aimed at understanding how large-scale structure drives galaxy and black hole evolution. Here, we introduce optical spectroscopy of over 800 H$\alpha$ emission-line galaxies, selected in 14 z~0.15-0.31 galaxy clusters, spanning a range of masses and dynamical states. We investigate the nature of the emission lines in relation to the host galaxy properties, its location within the cluster, and the properties of the parent cluster. We uncover remarkable differences between mergers and relaxed clusters. The majority of H$\alpha$ emission-line galaxies in merging cluster fields are located within 3 Mpc of their center. A large fraction of these line-emitters in merging clusters are powered by star formation irrespective of cluster-centric radius, while the rest are powered by active galactic nuclei. Star-forming galaxies are rare within 3 Mpc of relaxed clusters and active galactic nuclei are most abundant at their outskirts (~1.5-3 Mpc). We discover a population of star-forming galaxies with large equivalent widths and blue UV-optical colors, found exclusively in the merging clusters in our sample. The widespread emission-line activity in merging clusters is likely supported by triggered activity in recently-accreted, gas-rich galaxies. By contrast, our observations for relaxed clusters match established models, in which black hole activity is enhanced at the virial radius and star-formation is quenched within the infall region. We conclude that emission-line galaxies experience distinct evolutionary paths in merging and relaxed clusters., Comment: Accepted for publication in the Astrophysical Journal. 30 pages, 22 figures, 3 tables. Main results can be found in Figures 14, 15, 16, 17 and 19

We present the first results from the X-SHOOTER Lyman-$\alpha$ survey at $z=2$ (XLS-$z2$). XLS-$z2$ is a deep spectroscopic survey of 35 Lyman-$\alpha$ emitters (LAEs) utilising $\approx90$ hours of exposure time with VLT/X-SHOOTER and covers rest-frame Ly$\alpha$ to H$\alpha$ emission with R$\approx4000$. We present the sample selection, the observations and the data reduction. Systemic redshifts are measured from rest-frame optical lines for 33/35 sources. In the stacked spectrum, our LAEs are characterised by an interstellar medium with little dust, a low metallicity and a high ionisation state. The ionising sources are young hot stars that power strong emission-lines in the optical and high ionisation lines in the UV. The LAEs exhibit clumpy UV morphologies and have outflowing kinematics with blue-shifted SiII absorption, a broad [OIII] component and a red-skewed Ly$\alpha$ line. Typically 30 % of the Ly$\alpha$ photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly$\alpha$ photons escapes directly at the systemic suggesting clear channels enabling a $\approx10$ % escape of ionising photons, consistent with an inference based on MgII. A combination of a low effective HI column density, a low dust content and young star-burst determine whether a star forming galaxy is observed as a LAE. The first is possibly related to outflows and/or a fortunate viewing angle, while we find that the latter two in LAEs are typical for their stellar mass of 10$^9$ M$_{\odot}$., Comment: Main text 26 pages. Accepted for publication in MNRAS

We study the mass-metallicity relation for 19 members of a spectroscopically-confirmed protocluster in the COSMOS field at $z=2.2$ (CC2.2), and compare it with that of 24 similarly selected field galaxies at the same redshift. Both samples are $\rm H\alpha$ emitting sources, chosen from the HiZELS narrow-band survey, with metallicities derived from $\rm N2\ (\frac{\rm [NII] \lambda 6584}{\rm H \alpha})$ line ratio. For the mass-matched samples of protocluster and field galaxies, we find that protocluster galaxies with $10^{9.9} \rm M_\odot \leq M_* \leq 10^{10.9} \rm M_\odot$ are metal deficient by $0.10 \pm 0.04$ dex ($2.5\sigma$ significance) compared to their coeval field galaxies. This metal deficiency is absent for low mass galaxies, $\rm M_* < 10^{9.9} \rm M_\odot$. Moreover, relying on both SED-derived and $\rm {H\alpha}$ (corrected for dust extinction based on $\rm {M_*}$) SFRs, we find no strong environmental dependence of SFR-$\rm {M_*}$ relation, however, we are not able to rule out the existence of small dependence due to inherent uncertainties in both SFR estimators. The existence of $2.5\sigma$ significant metal deficiency for massive protocluster galaxies favors a model in which funneling of the primordial cold gas through filaments dilutes the metal content of protoclusters at high redshifts ($z \gtrsim 2$). At these redshifts, gas reservoirs in filaments are dense enough to cool down rapidly and fall into the potential well of the protocluster to lower the gas-phase metallicity of galaxies. Moreover, part of this metal deficiency could be originated from galaxy interactions which are more prevalent in dense environments., Comment: 10 pages, 5 figures, 2 tables, accepted for publication in ApJ

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., Comment: 19 pages, 7 figures, 2 tables, Accepted for publication in ApJ

Galaxy clusters grow by merging with other clusters, giving rise to Mpc-wide shock waves that travel at 1000-2500 km/s through the intra-cluster medium. To study the effects of merger shocks on the properties of cluster galaxies, we present the first spatially resolved spectroscopic view of 5 H$\alpha$ emitting galaxies located in the wake of shock fronts in the low redshift (z~0.2), massive (~2$\times10^{15}$ M$_\odot$), post-core passage merging cluster, CIZA J2242.8+5301 (nicknamed the `Sausage'). Our Gemini/GMOS-N integral field unit (IFU) observations, designed to capture H$\alpha$ and [NII] emission, reveal the nebular gas distribution, kinematics and metallicities in the galaxies over >16 kpc scales. While the galaxies show evidence for rotational support, the flux and velocity maps have complex features like tails and gas outflows aligned with the merger axis of the cluster. With gradients incompatible with inside-out disk growth, the metallicity maps are consistent with sustained star formation (SF) throughout and outside of the galactic disks. In combination with previous results, these pilot observations provide further evidence of a likely connection between cluster mergers and SF triggering in cluster galaxies, a potentially fundamental discovery revealing the interaction of galaxies with their environment., Comment: Published in ApJ Letters. Main results can be found in Figure 3. Data behind figures are available together with the journal paper

We present Atacama Large Millimeter/Submillimeter Array (ALMA) continuum observations of a sample of nine star-forming galaxies at redshifts 1.47 and 2.23 selected from the High-$z$ Emission Line Survey (HiZELS). Four galaxies in our sample are detected at high significance by ALMA at a resolution of 0.25'' at rest-frame 355 $\mu$m. Together with the previously observed H$\alpha$ emission, from adaptive optics-assisted integral-field-unit spectroscopy (0.15'' resolution), and F606W and F140W imaging from the Hubble Space Telescope (0.2'' resolution), we study the star-formation activity, stellar and dust mass in these high-redshift galaxies at $\sim$kpc-scale resolution. We find that ALMA detection rates are higher for more massive galaxies ($M_*>10^{10.5}$ M$_\odot$) and higher [N {\sc ii}]/H$\alpha$ ratios ($>0.25$, a proxy for gas-phase metallicity). The dust extends out to a radius of 8 kpc, with a smooth structure, even for those galaxies presenting clumpy H$\alpha$ morphologies. The half-light radii ($R_{\rm dust}$) derived for the detected galaxies are of the order $\sim$4.5 kpc, more than twice the size of submillimetre-selected galaxies at a similar redshift. Our global star-formation rate estimates -- from far-IR and extinction-corrected H$\alpha$ luminosities -- are in good agreement. However, the different morphologies of the different phases of the interstellar medium suggest complex extinction properties of the high-redshift normal galaxies., Comment: 13 pages, 7 figures, MNRAS accepted

We use the Fundamental Plane (FP) to measure the redshift evolution of the dynamical mass-to-light ratio ($M_{\mathrm{dyn}}/L$) and the dynamical-to-stellar mass ratio ($M_{\mathrm{dyn}}/M_*$). Although conventionally used to study the properties of early-type galaxies, we here obtain stellar kinematic measurements from the Large Early Galaxy Astrophysics Census (LEGA-C) Survey for a sample of $\sim1400$ massive ($\log( M_*/M_\odot) >10.5$) galaxies at $0.6

We present a novel approach to measure the attenuation curves of 485 individual star-forming galaxies with M$_*$ $>$ 10$^{10}$ M$_{\odot}$ based on deep optical spectra from the VLT/VIMOS LEGA-C survey and multi-band photometry in the COSMOS field. Most importantly, we find that the attenuation curves in the rest-frame $3000-4500$A range are typically almost twice as steep as the Milky Way, LMC, SMC, and Calzetti attenuation curves, which is in agreement with recent studies of the integrated light of present-day galaxies. The attenuation at $4500$A and the slope strongly correlate with the galaxy inclination: face-on galaxies show less attenuation and steeper curves compared to edge-on galaxies, suggesting that geometric effects dominate observed variations in attenuation. Our new method produces $2175$A UV bump detections for 260 individual galaxies. Even though obvious correlations between UV bump strength and global galaxy properties are absent, strong UV bumps are most often seen in face-on, lower-mass galaxies (10 $<$ log$_{10}$(M$_*$/M$_{\odot}$) $<$ 10.5) with low overall attenuation. Finally, we produce a typical attenuation curve for star-forming galaxies at $z\sim0.8$; this prescription represents the effect of dust on the integrated spectral energy distributions of high-redshift galaxies more accurately than commonly used attenuation laws., Comment: 15 pages, 10 figures, 1 table. Accepted for publication in The Astrophysical Journal

CR7 is among the most luminous Lyman-$\alpha$ emitters (LAEs) known at $z = 6.6$ and consists of at least three UV components that are surrounded by Lyman-$\alpha$ (Ly$\alpha$) 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$\alpha$ halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly$\alpha$ emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly$\alpha$ emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly$\alpha$ equivalent width (EW) is $\approx100$ {\AA}. This is a factor four higher than the EW measured in low-redshift analogues with carefully matched Ly$\alpha$ 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$\alpha$ spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly$\alpha$ blobs. CR7's Ly$\alpha$ 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., Comment: Main text 16 pages, 9 figures. Accepted for publication by MNRAS

We propose a new mission called Space Project for Astrophysical and Cosmological Exploration (SPACE) as part on the ESA long term planning Voyage 2050 programme. SPACE will study galaxy evolution at the earliest times, with the key goals of charting the formation of the heavy elements, measuring the evolution of the galaxy luminosity function, tracing the build-up of stellar mass in galaxies over cosmic time, and finding the first super-massive black holes (SMBHs) to form. The mission will exploit a unique region of the parameter space, between the narrow ultra-deep surveys with HST and JWST, and shallow wide-field surveys such as Roman Space Telescope and EUCLID, and should yield by far the largest sample of any current or planned mission of very high redshift galaxies at z > 10 which are sufficiently bright for detailed follow-up spectroscopy. Crucially, we propose a wide-field spectroscopic near-IR + mid-IR capability which will greatly enhance our understanding of the first galaxies by detecting and identifying a statistical sample of the first galaxies and the first SMBH, and to chart the metal enrichment history of galaxies in the early Universe - potentially finding signatures of the very first stars to form from metal-free primordial gas. The wide-field and wavelength range of SPACE will also provide us a unique opportunity to study star formation by performing a wide survey of the Milky Way in the near-IR + mid-IR. This science project can be enabled either by a stand-alone ESA-led M mission or by an instrument for an L mission (with ESA and/or NASA, JAXA and other international space agencies) with a wide-field (sub-)millimetre capability at wavelength > 500 microns., Comment: Submitted to Experimental Astronomy as one of the proposals for the European Space Agency Voyage 2050 Programme

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{\delta}_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $H{\delta}_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $H{\delta}_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ \r{A}/$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{A}$/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., Comment: 18 pages, 9 figures

We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.210^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.22 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s., Comment: 27 pages (6 in appendix), 21 figures (5 in appendix), submitted to A&A

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$\alpha$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $\mu$m) and $K$ ($\sim$ 1.92-2.40 $\mu$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 $\sigma_{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 $\delta_{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, $\delta_{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.

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

We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(H\delta) \geq 4$\AA) 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$., Comment: 16 pages, 8 figures, 1 table, accepted by ApJ

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., Comment: 5 pages, 3 figures. Proceeding of IAU Symposium #352, 2019, "Uncovering early galaxy evolution in the ALMA and JWST era"

[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]$\lambda$3727, H$\delta$, 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$ {\AA} to $2.0_{-0.4}^{+0.5}$ {\AA} 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., Comment: 20 pages, 13 figures, accepted for publication in A&A

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., Comment: 23 pages, 17 figures, 4 tables. Accepted for publication in ApJ

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 $\sigma$), 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., Comment: Accepted for publication in MNRAS. 23 pages, 11 figures, 3 tables + appendices. The full SC4K catalogue of LAEs with PSF photometry, stellar masses, SFRs and other properties for individual LAEs, is available at https://goo.gl/q9yfKo

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

The observed properties of the Lyman-$\alpha$ (Ly$\alpha$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$\alpha$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$\alpha$ equivalent width (EW$_0\approx38$ {\AA}). These data are combined with deep resolved [CII]$_{\rm 158 \mu m}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$\alpha$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$\alpha$ 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\alpha, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$\alpha$ 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$\alpha$ 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$\alpha$ properties in VR7 resemble the conditions in post-re-ionisation galaxies., Comment: 13 pages, 10 figures. Accepted for publication in MNRAS

We study the impact of local density and stellar mass on the morphology of ~500 quiescent (Q) and SFGs from the VIS3COS survey. We perform B/D decomposition of the SBPs and find ~41% of 10^10 Msun galaxies to be best fitted with 2 components. We complement our analysis with non-parametric measurements and qualitative visual classifications. We find that galaxy morphology depends on stellar mass and environment for our sample as a whole. We only find an impact of the environment on galaxy sizes for galaxies more massive than 10^11 Msun. We find higher n and B/T in high-density regions when compared to low-density counterparts at similar stellar masses. We also find that higher stellar mass galaxies have steeper light profiles compared to the lower ones. Using visual classifications, we find a morphology-density relation at z~0.84 for galaxies more massive than 10^10 Msun, with elliptical galaxies being dominant at high-density regions and disks more common in low-density regions. However, when splitting the sample into colour-colour selected SF and Q sub-populations, there are no statistically significant differences between low- and high-density regions. We find that Q galaxies are smaller, have higher n, and higher B/T when compared to SF counterparts. We confirm these trends with non-parametric quantities, finding Q galaxies to be smoother and having most of their light over smaller areas than SFGs. Overall, we find that the dependence of morphology on stellar mass is stronger than with local density and these relations are strongly correlated with the quenching fraction. The change in average morphology corresponds to a change in the relative fractions of blue disk-like galaxies and red elliptical galaxies with stellar mass and environment. We hypothesize that the processes responsible for the quenching of SF must also affect the galaxy morphology on similar timescales., Comment: 25 pages, 14 figures, change in the series number, to publish in Astronomy & Astrophysics

We present an analysis of the gas outflow energetics of 529 main-sequence star-forming galaxies at z~1 using KMOS observations of the broad, underlying H-alpha and forbidden lines of [N II] and [S II]. Based on the stacked spectra for a sample with median star-formation rates and stellar masses of SFR ~ 7 Mo/yr and M* = (1.0+/-0.1)x10^10 Mo respectively, we derive a typical mass outflow rate of dM/dt = 1-4 Mo/yr and a mass loading of dM/dt/SFR = 0.2--0.4. The mass loading of the wind does not show a strong trend with star-formation rate over the range SFR ~ 2--20 Mo/yr, although we identify a trend with stellar mass such that dM/dt/SFR ~ M*^(0.26+/-0.07). Finally, we find that the line width of the broad H-alpha increases with disk circular velocity with a sub-linear scaling relation FWHM_broad ~ v^(0.21+/-0.05). As a result of this behavior, in the lowest mass galaxies (M* < 10^10 Mo), a significant fraction of the outflowing gas should have sufficient velocity to escape the gravitational potential of the halo whilst in the highest mass galaxies (M* > 10^10 Mo) most of the gas will be retained, flowing back on to the galaxy disk at later times., Comment: Accepted for publication in MNRAS (14 pages, 8 figures). Data release: Datacubes and data products from KROSS are now available at http://astro.dur.ac.uk/KROSS/data.html

We use reconstructed star-formation histories (SFHs) of quiescent galaxies at $z=0.6-1$ in the LEGA-C survey to identify secondary star-formation episodes that, after an initial period of quiescence, moved the galaxies back to the star-forming main sequence (blue cloud). $16\pm3$\% of the $z\sim0.8$ quiescent population has experienced such rejuvenation events in the redshift range $0.7

We present a comparison of the observed, spatially integrated stellar and ionized gas velocity dispersions of $\sim1000$ massive ($\log M_{\star}/M_{\odot}\gtrsim\,10.3$) galaxies in the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $0.6\lesssim\,z\lesssim1.0$. The high $S/N\sim20{\rm\AA^{-1}}$ afforded by 20 hour VLT/VIMOS spectra allows for joint modeling of the stellar continuum and emission lines in all galaxies, spanning the full range of galaxy colors and morphologies. These observed integrated velocity dispersions (denoted as $\sigma'_{g, int}$ and $\sigma'_{\star, int}$) are related to the intrinsic velocity dispersions of ionized gas or stars, but also include rotational motions through beam smearing and spectral extraction. We find good average agreement between observed velocity dispersions, with $\langle\log(\sigma'_{g, int}/\sigma'_{\star, int})\rangle=-0.003$. This result does not depend strongly on stellar population, structural properties, or alignment with respect to the slit. However, in all regimes we find significant scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$, with an overall scatter of 0.13 dex of which 0.05 dex is due to observational uncertainties. For an individual galaxy, the scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$ translates to an additional uncertainty of $\sim0.24\rm{dex}$ on dynamical mass derived from $\sigma'_{g, int}$, on top of measurement errors and uncertainties from Virial constant or size estimates. We measure the $z\sim0.8$ stellar mass Faber-Jackson relation and demonstrate that emission line widths can be used to measure scaling relations. However, these relations will exhibit increased scatter and slopes that are artificially steepened by selecting on subsets of galaxies with progressively brighter emission lines., Comment: 9 pages, 5 figures, Accepted for publication in ApJ Letters

We investigate the clustering and halo properties of $\sim 5000$ Ly$\alpha$-selected emission line galaxies (LAEs) from the Slicing COSMOS 4K (SC4K) and from archival NB497 imaging of SA22 split in 15 discrete redshift slices between $z \sim 2.5 - 6$. We measure clustering lengths of $r_0 \sim 3 - 6\ h^{-1}$ Mpc and typical halo masses of $\sim 10^{11}$ M$_\odot$ for our narrowband-selected LAEs with typical $L_{\rm{Ly}\alpha} \sim 10^{42 - 43}$ erg s$^{-1}$. The intermediate band-selected LAEs are observed to have $r_0 \sim 3.5 - 15\ h^{-1}$ Mpc with typical halo masses of $\sim 10^{11 - 12}$ M$_\odot$ and typical $L_{\rm{Ly}\alpha} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$\alpha$ luminosity normalized by the characteristic Ly$\alpha$ luminosity, $L^\star(z)$. The faintest LAEs ($L \sim 0.1\ L^\star(z)$) typically identified by deep narrowband surveys are found in $10^{10}$ M$_\odot$ halos and the brightest LAEs ($L \sim 7\ L^\star(z)$) are found in $\sim 5 \times 10^{12}$ M$_\odot$ halos. A dependency on the rest-frame 1500 \AA~UV luminosity, M$_\rm{UV}$, is also observed where the halo masses increase from $10^{11}$ to $10^{13}$ M$_\odot$ for M$_\rm{UV} \sim -19$ to $-23.5$ mag. Halo mass is also observed to increase from $10^{9.8}$ to $10^{12.3}$ M$_\odot$ for dust-corrected UV star formation rates from $\sim 0.6$ to $10$ M$_\odot$ yr$^{-1}$ and continues to increase up to $10^{13.5}$ M$_\odot$ in halo mass, where the majority of those sources are AGN. All the trends we observe are found to be redshift-independent. Our results reveal that LAEs are the likely progenitors of a wide range of galaxies depending on their luminosity, from dwarf-like, to Milky Way-type, to bright cluster galaxies. LAEs therefore provide unique insight into the early formation and evolution of the galaxies we observe in the local Universe., Comment: 19 pages, 13 figures, 3 tables, submitted to MNRAS. Comments are welcomed

We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectroscopic redshifts, emission line fluxes, Lick/IDS indices, and observed stellar and gas velocity dispersions that are spatially integrated quantities including both rotational motions and genuine dispersion. To illustrate the new parameter space in the intermediate redshift regime probed by LEGA-C we explore relationships between dynamical and stellar population properties. The star-forming galaxies typically have observed stellar velocity dispersions of ~150 km/s and strong Hdelta absorption (Hd_A ~ 5 Angstrom), while passive galaxies have higher observed stellar velocity dispersions ~200 km/s and weak Hdelta absortion (Hd_A ~ 0 Angstrom). Strong [O III]5007 / Hbeta ratios tend to occur mostly for galaxies with weak Hd_A or galaxies with higher observed velocity dispersion. Beyond these broad trends, we find a large diversity of possible combinations of rest-frame colors, absorption line strengths and emission line detections, illustrating the utility of spectroscopic measurements to more accurately understand galaxy evolution. By making the spectra and value-added catalogs publicly available we encourage the community to take advantage of this very substantial investment in telescope time provided by ESO., Comment: accepted for publication by ApJ; LEGA-C data-release http://www.mpia.de/home/legac/

We analyze stellar age indicators (D$_n$4000 and EW(H$\delta$)) and sizes of 467 quiescent galaxies with $M_\ast \geq 10^{10} M_\odot$ at $z\sim0.7$ drawn from DR2 of the LEGA-C survey. Interpreting index variations in terms of equivalent single stellar population age, we find that the median stellar population is younger for larger galaxies at fixed stellar mass. The effect is significant, yet small; the ages of the larger and the smaller subsets differ by only $<500$ Myr, much less than the age variation among individual galaxies ($\sim1.5$ Gyr). At the same time, quiescent galaxies with the strongest H$\delta$ absorption --- those experienced recent and rapid quenching events --- tend to be smaller than the average. These co-existing trends unify seemingly contradictory results in the literature; the complex correlations between size and age indicators revealed by our large sample of galaxies with high-quality spectra suggest that there are multiple evolutionary pathways to quiescence. Regardless of the specific physical mechanisms responsible for the cessation of star formation in massive galaxies, the large scatter in D$_n$4000 and EW(H$\delta$) immediately implies that galaxies follow a large variety in evolutionary pathways. On the one hand, we see evidence for a process that slowly shuts off star-formation and transforms star-forming galaxies to quiescent galaxies without necessarily changing their structures. On the other hand, there is likely a mechanism that rapidly quenches galaxies, an event that coincides with dramatic structural changes, producing small post-starburst galaxies., Comment: Main part 16 pages, 13 figures. Appendix: 8 figures. Accepted by ApJ

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe., Comment: 71 pages, 42 figures, Submitted to ApJ

Distant luminous Lyman-alpha emitters (LAEs) are excellent targets for spectroscopic observations of galaxies in the epoch of reionisation (EoR). We present deep high-resolution (R=5000) VLT/X-SHOOTER observations, along with an extensive collection of photometric data of `COLA1', a proposed double peaked LAE at z=6.6 (Hu et al. 2016). We rule out that COLA1's emission line is an [OII] doublet at z=1.475 on the basis of i) the asymmetric red line-profile and flux ratio of the peaks (blue/red=$0.31\pm0.03$) and ii) an unphysical [OII]/Ha ratio ([OII]/Ha > 22). We show that COLA1's observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lya and [OIII] at z=2.142. We thus conclude that COLA1 is a real double-peaked LAE at z=6.593, the first discovered at z>6, confirming the result from Hu et al. (2016). COLA1 is UV luminous (M$_{1500}=-21.6\pm0.3$), has a high equivalent width (EW$_{0}$~120 \AA) and very compact Lya emission (r$_{50} = 0.3$ kpc). Relatively weak inferred Hb+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z<1/20 Z$_{\odot}$ or reduced strength of nebular lines due to high escape of ionising photons. The small Lya peak separation of $220\pm20$ km/s implies a low HI column density and an ionising photon escape fraction of ~15-30 %, providing the first direct evidence that such galaxies contribute actively to the reionisation of the Universe at z>6. Based on simple estimates, we find that COLA1 could have provided just enough photons to reionise its own ~0.3 pMpc bubble, allowing the blue Lya line to be observed. However, we also discuss alternative scenarios explaining the detected double peaked nature of COLA1. Our results show that future high-resolution observations of statistical samples of double peaked LAEs at z>5 are a promising probe of the occurrence of ionised regions around galaxies in the EoR., Comment: 14 pages, 10 figures. Accepted for publication by A&A. Data are fully public