Search

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

The baryon cycle is crucial for understanding galaxy formation, as gas inflows and outflows vary throughout a galaxy's lifetime and affect its star formation rate. Despite the necessity of accretion for galaxy growth at high redshifts, direct observations of inflowing gas have proven elusive especially at $z\gtrsim2$. We present spectroscopic analysis of a galaxy at redshift $z=2.45$ which exhibits signs of inflow in several ultraviolet interstellar absorption lines, with no clear outflow signatures. The absorption lines are redshifted by $\sim$250 km sec$^{-1}$ with respect to the systemic redshift, and C IV shows a prominent inverse P-Cygni profile. Simple stellar population models suggest that this galaxy has a low metallicity ($\sim$5\% solar), with a very young starburst of age $\sim$4 Myr dominating the ultraviolet luminosity. The gas inflow velocity and nebular velocity dispersion suggest an approximate halo mass of order $\sim 10^{11}M_{\odot}$, a regime in which simulations predict that bursty star formation is common at this redshift. We conclude that this system is likely in the beginning of a cycle of bursty star formation, where inflow and star formation rates are high, but where supernovae and other feedback processes have not yet launched strong outflows. In this scenario, we expect the inflow-dominated phase to be observable (e.g., with net redshifted ISM absorption) for only a short timescale after a starburst onset. This result represents a promising avenue for probing the full baryon cycle, including inflows, during the formative phases of low-mass galaxies at high redshifts., Comment: 12 pages, 4 figures

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

In this paper, we describe the wide-field spectroscopic survey telescope (WST) project. WST is a 12-metre wide-field spectroscopic survey telescope with simultaneous operation of a large field-of-view (3 sq. degree), high-multiplex (20,000) multi-object spectrograph (MOS), with both a low and high-resolution modes, and a giant 3x3 arcmin2 integral field spectrograph (IFS). In scientific capability, these specifications place WST far ahead of existing and planned facilities. In only 5 years of operation, the MOS would target 250 million galaxies and 25 million stars at low spectral resolution, plus 2 million stars at high resolution. Without need for pre-imaged targets, the IFS would deliver 4 billion spectra offering many serendipitous discoveries. 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 in synergy with future ground and space-based facilities. We show how it can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; the origin of stars and planets; and time domain and multi-messenger astrophysics. WST's uniquely rich dataset may yield unforeseen discoveries in many of these areas. The telescope and instruments are designed as an integrated system and will mostly use existing technology, with the aim to minimise the carbon footprint and environmental impact. We will propose WST as the next European Southern Observatory (ESO) project after completion of the 39-metre ELT., Comment: 15 pages, 15 figures

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