Lyman Break Galaxies Research Articles

ALPACA: a new semi-analytical model for metal absorption lines emerging from clumpy galactic environments

ABSTRACT We present a new semi-analytical formalism for modelling metal absorption lines that emerge from a clumpy galactic environment, ALPACA. We predict the “down-the-barrel” (DTB) metal absorption line profiles and the equivalent width (EW) of absorption at different impact parameters (b) as a function of the clump properties, including clump kinematics, clump volume filling factor, clump number density profile, and clump ion column densities. With ALPACA, we jointly model the stacked DTB C ii λ1334 spectrum of a sample of z ∼ 3 Lyman break galaxies and the EW versus b profile of a sample of z ∼ 2 star-forming galaxy–galaxy pairs. ALPACA successfully reproduced two data sets simultaneously, and the best fit prefers a low clump volume filling factor (∼3 × 10−3). The radial velocities of the clumps are a superposition of a rapidly accelerated outflow with a maximum velocity of $\sim 400 \, {\mathrm{km}\, \mathrm{s}^{-1}}$ and a velocity dispersion of $\sigma \sim 120 \, {\mathrm{km}\, \mathrm{s}^{-1}}$. The joint modelling reveals a physical scenario where the absorption observed at a particular velocity is contributed by the clumps distributed over a fairly broad range of radii. We also find that the commonly adopted Sobolev approximation is at best only applicable within a narrow range of radii where the clumps are undergoing rapid acceleration in a non-volume-filling clumpy medium. Lastly, we find that the clump radial velocity profile may not be fully constrained by the joint modelling and spatially resolved Ly α emission modelling may help break the degeneracy.

Galaxy clustering measurements out to redshift z ˜ 8 from Hubble Legacy Fields

ABSTRACT We present a novel approach for measuring the two-point correlation function of galaxies in narrow pencil beam surveys with varying depths. Our methodology is utilized to expand high-redshift galaxy clustering investigations up to z ∼ 8 by analysing a comprehensive sample consisting of Ng = 160 Lyman break galaxy candidates obtained through optical and near-infrared photometric data within the CANDELS GOODS data sets from the Hubble Space Telescope Legacy Fields. For bright sources with MUV < −19.8, we determine a galaxy bias of b = 9.33 ± 4.90 at $\overline{z} = 7.7$ and a correlation length of r0 = 10.74 ± 7.06 $h^{-1}\, \mathrm{Mpc}$. We obtain similar results for the XDF, with a galaxy bias measurement of b = 8.26 ± 3.41 at the same redshift for a slightly fainter sample with a median luminosity of MUV = −18.4. By comparing with dark-matter halo bias and employing abundance matching, we deduce a characteristic halo mass of Mh ∼ 1011.5 M⊙ and a duty cycle close to unity. To validate our approach for variable-depth data sets, we replicate the analysis in a region with near-uniform depth using a standard two-point correlation function estimator, yielding consistent outcomes. Our study not only provides a valuable tool for future utilization in JWST data sets but also suggests that the clustering of early galaxies continues to increase with redshift beyond z ≳ 8, potentially contributing to the existence of protocluster structures observed in early JWST imaging and spectroscopic surveys at z ≳ 8.

Lyman-α at Cosmic Noon II: The relationship between kinematics and Lyman-α in z ∼ 2–3 Lyman break galaxies

Abstract We report for the first time a relationship between galaxy kinematics and net Lyman- $\alpha$ equivalent width (net Ly $\alpha$ EW) in star-forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Ly $\alpha$ -emitting and Ly $\alpha$ -absorbing Lyman break galaxies (LBGs) at $z\sim2$ (Paper I in this series) and previously at $z\sim3$ , we use the Ly $\alpha$ spectral type classification method to study the relationship between net Ly $\alpha$ EW and nebular emission-line kinematics in samples of $z\sim2$ and $z\sim3$ LBGs drawn from the literature for which matching rest-frame UV photometry, consistently measured net Ly $\alpha$ EWs, and kinematic classifications from integral field unit spectroscopy are available. We show that $z\sim2$ and $z\sim3$ LBGs segregate in colour-magnitude space according to their kinematic properties and Lyman- $\alpha$ spectral type and conclude that LBGs with Ly $\alpha$ dominant in absorption (aLBGs) are almost exclusively rotation-dominated (presumably disc-like) systems, and LBGs with Ly $\alpha$ dominant in emission (eLBGs) characteristically have dispersion-dominated kinematics. We quantify the relationship between the strength of rotational dynamic support (as measured using ${v}_{\mathrm{obs}}/2{\sigma }_{\mathrm{int}}$ and ${v}_{\mathrm{rot}}/{\sigma}_{\mathrm{0}}$ ) and net Ly $\alpha$ EW for subsets of our kinematic sample where these data are available, and demonstrate the consistency of our result with other properties that scale with net Ly $\alpha$ EW and kinematics. Based on these findings, we suggest a method by which large samples of rotation- and dispersion-dominated galaxies might be selected using broadband imaging in as few as three filters and/or net Ly $\alpha$ EW alone. If confirmed with larger samples, application of this method will enable an understanding of galaxy kinematic behaviour over large scales in datasets from current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs in redshift ranges from $z\sim2-6$ across many hundreds to thousands of Mpc. Finally, we speculate that the combination of our result linking net Ly $\alpha$ EW and nebular emission-line kinematics with the known large-scale clustering behaviour of Ly $\alpha$ -absorbing and Ly $\alpha$ -emitting LBGs is evocative of an emergent bimodality of early galaxies that is consistent with a nascent morphology-density relation at $z\sim2-3$ .

The ALMA REBELS survey: obscured star formation in massive Lyman-break galaxies at z = 4–8 revealed by the IRX–β and M⋆ relations

ABSTRACT We investigate the degree of dust obscured star formation in 49 massive (log10(M⋆/M⊙) > 9) Lyman-break galaxies (LBGs) at z = 6.5–8 observed as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame ultraviolet (UV), optical, and far-infrared (FIR) properties which reveal a significant fraction (fobs = 0.4–0.7) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer (β ≃ −2.2) colours than expected from an extrapolation of the colour–magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\, {\rm K}$ and βd = 2.0), we find that the REBELS sources are in agreement with the local ‘Calzetti-like’ starburst Infrared-excess (IRX)–β relation. By re-analysing the data available for 108 galaxies at z ≃ 4–6 from the ALMA Large Program to Investigate C+ at Early Times (ALPINE) using a consistent methodology and assumed FIR SED, we show that from z ≃ 4–8, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX–β relation. When comparing the IRX–M⋆ relation derived from the combined ALPINE and REBELS sample to relations established at z < 4, we find a deficit in the IRX, indicating that at z > 4 the proportion of obscured star formation is lower by a factor of ≳ 3 at a given a M⋆. Our IRX–β results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for z ≃ 7 galaxies with similar stellar masses and star formation rates.

Hidden Giants in JWST's PEARLS: An Ultramassive z = 4.26 Submillimeter Galaxy that Is Invisible to HST

We present a multiwavelength analysis using the Submillimeter Array (SMA), James Clerk Maxwell Telescope, NOEMA, JWST, the Hubble Space Telescope (HST), and the Spitzer Space Telescope of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A 1489. These SMA-located sources both lie at z = 4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at ≲ 2 μm, even with deep JWST/NIRCam observations. We investigate their stellar, interstellar medium, and dynamical properties, including a pixel-level spectral energy distribution analysis to derive subkiloparsec-resolution stellar-mass and A V maps. We find that 850.1 is one of the most massive and highly obscured, A V ∼ 5, galaxies known at z > 4 with M * ∼1011.8 M ⊙ (likely forming at z > 6), and 850.2 is one of the least massive and least obscured, A V ∼ 1, members of the z > 4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z ≳ 3–4 based on imaging at ≲ 2 μm, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar-mass distribution, ∼1 kpc, but its expected evolution means that it matches both the properties of massive, quiescent galaxies at z ∼ 1.5 and ultramassive early-type galaxies at z ∼ 0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature that we link to the active dynamical environment it resides within.

Spatially resolved Kennicutt–Schmidt relation at z ≈ 7 and its connection with the interstellar medium properties

ABSTRACT We exploit moderately resolved [O $\scriptstyle \rm III$], [C $\scriptstyle \rm II$] and dust continuum ALMA observations to derive the gas density (n), the gas-phase metallicity (Z), and the deviation from the Kennicutt–Schmidt (KS) relation (κs) on $\approx \, \rm sub-kpc$ scales in the interstellar medium (ISM) of five bright Lyman Break Galaxies at the Epoch of Reionization (z ≈ 7). To do so, we use GLAM, a state-of-art, physically motivated Bayesian model that links the [C $\scriptstyle \rm II$]and [O $\scriptstyle \rm III$]surface brightness (Σ[CII], Σ[OIII]) and the SFR surface density (ΣSFR) to n, κs, and Z. All five sources are characterized by a central starbursting region, where the Σgas versus ΣSFR align ≈10 × above the KS relation (κs ≈ 10). This translates into gas depletion times in the range tdep ≈ 80 − 250 Myr. The inner starbursting centres are characterized by higher gas density (log (n/cm−3) ≈ 2.5–3.0) and higher metallicity (log (Z/Z⊙) ≈ −0.5) than the galaxy outskirts. We derive marginally negative radial metallicity gradients (∇log Z ≈ −0.03 ± 0.07 dex/kpc), and a dust temperature (Td ≈ 32 − 38 K) that anticorrelates with the gas depletion time.

Evolution of the Lyman-α-emitting fraction and UV properties of lensed star-forming galaxies in the range 2.9 < z < 6.7

Context. Faint galaxies are theorised to have played a major role, perhaps the dominant role, in reionising the Universe. Their properties, as well as the Lyman-α emitter (LAE) fraction, XLAE, could provide useful insights into this epoch. Aims. We used four clusters of galaxies from the Lensed Lyman-alpha MUSE Arcs Sample (LLAMAS) that also have deep HST photometry to select a population of intrinsically faint Lyman break galaxies (LBGs) and LAEs. We study the interrelation between these two populations, their properties, and the fraction of LBGs that display Lyman-α emission. Methods. The use of lensing clusters allows us to access an intrinsically faint population of galaxies, the largest such sample collected for this purpose: 263 LAEs and 972 LBGs with redshifts between 2.9 and 6.7, Lyman-α luminosities in the range 39.5 ≲ log(LLyα)(erg s−1)≲42, and absolute UV magnitudes in the range −22 ≲ M1500 ≲ −12. In addition to matching LAEs and LBGs, we define an LAE+continuum sample for the LAEs that match with a continuum object that is not selected as an LBG. Additionally, with the use of MUSE integral field spectroscopy, we detect a population of LAEs completely undetected in the continuum. Results. We find a redshift evolution of XLAE in line with literature results, with diminished values above z = 6. In line with past studies, we take this as signifying an increasingly neutral intervening intergalactic medium. When inspecting this redshift evolution with different limits on EWLyα and M1500, we find that the XLAE for the UV-brighter half of our sample is higher than the XLAE for the UV-fainter half, a difference that increases at higher redshifts. This is a surprising result and can be interpreted as the presence of a population of low Lyman-α equivalent width (EWLyα), UV-bright galaxies situated in reionised bubbles and overdensities. This result is especially interesting in the context of similar, UV-bright, low EWLyα objects recently detected during and around the epoch of reionisation. For intrinsically fainter objects, we confirm the previously observed trend of LAEs among LBGs as galaxies with high star formation rates and low dust content, as well as the trend of the strongest LAEs having, in general, fainter M1500 and steeper UV slopes.

Bright submillimeter galaxies do trace galaxy protoclusters

There is controversy in the literature regarding whether distant, massive, and dusty starbursts selected at (sub)millimeter wavelengths can trace galaxy overdensities. We thus performed the first systematic search for distant protoclusters around a homogeneously selected sample of 12 spectroscopically confirmed submillimeter galaxies (SMGs) at z ∼ 1.2 − 5.3, which we selected from the GOODS-N field. We applied the well-established Poisson probability method (PPM) to search for megaparsec-scale overdensities around these SMGs, using three different photometric redshift catalogs. We robustly detect galaxy overdensities for 11 out of the 12 SMGs (i.e., 92%±8%), distributed over eight large-scale protoclusters. We confirm all three previously discovered protoclusters, and we detect five new ones around the SMGs SMM J123634 (z = 1.225), ID.19 (z = 2.047), SMM J123607 (z = 2.487), SMM J123606 (z = 2.505), and GN10 (z = 5.303). A wavelet-based analysis of the protocluster fields shows that the SMGs are located in protocluster cores with a complex morphology (compact, filamentary, or clumpy) and an average size of ∼(0.4 − 1) Mpc. By comparing the PPM results obtained using the three redshift catalogs independently, each of which trace different galaxy populations and redshift ranges, we speculate that we are possibly witnessing a transitioning phase at z ≳ 4 for the galaxy population of protoclusters. While z ≲ 4 protoclusters appear to be populated by dusty galaxies, those at the highest redshifts, z ∼ 5, are detected as overdensities of Lyman α emitters or Lyman break galaxies. Further investigation with larger samples is required to reach a definitive conclusion. We also find a good correlation between the molecular (H2) gas mass of the SMGs and the significance of the associated overdensity. To explain the overall phenomenology, we suggest that galaxy interactions in dense environments likely triggered the starburst and gas-rich phase of the SMGs. Altogether, our findings support the scenario that SMGs are excellent tracers of distant protoclusters. The ones presented in this work are excellent targets for the James Webb Space Telescope. Similarly, future surveys with forthcoming facilities (e.g., Euclid and LSST) can be tuned to detect even larger samples of distant protoclusters.

An overdensity of Lyman break galaxies around the hot dust-obscured galaxy WISE J224607.56−052634.9

We report the identification of Lyman break galaxy (LBG) candidates around the most luminous hot dust-Obscured galaxy (Hot DOG) known, WISE J224607.56−052634.9 (W2246−0526) at z = 4.601, using deep r-, i-, and z-band imaging from the Gemini Multi-Object Spectrograph South (GMOS-S). We used the surface density of LBGs to probe the megaparsec-scale environment of W2246−0526 to characterize its richness and evolutionary state. We identified LBG candidates in the vicinity of W2246−0526 using the selection criteria developed in the Subaru Deep Field and in the Subaru XMM-Newton Deep Field, slightly modified to account for the difference between the filters used, and we find 37 and 55 LBG candidates, respectively. Matching to the z-band depths of those studies, this corresponds to δ = 5.8−1.9+2.4 times the surface density of LBGs expected in the field. Interestingly, the Hot DOG itself, as well as a confirmed neighbor, do not satisfy either LBG selection criteria, suggesting we may be missing a large number of companion galaxies. Our analysis shows that we have mostly likely only found those with a higher than average intergalactic medium (IGM) optical depth or moderately high dust obscuration. The number density of LBG candidates is not concentrated around W2246−0526, suggesting an early evolutionary stage for the proto-cluster, that the Hot DOG may not be the most massive galaxy, or that the Hot DOG may be affecting the IGM transparency in its vicinity. The overdensity around W2246−0526 is comparable to overdensities found around other Hot DOGs and is somewhat higher than what is typically found for radio galaxies and luminous quasars at a similar redshift.

SCUBA-2 High Redshift Bright Quasar Survey. II. The Environment of z ∼ 6 Quasars at Submillimeter Band

The formation of the first supermassive black holes is expected to have occurred in some most pronounced matter and galaxy overdensities in the early universe. We have conducted a submillimeter wavelength continuum survey of 54 z ∼ 6 quasars using the Submillimeter Common-User Bolometre Array-2 on the James Clerk Maxwell Telescope to study the environments around z ∼ 6 quasars. We identified 170 submillimeter galaxies (SMGs) with above 3.5σ detections in 450 or 850 μm maps. Their far-IR luminosities are (2.2–6.4) × 1012 L ⊙, and their star formation rates are ∼400–1200 M ⊙ yr−1. We also calculated the SMGs’ differential and cumulative number counts in a combined area of ∼620 arcmin2. To a 4σ detection (at ∼5.5 mJy), SMGs’ overdensity is (±0.19), exceeding the blank-field source counts by a factor of 1.68. We find that 13/54 quasars show overdensities (at ∼5.5 mJy) of δ SMG ∼ 1.5–5.4. The combined area of these 13 quasars exceeds the blank-field counts with the overdensity to 5.5 mJy of δ SMG ∼ (±0.25) in the regions of ∼150 arcmin2. However, the excess is insignificant on the bright end (e.g., 7.5 mJy). We also compare results with previous environmental studies of Lyα emitters and Lyman break galaxies on a similar scale. Our survey presents the first systematic study of the environment of quasars at z ∼ 6. The newly discovered SMGs provide essential candidates for follow-up spectroscopic observations to test whether they reside in the same large-scale structures as the quasars and search for protoclusters at an early epoch.

JADES NIRSpec Spectroscopy of GN-z11: Lyman-α emission and possible enhanced nitrogen abundance in a z = 10.60 luminous galaxy

We present JADES JWST/NIRSpec spectroscopy of GN-z11, the most luminous candidate z > 10 Lyman break galaxy in the GOODS-North field with MUV = −21.5. We derive a redshift of z = 10.603 (lower than previous determinations) based on multiple emission lines in our low and medium resolution spectra over 0.7 − 5.3 μm. We significantly detect the continuum and measure a blue rest-UV spectral slope of β = −2.4. Remarkably, we see spatially extended Lyman-α in emission (despite the highly neutral intergalactic medium expected at this early epoch), offset 555 km s−1 redwards of the systemic redshift. From our measurements of collisionally excited lines of both low and high ionisation (including [O II] λ3727, [Ne III] λ3869, and C III] λ1909), we infer a high ionisation parameter (log U ∼ −2). We detect the rarely seen N IV] λ1486 and N III] λ1748 lines in both our low and medium resolution spectra, with other high ionisation lines seen in the low resolution spectrum, such as He II (blended with O III]) and C IV (with a possible P-Cygni profile). Based on the observed rest-UV line ratios, we cannot conclusively rule out photoionisation from an active galactic nucleus (AGN), although the high C III]/He II and N III]/He II ratios are compatible with a star formation explanation. If the observed emission lines are powered by star formation, then the strong N III] λ1748 observed may imply an unusually high N/O abundance. Balmer emission lines (Hγ, Hδ) are also detected, and if powered by star formation rather than an AGN, we infer a star formation rate of ∼20 − 30 M⊙ yr−1 (depending on the initial mass function) and low dust attenuation. Our NIRSpec spectroscopy confirms that GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr after the Big Bang.

ALMA as a Redshift Machine: Using [C ii] to Efficiently Confirm Galaxies in the Epoch of Reionization

The [C ii]158μm line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper, we present the results of new ALMA observations spectral scanning for [C ii] in six particularly luminous Lyman-break galaxies at z ∼ 7. The six sources were drawn from a sample of bright z ∼ 7 galaxies identified using the wide-area optical, near-IR, and Spitzer/IRAC data over the COSMOS/UltraVISTA field and were targeted on the basis of tight constraints on their redshifts from their IRAC [3.6]–[4.5] colors. We detect significant (>9σ) [C ii] lines in three of our six targets (50%) cospatial with the rest-UV emission from the ground/space-based near-IR imaging. The luminosities of the [C ii] lines lie in the range 5.6–8.8 × 108 L ⊙, consistent with the local [C ii]–SFR relation. Meanwhile, their [C ii]/L IR ∼ 1–3 × 10−3 ratios are slightly elevated compared to local (U)LIRGS. This could be due to lower dust-to-gas or dust-to-metal ratios. We also find that our sources display a large kinematic diversity, with one source showing signs of rotation, one source a likely major merger, and one dispersion-dominated source that might contain a bright star-forming clump. Our results highlight the effectiveness of spectral scans with ALMA in spectroscopically confirming luminous galaxies in the epoch of reionization, something that is being be applied on a significantly larger sample in the ongoing REBELS large program.

The JWST FRESCO survey: legacy NIRCam/grism spectroscopy and imaging in the two GOODS fields

ABSTRACT We present the JWST cycle 1 53.8 h medium program FRESCO, short for ‘First Reionization Epoch Spectroscopically Complete Observations’. FRESCO covers 62 arcmin2 in each of the two GOODS/CANDELS fields for a total area of 124 arcmin2 exploiting JWST’s powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ∼2 h deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R ∼ 1600 covering 3.8–5.0 µm for most galaxies in the NIRCam field of view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z ∼ 0.2–0.5, to Pa α and Pa β at z ∼ 1–3, He i and [S iii] at z ∼ 2.5–4.5, H α and [N ii] at z ∼ 5–6.5, up to [O iii] and H β for z ∼ 7–9 galaxies. FRESCO’s grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ∼28.2 mag (5σ in 0${_{.}^{\prime\prime}}$32 diameter apertures). Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 yr ago, the discovery of broad-line active galactic nuclei at z > 4, and resolved Pa α maps of galaxies at z ∼ 1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations.

A Mixture of LBG Overdensities in the Fields of Three 6 < z < 7 Quasars: Implications for the Robustness of Photometric Selection

The most luminous quasars at z > 6 are suspected to be both highly clustered and to reside in the most massive dark matter halos in the early universe, making them prime targets to search for galaxy overdensities and/or protoclusters. We search for Lyman-break dropout-selected galaxies using HST WFC3/ACS broadband imaging in the fields of three 6 < z < 7 quasars, as well as for their simultaneously observed coordinated parallel fields, and constrain their photometric redshifts using EAZY. One field, J0305-3150, shows a volume density 10× higher than the blank-field UV luminosity function (UVLF) at M UV < −20, with tentative evidence of a 3σ overdensity in its parallel field located 15 cMpc away. Another field, J2054–0005, shows an angular overdensity within 500 ckpc from the quasar, but this is still consistent with UVLF predictions within 3σ, while the last field, J2348–3054, shows no enhancement. We discuss methods for reducing uncertainty in overdensity measurements when using photometric selection and show that we can robustly select Lyman-break galaxies (LBGs) consistent with being physically associated with the quasar, corroborated by existing JWST/NIRCam WFSS data in the J0305 field. Even accounting for incompleteness, the overdensities in J0305 and J2054 are higher for brighter galaxies at short angular separations, suggesting preferential enhancement of more massive galaxies in the quasar’s immediate vicinity. Finally, we compare the population of LBGs with previously identified [C ii] and millimeter-continuum companions. The LBG overdensities are not accompanied by an enhanced number of dusty galaxies, suggesting that the overdense quasar fields are not in the bursty star-forming phase sometimes seen in high-redshift protoclusters.

The bright end of the galaxy luminosity function at z ≃ 7 from the VISTA VIDEO survey

ABSTRACT We have conducted a search for z ≃ 7 Lyman-break galaxies over 8.2 deg2 of near-infrared imaging from the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey in the XMM–Newton-Large Scale Structure (XMM-LSS) and the Extended Chandra Deep Field-South (ECDF-S) fields. Candidate galaxies were selected from a full photometric redshift analysis down to a Y + J depth of 25.3 (5σ), utilizing deep auxiliary optical and Spitzer/Infrared Array Camera (IRAC) data to remove brown dwarf and red interloper galaxy contaminants. Our final sample consists of 28 candidate galaxies at 6.5 ≤ z ≤ 7.5 with −23.5 ≤ MUV ≤ −21.6. We derive stellar masses of 9.1 ≤ log10(M⋆/M⊙) ≤ 10.9 for the sample, suggesting that these candidates represent some of the most massive galaxies known at this epoch. We measure the rest-frame ultraviolet (UV) luminosity function (LF) at z ≃ 7, confirming previous findings of a gradual decline in number density at the bright end (MUV &amp;lt; −22) that is well described by a double power law (DPL). We show that quasar contamination in this magnitude range is expected to be minimal, in contrast to conclusions from recent pure-parallel Hubble studies. Our results are up to a factor of 10 lower than previous determinations from optical-only ground-based studies at MUV ≲ −23. We find that the inclusion of YJHKs photometry is vital for removing brown dwarf contaminants, and z ≃ 7 samples based on red optical data alone could be highly contaminated (≳50 per cent). In comparison with other robust z &amp;gt; 5 samples, our results further support little evolution in the very bright end of the rest-frame UV LF from z = 5–10, potentially signalling a lack of mass quenching and/or dust obscuration in the most massive galaxies in the first Gyr.

Big Three Dragons: Molecular Gas in a Bright Lyman-break Galaxy at z = 7.15

We report Atacama Large Millimeter/submillimeter Array Band 3 observations of CO(6−5), CO(7−6), and [C i](2−1) in B14-65666 (“Big Three Dragons”), one of the brightest Lyman-break galaxies at z > 7 in the rest-frame ultraviolet continuum, far-infrared continuum, and emission lines of [O iii] 88 μm and [C ii] 158 μm. CO(6−5), CO(7−6), and [C i](2−1), whose 3σ upper limits on the luminosities are approximately 40 times fainter than the [C ii] luminosity, are all not detected. The L [C II]/L CO(6–5) and L [C II]/L CO(7–6) ratios are higher than the typical ratios obtained in dusty star-forming galaxies or quasar host galaxies at similar redshifts, and they may suggest a lower gas density in the photodissociated region in B14-65666. By using the (1) [C ii] luminosity, (2) dust mass-to-gas mass ratio, and (3) a dynamical mass estimate, we find that the molecular gas mass (M mol) is (0.05–11) × 1010 M ⊙. This value is consistent with the upper limit inferred from the nondetection of mid-J CO and [C i](2−1). Despite the large uncertainty in M mol, we estimate a molecular gas-to-stellar mass ratio (μ gas) of 0.65–140 and a gas depletion time (τ dep) of 2.5–550 Myr; these values are broadly consistent with those of other high-redshift galaxies. B14-65666 could be an ancestor of a passive galaxy at z ≳ 4 if no gas is fueled from outside the galaxy.

On the Metallicities and Kinematics of the Circumgalactic Media of Damped Lyα Systems at z ∼ 2.5* * The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by

We use medium- and high-resolution spectroscopy of close pairs of quasars to analyze the circumgalactic medium (CGM) surrounding 32 damped Lyα absorption systems (DLAs). The primary quasar sightline in each pair probes an intervening DLA in the redshift range 1.6 < z abs < 3.5, such that the secondary sightline probes absorption from Lyα and a large suite of metal-line transitions (including O i, C ii, C iv, Si ii, and Si iv) in the DLA host galaxy’s CGM at transverse distances 24 kpc ≤ R ⊥ ≤ 284 kpc. Analysis of Lyα in the CGM sightlines shows an anticorrelation between R ⊥ and H i column density (N HI) with 99.8% confidence, similar to that observed around luminous galaxies. The incidences of C ii and Si ii with N > 1013 cm−2 within 100 kpc of DLAs are larger by 2σ than those measured in the CGM of Lyman break galaxies (C f (N CII ) > 0.89 and ). Metallicity constraints derived from ionic ratios for nine CGM systems with negligible ionization corrections and N HI > 1018.5 cm−2 show a significant degree of scatter (with metallicities/limits across the range ), suggesting inhomogeneity in the metal distribution in these environments. Velocity widths of C iv λ1548 and low-ionization metal species in the DLA versus CGM sightlines are strongly (>2σ) correlated, suggesting that they trace the potential well of the host halo over R ⊥ ≲ 300 kpc scales. At the same time, velocity centroids for C iv λ1548 differ in DLA versus CGM sightlines by >100 km s−1 for ∼50% of velocity components, but few components have velocities that would exceed the escape velocity assuming dark matter host halos of ≥1012 M ⊙.

Broad Emission Lines in Optical Spectra of Hot, Dust-obscured Galaxies Can Contribute Significantly to JWST/NIRCam Photometry

Selecting the first galaxies at z > 7 − 10 from JWST surveys is complicated by z < 6 contaminants with degenerate photometry. For example, strong optical nebular emission lines at z < 6 may mimic JWST/NIRCam photometry of z > 7–10 Lyman-break galaxies (LBGs). Dust-obscured 3 < z < 6 galaxies in particular are potentially important contaminants, and their faint rest-optical spectra have been historically difficult to observe. A lack of optical emission line and continuum measures for 3 < z < 6 dusty galaxies now makes it difficult to test their expected JWST/NIRCam photometry for degenerate solutions with NIRCam dropouts. Toward this end, we quantify the contribution by strong emission lines to NIRCam photometry in a physically motivated manner by stacking 21 Keck II/NIRES spectra of hot, dust-obscured, massive () and infrared (IR) luminous galaxies at z ∼ 1–4. We derive an average spectrum and measure strong narrow (broad) [O iii]5007 and Hα features with equivalent widths of 130 ± 20 Å (150 ± 50 Å) and 220 ± 30 Å (540 ± 80 Å), respectively. These features can increase broadband NIRCam fluxes by factors of 1.2 − 1.7 (0.2–0.6 mag). Due to significant dust attenuation (A V ∼ 6), we find Hα+[N ii] to be significantly brighter than [O iii]+Hβ and therefore find that emission-line dominated contaminants of high −z galaxy searches can only reproduce moderately blue perceived UV continua of S λ ∝ λ β with β > − 1.5 and z > 4. While there are some redshifts (z ∼ 3.75) where our stack is more degenerate with the photometry of z > 10 LBGs at λ rest ∼ 0.3–0.8 μm , redder filter coverage beyond λ obs > 3.5 μm and far-IR/submillimeter follow-up may be useful for breaking the degeneracy and making a crucial separation between two fairly unconstrained populations, dust-obscured galaxies at z ∼ 3–6 and LBGs at z > 10.

The Hα Luminosity Function of Galaxies at z ∼ 4.5

We present the Hα luminosity function (LF) derived from a large sample of Lyman break galaxies at z ∼ 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The Hα flux is derived from the offset between the continuum flux estimated from the best-fit spectral energy distribution, and the observed photometry in IRAC [3.6]. From these measurements, we build the Hα LF and study its evolution providing the best constraints of this property at high redshift, where spectroscopy of Hα is not yet available. Schechter parameterizations of the Hα LF show a decreasing evolution of Φ∗ with redshift, increasing evolution in L ∗, and no significant evolution in the faint-end slope at high z. We find that star formation rates (SFRs) derived from Hα are higher than those derived from the rest-frame UV for low SFR galaxies but the opposite happens for the highest SFRs. This can be explained by lower mass galaxies (also lower SFR) having, on average, rising star formation histories (SFHs), while at the highest masses the SFHs may be declining. The SFR function is steeper, and because of the excess SFR(Hα) compared to SFR(UV) at low SFRs, the SFR density estimated from Hα is higher than the previous estimates based on UV luminosities.

The hidden side of cosmic star formation at z &gt; 3

Our current understanding of the cosmic star formation history at z &gt; 3 is primarily based on UV-selected galaxies (Lyman-break galaxies, i.e., LBGs). Recent studies of H-dropouts (HST-dark galaxies) have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z &gt; 3. In this work, we extend the H-dropout criterion to lower masses to select optically dark or faint galaxies (OFGs) at high redshifts in order to complete the census between LBGs and H-dropouts. Our criterion (H &gt; 26.5 mag &amp; [4.5] &lt; 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies (typically E(B − V) &gt; 0.4) with lower stellar masses at high redshifts. In addition, with this criterion, our sample is not contaminated by massive passive or old galaxies. In total, we identified 27 OFGs at zphot &gt; 3 (with a median of zmed = 4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log(M⋆/M⊙) = 9.4 − 11.1 (with a median of log(M⋆med/M⊙) = 10.3). We find that up to 75% of the OFGs with log(M⋆/M⊙) = 9.5 − 10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing an optical-to-millimeter stacking analysis of the OFGs, we find that rather than being limited to a rare population of extreme starbursts, these OFGs represent a normal population of dusty star-forming galaxies at z &gt; 3. The OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than the scaling relation of typical star-forming galaxies. Additionally, the total star formation rate (SFRtot = SFRIR + SFRUV) of the stacked OFGs is much higher than the SFRUVcorr (SFRUV corrected for dust extinction), with an average SFRtot/SFRUVcorr = 8 ± 1, which lies above (∼0.3 dex) the 16–84th percentile range of typical star-forming galaxies at 3 ≤ z ≤ 6. All of the above suggests the presence of hidden dust regions in the OFGs that absorb all UV photons, which cannot be reproduced with dust extinction corrections. The effective radius of the average dust size measured by a circular Gaussian model fit in the uv plane is Re(1.13 mm) = 1.01 ± 0.05 kpc. After excluding the five LBGs in the OFG sample, we investigated their contributions to the cosmic star formation rate density (SFRD). We found that the SFRD at z &gt; 3 contributed by massive OFGs (log(M⋆/M⊙) &gt; 10.3) is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculated the combined contribution of OFGs and LBGs to the cosmic SFRD at z = 4 − 5 to be 4 × 10−2 M⊙ yr−1 Mpc−3, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift. This value could be even larger, as our calculations were performed in a very conservative way.