Search

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

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

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

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

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

The 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 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

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

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

We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products., Comment: Accepted for publication in MNRAS (11 pages, 8 figures)

We present optical spectra of nine Galactic H II regions observed with the 10.4 m Gran Telescopio Canarias telescope and located at Galactocentric distances (RG) from 4 to 8 kpc. The distances of the objects have been revised using Gaia DR2 parallaxes. We determine the electron temperature for all the nebulae, which allows a precise computation of their ionic abundances. We have included published data of an additional sample of Galactic H II regions, providing a final data set of 42 objects. The shape of the radial gradients of O/H and N/H is linear and constant, discarding any substantial change of the slope, at least for RG between 4 and 17 kpc. The small dispersion of the O/H and N/H values with respect to the computed gradients implies the absence of significant azimuthal variations of the chemical abundances, at least in the quadrant covered by our observations. We find an almost flat N/O versus O/H diagram relation. This result is not observed in other nearby spiral galaxy except M31. Finally, we compare our computed gradients with those obtained using far infrared (FIR) spectra. We confirm the significant offset in the N/O distribution between the optical and FIR observations. Possible explanations involve ionization correction factors and the strong dependence on density of the abundance determinations based on FIR lines., Comment: 18 pages, 11 figures, 5 tables. Accepted for publication in MNRAS

We present a reassessment of the radial abundance gradients of C, N, O, Ne, S, Cl and Ar in the Milky Way using deep spectra of 33 H II regions gathered from the literature, covering Galactocentric distances from 6 to 17 kpc. The distances of the objects have been revised using Gaia parallaxes. We recalculate the physical conditions and ionic abundances in an homogeneous way using updated atomic data. All the objects have direct determination of the electron temperature, permitting to derive their precise ionic abundances. We analyze and compare different ICF schemes for each element in order to obtain the most confident total abundances. Due to the revised distances, our results do not support previous claims about a possible flattening of the O/H gradient in the inner Galactic disk. We find that the Galactic N/O gradient is rather flat, in contrast to what has been found in other spiral galaxies. The slope of the gradients of some elements is sensitive to the ICF scheme used, especially in the case of Ne. The dispersion around the fit for the gradients of C, N, O, S, Cl and Ar is of the order of the typical uncertainties in the determination of the abundances, implying the absence of significant inhomogeneities in the chemical composition of the ionized gas phase of the ISM. We find flat gradients of log(S/O) and log(Cl/O) and very shallow or flat ones for log(Ne/O) and log(Ar/O), consistent with a lockstep evolution of Ne, S, Cl and Ar with respect to O., Comment: 26 pages, 12 figures. Accepted for publication in MNRAS