Search

JWST has revealed a large population of compact, red galaxies at $z>4$ known as Little Red Dots (LRDs). We analyze the spectral energy distributions (SEDs) of 95 LRDs from the JWST PRIMER survey with complete photometric coverage from $1-18\ \mu$m using NIRCam and MIRI imaging, representing the most extensive SED analysis on a large LRD sample with long-wavelength MIRI data. We examine SED models in which either galaxy or active galactic nucleus (AGN) emission dominates the rest-frame UV or optical continuum, extracting physical properties to explore each scenario's implications. In the galaxy-only model, we find massive, dusty stellar populations alongside unobscured, low-mass components, hinting at inhomogeneous obscuration. The AGN-only model indicates dusty, luminous AGNs with low hot dust fractions compared to typical quasars. A hybrid AGN and galaxy model suggests low-mass, unobscured galaxies in the UV, with stellar mass estimates spanning $\sim$2 dex across the different models, underscoring the need for caution in interpreting LRD stellar masses. With MIRI photometry, the galaxy-only model produces stellar masses within cosmological limits, but extremely high stellar mass densities are inferred. The hybrid model infers highly overmassive black holes exceeding those in recently reported high-redshift AGNs, hinting at a partial AGN contribution to the rest-optical continuum or widespread super-Eddington accretion. Our findings highlight the extreme conditions required for both AGN or galaxy dominated scenarios in LRDs, supporting a mixed contribution to the red continuum, or novel scenarios to explain the observed emission., Comment: 22 pages, 10 figures, submitted to ApJ. Machine-readable form of Table 2 available at: https://github.com/geneckleung/lrd_primer_miri

We present the properties of a massive, large, dusty, metal-rich, star-forming galaxy at z_spec=6.73. GOODSN-100182 was observed with JWST/NIRSpec as part of the AURORA survey, and is also covered by public multi-wavelength HST and JWST imaging. While the large mass of GOODSN-100182 (~10^10 M_sun) was indicated prior to JWST, NIRCam rest-optical imaging now reveals the presence of an extended disk (r_eff~1.5 kpc). In addition, the NIRSpec R~1000 spectrum of GOODSN-100182 includes the detection of a large suite of rest-optical nebular emission lines ranging in wavelength from [OII]3727 up to [NII]6583. The ratios of Balmer lines suggest significant dust attenuation (E(B-V)_gas=0.40+0.10/-0.09), consistent with the red rest-UV slope inferred for GOODSN-100182 (beta=-0.50+/-0.09). The star-formation rate based on dust-corrected H-alpha emission is log(SFR(H-alpha)/ M_sun/yr)=2.02+0.13/-0.14, well above the z~7 star-forming main sequence in terms of specific SFR. Strikingly, the ratio of [NII]6583/H-alpha emission suggests almost solar metallicity, as does the ratio ([OIII]5007/H-beta)/([NII]6583/H-alpha) and the detection of the faint [FeII]4360 emission feature, whereas the [OIII]5007/[OII]3727 ratio suggests roughly 50% solar metallicity. Overall, the excitation and ionization properties of GOODSN-100182 more closely resemble those of typical star-forming galaxies at z~2-3 rather than z~7. Based on public spectroscopy of the GOODS-N field, we find that GOODSN-100182 resides within a significant galaxy overdensity, and is accompanied by a spectroscopically-confirmed neighbor galaxy. GOODSN-100182 demonstrates the existence of mature, chemically-enriched galaxies within the first billion years of cosmic time, whose properties must be explained by galaxy formation models., Comment: 16 pages, 13 figures, submitted to ApJ

We present a sample of 50 H-alpha detected broad-line active galactic nuclei (BLAGN) at redshifts 3.50), independent of the contributions of emission lines to the broadband photometry. We construct the black hole (BH) mass function at 3.5

We use JWST/NIRSpec observations from the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) survey to constrain the shape of the nebular attenuation curve of a star-forming galaxy at z=4.41, GOODSN-17940. We utilize 11 unblended HI recombination lines to derive the attenuation curve spanning optical to near-infrared wavelengths (3751-9550 \r{A}). We then leverage a high-S/N spectroscopic detection of the rest-frame ultraviolet continuum in combination with rest-UV photometric measurements to constrain the shape of the curve at ultraviolet wavelengths. While this UV constraint is predominantly based on stellar emission, the large measured equivalent widths of H$\alpha$ and H$\beta$ indicate that GOODSN-17940 is dominated by an extremely young stellar population <10 Myr in age such that the UV stellar continuum experiences the same attenuation as the nebular emission. The resulting combined nebular attenuation curve spans 1400-9550 \r{A} and has a shape that deviates significantly from commonly assumed dust curves in high-redshift studies. Relative to the Milky Way, SMC, and Calzetti curves, the new curve has a steeper slope at long wavelengths ($\lambda>5000$ \r{A}) while displaying a similar slope across blue-optical wavelengths ($\lambda=3750-5000$ \r{A}). In the ultraviolet, the new curve is shallower than the SMC and Calzetti curves and displays no significant 2175 \r{A} bump. This work demonstrates that the most commonly assumed dust curves are not appropriate for all high-redshift galaxies. These results highlight the ability to derive nebular attenuation curves for individual high-redshift sources with deep JWST/NIRSpec spectroscopy, thereby improving the accuracy of physical properties inferred from nebular emission lines., Comment: 21 pages, 6 figures, 1 table, submitted to ApJ

We study, from both a theoretical and observational perspective, the physical origin and spectroscopic impact of extreme nebular emission in high-redshift galaxies. The nebular continuum, which can appear during extreme starbursts, is of particular importance as it tends to redden UV slopes and has a significant contribution to the UV luminosities of galaxies. Furthermore, its shape can be used to infer the gas density and temperature of the ISM. First, we provide a theoretical background, showing how different stellar populations (SPS models, IMFs, and stellar temperatures) and nebular conditions impact observed galaxy spectra. We demonstrate that, for systems with strong nebular continuum emission, 1) UV fluxes can increase by up to 0.7~magnitudes (or more in the case of hot/massive stars) above the stellar continuum, which may help reconcile the surprising abundance of bright high-redshift galaxies and the elevated UV luminosity density at $z>10$, 2) at high gas densities, UV slopes can redden from $\beta\lesssim-2.5$ to $\beta\sim-1$, 3) observational measurements of $\xi_{ion}$ are grossly underestimated, and 4) UV downturns from two-photon emission can masquerade as DLAs. Second, we present a dataset of 58 galaxies observed with NIRSpec on JWST at $2.510$ galaxies, finding that UV slopes and UV downturns are in some cases redder or steeper than expected from SPS models, which may hint at more exotic (e.g. hotter/more massive stars or AGN) ionizing sources., Comment: 26 pages, 18 figures, submitted to The Open Journal of Astrophysics

We present an analysis of JWST Ly$\alpha$ spectroscopy of $z\gtrsim6.5$ galaxies, using observations in the public archive covering galaxies in four independent fields (GOODS-N, GOODS-S, Abell 2744, EGS). We measure Ly$\alpha$ emission line properties for a sample of $210$ $z\simeq6.5-13$ galaxies, with redshifts confirmed independently of Ly$\alpha$ in all cases. We present $3$ new detections of Ly$\alpha$ emission in JWST spectra, including a large equivalent width (EW $=143\ \r{A}$) Ly$\alpha$ emitter with strong CIV emission (EW $=21\ \r{A}$) at $z=7.1$ in GOODS-N. We measure the redshift-dependent Ly$\alpha$ EW distribution across our sample. We find that strong Ly$\alpha$ emission (EW $>25\ \r{A}$) becomes increasingly rare at earlier epochs, suggesting that the transmission of Ly$\alpha$ photons decreases by $4\times$ between $z\simeq5$ and $z\simeq9$. We describe potential implications for the IGM neutral fraction. There is significant field to field variance in the Ly$\alpha$ emitter fraction. In contrast to the three other fields, the EGS shows no evidence for reduced transmission of Ly$\alpha$ photons at $z\simeq7-8$, suggesting a significantly ionized sightline may be present in the field. We use available NIRCam grism observations from the FRESCO survey to characterize overdensities on large scales around known Ly$\alpha$ emitters in the GOODS fields. The strongest overdensities appear linked with extremely strong Ly$\alpha$ detections (EW $>50\ \r{A}$) in most cases. Future Ly$\alpha$ spectroscopy with JWST has the potential to constrain the size of ionized regions around early galaxy overdensities, providing a new probe of the reionization process., Comment: 33 pages, 11 figures, accepted by ApJ

Strong gravitational lensing magnifies the light from a background source, allowing us to study these sources in detail. Here, we study the spectra of a $z = 1.95$ lensed Type Ia supernova SN~Encore for its brightest Image A, taken 39 days apart. We infer the spectral age with template matching using the supernova identification (SNID) software and find the spectra to be at 29.0 $\pm 5.0$ and 37.4 $\pm 2.8$ rest-frame days post maximum respectively, consistent with separation in the observer frame after accounting for time-dilation. Since SNe~Ia measure dark energy properties by providing relative distances between low- and high-$z$ SNe, it is important to test for evolution of spectroscopic properties. Comparing the spectra to composite low-$z$ SN~Ia spectra, we find strong evidence for similarity between the local sample of SN~Encore. The line velocities of common SN~Ia spectral lines, Si II 6355 and Ca II NIR triplet are consistent with the distribution for the low-$z$ sample as well as other lensed SNe~Ia, e.g. iPTF16geu ($z = 0.409$) and SN~H0pe ($z = 1.78$). The consistency in SN~Ia spectra across cosmic time demonstrates the utility of using SNe~Ia in the very high-$z$ universe for dark energy inference. We also find that the spectra of SN~Encore match the predictions for explosion models very well. With future large samples of lensed SNe~Ia, spectra at such late phases will be important to distinguish between different explosion scenarios., Comment: 9 pages, 8 figures, submitted to MNRAS. Comments welcome

While JWST has observed galaxies assembling as early as $z\sim14$, evidence of galaxies with significant old stellar populations in the Epoch of Reionisation (EoR) -- the descendants of these earliest galaxies -- are few and far between. Bursty star-formation histories (SFHs) have been invoked to explain the detectability of the earliest UV-bright galaxies, but also to interpret galaxies showing Balmer breaks without nebular emission lines. We present the first spectroscopic evidence of a $z\sim7.9$ galaxy, A2744-YD4, which shows a Balmer break and emission lines, indicating the presence of both a mature and young stellar population. The spectrum of A2744-YD4 shows peculiar emission line ratios suggesting a relatively low ionisation parameter and high gas-phase metallicity. A median stack of galaxies with similar emission line ratios reveals a clear Balmer break in their stacked spectrum. This suggests that a mature stellar population ($\sim 80$ Myr old) has produced a chemically enriched, disrupted interstellar medium. Based on SED-fitting and comparison to simulations, we conclude that the observed young stellar population is in fact the result of a rejuvenation event following a lull in star formation lasting $\sim 20$ Myr, making A2744-YD4 and our stack the first spectroscopic confirmation of galaxies that have rejuvenated following a mini-quenched phase. These rejuvenating galaxies appear to be in an exceptional evolutionary moment where they can be identified. Our analysis shows that a young stellar population of just $\sim 30 \%$ of the total stellar mass would erase the Balmer break. Hence, 'outshining' through bursty SFHs of galaxies in the early Universe is likely plaguing attempts to measure their stellar ages and masses accurately., Comment: 14 pages, 12 figures, 1 table. Submitted to MNRAS. Comments are welcome

We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM., Comment: 17 pages, 8 figures, submitted to ApJ

We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$\lambda$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance., Comment: 32 pages, 22 figures, submitted to ApJ

Recent stellar chemical abundance measurements of a handful of $z\sim2$ quiescent galaxies have suggested these galaxies exhibit a remarkably strong $\alpha$-enhancement compared to their local and intermediate redshift counterparts. This apparent chemical evolution following quenching suggests that even the innermost regions of massive early-type galaxies may have experienced substantial mixing of stars in mergers, challenging a purely inside-out growth model. However, larger samples are needed to determine whether a high $\alpha$-enhancement ([Mg/Fe] $\approx 0.5$) is common in $z \sim 2$ quiescent galaxies, and a comparative analysis is needed to determine whether it is consistently inferred using different stellar population synthesis models. We report age and stellar chemical abundance measurements for a sample of four gravitationally lensed quiescent galaxies at $z\sim2.1-2.65$ based on Magellan/FIRE spectroscopy. For three of these galaxies we constrain the $\alpha$-enhancement, and in two cases we measure high values comparable to earlier results when the spectra are analyzed consistently. We also find that the choice of modeling approach can exert a significant effect on the measured abundances. This model dependence can be partly, but not entirely, explained by the complex abundance patterns of $\alpha$ elements in galaxies, which has been observed at lower redshifts and in one $z \sim 2$ quiescent galaxy. Our investigation highlights the importance of independently varying abundance of $\alpha$ elements when fitting the spectra of such galaxies. Observations with JWST will soon deliver precise and spatially resolved abundances of these and other quiescent galaxies at cosmic noon, opening a new window into their evolution., Comment: 10 pages, submitted to The Astrophysical Journal

JWST observations have recently begun delivering the first samples of Ly$\alpha$ velocity profile measurements at $z>6$, opening a new window on the reionization process. Interpretation of $z\gtrsim6$ line profiles is currently stunted by limitations in our knowledge of the intrinsic Ly$\alpha$ profile (before encountering the IGM) of the galaxies that are common at $z\gtrsim6$. To overcome this shortcoming, we have obtained resolved ($R\sim3900$) Ly$\alpha$ spectroscopy of $42$ galaxies at $z=2.1-3.4$ with similar properties as are seen at $z>6$. We quantify a variety of Ly$\alpha$ profile statistics as a function of [OIII]+H$\beta$ EW. Our spectra reveal a new population of $z\simeq 2-3$ galaxies with large [OIII]+H$\beta$ EWs ($>1200\ \r{A}$) and a large fraction of Ly$\alpha$ flux emerging near the systemic redshift (peak velocity $\simeq0$ km s$^{-1}$). These spectra indicate that low density neutral hydrogen channels are able to form in a subset of low mass galaxies ($\lesssim1\times10^8\ M_{\odot}$) that experience a burst of star formation (sSFR $>100$ Gyr$^{-1}$). Other extreme [OIII] emitters show weaker Ly$\alpha$ that is shifted to higher velocities ($\simeq240$ km s$^{-1}$) with little emission near line center. We investigate the impact the IGM is likely to have on these intrinsic line profiles in the reionization era, finding that the centrally peaked Ly$\alpha$ emitters should be strongly attenuated at $z\gtrsim5$. We show that these line profiles are particularly sensitive to the impact of resonant scattering from infalling IGM and can be strongly attenuated even when the IGM is highly ionized at $z\simeq 5$. We compare these expectations against a new database of $z\gtrsim6.5$ galaxies with robust velocity profiles measured with JWST/NIRSpec., Comment: 28 pages, 11 figures, accepted by ApJ

We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN., Comment: 23 pages, 17 figures, submitted to ApJ

We report first results from a deep near infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra-Deep Field (HUDF), 10-15 years after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability, and measure (or constrain) the comoving number density of supermassive black holes (SMBHs), n_{SMBH}, at early times. In this Letter we present an overview of the program and preliminary results concerning eight objects. Three variables are supernovae, two of which are apparently hostless with indeterminable redshifts, although one has previously been recorded at a z\approx 6 object precisely because of its transient nature. Two further objects are clear AGN at z= 2.0 and 3.2, based on morphology and/or infrared spectroscopy from JWST. Three variable targets are identified at z = 6-7, which are also likely AGN candidates. These sources provide a first measure of n_{SMBH} in the reionization epoch by photometric variability, which places a firm lower limit of 3 \times 10^{-4} cMpc^{-3}. After accounting for variability and luminosity incompleteness, we estimate n_{SMBH} \gtrsim 8 \times 10{-3} cMpc{-3}, which is the largest value so far reported at these redshifts. This SMBH abundance is also strikingly similar to estimates of n_{SMBH} in the local Universe. We discuss how these results test various theories for SMBH formation., Comment: In press at ApJL

Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim500\,{\rm arcmin}^2$, we obtain a sample of $>$30,000 galaxies at $z_{\rm phot}\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a new population of UV-red galaxies that was missing from previous HST-based Lyman Break Galaxy (LBG) samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ down to limiting masses of $10^{7.5}-10^{8.5}\,{\rm M_\odot}$, finding steep low mass slopes over the entire redshift range, reaching values of $\alpha\approx-2$ at $z\gtrsim6$. At the high-mass end, UV-red galaxies dominate at least out to $z\sim6$. The implied redshift evolution of the SMF suggests a rapid build-up of massive dust-obscured or quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history., Comment: 31 pages, 17 figures, published in MNRAS

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

JWST has recently sparked a new era of Ly$\alpha$ spectroscopy, delivering the first measurements of the Ly$\alpha$ escape fraction and velocity profile in typical galaxies at $z\simeq6-10$. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Ly$\alpha$ properties in galaxies at $z\simeq5-6$ when the IGM is mostly ionized. Here we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Ly$\alpha$ velocity offsets, escape fractions, and EWs in $z\simeq5-6$ galaxies. We find that galaxies with large Ly$\alpha$ escape fractions (> 0.2) are common at $z\simeq5-6$, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Ly$\alpha$ becomes more prevalent in the galaxy population toward higher redshift from $z\sim3$ to $z\sim6$, although we find that this evolution slows considerably between $z\sim5$ and $z\sim6$, consistent with modest attenuation from residual HI in the mostly ionized IGM at $z\simeq5-6$. We find significant evolution in Ly$\alpha$ velocity profiles between $z\simeq2-3$ and $z\simeq5-6$, likely reflecting the influence of resonant scattering from residual intergalactic HI on the escape of Ly$\alpha$ emission near line center. This effect will make it challenging to use Ly$\alpha$ peak offsets as a probe of Lyman continuum leakage at $z\simeq5-6$. We use our $z\simeq5-6$ Ly$\alpha$ distributions to make predictions for typical Ly$\alpha$ properties at $z\gtrsim8$ and discuss implications of a recently-discovered Ly$\alpha$ emitter at $z\simeq8.5$ with a small peak velocity offset (156 km s$^{-1}$)., Comment: 31 pages, 12 figures, accepted by MNRAS

We study the spatially resolved outflow properties of CSWA13, an intermediate mass ($M_*=10^{9}~\mathrm{M}_{\odot}$), gravitationally lensed star-forming galaxy at $z=1.87$. We use Keck/KCWI to map outflows in multiple rest-frame ultraviolet ISM absorption lines, along with fluorescent Si II$^*$ emission, and nebular emission from C III] tracing the local systemic velocity. The spatial structure of outflow velocity mirrors that of the nebular kinematics, which we interpret to be a signature of a young galactic wind that is pressurizing the ISM of the galaxy but is yet to burst out. From the radial extent of Si II$^*$ emission, we estimate that the outflow is largely encapsulated within $3.5$ kpc. We explore the geometry (e.g., patchiness) of the outflow by measuring the covering fraction at different velocities, finding that the maximum covering fraction is at velocities $v\simeq-150$ km$\,$s$^{-1}$. Using the outflow velocity ($v_{out}$), radius ($R$), column density ($N$), and solid angle ($\Omega$) based on the covering fraction, we measure the mass loss rate $\log\dot{m}_{out}/(\mathrm{M}_{\odot}\text{yr}^{-1}) = 1.73\pm0.23$ and mass loading factor $\log\eta = 0.04\pm0.34$ for the low-ionization outflowing gas in this galaxy. These values are relatively large and the bulk of the outflowing gas is moving with speeds less than the escape velocity of the galaxy halo, suggesting that the majority of outflowing mass will remain in the circumgalactic medium and/or recycle back into the galaxy. The results support a picture of high outflow rates transporting mass and metals into the inner circumgalactic medium, providing the gas reservoir for future star formation., Comment: Submitted to ApJ. Comments welcome (17 pages, 10 figures, 1 table)

The velocity offset of Ly$\alpha$ emission from a galaxy's systemic redshift is an excellent tracer of conditions that enable the escape of Ly$\alpha$ photons from the galaxy, and potentially hydrogen ionizing Lyman continuum photons. However at $z\geq6$, Ly$\alpha$ is often heavily attenuated by the neutral intergalactic medium. Here we investigate the utility of Ciii],$\lambda\lambda1907,1909$ emission, usually the brightest UV line after Ly$\alpha$, as a proxy estimating the Ly$\alpha$ velocity offset ($\Delta v_{\rm{Ly}\alpha}$). To do so, we use analogues of reionization era galaxies based upon 52 star-forming galaxies with robust Ciii] detections drawn from the VANDELS survey. Our sample spans a broad UV magnitude range of $-18.5 < M_{\rm{UV}} < -22.0$, with a sample average value of EW(Ciii]) $=5.3$\,{\AA}. We find a slight increase of EW(Ciii]) with increasing EW(Ly$\alpha$), but find a large range of EW(Ciii]) $\sim$1-13\,{\AA} particularly at EW(Ly$\alpha$) $<10$\,{\AA}. Using the Ciii] line peak as the systemic redshift, we calculate $\Delta v_{\rm{Ly}\alpha}$ and recover the previously reported trend of decreasing $\Delta v_{\rm{Ly}\alpha}$ with increasing EW(Ly$\alpha$). Interestingly, we find an anti-correlation between $\Delta v_{\rm{Ly}\alpha}$ and EW(Ciii]), which also displays a dependence on the UV absolute magnitude. We derive a multi-variate fit to obtain $\Delta v_{\rm{Ly}\alpha}$ using both EW(Ciii]) and $M_{\rm{UV}}$, finding that $\Delta v_{\rm{Ly}\alpha}$ is more strongly dependent on EW(Ciii]), with a weaker but non-negligible dependence on $M_{\rm{UV}}$. We find that for a fixed EW(Ciii]), UV-bright Ly$\alpha$ emitting galaxies show smaller values of $\Delta v_{\rm{Ly}\alpha}$,which suggests that such galaxies may be undergoing more bursty star-formation compared to the UV-fainter ones, akin to a population of extremely UV-bright galaxies identified at $z>10$., Comment: 11 pages, 7 figures

During the epoch of reionisation the first galaxies were enshrouded in pristine neutral gas, with one of the brightest emission lines in star-forming galaxies, Lyman-$\alpha$ (Ly$\alpha$), expected to remain undetected until the Universe became ionised. Providing an explanation for the surprising detection of Ly$\alpha$ in these early galaxies is a major challenge for extra-galactic studies. Recent JWST observations have reignited the debate on whether residence in an overdensity of galaxies is a it sufficient and necessary condition for Ly$\alpha$ to escape. Here, we take unique advantage of both high-resolution and high-sensitivity images from the JWST instrument NIRCam to reveal that all galaxies in a sample of z>7 Ly$\alpha$ emitters have close companions. We exploit novel on-the-fly radiative transfer magnetohydrodynamical simulations with cosmic ray feedback to show that galaxies with frequent mergers have very bursty star formation which drives episodes of high intrinsic Ly$\alpha$ emission and facilitates the escape of Ly$\alpha$ photons along channels cleared of neutral gas. We conclude that the rapid build up of stellar mass through mergers presents a compelling solution to the long-standing puzzle of the detection of Ly$\alpha$ emission deep into the epoch of reionisation., Comment: Published in Nature Astronomy. 45 pages, 9 figures, 2 tables

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

We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distributions exploiting the extensive multi-wavelength data for the HDF-N. We identify the 7 highest-fidelity sources as CO emitters at 1=5.0. The relatively large area and depth of the NOEMA HDF-N survey extends the existing luminosity functions at 1

We present a novel technique called "photometric IGM tomography" to map the intergalactic medium (IGM) at $z\simeq4.9$ in the COSMOS field. It utilizes deep narrow-band (NB) imaging to photometrically detect faint Ly$\alpha$ forest transmission in background galaxies across the Subaru/Hyper-Suprime Cam (HSC)'s $1.8\rm\,sq.\,deg$ field of view and locate Ly$\alpha$ emitters (LAEs) in the same cosmic volume. Using ultra-deep HSC images and Bayesian spectral energy distribution fitting, we measure the Ly$\alpha$ forest transmission at $z\simeq4.9$ along a large number ($140$) of background galaxies selected from the DEIMOS10k spectroscopic catalogue at $4.98

We use the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation to study how Lyman Continuum (LyC) photons escape from galaxies and the observational signatures of this escape. We define two classes of LyC leaker: Bursty Leakers and Remnant Leakers, based on their star formation rates (SFRs) that are averaged over 10 Myr (SFR$_{10}$) or 100 Myr (SFR$_{100}$). Both have $f_{\rm esc}>20\%$ and experienced an extreme burst of star formation, but Bursty Leakers have ${\rm SFR_{10}>SFR_{100}}$, while Remnant Leakers have ${\rm SFR_{10}20\%$ dominate the ionizing photon budget at $z\gtrsim7.5$ but the contribution from galaxies with $f_{\rm esc}<5\%$ becomes significant at the tail-end of reionization., Comment: 16 pages, 14 figures, accepted by MNRAS

We study the kinematics of the interstellar medium (ISM) viewed "down the barrel" in 20 gravitationally lensed galaxies during Cosmic Noon ($z=1.5 - 3.5$). We use moderate-resolution spectra ($R\sim4000$) from Keck/ESI and Magellan/MagE to spectrally resolve the ISM absorption in these galaxies into $\sim$10 independent elements and use double Gaussian fits to quantify the velocity structure of the gas. We find that the bulk motion of gas in this galaxy sample is outflowing, with average velocity centroid $\left=-141$ km$\,$s$^{-1}$ ($\pm111$ km$\,$s$^{-1}$ scatter) measured with respect to the systemic redshift. 16 out of the 20 galaxies exhibit a clear positive skewness, with a blueshifted tail extending to $\sim -500$ km$\,$s$^{-1}$. We examine scaling relations in outflow velocities with galaxy stellar mass and star formation rate (SFR), finding correlations consistent with a momentum-driven wind scenario. Our measured outflow velocities are also comparable to those reported for FIRE-2 and TNG50 cosmological simulations at similar redshift and galaxy properties. We also consider implications for interpreting results from lower-resolution spectra. We demonstrate that while velocity centroids are accurately recovered, the skewness, velocity width, and probes of high velocity gas (e.g., $v_{95}$) are subject to large scatter and biases at lower resolution. We find that $R\gtrsim1700$ is required for accurate results for the gas kinematics of our sample. This work represents the largest available sample of well-resolved outflow velocity structure at $z>2$, and highlights the need for good spectral resolution to recover accurate properties., Comment: 42 pages, 37 figures (including appendix), Accepted for publication, ApJ

The recent launch of JWST has enabled the exciting prospect of detecting the first generation of metal-free, Population III (Pop. III) stars. Determining the emission line signatures that robustly signify the presence of Pop. III stars against other possible contaminants represents a key challenge for interpreting JWST data. To this end, we run high-resolution (sub-pc) cosmological radiation hydrodynamics simulations of the region around a dwarf galaxy at $z\geq10$ to predict the emission line signatures of the Pop. III/Pop. II transition. We show that the absence of metal emission lines is a poor diagnostic of Pop. III stars because metal-enriched galaxies in our simulation can maintain low [OIII] 5007${\rm \r{A}}$ emission that may be undetectable due to sensitivity limits. Combining spectral hardness probes (e.g. HeII 1640${\rm \r{A}}$/H$\alpha$) with metallicity diagnostics is more likely to probe the presence of metal-free stars, although contamination from Wolf-Rayet stars, X-ray binaries, or black holes may be important. The hard emission from Pop. III galaxies fades fast due to the short stellar lifetimes of massive Pop. III stars, which could further inhibit detection. Similarly, Pop. III stars may be detectable after they evolve off the main-sequence due to the cooling radiation from nebular gas or a supernova remnant; however, these signatures are also short-lived (i.e. few Myr), and contaminants such as flickering black holes might confuse this diagnostic. While JWST will provide a unique opportunity to spectroscopically probe the nature of the earliest galaxies, both the short timescales associated with pristine systems and ambiguities in interpreting key diagnostic emission lines may hinder progress. Special care will be needed before claiming the discovery of systems with pure Pop. III stars., Comment: 9 pages, 8 figures, submitted to MNRAS

The 100% detection rate of Ly$\alpha$ emission in a sample of four luminous $z\sim8$ galaxies with red Spitzer/IRAC colors suggests objects with unusual ionizing capabilities that created early ionized bubbles in a neutral era. Whether such bubbles reflect enhanced ionizing properties (nature) or an overdense environment (nurture), however, remains unclear. Here we aim to distinguish between these hypotheses via a search for Ly$\alpha$ emission in five fainter galaxies drawn from the CANDELS-GOODS fields using a similar IRAC excess and UV magnitudes that should reflect reduced clustering effects. Using Keck/MOSFIRE we tentatively detect $>4\sigma$ line emission in only two targets at redshifts $z_{\rm Ly\alpha}$=7.1081 and $z_{\rm Ly\alpha}$=7.9622 with rest-frame EWs of 16-17 A, $\sim$1.5$\times$ weaker compared to their brighter counterparts. Thus we find a reduced rate for Ly$\alpha$ emission of $0.40^{+0.30}_{-0.25}$ compared to $1.00^{+0.00}_{-0.44}$ for more luminous examples. The lower rate agrees with predictions from simulations of a mostly neutral IGM and an intrinsic EW$_{0,\rm Ly\alpha}$ distribution for $z\sim6$ galaxies. However, even with an extreme EW$_{0,\rm Ly\alpha}$ model, it is challenging to match the detection rate for the luminous objects. SED-fitting of our fainter sample indicates young and star-forming systems, albeit with less extreme SFRs and ionization parameters compared to their luminous counterparts. The enhanced Ly$\alpha$ rate in luminous galaxies is thus likely a byproduct of both extreme ionizing properties as well as environment effects. Further studies with JWST may be required to resolve the physical nature of this puzzling population., Comment: 15 pages, 4 figures, 2 tables. Accepted for publication in ApJ (in press), uploaded to match accepted version

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$\mu$m emission is detected at 10$\sigma$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $\Delta V_{\rm obs}/2\sigma_{\rm tot}=0.84\pm0.23$, where $\Delta V_{\rm obs}$ is the observed maximum velocity difference and $\sigma_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/\sigma_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $\sigma_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$., Comment: ApJL accepted

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 $z<10$, consistent with a claimed rapid rise in the star formation density at that time. Here we study a spectroscopically-confirmed sample of extreme [OIII] emitters at $z=1.3-3.7$, using both dynamical masses estimated from [OIII] line widths and rest-frame UV to near-infrared photometry to illustrate the dangers of assuming such systems are genuinely young. For the most extreme of our intermediate redshift line emitters, we find dynamical masses $10-100$ times that associated with a young stellar population mass, which are difficult to explain solely by the presence of additional dark matter or gaseous reservoirs. Adopting nonparametric star formation histories, we show how the near-infrared photometry of a subset of our sample reveals an underlying old ($>100$ 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$., Comment: 14 pages, 7 figure, accepted by MNRAS

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<-2$), reflecting minimal dust obscuration, young ages, and low metallicities. This contrasts with narrow-lined AGN at $z\sim2-3$, whose UV slopes are significantly redder ($\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 report the discovery of a double-peaked Lyman-alpha profile in a galaxy at $z=6.803$, A370p\_z1, in the parallel Frontier Field of Abell 370. The velocity separation between the blue and red peaks of the Lyman-$\alpha$ profile ($\Delta v=101_{-19}^{+38} (\pm48)\,\text{km s}^{-1}$) suggests an extremely high escape fraction of ionising photons $> 59(51)\% (2\sigma)$. The spectral energy distribution indicates a young ($50$ Myr), star-forming ($12\pm 6 \, \rm{M}_\odot \rm{yr}^{-1}$) galaxy with an IRAC excess implying strong [OIII]+H$\beta$ emission. On the basis of the high escape fraction measured, we demonstrate that A370p\_z1 was solely capable of creating an ionised bubble sufficiently large to account for the blue component of its Lyman-alpha profile. We discuss whether A370p\_z1 may be representative of a larger population of luminous $z\simeq7$ double-peaked Lyman-alpha emitting sources with high escape fractions that self-ionised their surroundings without contributions from associated UV-fainter sources., Comment: Accepted for publication in MNRAS. 7 pages, 4 figures, 3 tables

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

Measuring the chemical composition of galaxies is crucial to our understanding of galaxy formation and evolution models. However, such measurements are extremely challenging for quiescent galaxies at high redshifts, which have faint stellar continua and compact sizes, making it difficult to detect absorption lines and nearly impossible to spatially resolve them. Gravitational lensing offers the opportunity to study these galaxies with detailed spectroscopy that can be spatially resolved. In this work, we analyze deep spectra of MRG-M0138, a lensed quiescent galaxy at z = 1.98 which is the brightest of its kind, with an H-band magnitude of 17.1. Taking advantage of full spectral fitting, we measure $[{\rm Mg/Fe}]=0.51\pm0.05$, $[\rm{Fe/H}]=0.26\pm0.04$, and, for the first time, the stellar abundances of 6 other elements in this galaxy. We further constrained, also for the first time in a $z\sim2$ galaxy, radial gradients in stellar age, [Fe/H], and [Mg/Fe]. We detect no gradient in age or [Mg/Fe] and a slightly negative gradient in [Fe/H], which has a slope comparable to that seen in local early-type galaxies. Our measurements show that not only is MRG-M0138 very Mg-enhanced compared to the centers of local massive early-type galaxies, it is also very iron rich. These dissimilar abundances suggest that even the inner regions of massive galaxies have experienced significant mixing of stars in mergers, in contrast to a purely inside-out growth model. The abundance pattern observed in MRG-M0138 challenges simple galactic chemical evolution models that vary only the star formation timescale and shows the need for more elaborate models., Comment: Accepted for publication in ApJ Letters

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions., Comment: 21 pages, 13 figures, updated to match MNRAS version

Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$\mu$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$\mu$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$\mu$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch., Comment: 12 pages, 3 figures, 4 tables. Accepted for publication 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 examine the absorption line spectra of a sample of 31 luminous (M_UV=-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our lower resolution data to measure the depth of various interstellar and stellar absorption lines to probe the covering fraction of low ionization gas and the gas-phase and stellar metallicities near the end of the era of cosmic reionization. From maximum absorption line depths of SiII1260 and CII1334, we infer a mean covering fraction of >0.85+/-0.16 for our sample. This is larger than that determined using similar methods for lower luminosity galaxies at slightly lower redshifts, suggesting that the most luminous galaxies appear to have a lower escape fraction than fainter galaxies, and therefore may not play a prominent role in concluding reionization. Using various interstellar absorption lines we deduce gas-phase metallicities close to solar indicative of substantial early enrichment. Using selected stellar absorption lines, we model our spectra with a range of metallicities using techniques successfully employed at lower redshift and deduce a stellar metallicity of 0.4 +0.3/-0.1 solar, consistent with the stellar mass - stellar metallicity relation recently found at z~3-5. We discuss the implications of these metallicity estimates for the typical ages of our luminous galaxies and conclude our results imply initial star formation at redshifts z~10, consistent with independent analyses of earlier objects., Comment: 13 pages, 10 figures, 3 tables, Accepted for Publication in ApJ

Prior to the launch of JWST, Spitzer/IRAC photometry offers the only means of studying the rest-frame optical properties of z>7 galaxies. Many such high redshift galaxies display a red [3.6] - [4.5] micron colour, often referred to as the "IRAC excess", which has conventionally been interpreted as arising from intense [OIII]+Hbeta emission within the [4.5] micron bandpass. An appealing aspect of this interpretation is similarly intense line emission seen in star-forming galaxies at lower redshift as well as the redshift-dependent behaviour of the IRAC colours beyond z~7 modelled as the various nebular lines move through the two bandpasses. In this paper we demonstrate that, given the photometric uncertainties, established stellar populations with Balmer (4000 A, rest-frame) breaks, such as those inferred at z>9 where line emission does not contaminate the IRAC bands, can equally well explain the redshift-dependent behaviour of the IRAC colours in 7

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$\alpha$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$\alpha$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium (IGM) transmission at the end of reionisation. We have gathered a sample of $13$ spectroscopically confirmed Lyman-break galaxies (LBGs) and $21$ Lyman-$\alpha$ emitters (LAEs) at angular separations $20'' \lesssim \theta \lesssim 10'$ ($\sim 0.1-4$ pMpc at $z\sim 6$) from the sightlines to $8$ background $z\gtrsim 6$ quasars. We report for the first time the detection of an excess of Lyman-$\alpha$ transmission spikes at $\sim 10-60$ cMpc from LAEs ($3.2\sigma$) and LBGs ($1.9\sigma$). We interpret the data with an improved model of the galaxy-Lyman-$\alpha$ transmission and two-point cross-correlations which includes the enhanced photoionisation due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)-Lyman-$\alpha$ transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionisation to $\langle f_{\rm esc} \rangle_{M_{\rm UV}<-12} = 0.14_{-0.05}^{+0.28}\,(0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionisation. Our results show that a contribution from faint galaxies ($M_{\rm UV} > -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionisation might be driven by different sub-populations around LBGs and LAEs at $z\sim 6$., Comment: Accepted for publication in MNRAS. 17 pages, 18 figures, 3 tables + appendices. Key results: Figs. 6,13,14. Updated version matches the published one ('online' Appendices are directly in the arXiv pdf). Conclusions unchanged from earlier preprint

Dark matter-only simulations predict that dark matter halos have steep, cuspy inner density profiles, while observations of dwarf galaxies find a range of inner slopes that are often much shallower. There is debate whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. In Paper 1 of this series, we obtained high-resolution integral field H$\alpha$ observations for 26 dwarf galaxies with $M_*=10^{8.1}-10^{9.7}\textrm{M}_\odot$. We derived rotation curves from our observations, which we use here to construct mass models. We model the total mass distribution as the sum of a generalized Navarro-Frenk-White (NFW) dark matter halo and the stellar and gaseous components. Our analysis of the slope of the dark matter density profile focuses on the inner 300-800 pc, chosen based on the resolution of our data and the region resolved by modern hydrodynamical simulations. The inner slope measured using ionized and molecular gas tracers is consistent, and it is additionally robust to the choice of stellar mass-to-light ratio. We find a range of dark matter profiles, including both cored and cuspy slopes, with an average of $\rho_{\rm DM}\sim r^{-0.74\pm 0.07}$, shallower than the NFW profile, but steeper than those typically observed for lower-mass galaxies with $M_*\sim 10^{7.5}\textrm{M}_\odot$. Simulations that reproduce the observed slopes in those lower-mass galaxies also produce slopes that are too shallow for galaxies in our mass range. We therefore conclude that supernova feedback models do not yet provide a fully satisfactory explanation for the observed trend in dark matter slopes.

We present the stellar mass-[Fe/H] and mass-[Mg/H] relation of quiescent galaxies in two galaxy clusters at $z\sim0.39$ and $z\sim0.54$. We derive the age, [Fe/H], and [Mg/Fe] for each individual galaxy using a full-spectrum fitting technique. By comparing with the relations for $z\sim0$ SDSS galaxies, we confirm our previous finding that the mass-[Fe/H] relation evolves with redshift. The mass-[Fe/H] relation at higher redshift has lower normalization and possibly steeper slope. However, based on our sample, the mass-[Mg/H] relation does not evolve over the observed redshift range. We use a simple analytic chemical evolution model to constrain average outflow that these galaxies experience over their lifetime, via the calculation of mass-loading factor. We find that the average mass-loading factor $\eta$ is a power-law function of galaxy stellar mass, $\eta \propto M_*^{-0.21\pm0.09}$. The measured mass-loading factors are consistent with the results of other observational methods for outflow measurements and with the predictions where outflow is caused by star formation feedback in turbulent disks., Comment: 22 pages, 11 figures, Accepted for publication in ApJ

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., Comment: 15 pages, 6 figures. Accepted for publication in ApJ

We present the results of a model-independent investigation of the rest-frame UV spectra from a comprehensive sample of $394$ quasars in the redshift range $1.5\leq z \leq 7.5$. We fit the main Broad Emission Lines (BELs) in the rest-frame range $1280 \text{\AA} \leq \lambda \leq 3000 \text{\AA}$ (O I, C II, Si IV, C III, C IV and Mg II) with a lightly-supervised spline fitting technique. Redshifts are derived from the peaks of each fitted BEL and used to compute relative velocity shifts. We show that our method gives unbiased velocity shifts and is insensitive to spectral resolution and instrumental parameters. It is found that the average blueshift of the \cfour\, line with respect to several low-ionisation lines in luminosity-matched samples does not significantly evolve over $1.5\gtrsim z\gtrsim6$. However, the average blueshift increases significantly by a factor $\sim 2.5$ at $z\gtrsim 6$. We propose that this redshift evolution can be explained by \cfour\, winds launched perpendicularly to an accretion disk with increased torus opacity at high-redshift, coupled with a potential orientation-driven selection bias. Our results open new exciting avenues of investigation into young quasars in the reionisation epoch., Comment: Accepted for publication in MNRAS. 14 pages, 14 Figures, 2 Tables. Revised version to match the accepted one - results unchanged. Supplementary material for appendix A available at : www.homepages.ucl.ac.uk/~ucapram/QSOBEL_shifts_supp_material.tar.gz

We report the discovery of the first spectroscopically resolved C II /C II* 1334, 1335A doublet in the Lyman-break galaxy J0215-0555 at z = 5.754. The separation of the resonant and fluorescent emission channels was possible thanks to the large redshift of the source and long integration time, as well as the small velocity width of the feature, 0.6 +- 0.2A. We model this emission and find that at least two components are required to reproduce the combination of morphologies of C II* emission, C II absorption and emission, and Lyman-alpha emission from the object. We suggest that the close alignment between the fluorescence and Lyman-alpha emission could indicate an ionisation escape channel within the object. While the faintness of such a C II /C II* doublet makes it prohibitively difficult to pursue for similar systems with current facilities, we suggest it can become a valuable porosity diagnostic in the era of JWST and the upcoming generations of ELTs., Comment: 5 pages, 3 figures; accepted for publication in MNRAS Letters

We investigate the stellar populations for a sample of 24 quiescent galaxies at 1.5 < z < 2.5 using deep rest-frame optical spectra obtained with Keck MOSFIRE. By fitting templates simultaneously to the spectroscopic and photometric data, and exploring a variety of star formation histories, we obtain robust measurements of median stellar ages and residual levels of star formation. After subtracting the stellar templates, the stacked spectrum reveals the Halpha and [NII] emission lines, providing an upper limit on the ongoing star formation rate of 0.9 +/- 0.1 Msun/yr. By combining the MOSFIRE data to our sample of Keck LRIS spectra at lower redshift, we analyze in a consistent manner the quiescent population at 1 < z < 2.5. We find a tight relation (with a scatter of 0.13 dex) between the stellar age and the rest-frame U-V and V-J colors, which can be used to estimate the age of quiescent galaxies given their colors. Applying this age--color relation to large, photometric samples, we are able to model the number density evolution for quiescent galaxies of various ages. We find evidence for two distinct quenching paths: a fast quenching that produces compact post-starburst systems, and a slow quenching of larger galaxies. Fast quenching accounts for about a fifth of the growth of the red sequence at z~1.4, and half at z~2.2. We conclude that fast quenching is triggered by dramatic events such as gas-rich mergers, while slow quenching is likely caused by a different physical mechanism., Comment: 28 pages, 15 figures, accepted in ApJ

The presence of spectroscopically confirmed Balmer breaks in galaxy spectral energy distributions (SEDs) at $z>9$ provides one of the best probes of the assembly history of the first generations of stars in our Universe. Recent observations of the gravitationally lensed source, MACS 1149_JD1 (JD1), indicate that significant amounts of star formation likely occurred at redshifts as high as $z\simeq15$. The inferred stellar mass, dust mass, and assembly history of JD1, or any other galaxy at these redshifts that exhibits a strong Balmer break, can provide a strong test of our best theoretical models from high-resolution cosmological simulations. In this work, we present the results from a cosmological radiation-hydrodynamics simulation of the region surrounding a massive Lyman-break galaxy. For two of our most massive systems, we show that dust preferentially resides in the vicinity of the young stars thereby increasing the strength of the measured Balmer break such that the simulated SEDs are consistent with the photometry of JD1 and two other $z>9$ systems (GN-z10-3 and GN-z9-1) that have proposed Balmer breaks at high redshift. We find strong variations in the shape and luminosity of the SEDs of galaxies with nearly identical stellar and halo masses, indicating the importance of morphology, assembly history, and dust distribution in making inferences on the properties of individual galaxies at high redshifts. Our results stress the importance that dust may play in modulating the observable properties of galaxies, even at the extreme redshifts of $z>9$., Comment: 16 pages, 13 Figures, Accepted to MNRAS