Showing 1–50 of 195 results for author: Kewley, L J

Does the Fundamental Metallicity Relation Evolve with Redshift? II: The Evolution in Normalisation of the Mass-Metallicity Relation

Authors: Alex M. Garcia, Paul Torrey, Sara L. Ellison, Kathryn Grasha, Qian-Hui Chen, Z. S. Hemler, Dhruv T. Zimmerman, Ruby J. Wright, Henry R. M. Zovaro, Erica J. Nelson, Ryan L. Sanders, Lisa J. Kewley, Lars Hernquist

Abstract: 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 JWST data indicate deviations from this model. In this… ▽ More 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 JWST data indicate deviations from this model. 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 Illustris and SIMBA and dynamic FMRs in 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 be dynamic as well. Overall, our findings imply that the current FMR framework neglects important variations in the baryon cycle through cosmic time. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856). Comments welcome

The significant contribution of supersoft X-ray Sources to the nebular HeII line emission

Authors: Dian P. Triani, Rosanne Di Stefano, Lisa J. Kewley

Abstract: 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… ▽ More 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. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 17 pages, 7 figures, comments are welcomed

The mass-metallicity and fundamental metallicity relations in non-AGN and AGN-host galaxies

Authors: Song-Lin Li, Kathryn Grasha, Mark R. Krumholz, Emily Wisnioski, Ralph S. Sutherland, Lisa J. Kewley, Yan-Mei Chen, Zefeng Li

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

Does the Fundamental Metallicity Relation Evolve with Redshift? I: The Correlation Between Offsets from the Mass-Metallicity Relation and Star Formation Rate

Authors: Alex M. Garcia, Paul Torrey, Sara Ellison, Kathryn Grasha, Lars Hernquist, Henry R. M. Zovaro, Qian-Hui Chen, Z. S. Hemler, Lisa J. Kewley, Erica J. Nelson, Ruby J. Wright

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures (Figure 5 is summary). Accepted to MNRAS

Interplay of Stellar and Gas-Phase Metallicities: Unveiling Insights for Stellar Feedback Modeling with Illustris, IllustrisTNG, and EAGLE

Authors: Alex M. Garcia, Paul Torrey, Kathryn Grasha, Lars Hernquist, Sara Ellison, Henry R. M. Zovaro, Z. S. Hemler, Erica J. Nelson, Lisa J. Kewley

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 11 March, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures, Accepted to MNRAS

The Next Generation Deep Extragalactic Exploratory Public Near-Infrared Slitless Survey Epoch 1 (NGDEEP-NISS1): Extra-Galactic Star-formation and Active Galactic Nuclei at 0.5 < z < 3.6

Authors: Nor Pirzkal, Barry Rothberg, Casey Papovich, Lu Shen, Gene C. K. Leung, Micaela B. Bagley, Steven L. Finkelstein, Brittany N. Vanderhoof, Jennifer M. Lotz, Anton M. Koekemoer, Nimish P. Hathi, Yingjie Cheng, Nikko J. Cleri, Norman A. Grogin, L. Y. Aaron Yung, Mark Dickinson, Henry C. Ferguson, Jonathan P. Gardner, Intae Jung, Jeyhan S. Kartaltepe, Russell Ryan, Raymond C. Simons, Swara Ravindranath, Danielle A. Berg, Bren E. Backhaus , et al. (26 additional authors not shown)

Abstract: 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 firs… ▽ More 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$σ$ 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$β$), and dust content, of galaxies at $1<z<3.5$. The latter showing intriguing initial results on the applicability and assumptions made regarding the use of Case B recombination. Further, we identify the presence of active galactic nuclei (AGN) and infer the mass of their supermassive black holes (SMBHs) using broadened restframe MgII and H$β$ emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution (SED) models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts. △ Less

Submitted 20 April, 2024; v1 submitted 15 December, 2023; originally announced December 2023.

Comments: 33 pages, 28 Figures, Accepted (ApJ)

The SAMI Galaxy Survey: $Σ_{\rm SFR}$ drives the presence of complex emission line profiles in star-forming galaxies

Authors: Henry R. M. Zovaro, J. Trevor Mendel, Brent Groves, Lisa J. Kewley, Matthew Colless, Andrei Ristea, Luca Cortese, Sree Oh, Francesco D'Eugenio, Scott M. Croom, Ángel R. López-Sánchez, Jesse van de Sande, Sarah Brough, Anne M. Medling, Joss Bland-Hawthorn, Julia J. Bryant

Abstract: 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$α$ components, and distinguis… ▽ More 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$α$ 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}$ ($Σ_{{\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 $Σ_{\rm SFR}$, but by the $Σ_{{\rm SFR},R_{\rm e}}$ of the host galaxy. We also parametrise the clumpiness of the SFR within individual galaxies, and find that $Σ_{{\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. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 21 pages, 17 figures, accepted for publication in MNRAS

Galaxies Going Bananas: Inferring the 3D Geometry of High-Redshift Galaxies with JWST-CEERS

Authors: Viraj Pandya, Haowen Zhang, Marc Huertas-Company, Kartheik G. Iyer, Elizabeth McGrath, Guillermo Barro, Steven L. Finkelstein, Martin Kuemmel, William G. Hartley, Henry C. Ferguson, Jeyhan S. Kartaltepe, Joel Primack, Avishai Dekel, Sandra M. Faber, David C. Koo, Greg L. Bryan, Rachel S. Somerville, Ricardo O. Amorin, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Emmanuel Bertin, Luca Costantin, Romeel Dave, Mark Dickinson , et al. (31 additional authors not shown)

Abstract: 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 e… ▽ More 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érsic 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. △ Less

Submitted 15 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

CEERS Key Paper VII: JWST/MIRI Reveals a Faint Population of Galaxies at Cosmic Noon Unseen by Spitzer

Authors: Allison Kirkpatrick, Guang Yang, Aurelien Le Bail, Greg Troiani, Eric F. Bell, Nikko J. Cleri, David Elbaz, Steven L. Finkelstein, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Dale D. Kocevski, Ray A. Lucas, Jed McKinney, Casey Papovich, Pablo G. Perez-Gonzalez, Alexander de la Vega, Micaela B. Bagley, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Adriano Fontana, Andrea Grazian, Norman A. Grogin, Pablo Arrabal Haro , et al. (11 additional authors not shown)

Abstract: 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$μ$m po… ▽ More 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$μ$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$μ$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. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 21 pages, 10 figures. Resubmitted to ApJS after revision

CEERS Key Paper VIII: Emission Line Ratios from NIRSpec and NIRCam Wide-Field Slitless Spectroscopy at z>2

Authors: Bren E. Backhaus, Jonathan R. Trump, Nor Pirzkal, Guillermo Barro, Steven L. Finkelstein, Pablo Arrabal Haro, Raymond C. Simons, Jessica Wessner, Nikko J. Cleri, Michaela Hirschmann, Micaela B. Bagley, David C. Nicholls, Mark Dickinson, Jeyhan S. Kartaltepe, Casey Papovich, Dale D. Kocevski, Anton M. Koekemoer, Laura Bisigello, Anne E. Jaskot, Ray A. Lucas, Intae Jung, Stephen M. Wilkins, L. Y. Aaron Yung, Henry C. Ferguson, Adriano Fontana , et al. (15 additional authors not shown)

Abstract: 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.… ▽ More 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$σ$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$σ$) 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<z<9. Comparison with MAPPINGS~V models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We additionally compare to IllustriousTNG predictions and find that they effectively describe the highest [OIII]/Hb ratios observed in our sample, without the need to invoke MAPPINGS models with significant shock ionizionation components. △ Less

Submitted 7 September, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 16 pages, 11 figures

CLASSY VIII: Exploring the Source of Ionization with UV ISM diagnostics in local High-$z$ Analogs

Authors: Matilde Mingozzi, Bethan L. James, Danielle A. Berg, Karla Z. Arellano-Córdova, Adele Plat, Claudia Scarlata, Alessandra Aloisi, Ricardo O. Amorín, Jarle Brinchmann, Stéphane Charlot, John Chisholm, Anna Feltre, Simon Gazagnes, Matthew Hayes, Timothy Heckman, Svea Hernandez, Lisa J. Kewley, Nimisha Kumari, Claus Leitherer, Crystal L. Martin, Michael Maseda, Themiya Nanayakkara, Swara Ravindranath, Jane R. Rigby, Peter Senchyna , et al. (5 additional authors not shown)

Abstract: 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… ▽ More 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] $λλ$1907,9 He II $\lambda1640$ and O III] $λ$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] $λ$1666 still allow us to define a SF-locus with some caveats. Diagrams including C IV $λλ$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. △ Less

Submitted 3 December, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 36 pages, 15 figures, accepted for publication in ApJ

A New Photoionization Model of the Narrow Line Region in Active Galactic Nuclei

Authors: Peixin Zhu, Lisa J. Kewley, Ralph S. Sutherland

Abstract: 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 observation… ▽ More 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. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Comments: 23 pages, 11 figures, Accepted for publication in ApJ

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey

Authors: Micaela B. Bagley, Nor Pirzkal, Steven L. Finkelstein, Casey Papovich, Danielle A. Berg, Jennifer M. Lotz, Gene C. K. Leung, Henry C. Ferguson, Anton M. Koekemoer, Mark Dickinson, Jeyhan S. Kartaltepe, Dale D. Kocevski, Rachel S. Somerville, L. Y. Aaron Yung, Bren E. Backhaus, Caitlin M. Casey, Marco Castellano, Óscar A. Chávez Ortiz, Katherine Chworowsky, Isabella G. Cox, Romeel Davé, Kelcey Davis, Vicente Estrada-Carpenter, Adriano Fontana, Seiji Fujimoto , et al. (23 additional authors not shown)

Abstract: 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 r… ▽ More 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.5<z<4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam (m=30.6-30.9, 5sigma) to discover galaxies to z>12, 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. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 21 pages, 9 figures, submitted to AAS Journals

CEERS Spectroscopic Confirmation of NIRCam-Selected z > 8 Galaxy Candidates with JWST/NIRSpec: Initial Characterization of their Properties

Authors: Seiji Fujimoto, Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Rebecca L. Larson, Denis Burgarella, Micaela B. Bagley, Peter Behroozi, Katherine Chworowsky, Michaela Hirschmann, Jonathan R. Trump, Stephen M. Wilkins, L. Y. Aaron Yung, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Henry C. Ferguson, Adriano Fontana, Norman A. Grogin, Andrea Grazian, Lisa J. Kewley, Dale D. Kocevski, Jennifer M. Lotz, Laura Pentericci , et al. (16 additional authors not shown)

Abstract: 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μ$m either with the NIRSpec prism or its three medium resolution gratings… ▽ More 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μ$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$β$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$α$ 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$β$ equivalent widths ($\simeq$1100$Å$), and elevated production efficiency of ionizing photons ($\log(ξ_{\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$β$ 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. △ Less

Submitted 26 May, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

CEERS Key Paper IV: Galaxies at $4 < z < 9$ are Bluer than They Appear -- Characterizing Galaxy Stellar Populations from Rest-Frame $\sim 1$ micron Imaging

Authors: Casey Papovich, Justin Cole, Guang Yang, Steven L. Finkelstein, Guillermo Barro, Véronique Buat, Denis Burgarella, Pablo G. Pérez-González, Paola Santini, Lise-Marie Seillé, Lu Shen, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, Caitlin M. Casey, Marco Castellano, Katherine Chworowsky, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Mark Dickinson, Henry C. Ferguson, Adriano Fontana , et al. (24 additional authors not shown)

Abstract: We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame… ▽ More We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame spectral-energy distribution (SED) to nearly 1 micron for galaxies in this redshift range. By modeling the galaxies' SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å) $-$ $I$-band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower (stellar)-mass-to-light ratios. The MIRI data reduce the stellar masses by $\langle Δ\log M_\ast\rangle=0.25$ dex at $4<z<6$ (a factor of 1.8) and 0.37 dex at $6<z<9$ (a factor of 2.3). This also reduces the star-formation rates (SFRs) by $\langle Δ\log\mathrm{SFR} \rangle=0.14$ dex at $4<z<6$ and 0.27 dex at $6<z<9$. The MIRI data also improve constraints on the allowable stellar mass formed in early star-formation. We model this using a star-formation history that includes both a "burst' at $z_f=100$ and a slowly varying ("delayed-$τ$") model. The MIRI data reduce the allowable stellar mass by 0.6 dex at $4<z< 6$ and by $\approx$1 dex at $6<z<9$. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early universe ($z\approx9$). Therefore, observations of rest-frame $\gtrsim$1 $μ$m are paramount for constraining the stellar-mass build-up in galaxies at very high-redshifts. △ Less

Submitted 25 March, 2023; v1 submitted 30 December, 2022; originally announced January 2023.

Comments: Updated with accepted ApJ version. Part of the CEERS Focus Issue. 27 pages, many figures (4 Figure Sets, available upon reasonable request)

Metallicity Gradient of Barred Galaxies with TYPHOON

Authors: Qian-Hui Chen, Kathryn Grasha, Andrew J. Battisti, Lisa J. Kewley, Barry F. Madore, Mark Seibert, Jeff A. Rich, Rachael L. Beaton

Abstract: 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 g… ▽ More 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. △ Less

Submitted 24 January, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 17 pages, 6 figures, 2 tables, published on MNRAS

Can pre-supernova winds from massive stars enrich the interstellar medium with nitrogen at high redshift?

Authors: Arpita Roy, Mark R. Krumholz, Michael A. Dopita, Ralph S. Sutherland, Lisa J. Kewley, Alexander Heger

Abstract: 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… ▽ More 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. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 6 pages, 1 figure, published in Proceedings of the International Astronomical Union No. 366, 2022 titled "The Origin of Outflows in Evolved Stars"

Journal ref: Proceedings of the International Astronomical Union, 16(S366), 33-38, 2022

Gas-phase metallicity break radii of star-forming galaxies in IllustrisTNG

Authors: Alex M. Garcia, Paul Torrey, Z. S. Hemler, Lars Hernquist, Lisa J. Kewley, Erica J. Nelson, Kathryn Grasha, Henry R. M. Zovaro, Qian-Hui Chen

Abstract: 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 radi… ▽ More 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. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 17 pages, 11 figures, submitted to MNRAS

MOSEL Survey: Extremely weak outflows in EoR analogues at z=3-4

Authors: Anshu Gupta, Kim-Vy Tran, Trevor Mendel, Anishya Harshan, Ben Forrest, R. L. Davies, Emily Wisnioski, Themiya Nanayakkara, Glenn G. Kacprzak, Lisa J. Kewley

Abstract: 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 fin… ▽ More 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. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 14 pages, 9 figures, accepted for publication in MNRAS

CEERS Key Paper I: An Early Look into the First 500 Myr of Galaxy Formation with JWST

Authors: Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Stephen M. Wilkins, Jeyhan S. Kartaltepe, Casey Papovich, L. Y. Aaron Yung, Pablo Arrabal Haro, Peter Behroozi, Mark Dickinson, Dale D. Kocevski, Anton M. Koekemoer, Rebecca L. Larson, Aurelien Le Bail, Alexa M. Morales, Pablo G. Perez-Gonzalez, Denis Burgarella, Romeel Dave, Michaela Hirschmann, Rachel S. Somerville, Stijn Wuyts, Volker Bromm, Caitlin M. Casey, Adriano Fontana, Seiji Fujimoto , et al. (42 additional authors not shown)

Abstract: 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 det… ▽ More 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. △ Less

Submitted 4 November, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: Replaced with published version

CEERS Epoch 1 NIRCam Imaging: Reduction Methods and Simulations Enabling Early JWST Science Results

Authors: Micaela B. Bagley, Steven L. Finkelstein, Anton M. Koekemoer, Henry C. Ferguson, Pablo Arrabal Haro, Mark Dickinson, Jeyhan S. Kartaltepe, Casey Papovich, Pablo G. Pérez-González, Nor Pirzkal, Rachel S. Somerville, Christopher N. A. Willmer, Guang Yang, L. Y. Aaron Yung, Adriano Fontana, Andrea Grazian, Norman A. Grogin, Michaela Hirschmann, Lisa J. Kewley, Allison Kirkpatrick, Dale D. Kocevski, Jennifer M. Lotz, Aubrey Medrano, Alexa M. Morales, Laura Pentericci , et al. (12 additional authors not shown)

Abstract: 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). W… ▽ More 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. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: 27 pages, 14 figures, submitted to ApJ. Accompanying CEERS public Data Release 0.5 available at ceers.github.io/releases.html

CEERS Key Paper V: A triality on the nature of HST-dark galaxies

Authors: Pablo G. Pérez-González, Guillermo Barro, Marianna Annunziatella, Luca Costantin, Ángela García-Argumánez, Elizabeth J. McGrath, Rosa M. Mérida, Jorge A. Zavala, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Peter Behroozi, Eric F. Bell, Laura Bisigello, Véronique Buat, Antonello Calabrò, Caitlin M. Casey, Nikko J. Cleri, Rosemary T. Coogan, M. C. Cooper, Asantha R. Cooray, Avishai Dekel, Mark Dickinson, David Elbaz, Henry C. Ferguson , et al. (30 additional authors not shown)

Abstract: 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 photometr… ▽ More 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 2<z<6 with masses 9<log M/M_sun<11 and a variety of specific SFRs (<1 to >100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3<z<5, masses log M/M_sun~10 and post-starburst stellar mass-weighted ages (0.5-1 Gyr); and (3) 11% are strong young starbursts with indications of high-EW emission lines (typically, [OIII]+Hbeta) at 6<z<7 and log M/M_sun~9.5. The sample is dominated by disk-like galaxies with a remarkable compactness for XELG-z6 (effective radii smaller than 0.4 kpc). Large attenuations in SFGs, 2<A(V)<5 mag, are found within 1.5 times the effective radius, approximately 2 kpc, while QGs present A(V)~0.2 mag. Our SED-fitting technique reproduces the expected dust emission luminosities of IR-bright and sub-millimeter galaxies. This study implies high levels of star formation activity between z~20 and z~10, where virtually 100% of our galaxies had already formed 10^8 M_sun of their stellar content, 60% of them had assembled 10^9 M_sun, and 10% up to 10^10 M_sun (in situ or ex situ). (abridged) △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Published in CEERS ApJL Focus Issue, ApJL 946, L16

Journal ref: ApJL 946, L16 https://iopscience.iop.org/article/10.3847/2041-8213/acb3a5

CEERS Key Paper III: The Diversity of Galaxy Structure and Morphology at z=3-9 with JWST

Authors: Jeyhan S. Kartaltepe, Caitlin Rose, Brittany N. Vanderhoof, Elizabeth J. McGrath, Luca Costantin, Isabella G. Cox, L. Y. Aaron Yung, Dale D. Kocevski, Stijn Wuyts, Henry C. Ferguson Brett H. Andrews, Micaela B. Bagley, Steven L. Finkelstein, Ricardo O. Amorin, Pablo Arrabal Haro, Bren E. Backhaus, Peter Behroozi, Laura Bisigello, Antonello Calabro, Caitlin M. Casey, Rosemary T. Coogan, Darren Croton, Alexander de la Vega, Mark Dickinson, M. C. Cooper, Adriano Fontana , et al. (36 additional authors not shown)

Abstract: 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, wit… ▽ More 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érsic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher Sérsic 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. △ Less

Submitted 13 January, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: Accepted for publication to ApJL, 24 pages, 14 figures

CEERS Key Paper III: The Resolved Host Properties of AGN at 3 < z < 5 with JWST

Authors: Dale D. Kocevski, Guillermo Barro, Elizabeth J. McGrath, Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Shardha Jogee, Guang Yang, Mark Dickinson, Nimish P. Hathi, Bren E. Backhaus, Eric F. Bell, Laura Bisigello, Véronique Buat, Denis Burgarella, Caitlin M. Casey, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Darren Croton, Emanuele Daddi, Adriano Fontana, Seiji Fujimoto, Jonathan P. Gardner, Eric Gawiser , et al. (34 additional authors not shown)

Abstract: 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 redshif… ▽ More 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. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ; 12 Pages, 6 Figures

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: 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 als… ▽ More 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. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: 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 f… ▽ More 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. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 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

The Physical Conditions of Emission-Line Galaxies at Cosmic Dawn from JWST/NIRSpec Spectroscopy in the SMACS 0723 Early Release Observations

Authors: Jonathan R. Trump, Pablo Arrabal Haro, Raymond C. Simons, Bren E. Backhaus, Ricardo O. Amorín, Mark Dickinson, Vital Fernández, Casey Papovich, David C. Nicholls, Lisa J. Kewley, Samantha W. Brunker, John J. Salzer, Stephen M. Wilkins, Omar Almaini, Micaela B. Bagley, Danielle A. Berg, Rachana Bhatawdekar, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Laure Ciesla, Nikko J. Cleri, Justin W. Cole , et al. (39 additional authors not shown)

Abstract: 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 wav… ▽ More 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]$λ$5008/H$β$ ratios, similar [OIII]$λ$4364/H$γ$ ratios, and higher ($\sim$0.5 dex) [NeIII]$λ$3870/[OII]$λ$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$λ$3870/[OII]$λ$3728, [OIII]$λ$4364/H$γ$, and [OIII]$λ$5008/H$β$ 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]$λ$4364/H$γ$ and [OIII]$λ$(4960+5008)/H$β$ 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. △ Less

Submitted 19 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

Theoretically Modelling Photoionized Regions with Fractal Geometry in Three Dimensions

Authors: Yifei Jin, Lisa J. Kewley, Ralph S. Sutherland

Abstract: 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 ne… ▽ More 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. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: Accepted for Publication in ApJL

Spatial metallicity distribution statistics at $\lesssim 100$ pc scales in the AMUSING++ nearby galaxy sample

Authors: Zefeng Li, Emily Wisnioski, J. Trevor Mendel, Mark R. Krumholz, Lisa J. Kewley, Carlos López-Cobá, Sebastian F. Sánchez, Joseph P. Anderson, Lluís Galbany

Abstract: 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 r… ▽ More 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. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 18 pages, 18 figures, 1 table, submitted to MNRAS. Comments are welcome

CLASSY II: A technical Overview of the COS Legacy Archive Spectroscopic SurveY

Authors: Bethan L. James, Danielle A. Berg, Teagan King, David J. Sahnow, Matilde Mingozzi, John Chisholm, Timothy Heckman, Crystal L. Martin, Dan P. Stark, The Classy Team, :, Alessandra Aloisi, Ricardo O. Amorín, Karla Z. Arellano-Córdova, Matthew Bayliss, Rongmon Bordoloi, Jarle Brinchmann, Stéphane Charlot, Jacopo Chevallard, Ilyse Clark, Dawn K. Erb, Anna Feltre, Matthew Hayes, Alaina Henry, Svea Hernandez , et al. (23 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: 31 pages, 17 figures, accepted for publication in ApJS

The AGEL Survey: Spectroscopic Confirmation of Strong Gravitational Lenses in the DES and DECaLS Fields Selected Using Convolutional Neural Networks

Authors: Kim-Vy H. Tran, Anishya Harshan, Karl Glazebrook, G. C. Keerthi Vasan, Tucker Jones, Colin Jacobs, Glenn G. Kacprzak, Tania M. Barone, Thomas E. Collett, Anshu Gupta, Astrid Henderson, Lisa J. Kewley, Sebastian Lopez, Themiya Nanayakkara, Ryan L. Sanders, Sarah M. Sweet

Abstract: 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 t… ▽ More 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.21<z_defl<0.89) and (0.88<z_src<3.55). The AGEL systems include 41 deflectors at zdefl>0.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. △ Less

Submitted 26 September, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Updated final version of manuscript published by the Astronomical Journal

CLASSY III: The Properties of Starburst-Driven Warm Ionized Outflows

Authors: Xinfeng Xu, Timothy Heckman, Alaina Henry, Danielle A. Berg, John Chisholm, Bethan L. James, Crystal L. Martin, Daniel P. Stark, Alessandra Aloisi, Ricardo O. Amorín, Karla Z. Arellano-Córdova, Rongmon Bordoloi, Stéphane Charlot, Zuyi Chen, Matthew Hayes, Matilde Mingozzi, Yuma Sugahara, Lisa J. Kewley, Masami Ouchi, Claudia Scarlata, Charles C. Steidel

Abstract: 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 en… ▽ More 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. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 34 pages, 16 figures, 6 tables, submitted to ApJ

The COS Legacy Archive Spectroscopy SurveY (CLASSY) Treasury Atlas

Authors: Danielle A. Berg, Bethan L. James, Teagan King, Meaghan Mcdonald, Zuyi Chen, John Chisholm, Timothy Heckman, Crystal L. Martin, Dan P. Stark, The Classy Team, :, Alessandra Aloisi, Ricardo O. AmorÍn, Karla Z. Arellano-CÓrdova, Matthew Bayliss, Rongmon Bordoloi, Jarle Brinchmann, StÉphane Charlot, Jacopo Chevallard, Ilyse Clark, Dawn K. Erb, Anna Feltre, Matthew Hayes, Alaina Henry, Svea Hernandez , et al. (24 additional authors not shown)

Abstract: 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 comprehensi… ▽ More 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. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Metallicity, ionization parameter, and pressure variations of HII regions in the TYPHOON spiral galaxies

Authors: K. Grasha, Q. H. Chen, A. J. Battisti, A. Acharyya, S. Ridolfo, E. Poehler, S. Mably, A. A. Verma, K. L. Hayward, A. Kharbanda, H. Poetrodjojo, M. Seibert, J. A. Rich, B. F. Madore, L. J. Kewley

Abstract: We present a spatially-resolved HII region study of the gas-phase metallicity, ionization parameter, and ISM pressure maps of 6 local star-forming and face-on spiral galaxies from the TYPHOON program. Self-consistent metallicity, ionization parameter, and pressure maps are calculated simultaneously through an iterative process to provide useful measures of the local chemical abundance and its rela… ▽ More We present a spatially-resolved HII region study of the gas-phase metallicity, ionization parameter, and ISM pressure maps of 6 local star-forming and face-on spiral galaxies from the TYPHOON program. Self-consistent metallicity, ionization parameter, and pressure maps are calculated simultaneously through an iterative process to provide useful measures of the local chemical abundance and its relation to localized ISM properties. We constrain the presence of azimuthal variations in metallicity by measuring the residual metallicity offset $Δ$(O/H) after subtracting the linear fits to the radial metallicity profiles. We however find weak evidence of azimuthal variations in most of the galaxies, with small (mean 0.03 dex) scatter. The galaxies instead reveal that HII regions with enhanced and reduced abundances are found distributed throughout the disk. While the spiral pattern plays a role in organizing the ISM, it alone does not establish the relatively uniform azimuthal variations we observe. Differences in the metal abundances are more likely driven by the strong correlations with the local physical conditions. We find a strong and positive correlation between the ionization parameter and the local abundances as measured by the relative metallicity offset $Δ$(O/H), indicating a tight relationship between local physical conditions and their localized enrichment of the ISM. Additionally, we demonstrate the impact of unresolved observations on the measured ISM properties by rebinning the datacubes to simulate low-resolution (1 kpc) observations, typical of large IFU surveys. We find that the ionization parameter and ISM pressure diagnostics are impacted by the loss of resolution such that their measured values are larger relative to the measured values on sub-HII region scales. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 31 pages, 14 figures. Accepted for publication in ApJ

The SAMI Galaxy Survey: the difference between ionised gas and stellar velocity dispersions

Authors: Sree Oh, Matthew Colless, Francesco D'Eugenio, Scott M. Croom, Luca Cortese, Brent Groves, Lisa J. Kewley, Jesse van de Sande, Henry Zovaro, Mathew R. Varidel, Stefania Barsanti, Sarah Brough, Julia J. Bryant, Sarah Casura, Jon S. Lawrence, Nuria P. F. Lorente, Anne M. Medling, Matt S. Owers, Sukyoung K. Yi

Abstract: We investigate the mean locally-measured velocity dispersions of ionised gas ($σ_{\rm gas}$) and stars ($σ_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $σ_*$ tends to be larger than $σ_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference betwe… ▽ More We investigate the mean locally-measured velocity dispersions of ionised gas ($σ_{\rm gas}$) and stars ($σ_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $σ_*$ tends to be larger than $σ_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference between $σ_{\rm gas}$ and $σ_*$ ($Δσ$) correlates with various galaxy properties. We establish that the strongest correlation of $Δσ$ is with beam smearing, which inflates $σ_{\rm gas}$ more than $σ_*$, introducing a dependence of $Δσ$ 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 ($σ_{\rm aper}$) results in less correlation between the aperture-based $Δσ$ ($Δσ_{\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 $Δσ_{\rm aper}$ is much smaller than that of $Δσ$, a correlation between $Δσ_{\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. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 18 pages, 13 figures

The Black Hole-Galaxy Connection: Interplay between Feedback, Obscuration, and Host Galaxy Substructure

Authors: Stéphanie Juneau, Andy D. Goulding, Julie Banfield, Stefano Bianchi, Pierre-Alain Duc, I-Ting Ho, Michael A. Dopita, Julia Scharwächter, Franz E. Bauer, Brent Groves, David M. Alexander, Rebecca L. Davies, David Elbaz, Emily Freeland, Elise Hampton, Lisa J. Kewley, Robert Nikutta, Prajval Shastri, Xinwen Shu, Frédéric P. A. Vogt, Tao Wang, O. Ivy Wong, Jong-Hak Woo

Abstract: There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple… ▽ More There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN. △ Less

Submitted 10 January, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 31 pages, 16 figures (including two animations), accepted for publication in the Astrophysical Journal. For a preview of the animations, see https://noirlab.edu/science/videos/rotating-gas-disk-of-the-ngc-7582-galaxy for Figure 15, and see https://noirlab.edu/science/videos/black-hole-wind-from-the-ngc-7582-galaxy for Figure 16

Tracing the Ionization Structure of the Shocked Filaments of NGC 6240

Authors: Anne M. Medling, Lisa J. Kewley, Daniela Calzetti, George C. Privon, Kirsten Larson, Jeffrey A. Rich, Lee Armus, Mark G. Allen, Geoffrey V. Bicknell, Tanio Díaz-Santos, Timothy M. Heckman, Claus Leitherer, Claire E. Max, David S. N. Rupke, Ezequiel Treister, Hugo Messias, Alexander Y. Wagner

Abstract: We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and fi… ▽ More We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and filaments seen in shock tracers from warm molecular gas (H$_2$ 2.12$μ$m) to optical ionized gas ([O III], [N II], [S II], [O I]) and hot plasma (Fe XXV). In the most distinct bubble, we see a clear shock front traced by high [O III]/H$β$ and [O III]/[O I]. Cool molecular gas (CO(2-1)) is only present near the base of the bubble, towards the nuclei launching the outflow. We interpret the lack of molecular gas outside the bubble to mean that the shock front is not responsible for dissociating molecular gas, and conclude that the molecular clouds are partly shielded and either entrained briefly in the outflow, or left undisturbed while the hot wind flows around them. Elsewhere in the galaxy, shock-excited H$_2$ extends at least $\sim$4 kpc from the nuclei, tracing molecular gas even warmer than that between the nuclei, where the two galaxies' interstellar media are colliding. A ridgeline of high [O III]/H$β$ emission along the eastern arm aligns with the south nucleus' stellar disk minor axis; optical integral field spectroscopy from WiFeS suggests this highly ionized gas is centered at systemic velocity and likely photoionized by direct line-of-sight to the south AGN. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 27 pages, 18 figures; accepted for publication in ApJ

Revisiting the Giant Radio Galaxy ESO 422-G028: Part I. Discovery of a neutral inflow and recent star formation in a restarted giant

Authors: Henry R. M. Zovaro, Chris J. Riseley, Philip Taylor, Nicole P. H. Nesvadba, Tim J. Galvin, Umang Malik, Lisa J. Kewley

Abstract: Giant radio galaxies provide important clues into the life cycles and triggering mechanisms of radio jets. With large-scale jets spanning 1.8 Mpc, ESO 422-G028 ($z = 0.038$) is a giant radio galaxy that also exhibits signs of restarted jet activity in the form of pc-scale jets. We present a study of the spatially-resolved stellar and gas properties of ESO 422-G028 using optical integral field spec… ▽ More Giant radio galaxies provide important clues into the life cycles and triggering mechanisms of radio jets. With large-scale jets spanning 1.8 Mpc, ESO 422-G028 ($z = 0.038$) is a giant radio galaxy that also exhibits signs of restarted jet activity in the form of pc-scale jets. We present a study of the spatially-resolved stellar and gas properties of ESO 422-G028 using optical integral field spectroscopy from the WiFeS spectrograph. In addition to the majority $\sim 13\,\rm Gyr$ old stellar population, ESO 422-G028 exhibits a much younger ($\lesssim 10\,\rm Myr$ old) component with an estimated mass of $ 10^{7.6}\,\rm M_\odot$ which is predominantly located in the North-West region of the galaxy. Unusually, the ionised gas kinematics reveal two distinct disks traced by narrow ($σ_{\rm Hα} < 100 \,\rm km\,s^{-1}$) and broad ($σ_{\rm Hα} > 150 \,\rm km\,s^{-1}$) H$α$ emission respectively. Both ionised gas disks are misaligned with the axis of stellar rotation, suggesting an external origin. This is consistent with the prominent interstellar Na D absorption, which traces a $1 - 3 \,\rm M_\odot \, yr^{-1}$ inflow of neutral gas from the North. We posit that an inflow of gas - either from an accretion event or a gas-rich merger - has triggered both the starburst and the restarted jet activity, and that ESO 422-G028 is potentially on the brink of an epoch of powerful AGN activity. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 21 pages, 19 figures; submitted to the Monthly Notices of the Royal Astronomical Society

Detection of metallicity correlations in 100 nearby galaxies

Authors: Zefeng Li, Mark R. Krumholz, Emily Wisnioski, J. Trevor Mendel, Lisa J. Kewley, Sebastian F. Sánchez, Lluís Galbany

Abstract: In this paper we analyse the statistics of the 2D gas-phase oxygen abundance distributions of 100 nearby galaxies drawn from the Calar Alto Legacy Integral Field spectroscopy Area survey. After removing the large-scale radial metallicity gradient, we compute the two-point correlation functions of the resulting metallicity fluctuation maps. We detect correlations in the majority of our targets, whi… ▽ More In this paper we analyse the statistics of the 2D gas-phase oxygen abundance distributions of 100 nearby galaxies drawn from the Calar Alto Legacy Integral Field spectroscopy Area survey. After removing the large-scale radial metallicity gradient, we compute the two-point correlation functions of the resulting metallicity fluctuation maps. We detect correlations in the majority of our targets, which we show are significantly in excess of what is expected due to beam-smearing, and are robust against the choice of metallicity diagnostic. We show that the correlation functions are generally well-fit by the predictions of a simple model for stochastic metal injection coupled with diffusion, and from the model we show that, after accounting for the effects of both beam smearing and noise, the galaxies in our sample have characteristic correlation lengths of $\sim1$ kpc. Correlation lengths increase with both stellar mass and star formation rate, but show no significant variation with Hubble type, barredness, or merging state. We also find no evidence for a theoretically-predicted relationship between gas velocity dispersion and correlation length, though this may be due to the small dynamic range in gas velocity dispersion across our sample. Our results suggest that measurements of 2D metallicity correlation functions can be a powerful tool for studying galaxy evolution. △ Less

Submitted 30 April, 2021; originally announced April 2021.

Comments: 17 pages, 16 figures, 1 table, accepted for publication in MNRAS

Stromlo Stellar Tracks: non-solar scaled abundances for massive stars

Authors: K. Grasha, A. Roy, R. S. Sutherland, L. J. Kewley

Abstract: We present the Stromlo Stellar Tracks, a set of stellar evolutionary tracks, computed by modifying the Modules for Experiments in Stellar Astrophysics (MESA) 1D stellar evolution package, to fit the Galactic Concordance abundances for hot ($T > 8000$ K) massive ($\geq 10M_\odot$) Main-Sequence (MS) stars. Until now, all stellar evolution tracks are computed at solar, scaled-solar, or alpha-element… ▽ More We present the Stromlo Stellar Tracks, a set of stellar evolutionary tracks, computed by modifying the Modules for Experiments in Stellar Astrophysics (MESA) 1D stellar evolution package, to fit the Galactic Concordance abundances for hot ($T > 8000$ K) massive ($\geq 10M_\odot$) Main-Sequence (MS) stars. Until now, all stellar evolution tracks are computed at solar, scaled-solar, or alpha-element enhanced abundances, and none of these models correctly represent the Galactic Concordance abundances at different metallicities. This paper is the first implementation of Galactic Concordance abundances to the stellar evolution models. The Stromlo tracks cover massive stars ($10\leq M/M_\odot \leq 300$) with varying rotations ($v/v_{\rm crit} = 0.0, 0.2, 0.4$) and a finely sampled grid of metallicities ($-2.0 \leq {\rm [Z/H]} \leq +0.5$; $Δ{\rm [Z/H]} = 0.1$) evolved from the pre-main sequence to the end of $^{12}$Carbon burning. We find that the implementation of Galactic Concordance abundances is critical for the evolution of main-sequence, massive hot stars in order to estimate accurate stellar outputs (L, T$_{\rm eff}$, $g$), which, in turn, have a significant impact on determining the ionizing photon luminosity budgets. We additionally support prior findings of the importance that rotation plays on the evolution of massive stars and their ionizing budget. The evolutionary tracks for our Galactic Concordance abundance scaling provide a more empirically motivated approach than simple uniform abundance scaling with metallicity for the analysis of HII regions and have considerable implications in determining nebular emission lines and metallicity. Therefore, it is important to refine the existing stellar evolutionary models for comprehensive high-redshift extragalactic studies. The Stromlo tracks are publicly available to the astronomical community online. △ Less

Submitted 28 February, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 13 pages, 12 figures. Accepted for publication in ApJ. The Stromlo tracks are available online at https://sites.google.com/view/stromlotracks

Journal ref: ApJ, 908 (2021), 241

Spatially resolved direct method metallicity in a high-redshift analogue local galaxy: temperature structure impact on metallicity gradients

Authors: Alex J. Cameron, Tiantian Yuan, Michele Trenti, David C. Nicholls, Lisa J. Kewley

Abstract: We investigate how HII region temperature structure assumptions affect "direct-method" spatially-resolved metallicity observations using multispecies auroral lines in a galaxy from the SAMI Galaxy Survey. SAMI609396B, at redshift $z=0.018$, is a low-mass galaxy in a minor merger with intense star formation, analogous to conditions at high redshifts. We use three methods to derive direct metallicit… ▽ More We investigate how HII region temperature structure assumptions affect "direct-method" spatially-resolved metallicity observations using multispecies auroral lines in a galaxy from the SAMI Galaxy Survey. SAMI609396B, at redshift $z=0.018$, is a low-mass galaxy in a minor merger with intense star formation, analogous to conditions at high redshifts. We use three methods to derive direct metallicities and compare with strong-line diagnostics. The spatial metallicity trends show significant differences among the three direct methods. Our first method is based on the commonly used electron temperature $T_e$([OIII]) from the [OIII]$λ$4363 auroral line and a traditional $T_e$([OII]) -- $T_e$([OIII]) calibration. The second method applies a recent empirical correction to the O$^+$ abundance from the [OIII]/[OII] strong-line ratio. The third method infers the $T_e$([OII]) from the [SII]$λλ$4069,76 auroral lines. The first method favours a positive metallicity gradient along SAMI609396B, whereas the second and third methods yield flattened gradients. Strong-line diagnostics produce mostly flat gradients, albeit with unquantified contamination from shocked regions. We conclude that overlooked assumptions about the internal temperature structure of HII regions in the direct method can lead to large discrepancies in metallicity gradient studies. Our detailed analysis of SAMI609396B underlines that high-accuracy metallicity gradient measurements require a wide array of emission lines and improved spatial resolutions in order to properly constrain excitation sources, physical conditions, and temperature structures of the emitting gas. Integral-field spectroscopic studies with future facilities such as JWST/NIRSpec and ground-based ELTs will be crucial in minimising systematic effects on measured gradients in distant galaxies. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 21 pages, 14 figures, Accepted for Publication in MNRAS

Gas-phase metallicity gradients of TNG50 star-forming galaxies

Authors: Z. S. Hemler, Paul Torrey, Jia Qi, Lars Hernquist, Mark Vogelsberger, Xiangcheng Ma, Lisa J. Kewley, Dylan Nelson, Annalisa Pillepich, Rüdiger Pakmor, Federico Marinacci

Abstract: We present the radial gas-phase, mass-weighted metallicity profiles and gradients of the TNG50 star-forming galaxy population measured at redshifts $z=$ 0--3. We investigate the redshift evolution of gradients and examine relations between gradient steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity am… ▽ More We present the radial gas-phase, mass-weighted metallicity profiles and gradients of the TNG50 star-forming galaxy population measured at redshifts $z=$ 0--3. We investigate the redshift evolution of gradients and examine relations between gradient steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity among these negative gradients. We determine that the gradient steepness of all galaxies increases approximately monotonically with redshift at a roughly constant rate. This rate does not vary significantly with galaxy mass. We observe a weak negative correlation between gradient steepness and galaxy stellar mass at redshifts $z\leq2$. However, when we normalize gradients by a characteristic radius defined by the galactic star formation distribution, we find that these normalized gradients remain invariant with both stellar mass and redshift. We place our results in the context of previous simulations and show that TNG50 high-redshift gradients are steeper than those of models featuring burstier feedback, which may further highlight high-redshift gradients as important discriminators of galaxy formation models. We also find that redshift $z=0$ and $z=0.5$ TNG50 gradients are consistent with the gradients observed in galaxies at these redshifts, although the preference for flat gradients observed in redshift $z\gtrsim1$ galaxies is not present in TNG50. If future JWST and ELT observations validate these flat gradients, it may indicate a need for simulation models to implement more powerful radial gas mixing within the ISM, possibly via turbulence and/or stronger winds △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 24 pages, 14 figures, submitted to MNRAS

Oxygen Loss from Simulated Galaxies and the Metal Flow Main Sequence: Predicting the Dependence on Mass and Environment

Authors: Philip Taylor, Chiaki Kobayashi, Lisa J. Kewley

Abstract: We predict the mass fraction of oxygen lost from galaxies in a cosmological simulation as a function of stellar mass and environment at the present day. The distribution with stellar mass is bimodal, separating star-forming and quenched galaxies. The metallicity of gas and stars is self-consistently calculated using a chemical evolution model that includes supernovae type II and Ia, hypernovae, an… ▽ More We predict the mass fraction of oxygen lost from galaxies in a cosmological simulation as a function of stellar mass and environment at the present day. The distribution with stellar mass is bimodal, separating star-forming and quenched galaxies. The metallicity of gas and stars is self-consistently calculated using a chemical evolution model that includes supernovae type II and Ia, hypernovae, and asymptotic giant branch stars. The mass of oxygen lost from each galaxy is calculated by comparing the existing oxygen in gas and stars in the galaxy to the oxygen that should have been produced by the present-day population of stars. More massive galaxies are able to retain a greater fraction of their metals ($\sim 100$ per cent) than low-mass galaxies ($\sim 40 - 70$ per cent). As in the star formation main sequence, star-forming galaxies follow a tight relationship also in terms of oxygen mass lost -- a metal flow main sequence, ZFMS -- whereas massive quenched galaxies tend to have lost a greater fraction of oxygen (up to 20 per cent), due to AGN-driven winds. The amount of oxygen lost by satellite galaxies depends on the details of their interaction history, and those in richer groups tend to have lost a greater fraction of their oxygen. Observational estimates of metal retention in galaxies will provide a strong constraint on models of galaxy evolution. △ Less

Submitted 27 June, 2020; originally announced June 2020.

Comments: 10 pages, 7 figures, accepted for publication in MNRAS

On the origin of nitrogen at low metallicity

Authors: Arpita Roy, Michael A. Dopita, Mark R. Krumholz, Lisa J. Kewley, Ralph S. Sutherland, Alexander Heger

Abstract: Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment… ▽ More Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds. △ Less

Submitted 7 February, 2021; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: 19 pages, 14 figures, 4 tables, accepted for publication in MNRAS

Quantifying the effects of spatial resolution and noise on galaxy metallicity gradients

Authors: Ayan Acharyya, Mark R. Krumholz, Christoph Federrath, Lisa J. Kewley, Nathan J. Goldbaum, Rob Sharp

Abstract: Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galax… ▽ More Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 21 pages, 11 figures, 20 pages Supplementary Online Material provided with 10 additional figures, accepted for publication in MNRAS

Spinning Bar and a Star-formation Inefficient Repertoire: Turbulence in Hickson Compact Group NGC7674

Authors: Diane M. Salim, Katherine Alatalo, Christoph Federrath, Brent Groves, Lisa J. Kewley

Abstract: The physics regulating star formation (SF) in Hickson Compact Groups (HCG) has thus far been difficult to describe, due to their unique kinematic properties. In this study we expand upon previous works to devise a more physically meaningful SF relation able to better encompass the physics of these unique systems. We combine CO(1--0) data from the Combined Array from Research in Millimeter Astronom… ▽ More The physics regulating star formation (SF) in Hickson Compact Groups (HCG) has thus far been difficult to describe, due to their unique kinematic properties. In this study we expand upon previous works to devise a more physically meaningful SF relation able to better encompass the physics of these unique systems. We combine CO(1--0) data from the Combined Array from Research in Millimeter Astronomy (CARMA) to trace the column density of molecular gas $Σ_\mathrm{gas}$ and deep H$α$ imaging taken on the Southern Astrophysical Research (SOAR) Telescope tracing $Σ_\mathrm{SFR}$ to investigate star formation efficiency across face-on HCG, NGC7674. We find a lack of universality in star formation, with two distinct sequences present in the $Σ_\mathrm{gas}-Σ_\mathrm{SFR}$ plane; one for inside and one for outside the nucleus. We devise a SF relation based on the multi-freefall nature of gas and the critical density, which itself is dependent on the virial parameter $α_{\mathrm{vir}}$, the ratio of turbulent to gravitational energy. We find that our modified SF relation fits the data and describes the physics of this system well with the introduction of a virial parameter of about 5--10 across the galaxy. This $α_{\mathrm{vir}}$ leads to an order-of-magnitude reduction in SFR compared to $α_{\mathrm{vir}}\approx 1$ systems. △ Less

Submitted 23 March, 2020; originally announced March 2020.

Comments: 16 pages, 7 figures

A comparison of UV and optical metallicities in star-forming galaxies

Authors: Nell Byler, Lisa J Kewley, Jane R Rigby, Ayan Acharyya, Danielle A Berg, Matthew Bayliss, Keren Sharon, .

Abstract: Our ability to study the properties of the interstellar medium (ISM) in the earliest galaxies will rely on emission line diagnostics at rest-frame ultraviolet (UV) wavelengths. In this work, we identify metallicity-sensitive diagnostics using UV emission lines. We compare UV-derived metallicities with standard, well-established optical metallicities using a sample of galaxies with rest-frame UV an… ▽ More Our ability to study the properties of the interstellar medium (ISM) in the earliest galaxies will rely on emission line diagnostics at rest-frame ultraviolet (UV) wavelengths. In this work, we identify metallicity-sensitive diagnostics using UV emission lines. We compare UV-derived metallicities with standard, well-established optical metallicities using a sample of galaxies with rest-frame UV and optical spectroscopy. We find that the He2-O3C3 diagnostic (He II 1640 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, particularly at low metallicity (12+log(O/H) < 8), where stellar contributions are minimal. We find that the Si3-O3C3 diagnostic (Si III 1883 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, though with large scatter (0.2-0.3 dex), which we suggest is driven by variations in gas-phase abundances. We find that the C4-O3C3 diagnostic (C IV 1548,1550 / O III 1666 vs. O III 1666 / C III 1906,1909) correlates poorly with optically-derived metallicities. We discuss possible explanations for these discrepant metallicity determinations, including the hardness of the ionizing spectrum, contribution from stellar wind emission, and non-solar-scaled gas-phase abundances. Finally, we provide two new UV oxygen abundance diagnostics, calculated from polynomial fits to the model grid surface in the He2-O3C3 and Si3-O3C3 diagrams. △ Less

Submitted 15 March, 2020; originally announced March 2020.

Comments: ApJ accepted

MOSEL Survey: Tracking the Growth of Massive Galaxies at 2<z<4 using Kinematics and the IllustrisTNG Simulation

Authors: Anshu Gupta, Kim-Vy Tran, Jonathan Cohn, Leo Y. Alcorn, Tiantian Yuan, Vicente Rodriguez-Gomez, Anishya Harshan, Ben Forrest, Lisa J. Kewley, Karl Glazebrook, Caroline M. Straatman, Glenn G. Kacprzak, Themiya Nanayakkara, Ivo Labbé, Casey Papovich, Michael Cowley

Abstract: We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mas… ▽ More We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 19 pages, 9 figures, Accepted for publication in the Astrophysical Journal

ZFIRE: Measuring Electron Density with [OII] as a function of environment at z = 1.62

Authors: Anishya Harshan, Anshu Gupta, Kim-Vy Tran, Leo Y. Alcorn, Tiantian Yuan, Glenn G. Kacprzak, Themiya Nanayakkara, Karl Glazebrook, Lisa J. Kewley, Ivo Labb'e, Casey Papovich

Abstract: The global star formation rates (SFR) of galaxies at fixed stellar masses increase with redshift and are known to vary with environment unto z~2. We explore here whether the changes in the star formation rates can also apply to the electron densities of the inter-stellar medium (ISM) by measuring the [OII] (3727A/3729A) ratio for cluster and field galaxies at z~2. We measure a median electron dens… ▽ More The global star formation rates (SFR) of galaxies at fixed stellar masses increase with redshift and are known to vary with environment unto z~2. We explore here whether the changes in the star formation rates can also apply to the electron densities of the inter-stellar medium (ISM) by measuring the [OII] (3727A/3729A) ratio for cluster and field galaxies at z~2. We measure a median electron density of ne = 366+/-84 cm-3 for six galaxies (with 1-sigma scatter = 163 cm-3) in the UDS proton-cluster at z=1.62. We find that the median electron density of galaxies in the UDS photo-cluster environment is three times higher compared to the median electron density of field galaxies (ne = 113+/- 63 cm-3 and 1-sigma scatter = 79 cm-3) at comparable redshifts, stellar mass and SFR. However, we note that a sample of six photo-cluster galaxies is insufficient to reliably measure the electron density in the average porto-cluster environment at z~2. We conclude that the electron density increases with redshift in both cluster and field environments up to z~2 (ne = 30 +/- 1 cm-3 for z ~ 0 to ne =254+/- 76 cm-3 for z~1.5). We find tentative evidence (~2.6 sigma ) for a possible dependence of electron density on environment, but the results require confirmation with larger sample sizes. △ Less

Submitted 25 February, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: 18 Pages, 6 figures, 3 tables, Accepted for publication in Astrophysical Journal

What Drives the Redshift Evolution of Strong Emission Line Ratios?

Authors: Fuyan Bian, Lisa J. Kewley, Brent Groves, Michael A. Dopita

Abstract: We study the physical mechanisms that cause the offset between low-redshift and high-redshift galaxies on the [OIII]/H$β$ versus [NII]/H$α$ ``Baldwin, Phillips & Terlevich'' (BPT) diagram using a sample of local analogues of high-redshift galaxies. These high-redshift analogue galaxies are selected from the Sloan Digital Sky Survey. Located in the same region on the BPT diagram as the ultra-violet… ▽ More We study the physical mechanisms that cause the offset between low-redshift and high-redshift galaxies on the [OIII]/H$β$ versus [NII]/H$α$ ``Baldwin, Phillips & Terlevich'' (BPT) diagram using a sample of local analogues of high-redshift galaxies. These high-redshift analogue galaxies are selected from the Sloan Digital Sky Survey. Located in the same region on the BPT diagram as the ultra-violet selected galaxies at $z\sim2$, these high-redshift analogue galaxies provide an ideal local benchmark to study the offset between the local and high-redshift galaxies on the BPT diagram. We compare the nitrogen-to-oxygen ratio (N/O), the shape of the ionising radiation field, and ionisation parameters between the high-redshift analogues and a sample of local reference galaxies. The higher ionisation parameter in the high-redshift analogues is the dominant physical mechanism driving the BPT offset from low- to high-redshift, particularly at high {\nii/\ha}. Furthermore, the N/O ratio enhancement also plays a minor role to cause the BPT offset. However, the shape of the ionising radiation field is unlikely to cause the BPT offset because the high-redshift analogues have a similar hard ionising radiation field as local reference galaxies. This hard radiation field cannot be produced by the current standard stellar synthesis models. The stellar rotation and binarity may help solve the discrepancy. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Comments: 7 pages, 5 figures. accepted for publication in MNRAS