Search

We present the Cosmic Evolution Early Release Science (CEERS) Survey, a 77.2 hour Director's Discretionary Early Release Science Program. CEERS demonstrates, tests, and validates efficient extragalactic surveys using coordinated, overlapping parallel observations with the JWST instrument suite, including NIRCam and MIRI imaging, NIRSpec low (R~100) and medium (R~1000) resolution spectroscopy, and NIRCam slitless grism (R~1500) spectroscopy. CEERS targets the Hubble Space Telescope-observed region of the Extended Groth Strip (EGS) field, supported by a rich set of multiwavelength data. CEERS facilitated immediate community science in both of the extragalactic core JWST science drivers ``First Light" and ``Galaxy Assembly," including: 1) The discovery and characterization of large samples of galaxies at z >~ 10 from ~90 arcmin^2 of NIRCam imaging, constraining their abundance and physical nature; 2) Deep spectra of >1000 galaxies, including dozens of galaxies at 63; and 4) Characterizing galaxy mid-IR emission with MIRI to study dust-obscured star-formation and supermassive black hole growth at z~1-3. As a legacy product for the community, the CEERS team has provided several data releases, accompanied by detailed notes on the data reduction procedures and notebooks to aid in reproducibility. In addition to an overview of the survey and quality of the data, we provide science highlights from the first two years with CEERS data., Comment: 38 pages, 13 figures, 6 tables

Gas metallicity, ionization parameter, and gas pressure can affect the observed ratios of specific strong emission lines within galaxies. While the theoretical strong lines diagnostics for gas metallicity, ionization parameters, and gas pressure in star-forming regions are well-established, theoretical diagnostics for active galactic nuclei (AGNs) narrow line regions are still lacking. In Zhu et al. (2023), we presented a new AGN model that provides the best predictions for observations spanning the UV, optical, and infrared wavelengths. This paper presents a suite of theoretical diagnostics for the gas metallicity, ionization parameter, gas pressure, and the peak energy in AGN ionizing radiation field $E_{peak}$ for AGN narrow-line regions spanning the UV and optical wavelengths. We investigate the model dependency on the ionization parameter, gas pressure, $E_{peak}$, and the nitrogen scaling relation and make recommendations on metallicity diagnostics that are most robust against these parameters. We test our new AGN metallicity diagnostics using optical galaxy spectra from Sloan Digital Sky Survey DR16. These tests show that the metallicities measured from different diagnostics in this paper are consistent within $\sim0.3$ dex. We compare consistent HII and AGN diagnostics and demonstrate that HII and AGN diagnostics should not be used interchangeably. With a wide wavelength coverage, we anticipate that these AGN diagnostics will enable new metallicity studies of galaxies dominated by AGN., Comment: 33 pages, 22 figures, 7 tables, Accepted for publication in ApJ

We investigate the coevolution of metals and dust for 173 galaxies at 4.010 by JWST. Besides, we observe that the metallicity of galaxies at z>8 presents a metal-to-stellar mass ratio larger than a few 10^-3, above a floor. This suggests a very fast rise of metals at high redshift, impacting the tentative detections of population III objects., Comment: This paper is submitted. It contains a main paper with 5 figures and 2 tables, plus supplementary materials

We present the spatially resolved relationship between the dust-to-gas mass ratio (DGR) and gas-phase metallicity (Zgas or 12+log(O/H)) (i.e., DGR-Zgas relation) of 11 nearby galaxies with a large metallicity range (1.5 dex of 12+log(O/H)) at (sub-)kpc scales. We used the large field-of-view (> 3') optical pseudo-Integral Field Spectroscopy data taken by the TYPHOON/PrISM survey, covering the optical size of galaxies, combining them with multi-wavelength data (far-UV to far-IR, CO, and HI 21 cm radio). A large scatter of DGR in the intermediate metallicity galaxies (8.0 < 12+log(O/H) < 8.3) is found, which is in line with dust evolution models, where grain growth begins to dominate the mechanism of dust mass accumulation. In the lowest metallicity galaxy of our sample, Sextans A (12+log(O/H) < 7.6), the star-forming regions have significantly higher DGR values (by 0.5-2 dex) than the global estimates from literature at the same metallicity but aligns with the DGR values from metal depletion method from Damped Lyman Alpha systems and high hydrogen gas density regions of Sextans A. Using dust evolution models with a Bayesian MCMC approach suggests: 1) a high SN dust yield and 2) a negligible amount of photofragmentation by UV radiation, although we note that our sample in the low-metallicity regime is limited to Sextans A. On the other hand, it is also possible that while metallicity influences DGR, gas density also plays a role, indicating an early onset of dust grain growth in the dust mass build-up process despite its low metallicity., Comment: 20 pages, 10 figures, 7 tables, accepted for publication in MNRAS

Most star formation in the local Universe occurs in spiral galaxies, but their origin remains an unanswered question. Various theories have been proposed to explain the development of spiral arms, each predicting different spatial distributions of the interstellar medium. This study maps the star formation rate (SFR) and gas-phase metallicity of nine spiral galaxies with the TYPHOON survey to test two dominating theories: density wave theory and dynamic spiral theory. We discuss the environmental effects on our galaxies, considering reported environments and merging events. Taking advantage of the large field of view covering the entire optical disk, we quantify the fluctuation of SFR and metallicity relative to the azimuthal distance from the spiral arms. We find higher SFR and metallicity in the trailing edge of NGC~1365 (by 0.117~dex and 0.068~dex, respectively) and NGC~1566 (by 0.119~dex and 0.037~dex, respectively), which is in line with density wave theory. NGC~2442 shows a different result with higher metallicity (0.093~dex) in the leading edge, possibly attributed to an ongoing merging. The other six spiral galaxies show no statistically significant offset in SFR or metallicity, consistent with dynamic spiral theory. We also compare the behaviour of metallicity inside and outside the co-rotation radius (CR) of NGC~1365 and NGC~1566. We find comparable metallicity fluctuations near and beyond the CR of NGC~1365, indicating gravitational perturbation. NGC~1566 shows the greatest fluctuation near the CR, in line with the analytic spiral arms. Our work highlights that a combination of mechanisms explains the origin of spiral features in the local Universe., Comment: 22 pages, 10 figures, 4 tables, accepted for publication in MNRAS

A crucial yet challenging task in galaxy evolution studies is the identification of distant merging galaxies, a task which suffers from a variety of issues ranging from telescope sensitivities and limitations to the inherently chaotic morphologies of young galaxies. In this paper, we use random forests and convolutional neural networks to identify high-redshift JWST CEERS galaxy mergers. We train these algorithms on simulated $3

The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent high-$z$ JWST data indicate deviations from the FMR established at low-$z$. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a $z = 0$ calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a ``static'' FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a ``dynamic'' FMR, where the role of SFR evolves over time. We find static FMRs in SIMBA and dynamic FMRs in Illustris, IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at $z > 4$ by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may have a similar dynamic trend as these simulations. Overall, our findings imply that the current FMR framework neglects important time variations of these simulations' baryon cycles., Comment: Accepted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856)

Nebular spectral lines provide insight into the properties of the interstellar medium and the ionizing radiation within galaxies. The presence of high-energy ionization lines, such as He II, indicates the existence of hard ionizing photons, because the second ionization energy of helium is high (54 eV). The enigma surrounding the origin of these lines observed in star-forming galaxies persists, as stellar ionization cannot account for such high energy emission. This paper proposes that supersoft X-ray sources (SSS) may produce these high energy ionization lines in star-forming galaxies. We model the spectra of such sources using blackbody radiation and then add the blackbody to stellar population spectra to represent the overall spectra of galaxies. We then use a photoionization model to predict the resulting emission lines and compare the contribution of SSS to the observation of highly ionized lines in star-forming galaxies, both at low and high redshifts. We find that incorporating a blackbody with temperatures between kT = 20-100eV can boost the He II emission line ratio to the observed level. The blackbody temperature range aligns with the observed temperatures of the SSSs. The number of SSSs in spiral galaxies listed in Chandra catalogues, and our estimates of the total population, confirms that SSSs are promising candidates for the source of the He II ionization., Comment: 17 pages, 7 figures, comments are welcomed

Galaxies' stellar masses, gas-phase oxygen abundances (metallicity), and star formation rates (SFRs) obey a series of empirical correlations, most notably the mass-metallicity relation (MZR) and fundamental metallicity relation (FZR), which relates oxygen abundance to a combination of stellar mass and SFR. However, due to the difficulty of measuring oxygen abundances and SFRs in galaxies that host powerful active galactic nuclei (AGN), to date it is unknown to what extent AGN-host galaxies also follow these correlations. In this work, we apply Bayesian methods to the MaNGA integral field spectrographic (IFS) survey that allow us to measure oxygen abundances and SFRs in AGN hosts, and use these measurements to explore how the MZR and FZR differ between galaxies that do and do not host AGN. We find similar MZRs at stellar masses above $10^{10.5} \mathrm{M}_\odot$, but that at lower stellar masses AGN hosts show up to $\sim 0.2$ dex higher oxygen abundances. The offset in the FZR is significantly smaller, suggesting that the larger deviation in the MZR is a result of AGN-host galaxies having systematically lower SFRs at fixed stellar mass. However, within the AGN-host sample there is little correlation between SFR and oxygen abundance. These findings support a scenario in which an AGN can halt efficient gas accretion, which drives non-AGN host galaxies to both higher SFR and lower oxygen abundance, resulting in the galaxy evolving off the star-forming main sequence (SFMS). As a consequence, as the SFR declines for an individual system its metallicity remains mostly unchanged., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i.e., redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and EAGLE all predict MZRs that exhibit scatter with a secondary dependence on SFR up to $z=8$. We introduce the concept of a "strong" FMR, where the strength of correlated scatter does not evolve with time, and a "weak" FMR, where there is some time evolution. We find that each simulation analysed has a weak FMR -- there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ~10-40% at $z\gtrsim3$ when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR., Comment: 10 pages, 5 figures (Figure 5 is summary). Accepted to MNRAS

The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (${\rm MZ}_*{\rm R}$) in the hydrodynamic simulations Illustris, TNG, and EAGLE to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the ${\rm MZ}_*{\rm R}$ across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g., SFR, gas fraction, atomic gas mass). Since we find this correlation exists in all three simulations, each employing a sub-grid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ``burstiness'' of feedback within galaxies., Comment: 15 pages, 9 figures, Accepted to MNRAS

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$\sigma$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$\beta$), and dust content, of galaxies at $1

Galactic fountains driven by star formation result in a variety of kinematic structures such as ionised winds and thick gas disks, both of which manifest as complex emission line profiles that can be parametrised by multiple Gaussian components. We use integral field spectroscopy (IFS) from the SAMI Galaxy Survey to spectrally resolve these features, traced by broad H$\alpha$ components, and distinguish them from the star-forming thin disk, traced by narrow components, in 3068 galaxies in the local Universe. Using a matched sample analysis technique, we demonstrate that the presence of complex emission line profiles in star-forming galaxies is most strongly correlated with the global star formation rate (SFR) surface density of the host galaxy measured within $1R_{\rm e}$ ($\Sigma_{{\rm SFR},R_{\rm e}}$), even when controlling for both observational biases, including inclination, amplitude-to-noise and angular scale, and sample biases in parameters such as stellar mass and SFR. Leveraging the spatially resolved nature of the dataset, we determine that the presence of complex emission line profiles within individual spaxels is driven not only by the local $\Sigma_{\rm SFR}$, but by the $\Sigma_{{\rm SFR},R_{\rm e}}$ of the host galaxy. We also parametrise the clumpiness of the SFR within individual galaxies, and find that $\Sigma_{{\rm SFR},R_{\rm e}}$ is a stronger predictor of the presence of complex emission line profiles than clumpiness. We conclude that, with a careful treatment of observational effects, it is possible to identify structures traced by complex emission line profiles, including winds and thick ionised gas disks, at the spatial and spectral resolution of SAMI using the Gaussian decomposition technique., Comment: 21 pages, 17 figures, accepted for publication in MNRAS

Spiral structures are important drivers of the secular evolution of disc galaxies, however, the origin of spiral arms and their effects on the development of galaxies remain mysterious. In this work, we present two three-armed spiral galaxies at z~0.3 in the Middle Age Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey. Taking advantage of the high spatial resolution (~0.6'') of the Multi-Unit Spectroscopic Unit (MUSE), we investigate the two-dimensional distributions of different spectral parameters: Halpha, gas-phase metallicity, and D4000. We notice significant offsets in Halpha (~0.2 dex) as well as gas-phase metallicities (~0.05 dex) among the spiral arms, downstream and upstream of MAGPI1202197197 (SG1202). This observational signature suggests the spiral structure in SG1202 is consistent with arising from density wave theory. No azimuthal variation in Halpha or gas-phase metallicities is observed in MAGPI1204198199 (SG1204), which can be attributed to the tighter spiral arms in SG1204 than SG1202, coming with stronger mixing effects in the disc. The absence of azimuthal D4000 variation in both galaxies suggests the stars at different ages are well-mixed between the spiral arms and distributed around the disc regions. The different azimuthal distributions in Halpha and D4000 highlight the importance of time scales traced by various spectral parameters when studying 2D distributions in spiral galaxies. This work demonstrates the feasibility of constraining spiral structures by tracing interstellar medium (ISM) and stellar population at z~0.3, with a plan to expand the study to the full MAGPI survey., Comment: 15 pages, 11 figures, 2 tables, accepted for publication in MNRAS

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

Photoionization models frequently assume constant temperature or density within HII regions. We investigate this assumption by measuring the detailed temperature and density structures of four HII regions in the Large Magellanic Cloud and the Small Magellanic Cloud, using integral-field spectroscopic data from the Wide-Field Spectrograph on the ANU 2.3m telescope. We analyse the distribution of emission-lines of low-ionization species, intermediate-ionization species and high-ionization species. We present the complex electron temperature and density structures within HII regions. All four nebulae present a negative gradient in the electron density profile. Both positive and negative temperature gradients are observed in the nebulae. We create a series of nebula models with a constant ISM pressure and varying temperature and density distributions. Comparison of the line ratios between our HII regions and models suggests that none of the simple nebula models can reproduce the observed temperature and density structures. Comparison between the models and the data suggests that the ISM pressure of nebulae in LMC and SMC is between log(P/k)=6-7.5. Complex internal structures of the nebulae highlight the importance of future Monte-Carlo photoionization codes for accurate nebula modeling, which include a comprehensive consideration of arbitrary geometries of HII regions., Comment: 24 pages, 12 figures, 11 tables, Accepted for publication in ApJ

The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24$\mu$m population in the EGS field. We find that MIRI can observe an order of magnitude deeper than MIPS in significantly shorter integration times, attributable to JWST's much larger aperture and MIRI's improved sensitivity. MIRI is exceptionally good at finding faint ($L_{\rm IR}<10^{10} L_\odot$) galaxies at $z\sim1-2$. We find that a significant portion of MIRI galaxies are "mid-IR weak"--they have strong near-IR emission and relatively weaker mid-IR emission, and most of the star formation is unobscured. We present new IR templates that capture how the mid-IR to near-IR emission changes with increasing infrared luminosity. We present two color-color diagrams to separate mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find that these color diagrams are most effective when used in conjunction with each other. We present the first number counts of 10$\mu$m sources and find that there are $\lesssim10$ IR AGN per MIRI pointing, possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to low metallicities or low dust content). We conclude that MIRI is most effective at observing moderate luminosity ($L_{\rm IR}=10^9-10^{10}L_\odot$) galaxies at $z=1-2$, and that photometry alone is not effective at identifying AGN within this faint population., Comment: 21 pages, 10 figures. Resubmitted to ApJS after revision

The Cosmic Evolution Early Release Science Survey (CEERS), targeting the Extended Groth Strip extragalactic field, is one of the JWST Director's Discretionary Early Release Science programs. To date, all observations have been executed and include NIRCam/MIRI imaging and NIRSpec/NIRCam spectroscopic exposures. Here, we discuss the MIRI imaging, which includes eight pointings, four of which provide deep imaging with the bluer bands (F560W, F770W) and four with contiguous wavelength coverage in F1000W, F1280W, F1500W, and F1800W, where two of these also include coverage in F770W and F2100W. We present a summary of the data, the data quality, and data reduction. The data reduction is based on the JWST Calibration Pipeline combined with custom modifications and additional steps designed to enhance the output quality, including improvements in astrometry and the removal of detector artifacts. We estimate the image depth of the reduced mosaics, and show that these generally agree with expectations from the Exposure Time Calculator. We compare the MIRI F560W and F770W flux densities for bright sources to measurements from Spitzer/IRAC Ch3 (5.8 $\mu$m) and Ch4 (8.0 $\mu$m), and we find that they agree with systematic differences of $<0.1$ mag. For the redder MIRI bands, we assess their quality by studying the spectral energy distributions (SEDs) of Galactic stars. The SEDs are consistent with the expected Rayleigh-Jeans law with a deviation $\sim 0.03$ mag, indicating that the MIRI colors are reliable. We also discuss all publicly released data products (images and source catalogs), which are available on the CEERS website (https://ceers.github.io/)., Comment: 17 pages, 11 figures, and 4 tables. ApJL in press

Current methods of identifying the ionizing source of nebular emission in galaxies are well defined for the era of single fiber spectroscopy, but still struggle to differentiate the complex and overlapping ionization sources in some galaxies. With the advent of integral field spectroscopy, the limits of these previous classification schemes are more apparent. We propose a new method for distinguishing the ionizing source in resolved galaxy spectra by use of a multi-dimensional diagnostic diagram that compares emission line ratios with velocity dispersion on a spaxel by spaxel basis within a galaxy. This new method is tested using the SAMI Galaxy Survey Data Release 3, which contains 3068 galaxies at z $<$ 0.12. Our results are released as ionization maps available alongside the SAMI DR3 public data. Our method accounts for a more diverse range of ionization sources than the standard suite of emission line diagnostics; we find 1433 galaxies with significant contribution from non-star-forming ionization using our improved method as compared to 316 galaxies identified using only emission line ratio diagnostics. Within these galaxies, we further identify 886 galaxies hosting unique signatures inconsistent with standard ionization by H2 regions, AGN, or shocks. These galaxies span a wide range of masses and morphological types and comprise a sizable portion of the galaxies used in our sample. With our revised method, we show that emission line diagnostics alone do not adequately differentiate the multiple ways to ionize gas within a galaxy.

We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey (CEERS) to measure rest-frame optical emission-line of 155 galaxies at z>2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the Ha, [OIII]/Hb, and [NeIII]/[OII] emission lines in terms of redshift by comparing to lower redshift SDSS and CLEAR samples. We find a significant (>3$\sigma$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$\sigma$) correlation with redshift. We compare [OIII]/Hb and [NeIII]/[OII] to stellar mass and Hb SFR. We find that both emission-line ratios have a correlation with Hb SFR and an anti-correlation with stellar mass across the redshifts 0

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

The photoionization model of narrow-line regions (NLRs) in active galactic nuclei (AGNs) has been investigated for decades. Many published models are restricted to simple linear scaling abundance relations, dust-free assumption, uniform AGN radiation field, and using one specific photoionization code, which restricts them from providing a satisfactory prediction on a broad range of AGN observations. Through a comprehensive investigation, here we present how the choice of abundance scaling relations, dust inclusion, AGN radiation fields, and different photoionization codes CLOUDY and MAPPINGS affect the predictions on the strength of strong UV, optical, and infrared emission lines. We find the dust-depleted radiation pressure-dominated AGN model built with the latest non-linear abundance sets and photoionization code MAPPINGS V are consistent with AGN observations across a broad range of wavelengths. We also assess new potential HII-AGN separation diagrams in the optical and UV wavelengths., Comment: 23 pages, 11 figures, Accepted for publication in ApJ

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.512, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products., Comment: 21 pages, 9 figures, submitted to AAS Journals

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

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4

Bars play an important role in mixing material in the inner regions of galaxies and stimulating radial migration. Previous observations have found evidence for the impact of a bar on metallicity gradients but the effect is still inconclusive. We use the TYPHOON/PrISM survey to investigate the metallicity gradients along and beyond the bar region across the entire star-forming disk of five nearby galaxies. Using emission line diagrams to identify star-forming spaxels, we recover the global metallicity gradients ranging from -0.0162 to -0.073 dex/kpc with evidence that the galactic bars act as an agent in affecting in-situ star formation as well as the motions of gas and stars. We observe cases with a `shallow-steep' metallicity radial profile, with evidence of the bar flattening the metallicity gradients inside the bar region (NGC~5068 and NGC~1566) and also note instances where the bar appears to drive a steeper metallicity gradient producing `steep-shallow' metallicity profiles (NGC~1365 and NGC~1744). For NGC~2835, a `steep-shallow' metallicity gradient break occurs at a distance $\sim$ 4 times the bar radius, which is more likely driven by gas accretion to the outskirt of the galaxy instead of the bar. The variation of metallicity gradients around the bar region traces the fluctuations of star formation rate surface density in NGC~1365, NGC~1566 and NGC~1744. A larger sample combined with hydrodynamical simulations is required to further explore the diversity and the relative importance of different ISM mixing mechanisms on the gas-phase metallicity gradients in local galaxies., Comment: 17 pages, 6 figures, 2 tables, published on MNRAS

Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M$\gtrsim 10 \, \mathrm{M}_\odot$, $v/v_{\rm{crit}} \gtrsim 0.4$). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances., Comment: 6 pages, 1 figure, published in Proceedings of the International Astronomical Union No. 366, 2022 titled "The Origin of Outflows in Evolved Stars"

We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them., Comment: 17 pages, 11 figures, submitted to MNRAS

This paper presents deep K-band spectroscopic observations of galaxies at z=3-4 with composite photometric rest-frame Hb+[OIII] equivalent widths EW_0>600A, comparable to the EW of galaxies observed during the epoch of reionisation (EoR, z>6). The typical spectroscopic [OIII] EW_0 and stellar mass of our targets is ~ 700A and log(M_star/M_sun)=8.98. By stacking the [OIII] emission profiles, we find evidence of a weak broad component with F_broad/F_narrow ~ 0.2 and velocity width sigma_{broad} ~ 170 km/s. The strength and velocity width of the broad component does not change significantly with stellar mass and [OIII] EW_0 of the stacked sample. Assuming similar broad component profiles for [OIII] and Halpha emission, we estimate a mass loading factor ~0.2, similar to low stellar mass galaxies at z>1 even if the star formation rates of our sample is 10 times higher. We hypothesize that either the multi-phase nature of supernovae driven outflows or the suppression of winds in the extreme star-forming regime is responsible for the weak signature of outflows in the EoR analogues., Comment: 14 pages, 9 figures, accepted for publication in MNRAS

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations (Yang et al. in prep). We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known pre-launch issues such as 1/f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with pre-launch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated dataset in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public datasets released from the Directors Discretionary Early Release Science (DD-ERS) program., Comment: 27 pages, 14 figures, submitted to ApJ. Accompanying CEERS public Data Release 0.5 available at ceers.github.io/releases.html

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3

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

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

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state., Comment: Submitted to ApJ; 12 Pages, 6 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

We create a photoionization model embedded in the turbulent ISM by using the state-of-the-art Messenger Monte-Carlo MAPPINGS~V code (M$^3$) in conjunction with the CMFGEN stellar atmosphere model. We show that the turbulent ISM causes the inhomogeneity of electron temperature and density within the nebula. The fluctuation in the turbulent ISM creates complex ionization structures seen in nearby nebulae. The inhomogeneous density distribution within the nebula creates a significant scatter on the spatially-resolved standard optical diagnostic diagrams, which cannot be represented by the spherical constant density photoionization model. We analyze the dependence of different optical emission lines on the complexity of nebular geometry, finding that the emission-lines residing on the nebular boundary are highly sensitive to the complexity of nebular geometry, while the emission-lines produced throughout the nebula are sensitive to the density distribution of the ISM within the nebula. Our fractal photoionization model demonstrates that a complex nebular geometry is required for accurate modeling of HII regions and emission-line galaxies, especially for the high-redshift galaxies, where the ISM is highly turbulent based on the increasing observational evidence., Comment: Accepted for Publication in ApJL

We analyse the spatial statistics of the 2D gas-phase oxygen abundance distributions in a sample of 219 local galaxies. We introduce a new adaptive binning technique to enhance the signal-to-noise ratio of weak lines, which we use to produce well-filled metallicity maps for these galaxies. We show that the two-point correlation functions computed from the metallicity distributions after removing radial gradients are in most cases well described by a simple injection-diffusion model. Fitting the data to this model yields the correlation length $l_{\rm corr}$, which describes the characteristic interstellar medium mixing length scale. We find typical correlation lengths $l_{\rm corr} \sim 1$ kpc, with a strong correlation between $l_{\rm corr}$ and stellar mass, star formation rate, and effective radius, a weak correlation with Hubble type, and significantly elevated values of $l_{\rm corr}$ in interacting or merging galaxies. We show that the trend with star formation rate can be reproduced by a simple transport+feedback model of interstellar medium turbulence at high star formation rate, and plausibly also at low star formation rate if dwarf galaxy winds have large mass-loading factors. We also report the first measurements of the injection width that describes the initial radii over which supernova remnants deposit metals. Inside this radius the metallicity correlation function is not purely the product of a competition between injection and diffusion. We show that this size scale is generally smaller than 60 pc., Comment: 18 pages, 18 figures, 1 table, submitted to MNRAS. Comments are welcome

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 present spectroscopic confirmation of candidate strong gravitational lenses using the Keck Observatory and Very Large Telescope as part of our ASTRO 3D Galaxy Evolution with Lenses (AGEL) survey. We confirm that 1) search methods using Convolutional Neural Networks (CNN) with visual inspection successfully identify strong gravitational lenses and 2) the lenses are at higher redshifts relative to existing surveys due to the combination of deeper and higher resolution imaging from DECam and spectroscopy spanning optical to near-infrared wavelengths. We measure 104 redshifts in 77 systems selected from a catalog in the DES and DECaLS imaging fields (r<22 mag). Combining our results with published redshifts, we present redshifts for 68 lenses and establish that CNN-based searches are highly effective for use in future imaging surveys with a success rate of 88% (defined as 68/77). We report 53 strong lenses with spectroscopic redshifts for both the deflector and source (z_src>z_defl), and 15 lenses with a spectroscopic redshift for either the deflector (z_defl>0.21) or source (z_src>1.34). For the 68 lenses, the deflectors and sources have average redshifts and standard deviations of 0.58+/-0.14 and 1.92+/-0.59 respectively, and corresponding redshift ranges of (0.210.5 that are ideal for follow-up studies to track how mass density profiles evolve with redshift. Our goal with AGEL is to spectroscopically confirm ~100 strong gravitational lenses that can be observed from both hemispheres throughout the year. The AGEL survey is a resource for refining automated all-sky searches and addressing a range of questions in astrophysics and cosmology., Comment: Updated final version of manuscript published by the Astronomical Journal

We use deep spectroscopy from the Hubble Space Telescope (HST) Wide-Field-Camera 3 (WFC3) IR grisms combined with broad-band photometry to study the stellar populations, gas ionization and chemical abundances in star-forming galaxies at $z\sim 1.1-2.3$. The data stem from the CANDELS Lyman-$\alpha$ Emission At Reionization (CLEAR) survey. At these redshifts the grism spectroscopy measure the [OII] 3727, 3729, [OIII] 4959, 5008, H-$\beta$ strong emission features, which constrain the ionization parameter and oxygen abundance of the nebular gas. We compare the line flux measurements to predictions from updated photoionization models (MAPPINGS (Kewley et al. 2019), which include an updated treatment of nebular gas pressure, log P/k = $n_e T_e$. Compared to low-redshift samples ($z\sim 0.2$) at fixed stellar mass, llog M / M$_\odot$ = 9.4-9.8, the CLEAR galaxies at z=1.35 (z=1.90) have lower gas-phase metallicity, $\Delta$(log Z) = 0.25 (0.35) dex, and higher ionization parameters, $\Delta$(log q) = 0.25 (0.35) dex, where U = q/c. We provide updated analytic calibrations between the [OIII], [OII], and H-$\beta$ emission line ratios, metallicity, and ionization parameter. The CLEAR galaxies show that at fixed stellar mass, the gas ionization parameter is correlated with the galaxy specific star-formation rates (sSFRs), where $\Delta$ log q = 0.4 $\Delta$(log sSFR), derived from changes in the strength of galaxy H-$\beta$ equivalent width. We interpret this as a consequence of higher gas densities, lower gas covering fractions, combined with higher escape fraction of H-ionizing photons. We discuss both tests to confirm these assertions and implications this has for future observations of galaxies at higher redshifts., Comment: 37 pages, plethora of figures, accepted for publication in ApJ. Updated with accepted version

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

The Messenger Interface Monte-Carlo Mappings V (M^3) is a photoionization code adopting the fully self-consistent Monte-Carlo radiative transfer technique, which presents a major advance over previous photoionization models with simple geometries. M^3 is designed for modeling nebulae in arbitrary three-dimensional geometries. In this paper, we describe the Monte-Carlo radiative transfer technique and the microphysics implemented in M^3, including the photoionization, collisional ionization, the free-free and free-bound recombination, and two-photon radiation. We put M^3 through the Lexington/Meudon benchmarks to test the reliability of the new code. We apply M^3 to three HII region models with fiducial geometries, demonstrating that M^3 is capable of dealing with nebulae with complex geometries. M^3 is a promising tool for understanding emission-line behavior in the era of SDSS-V/LVM and JWST, which will provide high-quality data of spatially-resolved nearby HII regions and highly turbulent local and high-redshift HII regions., Comment: Published on ApJ, 28 pages, 16 figures, 4 tables

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 investigate the mean locally-measured velocity dispersions of ionised gas ($\sigma_{\rm gas}$) and stars ($\sigma_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $\sigma_*$ tends to be larger than $\sigma_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference between $\sigma_{\rm gas}$ and $\sigma_*$ ($\Delta\sigma$) correlates with various galaxy properties. We establish that the strongest correlation of $\Delta\sigma$ is with beam smearing, which inflates $\sigma_{\rm gas}$ more than $\sigma_*$, introducing a dependence of $\Delta\sigma$ on both the effective radius relative to the point spread function and velocity gradients. The second-strongest correlation is with the contribution of active galactic nuclei (AGN) (or evolved stars) to the ionised gas emission, implying the gas velocity dispersion is strongly affected by the power source. In contrast, using the velocity dispersion measured from integrated spectra ($\sigma_{\rm aper}$) results in less correlation between the aperture-based $\Delta\sigma$ ($\Delta\sigma_{\rm aper}$) and the power source. This suggests that the AGN (or old stars) dynamically heat the gas without causing significant deviations from dynamical equilibrium. Although the variation of $\Delta\sigma_{\rm aper}$ is much smaller than that of $\Delta\sigma$, a correlation between $\Delta\sigma_{\rm aper}$ and gas velocity gradient is still detected, implying there is a small bias in dynamical masses derived from stellar and ionised gas velocity dispersions., Comment: 18 pages, 13 figures