ABSTRACT
We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* \gt 10^{10.3} \, \mathrm{M_\odot }$) UVJ-selected galaxies at redshifts of 1.0 < z < 1.3, combining extremely deep ...rest-frame ultraviolet spectroscopy from VANDELS with multiwavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We first constrain the stellar mass versus stellar age relationship, finding a steep trend towards earlier average formation time with increasing stellar mass (downsizing) of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* \gt 10^{11} \, \mathrm{M_\odot }$. We show that this is consistent with other spectroscopic studies from 0 < z < 2. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by simba and illustristng agree well with observations at z = 0.1, they are too shallow at z = 1, predicting an evolution of ≲0.5 Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst, and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} \, \mathrm{M_\odot }$) are consistent with formation in recent (z < 2), intense starburst events, with time-scales of ≲500 Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by z = 1, they formed at z ∼ 5 and quenched by z = 3. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity.
Context. The form and evolution of the galaxy stellar mass function (GSMF) at high redshifts provide crucial information on star formation history and mass assembly in the young Universe, close or ...even prior to the epoch of reionization. Aims. We used the unique combination of deep optical/near-infrared/mid-infrared imaging provided by HST, Spitzer, and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5 ≤ z ≤ 7.5. Methods. We used the HST WFC3/IR near-infrared imaging from CANDELS and HUDF09, reaching H ≃ 27 − 28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5 ≤ z ≤ 7.5, resulting in an overall precision of σz/ (1 + z) ~ 0.037. With this database we have determined the low-mass end of the high-redshift GSMF with unprecedented precision, reaching down to masses as low as M∗ ~ 109 M⊙ at z = 4 and ~6 × 109 M⊙ at z = 7. Results. We find that the GSMF at 3.5 ≤ z ≤ 7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photometric redshifts, the star formation histories, the nebular contribution, or the presence of AGN in the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Meanwhile the high-mass end of the GSMF appears to evolve primarily in density, although there is also some evidence of evolution in characteristic mass. Our results are very different from previous mass function estimates based on converting UV galaxy luminosity functions into mass functions via tight mass-to-light relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the star formation rate density over cosmic time at z> 4. Conclusions. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts, and reaffirm the importance of space-based infrared selection for the unbiased measurement of the evolving GSMF in the young Universe.
We unveil the complex shape of a proto-supercluster at z ∼ 2.45 in the COSMOS field exploiting the synergy of both spectroscopic and photometric redshifts. Thanks to the spectroscopic redshifts of ...the VIMOS Ultra-Deep Survey (VUDS), complemented by the zCOSMOS-Deep spectroscopic sample and high-quality photometric redshifts, we compute the three-dimensional (3D) overdensity field in a volume of ∼100 × 100 × 250 comoving Mpc3 in the central region of the COSMOS field, centred at z ∼ 2.45 along the line of sight. The method relies on a two-dimensional (2D) Voronoi tessellation in overlapping redshift slices that is converted into a 3D density field, where the galaxy distribution in each slice is constructed using a statistical treatment of both spectroscopic and photometric redshifts. In this volume, we identify a proto-supercluster, dubbed “Hyperion” for its immense size and mass, which extends over a volume of ∼60 × 60 × 150 comoving Mpc3 and has an estimated total mass of ∼4.8 × 1015 M⊙. This immensely complex structure contains at least seven density peaks within 2.4 ≲ z ≲ 2.5 connected by filaments that exceed the average density of the volume. We estimate the total mass of the individual peaks, Mtot, based on their inferred average matter density, and find a range of masses from ∼0.1 × 1014 M⊙ to ∼2.7 × 1014 M⊙. By using spectroscopic members of each peak, we obtain the velocity dispersion of the galaxies in the peaks, and then their virial mass Mvir (under the strong assumption that they are virialised). The agreement between Mvir and Mtot is surprisingly good, at less than 1−2σ, considering that (almost all) the peaks are probably not yet virialised. According to the spherical collapse model, these peaks have already started or are about to start collapsing, and they are all predicted to be virialised by redshift z ∼ 0.8−1.6. We finally perform a careful comparison with the literature, given that smaller components of this proto-supercluster had previously been identified using either heterogeneous galaxy samples (Lyα emitters, sub-mm starbursting galaxies, CO emitting galaxies) or 3D Lyα forest tomography on a smaller area. With VUDS, we obtain, for the first time across the central ∼1 deg2 of the COSMOS field, a panoramic view of this large structure, that encompasses, connects, and considerably expands in a homogeneous way on all previous detections of the various sub-components. The characteristics of this exceptional proto-supercluster, its redshift, its richness over a large volume, the clear detection of its sub-components, together with the extensive multi-wavelength imaging and spectroscopy granted by the COSMOS field, provide us the unique possibility to study a rich supercluster in formation.
Context. Determining the average fraction of Lyman continuum (LyC) photons escaping high redshift galaxies is essential for understanding how reionization proceeded in the z> 6 Universe. Aims. We ...want to measure the LyC signal from a sample of sources in the Chandra Deep Field South (CDFS) and COSMOS fields for which ultra-deep VIMOS spectroscopy as well as multi-wavelength Hubble Space Telescope (HST) imaging are available. Methods. We select a sample of 46 galaxies at z ~ 4 from the VIMOS Ultra Deep Survey (VUDS) database, such that the VUDS spectra contain the LyC part, that is, the rest-frame range 880−910 Å. Taking advantage of the HST imaging, we apply a careful cleaning procedure and reject all the sources showing nearby clumps with different colours, that could potentially be lower-redshift interlopers. After this procedure, the sample is reduced to 33 galaxies. We measure the ratio between ionizing flux (LyC at 895 Å) and non-ionizing emission (at ~ 1500 Å) for all individual sources. We also produce a normalized stacked spectrum of all sources. Results. Assuming an intrinsic average Lν(1470) /Lν(895) of 3, we estimate the individual and average relative escape fraction. We do not detect ionizing radiation from any individual source, although we identify a possible LyC emitter with very high Lyα equivalent width (EW). From the stacked spectrum and assuming a mean transmissivity for the sample, we measure a relative escape fraction f escrel = 0.09 ± 0.04. We also look for correlations between the limits in the LyC flux and source properties and find a tentative correlation between LyC flux and the EW of the Lyα emission line. Conclusions. Our results imply that the LyC flux emitted by V = 25−26 star-forming galaxies at z ~ 4 is at most very modest, in agreement with previous upper limits from studies based on broad and narrow band imaging.
Abstract
We constrain the distribution of spatially offset Lyman-alpha emission (Ly α) relative to rest-frame ultraviolet emission in ∼300 high redshift (3 < z < 5.5) Lyman-break galaxies (LBGs) ...exhibiting Ly α emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (≃0.2 deg2 total). Because slit spectroscopy only provides one spatial dimension, we use Bayesian inference to recover the underlying two-dimensional Ly α spatial offset distribution. We model the distribution using a two-dimensional circular Gaussian, defined by a single parameter σr,Ly α, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample (3 < z < 5.5), we find $\sigma _{r,\mathrm{Ly}\,\alpha }=1.70^{+0.09}_{-0.08}$ kpc ($68\hbox{ per cent}$ conf.), corresponding to ∼0${^{\prime\prime}_{.}}$25 at 〈z〉 = 4.5. We also find that σr,Ly α decreases significantly with redshift. Because Ly α spatial offsets can cause slit losses, the decrease in σr,Ly α with redshift can partially explain the increase in the fraction of Ly α emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If σr,Ly α continues to decrease into the reionization epoch, then the decrease in Ly α transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the intergalactic medium, slit losses may explain some of the drop in Ly α transmission observed at z > 6. Spatially resolved observations of Ly α and UV continuum at 6 < z < 8 are needed to settle the issue.
Context. The Frontier Fields survey is a pioneering observational program aimed at collecting photometric data, both from space (Hubble Space Telescope and Spitzer Space Telescope) and from ...ground-based facilities (VLT Hawk-I), for six deep fields pointing at clusters of galaxies and six nearby deep parallel fields, in a wide range of passbands. The analysis of these data is a natural outcome of the Astrodeep project, an EU collaboration aimed at developing methods and tools for extragalactic photometry and creating valuable public photometric catalogues. Aims. We produce multiwavelength photometric catalogues (from B to 4.5 mu m) for the first two of the Frontier Fields, Abell-2744 and MACS-J0416 (plus their parallel fields). Methods. To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members, and the intra-cluster light. We perform the detection on the processed HST H 160 image to obtain a pure H-selected sample, which is the primary catalogue that we publish. We also add a sample of sources which are undetected in the H 160 image but appear on a stacked infrared image. Photometry on the other HST bands is obtained using SExtractor, again on processed images after the procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands is obtained using our PSF-matching deconfusion code t-phot. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. Results. The procedure of foreground light subtraction allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of the complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K-band from Hawk-I, and 3.6?4.5 mu m from Spitzer. full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. Conclusions. The multiwavelength photometric catalogues are available publicly and are ready to be used for scientific purposes. Our procedures allows for the detection of outshone objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in Paper II.
ABSTRACT We present a model using both strong and weak gravitational lensing of the galaxy cluster MACS J0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space ...(GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 30 multiple images belonging to 15 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for five of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and nine multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting map of total mass density is compared with a map of stellar mass density obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, f⋆, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is (stat.), but with a systematic error as large as 0.004-0.005, dominated by the choice of the initial mass function. We find agreement with several recent measurements of f⋆ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.
The Atacama Large Millimeter Array (ALMA) Large Program to INvestigate CII at Early times (ALPINE) targets the CII 158
μ
m line and the far-infrared continuum in 118 spectroscopically confirmed ...star-forming galaxies between
z
= 4.4 and
z
= 5.9. It represents the first large CII statistical sample built in this redshift range. We present details regarding the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates (SFRs), we measured the conversion factor from the ALMA 158
μ
m rest-frame dust continuum luminosity to the total infrared luminosity (
L
IR
) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median
L
IR
of 4.4 × 10
11
L
⊙
. We also detected 57 additional continuum sources in our ALMA pointings. They are at a lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5 ± 0.2. We measured the 850
μ
m number counts between 0.35 and 3.5 mJy, thus improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40% of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the CII line in 75 of our targets. Their median CII luminosity is 4.8 × 10
8
L
⊙
and their median full width at half maximum is 252 km s
−1
. After measuring the mean obscured SFR in various CII luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR–
L
CII
relations.
ABSTRACT
We present a sample of 151 massive (M* > 1010 M⊙) quiescent galaxies at 2 < z < 5, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and ...GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at z > 3, of which 2 are at z > 4. We report formation redshifts, demonstrating that the oldest objects formed at z > 6; however, individual ages from our photometric data have significant uncertainties, typically ∼0.5 Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at z > 3, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this, we construct a model for the time evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at z ∼ 6–7 and quenching at z ∼ 5. We report spectroscopic redshifts for two of our objects at z = 3.440 and 3.396, which exhibit extremely weak Ly α emission in ultra-deep VANDELS spectra. We calculate star formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly α escape fractions are >3 and >10 per cent, respectively. We finally report that our highest redshift robust object exhibits a continuum break at λ ∼ 7000 Å in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent, this object would be the highest redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to z ≳ 3 quiescent objects.