Search

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

Isolated, low-mass galaxies provide the opportunity to assess the impact of reionization on their star formation histories (SFHs) without the ambiguity of environmental processes associated with massive host galaxies. There are very few isolated, low-mass galaxies that are close enough to determine their SFHs from resolved star photometry reaching below the oldest main sequence turnoff. JWST has increased the volume for which this is possible, and here we report on JWST observations of the low-mass, isolated galaxy Leo P. From NIRCam imaging in F090W, F150W, and F277W, we derive a SFH which shows early star formation followed by a pause subsequent to the epoch of reionization which is then later followed by a re-ignition of star formation. This is very similar to the SFHs from previous studies of other dwarf galaxies in the ``transition zone'' between quenched very low-mass galaxies and the more massive galaxies which show no evidence of the impact of reionization on their SFHs; this pattern is rarely produced in simulations of SFHs. The lifetime SFH reveals that Leo P's stellar mass at the epoch of reionization was in the range that is normally associated with being totally quenched. The extended pause in star formation from z~5-1 has important implications for the contribution of low-mass galaxies to the UV photon budget at intermediate redshifts. We also demonstrate that, due to higher sensitivity and angular resolution, observing in two NIRCam short wavelength filters is superior to observing in a combination of a short and a long wavelength filter., Comment: 24 pages, 9 figures, 3 tables

We present $\sim300$ stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI ($M_V = -7.5$) and Andromeda XXVIII ($M_V = -8.8$) derived using metallicity-sensitive Calcium H & K narrow-band Hubble Space Telescope imaging. These are the first individual stellar metallicities in And~XVI (95 stars). Our And~XXVIII sample (191 stars) is a factor of $\sim15$ increase over literature metallicities. For And~XVI, we measure $\langle \mbox{[Fe/H]}\rangle = -2.17^{+0.05}_{-0.05}$, $\sigma_{\mbox{[Fe/H]}}=0.33^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.23\pm0.15$ dex $R_e^{-1}$. We find that And XVI is more metal-rich than MW UFDs of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure $\langle \mbox{[Fe/H]}\rangle = -1.95^{+0.04}_{-0.04}$, $\sigma_{\mbox{[Fe/H]}}=0.34^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.46 \pm 0.10$~dex~$R_e^{-1}$, placing it on the dwarf galaxy mass-metallicity relation. Neither galaxy has a metallicity distribution function with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint ($L \lesssim 10^6~L_{\odot}$) galaxies outside the Milky Way halo. Both galaxies' gradients are consistent with predictions from the FIRE simulations, where an age-gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal poor ($\mbox{[Fe/H]} < -3.0$) star candidates, which make up 7% of And~XVI's MDF and 6% of And~XXVIII's., Comment: 17 pages, 5 figures, 4 tables, ApJ submitted; comments welcome!

We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $\Gamma$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overline{\Gamma} = 1.49\pm0.18$, with an intrinsic scatter of $\sigma^{2}_{\Gamma} = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overline{\Gamma}$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($\Gamma = 1.30$) IMF, as well as the Salpeter ($\Gamma = 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields., Comment: Accepted for publication in ApJ. 9 Figures, 1 Table

The tip of the red giant (TRGB) is a standardizable candle and is identifiable as the discontinuity at the bright extreme of the red giant branch (RGB) stars in color-magnitude diagram (CMD) space. The TRGB-based distance method has been calibrated and used to measured distances to galaxies out to $D\leq20$ Mpc with the $I$-band equivalent Hubble Space Telescope ($HST$) $F814W$ filter, and as an important rung in the distance ladder to measure the Hubble constant, $H_0$. In the infrared (IR), the TRGB apparent magnitude ranges from $1-2$ magnitudes brighter than in the optical, and now with the IR James Webb Space Telescope ($JWST$) observatory the feasible distance range of the TRGB method can be extended to $\sim50$ Mpc. However, in the IR the TRGB luminosity depends to varying degrees on stellar metallicity and age. In this study we standardize the TRGB luminosity using stellar colors as a proxy for metallicity/age to derive color-based corrections for the $JWST$ Near-Infrared Camera (NIRCam) short wavelength (SW) filters $F090W$, $F115W$, $F150W$ and the long wavelength (LW) filters $F277W$, $F356W,$ and $F444W$. We provide recommended filter combinations for distance measurements depending on the requisite precision. For science requiring high precision ($\leq1\%$ in distance) and robustness we recommend measuring the TRGB in $F090W$ vs $F090W-F150W$ or $F115W$ vs. $F115W-F277W$ with the caveat that even with $JWST$ long integration times will be necessary at further distances. If lower precision ($>1.5\%$ in distance) can be tolerated, or if shorter integration times are desirable, we recommend measuring the TRGB in either $F115W$ or $F150W$ paired with $F356W$. We do not recommend $F444W$ for precision TRGB measurements due to its lower angular resolution., Comment: 25 pages, 11 figures

We investigate the stellar and nebular properties of 9 H II regions in the spiral galaxy M101 with far-ultraviolet (FUV; ~900-2000 \r{A}) and optical (~3200-10000 \r{A}) spectra. We detect significant C III] 1907,1909 nebular emission in 7 regions, but O III] 1666 only in the lowest-metallicity region. We produce new analytic functions of the carbon ICF as a function of metallicity in order to perform a preliminary C/O abundance analysis. The FUV spectra also contain numerous stellar emission and P-Cygni features that we fit with luminosity-weighted combinations of single-burst Starburst99 and BPASS models. We find that the best-fit Starburst99 models closely match the observed very-high-ionization P-Cygni features, requiring very-hot, young (~< 3 Myr), metal-enriched massive stars. The youngest stellar populations are strongly correlated with broad He II emission, nitrogen Wolf-Rayet (WR) FUV and optical spectral features, and enhanced N/O gas abundances. Thus, the short-lived WR phase may be driving excess emission in several N P-Cygni wind features (955 \r{A}, 991 \r{A}, 1720 \r{A}) that bias the stellar continuum fits to higher metallicities relative to the gas-phase metallicities. Accurate characterization of these H II regions requires additional inclusion of WR stars in the stellar population synthesis models. Our FUV spectra demonstrate that the ~900-1200 \r{A} FUV can provide a strong test-bed for future WR atmosphere and evolution models., Comment: 24 ages, 12 figures, accepted to ApJ

We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $\mu$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $\mu$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies., Comment: 23 pages, 13 figures, accepted for publication in ApJ

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

The Tip of the Red Giant Branch (TRGB)-based distance method in the I band is one of the most efficient and precise techniques for measuring distances to nearby galaxies (D <= 15 Mpc). The TRGB in the near infrared (NIR) is 1 to 2 magnitudes brighter relative to the I band, and has the potential to expand the range over which distance measurements to nearby galaxies are feasible. Using Hubble Space Telescope (HST) imaging of 12 fields in 8 nearby galaxies, we determine color-based corrections and zero points of the TRGB in the Wide Field Camera 3 IR (WFC3/IR) F110W and F160W filters. First, we measure TRGB distances in the I band equivalent Advanced Camera System (ACS) F814W filter from resolved stellar populations with the HST. The TRGB in the ACS F814W filter is used for our distance anchor and to place the WFC3/IR magnitudes on an absolute scale. We then determine the color dependence (a proxy for metallicity/age) and zero point of the NIR TRGB from photometry of WFC3/IR fields which overlap with the ACS fields. The new calibration is accurate to ~1% in distance, relative to the F814W TRGB. Validating the accuracy of the calibrations, we find that the distance modulus for each field using the NIR TRGB calibration agrees with the distance modulus of the same fields as determined from the F814W TRGB. This is a JWST preparatory program and the work done here will directly inform our approach to calibrating the TRGB in JWST NIRCam and NIRISS photometric filters., Comment: 23 pages, 8 figures, accepted by the Astrophysical Journal (ApJ)

We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data., Comment: 30 pages, 17 figures. Accepted to ApJS. Data products to be hosted on MAST. For DOLPHOT/JWST tutorials, see https://dolphot-jwst.readthedocs.io/en/latest/ . For more program and DOLPHOT info, see https://ers-stars.github.io

We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas., Comment: Submitted to The Astrophysical Journal

We present an investigation of evolved stars in the nearby star-forming galaxy WLM, using NIRCam imaging from the JWST resolved stellar populations early-release science (ERS) program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants (RSGs) and thermally-pulsing asymptotic giant branch (TP-AGB) stars from one another, while also providing a reasonable separation of the primary TP-AGB subtypes: carbon-rich C-type stars and oxygen-rich M-type stars. The classification scheme we present here agrees very well with the well-established Hubble Space Telescope (HST) medium-band filter technique. The ratio of C to M-type stars (C/M) is 0.8$\pm$0.1 for both the new JWST and the HST classifications, which is within one sigma of empirical predictions from optical narrow-band CN and TiO filters. The evolved star colors show good agreement with the predictions from the PARSEC$+$COLIBRI stellar evolutionary models, and the models indicate a strong metallicity dependence that makes stellar identification even more effective at higher metallicity. However, the models also indicate that evolved star identification with NIRCam may be more difficult at lower metallicies. We test every combination of NIRCam filters using the models and present additional filters that are also useful for evolved star studies. We also find that $\approx$90\% of the dusty evolved stars are carbon-rich, suggesting that carbonaceous dust dominates the present-day dust production in WLM, similar to the findings in the Magellanic Clouds. These results demonstrate the usefulness of NIRCam in identifying and classifying dust-producing stars without the need for mid-infrared data., Comment: 15 page, 12 figures, submitted to AAS Journals

We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering a SFH and AMR with the NIRCam F090W-F150W filter combination and provide validation of the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation re-ignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo~A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ~5 Gyr; the stellar mass then rapidly doubled in ~2.5 Gyr, followed by constant star formation over the last ~5 Gyr., Comment: 27 pages, 10 figures, 1 table

Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies., Comment: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use them

We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$ completely distinct from what is seen in broadband optical imaging, and the first-ever confirmation of a bar in M33. We estimate a bar extent of $\sim$1 kpc. The two spiral arms are asymmetric in orientation and strength, and likely represent the innermost impact of the recent tidal interaction responsible for M33's warp at larger scales. The flocculent multi-armed morphology for which M33 is known is only visible in the young upper main sequence population, which closely tracks the morphology of the ISM. We investigate the stability of M33's disk, finding $Q{\sim}1$ over the majority of the disk. We fit multiple components to the old stellar density distribution and find that, when considering recent stellar kinematics, M33's bulk structure favors the inclusion of an accreted halo component, modeled as a broken power-law. The best-fit halo model has an outer power-law index of $-$3 and accurately describes observational evidence of M33's stellar halo from both resolved stellar spectroscopy in the disk and its stellar populations at large radius. Integrating this profile yields a total halo stellar mass of ${\sim}5{\times}10^8\ M_{\odot}$, giving a total stellar halo mass fraction of 16%, most of which resides in the innermost 2.5 kpc., Comment: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

Studying the galaxies responsible for reionization is often conducted through local reionization-era analogs; however, many of these local analogs are too massive to be representative of the low-mass star-forming galaxies that are thought to play a dominant role in reionization. The local, low-mass dwarf starburst galaxy Pox 186 is one such system with physical conditions representative of a reionization-era starburst galaxy. We present deep ultraviolet (UV) spectroscopy of Pox 186 to study its stellar population and ionization conditions and to compare these conditions to other local starburst galaxies. The new Cosmic Origins Spectrograph data are combined with archival observations to cover $\sim$1150-2000 A and allow for an assessment of Pox 186's stellar population, the relative enrichment of C and O, and the escape of ionizing photons. We detect significant Ly$\alpha$ and low-ionization state absorption features, indicative of previously undetected neutral gas in Pox 186. The C/O relative abundance, log(C/O) = -0.62$\pm$0.02, is consistent with other low-metallicity dwarf galaxies and suggests a comparable star formation history in these systems. We compare UV line ratios in Pox 186 to those of dwarf galaxies and photoionization models, and we find excellent agreement for the ratios utilizing the intense C III], O III], and double-peaked C IV lines. However, the UV and optical He II emission is faint and distinguishes Pox 186 from other local starburst dwarf galaxies. We explore mechanisms that could produce faint He II, which have implications for the low-mass reionization-era galaxies which may have similar ionization conditions., Comment: 22 pages, 9 figures, accepted for publication in The Astrophysical Journal

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

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 provide one of the most comprehensive metallicity studies at z $\sim$ 4 by analyzing the UV/optical HST photometry, and rest-frame VLT-FORS2 ultraviolet and VLT-XSHOOTER optical spectra of J0332-3557, a gravitationally lensed galaxy magnified by a factor of 20. With a 5$\sigma$ detection of the auroral O III]1666 line, we are able to derive a direct gas metallicity estimate for our target. We find Z_gas = 12 + log(O/H) = 8.26 $\pm$ 0.06, which is compatible with an increasing of both the gas fraction and the outflow metal loading factor from z $\sim$ 0 to z $\sim$ 4. J0332-3557 is the most metal-rich individual galaxy at z $\sim$ 4 for which the C/O ratio has been measured. We derive a low log(C/O) = - 1.02 $\pm$ 0.2, which suggests that J0332-3557 is in the early stages of ISM carbon enrichment driven mostly by massive stars. The low C/O abundance also indicates that J0332-3557 is characterized by a low star formation efficiency, higher yields of oxygen, and longer burst duration. We find that EW(C III])1906,9 is as low as $\sim$ 3 \AA. The main drivers of the low EW(C III])1906,9 are the higher gas metallicity and the low C/O abundance. J0332-3557 is characterized by one diffuse and two more compact regions $\sim$ 1 kpc in size. We find that the carbon emission mostly originates in the compact knots. Our study on J0332-3557 serves as an anchor for studies investigating the evolution of metallicity and C/O abundance across different redshifts., Comment: 27 pages, 17 figures. Accepted for publication in ApJ

We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM., Comment: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come., Comment: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

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 dispersion in chemical abundances provides a very strong constraint on the processes that drive the chemical enrichment of galaxies. Due to its proximity, the spiral galaxy M33 has been the focus of numerous chemical abundance surveys to study the chemical enrichment and dispersion in abundances over large spatial scales. The CHemical Abundances Of Spirals (CHAOS) project has observed $\sim$100 H II regions in M33 with the Large Binocular Telescope (LBT), producing the largest homogeneous sample of electron temperatures (T$_e$) and direct abundances in this galaxy. Our LBT observations produce a robust oxygen abundance gradient of $-$0.037 $\pm$ 0.007 dex/kpc and indicate a relatively small (0.043 $\pm$ 0.015 dex) intrinsic dispersion in oxygen abundance relative to this gradient. The dispersions in N/H and N/O are similarly small and the abundances of Ne, S, Cl, and Ar relative to O are consistent with the solar ratio as expected for $\alpha$-process or $\alpha$-process-dependent elements. Taken together, the ISM in M33 is chemically well-mixed and homogeneously enriched from inside-out with no evidence of significant abundance variations at a given radius in the galaxy. Our results are compared to those of the numerous studies in the literature, and we discuss possible contaminating sources that can inflate abundance dispersion measurements. Importantly, if abundances are derived from a single T$_e$ measurement and T$_e$-T$_e$ relationships are relied on for inferring the temperature in the unmeasured ionization zone, this can lead to systematic biases which increase the measured dispersion up to 0.11 dex., Comment: 52 pages, 13 figures, accepted for publication in The Astrophysical Journal

We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improved zeropoints for the ERS filters and show that alternate zeropoints derived by the community also improve the calibration significantly. We also find that the detector offsets appear to be time variable by up to at least 0.1 mag., Comment: Accepted for publication in RNAAS

We analyze the rest-frame near-UV and optical nebular spectra of three $z > 7$ galaxies from the Early Release Observations taken with the Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST). These three high-z galaxies show the detection of several strong-emission nebular lines, including the temperature-sensitive [O III] $\lambda$4363 line, allowing us to directly determine the nebular conditions and abundances for O/H, C/O, and Ne/O. We derive O/H abundances and ionization parameters that are generally consistent with other recent analyses. We analyze the mass-metallicity relationship (i.e., slope) and its redshift evolution by comparing between the three z > 7 galaxies and local star-forming galaxies. We also detect the C III] $\lambda\lambda$1907,1909 emission in a z > 8 galaxy from which we determine the most distant C/O abundance to date. This valuable detection of log(C/O) = $-0.83\pm0.38$ provides the first test of C/O redshift evolution out to high-redshift. For neon, we use the high-ionization [Ne III] $\lambda$3869 line to measure the first Ne/O abundances at z>7, finding no evolution in this $\alpha$-element ratio. We explore the tentative detection of [Fe II] and [Fe III] lines in a z>8 galaxy, which would indicate a rapid build up of metals. Importantly, we demonstrate that properly flux-calibrated and higher S/N spectra are crucial to robustly determine the abundance pattern in z>7 galaxies with NIRSpec/JWST., Comment: Accepted for publication on ApJL

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

We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system., Comment: 32 pages, 15 figures, accepted for publication on 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

The ratio of baryonic-to-dark matter in present-day galaxies constrains galaxy formation theories and can be determined empirically via the baryonic Tully-Fisher relation (BTFR), which compares a galaxy's baryonic mass (Mbary) to its maximum rotation velocity (Vmax). The BTFR is well-determined at Mbary >10^8 Msun, but poorly constrained at lower masses due to small samples and the challenges of measuring rotation velocities in this regime. For 25 galaxies with high-quality data and Mbary <~10^8 Msun, we estimate Mbary from infrared and HI observations and Vmax from the HI gas rotation. Many of the Vmax values are lower limits because the velocities are still rising at the edge of the detected HI disks (Rmax); consequently, most of our sample has lower velocities than expected from extrapolations of the BTFR at higher masses. To estimate Vmax, we map each galaxy to a dark matter halo assuming density profiles with and without cores. In contrast to non-cored profiles, we find the cored profile rotation curves are still rising at Rmax values, similar to the data. When we compare the Vmax values derived from the cored density profiles to our Mbary measurements, we find a turndown of the BTFR at low masses that is consistent with LCDM predictions and implying baryon fractions of 1-10% of the cosmic value. Although we are limited by the sample size and assumptions inherent in mapping measured rotational velocities to theoretical rotation curves, our results suggest that galaxy formation efficiency drops at masses below Mbary~10^8 Msun, corresponding to Mhalo~10^10 Msun., Comment: 39 pages, 4 tables, 24 figures

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