Search

We present JWST/NIRSpec PRISM spectroscopic characterization of GHZ9 at z= 10.145 $\pm$ 0.010, currently the most distant source detected by the Chandra X-ray Observatory. The spectrum reveals several UV high-ionization lines, including CII, SiIV, [NIV], CIV, HeII, OIII], NIII], and CIII]. The prominent rest-frame equivalent widths (EW(CIV)$\simeq$65A, EW(HeII)$\simeq$18A, EW(CIII])$\simeq$48A) show the presence of a hard radiation field, while the analysis of line ratio diagnostics suggest this galaxy hosts both AGN and star-formation activity. GHZ9 is nitrogen-enriched (6--9.5 times solar), carbon-poor (0.2--0.65 times solar), metal-poor (Z = 0.01--0.1 Z$_{\odot}$), and compact ($<$ 106 pc), similarly to GNz11, GHZ2, and recently discovered N-enhanced high redshift objects. We exploited the newly available JWST/NIRSpec and NIRCam dataset to perform an independent analysis of the Chandra data confirming that GHZ9 is the most likely JWST source associated to X-ray emission at 0.5-7 keV. Assuming a spectral index $\Gamma$ = 2.3 (1.8), we estimate a black hole (BH) mass of 1.60 $\pm$ 0.31 (0.48 $\pm$ 0.09) $\times$ 10$^8$M$_{\odot}$, which is consistent either with Eddington-accretion onto heavy ($\geq$ 10$^6$ M$_{\odot}$) BH seeds formed at z=18, or super-Eddington accretion onto a light seed of $\sim$ 10$^2-10^4$ M$_{\odot}$ at z = 25. The corresponding BH-to-stellar mass ratio M$_{BH}$/M$_{star}$= 0.33$\pm$0.22 (0.10$\pm$0.07), with a stringent limit $>$0.02, implies an accelerated growth of the BH mass with respect to the stellar mass. GHZ9 is the ideal target to constrain the early phases of AGN-galaxy coevolution with future multi-frequency observations., Comment: Submitted to ApJ

The origin of Lyman Continuum (LyC) photons responsible for reionizing the universe remains a mystery, with the fraction of escaping LyC photons from galaxies at z$\sim$ 6 to 12 being highly uncertain. While direct detection of LyC photons from this epoch is hindered by absorption from the intergalactic medium, lower redshift analogs offer a promising avenue to study LyC leakage. We present Hubble Space Telescope Cosmic Origins Spectrograph (HST COS) observations of five low redshift (z$\sim$ 0.3) massive starburst galaxies, selected for their high stellar mass and weak [SII] nebular emission - an indirect tracer of LyC escape. Three of the five galaxies show LyC leakage, highlighting the reliability of weak [SII] as a tracer, especially in light of recent JWST discoveries of z $>$ 5 galaxies with similarly weak [SII] emission. The dust corrected LyC escape fractions, which represent the LyC photons that would escape in the absence of dust, range from 33% to 84%. However, the absolute escape fractions, which show the LyC photons escaping after passing through both neutral hydrogen absorption and dust attenuation, are significantly lower, ranging between 1% and 3%. This suggests that while the galaxies are nearly optically thin to HI, their high dust content significantly suppresses LyC photon escape. These [SII] weak, massive leakers are distinct from typical low-redshift LyC emitters, showing higher metallicity, lower ionization states, more dust extinction and higher star formation surface densities. This suggests that these galaxies constitute a distinct population, likely governed by a different mechanism facilitating LyC photon escape. We propose that the feedback-driven winds in these compact starbursts create ionized channels through which LyC photons escape, aligning with a picket-fence model., Comment: 16 pages, 10 figures, 1 Table

We conducted systematic radiative transfer (RT) modeling of the Mg II doublet line profiles for 33 low-redshift Lyman continuum (LyC) leakers, and Ly$\alpha$ modeling for a subset of six objects, using a multiphase, clumpy circumgalactic medium (CGM) model. Our RT models successfully reproduced the Mg II line profiles for all 33 galaxies, revealing a necessary condition for strong LyC leakage: high maximum clump outflow velocity ($v_{\rm MgII,\,max} \gtrsim 390\,\rm km\,s^{-1}$) and low total Mg II column density ($N_{\rm MgII,\,tot} \lesssim 10^{14.3}\,\rm cm^{-2}$). We found that the clump outflow velocity and total Mg II column density have the most significant impact on Mg II spectra and emphasized the need for full RT modeling to accurately extract the CGM gas properties. In addition, using archival HST COS/G160M data, we modeled Ly$\alpha$ profiles for six objects and found that their spectral properties do not fully align with the conventional LyC leakage criteria, yet no clear correlation was identified between the modeled parameters and observed LyC escape fractions. We inferred LyC escape fractions based on HI properties from Ly$\alpha$ RT modeling and found that LyC leakage is primarily governed by the number of optically thick HI clumps per sightline ($f_{\rm cl}$). Intriguingly, two galaxies with relatively low observed LyC leakage exhibited the highest RT-inferred LyC escape fractions due to their lowest $f_{\rm cl}$ values, driven by the strong blue peaks of their Ly$\alpha$ emission. Future high-resolution, spatially resolved observations are crucial for resolving this puzzle. Overall, our results support a "picket fence" geometry over a "density-bounded" scenario for the CGM, where a combination of high Mg II outflow velocities and low Mg II column densities may be correlated with the presence of more low-density HI channels that facilitate LyC escape., Comment: 22 pages, 16 figures, comments are welcome

Rest-frame far-ultraviolet (FUV) observations from JWST are revolutionizing our understanding of the high-z galaxies that drove reionization and the mechanisms by which they accomplished it. To fully interpret these observations, we must be able to diagnose how properties of the interstellar medium (ISM; e.g., column density, covering fraction, outflow velocity) directly relate to the absorption features produced. Using the high-S/N and high-resolution FUV spectra of 45 nearby star-forming galaxies from CLASSY, we present the largest uniform, simultaneous characterization of neutral and low-ionization state (LIS) interstellar UV absorption lines (OI, SiII, SII, CII, AlII) across a wide range of galaxy properties. We also present 21-cm HI observations for 35 galaxies, multiple of which are gas-poor or non-detected, possibly indicating the onset of a post-starburst phase. We find that our simultaneous 1-component Voigt profile fits are capable of accurately modeling the LIS absorption for ~75% of galaxies, mitigating challenges associated with saturation, infilling, and degeneracies. While the most massive galaxies require additional components, our 1-component fits return average properties of the absorbing gas and follow the scaling relations described by a single gas cloud. We explore connections between LIS absorption and direct tracers of the neutral ISM (OI, Ly-alpha, HI 21-cm), finding that CII most closely traces the neutral gas trends while other ions exhibit weaker correlations. Given the challenges with directly observing HI at higher-z, we demonstrate that LIS absorption can be a powerful means to study the neutral ISM and present empirical relationships for predicting neutral gas properties., Comment: 25 pages with 16 figures and 2 tables. Long appendix with figure sets and tables. Accepted to ApJ

Galactic outflows provide important feedback effects to regulate the evolution of the host galaxies. Two primary diagnostics of galactic outflows are broad and/or blueshifted emission and absorption lines. Even though well-established methods exist to analyze these outflow signatures, connections between them are rarely studied and largely unknown. In this paper, we present the first detailed comparisons of the outflow properties measured independently from the two outflow diagnostics for a sample of 33 low-redshift star-forming galaxies. Their UV absorption lines are detected by the Hubble Space Telescope/Cosmic Origin Spectrograph, and optical emission lines are observed by the Keck/Echellette Spectrograph and Imager. We find that several outflow properties derived from emission and absorption lines are tightly correlated. These include outflow maximum velocity, line width, and sizes. Specifically, in a given galaxy, outflows seen in emission lines have smaller maximum velocities, narrower line widths, and smaller sizes than those measured from the absorption lines. These findings can be interpreted by the fact that emission line luminosity is weighted by density squared, while absorption line depth is weighted by density. We then test both spherical and bi-conical outflow models, and find the same outflow velocity and density distributions can explain the observed outflow features in emission and absorption lines for individual galaxies. These results provide novel calibration between galactic outflow properties measured from the two diagnostics and provide valuable insights for future models of galactic outflows by potentially doubling the number of observational constraints., Comment: 3 tables, 17 figures, 30 pages, submitted to ApJ

One of the fundamental questions of cosmology is the origin and mechanism(s) responsible for the reionization of the Universe beyond $z\sim6$. To address this question, many studies over the past decade have focused on local ($z\sim0.3$) galaxies which leak ionizing radiation (Lyman continuum or LyC). However, line-of-sight effects and data quality have prohibited deeper insight into the nature of LyC escape. To circumvent these limitations, we analyze stacks of a consolidated sample of {\it HST}/COS observations of the LyC in 89 galaxies at $z\sim0.3$. From fitting of the continuum, we obtain information about the underlying stellar populations and neutral ISM geometry. We find that most LyC non-detections are not leaking appreciable LyC ($f_{esc}^{\rm LyC}<1$\%) but also that exceptional cases point to spatial variations in the LyC escape fraction $f_{esc}^{\rm LyC}$. Stellar populations younger than 3 Myr lead to an increase in ionizing feedback, which in turn increases the isotropy of LyC escape. Moreover, mechanical feedback from supernovae in 8-10 Myr stellar populations is important for anisotropic gas distributions needed for LyC escape. While mechanical feedback is necessary for any LyC escape, high $f_{esc}^{\rm LyC}$ ($>5$\%) also requires a confluence of young stars and ionizing feedback. A two-stage burst of star formation could facilitate this optimal LyC escape scenario., Comment: Submitted to AAS publications, 42 pages, 25 figures

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg II 2796\r{A}, 2804\r{A} absorption + emission lines in 29 galaxies taken from the Low-z LyC Survey (LzLCS) to investigate the impact of (radiation + mechanical) feedback on LyC escape. Using constraints on Mg$^+$ and photoionization models, we map the outflows' neutral hydrogen content and predict $f_{esc}^{LyC}$ with a multiphase wind model. We measure mass, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in SFR-based deposition rates. We use SED template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high star formation rate surface density and compact UV half-light radii. Our findings indicate that the strongest leakers are radiation feedback dominant, lack Mg II outflows, but have extended broad components in higher ionization lines like [O III] 5007\r{A}, as observed by Amor\'in et al. (2024). In contrast, galaxies experiencing supernovae feedback typically exhibit weaker $f_{esc}^{LyC}$ and show evidence of outflows in both Mg II and higher ionization lines. We attribute these findings to rapid or "catastrophic" cooling in the radiation-dominant systems, which, given the low metallicities in our sample, are likely experiencing delayed supernovae., Comment: 34 pages, 16 figures, 7 tables

The radial gradient of gas-phase metallicity is a powerful probe of the chemical and structural evolution of star-forming galaxies, closely tied to disk formation and gas kinematics in the early universe. We present spatially resolved chemical and dynamical properties for a sample of 26 galaxies at $0.5 \lesssim z \lesssim 1.7$ from the MSA-3D survey. These innovative observations provide 3D spectroscopy of galaxies at a spatial resolution approaching JWST's diffraction limit and a high spectral resolution of $R\simeq2700$. The metallicity gradients measured in our galaxy sample range from $-$0.05 to 0.02 dex~kpc$^{-1}$. Most galaxies exhibit negative or flat radial gradients, indicating lower metallicity in the outskirts or uniform metallicity throughout the entire galaxy. We confirm a tight relationship between stellar mass and metallicity gradient at $z\sim1$ with small intrinsic scatter of 0.02 dex~kpc$^{-1}$. Our results indicate that metallicity gradients become increasingly negative as stellar mass increases, likely because the more massive galaxies tend to be more "disky". This relationship is consistent with the predictions from cosmological hydrodynamic zoom-in simulations with strong stellar feedback. This work presents the effort to harness the multiplexing capability of JWST NIRSpec/MSA in slit-stepping mode to map the chemical and kinematic profiles of high-redshift galaxies in large samples and at high spatial and spectral resolution., Comment: 23 pages, 6 figures, 1 table; submitted to ApJL

Integral field spectroscopy (IFS) is a powerful tool for understanding the formation of galaxies across cosmic history. We present the observing strategy and first results of MSA-3D, a novel JWST program using multi-object spectroscopy in a slit-stepping strategy to produce IFS data cubes. The program observed 43 normal star-forming galaxies at redshifts $0.5 \lesssim z \lesssim 1.5$, corresponding to the epoch when spiral thin-disk galaxies of the modern Hubble sequence are thought to emerge, obtaining kpc-scale maps of rest-frame optical nebular emission lines with spectral resolution $R\simeq2700$. Here we describe the multiplexed slit-stepping method which is $>15$ times more efficient than the NIRSpec IFS mode for our program. As an example of the data quality, we present a case study of an individual galaxy at $z=1.104$ (stellar mass $M_{*} = 10^{10.3}~M_{\odot}$, star formation rate~$=3~M_{\odot}$ yr$^{-1}$) with prominent face-on spiral structure. We show that the galaxy exhibits a rotationally supported disk with moderate velocity dispersion ($\sigma = 36^{+5}_{-4}$~\kms), a negative radial metallicity gradient ($-0.020\pm0.002$~dex\,kpc$^{-1}$), a dust attenuation gradient, and an exponential star formation rate density profile which closely matches the stellar continuum. These properties are characteristic of local spirals, indicating that mature galaxies are in place at $z\sim1$. We also describe the customized data reduction and original cube-building software pipelines which we have developed to exploit the powerful slit-stepping technique. Our results demonstrate the ability of JWST slit-stepping to study galaxy populations at intermediate to high redshifts, with data quality similar to current surveys of the $z\sim0.1$ universe., Comment: 24 pages, 12 figures, 1 table; submitted to ApJ, comments welcome. Our custom designed data reduction pipeline (MSA3D) and example processed data-cube will be publicly released before the JWST Cycle 4 deadline. MSA3D package on GitHub: https://github.com/barisiciv/msa3d.git

JWST is uncovering the properties of ever increasing numbers of galaxies at z>6, during the epoch of reionization. Connecting these observed populations to the process of reionization requires understanding how efficiently they produce Lyman continuum (LyC) photons and what fraction (fesc) of these photons escape into the intergalactic medium. By applying the Cox proportional hazards model, a survival analysis technique, to the Low-redshift Lyman Continuum Survey (LzLCS), we develop new, empirical, multivariate predictions for fesc. The models developed from the LzLCS reproduce the observed fesc for z~3 samples, which suggests that LyC emitters may share similar properties at low and high redshift. Our best-performing models for the z~3 galaxies include information about dust attenuation, ionization, and/or morphology. We then apply these models to z$\gtrsim$6 galaxies. For large photometric samples, we find a median predicted fesc=0.047-0.14. For smaller spectroscopic samples, which may include stronger emission line galaxies, we find that $\geq$33% of the galaxies have fesc >0.2, and we identify several candidate extreme leakers with fesc $\geq$0.5. The current samples show no strong trend between predicted fesc and UV magnitude, but limited spectroscopic information makes this result uncertain. Multivariate predictions can give significantly different results from single variable predictions, and the predicted fesc for high-redshift galaxies can differ significantly depending on whether star formation rate surface density or radius is used as a measure of galaxy morphology. We provide all parameters necessary to predict fesc for additional samples of high-redshift galaxies using these models., Comment: Accepted for publication in ApJ. 33 pages, 9 figures, 10 tables, plus appendix

To understand how galaxies reionized the universe, we must determine how the escape fraction of Lyman Continuum (LyC) photons (fesc) depends on galaxy properties. Using the z~0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of fesc. These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detections and upper limits. Our best model predicts the LzLCS fesc detections with a root-mean-square (RMS) scatter of 0.31 dex, better than single-variable correlations. According to ranking techniques, the most important predictors of fesc are the equivalent width (EW) of Lyman-series absorption lines and the UV dust attenuation, which track line-of-sight absorption due to HI and dust. The HI absorption EW is uniquely crucial for predicting fesc for the strongest LyC emitters, which show properties similar to weaker LyC emitters and whose high fesc may therefore result from favorable orientation. In the absence of HI information, star formation rate surface density ($\Sigma_{\rm SFR}$) and [O III]/[O II] ratio are the most predictive variables and highlight the connection between feedback and fesc. We generate a model suitable for z>6, which uses only the UV slope, $\Sigma_{\rm SFR}$, and [O III]/[O II]. We find that $\Sigma_{\rm SFR}$ is more important in predicting fesc at higher stellar masses, whereas [O III]/[O II] plays a greater role at lower masses. We also analyze predictions for other parameters, such as the ionizing-to-non ionizing flux ratio and Ly=alpha escape fraction. These multivariate models represent a promising tool for predicting fesc at high redshift., Comment: Accepted for publication in ApJ. 34 pages + appendix, 12 figures

The escape fraction of LyC ionizing radiation $f_{LyC}^{esc}$ is crucial for understanding reionization, yet impossible to measure at z $\gtrsim$ 5.3. Recently, studies have focused on calibrating indirect indicators of $f_{LyC}^{esc}$ at z $\sim$ 0.3, finding that Ly$\alpha$ is closely linked to it. What is still unclear is whether the LyC - Ly$\alpha$ relation evolves with redshift, and if Ly$\alpha$ is truly applicable as an $f_{LyC}^{esc}$ indicator during the reionization epoch. In this study, we investigate seven $-21 \lesssim M_{UV} \lesssim -19$ gravitationally lensed galaxies from the BELLS GALLERY survey at z $\sim$ 2.3. Our targets have rest-frame Ly$\alpha$ equivalent widths ranging from 10 \r{A} to 100 \r{A} and low dust content ($-2.5 \lesssim \beta \lesssim -1.9$), both indicative of high LyC escape. Surprisingly, direct estimates of $f_{LyC}^{esc}$ using Hubble Space Telescope imaging with F275W and F225W reveal that our targets are not LyC emitters, with absolute $f_{LyC}^{esc}$ $\lesssim$ 6.5% at 3$\sigma$ significance (with two sources having absolute $f_{LyC}^{esc}$(3$\sigma$) $\lesssim$ 10% and $\lesssim$ 16%). The low $f_{LyC}^{esc}$, coupled with the high Ly$\alpha$ escape fraction and equivalent width could potentially be attributed to the redshift evolution of the neutral hydrogen column density and dust content. Our results challenge previous studies based on local samples, suggesting that the extrapolation of z ~ 0 Ly$\alpha$-based LyC indirect estimators into the reionization epoch might not be correct., Comment: 19 Pages, 8 Figures. Submitted to ApJ

Using more than 100 galaxies in the MUSE Ultra Deep Field with spectroscopy from the Hubble Space Telescope's Wide Field Camera 3 and the Very Large Telescope's Multi Unit Spectroscopic Explorer, we extend the gas-phase mass-metallicity relation (MZR) at $z\approx\,$1$\,$-$\,$2 down to stellar masses of M$_{\star}$ $\approx$ 10$^{7.5}$ M$_{\odot}$. The sample reaches six times lower in stellar mass and star formation rate (SFR) than previous HST studies at these redshifts, and we find that galaxy metallicities decrease to log(O/H) + 12 $\approx$ 7.8 $\pm$ 0.1 (15% solar) at log(M$_{\star}$/M$_{\odot}$) $\approx$ 7.5, without evidence of a turnover in the shape of the MZR at low masses. We validate our strong-line metallicities using the direct method for sources with [O III] $\lambda$4363 and [O III] $\lambda$1666 detections, and find excellent agreement between the techniques. The [O III] $\lambda$1666-based metallicities double existing measurements with S/N $\geq$ 5 for unlensed sources at $z~>$ 1, validating the strong-line calibrations up to $z \sim$2.5. We confirm that the MZR resides $\sim$0.3 dex lower in metallicity than local galaxies and is consistent with the fundamental metallicity relation (FMR) if the low mass slope varies with SFR. At lower redshifts ($z\sim$0.5) our sample reaches $\sim$0.5 dex lower in SFR than current calibrations and we find enhanced metallicities that are consistent with extrapolating the MZR to lower SFRs. Finally, we detect only a $\sim$0.1 dex difference in the metallicities of galaxies in groups versus isolated environments. These results are based on robust calibrations and reach the lowest masses and SFRs that are accessible with HST, providing a critical foundation for studies with the Webb and Roman Space Telescopes., Comment: Accepted for publication in ApJ on March 23, 2024. The paper has 29 pages, 12 figures, and 6 tables. The calibrated data are available through MAST at: https://archive.stsci.edu/hlsp/mudf

Low-mass galaxies can significantly contribute to reionization due to their potentially high Lyman continuum (LyC) escape fraction and relatively high space density. We present a constraint on the LyC escape fraction from low-mass galaxies at z = 1.3 - 3.0. We obtained rest-frame UV continuum imaging with the ACS/SBC and the WFC3/UVIS from the Hubble Space Telescope for eight strongly-lensed galaxies that were identified in the Sloan Giant Arc Survey (SGAS) and the Cluster Lensing And Supernova survey with Hubble (CLASH). The targeted galaxies were selected to be spectroscopically confirmed, highly magnified, and blue in their UV spectral shapes ($\beta<-1.7$). Our targets include intrinsically low luminosity galaxies down to a magnification-corrected absolute UV magnitude of $M_{\rm UV}\sim-14$. We perform custom-defined aperture photometry to place the most reliable upper limits of LyC escape from our sample. From our observations, we report no significant ($>$$2\sigma$) detections of LyC fluxes, placing 1$\sigma$ upper limits on the absolute LyC escape fractions of 3 - 15%. Our observations do not support the expected increased escape fractions of LyC photons from intrinsically UV faint sources. Considering the highly anisotropic geometry of LyC escape, increasing the sample size of faint galaxies in future LyC observations is crucial., Comment: Submitted to ApJ (16 pages, 6 figures, 5 tables)

The geometry of the neutral gas in and around galaxies is a key regulator of the escape of ionizing photons. We present the first statistical study aiming at linking the neutral and ionized gas distributions to the Lyman continuum (LyC) escape fraction (fesc(LyC)) in a sample of 22 confirmed LyC leakers and non-leakers at z~0.35 using the Keck Cosmic Web Imager (Keck/KCWI) and the Low Resolution Spectrograph 2 (HET/LRS2). Our integral field unit data enable the detection of neutral and low-ionization gas, as traced by Mg II, and ionized gas, as traced by [O II], extending beyond the stellar continuum for 7 and 10 objects, respectively. All but one object with extended Mg II emission also shows extended [O II] emission; in this case, Mg II emission is always more extended than [O II] by a factor 1.3 on average. Most of the galaxies with extended emission are non or weak LyC leakers (fesc(LyC) < 5%), but we find a large diversity of neutral gas configurations around these weakly LyC-emitting galaxies. Conversely, the strongest leakers (fesc(LyC) > 10%) appear uniformly compact in both Mg II and [O II] with exponential scale lengths <1 kpc. We also find a trend between fesc(LyC) and the spatial offsets of the nebular gas and the stellar continuum emission. Moreover, we find significant anti-correlations between the spatial extent of the neutral gas and the [O III]/[O II] ratio, and H$\beta$ equivalent width, as well as positive correlations with metallicity and UV size, suggesting that galaxies with more compact neutral gas sizes are more highly ionized. The observations suggest that strong LyC emitters do not have extended neutral gas halos and ionizing photons may be emitted in many directions. Combined with high ionization diagnostics, we propose the Mg II, and potentially [O II], spatial compactness are indirect indicators of LyC emitting galaxies at high-redshift., Comment: 31 pages, 24 figures, 4 tables, submitted to A&A

Galactic winds play essential roles in the evolution of galaxies through the feedback they provide. Despite intensive studies of winds, the radial distributions of their properties and feedback are rarely observable. Here we present such measurements for the prototypical starburst galaxy, M 82, based on observations by Subaru telescope. We determine the radial distribution of outflow densities ($n_e$) from the spatially-resolved [S II] $\lambda\lambda$ 6717, 6731 emission-lines. We find $n_e$ drops from 200 to 40 cm$^{-3}$ with radius ($r$) between 0.5 and 2.2 kpc with a best-fit power-law index of $r^{-1.2}$. Combined with resolved H$\alpha$ lines, we derive mass, momentum, and energy outflow rates, which drop quite slowly (almost unchanged within error bars) over this range of $r$. This suggests that the galactic wind in M 82 can carry mass, momentum, and energy from the central regions to a few kpc with minimal losses. We further derive outflow cloud properties, including size and column densities. The clouds we measure have pressures and densities that are too high to match those from recent theoretical models and numerical simulations of winds. By comparing with a sample of outflows in local star-forming galaxies studied with UV absorption-lines, the above-derived properties for M 82 outflows match well with the published scaling relationships. These matches suggest that the ionized gas clouds traced in emission and absorption are strongly related. Our measurements motivate future spatially resolved studies of galactic winds, which is the only way to map the structure of their feedback effects., Comment: 15 pages, 9 figures, 3 tables; Accepted to ApJ

Lyman-alpha line profiles are a powerful probe of ISM structure, outflow speed, and Lyman continuum escape fraction. In this paper, we present the Ly$\alpha$ line profiles of the COS Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Ly$\alpha$ emission profile in the bottom of a damped, Ly$\alpha$ absorption trough. Such profiles reveal an inhomogeneous interstellar medium (ISM). We successfully fit the damped Ly$\alpha$ absorption (DLA) and the Ly$\alpha$ emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Ly$\alpha$ exchange between high-$N_\mathrm{HI}$ and low-$N_\mathrm{HI}$ paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Ly$\alpha$ peak separation and the [O III]/[O II] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Ly$\alpha$ peak separation decreases. We combine measurements of Ly$\alpha$ peak separation and Ly$\alpha$ red peak asymmetry in a diagnostic diagram which identifies six Lyman continuum leakers in the CLASSY sample. We find a strong correlation between the Ly$\alpha$ trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Ly$\alpha$ peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Ly$\alpha$ photons outside the spectroscopic aperture reshapes Ly$\alpha$ profiles as the distances to these compact starbursts span a large range., Comment: 40 pages, 19 figures, 5 tables, accepted for publication in ApJ, comments welcome

Observations of low-ionization state (LIS) metal lines provide crucial insights into the interstellar medium of galaxies, yet, disentangling the physical processes responsible for the emerging line profiles is difficult. This work investigates how mock spectra generated using a single galaxy in a radiation-hydrodynamical simulation can help us interpret observations of a real galaxy. We create 22,500 C II and Si II spectra from the virtual galaxy at different times and through multiple lines of sight and compare them with the 45 observations of low-redshift star-forming galaxies from the COS Legacy Spectroscopic SurveY (CLASSY). We find that the mock profiles provide accurate replicates to the observations of 38 galaxies with a broad range of stellar masses ($10^6$ to $10^9$ $M_\odot$) and metallicities (0.02 to 0.55 $Z_\odot$). Additionally, we highlight that aperture losses explain the weakness of the fluorescent emission in several CLASSY spectra and must be accounted for when comparing simulations to observations. Overall, we show that the evolution of a single simulated galaxy can produce a large diversity of spectra whose properties are representative of galaxies of comparable or smaller masses. Building upon these results, we explore the origin of the continuum, residual flux, and fluorescent emission in the simulation. We find that these different spectral features all emerge from distinct regions in the galaxy's ISM, and their characteristics can vary as a function of the viewing angle. While these outcomes challenge simplified interpretations of down-the-barrel spectra, our results indicate that high-resolution simulations provide an optimal framework to interpret these observations., Comment: Accepted for publication in Apj

We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $\lambda\lambda$1907,1909, [OII] $\lambda$3727, [NeIII] $\lambda$3869, [NeIII] $\lambda$3968, H$\delta$ $\lambda$4101, H$\gamma$ $\lambda$4340, and [OIII] $\lambda$4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity $Z =$ log(O/H) = $7.5 - 8.0$, or $(0.06 - 0.2)$ $Z_\odot$, ionization parameter log($U$) $\sim -1.9\pm0.2$, and an ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})=25.2\pm0.2\,$erg$^{-1}$ Hz. The spectrum has a softened Lyman-$\alpha$ break, evidence for a strong Ly$\alpha$ damping wing, suggesting that MACS0647-JD was unable to ionize its surroundings beyond its immediate vicinity ($R_{\text{HII}} \ll 1$ pMpc). The Ly$\alpha$ damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift $z = 10.6 \pm 0.3$. MACS0647-JD has a stellar mass log($M/M_\odot$) = $8.1 \pm 0.3$, including $\sim$ 6$\times 10^7 M_\odot$ in component A, most of which formed recently (within $\sim$ 20 Myr) with a star formation rate $2\pm1 M_\odot$ / yr, all within an effective radius $70\pm24\,$pc. The smaller component B ($r \sim 20$) pc is likely older ($\sim$100 Myr) with more dust ($A_V \sim 0.1$ mag), as found previously. Spectroscopy of a fainter companion galaxy C separated by a distance of \about\ 3$\,$kpc reveals a Lyman break consistent with $z = 10.17$. MACS0647-JD is likely the most distant galaxy merger known., Comment: 22 pages, 6 figures, accepted for publication in ApJ

Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit., Comment: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Lyman-alpha (Ly$\alpha$) emission from galaxies can be used to trace neutral hydrogen in the epoch of reionization, however, there is a degeneracy between the attenuation of Ly$\alpha$ in the intergalactic medium (IGM) and the line profile emitted from the galaxy. Large shifts of Ly$\alpha$ redward of systemic due to scattering in the interstellar medium can boost Ly$\alpha$ transmission in the IGM during reionization. The relationship between Ly$\alpha$ velocity offset from systemic and other galaxy properties is not well-established at high-redshift or low luminosities, due to the difficulty of observing emission lines which trace systemic redshift. Rest-frame optical spectroscopy with JWST/NIRSpec has opened a new window into understanding of Ly$\alpha$ at z>3. We present a sample of 12 UV-faint galaxies ($-20 \lesssim$ MUV $\lesssim -16$) at $3 \lesssim z \lesssim 6$, with Ly$\alpha$ velocity offsets, $\Delta v_{\mathrm{Ly}\alpha}$, measured from VLT/MUSE and JWST/NIRSpec from the GLASS-JWST Early Release Program. We find median $\Delta v_{\mathrm{Ly}\alpha}$ of 205 km s$^{-1}$ and standard deviation 75 km s$^{-1}$, compared to 320 and 170km s$^{-1}$ for MUV < -20 galaxies in the literature. Our new sample demonstrates the previously observed trend of decreasing Ly$\alpha$ velocity offset with decreasing UV luminosity and optical line velocity dispersion, extends to MUV $\gtrsim$ -20, consistent with a picture where the Ly$\alpha$ profile is shaped by gas close to the systemic redshift. Our results imply that during reionization Ly$\alpha$ from UV-faint galaxies will be preferentially attenuated, but that detecting Ly$\alpha$ with low $\Delta v_{\mathrm{Ly}\alpha}$ can be an indicator of large ionized bubbles., Comment: 13 pages, 4 figures

Ly$\alpha$ emission is possibly the best indirect diagnostic of Lyman continuum (LyC) escape since the conditions that favor the escape of Ly$\alpha$ photons are often the same that allows for the escape of LyC photons. In this work, we present the rest UV-optical spectral characteristics of 11 Ly$\alpha$ emitting galaxies at 3 $<$ z $<$ 6 - the optimal redshift range chosen to avoid the extreme IGM attenuation while simultaneously studying galaxies close enough to the epoch of reionization. From a combined analysis of JWST/NIRSpec and MUSE data, we present the Ly$\alpha$ escape fraction and study their correlations with other physical properties of galaxies that might facilitate Ly$\alpha$ escape. We find that our galaxies have low masses (80\% of the sample with $\rm log_{10} \ M_{\star} < 9.5\ M_{\odot}$), compact sizes (median $\rm R_e \sim 0.7 \ kpc $), low dust content, moderate [OIII]/[OII] flux ratios (mean $\sim$ 6.8 $\pm$ 1.2), and moderate Ly$\alpha$ escape fraction (mean $\rm f_{esc}^{Ly\alpha} \ \sim$ 0.11). Our sample show characteristics that are broadly consistent with the low redshift galaxies with Ly$\alpha$ emission, which are termed as "analogs" of high redshift population. We predict the Lyman continuum escape fraction in our sample to be low (0.03-0.07), although larger samples in the post-reionization epoch are needed to confirm these trends., Comment: 15 pages, 4 figures, 1 table; Submitted to ApJL

We present results of [CII]$\,158\,\rm{\mu m}$ emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed ($\mu\sim20$) star-forming galaxy, MACS0308-zD1 at $z=6.2078\pm0.0002$. The [CII] emission line is detected with a signal-to-noise ratio $>6$ within the rest-frame UV bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the narrow [CII] has a velocity full-width-half-maximum (FWHM) of $110\pm20\,\rm{km/s}$, while broader [CII] is seen with an FWHM of $230\pm20\,\rm{km/s}$. The broader [CII] component is blueshifted ($-80\pm20\,\rm{km/s}$) with respect to the narrow [CII] component, and has a morphology which extends beyond the UV-bright clump. We find that while the narrow [CII] emission is most likely associated with zD1.1, the broader component is possibly associated with outflowing gas. Based on the non-detection of $\lambda_{\rm 158\,\mu m}$ dust continuum, we find that MACS0308-zD1's star-formation activity occurs in a dust-free environment with the stringent upper limit of infrared luminosity $\lesssim9\times10^{8}\,{\rm L_{\odot}}$. Targeting this strongly lensed faint galaxy for follow-up ALMA and JWST observations will be crucial to characterize the details of typical galaxy growth in the early Universe., Comment: 10 pages, 6 figures, 1 table; submitted to ApJ

JWST observations confirm the existence of galaxies as early as 300Myr and at a higher number density than expected based on galaxy formation models and HST observations. Yet, sources confirmed spectroscopically in the first 500Myr have estimated stellar masses $<5\times10^8M_\odot$, limiting the signal to noise ratio (SNR) for investigating substructure. We present a high-resolution spectroscopic and spatially resolved study of a rare bright galaxy at $z=9.3127\pm0.0002$ with a stellar mass of $(2.5^{+0.7}_{-0.5})\times10^9M_\odot$, forming $25^{+3}_{-4}M_\odot/yr$ and with a metallicity of $\sim0.1Z_\odot$- lower than in the local universe for the stellar mass but in line with expectations of chemical enrichment in galaxies 1-2Gyr after the Big Bang. The system has a morphology typically associated to two interacting galaxies, with a two-component main clump of very young stars (age$<10$Myr) surrounded by an extended stellar population ($130\pm20$Myr old, identified by modeling the NIRSpec spectrum) and an elongated clumpy tidal tail. The spectroscopic observations identify O, Ne and H emission lines, and the Lyman break, where there is evidence of substantial Ly$\alpha$ absorption. The [OII] doublet is resolved spectrally, enabling an estimate of the electron number density and ionization parameter of the interstellar medium and showing higher densities and ionization than in lower redshift analogs. For the first time at $z>8$, we identify evidence of absorption lines (Si, C and Fe), with low confidence individual detections but SNR$>6$ when stacked. The absorption features suggest that Ly$\alpha$ is damped by the interstellar and circumgalactic medium. Our observations provide evidence of rapid efficient build-up of mass and metals in the immediate aftermath of the Big Bang through mergers, demonstrating that massive galaxies with several billion stars exist earlier than expected., Comment: 52 pages, 10 figures This version of the article has been accepted for publication, after peer review and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at: http://dx.doi.org/10.1038/s41550-024-02218-7

We present extremely deep Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) observations of the MUSE Ultra Deep Field (MUDF). This unique region of the sky contains two quasars at $z \approx$ 3.22 that are separated by only $\sim$500 kpc, providing a stereoscopic view of gas and galaxies in emission and absorption across $\sim$10 billion years of cosmic time. We have obtained 90 orbits of HST WFC3 G141 near-infrared grism spectroscopy of this field in a single pointing, as well as 142 hours of optical spectroscopy with the Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE). The WFC3 (F140W, F125W, and F336W) and archival WFPC2 (F702W and F450W) imaging provides five-filter photometry that we use to detect 3,375 sources between $z \approx$ 0 - 6, including 1,536 objects in a deep central pointing with both spectroscopic and photometric coverage. The F140W and F336W mosaics reach exceptional depths of $m_\mathrm{AB}\approx$ 28 and 29, respectively, providing near-infrared and rest-frame ultraviolet information for 1,580 sources, and we reach 5$\sigma$ continuum detections for objects as faint as $m_\mathrm{AB}\approx$ 27 in the grism spectra. The extensive wavelength coverage of MUSE and WFC3 allows us to measure spectroscopic redshifts for 419 sources, down to galaxy stellar masses of log(M/M$_{\odot}$) $\approx$ 7 at $z \approx$ 1 - 2. In this publication, we provide the calibrated HST data and source catalogs as High Level Science Products for use by the community, which includes photometry, morphology, and redshift measurements that enable a variety of studies aimed at advancing our models of galaxy formation and evolution in different environments., Comment: Accepted for publication in ApJS on February 1, 2023. The paper has 25 pages, 13 figures, and 7 tables. Version two includes minor corrections to match the journal publication. The calibrated data are available through MAST at: https://archive.stsci.edu/hlsp/mudf

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

The Mg II 2796, 2803 doublet has been suggested to be a useful indirect indicator for the escape of Ly-alpha and Lyman continuum (LyC) photons in local star-forming galaxies. However, studies to date have focused on small samples of galaxies with strong Mg II or strong LyC emission. Here we present the first study of Mg II probing a large dynamic range of galaxy properties, using newly obtained high signal-to-noise, moderate-resolution spectra of Mg II for a sample of 34 galaxies selected from the Low-redshift Lyman Continuum Survey. We show that the galaxies in our sample have Mg II profiles ranging from strong emission to P-Cygni profiles, and to pure absorption. We find there is a significant trend (with a possibility of spurious correlations of ~ 2%) that galaxies detected as strong LyC Emitters (LCEs) also show larger equivalent widths of Mg II emission, and non-LCEs tend to show evidence of more scattering and absorption features in Mg II We then find Mg II strongly correlates with Ly-alpha in both equivalent width and escape fraction, regardless of whether the emission or absorption dominates the Mg II profiles. Furthermore, we present that, for galaxies categorized as Mg II emitters (MgE), one can adopt the information of Mg II, metallicity, and dust to estimate the escape fraction of LyC within a factor of 3. These findings confirm that Mg II lines can be used as a tool to select galaxies as LCEs and to serve as an indirect indicator for the escape of Ly-alpha and LyC., Comment: 12 pages, 5 figures, 2 tables, accepted to ApJ

We present the spectroscopic confirmation of a protocluster at $z=7.88$ behind the galaxy cluster Abell2744 (hereafter A2744-z7p9OD). Using JWST NIRSpec, we find seven galaxies within a projected radius of 60kpc. Although the galaxies reside in an overdensity around $>20\times$ greater than a random volume, they do not show strong Lyman-alpha emission. We place 2-$\sigma$ upper limits on the rest-frame equivalent width $<16$-$28$AA. Based on the tight upper limits to the Lyman-alpha emission, we constrain the volume-averaged neutral fraction of hydrogen in the intergalactic medium to be $x_{\rm HI} > 0.45$ (68% CI). Using an empirical $M_{\rm UV}$-$M_{\rm halo}$ relation for individual galaxies, we estimate that the total halo mass of the system is $\gtrsim 4\times10^{11}\,M_\odot$. Likewise, the line of sight velocity dispersion is estimated to be $1100 \pm 200$km/s. Using an empirical relation, we estimate the present-day halo mass of A2744-z7p9OD to be $\sim2\times10^{15}\,M_\odot$, comparable to the Coma cluster. A2744-z7p9OD is the highest redshift spectroscopically confirmed protocluster to date, demonstrating the power of JWST to investigate the connection between dark-matter halo assembly and galaxy formation at very early times with medium-deep observations at $<20$hrs total exposure time. Follow-up spectroscopy of the remaining photometric candidates of the overdensity will further refine the features of this system and help characterize the role of such overdensities in cosmic reionization., Comment: The title has been updated to reflect the published numbering; a minor change has been made to Figure 1 with regard to the MSA shutters on the rgb stamp image. NASA press release article can be found at: https://www.nasa.gov/feature/goddard/2023/webb-reveals-early-universe-prequel-to-huge-galaxy-cluster

In the first billion years after the Big Bang, sources of ultraviolet (UV) photons are believed to have ionized intergalactic hydrogen, rendering the Universe transparent to UV radiation. Galaxies brighter than the characteristic luminosity $L^{*}$ do not provide enough ionizing photons to drive this cosmic reionization. Fainter galaxies are thought to dominate the photon budget; however they are surrounded by neutral gas that prevents the escape of the Lyman-$\alpha$ photons, which has been the dominant way to identify them so far. JD1 was previously identified as a triply-imaged galaxy with a magnification factor of 13 provided by the foreground cluster Abell 2744, and a photometric redshift of $z\sim10$. Here we report the spectroscopic confirmation of this very low luminosity ($\sim0.05 L^{*}$) galaxy at $z=9.79$, observed 480 Myr after the Big Bang, by means of the identification of the Lyman break and redward continuum, as well as multiple $\gtrsim4\sigma$ emission lines, with the Near-InfraRed Spectrograph (NIRSpec) and Near-InfraRed Camera (NIRCam) instruments. The combination of the James Webb Space Telescope (JWST) and gravitational lensing shows that this ultra-faint galaxy ($M_{\rm UV}=-17.35$) -- with a luminosity typical of the sources responsible for cosmic reionization -- has a compact ($\sim$150 pc) and complex morphology, low stellar mass (10$^{7.19}$ M$_\odot$), and subsolar ($\sim$0.6 $Z_{\odot}$) gas-phase metallicity., Comment: Matched to the published Nature version of the article. 19 pages, 4 main figures, 1 supplementary figure, 1 supplementary tables. This version includes an updated version of the NIRSpec spectrum and NIRCam photometry

MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A" is intrinsically very blue ($\beta\sim-2.6$), likely due to very recent star formation and no dust, and is spatially extended with an effective radius $\sim70\,\rm pc$. The smaller component "B" appears redder ($\beta\sim-2$), likely because it is older ($100-200\,\rm Myr$) with mild dust extinction ($A_V\sim0.1\,\rm mag$), and a smaller radius $\sim20\,\rm pc$. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be out of phase. With an estimated stellar mass ratio of roughly 2:1 and physical projected separation $\sim400\,\rm pc$, we may be witnessing a galaxy merger 400 million years after the Big Bang. We also identify a candidate companion galaxy C $\sim3\,{\rm kpc}$ away, likely destined to merge with galaxies A and B. The combined light from galaxies A+B is magnified by factors of $\sim$8, 5, and 2 in three lensed images JD1, 2, and 3 with F356W fluxes $\sim322$, $203$, $86\,\rm nJy$ (AB mag 25.1, 25.6, 26.6). MACS0647$-$JD is significantly brighter than other galaxies recently discovered at similar redshifts with JWST. Without magnification, it would have AB mag 27.3 ($M_{UV}=-20.4$). With a high confidence level, we obtain a photometric redshift of $z=10.6\pm0.3$ based on photometry measured in 6 NIRCam filters spanning $1-5\rm\mu m$, out to $4300\,\r{A}$ rest-frame. JWST NIRSpec observations planned for January 2023 will deliver a spectroscopic redshift and a more detailed study of the physical properties of MACS0647$-$JD., Comment: 27 pages, 14 figures, accepted for publication in ApJL

Absorption line spectroscopy offers one of the best opportunities to constrain the properties of galactic outflows and the environment of the circumgalactic medium. Extracting physical information from line profiles is difficult, however, for the physics governing the underlying radiation transfer is complicated and depends on many different parameters. Idealized analytical models are necessary to constrain the large parameter spaces efficiently, but are typically plagued by model degeneracy and systematic errors. Comparison tests with idealized numerical radiation transfer codes offer an excellent opportunity to confront both of these issues. In this paper, we present a detailed comparison between SALT, an analytical radiation transfer model for predicting UV spectra of galactic outflows, with the numerical radiation transfer software, RASCAS. Our analysis has lead to upgrades to both models including an improved derivation of SALT and a customizable adaptive mesh refinement routine for RASCAS. We explore how well SALT, when paired with a Monte Carlo fitting procedure, can recover flow parameters from non-turbulent and turbulent flows. When the velocity and density gradients are excluded, we find that flow parameters are well recovered from high resolution (20 $\rm{km}$ $\rm{s}^{-1}$) data and moderately well from medium resolution (100 $\rm{km}$ $\rm{s}^{-1}$) data without turbulence at a S/N = 10, while derived quantities (e.g., mass outflow rates, column density, etc.) are well recovered at all resolutions. In the turbulent case, biased errors emerge in the recovery of individual parameters, but derived quantities are still well recovered.

The COS Legacy Archive Spectroscopic SurveY (CLASSY) HST/COS treasury program provides the first high-resolution spectral catalogue of 45 local high-z analogues in the UV (1200-2000\AA) to investigate their stellar and gas properties. We present a toolkit of UV interstellar medium (ISM) diagnostics, analyzing the main emission lines of CLASSY spectra (i.e., NIV]$\lambda\lambda$1483,87, CIV$\lambda\lambda$1548,51, HeII$\lambda$1640, OIII]$\lambda\lambda$1661,6, SiIII]$\lambda\lambda$1883,92, CIII]$\lambda\lambda$1907,9). Specifically, we focus our investigation on providing accurate diagnostics for reddening, electron density and temperature, gas-phase metallicity and ionization parameter, taking into account the different ionization zones of the ISM. We calibrate our UV toolkit using well-known optical diagnostics, analyzing archival optical spectra for all the CLASSY targets. We find that UV density diagnostics estimate ne values that are ~1-2 dex higher (e.g., ne(CIII]$\lambda\lambda$}1907,9)~10$^4$cm$^{-3}$) than those inferred from their optical counterparts (e.g., ne([SII]$\lambda\lambda$6717,31)~10$^2$cm$^{-3}$). Te derived from the hybrid ratio OIII]$\lambda$1666/[OIII]$\lambda$}5007 proves to be a reliable Te diagnostic, with differences in 12+log(O/H) within ~$\pm$0.3dex. We also investigate the relation between the stellar and gas E(B-V), finding consistent values at high specific star formation rates, while at low sSFR we confirm an excess of dust attenuation in the gas. Finally, we investigate UV line ratios and equivalent widths to provide correlations with 12+log(O/H) and log(U), but note there are degeneracies between the two. With this suite of UV-based diagnostics, we illustrate the pivotal role CLASSY plays in understanding the chemical and physical properties of high-z systems that JWST can observe in the rest-frame UV., Comment: Accepted for publication on ApJ

The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$\mu$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $\mu > 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe., Comment: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

We report the first gas-phase metallicity map of a distant galaxy measured with the James Webb Space Telescope (JWST). We use the NIRISS slitless spectroscopy acquired by the GLASS Early Release Science program to spatially resolve the rest-frame optical nebular emission lines in a gravitationally lensed galaxy at $z=3.06$ behind the Abell 2744 galaxy cluster. This galaxy (dubbed GLASS-Zgrad1) has stellar mass $\sim10^{8.6} M_\odot$, instantaneous star formation rate $\sim8.6$ $M_\odot$/yr (both corrected for lensing magnification), and global metallicity one-fourth solar. From its emission line maps ([O III], H$\beta$, H$\gamma$, [Ne III], and [O II]) we derive its spatial distribution of gas-phase metallicity using a well-established forward-modeling Bayesian inference method. The exquisite resolution and sensitivity of JWST/NIRISS, combined with lensing magnification, enable us to resolve this $z\sim3$ dwarf galaxy in $\gtrsim$50 resolution elements with sufficient signal, an analysis hitherto not possible. We find that the radial metallicity gradient of GLASS-Zgrad1 is strongly inverted (i.e. positive): $\Delta\log({\rm O/H})/\Delta r$ = $0.165\pm0.023$ $\mathrm{dex~kpc^{-1}}$. This inverted gradient may be due to tidal torques induced by a massive nearby ($\sim$15 kpc projected) galaxy, which can cause inflows of metal-poor gas into the central regions of GLASS-Zgrad1. These first results showcase the power of JWST wide-field slitless spectroscopic modes to resolve the mass assembly and chemical enrichment of low-mass galaxies in and beyond the peak epoch of cosmic star formation ($z\gtrsim2$). Reaching masses $\lesssim 10^9~M_\odot$ at these redshifts is especially valuable to constrain the effects of galactic feedback and environment, and is possible only with JWST's new capabilities., Comment: Accepted to ApJL. 11 pages, 4 figures, and 1 table

The James Webb Space Telescope (JWST) Early Release Observations (EROs) is a set of public outreach products created to mark the end of commissioning and the beginning of science operations for JWST. Colloquially known as the "Webb First Images and Spectra", these products were intended to demonstrate to the worldwide public that JWST is ready for science, and is capable of producing spectacular results. The package was released on July 12, 2022, and included images and spectra of the galaxy cluster SMACS~J0723.3-7327 and distant lensed galaxies, the interacting galaxy group Stephan's Quintet, NGC 3324 in the Carina star-forming complex, the Southern Ring planetary nebula NGC 3132, and the transiting hot Jupiter WASP 96b. This paper describes the ERO technical design, observations, and scientific processing of data underlying the colorful outreach products., Comment: 15 pages, accepted by ApJ Letters

We present the first search for $z\geqslant7$, continuum-confirmed sources with NIRISS/WFS spectroscopy over the Abell 2744 Frontier Fields cluster, as part of the GLASS-JWST ERS survey. With $\sim15$ hrs of pre-imaging and multi-angle grism exposures in the F115W, F150W, and F200W filters, we describe the general data handling (i.e., reduction, cleaning, modeling, and extraction processes) and analysis for the GLASS-JWST survey. We showcase the power of JWST to peer deep into reionization, when most intergalactic hydrogen is neutral, by confirming two galaxies at $z=8.04\pm0.15$ and $z=7.90\pm0.13$ by means of their Lyman breaks. Fainter continuum spectra are observed in both the F150W and F200W bands, indicative of blue ($-1.69$ and $-1.33$) UV slopes and moderately-bright absolute magnitudes ($-20.37$ and $-19.68$ mag). We do not detect strong Ly$\alpha$ in either galaxy, but do observe tentative ($\sim2.7-3.8\sigma$) HeII$\lambda$1640 A, OIII]$\lambda\lambda$1661,1666 A, and NIII]$\lambda\lambda$1747,1749 A line emission in one, suggestive of low metallicity, star-forming systems with possible non-thermal contributions. These novel observations provide a first look at the extraordinary potential of JWST/NIRISS for confirming representative samples of bright $z\geqslant7$ sources in the absence of strong emission lines, and gain unprecedented insight into their contributions towards cosmic reionization., Comment: Accepted for publication in ApJL on 1st September 2022, after minor modifications based on the referee report. 10 pages, 3 figures, 1 table

This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies., Comment: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

The COS Legacy Archive Spectroscopic SurveY (CLASSY) is designed to provide the community with a spectral atlas of 45 nearby star-forming galaxies which were chosen to cover similar properties as those seen at high-z (z>6). The prime high level science product of CLASSY is accurately coadded UV spectra, ranging from ~1000-2000A, derived from a combination of archival and new data obtained with HST's Cosmic Origins Spectrograph (COS). This paper details the multi-stage technical processes of creating this prime data product, and the methodologies involved in extracting, reducing, aligning, and coadding far-ultraviolet (FUV) and near-ultraviolet (NUV) spectra. We provide guidelines on how to successfully utilize COS observations of extended sources, despite COS being optimized for point sources, and best-practice recommendations for the coaddition of UV spectra in general. Moreover, we discuss the effects of our reduction and coaddition techniques in the scientific application of the CLASSY data. In particular, we find that accurately accounting for flux calibration offsets can affect the derived properties of the stellar populations, while customized extractions of NUV spectra for extended sources are essential for correctly diagnosing the metallicity of galaxies via CIII] nebular emission. Despite changes in spectral resolution of up to ~25% between individual datasets (due to changes in the COS line spread function), no adverse affects were observed on the difference in velocity width and outflow velocities of isolated absorption lines when measured in the final combined data products, owing in-part to our signal-to-noise regime of S/N<20., Comment: 31 pages, 17 figures, accepted for publication in ApJS

Star-forming galaxies are considered the likeliest source of the H I ionizing Lyman Continuum (LyC) photons that reionized the intergalactic medium at high redshifts. However, above z >~ 6, the neutral intergalactic medium prevents direct observations of LyC. Therefore, recent years have seen the development of indirect indicators for LyC that can be calibrated at lower redshifts and applied in the Epoch of Reionization. Emission from Mg II \ly\ly 2796, 2803 doublet has been proposed as a promising LyC proxy. In this paper, we present new Hubble Space Telescope/Cosmic Origins Spectrograph observations for 8 LyC emitter candidates, selected to have strong Mg II emission lines. We securely detect LyC emission in 50% (4/8) galaxies with 2$\sigma$ significance. This high detection rate suggests that strong Mg II emitters might be more likely to leak LyC than similar galaxies without strong Mg II. Using photoionization models, we constrain the escape fraction of Mg II as ~ 15 -- 60%. We confirm that the escape fraction of Mg II correlates tightly with that of Lyman-alpha (LyA), which we interpret as an indication that the escape fraction of both species is controlled by resonant scattering in the same low column density gas. Furthermore, we show that the combination of the Mg II emission and dust attenuation can be used to estimate the escape fraction of LyC statistically. These findings confirm that Mg II emission can be adopted to estimate the escape fraction of LyA and LyC in local star-forming galaxies and may serve as a useful indirect indicator at the Epoch of Reionization., Comment: 24 pages, 15 figures, 7 tables, accepted to ApJ

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

We report the first spatially resolved measurements of gas-phase metallicity radial gradients in star-forming galaxies in overdense environments at $z\gtrsim2$. The spectroscopic data are acquired by the \mg\ survey, a Hubble Space Telescope (HST) cycle-28 medium program. This program is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters (BOSS 1244, BOSS 1542 and BOSS 1441) at $z=2-3$. Our sample in the BOSS 1244 field consists of 20 galaxies with stellar-mass ranging from $10^{9.0}$ to $10^{10.3}$ \Msun\ , star formation rate (SFR) from 10 to 240 \Msun\,yr$^{-1}$, and global gas-phase metallicity (\oh) from 8.2 to 8.6. At $1\sigma$ confidence level, 2/20 galaxies in our sample show positive (inverted) gradients -- the relative abundance of oxygen increasing with galactocentric radius, opposite the usual trend. Furthermore, 1/20 shows negative gradients and 17/20 are consistent with flat gradients. This high fraction of flat/inverted gradients is uncommon in simulations and previous observations conducted in blank fields at similar redshifts. To understand this, we investigate the correlations among various observed properties of our sample galaxies. We find an anticorrelation between metallicity gradient and global metallicity of our galaxies residing in extreme overdensities, and a marked deficiency of metallicity in our massive galaxies as compared to their coeval field counterparts. We conclude that the cold-mode gas accretion plays an active role in shaping the chemical evolution of galaxies in the protocluster environments, diluting their central chemical abundance, and flattening/inverting their metallicity gradients., Comment: 13 pages, 6 figures, 1 table. Accepted for publication in ApJL

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Ly{\alpha}, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape., Comment: ApJ, accepted. 31 pages, 26 figures, 1 table

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time., Comment: Accepted for publication in ApJ

Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV attenuation, by fitting the stellar continuum with a combination of simple stellar population models. The models, together with accurate LyC flux measurements, allow to determine the absolute LyC photon escape fraction for each galaxy ($f_{\rm esc}^{\rm abs}$). We measure the equivalent widths and residual fluxes of multiple HI and low-ionization state (LIS) lines, and the geometrical covering fraction adopting the picket-fence model. The $f_{\rm esc}^{\rm abs}$ spans a wide range, with a median (0.16, 0.84 quantiles) of 0.04 (0.02, 0.20), and 50 out of the 89 galaxies detected in the LyC. The HI and LIS line equivalent widths scale with the UV luminosity and attenuation, and inversely with the residual flux of the lines. The HI and LIS residual fluxes are correlated, indicating that the neutral gas is spatially traced by the LIS transitions. We find the observed trends of the absorption lines and the UV attenuation are primarily driven by the covering fraction. The non-uniform gas coverage demonstrates that LyC photons escape through low-column density channels in the ISM. The equivalent widths and residual fluxes of the UV lines strongly correlate with $f_{\rm esc}^{\rm abs}$: strong LyC leakers show weak absorption lines, low UV attenuation, and large Ly$\alpha$ equivalent widths. We finally show that simultaneous UV absorption line and dust attenuation measurements can predict, on average, the escape fraction of galaxies and the method can be applied to galaxies across a wide redshift range., Comment: 30 pages, 16 figures, 3 tables; accepted for publication in Astronomy and Astrophysics on December 16, 2021. Tables A1 to A4 are part of the LzLCS science products and will be publicly available in a dedicated website