Search

To study the chemical evolution across cosmic epochs, we investigate Ne, S, Cl, and Ar abundance patterns in the COS Legacy Archive Spectroscopic SurveY (CLASSY). CLASSY comprises local star-forming galaxies (0.02 < z < 0.18) with enhanced star-formation rates, making them strong analogues to high-z star-forming galaxies. With direct measurements of electron temperature, we derive accurate ionic abundances for all elements and assess ionization correction factors (ICFs) to account for unseen ions and derive total abundances. We find Ne/O, S/O, Cl/O, and Ar/O exhibit constant trends with gas-phase metallicity for 12+log(O/H) < 8.5 but significant correlation for Ne/O and Ar/O with metallicity for 12+log(O/H) > 8.5, likely due to ICFs. Thus, applicability of the ICFs to integrated spectra of galaxies could bias results, underestimating true abundance ratios. Using CLASSY as a local reference, we assess the evolution of Ne/O, S/O, and Ar/O in galaxies at z>3, finding no cosmic evolution of Ne/O, while the lack of direct abundance determinations for S/O and Ar/O can bias the interpretation of the evolution of these elements. We determine the fundamental metallicity relationship (FMR) for CLASSY and compare to the high-redshift FMR, finding no evolution. Finally, we perform the first mass-neon relationship analysis across cosmic epochs, finding a slight evolution to high Ne at later epochs. The robust abundance patterns of CLASSY galaxies and their broad range of physical properties provide essential benchmarks for interpreting the chemical enrichment of the early galaxies observed with the JWST., Comment: Accepted for publication in The Astrophysical Journal. 20 pages (main body), 10 figures, 6 Tables

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

Ultraviolet (UV; rest-frame $\sim1200-2000$ A) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at $5.6

We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $\sigma(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $\sigma(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$., Comment: 28 pages, 12 figures, 2 Appendix figures

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 45000 AA. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: all spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the SED-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hb + [OIII] and Ha in the neighboring broad-band filters. We present examples of EELGs that could be incorrectly classified at ultra-high redshift (z>12) as a result of extreme Hb + [OIII] emission blended across the reddest photometric filters. We compare the EELGs to the broader (sub-extreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high equivalent width tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. Photometrically inferred emission-line ratios are consistent with ISM conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies., Comment: 24 pages, 25 figures, submitted to ApJ

The wavelength-coverage and sensitivity of JWST now enables us to probe the rest-frame UV - optical spectral energy distributions (SEDs) of galaxies at high-redshift ($z>4$). From these SEDs it is, in principle, through SED fitting possible to infer key physical properties, including stellar masses, star formation rates, and dust attenuation. These in turn can be compared with the predictions of galaxy formation simulations allowing us to validate and refine the incorporated physics. However, the inference of physical properties, particularly from photometry alone, can lead to large uncertainties and potential biases. Instead, it is now possible, and common, for simulations to be \emph{forward-modelled} to yield synthetic observations that can be compared directly to real observations. In this work, we measure the JWST broadband fluxes and colours of a robust sample of $58$ the distributions differ somewhat, though our observed sample size is small and thus susceptible to statistical fluctuations. Likewise, the predicted and observed colour evolution show broad agreement, at least at $58$, though, again, the sample size is small here., Comment: 11 pages, 10 figures, submitted to MNRAS

Sources that leak Lyman-continuum (LyC) photons and lead to the reionisation of the universe are intensely studied using multiple observing facilities. Recently, the Low-redshift LyC Survey (LzLCS) has found the first large sample of LyC emitting galaxies at low redshift ($z\sim 0.3$) with the Hubble Space Telescope/Cosmic Origins Spectrograph. The LzLCS sample contains a robust estimate of the LyC escape fraction ($f_\mathrm{esc}^\mathrm{LyC}$) for 66 galaxies spanning a wide range of $f_\mathrm{esc}^\mathrm{LyC}$. Here we, for the first time, aim to study the radio continuum (RC) properties of LzLCS sources and their dependence on $f_\mathrm{esc}^\mathrm{LyC}$. We present Karl G. Jansky Very Large Array RC observations at C (4-8 GHz), S (2-4 GHz) and L (1-2 GHz) bands for a sub-sample of the LzLCS sources. The radio spectral index ($\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$) spans a wide range from being flat ( $\geq -0.1$) to very steep ($\leq -1.0$). The strongest leakers in our sample show flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$, weak leakers have $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ close to normal star-forming galaxies, and non-leakers are characterized by steep $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$. We argue that a combination of young ages, free-free absorption, and a flat cosmic-ray energy spectrum can altogether lead to a flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ for strong leakers. Non-leakers are characterized by steep spectra which can arise due to break/cutoff at high frequencies. Such a cutoff in the spectrum can arise in a single injection model of CRs characteristic of galaxies which have recently stopped star formation. Such a relation between $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ and $f_\mathrm{esc}^\mathrm{LyC}$ hints at the interesting role of supernovae, CRs, and magnetic fields in facilitating the escape ( and/or the lack) of LyC photons. (Abridged), Comment: 20 pages, 17 figures, 4 tables, 8 appendix pages, Astronomy & Astrophysics in press

The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar S\'ersic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications., Comment: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

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

In the current JWST era, rest-frame UV spectra play a crucial role in enhancing our understanding of the interstellar medium (ISM) and stellar properties of the first galaxies in the epoch of reionization (EoR, $z>6$). Here, we compare well-known and reliable optical diagrams sensitive to the main ionization source (i.e., star formation, SF; active galactic nuclei, AGN; shocks) to UV counterparts proposed in the literature - the so-called ``UV-BPT diagrams'' - using the HST COS Legacy Archive Spectroscopic SurveY (CLASSY), the largest high-quality, high-resolution and broad-wavelength range atlas of far-UV spectra for 45 local star-forming galaxies. In particular, we explore where CLASSY UV line ratios are located in the different UV diagnostic plots, taking into account state-of-the-art photoionization and shock models and, for the first time, the measured ISM and stellar properties (e.g., gas-phase metallicity, ionization parameter, carbon abundance, stellar age). We find that the combination of C III] $\lambda\lambda$1907,9 He II $\lambda1640$ and O III] $\lambda$1666 can be a powerful tool to separate between SF, shocks and AGN at sub-solar metallicities. We also confirm that alternative diagrams without O III] $\lambda$1666 still allow us to define a SF-locus with some caveats. Diagrams including C IV $\lambda\lambda$1548,51 should be taken with caution given the complexity of this doublet profile. Finally, we present a discussion detailing the ISM conditions required to detect UV emission lines, visible only in low gas-phase metallicity (12+log(O/H) $\lesssim8.3$) and high ionization parameter (log($U$) $\gtrsim-2.5$) environments. Overall, CLASSY and our UV toolkit will be crucial in interpreting the spectra of the earliest galaxies that JWST is currently revealing., Comment: 36 pages, 15 figures, accepted for publication in ApJ

We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$\alpha$-emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) survey. From NIRSpec data, we confirm the systemic redshifts of these Ly$\alpha$ emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal poor. We investigate Ly$\alpha$ line properties, including the line flux, velocity offset, and spatial extension. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements ($\sim5\times$ greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly$\alpha$ emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly$\alpha$ in one galaxy is significantly more extended than the non-resonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly$\alpha$ sources, discussing viable scenarios for the creation of sizable ionized bubbles ($>$1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, requiring additional sources of ionizing photons. Therefore, the fact that we detect Ly$\alpha$ from these galaxies suggests diverse scenarios on escape of Ly$\alpha$ during the epoch of reionization. High spectral resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization., Comment: Submitted to ApJ (18 pages, 7 figures, 2 tables)

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV<-20), with blue continua (-2.38 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8., Comment: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times., Comment: Accelerated Article Preview (AAP) version published in Nature

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

We report on the discovery of two low-luminosity, broad-line AGN at $z>5$ identified using JWST NIRSpec spectroscopy from the CEERS Survey. We detect broad H$\alpha$ emission from both sources, with FWHM of $2038\pm286$ and $1807\pm207$ km s$^{-1}$, resulting in black hole (BH) masses that are 1-2 dex below that of existing samples of luminous quasars at $z>5$. The first source, CEERS 1670 at $z=5.242$, is 2-3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its rest-frame optical SED. We measure a BH mass of $M_{\rm BH}=1.3\pm0.4\times 10^{7}~M_{\odot}$, confirming that this AGN is powered by the least-massive BH known in the universe at the end of cosmic reionization. The second source, CEERS 3210 at $z=5.624$, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be $M_{\rm BH}\simeq 0.9-4.7 \times 10^{7}~M_{\odot}$, depending on the level of dust obscuration assumed. We derive host stellar masses, $M_\star$, allowing us to place constraints on the BH-galaxy mass relationship in the lowest mass range yet probed in the early universe. The $M_{\rm BH}/M_\star$ ratio for CEERS 1670, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at $z=0$. We examine the emission-line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for low-metallicity AGN with $Z/Z_\odot \simeq 0.2-0.4$. The spectroscopic identification of low-luminosity, broad-line AGN at $z>5$ with $M_{\rm BH}\simeq 10^{7}~M_{\odot}$ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH-galaxy assembly., Comment: 16 pages, 8 figures, submitted to ApJL

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5\mu$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$\beta$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$\alpha$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$\beta$ equivalent widths ($\simeq$1100$\r{A}$), and elevated production efficiency of ionizing photons ($\log(\xi_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$\beta$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models., Comment: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

Spectroscopic studies of extreme-ionization galaxies (EIGs) are critical to our understanding of exotic systems throughout cosmic time. These EIGs exhibit spectral features requiring >54.42 eV photons: the energy needed to fully ionize helium into He2+ and emit He II recombination lines. They are likely key contributors to reionization, and they can also probe exotic stellar populations or accretion onto massive black holes. To facilitate the use of EIGs as probes of high ionization, we focus on ratios constructed from strong rest-frame UV/optical emission lines, specifically [O III] 5008, H-beta, [Ne III] 3870, [O II] 3727,3729, and [Ne V] 3427. These lines probe the relative intensity at energies of 35.12, 13.62, 40.96, 13.62 eV, and 97.12, respectively, covering a wider range of ionization than traced by other common rest-frame UV/optical techniques. We use ratios of these lines ([Ne V]/[Ne III] = Ne53 and [Ne III]/[O II]), which are closely separated in wavelength, and mitigates effects of dust attenuation and uncertainties in flux calibration. We make predictions from photoionization models constructed from Cloudy that use a broad range of stellar populations and black hole accretion models to explore the sensitivity of these line ratios to changes in the ionizing spectrum. We compare our models to observations from the Hubble Space Telescope and James Webb Space Telescope of galaxies with strong high-ionization emission lines at z ~ 0, z ~ 2, and z ~ 7. We show that the Ne53 ratio can separate galaxies with ionization from 'normal' stellar populations from those with AGN and even 'exotic' Population III models. We introduce new selection methods to identify galaxies with photoionization driven by Population III stars or intermediate-mass black hole accretion disks that could be identified in upcoming high-redshift spectroscopic surveys., Comment: 16 pages, 5 figures, 1 table. Accepted in 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 ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$\mu$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$~GHz exactly coinciding with the JWST source position, with a 2\% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11--13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11--17 is favored in every candidate as a very blue (UV continuum slope of $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$--5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in survey volumes used for the identification of the $z\sim$11--17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search., Comment: 25 pages, 8 figures, 7 tables. ApJ in press

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of S\'ersic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher S\'ersic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed., Comment: Accepted for publication to ApJL, 24 pages, 14 figures

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$\lambda$5008/H$\beta$ ratios, similar [OIII]$\lambda$4364/H$\gamma$ ratios, and higher ($\sim$0.5 dex) [NeIII]$\lambda$3870/[OII]$\lambda$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$\lambda$3870/[OII]$\lambda$3728, [OIII]$\lambda$4364/H$\gamma$, and [OIII]$\lambda$5008/H$\beta$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$\lambda$4364/H$\gamma$ and [OIII]$\lambda$(4960+5008)/H$\beta$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn., Comment: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

Strong nebular emission lines are an important diagnostic tool for tracing the evolution of star-forming galaxies across cosmic time. However, different observational setups can affect these lines, and the derivation of the physical nebular properties. We analyze 12 local star-forming galaxies from the COS Legacy Spectroscopy SurveY (CLASSY) to assess the impact of using different aperture combinations on the determination of the physical conditions and gas-phase metallicity. We compare optical spectra observed with the SDSS aperture, which has a 3" of diameter similar to COS, to IFU and longslit spectra, including new LBT/MODS observations of five CLASSY galaxies. We calculate the reddening, electron densities and temperatures, metallicities, star formation rates, and equivalent widths (EWs). We find that measurements of the electron densities and temperatures, and metallicity remained roughly constant with aperture size, indicating that the gas conditions are relatively uniform for this sample. However, using the IFU observations of 3 galaxies, we find that the E(B-V) values derived from the Balmer ratios decrease ( by up to 53%) with increasing aperture size. The values change most significantly in the center of the galaxies, and level out near the COS aperture diameter of 2.5". We examine the relative contributions from the gas and stars using the H$\alpha$ and [OIII] $\lambda$5007 EWs as a function of aperture light fraction, but find little to no variations within a given galaxy. These results imply that the optical spectra provide nebular properties appropriate for the FUV CLASSY spectra, even when narrow 1.0" long-slit observations are used., Comment: Accepted for publication in ApJ

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

Extreme emission line galaxies (EELGs) are a notable galaxy genus, ultimately being regarded as local prototypes of early galaxies at the cosmic noon. Robust characterisation of their stellar content, however, is hindered by the exceptionally high nebular emission present in their optical spectroscopic data. This study is dedicated into recovering the stellar properties of a sample of 414 EELGs as observed by the SDSS Survey. Such is achieved by means of the spectral synthesis code FADO, which self-consistently considers the stellar and nebular emission in an optical spectrum. Additionally, a comparative analysis was carried on, by further processing the EELGs sample with the purely stellar spectral synthesis code Starlight, and by extending the analysis to a sample of 697 normal star-forming galaxies, expected to be less affected by nebular contribution. We find that, for both galaxy samples, stellar mass and mean age estimates by Starlight are systematically biased towards higher values, and that an adequate determination of the physical and evolutionary properties of EELGs via spectral synthesis is only possible when nebular continuum emission is taken into account. Moreover, the differences between the two population synthesis codes can be ascribed to the degree of star-formation activity through the specific star-formation rate and the sum of the flux of the most prominent emission lines. As expected, on the basis of the theoretical framework, our results emphasise the importance of considering the nebular emission while performing spectral synthesis, even for galaxies hosting typical levels of star-formation activity., Comment: 14 pages, 13 figures, accepted for publication in A&A

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

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

We have investigated the oxygen and nitrogen chemical abundances in extremely compact star-forming galaxies with redshifts between $\sim$0.11-0.35, popularly referred to as "green peas". Direct and strong-line methods sensitive to the N/O ratio applied to their SDSS spectra reveals that these systems are genuine metal-poor galaxies, with mean oxygen abundances 20% solar. At a given metallicity these galaxies display systematically large N/O ratios compared to normal galaxies, which can explain the strong difference between our metallicities measurements and previous ones. While their N/O ratios follow the relation with stellar mass of local star-forming galaxies in the SDSS, we find that the mass--metallicity relation of the "green peas" is offset $\ga$0.3 dex to lower metallicities. We argue that recent interaction-induced inflow of gas, possibly coupled with a selective metal-rich gas loss, driven by supernova winds, may explain our findings and the known galaxy properties, namely high specific star formation rates, extreme compactness, and disturbed optical morphologies. The "green pea" galaxy properties seem to be not common in the nearby Universe, suggesting a short and extreme stage of their evolution. Therefore, these galaxies may allow us to study in great detail many processes, such as starburst activity and chemical enrichment, under physical conditions approaching those in galaxies at higher redshifts., Comment: 8 pages, 3 Figures, ApJ Letters in press. Proof version with corrections

The structural properties of the low surface brightness stellar host in blue compact dwarf galaxies are often studied by fitting r^{1/n} models to the outer regions of their radial profiles. The limitations imposed by the presence of a large starburst emission overlapping the underlying component makes this kind of analysis a difficult task. We propose a two-dimensional fitting methodology in order to improve the extraction of the structural parameters of the LSB host. We discuss its advantages and weaknesses by using a set of simulated galaxies and compare the results for a sample of eight objects with those already obtained using a one-dimensional technique. We fit a PSF convolved Sersic model to synthetic galaxies, and to real galaxy images in the B, V, R filters. We restrict the fit to the stellar host by masking out the starburst region and take special care to minimize the sky-subtraction uncertainties. In order to test the robustness and flexibility of the method, we carry out a set of fits with synthetic galaxies. Furthermore consistency checks are performed to assess the reliability and accuracy of the derived structural parameters. The more accurate isolation of the starburst emission is the most important advantage and strength of the method. Thus, we fit the host galaxy in a range of surface brightness and in a portion of area larger than in previous published 1D fits with the same dataset. We obtain robust fits for all the sample galaxies, all of which, except one, show Sersic indices n very close to 1, with good agreement in the three bands. These findings suggest that the stellar hosts in BCDs have near-exponential profiles, a result that will help us to understand the mechanisms that form and shape BCD galaxies, and how they relate to the other dwarf galaxy classes., Comment: 22 pages, 15 figures (low resolution), accepted for publication in A&A. A higher resolution version of the figures can be provided upon request