Search

We present new observations of sixteen bright ($r=19-21$) gravitationally lensed galaxies at $z\simeq 1-3$ selected from the CASSOWARY survey. Included in our sample is the $z=1.42$ galaxy CSWA-141, one of the brightest known reionization-era analogs at high redshift (g=20.5), with a large sSFR (31.2 Gyr$^{-1}$) and an [OIII]+H$\beta$ equivalent width (EW$_{\rm{[OIII]+H\beta}}$=730~\r{A}) that is nearly identical to the average value expected at $z\simeq 7-8$. In this paper, we investigate the rest-frame UV nebular line emission in our sample with the goal of understanding the factors that regulate strong CIII] emission. Whereas most of the sources in our sample show weak UV line emission, we find elevated CIII] in the spectrum of CSWA-141 (EW$_{\rm{CIII]}}$=4.6$\pm1.9$~\r{A}) together with detections of other prominent emission lines (OIII], Si III], Fe II$^\star$, Mg II). We compare the rest-optical line properties of high redshift galaxies with strong and weak CIII] emission, and find that systems with the strongest UV line emission tend to have young stellar populations and nebular gas that is moderately metal-poor and highly ionized, consistent with trends seen at low and high redshift. The brightness of CSWA-141 enables detailed investigation of the extreme emission line galaxies which become common at $z>6$. We find that gas traced by the CIII] doublet likely probes higher densities than that traced by [OII] and [SII]. Characterisation of the spectrally resolved Mg II emission line and several low ionization absorption lines suggests neutral gas around the young stars is likely optically thin, potentially facilitating the escape of ionizing radiation., Comment: 20 pages, 9 figures, Accepted for publication in MNRAS

With JWST, new opportunities to study the evolution of galaxies in the early Universe are emerging. Spitzer constraints on rest-optical properties of $z\gtrsim7$ galaxies demonstrated the power of using galaxy stellar masses and star formation histories (SFHs) to indirectly infer the cosmic star formation history. However, only the brightest individual $z\gtrsim8$ objects could be detected with Spitzer, making it difficult to robustly constrain activity at $z\gtrsim10$. Here, we leverage the greatly improved rest-optical sensitivity of JWST at $z\gtrsim8$ to constrain the ages of seven UV-bright ($M_{UV}\lesssim-19.5$) galaxies selected to lie at $z\sim8.5-11$, then investigate implications for $z\gtrsim15$ star formation. We infer the properties of individual objects with two spectral energy distribution modelling codes, then infer a distribution of ages for bright $z\sim8.5-11$ galaxies. We find a median age of $\sim20$ Myr, younger than that inferred at $z\sim7$ with a similar analysis, consistent with an evolution towards larger specific star formation rates at early times. The age distribution suggests that only $\sim3$ percent of bright $z\sim8.5-11$ galaxies would be similarly luminous at $z\gtrsim15$, implying that the number density of bright galaxies declines by at least an order of magnitude between $z\sim8.5-11$ and $z\sim15$. This evolution is challenging to reconcile with some early JWST results suggesting the abundance of bright galaxies does not significantly decrease towards very early times, but we suggest this tension may be eased if young stellar populations form on top of older stellar components, or if bright $z\sim15$ galaxies are observed during a burst of star formation., 16 pages, 10 figures, 1 appendix, accepted to MNRAS

The onset of the {\it JWST}-era provides a much-improved opportunity to characterize the resolved structure of early star forming systems. Previous {\it Spitzer} observations of $z\gtrsim 6$ galaxies revealed the presence of old stars and luminous HII regions (via [OIII]+H$\beta$ emission), but the poor resolution stunted our ability to map their locations with respect to the star forming regions identified in the rest-UV. In this paper, we investigate the internal structure of 12 of the most luminous $z\simeq 6-8$ galaxies in the EGS field observed with recent {\it JWST}/NIRCam imaging. The systems appear clumpy in the rest-UV, with more than half of the light coming from $\simeq 10^7$ to 10$^{9}$ M$_\odot$ star forming complexes that are $\simeq 150$ - 480 pc in size. The clumps tend to be dominated by young stars (median = 36 Myr), but we also find large variations in clump ages within individual galaxies. The [OIII]+H$\beta$ EW varies significantly across individual galaxies (reflecting differences in stellar and gas properties), but the HII regions largely track the UV-bright complexes. Perhaps surprisingly, the rest-optical continuum is just as clumpy as the UV, and we do not find older (and redder) nuclear stellar components that were previously undetected or faint in the UV. The majority of the stellar mass in bright $6, Comment: main text 12 pages, 7 figures; accepted for publication by MNRAS

We report the discovery of a double-peaked Lyman-$\alpha$ (Ly$\alpha$) emitter (LAE) at $z=3.2177\pm0.0001$ in VLT/MUSE data. The galaxy is strongly lensed by the galaxy cluster RXC~J0018.5+1626 recently observed in the RELICS survey, and the double-peaked Ly$\alpha$ emission is clearly detected in the two counter images in the MUSE field-of-view. We measure a relatively high Ly$\alpha$ rest-frame equivalent width (EW) of $\mathrm{EW}_{\mathrm{Ly}\alpha,0}=(63\pm2)\,\mathring{\mathrm{A}}$. Additional near-infrared (NIR) spectroscopy allows us to measure the H$\beta$, [OIII]$\lambda4959\,\mathring{\mathrm{A}}$ and [OIII]$\lambda5007\,\mathring{\mathrm{A}}$ emission lines, which show moderate rest-frame EWs of the order of a few $\sim10-100\,\mathring{\mathrm{A}}$, an [OIII]$\lambda5007\,\mathring{\mathrm{A}}$/H$\beta$ ratio of $4.8\pm0.7$, and a lower limit on the [OIII]/[OII] ratio of $>5.6$. The galaxy has very blue UV-continuum slopes of $\beta_{\mathrm{FUV}}=-2.23\pm0.06$ and $\beta_{\mathrm{NUV}}=-3.0\pm0.2$, and is magnified by factors $\mu\sim7-10$ in each of the two images, thus enabling a view into a low-mass ($M_{\star}\simeq10^{7.5}\,\mathrm{M}_{\odot}$) high-redshift galaxy analog. Notably, the blue peak of the Ly$\alpha$ profile is significantly stronger than the red peak, which suggests an inflow of matter and possibly very low HI column densities in its circumgalactic gas. Combined with the high lensing magnification and image multiplicity, these properties make this galaxy a prime candidate for follow-up observations to search for LyC emission and constrain the LyC photon escape fraction., Comment: Updated to match published version

We present a Keck/MOSFIRE, rest-optical, composite spectrum of 16 typical, gravitationally-lensed, star-forming, dwarf galaxies at $1.7 \lesssim z \lesssim 2.6$ ($z_{\rm{mean}}=2.30$), all chosen independent of emission-line strength. These galaxies have a median stellar mass of $\log(M_\ast/\rm{M_\odot})_{\rm{med}} = 8.29^{+0.51}_{-0.43}$ and a median star formation rate of $\rm{SFR_{H\alpha}^{med} = 2.25^{+2.15}_{-1.26}\ M_\odot\ yr^{-1}}$. We measure the faint, electron-temperature-sensitive, [O III] $\lambda$4363 emission line at $2.5\sigma$ ($4.1\sigma$) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of $12+\log(\rm{O/H})_{\rm{direct}}=7.88^{+0.25}_{-0.22}$ ($0.15^{+0.12}_{-0.06}\ \rm{Z_\odot}$). We investigate the applicability at high-$z$ of locally-calibrated, oxygen-based, strong-line metallicity relations, finding that the local reference calibrations of arXiv:1805.08224 best reproduce ($\lesssim 0.12$ dex) our composite metallicity at fixed strong-line ratio. At fixed $M_\ast$, our composite is well-represented by the $z \sim 2.3$ direct-method stellar mass$\,-\,$gas-phase metallicity relation (MZR) of arXiv:1907.00013. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, at fixed $M_\ast$ and SFR, of the locally-defined fundamental metallicity relation. We measure the doublet ratio [O II] $\lambda$3729/[O II] $\lambda3726 = 1.56 \pm 0.32$ ($1.51 \pm 0.12$) and a corresponding electron density of $n_e = 1^{+215}_{-0}\ \rm{cm^{-3}}$ ($n_e = 1^{+74}_{-0}\ \rm{cm^{-3}}$) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies at $z \sim 2$., Comment: 30 pages, 6 figures, 4 tables, abstract slightly abridged, use of PyNeb adopted (replacing use of IRAF) to calculate electron temperatures and densities; change in results is negligible, accepted for publication in ApJ

The JWST observations of high-redshift galaxies are used to measure their star formation histories—the buildup of stellar mass in the earliest galaxies. Here we use a novel analysis program, SEDz*, to compare near-IR spectral energy distributions for galaxies with redshifts 5 < z < 7 to combinations of stellar population templates evolved from z = 12. We exploit NIRCam imaging in seven wide bands covering 1–5 μm taken in the context of the GLASS-JWST-ERS program and use SEDz* to solve for well-constrained star formation histories for 24 exemplary galaxies. In this first look, we find a variety of histories, from long, continuous star formation over 5 < z < 12 to short but intense starbursts, sometimes repeating, and, most commonly, contiguous mass buildup lasting ∼0.5 Myr, possibly the seeds of today’s typical M* galaxies.

Support by grant 2020750 from the United States- Israel Binational Science Foundation (BSF) and grant 2109066 from the United States National Science Foundation (NSF), and by the Ministry of Science and Technology, Israel. J.C. acknowl- edges funding from the ‘FirstGalaxies’ Advanced Grant from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (Grant agreement No. 789056). E.C.L. acknowledges support of an STFC Webb Fellowship (ST/W001438/1). K.K. acknowledges the support by JSPS KAKENHI Grant Number JP17H06130 and the NAOJ ALMA Scientific Research Grant Number 2017–06B. D.E. acknowledges support from a Beatriz Galindo senior fellowship (BG20/00224) from the Spanish Ministry of Science and Innovation, projects PID2020-114414GB-100 and PID2020-113689GB-I00 financed by MCIN/AEI/10.13039/501100011033, project P20-00334 financed by the Junta de Andaluc´ıa, and project A-FQM-510-UGR20 of the FEDER/Junta de Andaluc´ıa-Consejer´ıa de Transformaci ´on Econ ´omica, Industria, Conocimiento y Universidades. G.E.M. ac- knowledges financial support from the Villum Young Investiga- tor grant 37440 and 13160 and the The Cosmic Dawn Center (DAWN), funded by the Danish National Research Foundation under grant No. 140. F.E.B. acknowledges support from ANID- Chile BASAL CATA FB210003, FONDECYT Regular 1200495 and 1190818, and Millennium Science Initiative Program–ICN12 009. K.K.K. acknowledges support from the Knut and Alice Wallenberg Foundation., We report the discovery of a triply imaged active galactic nucleus (AGN), lensed by the galaxy cluster MACS J0035.4−2015 (z d = 0.352). The object is detected in Hubble Space Telescope imaging taken for the RELICS program. It appears to have a quasi-stellar nucleus consistent with a point-source, with a de-magnified radius of re ≲ 100 pc. The object is spectroscopically confirmed to be an AGN at z spec = 2.063 ± 0.005 showing broad rest-frame UV emission lines, and detected in both X-ray observations with Chandra and in ALCS ALMA band 6 (1.2 mm) imaging. It has a relatively faint rest-frame UV luminosity for a quasar-like object, MUV, 1450 = −19.7 ± 0.2. The object adds to just a few quasars or other X-ray sources known to be multiply lensed by a galaxy cluster. Some diffuse emission from the host galaxy is faintly seen around the nucleus, and there is a faint object nearby sharing the same multiple-imaging symmetry and geometric redshift, possibly an interacting galaxy or a star-forming knot in the host. We present an accompanying lens model, calculate the magnifications and time delays, and infer the physical properties of the source. We find the rest-frame UV continuum and emission lines to be dominated by the AGN, and the optical emission to be dominated by the host galaxy of modest stellar mass M✶ ≃ 109.2 M⊙. We also observe some variation in the AGN emission with time, which may suggest that the AGN used to be more active. This object adds a low-redshift counterpart to several relatively faint AGN recently uncovered at high redshifts with HST and JWST., Horizon 2020 Framework Programme H2020, Consejería de Transformación Económica, Industria, Conocimiento y Universidades 13160, 37440, Science and Technology Facilities Council ST/W001438/1 STFC, European Research Council ERC, Fondo Nacional de Desarrollo Científico y Tecnológico 1190818, 1200495, ICN12_009 FONDECYT, Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033, P20-00334, PID2020-113689GB-I00, PID2020-114414GB-100 MICINN, Horizon 2020 789056, European Regional Development Fund ERDF, Junta de Andalucía A-FQM-510-UGR20

Reionisation-era galaxies often display intense nebular emission lines, both in rest-frame optical ([OIII]+H$\beta$) and ultraviolet (UV; CIII], CIV). How such strong nebular emission is powered remains unclear, with both active galactic nuclei (AGN) and hot stars considered equally viable. The UV continuum slopes of these early systems tend to be very blue ($\beta-1$) than typical star-forming systems in the reionisation era. To investigate the properties of AGN in the reionisation era, we have conducted a search for potential examples of rare analogues with blue continua at intermediate redshift ($z\sim2-3$). Our goals are to determine whether AGN with intense line emission and blue continua exist and thereby to establish the range of rest-frame UV and optical line ratios in this population. In this paper we report the detection of a X-ray luminous AGN at $z=3.21$ (UDS-24561) with extreme [OIII]+H$\beta$ line emission (EW $=1300$ \r{A}) and a blue UV continuum slope ($\beta=-2.34$). MMT/Binospec and Keck/MOSFIRE spectra indicate rest-frame UV line ratios consistent with AGN photoionisation models and rest-frame optical lines with both a narrow component (FWHM $=154$ km$/$s) and extended broad wings (FWHM $=977$ km$/$s), consistent with outflowing gas. We describe how such objects can be identified in future JWST emission line surveys in the reionisation era, thereby providing a valuable census of AGN activity at $z>6$ and understanding their contribution to cosmic reionisation., Comment: 11 pages, 7 figures, accepted by MNRAS

We present ultra-deep Keck/MOSFIRE rest-optical spectra of two star-forming galaxies at z=2.18 in the COSMOS field with bright emission lines, representing more than 20~hours of total integration. The fidelity of these spectra enabled the detection of more than 20 unique emission lines for each galaxy, including the first detection of the auroral [O II]$\lambda\lambda$7322,7332 lines at high redshift. We use these measurements to calculate the electron temperature in the low-ionization O$^+$ zone of the ionized ISM and derive abundance ratios of O/H, N/H, and N/O using the direct method. The N/O and $\alpha$/Fe abundance patterns of these galaxies are consistent with rapid formation timescales and ongoing strong starbursts, in accord with their high specific star-formation rates. These results demonstrate the feasibility of using auroral [O II] measurements for accurate metallicity studies at high redshift in a higher metallicity regime previously unexplored with the direct method in distant galaxies. These results also highlight the difficulty in obtaining the measurements required for direct-method metallicities from the ground. We emphasize the advantages that the JWST/NIRSpec instrument will bring to high-redshift metallicity studies, where the combination of increased sensitivity and uninterrupted wavelength coverage will yield more than an order of magnitude increase in efficiency for multiplexed auroral-line surveys relative to current ground-based facilities. Consequently, the advent of JWST promises to be the beginning of a new era of precision chemical abundance studies of the early universe at a level of detail rivaling that of local galaxy studies., Comment: 18 pages, 5 figures, 3 tables, submitted to ApJ

Galaxy formation and evolution are regulated by the feedback from galactic winds. Absorption lines provide the most widely available probe of winds. However, since most data only provide information integrated along the line of sight, they do not directly constrain the radial structure of the outflows. In this paper, we present a method to directly measure the gas electron density in outflows (n e ), which in turn yields estimates of outflow cloud properties (e.g., density, volume filling factor, and sizes/masses). We also estimate the distance (r n ) from the starburst at which the observed densities are found. We focus on 22 local star-forming galaxies primarily from the COS Legacy Archive Spectroscopic SurveY (CLASSY). In half of them, we detect absorption lines from fine-structure excited transitions of Si ii (i.e., Si ii*). We determine n e from relative column densities of Si ii and Si ii*, given Si ii* originates from collisional excitation by free electrons. We find that the derived n e correlates well with the galaxy’s star formation rate per unit area. From photoionization models or assuming the outflow is in pressure equilibrium with the wind fluid, we get r n ∼ 1–2r * or ∼5r *, respectively, where r * is the starburst radius. Based on comparisons to theoretical models of multiphase outflows, nearly all of the outflows have cloud sizes large enough for the clouds to survive their interaction with the hot wind fluid. Most of these measurements are the first ever for galactic winds detected in absorption lines and, thus, will provide important constraints for future models of galactic winds.

We present band 6 ALMA observations of a heavily-obscured radio-loud ($L_{1.4\ \mathrm{GHz}}=10^{25.4}$ W Hz$^{-1}$) AGN candidate at $z_\mathrm{phot}=6.83\pm0.06$ found in the 1.5 deg$^2$ COSMOS field. The ALMA data reveal detections of exceptionally strong [CII]158$\mu$m ($z_\mathrm{[CII]}=6.8532$) and underlying dust continuum emission from this object (COS-87259), where the [CII] line luminosity, line width, and 158$\mu$m continuum luminosity are comparable to that seen from $z\sim7$ sub-mm galaxies and quasar hosts. The 158$\mu$m continuum detection suggests a total infrared luminosity of $9\times10^{12}$ $L_\odot$ with corresponding very large obscured star formation rate (1300 $M_\odot$/yr) and dust mass ($2\times10^9$ $M_\odot$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization era galaxy with $M_\ast\approx1.7\times10^{11}$ $M_\odot$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbors an AGN that is heavily obscured ($\tau_{_{\mathrm{9.7\mu m}}}=2.3$) with a bolometric luminosity of approximately $5\times10^{13}$ $L_\odot$. Such a very high AGN luminosity suggests this object is powered by an $\approx$1.6 $\times$ 10$^9$ $M_\odot$ black hole if accreting near the Eddington limit, and is effectively a highly-obscured version of an extremely UV-luminous ($M_{1450}\approx-27.3$) $z\sim7$ quasar. Notably, these $z\sim7$ quasars are an exceedingly rare population ($\sim$0.001 deg$^{-2}$) while COS-87259 was identified over a relatively small field. Future very wide-area surveys with, e.g., Roman and Euclid have the potential to identify many more extremely red yet UV-bright $z\gtrsim7$ objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population., Comment: 12 pages, 7 figures, 1 table. Updated to accepted version (MNRAS)

Deep spectroscopy of galaxies in the reionization-era has revealed intense CIII] and CIV line emission (EW $>15-20$ \r{A}). In order to interpret the nebular emission emerging at $z>6$, we have begun targeting rest-frame UV emission lines in galaxies with large specific star formation rates (sSFRs) at $1.36$ only in large sSFR galaxies with [OIII]+H$\beta$ EW $>1500$ \r{A}. In contrast to previous studies, we find that many galaxies with intense [OIII] have weak CIII] emission (EW $=5-8$ \r{A}), suggesting that the radiation field associated with young stellar populations is not sufficient to power strong CIII]. Photoionization models demonstrate that the spread in CIII] among systems with large sSFRs ([OIII]+H$\beta$ EW $>1500$ \r{A}) is driven by variations in metallicity, a result of the extreme sensitivity of CIII] to electron temperature. We find that the strong CIII] emission seen at $z>6$ (EW $>15$ \r{A}) requires metal poor gas ($\simeq0.1\ Z_\odot$) whereas the weaker CIII] emission in our sample tends to be found at moderate metallicities ($\simeq0.3\ Z_\odot$). The luminosity distribution of the CIII] emitters in our $z\simeq1-3$ sample presents a consistent picture, with stronger emission generally linked to low luminosity systems ($M_{\rm{UV}}>-19.5$) where low metallicities are more likely. We quantify the fraction of strong CIII] and CIV emitters at $z\simeq1-3$, providing a baseline for comparison against $z>6$ samples. We suggest that the first UV line detections at $z>6$ can be explained if a significant fraction of the early galaxy population is found at large sSFR ($>200$ Gyr$^{-1}$) and low metallicity ($, Comment: 21 pages, 10 figures, accepted by MNRAS

Connecting the observed rest-ultraviolet (UV) luminosities of high-$z$ galaxies to their intrinsic luminosities (and thus star formation rates) requires correcting for the presence of dust. We bypass a common dust-correction approach that uses empirical relationships between infrared (IR) emission and UV colours, and instead augment a semi-empirical model for galaxy formation with a simple -- but self-consistent -- dust model and use it to jointly fit high-$z$ rest-UV luminosity functions (LFs) and colour-magnitude relations ($M_{\mathrm{UV}}$-$\beta$). In doing so, we find that UV colours evolve with redshift (at fixed UV magnitude), as suggested by observations, even in cases without underlying evolution in dust production, destruction, absorption, or geometry. The observed evolution in our model arises due to the reduction in the mean stellar age and rise in specific star formation rates with increasing $z$. The UV extinction, $A_{\mathrm{UV}}$, evolves similarly with redshift, though we find a systematically shallower relation between $A_{\mathrm{UV}}$ and $M_{\mathrm{UV}}$ than that predicted by IRX-$\beta$ relationships derived from $z \sim 3$ galaxy samples. Finally, assuming that high $1600 \r{A}$ transmission ($\gtrsim 0.6$) is a reliable LAE indicator, modest scatter in the effective dust surface density of galaxies can explain the evolution both in $M_{\mathrm{UV}}$-$\beta$ and LAE fractions. These predictions are readily testable by deep surveys with the James Webb Space Telescope., Comment: 14+4 pages, 11+5 figures, accepted for publication in MNRAS

Despite significant progress both observationally and theoretically, the origin of high-ionization nebular HeII emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing $\mathrm{He^+}$-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input SEDs containing the joint contribution of both stellar atmospheres and a multi-color disk model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power HeII, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of $10^{42}$ erg/s per unit star formation rate (SFR). To test this, we assemble a sample of eleven low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and HeII emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed HeII strengths, with typical X-ray luminosities or upper limits thereon of only $10^{40}$-$10^{41}$ erg/s per SFR. This indicates that HMXBs are not the dominant source of $\mathrm{He^+}$ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity., Comment: 18 pages, 9 figures, published in MNRAS

As the James Webb Space Telescope approaches scientific operation, there is much interest in exploring the redshift range beyond that accessible with Hubble Space Telescope imaging. Currently, the only means to gauge the presence of such early galaxies is to age-date the stellar population of systems in the reionisation era. As a significant fraction of $z\simeq7-8$ galaxies are inferred from Spitzer photometry to have extremely intense [OIII] emission lines, it is commonly believed these are genuinely young systems that formed at redshifts $z100$ Myr) population whose stellar mass is $\simeq40$ times that associated with the starburst responsible for the extreme line emission. Without adequate rest-frame near-infrared photometry we argue it may be premature to conclude that extreme line emitters in the reionisation era are low mass systems that formed at redshifts below $z\simeq10$., 14 pages, 7 figure, accepted by MNRAS

The variety of star formation histories (SFHs) of $z\gtrsim6$ galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many $z\sim6-9$ galaxies are dominated by $\gtrsim200$ Myr stellar populations, implying significant star formation at $z\gtrsim9$, while others find that most reionization era galaxies are $\lesssim10$ Myr, consistent with little $z\gtrsim9$ star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{UV}\lesssim-21$) galaxies colour-selected to lie at $z\simeq6.6-6.9$, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of $\sim30-70$ Myr depending on SED model. We quantify the fractions of $\leq10$ Myr and $\geq250$ Myr galaxies, finding that these systems comprise $\sim15-30$ percent and $\sim20-25$ percent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for $z\gtrsim9$ stellar mass assembly and discuss improvements expected from JWST., Comment: 24 pages, 20 figures, 2 appendices, accepted to MNRAS

The ultraviolet (UV) continuum slope ($\beta~$where$~f_\lambda\propto\lambda^\beta$) of galaxies is sensitive to a variety of properties, from the metallicity and age of the stellar population to dust attenuation throughout the galaxy. Considerable attention has focused on identifying reionization-era galaxies with very blue UV slopes ($\beta7$., Comment: 13 pages, 10 figures, 3 tables; Accepted to ApJ

Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly α transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly α profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at z ∼ 7 remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [C ii]-based redshifts from ALMA. Here, we report MMT/Binospec Ly α spectroscopy of eight UV-bright (MUV ∼ −22) galaxies at z ≃ 7 selected from the ALMA REBELS survey. We detect Ly α in four of eight galaxies and use the [C ii] systemic redshifts to investigate the Ly α velocity profiles. The Ly α lines are significantly redshifted from systemic (average velocity offset = 223 km s–1) and broad (FWHM ≈ 300–650 km s−1), with two sources showing emission extending to ≈750 km s−1. We find that the broadest Ly α profiles are associated with the largest [C ii] line widths, suggesting a potential link between the Ly α FWHM and the dynamical mass. Since Ly α photons at high velocities transmit efficiently through the z = 7 IGM, our data suggest that velocity profiles play a significant role in boosting the Ly α visibility of the most UV-luminous reionization-era galaxies.

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6, 18 pages, 10 figures, submitted to MNRAS

We report the identification of radio (0.144-3 GHz), mid-IR, far-IR, and sub-mm (24-850$\mu$m) emission at the position of one of 41 UV-bright (M$_\mathrm{UV}^{}\lesssim-21.25$) $z\simeq6.6-6.9$ Lyman-break galaxy candidates in the 1.5 deg$^2$ COSMOS field. This source, COS-87259, exhibits a sharp flux discontinuity (factor $>$3) between two narrow/intermediate bands at 9450 and 9700 Angstroms and is undetected in all nine bands blueward of 9600 Angstroms, as expected from a Lyman-alpha break at $z\simeq6.8$. The full multi-wavelength (X-ray through radio) data of COS-87529 can be self-consistently explained by a very massive (M$_{\ast}=10^{10.8}$ M$_{\odot}$) and extremely red (rest-UV slope $\beta=-0.59$) $z\simeq6.8$ galaxy with hyperluminous infrared emission (L$_{\mathrm{IR}}=10^{13.6}$ L$_{\odot}$) powered by both an intense burst of highly-obscured star formation (SFR$\approx$1800 M$_{\odot}$ yr$^{-1}$) and an obscured ($\tau_{\mathrm{9.7\mu m}}=7.7\pm2.5$) radio-loud (L$_{\mathrm{1.4\ GHz}}\approx10^{25.4}$ W Hz$^{-1}$) AGN. The radio emission is compact (1.04$\pm$0.12 arcsec) and exhibits an ultra-steep spectrum between 1.32-3 GHz ($\alpha=-1.57^{+0.22}_{-0.21}$) that flattens at lower frequencies ($\alpha=-0.86^{+0.22}_{-0.16}$ between 0.144-1.32 GHz), consistent with known $z>4$ radio galaxies. We also demonstrate that COS-87259 may reside in a significant (11$\times$) galaxy overdensity at $z\simeq6.6-6.9$, as common for systems hosting radio-loud AGN. Nonetheless, a spectroscopic redshift will ultimately be required to establish the true nature of COS-87259 as we cannot yet completely rule out low-redshift solutions. If confirmed to lie at $z\simeq6.8$, the properties of COS-87259 would be consistent with a picture wherein AGN and highly-obscured star formation activity are fairly common among very massive (M$_{\ast}>10^{10}$ M$_{\odot}$) reionization-era galaxies., Comment: 14 pages, 6 figures. Accepted in MNRAS with minor revisions. This accepted version considers very recent data from LOFAR and MeerKAT which improve our analyses on the radio slope and luminosity

Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization., Comment: MNRAS Accepted (01 Nov 2022)

Our understanding of reionization has advanced considerably over the past decade, with several results now demonstrating that the IGM transitioned from substantially neutral at $z=7$ to largely reionized at $z=6$. However, little remains known about the sizes of ionized bubbles at $z\gtrsim7$ as well as the galaxy overdensities which drive their growth. Fortunately, rest-UV spectroscopic observations offer a pathway towards characterizing these ionized bubbles thanks to the resonant nature of Lyman-alpha photons. In a previous work, we presented Ly$\alpha$ detections from three closely-separated Lyman-break galaxies at $z\simeq6.8$, suggesting the presence of a large ($R>1$ physical Mpc) ionized bubble in the 1.5 deg$^2$ COSMOS field. Here, we present new deep Ly$\alpha$ spectra of ten UV-bright ($\mathrm{M}_{\mathrm{UV}}^{} \leq -20.4$) $z\simeq6.6-6.9$ galaxies in the surrounding area, enabling us to better characterize this potential ionized bubble. We confidently detect (S/N$>$7) Ly$\alpha$ emission at $z=6.701-6.882$ in nine of ten observed galaxies, revealing that the large-scale volume spanned by these sources (characteristic radius $R = 3.2$ physical Mpc) traces a strong galaxy overdensity ($N/\langle N\rangle \gtrsim 3$). Our data additionally confirm that the Ly$\alpha$ emission of UV-bright galaxies in this volume is significantly enhanced, with 40% (4/10) showing strong Ly$\alpha$ emission (equivalent width$>$25 $\mathrm{\mathring{A}}$) compared to the 8$-$9% found on average at $z\sim7$. The median Ly$\alpha$ equivalent width of our observed galaxies is also $\approx$2$\times$ that typical at $z\sim7$, consistent with expectations if a very large ($R\sim3$ physical Mpc) ionized bubble is allowing the Ly$\alpha$ photons to cosmologically redshift far into the damping wing before encountering HI., Comment: 13 pages, 3 figures, 3 tables. Accepted in MNRAS with minor revisions

We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry., 35 pages, 18 figures

Metal-poor nearby galaxies hosting massive stars have a fundamental role to play in our understanding of both high-redshift galaxies and low metallicity stellar populations. But while much attention has been focused on their bright nebular gas emission, the massive stars that power it remain challenging to constrain. Here we present exceptionally deep Hubble Space Telescope ultraviolet spectra targeting six galaxies that power strong nebular C IV emission approaching that encountered at $z>6$. We find that the strength and spectral profile of the nebular C IV in these new spectra follow a sequence evocative of resonant scattering models, indicating that the hot circumgalactic medium likely plays a key role in regulating C IV escape locally. We constrain the metallicity of the massive stars in each galaxy by fitting the forest of photospheric absorption lines, reporting measurements driven by iron that lie uniformly below 10% solar. Comparison with the gas-phase oxygen abundances reveals evidence for enhancement in O/Fe above solar across the sample, robust to assumptions about the absolute gas-phase metallicity scale. This supports the idea that these local systems are more chemically-similar to their primordial high-redshift counterparts than to the bulk of nearby galaxies. Finally, we find significant tension between the strong stellar wind profiles observed and our population synthesis models constrained by the photospheric forest in our highest-quality spectra. This reinforces the need for caution in interpreting wind lines in isolation at high-redshift, but also suggests a unique path towards validating fundamental massive star physics at extremely low metallicity with integrated ultraviolet spectra., 34 pages, 14 figures, submitted to AAS Journals

As deep spectroscopic campaigns extend to higher redshifts and lower stellar masses, the interpretation of galaxy spectra depends increasingly upon models for very young stellar populations. Here we present new HST/COS ultraviolet spectroscopy of seven nearby ($, 25 pages, 9 figures; submitted to MNRAS

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with COS data. The outflows are traced by blueshifted absorption-lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static ISM from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow's mean velocity and velocity width, and find that both correlate in a highly significant way with the star-formation rate, galaxy mass, and circular velocity over ranges of four orders-of-magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried in the warm phase we observe. The outflows' mass-loading factor increases steeply and inversely with both circular and outflow velocity (log-log slope $\sim$ -1.6), and reaches $\sim 10$ for dwarf galaxies. We find that the outflows typically carry about 10 to 100% of the momentum injected by massive stars and about 1 to 20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds., Comment: 34 pages, 16 figures, 6 tables, submitted to ApJ

Spectroscopic observations of massive $z>7$ galaxies selected to have extremely large [OIII]+H$\beta$ equivalent width (EW $\sim1500$ \r{A}) have recently revealed large Ly$\alpha$ detection rates, in contrast to the weak emission seen in the general population. Why these systems are uniquely visible in Ly$\alpha$ at redshifts where the IGM is likely significantly neutral is not clear. With the goal of better understanding these results, we have begun a campaign with MMT and Magellan to measure Ly$\alpha$ in galaxies with similar [OIII]+H$\beta$ EWs at $z\simeq2-3$. At these redshifts, the IGM is highly ionized, allowing us to clearly disentangle how the Ly$\alpha$ properties depend on the [OIII]+H$\beta$ EW. Here we present Ly$\alpha$ EWs of $49$ galaxies at $z=2.2-3.7$ with intense [OIII]+H$\beta$ line emission (EW $=300-3000$ \r{A}). Our results demonstrate that strong Ly$\alpha$ emission (EW $>20$ \r{A}) becomes more common in galaxies with larger [OIII]+H$\beta$ EW, reflecting a combination of increasingly efficient ionizing photon production and enhanced transmission of Ly$\alpha$. Among the galaxies with the most extreme [OIII]+H$\beta$ emission (EW $\sim1500$ \r{A}), we find that strong Ly$\alpha$ emission is not ubiquitous, with only $50$ per cent of our population showing Ly$\alpha$ EW $>20$ \r{A}. Our data suggest that the range of Ly$\alpha$ strengths is related to the observed ellipticity, with those systems that appear edge-on or elongated having weaker Ly$\alpha$ emission. We use these results to interpret the anomalous Ly$\alpha$ properties seen in $z>7$ galaxies with extreme [OIII]+H$\beta$ emission and discuss implications for the escape of ionizing radiation from these extreme line emitting galaxies.

Reionization-era galaxies tend to exhibit weak Ly$\alpha$ emission, likely reflecting attenuation from an increasingly neutral IGM. Recent observations have begun to reveal exceptions to this picture, with strong Ly$\alpha$ emission now known in four of the most massive z=7$-$9 galaxies in the CANDELS fields, all of which also exhibit intense [OIII]$+$H$\beta$ emission (EW$>$800 $\mathrm{\mathring{A}}$). To better understand why Ly$\alpha$ is anonymously strong in a subset of massive z$\simeq$7$-$9 galaxies, we have initiated an MMT/Binospec survey targeting a larger sample (N=22) of similarly luminous ($\simeq$1$-$6 L$^{\ast}_{\mathrm{UV}}$) z$\simeq$7 galaxies selected over very wide-area fields ($\sim$3 deg$^2$). We confidently ($>$7$\sigma$) detect Ly$\alpha$ in 78% (7/9) of galaxies with strong [OIII]$+$H$\beta$ emission (EW$>$800 $\mathrm{\mathring{A}}$) as opposed to only 8% (1/12) of galaxies with more moderate (EW=200$-$800 $\mathrm{\mathring{A}}$) [OIII]$+$H$\beta$. We argue that the higher Ly$\alpha$ EWs of the strong [OIII]$+$H$\beta$ population likely reflect enhanced ionizing photon production efficiency owing to their large sSFRs ($\gtrsim$30 Gyr$^{-1}$). We also find evidence that Ly$\alpha$ transmission from massive galaxies declines less rapidly over $6, Comment: 20 pages, 11 figures. Accepted in MNRAS with minor revisions

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5�� allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar), 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$\sigma$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

We measure the ionizing photon production efficiency ($\xi_{ion}$) of low-mass galaxies ($10^{7.8}$-$10^{9.8}$ $M_{\odot}$) at $1.4, Comment: 16 pages, 6 figures, 1 table, Accepted in ApJ

Galaxy evolution is thought to be driven in large part by the flow of gas between galaxies and the circumgalactic medium (CGM), a halo of metal-enriched gas extending out to $\gtrsim100$ kpc from each galaxy. Studying the spatial structure of the CGM holds promise for understanding these gas flow mechanisms; however, the common method using background quasar sightlines provides minimal spatial information. Recent works have shown the utility of extended background sources such as giant gravitationally lensed arcs. Using background lensed arcs from the CSWA 38 lens system, we continuously probed, at a resolution element of about 15 kpc$^2$, the spatial and kinematic distribution of MgII absorption in a star-forming galaxy at $z=0.77$ (stellar mass $\approx 10^{9.7}$ M$_\odot$, star formation rate $\approx 10$ M$_\odot$ yr$^{-1}$) at impact parameters $D \simeq 5-30$ kpc. Our results present an anisotropic, optically thick medium whose absorption strength decreases with increasing impact parameter, in agreement with the statistics towards quasars and other gravitational arcs. Furthermore, we find generally low line-of-sight velocities in comparison to the relatively high velocity dispersion in the MgII gas (with typical $\sigma\approx 50$ km s$^{-1}$). While the galaxy itself exhibits a clear outflow (with MgII velocities up to $\sim 500$ km s$^{-1}$) in the down-the-barrel spectrum, the outflow component is sub-dominant and only weakly detected at larger impact parameters probed by the background arcs. Our results provide evidence of mainly dispersion-supported, metal-enriched gas recycling through the CGM., Comment: 27 pages, 10 figures, 5 tables

We quantify the distribution of [OIII]+H$\beta$ line strengths at z$\simeq$7 using a sample of 20 bright (M$_{\mathrm{UV}}$ $\lesssim$ $-$21) galaxies. We select these systems over wide-area fields (2.3 deg$^2$ total) using a new colour-selection which precisely selects galaxies at z$\simeq$6.63$-$6.83, a redshift range where blue Spitzer/IRAC [3.6]$-$[4.5] colours unambiguously indicate strong [OIII]$+$H$\beta$ emission. These 20 galaxies suggest a log-normal [OIII]$+$H$\beta$ EW distribution with median EW = 759$^{+112}_{-113}$ $\mathrm{\mathring{A}}$ and standard deviation = 0.26$^{+0.06}_{-0.05}$ dex. We find no evidence for strong variation in this EW distribution with UV luminosity. The typical [OIII]+H$\beta$ EW at z$\simeq$7 implied by our sample is considerably larger than that in massive star forming galaxies at z$\simeq$2, consistent with a shift toward larger average sSFR (4.4 Gyr$^{-1}$) and lower metallicities (0.16 Z$_\odot$). We also find evidence for the emergence of a population with yet more extreme nebular emission ([OIII]+H$\beta$ EW$>$1200 $\mathrm{\mathring{A}}$) that is rarely seen at lower redshifts. These objects have extremely large sSFR ($>$30 Gyr$^{-1}$), as would be expected for systems undergoing a burst or upturn in star formation. While this may be a short-lived phase, our results suggest that 20% of the z$\simeq$7 population has such extreme nebular emission, implying that galaxies likely undergo intense star formation episodes regularly at z$>$6. We argue that this population may be among the most effective ionizing agents in the reionization era, both in terms of photon production efficiency and escape fraction. We furthermore suggest that galaxies passing through this large sSFR phase are likely to be very efficient in forming bound star clusters., Comment: 20 pages, 11 figures. Accepted in MNRAS with minor revisions

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the $z\sim 6\mbox{--}8$ candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ≈200 arcmin2. These clusters were selected to be excellent lenses, and we find similar high-redshift sample sizes and magnitude distributions as the Cluster Lensing And Supernova survey with Hubble (CLASH). We discover 257, 57, and eight candidate galaxies at z ∼ 6, 7, and 8 respectively, (322 in total). The observed (lensed) magnitudes of the z ∼ 6 candidates are as bright as AB mag ∼23, making them among the brightest known at these redshifts, comparable with discoveries from much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Among the nearest most metal-poor starburst-dwarf galaxies known, SBS 0335-052E is the most luminous in integrated nebular He II {\lambda}4686 emission. This makes it a unique target to test spectral synthesis models and spectral interpretation tools of the kind that will be used to interpret future rest-frame UV observations of primeval galaxies. Previous attempts to reproduce its He II {\lambda}4686 emission luminosity found that X-ray sources, shocks, and single Wolf-Rayet stars are not main contributors to the He II-ionizing budget; and that only metal-free single rotating stars or binary stars with a top-heavy IMF and an unphysically-low metallicity can reproduce it. We present new UV (COS) and optical (MUSE) spectra which integrate the light of four super star clusters in SBS 0335-052E. Nebular He II, [C III], C III], C IV, and O III] UV emission lines with equivalent widths between 1.7 and 5 {\AA}, and a C IV {\lambda}{\lambda}1548, 1551 P-Cygni like profile are detected. Recent extremely-metal poor shock + precursor models and binary models fail to reproduce the observed optical emission-line ratios. We use different sets of UV and optical observables to test models of constant star formation with single non-rotating stars which account for very massive stars, as blueshifted O V {\lambda}1371 absorption is present. Simultaneously fitting the fluxes of all high-ionization UV lines requires an unphysically-low metallicity. Fitting the P-Cygni like + nebular components of C IV {\lambda}{\lambda}1548, 1551 does not constrain the stellar metallicity and time since the beginning of star formation. We obtain 12+log(O/H)=7.45\pm0.04 and log(C/O)=-0.45(+0.03)(-0.04) for the galaxy. Model-testing would benefit from higher spatial resolution UV and optical spectroscopy of the galaxy., Comment: 20 pages, 12 figures, accepted for publication in MNRAS

Escaping Lyman continuum photons from galaxies likely reionized the intergalactic medium at redshifts $z\gtrsim6$. However, the Lyman continuum is not directly observable at these redshifts and secondary indicators of Lyman continuum escape must be used to estimate the budget of ionizing photons. Observationally, at redshifts $z\sim2-3$ where the Lyman continuum is observationally accessible, surveys have established that many objects that show appreciable Lyman continuum escape fractions $f_{esc}$ also show enhanced [OIII]/[OII] (O$_{32}$) emission line ratios. Here, we use radiative transfer analyses of cosmological zoom-in simulations of galaxy formation to study the physical connection between $f_{esc}$ and O$_{32}$. Like the observations, we find that the largest $f_{esc}$ values occur at elevated O$_{32}\sim3-10$ and that the combination of high $f_{esc}$ and low O$_{32}$ is extremely rare. While high $f_{esc}$ and O$_{32}$ often are observable concurrently, the timescales of the physical origin for the processes are very different. Large O$_{32}$ values fluctuate on short ($\sim$1 Myr) timescales during the Wolf-Rayet-powered phase after the formation of star clusters, while channels of low absorption are established over tens of megayears by collections of supernovae. We find that while there is no direct causal relation between $f_{esc}$ and O$_{32}$, high $f_{esc}$ most often occurs after continuous input from star formation-related feedback events that have corresponding excursions to large O$_{32}$ emission. These calculations are in agreement with interpretations of observations that large $f_{esc}$ tends to occur when O$_{32}$ is large, but large O$_{32}$ does not necessarily imply efficient Lyman continuum escape., Comment: 12 pages, 5 figures, submitted to ApJL on September 15, 2020

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data., Comment: 39 pages, 27 figures. Submitted to ApJ. Comments welcome

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation., Comment: 21 pages, 16 figures, the version accepted in MNRAS

We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$\lambda\lambda1907,1909$ doublet for three $z\sim8$ galaxy candidates in the AEGIS field. Along with non-detections in two galaxies, we obtain one of the highest-redshift detections to-date of [CIII]$\lambda 1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$\lambda$1907 flux of $2.24\pm0.71\times10^{-18} \mbox{ erg}\mbox{ s}^{-1} \mbox{ cm}^{-2}$, corresponding to a rest-frame equivalent width of $20.3\pm6.5 \unicode{x212B}$ for the single line. Given the not very constraining upper limit for CIII]$\lambda 1909$ based on strong sky-line contamination, we assume a [CIII]$\lambda$1907/CIII]$\lambda 1909$ doublet ratio of 1.5 and infer a total [CIII],CIII]$\lambda\lambda1907,1909$ equivalent width of $33.7\pm 10.8 \unicode{x212B}$. We repeat the same reductions and analysis on multiple subsets of our raw data divided on the basis of time and observing conditions, verifying that the [CIII]$\lambda 1907$ emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [CIII]$\lambda 1907$ detection in different subsets of our data tracks that of brighter emission features detected on the same multi-slit mask. These multiple tests suggest that the observed emission line is real and associated with the $z\sim 8$ target. The strong observed [CIII],CIII]$\lambda\lambda1907,1909$ in AEGIS-33376 likely indicates ISM conditions of low metallicity, high ionization parameter, and a hard ionizing spectrum, although AGN contributions are possible. This single detection represents a sizable increase in the current sample [CIII],CIII]$\lambda\lambda1907,1909$ detections at $z>7$, while $\textit{JWST}$ will provide the first statistical samples of such measurements at these redshifts., Comment: 6 pages, 3 figures, 1 table; Accepted to ApJL

Strong [OIII]$\lambda\lambda$4959,5007+H$\beta$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$\alpha$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$\lambda\lambda$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$\alpha$ and CIII]$\lambda\lambda$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$\alpha$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$\alpha$ emission starts to emerge. Using stacked spectra, we find that Ly$\alpha$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$\lambda\lambda$1548,1550, HeII$\lambda$1640, and OIII]$\lambda\lambda$1661,1665 emission in addition to intense Ly$\alpha$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts., Comment: 22 pages, 9 figures, accepted for publication in ApJ

Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$., Comment: Accepted for publication in MNRAS

We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization., 15 pages, 6 figures. Accepted for publication in ApJ

We present Keck/MOSFIRE and Keck/LRIS spectroscopy of A1689-217, a lensed (magnification $\sim7.9$), star-forming (SFR $\sim$ 16 $\rm M_\odot\ yr^{-1}$), dwarf (log($M_\ast/M_\odot$) = $8.07-8.59$), Ly$\alpha$-emitter ($\rm EW_0\sim138$ A) at $z$ = 2.5918. Dwarf galaxies similar to A1689-217 are common at high redshift and likely responsible for reionization, yet few have been studied with detailed spectroscopy. We report a 4.2$\sigma$ detection of the electron-temperature-sensitive [O III]$\lambda$4363 emission line and use this line to directly measure an oxygen abundance of 12+$\log$(O/H) = 8.06 $\pm$ 0.12 ($\sim1/4\ \rm Z_\odot$). A1689-217 is the lowest mass galaxy at $z>$ 2 with an [O III]$\lambda$4363 detection. Using the rest-optical emission lines, we measure A1689-217's other nebular conditions including electron temperature ($T_e$([O III]) $\sim$ 14,000 K), electron density ($n_e\sim$ 220 $\rm cm^{-3}$) and reddening ($E(B-V)\sim$ 0.39). We study relations between strong-line ratios and direct metallicities with A1689-217 and other galaxies with [O III]$\lambda$4363 detections at $z\sim0-3.1$, showing that the locally-calibrated, oxygen-based, strong-line relations are consistent from $z\sim0-3.1$. We also show additional evidence that the $\rm O_{32}$ vs. $\rm R_{23}$ excitation diagram can be utilized as a redshift-invariant, direct-metallicity-based, oxygen abundance diagnostic out to $z\sim$ 3.1. From this excitation diagram and the strong-line ratio $-$ metallicity plots, we observe that the ionization parameter at fixed O/H is consistent with no redshift evolution. Although A1689-217 is metal-rich for its $M_\ast$ and SFR, we find it to be consistent within the large scatter of the low-mass end of the Fundamental Metallicity Relation., Comment: 23 pages, 9 figures, 3 tables, 1 appendix - v2: Revised manuscript submitted to ApJ (added discussion/comments/refs), Rest-frame EW of H$\alpha$ and [O III]$\lambda$5007 corrected (results unaffected)

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST)., 20 pages, 8 figures. Accepted to ApJ

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$\sigma$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization., Comment: 17 pages, 7 figures

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches., Comment: 16 pages, 7 figures, 1 table. Accepted to ApJL

Large surveys of galaxy clusters with the Hubble Space Telescope (HST) and Spitzer, including the Cluster Lensing And Supernova survey with Hubble and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey, described here. Our 188-orbit Hubble Treasury Program observed 41 clusters at 0.182 ≤ z ≤ 0.972 with Advanced Camera for Surveys (ACS) and WFC3/IR imaging spanning 0.4-1.7 μm. We selected 21 of the most massive clusters known based on Planck PSZ2 estimates and 20 additional clusters based on observed or inferred lensing strength. RELICS observed 46 WFC3/IR pointings (∼200 arcmin) each with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered over 300 z ∼ 6-10 candidates, including the brightest z ∼ 6 candidates known, and the most distant spatially resolved lensed arc known at z ∼ 10. Spitzer IRAC imaging (945 hr awarded, plus 100 archival, spanning 3.0-5.0 μm) has crucially enabled us to distinguish z ∼ 10 candidates from z ∼ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. Reduced HST images, catalogs, and lens models are available on MAST, and reduced Spitzer images are available on IRSA.© 2019. The American Astronomical Society. All rights reserved., We thank Lindsey Bleem for providing Magellan Megacam LDSS3 images of SPT0254-58 and AS295 prior to RELICS to inform our HST observations of these clusters. We thank Florian Pacaud and Matthias Klein for discussions regarding our AS295 HST pointings. We thank Dale Kocevski for an image of RXC 0142+44 obtained with the University of Hawaii 2.2 m telescope. And we thank Stella Seitz et al. for sharing their HST observations of PLCK G287+32 obtained during the same cycle. We thank the STScI and SSC directors and time allocation committees for enabling these large observing programs. We are grateful to our HST program coordinator William Januszewski for implementing the RELICS HST observations. We thank Jennifer Mack for expert mentoring of our HST image reduction gurus RA and SO. And we thank Gabriel Brammer for providing an updated WFC3/IR hot pixel mask derived from science observations from GO 14114. The RELICS Hubble Treasury Program (GO 14096) consists of observations obtained by the NASA/ESA Hubble Space Telescope (HST). HST is operated by the Association of Universities for Research in Astronomy, Inc. (AURA), under NASA contract NAS5-26555. Data from the NASA/ESA HST presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST), operated by the Space Telescope Science Institute (STScI). STScI is operated by the Association of Universities for Research in Astronomy, Inc..(AURA), under NASA contract NAS 5-26555. The HST Advanced Camera for Surveys (ACS) was developed under NASA contract NAS 5-32864. Spitzer Space Telescope data presented in this paper were obtained from the NASA/IPAC Infrared Science Archive (IRSA), operated by the Jet Propulsion Laboratory, California Institute of Technology. Spitzer and IRSA are operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA. We gratefully acknowledge support from JPL for the Spitzer analysis. M.B. and V.S. also acknowledge support by NASA through ADAP grant 80NSSC18K0945, NASA/HST through HST-GO-14096, HST-GO-13666, and two awards issued by Spitzer/JPL/Caltech associated with the SRELICS_DEEP and SRELICS programs. Part of this work by W.D. was performed under the auspices of the U.S..DOE by LLNL under contract DE-AC5207NA27344. K.U. acknowledges support from the Ministry of Science and Technology of Taiwan (grant MOST 106-2628M-001-003-MY3) and from Academia Sinica (grant AS-IA107-M01). O.G. is supported by an NSF Astronomy and Astrophysics Fellowship under award AST-1602595. S.A.R. was supported by NASA grant HST-GO-14208 from STScI, which is operated by Associated Universities for Research in Astronomy, Inc. (AURA), under NASA contract NAS 5-26555. A.M. acknowledges the financial support of the Brazilian funding agency FAPESP (Post-doc fellowshipprocess No. 2014/11806-9). J.H. was supported by a VILLUM FONDEN Investigator grant (project No. 16599).

Galaxies in the reionization-era have been shown to have prominent [OIII]+H$\beta$ emission. Little is known about the gas conditions and radiation field of this population, making it challenging to interpret the spectra emerging at $z\gtrsim6$. Motivated by this shortcoming, we have initiated a large MMT spectroscopic survey identifying rest-frame optical emission lines in 227 intense [OIII] emitting galaxies at $1.3, Comment: 25 pages, 11 figures, accepted by MNRAS

We present constraints on the physical properties (including stellar mass, age, and star formation rate) of 207 $6\lesssim z \lesssim8$ galaxy candidates from the Reionization Lensing Cluster Survey (RELICS) and companion Spitzer-RELICS surveys. We measure photometry using T-PHOT and perform spectral energy distribution fitting using EA$z$Y and BAGPIPES. Of the 207 candidates for which we could successfully measure Spitzer fluxes, 23 were demoted to likely low redshift ($z, Comment: 20 pages, 9 figures