Showing 1–50 of 107 results for author: Roman-Duval, J

JWST Whirlpool Galaxy Treasury: Mid-Infrared Emission in M51 and its Relation to Gas Column and Star Formation

Authors: Mansi Padave, Karin M. Sandstrom, Daniel A. Dale, Adam K. Leroy, Eric W. Koch, Tony D. Weinbeck, Angela Adamo, Jessica Sutter, Lindsey Hands, Torsten Böker, Martha L. Boyer, Ryan Chown, Bruce T. Draine, Ilse de Looze, Brandt A. L. Gaches, Simon C. O. Glover, Dario Colombo, Robert C. Kennicutt Jr., Hannah Koziol, Ralf S. Klessen, Sean T. Linden, Desika Narayanan, Alex Pedrini, Julia Roman-Duval, Eva Schinnerer , et al. (3 additional authors not shown)

Abstract: Using JWST/MIRI imaging of M51 in eight broadband filters, we investigate correlations of mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) and dust continuum with molecular, atomic, and ionized gas traced by CO(1-0), HI, and Pa-alpha, respectively. In molecular gas-dominated regions, PAH-dominated filters (F560W, F770W, F1130W, F1280W) exhibit near-linear correlations with CO(1-0… ▽ More Using JWST/MIRI imaging of M51 in eight broadband filters, we investigate correlations of mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) and dust continuum with molecular, atomic, and ionized gas traced by CO(1-0), HI, and Pa-alpha, respectively. In molecular gas-dominated regions, PAH-dominated filters (F560W, F770W, F1130W, F1280W) exhibit near-linear correlations with CO(1-0) at 40 pc scale, indicating that PAHs are well-mixed with gas and experience relatively constant radiation field intensities. The F1500W, F1800W, and F2100W dust continuum-dominated filters show shallower slopes with CO(1-0), reflecting contributions from star-forming regions with high radiation field intensities. This is reinforced by the near-linear scaling between F2100W and Pa-alpha. PAH-dominated bands do not show this linear trend with Pa-alpha, likely due to their destruction in ionized regions. F1000W behaves similarly to PAH bands in its correlations with CO(1-0) and Pa-alpha. Modeling mid-infrared emission with an empirical decomposition into gas- and star-formation-associated components shows that PAH-dominated filters receive comparable contributions from both, while the relative contribution associated with the Pa-alpha template increases toward longer wavelengths, reaching $\sim$75% in F2100W. These results demonstrate that mid-infrared simultaneously traces the gas column and star formation, but with a systematic wavelength-dependent shift in what drives the correlations: PAHs being more gas-tracing and dust-continuum reflecting star formation. Lastly, considering both HI and H$_2$ at 440 pc resolution, we find a tight, linear relation between $Σ_{HI+H_2}$ and PAH-dominated filters. Although most of our coverage is in H$_2$-dominated regions, we note similar observations with HI, suggesting that PAHs are also well-mixed with atomic gas. △ Less

Submitted 17 August, 2026; originally announced August 2026.

Comments: (Accepted for publication in ApJ)

Observing Co-Located Neutral and Ionized Gas-Phase Iron Depletion in the Magellanic Clouds

Authors: Yun Qi Li, Jessica K. Werk, Caleb R. Choban, Julia Roman-Duval, Kirill Tchernyshyov, J. Xavier Prochaska, Doyeon A. Kim, Arianna S. Long

Abstract: Depletion is the observed phenomenon where gas-phase elemental abundances are reduced through accretion onto dust grains. We measure neutral gas-phase elemental abundances (S, Fe) in the Magellanic Clouds along 33 sightlines using high-resolution UV spectroscopy (HST/COS and HST/STIS), and compare them to ionized gas-phase abundances (S, Fe) adopted from the literature for six co-located H\,\texts… ▽ More Depletion is the observed phenomenon where gas-phase elemental abundances are reduced through accretion onto dust grains. We measure neutral gas-phase elemental abundances (S, Fe) in the Magellanic Clouds along 33 sightlines using high-resolution UV spectroscopy (HST/COS and HST/STIS), and compare them to ionized gas-phase abundances (S, Fe) adopted from the literature for six co-located H\,\textsc{ii} regions (with the furthest separation of $\lesssim3'$, 50 pc). Comparing S abundances show that S is minimally depleted in the H\,\textsc{ii} regions and surrounding diffuse ISM. However, we find that the gas-phase Fe abundances in H\,\textsc{ii} regions can be lower than those of the neighboring neutral ISM by 0.3 to 2 dex. This difference is likely an offset in the amount of Fe depleted into dust grains. As accretion of gas-phase Fe is likely not effective at the temperatures of the H\,\textsc{ii} regions, Fe depletion into solid form would have occurred in the dense atomic or molecular clouds prior to star formation. Stronger depletion in the H\,\textsc{ii} regions shows that Fe-bearing grains survive destruction in the first few million years following ionization. Our observations highlight that Fe depletion in H\,\textsc{ii} regions can be a useful tracer of Fe depletion in dense molecular clouds, which are challenging to observe directly via UV absorption. △ Less

Submitted 12 August, 2026; originally announced August 2026.

GLOW II: A Census of Oxygen in Low-Mass Galaxies

Authors: Kristen B. W. McQuinn, O. Grace Telford, Andrey V. Kravtsov, Joseph Burchett, Roger Cohen, Laura Hunter, Evan D. Skillman, Julia Roman-Duval, Danielle A. Berg, John M. Cannon, John Chisholm, Julianne J. Dalcanton, Andrew E. Dolphin, Liese van Zee, Benjamin Williams

Abstract: Oxygen is forged by stars and redistributed through galaxies by feedback-driven outflows, leaving a record of star formation and the baryon cycle imprinted on its present-day abundance and distribution. The Galaxies Losing Oxygen via Winds (GLOW) project quantifies the production, distribution, and retention of oxygen in 37 low-mass, low-metallicity, gas-rich galaxies in the nearby universe (D<6 M… ▽ More Oxygen is forged by stars and redistributed through galaxies by feedback-driven outflows, leaving a record of star formation and the baryon cycle imprinted on its present-day abundance and distribution. The Galaxies Losing Oxygen via Winds (GLOW) project quantifies the production, distribution, and retention of oxygen in 37 low-mass, low-metallicity, gas-rich galaxies in the nearby universe (D<6 Mpc) spanning a critical stellar mass range (10^6.5<Mstar<10^9.5) where feedback is expected to strongly shape galaxy evolution. We incorporate literature results to extend the analysis over 10^5<Mstar<10^11.5. We couple star formation histories derived from resolved stars with measurements of ISM metallicity and gas mass, and compares the oxygen budgets with empirical constraints in the circumgalactic medium (CGM). We find low-mass galaxies retain between 4-25% (mean 11%) of their oxygen produced by stellar nucleosynthesis, and nearly all retained oxygen residing in the gas. Galaxies exhibit increased scatter to higher oxygen retention fractions when residing in denser environments (Theta_5>0.5) and within 1 Mpc of a massive galaxy (Mstar>10^9 Msun), indicating environment likely affects the amount of oxygen retained, recycled, or accreted to galaxies. Contrary to expectations, oxygen retention does not correlate with position on the mass-metallicity relation. Although ionized oxygen is detected in the CGM, it remains unclear whether most missing oxygen resides there or has been expelled entirely. Hydrodynamical simulations, despite successfully reproducing the mass-metallicity relation, predict much higher retention fractions than measured in low-mass galaxies. Simple modeling indicates that wind mass-loading and outflow metallicity govern oxygen retention, and that low-mass galaxies accrete less baryonic material relative to the cosmic baryon fraction. △ Less

Submitted 13 August, 2026; v1 submitted 7 August, 2026; originally announced August 2026.

Comments: 31 pages, 12 figures, 2 tables

Scylla: Observational Evidence for an Order of Magnitude in Dust Mass Opacity Evolution with ISM Density in the Large Magellanic Cloud

Authors: Christina W. Lindberg, Christopher J. R. Clark, Claire E. Murray, Julia Roman-Duval, Caroline Bot, Yumi Choi, Roger E. Cohen, Steven R. Goldman, Karl D. Gordon, Kristen B. W. McQuinn, Elizabeth Tarantino, Benjamin F. Williams, Petia Yanchulova Merica-Jones, Catherine Zucker

Abstract: The emissivity of dust is known to vary greatly with radiative environment, density, grain chemistry, and geometry. Discrepancies between dust mass surface densities derived from far-infrared (FIR) emission and visible extinction persist across and within galaxies in the local Universe. Here, we use new extinction and emission measurements towards the LMC to show that this discrepancy is driven by… ▽ More The emissivity of dust is known to vary greatly with radiative environment, density, grain chemistry, and geometry. Discrepancies between dust mass surface densities derived from far-infrared (FIR) emission and visible extinction persist across and within galaxies in the local Universe. Here, we use new extinction and emission measurements towards the LMC to show that this discrepancy is driven by the dust mass opacity evolving with the intrinsic density of the ISM, and that the ratio between FIR and optical dust mass opacity varies with gas surface density. These new findings imply that the dust mass opacity in the FIR could increase by nearly an order of magnitude (e.g., $κ_{160} = 0.3 - 6\ m^2\ kg^{-2}$) across over an order of magnitude of total hydrogen surface density $(Σ_H = 4 - 100 M_{\odot}\ pc^{-2})$, corroborating previous theoretical models for dust mass opacity evolution in the FIR, and providing new implications for emission-based dust mass estimates. △ Less

Submitted 5 August, 2026; originally announced August 2026.

Comments: accepted for publication by ApJ Letters

COS2035: Extending COS/FUV Operations Through the 2030s

Authors: Marc Rafelski, David Sahnow, Christian I. Johnson, Bethan James, Svea Hernandez, John Debes, Beverly Serrano, Kate Davis, Serge Dieterich, Van Dixon, Leonardo Dos Santos, Travis Fischer, Elaine Frazer, David French, Mark Giuliano, Joshua Goldberg, Sten Hasselquist, Jacqueline Hernandez, Nick Indriolo, Mike Kelly, Olivia Lupie, Lauren Miller, Diego Mundo, Anna Payne, Karla Peterson , et al. (7 additional authors not shown)

Abstract: The far-ultraviolet (FUV) detector of the Cosmic Origins Spectrograph (COS) accumulates gain sag where photons land, and without continued mitigation this degradation would render the most used modes unusable. To extend COS FUV operations through the 2030s, the COS team developed the COS2035 strategy, which builds on the existing COS2025 rules with four technical breakthroughs and two new usage po… ▽ More The far-ultraviolet (FUV) detector of the Cosmic Origins Spectrograph (COS) accumulates gain sag where photons land, and without continued mitigation this degradation would render the most used modes unusable. To extend COS FUV operations through the 2030s, the COS team developed the COS2035 strategy, which builds on the existing COS2025 rules with four technical breakthroughs and two new usage policies. First, SPLIT-wavecals decouple wavelength calibration from science exposures and open detector real estate above the Pt-Ne lamp light leak. Second, a hybrid lifetime position (LP) architecture allows different gratings to operate at different LPs simultaneously. Third, the LP-infinity framework removes the dependence on the eight-LP limit in the COS flight software, supported by a new table-based APT and TRANS rules architecture. Fourth, a revised gain-sag flagging method evaluates integrated column count loss against the maximum achievable signal-to-noise (S/N) per mode. The two new usage policies cap per-target S/N at the maximum achievable value set by fixed-pattern noise, and limit any single program to 2\% of the lifetime at any single LP. With LP7 and LP10 enabled in Cycle 33 and LP11 and LP12 in active commissioning for Cycles 34 and 35, the COS2035 strategy positions the FUV channel for continued high productivity through the 2030s. △ Less

Submitted 15 July, 2026; originally announced July 2026.

Comments: 23 pages, 5 figures

Report number: Instrument Science Report COS 2026-04

Dust in the very metal-poor galaxy Sextans A with JWST. I: Characterizing the evolved stellar population of Sextans A based on JWST observations and stellar evolution models

Authors: C. Gavetti, F. Dell'Agli, E. Tarantino, M. L. Boyer, I. McDonald, J. Th. Van Loon, D. A. Garcia-Hernandez, M. A. T. Groenewegen, A. Nanni, J. A. D. L. Blommaert, R. D. Gehrz, L. M. Gerlach, S. Goldman, M. Marengo, K. B. W. McQuinn, J. M. Oliveira, J. Roman-Duval, R. Sahai, E. D. Skillman, B. F. Williams, A. Javadi, O. C. Jones, F. Kemper, F. La Franca, G. C Sloan

Abstract: The nearby star-forming dwarf galaxy Sextans A offers a unique window into galaxy evolution in the early Universe, owing to its extremely low metallicity (about 1-7% Zsun). Recent JWST imaging of Sextans A spanning 1-21 micron enables a detailed characterization of its dusty stellar populations and interstellar medium. In this work, we compare the observed JWST color-magnitude distributions of evo… ▽ More The nearby star-forming dwarf galaxy Sextans A offers a unique window into galaxy evolution in the early Universe, owing to its extremely low metallicity (about 1-7% Zsun). Recent JWST imaging of Sextans A spanning 1-21 micron enables a detailed characterization of its dusty stellar populations and interstellar medium. In this work, we compare the observed JWST color-magnitude distributions of evolved stars with stellar evolution and dust-formation models to characterize the properties of the asymptotic giant branch (AGB) population, including progenitor mass, formation epoch, metallicity, and dust production. Evolutionary tracks for 0.8-7 Msun stars with metallicity Z=10^-3 provide good agreement with the overall distribution of AGB stars in Sextans A. More than 90% of the AGB population occupies a nearly vertical sequence in the color-magnitude diagrams, corresponding to stars spanning a wide range of masses and ages but exhibiting little or no circumstellar dust. This sequence appears to be dominated by oxygen-rich (M-type) AGB stars and reveals that the F444W flux is a robust luminosity diagnostic. A small subset of sources displays strong infrared excesses and is dominated by carbon stars descending from 1.25-1.5 Msun progenitors that formed about 2-3 Gyr ago and are currently in the final AGB phases. Their MIRI colors imply very low metallicities, consistent with estimates from the red giant branch morphology (about 1-2% Zsun). Finally, we show that the JWST/NIRCam F277W-F444W color serves as an effective proxy for the dust production rate, with models predicting rates up to 10^-7 Msun/yr for the reddest sources in Sextans A. △ Less

Submitted 5 June, 2026; originally announced June 2026.

Comments: 18 pages, 4 figures, accepted for publication in ApJ

The Panchromatic Hubble Andromeda Southern Treasury (PHAST). II. The Spatially Resolved Recent Star Formation History in M31

Authors: Tobin M. Wainer, Benjamin F. Williams, Zhuo Chen, Margaret Lazzarini, Julianne J. Dalcanton, Eric F. Bell, Kameron Goold, Andrew Dolphin, Meredith J. Durbin, Stefany L. Fabian Dubón, Karoline M. Gilbert, Puragra Guhathakurta, Francois Hammer, L. Clifton Johnson, Eric W. Koch, Kristen. B. W. McQuinn, Ekta Patel, Vaishnav V. Rao, Julia Roman-Duval, Adam Smercina, Debby Tran, Daniel R. Weisz

Abstract: We use Hubble Space Telescope optical imaging from the Panchromatic Hubble Andromeda Southern Treasury (PHAST) to measure the spatially resolved recent star formation history (SFH) across the southern disk of M31. We fit color-magnitude diagrams (CMDs) of over 6500 individual 0.01 kpc$^2$ regions to measure SFHs over the last $\sim$500 Myr. The resulting maps show coherent structure that traces th… ▽ More We use Hubble Space Telescope optical imaging from the Panchromatic Hubble Andromeda Southern Treasury (PHAST) to measure the spatially resolved recent star formation history (SFH) across the southern disk of M31. We fit color-magnitude diagrams (CMDs) of over 6500 individual 0.01 kpc$^2$ regions to measure SFHs over the last $\sim$500 Myr. The resulting maps show coherent structure that traces the ringed morphology of the disk. We find a clear global decline in the recent SFR, with a pronounced drop in the last $\sim$40 Myr that is most evident in the region closest to M32. Combining PHAST and PHAT measurements, we now cover two thirds of M31's star-forming disk with homogeneous SFHs, yielding the highest-resolution spatially resolved SFHs of M31. Inside the joint footprint, we measure mean SFRs of $0.445 \pm0.006$ M$_\odot$ yr$^{-1}$ over the last 100 Myr and $0.285 \pm 0.014$ M$_\odot$ yr$^{-1}$ over the last 20 Myr, implying total disk SFRs of $\sim$0.67 and $\sim$0.43 M$_\odot$ yr$^{-1}$, respectively. The observed decline is interpreted as the late stage of a multi-Gyr wind-down from a previously more active state. Because recent star formation in M31 is concentrated primarily in the rings, the global decline is driven mainly by decreasing activity within those features. We also compare the CMD-based SFR surface densities to those inferred from FUV+24 $μ$m prescriptions and find that the FUV-based calibration underestimates the CMD-based 100 Myr average by a factor of $\sim$2.1. However, the PHAST SFHs produce a synthetic GALEX FUV image that agrees well with observations, indicating that the CMD-derived SFHs provide an accurate description of recent star formation. The mismatch with the FUV+24 $μ$m estimates underscores that tracers implicitly averaged over $\sim$100 Myr are not reliable when the recent SFR is evolving. △ Less

Submitted 12 May, 2026; originally announced May 2026.

Comments: Accepted to ApJ. 38 Pages, 15 figures, and 1 table. 2 Animations are available with the journal version

Journal ref: Tobin M. Wainer et al 2026 ApJ 1006 195

Scylla VI: Parsec-Scale Dust Extinction Maps in the SMC and LMC

Authors: Christina W. Lindberg, Claire E. Murray, Christopher J. R. Clark, Caroline Bot, Clare Burhenne, Yumi Choi, Roger E. Cohen, Steven R. Goldman, Karl D. Gordon, Kristen B. W. McQuinn, Julia Roman-Duval, Karin M. Sandstrom, Edward F. Schlafly, Elizabeth Tarantino, Benjamin F. Williams, Petia Yanchulova Merica-Jones, Catherine Zucker

Abstract: We present a novel methodology for mapping dust extinction in nearby galaxies at parsec-scale resolution. We apply it to HST 68 fields within the Small and Large Magellanic Clouds (23 fields in the SMC and 45 fields in the LMC) using multi-band HST photometry from the Scylla and METAL surveys. Our technique leverages \textit{kriging}, a geostatistical interpolation method built on the principles o… ▽ More We present a novel methodology for mapping dust extinction in nearby galaxies at parsec-scale resolution. We apply it to HST 68 fields within the Small and Large Magellanic Clouds (23 fields in the SMC and 45 fields in the LMC) using multi-band HST photometry from the Scylla and METAL surveys. Our technique leverages \textit{kriging}, a geostatistical interpolation method built on the principles of Gaussian Process regression, combined with Gaussian mixture modeling to statistically isolate background stellar sources and account for line-of-sight depth effects. 3D dust simulations demonstrate the method's capability to recover column densities to an accuracy of $A_V \approx 0.1$ mag in fields with at least 1000 sources. The resulting $4^{\prime\prime}$ resolution ($\sim1$-pc) dust maps reveal detailed structure and strong spatial correlation with ancillary ISM tracers, especially in star-forming regions like 30 Doradus. Global extinction of total column densities follows log-normal profiles in both galaxies, with the SMC exhibiting slightly higher mean extinction ($e^μ=0.47$ mag) than the broader LMC ($e^μ=0.43$ mag), likely due to significant line-of-sight depths. We find systematic offsets between dust mass surface densities ($Σ_{D}$) derived from extinction versus FIR emission in both galaxies, with $Σ_{D, FIR}/Σ_{D, A_V}$ ratios ranging from $0.6-1.8$. This work provides the highest-resolution dust extinction maps in SMC and LMC to date, which offer a vital independent benchmark for constraining dust emissivity, $\text{CO}$-dark gas fractions, and the multi-scale structure of the ISM in low-metallicity environments. △ Less

Submitted 7 May, 2026; originally announced May 2026.

Comments: 28 pages, 14 figures, accepted to ApJ

The Hubble Advanced Spectral Product (HASP) Program

Authors: John Debes, Ravi Sankrit, Travis Fischer, Elaine Frazer, Alec Hirschauer, Kate Rowlands, Matthew Burger, Robert Swaters, Robert Jedrzejewski, Sierra Gomez, Leonardo Dos Santos, Svea Hernandez, Lauren Miller, Anna Payne, Marc Rafelski, Thomas Wevers, Sara Anderson, Tom Bair, Kathryn Bello, Joleen Carlberg, Brian Charlow, Andrew Cortese, Nadia Dencheva, Tracy Ellis, Ben Falk , et al. (17 additional authors not shown)

Abstract: The Hubble Advanced Spectral Products (HASP) program is designed to robustly coadd Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS) spectra within the Mikulski Archive for Space Telescopes (MAST) in an automated fashion such that coadds are available for new data or archival data with updated calibrations. For each target within a visit or program, HASP employs a m… ▽ More The Hubble Advanced Spectral Products (HASP) program is designed to robustly coadd Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS) spectra within the Mikulski Archive for Space Telescopes (MAST) in an automated fashion such that coadds are available for new data or archival data with updated calibrations. For each target within a visit or program, HASP employs a meticulous multi-stage filtering process to ensure data quality and creates coadded products for all central wavelengths (CENWAVEs) within specific gratings, as well as combined products using different gratings and instruments. The project also emphasizes making the code accessible to the user community for custom coaddition. As calibrations improve and new data are added to the archive, HASP products are re-created automatically so that they represent the best reduction of a given visit or program. Automated coadditions like those achieved by HASP can significantly enhance the combination of different CENWAVES, increase signal-to-noise ratios, and increase wavelength coverage. These properties make HASP a vital resource for astronomers using archival spectroscopic data from HST. △ Less

Submitted 7 May, 2026; originally announced May 2026.

Comments: 30 pages, 8 figures

Report number: Instrument Science Report COS 2024-01

Ashes of FIRE: Modeling Dust Grain Size Evolution in the Local Group with FIRE

Authors: Caleb R. Choban, Samir Salim, Dušan Kereš, Julia Roman-Duval, Karin M. Sandstrom

Abstract: We introduce a new, discretized grain size evolution model, incorporated into the GIZMO code and coupled with FIRE-3 stellar feedback and ISM physics, to investigate variations in dust abundance, chemical composition, and grain sizes observed in the Local Group. This model tracks the size evolution of specific dust species, and includes stellar production of dust, dust growth through gas-phase met… ▽ More We introduce a new, discretized grain size evolution model, incorporated into the GIZMO code and coupled with FIRE-3 stellar feedback and ISM physics, to investigate variations in dust abundance, chemical composition, and grain sizes observed in the Local Group. This model tracks the size evolution of specific dust species, and includes stellar production of dust, dust growth through gas-phase metal accretion, dust destruction by sputtering, SNe shocks, and astration, grain-grain collisional shattering and coagulation, and turbulent dust diffusion. Using idealized galaxy simulations, we test the dependence of MW dust properties on variations in each dust process and find that our model uniquely predicts a bimodal grain size distribution. This bimodality is due to our simulation's ability to resolve each dust process and where they occur in the ISM, unlike other works. We find that Local Group dust abundances are determined by dust growth and destruction, with little dependence on coagulation or shattering, explaining why models that do not include these processes can match abundance observations. We also find that variations in Local Group extinction curve slopes are determined by coagulation, with inefficient coagulation leading to steeper slopes. However, inefficient coagulation also results in stronger extinction curve bumps, which are not observed. We also do not predict a population of very small (${<}1$ nm) carbonaceous grains, required for MIR emission features, due to their rapid growth by accretion. These results highlight the possible necessity of ``top-down'' PAH formation from preexisting grains as a means to inhibit carbonaceous dust growth. △ Less

Submitted 9 March, 2026; originally announced March 2026.

Comments: 29 pages, 17 figures, submitted to MNRAS. Comments welcomed

Journal ref: Mon Not R Astron Soc (2026)

JWST Captures Growth of Aromatic Hydrocarbon Dust Particles in the Extremely Metal-poor Galaxy Sextans A

Authors: Elizabeth J. Tarantino, Julia Roman-Duval, Karin M. Sandstrom, J. -D. T. Smith, Cory M. Whitcomb, Bruce T. Draine, Martha L. Boyer, Jérémy Chastenet, Ryan Chown, Christopher J. R. Clark, Karl D. Gordon, Brandon S. Hensley, Thomas S. -Y. Lai, Christina W. Lindberg, Kristen B. W. McQuinn, Max J. B. Newman, O. Grace Telford, Dries Van De Putte, Benjamin F. Williams

Abstract: The mid-infrared spectrum of star-forming, high metallicity galaxies is dominated by emission features from aromatic and aliphatic bonds in small carbonaceous dust grains, often referred to as polycyclic aromatic hydrocarbons (PAHs). In metal-poor galaxies, the abundance of PAHs relative to the total dust sharply declines, but the origin of this deficit is unknown. We present JWST observations tha… ▽ More The mid-infrared spectrum of star-forming, high metallicity galaxies is dominated by emission features from aromatic and aliphatic bonds in small carbonaceous dust grains, often referred to as polycyclic aromatic hydrocarbons (PAHs). In metal-poor galaxies, the abundance of PAHs relative to the total dust sharply declines, but the origin of this deficit is unknown. We present JWST observations that detect and resolve emission from PAHs in the 7% Solar metallicity galaxy Sextans A, representing the lowest metallicity detection of PAH emission to date. In contrast to higher metallicity galaxies, the clumps of PAH emission are compact (0.5-1.5'' or 3-10 pc), which explains why PAH emission evaded detection by lower resolution instruments like Spitzer. Ratios between the 3.3, 7.7, and 11.3 $μ$m PAH features indicate that the PAH grains in Sextans A are small and neutral, with no evidence of significant processing from the hard radiation fields within the galaxy. These results favor inhibited grain growth over enhanced destruction as the origin of the low PAH abundance in Sextans A. The compact clumps of PAH emission are likely active sites of in-situ PAH growth within a dense, well-shielded phase of the interstellar medium. Our results show that PAHs can form and survive in extremely metal-poor environments common early in the evolution of the Universe. △ Less

Submitted 3 December, 2025; originally announced December 2025.

Comments: Under review

Ashes of Creation: JWST Uncovers Silicate Dust in Massive Star Clusters

Authors: Daniel Maschmann, Bradley C. Whitmore, David A. Thilker, Ivan Gerasimov, Simon C. O. Glover, B. T. Draine, Bret Lehmer, Varun Bajaj, Sumit Sarbadhicary, Médéric Boquien, G. C. Sloan, Tony D. Weinbeck, Daniel A. Dale, Kiana Henny, Kirsten L. Larson, M. Jimena Rodríguez, Robert Kennicutt, Amirnezam Amiri, Ashley. T. Barnes, Torsten Böker, Martha Boyer, Daizhong Liu, Oleg V. Egorov, Hwihyun Kim, Ralf S. Klessen , et al. (13 additional authors not shown)

Abstract: Dust production is a fundamental aspect of the baryonic cycle of star formation. It is known that dust is injected into the interstellar medium during early star formation by supernovae and later on by evolved stars. From individual objects, these mechanisms are well understood, but the overall dust production in star clusters at different evolutionary stages is still challenging to quantify. We p… ▽ More Dust production is a fundamental aspect of the baryonic cycle of star formation. It is known that dust is injected into the interstellar medium during early star formation by supernovae and later on by evolved stars. From individual objects, these mechanisms are well understood, but the overall dust production in star clusters at different evolutionary stages is still challenging to quantify. We present 22 massive (> 105M$_{\odot}$) extra galactic star clusters with ages between 3 and 100 Myr exhibiting a compact dust morphology seen with JWST-MIRI. We only find PAH features associated with one star cluster and nineteen have already cleared themselves from their natal dust. Their main characteristic is a significant enhancement at 10$μ$m, which is likely due to silicate emission and cannot be explained by ionized gas. We discuss several possible explanations including dust production from evolved stars such as red super giants, more exotic star types like yellow hypergiants and luminous blue variable stars. Stochastic dust injection from supernovae or a single supernova in dense gas can also create significant silicate emission. However, for this scenario secondary tracers such as a X-ray signal are expected which we only observe in three star clusters. We find the most luminous 10$μ$m emitter to be the three most massive star clusters (> 106M$_{\odot}$) which is at least a magnitude stronger than any known stellar sources indicating a rare mechanism that only appears at extreme masses and a short lifetime. △ Less

Submitted 25 November, 2025; originally announced November 2025.

Comments: 38 Pages, 13Figures, submitted to APJ

The Local Group L-band Survey: Probing Cold Atomic Gas in IC10 with Neutral Hydrogen Absorption

Authors: Ioana A. Stelea, Snezana Stanimirovic, Nickolas M. Pingel, Hongxing Chen, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Chang-Goo Kim, Alberto D. Bolatto, Julianne J. Dalcanton, Michael P. Busch, Harrisen Corbould, J. R. Dawson, Cosima Eibensteiner, Amanda Kepley, Melanie Krips, Claire E. Murray, Julia Roman-Duval, Daniel R. Rybarczyk, Evan D. Skillman, Elizabeth Tarantino, Vicente Villanueva, Thomas G. Williams

Abstract: We present the first localized detections of the cold neutral medium (CNM) in IC10, offering a rare view of dense atomic gas in a low-metallicity (0.27 solar metallicity) dwarf galaxy. As a low-metallicity starburst, IC10's interstellar medium conditions could reflect small-scale physical conditions that mirror those of early galaxies, providing a unique window into the heating and cooling process… ▽ More We present the first localized detections of the cold neutral medium (CNM) in IC10, offering a rare view of dense atomic gas in a low-metallicity (0.27 solar metallicity) dwarf galaxy. As a low-metallicity starburst, IC10's interstellar medium conditions could reflect small-scale physical conditions that mirror those of early galaxies, providing a unique window into the heating and cooling processes that shaped the interstellar medium in early-Universe environments. Leveraging the high angular (<5'' ~ 15pc) and spectral (0.4 km/s) resolution of the Local Group L-band Survey, we searched for HI absorption against nine continuum radio sources and detected absorption along three sightlines corresponding to internal radio emission sources within IC10. Using Gaussian decomposition and radiative transfer, we characterize the CNM, deriving spin temperatures of ~30-55 K, column densities of (0.6-3.0)x 1$0^{21}$ cm$^{-2}$, cold HI fractions of ~ 21-37%, and line widths of ~ 5.6-13.6 km/s. For each individual detection of HI absorption, we find corresponding molecular emission from 12CO (J=1-0), HCO+ (J=1-0), and HCN (J=1-0) at similar velocities and with comparable linewidths, indicating a well-mixed cold atomic and molecular medium. In IC10, the CNM shows a clear kinematic connection to the high-density ISM, implying a stronger dynamical coupling with molecular gas than in the Milky Way, in line with expectations for low-metallicity environments. At the ~ 15 pc scales probed by slightly extended HII regions in IC10, unresolved CNM clouds likely contribute to line blending, so the observed broad HI linewidths may partly reflect spatial and kinematic averaging. △ Less

Submitted 10 December, 2025; v1 submitted 13 November, 2025; originally announced November 2025.

Comments: 18 pages, 8 figures, accepted for publication in ApJ

Scylla V: Constraints on the spatial and temporal distribution of bursts and the interaction history of the Magellanic Clouds from their resolved stellar populations

Authors: Clare Burhenne, Kristen B. W. McQuinn, Roger E. Cohen, Claire E. Murray, Ekta Patel, Benjamin F. Williams, Christina W. Lindberg, Petia Yanchulova Merica-Jones, Karl D. Gordon, Yumi Choi, Andrew E. Dolphin, Julia C. Roman-Duval

Abstract: We measure the star formation histories (SFHs) from the Scylla survey in approximately 98,000 pc^2 and 75,000 pc^2 of the SMC and LMC, respectively, using deep Hubble Space Telescope imaging (80% complete to more than 1 mag below the ancient main-sequence turnoff, 25.1 and 26.0 mag in F475W and F814W) from 74 pointings. We group the fields into eight sub-regions in the SMC and seven in the LMC. We… ▽ More We measure the star formation histories (SFHs) from the Scylla survey in approximately 98,000 pc^2 and 75,000 pc^2 of the SMC and LMC, respectively, using deep Hubble Space Telescope imaging (80% complete to more than 1 mag below the ancient main-sequence turnoff, 25.1 and 26.0 mag in F475W and F814W) from 74 pointings. We group the fields into eight sub-regions in the SMC and seven in the LMC. We use the birth rate parameter to identify bursts of star formation and measure their properties in each sub-region. Our methodology provides a standardized framework for burst identification and reveals both broad and fine burst characteristics. We identify global and local bursts, defined as those occurring in at least half or less than half of a galaxy's sub-regions, respectively. In the SMC we find two global (about 5 and 1.5 Gyr ago) and one local burst (about 3 Gyr ago). In the LMC we find one global burst (about 3 Gyr ago). Comparing these findings with dynamical models of the LMC and SMC orbital histories, we find that when models predict a shared dynamical trigger for bursts across both galaxies, the burst begins earlier in the SMC with a greater enhancement in star formation rate than in the LMC. Finally, using age-metallicity relations (AMRs) and cumulative SFHs, we report that the Wing/Bridge region in the SMC resembles the southwestern LMC both chemically and in stellar mass assembly over the last about 7 Gyr, possibly due to stellar material stripped from the LMC during their last interaction. △ Less

Submitted 4 November, 2025; originally announced November 2025.

Comments: 32 pages, 15 figures

The Metallicity Dependence of PAH Emission in Galaxies II: Insights from JWST/NIRCam Imaging of the Smallest Dust Grains in M101

Authors: Cory M. Whitcomb, J. -D. T. Smith, Elizabeth Tarantino, Karin Sandstrom, Thomas S. -Y. Lai, Lee Armus, Alberto Bolatto, Martha Boyer, Daniel A. Dale, Bruce T. Draine, Brandon S. Hensley, Desika Narayanan, Julia Roman-Duval, Evan D. Skillman

Abstract: We explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity using JWST/NIRCam imaging of the nearby galaxy M101, covering regions from solar metallicity (Z$_{\odot}$) down to 0.4 Z$_{\odot}$. These maps are used to trace the radial evolution of the shortest-wavelength PAH feature at 3.3 $μ$m, which is emitted preferentially by the sma… ▽ More We explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity using JWST/NIRCam imaging of the nearby galaxy M101, covering regions from solar metallicity (Z$_{\odot}$) down to 0.4 Z$_{\odot}$. These maps are used to trace the radial evolution of the shortest-wavelength PAH feature at 3.3 $μ$m, which is emitted preferentially by the smallest PAHs ($<100$ carbon atoms). The fractional contribution of PAH 3.3 $μ$m to the total PAH luminosity ($Σ$PAH) increases by 3x as metallicity declines, rising from $\sim$1$\%$ to $\sim$3$\%$ over the observed range, consistent with prior predictions from the inhibited grain growth model based on Spitzer spectroscopy. We explore model refinements including photon effects and alternative size evolution prescriptions, and find that a modest amount of small grain photo-destruction remains possible, provided the grain size cutoff does not exceed $\sim55$ carbon atoms. The best-fit models predict 3.3 $μ$m/$Σ$PAH will rise to $\sim5.6-7.7\%$ at 10$\%$ Z$_{\odot}$. Surprisingly, even as $Σ$PAH drops significantly relative to the total infrared luminosity (TIR) as metallicity declines, 3.3 $μ$m/TIR alone rises, potentially indicating the mass fraction of the smallest PAH grains increases as the total dust content in galaxies drops. The current model cannot fully reproduce this trend even if the unusually strong effects of changing radiation field hardness on 3.3 $μ$m/TIR are included. This may be evidence that the smallest PAHs are uniquely robust against destruction and inhibited growth effects. These results highlight the pivotal role that short-wavelength PAH emission can play in studies of low-metallicity and high-redshift galaxies. △ Less

Submitted 23 February, 2026; v1 submitted 22 September, 2025; originally announced September 2025.

Comments: 19 pages, 6 figures, accepted for publication in ApJ

An Ultraviolet Study of CO Chemistry in the Magellanic Clouds

Authors: Kirill Tchernyshyov, Jessica K. Werk, Julia Roman-Duval

Abstract: How does molecular cloud chemistry change with metallicity? In this work, we study the relation between molecular hydrogen ($H_2$) and carbon monoxide (CO) at $1/2$ and $1/5$ solar metallicity using ultraviolet absorption spectroscopy obtained as part of the UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Hubble Space Telescope (HST) program. We determine CO column densities or u… ▽ More How does molecular cloud chemistry change with metallicity? In this work, we study the relation between molecular hydrogen ($H_2$) and carbon monoxide (CO) at $1/2$ and $1/5$ solar metallicity using ultraviolet absorption spectroscopy obtained as part of the UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Hubble Space Telescope (HST) program. We determine CO column densities or upper limits for a sample of 50 lines of sight through the Large and Small Magellanic Clouds (LMC and SMC). $^{12}$CO is detected along eight lines of sight and $^{13}$CO is detected along two. Combining our new CO column densities with $N_{H_2}$ measurements from the literature, we find that the evolution of $N_{CO}(N_{H_2})$ from the Milky Way to the LMC and SMC is a relatively shallow function of metallicity. Taking $N_{CO}>3\times10^{15}$ cm$^{-2}$ as a threshold value above which CO emission is likely to be detectable at the distance of the Magellanic Clouds, the $\log_{10} N_{H_2}$ at which a sightline has a 50% probability of having $N_{CO}$ above this threshold is 20.8 in the Milky Way, 20.9 in the LMC, and 21.1 in the SMC. This is an 0.3 dex change in threshold $\log_{10} N_{H_2}$ over an 0.7 dex change in metallicity. We compare our measurements with $N_{CO}(N_{H_2})$ relations from literature chemical models and find that the measured relations agree best with models in which the dynamical timescale is longer than the chemical timescale for $H_2$ but shorter than the chemical timescale for CO. △ Less

Submitted 14 August, 2025; originally announced August 2025.

Comments: Accepted for publication to AJ

Discovery of SiC and Iron Dust Around AGB Stars in the very Metal-Poor Sextans A Dwarf Galaxy with JWST: Implications for Dust Production at High Redshift

Authors: M. L. Boyer, G. C. Sloan, A. Nanni, E. Tarantino, I. McDonald, S. Goldman, J. A. D. L Blommaert, F. Dell'Agli, M. Di Criscienzo, D. A. Garcia-Hernandez, R. D. Gehrz, M. A. T. Groenewegen, A. Javadi, O. C. Jones, F. Kemper, M. Marengo, K. B. W. McQuinn, J. M. Oliveira, G. Pastorelli, J. Roman-Duval, R. Sahai, E. D. Skillman, S. Srinivasan, J. Th. van Loon, D. R. Weisz , et al. (1 additional authors not shown)

Abstract: Low-resolution infrared spectroscopy from JWST confirms the presence of SiC and likely metallic iron dust around asymptotic giant branch (AGB) stars in the Sextans A dwarf galaxy, which has a metallicity ~1%-7% Z_sun. While metal-poor carbon-rich AGB stars are known to produce copious amounts of amorphous carbon dust owing to the dredge up of newly synthesized carbon, this is the first time that S… ▽ More Low-resolution infrared spectroscopy from JWST confirms the presence of SiC and likely metallic iron dust around asymptotic giant branch (AGB) stars in the Sextans A dwarf galaxy, which has a metallicity ~1%-7% Z_sun. While metal-poor carbon-rich AGB stars are known to produce copious amounts of amorphous carbon dust owing to the dredge up of newly synthesized carbon, this is the first time that Si- and Fe-bearing dust has been detected at this extreme metallicity. Of the six AGB stars observed, one is an intermediate-mass (~1.2-4 M_sun) carbon star showing SiC dust, and another is an oxygen-rich M-type star with mass ~4-5 M_sun that is likely undergoing hot bottom burning. The infrared excess of the M-type star is strong, but featureless. We tested multiple dust species, and find that it is best fit with metallic iron dust. Assuming its dust-production rate stays constant over the final 2-3x10^4 yr of its evolution, this star will produce ~0.9-3.7 times the iron dust mass predicted by models, with the range depending on the adopted stellar mass. The implications for dust production in high-redshift galaxies are potentially significant, especially regarding the assumed dust species used in dust evolution models and the timescale of AGB dust formation. Stars on the upper end of the AGB mass range can begin producing dust as early as 30-50 Myr after they form, and they may therefore rival dust production by supernovae at high redshift. △ Less

Submitted 22 July, 2025; originally announced July 2025.

Comments: Accepted for publication in ApJ

Probing the variations of interstellar dust abundance and properties within and between galaxies with HWO UV spectroscopy in the Local Volume

Authors: Julia Roman-Duval, Yumi Choi, Mederic Boquien

Abstract: The cycle of metals between the gas and the dust phases in the neutral interstellar medium (ISM) is an integral part of the baryon cycle in galaxies. The resulting variations in the abundance and properties of interstellar dust have important implications for how accurately we can trace the chemical enrichment of the universe over cosmic time. Multi-object UV spectroscopy with HWO can provide the… ▽ More The cycle of metals between the gas and the dust phases in the neutral interstellar medium (ISM) is an integral part of the baryon cycle in galaxies. The resulting variations in the abundance and properties of interstellar dust have important implications for how accurately we can trace the chemical enrichment of the universe over cosmic time. Multi-object UV spectroscopy with HWO can provide the large samples of abundance and dust depletion measurements needed to understand how the abundance and properties of interstellar dust vary within and between galaxies, thereby observationally addressing important questions about chemical enrichment and galaxy evolution. Medium-resolution (R~50,000) spectroscopy in the full UV range (950-3150 A) toward massive stars in Local Volume galaxies (D < 10 Mpc) will enable gas- and dust-phase abundance measurements of key elements, such as Fe, Si, Mg, S, Zn. These measurements will provide an estimate of how the dust abundance varies with environment, in particular metallicity and gas density. However, measuring the carbon and oxygen contents of dust requires very high resolution (R > 100,000) and high signal-to-noise (S/N > 100) owing to the non-saturated UV transitions for those elements being extremely weak. Since carbon and oxygen in the neutral ISM contribute the largest metal mass reservoir for dust, it is critical that the HWO design include a grating similar to the HST STIS H gratings providing very high resolution, as well as FUV and NUV detectors capable of reaching very high S/N. △ Less

Submitted 9 September, 2025; v1 submitted 30 June, 2025; originally announced July 2025.

Comments: 12 pages; 6 figures; will be published in ASP conference proceedings of the HWO2025 conference

Investigating silicate, carbon, and water in the diffuse interstellar medium: the first shots from WISCI

Authors: S. T. Zeegers, Jonathan P. Marshall, Karl D. Gordon, Karl A. Misselt, G. P. P. L. Otten, Jeroen Bouwman, Jean Chiar, Marjorie Decleir, Thavisha Dharmawardena, F. Kemper, Aigen Li, Mayank Narang, Alexey Potapov, Manoj Puravankara, Peter Scicluna, Himanshu Tyagi, Eleonora Zari, ChuanYu Wei, Lex Kaper, Frank Backs, Stefan T. Bromley, Laurie Chu, Elisa Costantini, T. R. Geballe, Joel D. Green , et al. (11 additional authors not shown)

Abstract: The dusty interstellar medium (ISM) of the Milky Way is distributed in a complex, cloudy structure. It is fundamental to the radiation balance within the Milky Way, provides a reaction surface to form complex molecules, and is the feedstock for future generations of stars and planets. The life cycle of interstellar dust is not completely understood, and neither are its structure nor composition. T… ▽ More The dusty interstellar medium (ISM) of the Milky Way is distributed in a complex, cloudy structure. It is fundamental to the radiation balance within the Milky Way, provides a reaction surface to form complex molecules, and is the feedstock for future generations of stars and planets. The life cycle of interstellar dust is not completely understood, and neither are its structure nor composition. The abundance, composition, and structure of dust in the diffuse ISM can be determined by combining infrared, optical and ultraviolet spectroscopy. JWST enables measurement of the faint absorption of ISM dust grains against bright stars at kiloparsec distances across the infrared spectrum. Here we present an overview of the project `Webb Investigation of Silicates, Carbons, and Ices' (WISCI) along with interpretation of two targets, GSC 08152-02121 and CPD-59 5831. Observations of 12 WISCI target stars were taken by JWST, the Hubble Space Telescope, Himalayan Chandra Telescope, and the Very Large Telescope. We use these to characterize the targets' spectral types and calculate their line-of-sight extinction parameters, $A_{\rm V}$ and $R_{\rm V}$. We find absorption in the JWST spectra of GSC 08152-02121, and CPD-59 5831 associated with carbonaceous dust around 3.4 and 6.2 micron and amorphous silicates at 9.7 micron. In GSC 08152-02121 we also find indications of absorption by trapped water around 3 micron. This first look from WISCI demonstrates the line-of-sight variability within the sample, and the program's potential to identify and correlate features across ultraviolet to mid-infrared wavelengths. △ Less

Submitted 24 June, 2025; originally announced June 2025.

Comments: Accepted for publication in ApJ

The Karl G. Jansky Very Large Array Local Group L-band Survey (LGLBS)

Authors: Eric W. Koch, Adam K. Leroy, Erik W. Rosolowsky, Laura Chomiuk, Julianne J. Dalcanton, Nickolas M. Pingel, Sumit K. Sarbadhicary, Snežana Stanimirović, Fabian Walter, Haylee N. Archer, Alberto D. Bolatto, Michael P. Busch, Hongxing Chen, Ryan Chown, Harrisen Corbould, Serena A. Cronin, Jeremy Darling, Thomas Do, Jennifer Donovan Meyer, Cosima Eibensteiner, Deidre Hunter, Rémy Indebetouw, Preshanth Jagannathan, Amanda A. Kepley, Chang-Goo Kim , et al. (23 additional authors not shown)

Abstract: We present the Local Group L-Band Survey (LGLBS), a Karl G. Jansky Very Large Array (VLA) survey producing the highest quality 21-cm and 1-2 GHz radio continuum images to date for the six VLA-accessible, star-forming, Local Group galaxies. Leveraging the VLA's spectral multiplexing power, we simultaneously survey the 21-cm line at high 0.4 km/s velocity resolution, the 1-2 GHz polarized continuum,… ▽ More We present the Local Group L-Band Survey (LGLBS), a Karl G. Jansky Very Large Array (VLA) survey producing the highest quality 21-cm and 1-2 GHz radio continuum images to date for the six VLA-accessible, star-forming, Local Group galaxies. Leveraging the VLA's spectral multiplexing power, we simultaneously survey the 21-cm line at high 0.4 km/s velocity resolution, the 1-2 GHz polarized continuum, and four OH lines. For the massive spiral M31, the dwarf spiral M33, and the dwarf irregular galaxies NGC6822, IC10, IC1613, and the Wolf-Lundmark-Melotte Galaxy (WLM), we use all four VLA configurations and the Green Bank Telescope to reach angular resolutions of $< 5''$ ($10{-}20$~pc) for the 21-cm line with $<10^{20}$~cm$^{-2}$ column density sensitivity, and even sharper views ($< 2''$; $5{-}10$~pc) of the continuum. Targeting these nearby galaxies ($D\lesssim1$ Mpc) reveals a sharp, resolved view of the atomic gas, including 21-cm absorption, and continuum emission from supernova remnants and HII regions. These datasets can be used to test theories of the abundance and formation of cold clouds, the driving and dissipation of interstellar turbulence, and the impact of feedback from massive stars and supernovae. Here, we describe the survey design and execution, scientific motivation, data processing, and quality assurance. We provide a first look at and publicly release the wide-field 21-cm HI data products for M31, M33, and four dwarf irregular targets in the survey, which represent some of the highest physical resolution 21-cm observations of any external galaxies beyond the LMC and SMC. △ Less

Submitted 13 June, 2025; originally announced June 2025.

Comments: ApJS in press. LGLBS HI v1.0 data release is available here: https://www.canfar.net/storage/vault/list/LGLBS/RELEASES/LGLBS-HI-v1.0 (with permanent DOI to follow)

Measuring Interstellar Carbon Abundance via 158 um [CII] Absorption with SOFIA -- A Potential Detection, and Proof-of-Concept for Depletion Studies with Future Far-IR Facilities

Authors: Christopher J. R. Clark, Julia. C. Roman-Duval, Suzanne C. Madden Marc Mertens, Claire E. Murray, Jürgen Stutzki, Elizabeth Tarantino, and Kirill Tchernyshyov

Abstract: Carbon plays key roles in the InterStellar Medium (ISM) -- as a constituent of dust, as the carrier of the dominant far-infrared cooling line, and as a component of various important molecules. But despite this, there are very few measurements of the abundance and depletion of carbon in the diffuse ISM. As with other elements, these measurements are traditionally performed in the ultraviolet. But… ▽ More Carbon plays key roles in the InterStellar Medium (ISM) -- as a constituent of dust, as the carrier of the dominant far-infrared cooling line, and as a component of various important molecules. But despite this, there are very few measurements of the abundance and depletion of carbon in the diffuse ISM. As with other elements, these measurements are traditionally performed in the ultraviolet. But for carbon, such measurements are extremely difficult, and less than 20 have been reported in the literature to date. Here, we present a novel method of measuring the abundance and depletion of carbon in the diffuse ISM: by observing absorption of the 158 $μ$m [CII] line in the far-infrared. We present a catalog of 432 candidate sightlines that use bright nearby galaxies as background sources, and predict the [CII] absorption expected towards each. We conducted a pilot study using SOFIA, targeting sightlines towards the galaxies IC342 and Circinus. We report a potential detection of Galactic [CII] absorption along the IC342 sightline, although it requires disentangling [CII] emission from IC342 itself. The Circinus sightline had an insufficiently stable instrumental baseline to allow a detection. This SOFIA study informs the prospects for [CII] absorption measurements with future facilities. To that end, we explore the potential for four proposed future FIR telescopes -- PRIMA, FIRSST, SALTUS, and Origins -- to detect [CII] absorption. We find that all four facilities would be able to detect [CII] absorption along a significant number of sightlines. △ Less

Submitted 5 June, 2025; v1 submitted 5 May, 2025; originally announced May 2025.

Comments: Accepted for publication in the Astronomical Journal; machine-readable versions of data tables accessible with online open-access version of published article at the journal: https://iopscience.iop.org/article/10.3847/1538-3881/add40f

Journal ref: AJ 170 8 (2025)

The UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Large Director's Discretionary Program with Hubble. I. Goals, Design, and Initial Results

Authors: Julia Roman-Duval, William J. Fischer, Alexander W. Fullerton, Jo Taylor, Rachel Plesha, Charles Proffitt, TalaWanda Monroe, Travis C. Fischer, Alessandra Aloisi, Jean-Claude Bouret, Christopher Britt, Nuria Calvet, Joleen K. Carlberg, Paul A. Crowther, Gisella De Rosa, William V. Dixon, Catherine C. Espaillat, Christopher J. Evans, Andrew J. Fox, Kevin France, Miriam Garcia, Sott W. Fleming, Elaine M. Frazer, Ana I. Gómez De Castro, Gregory J. Herczeg , et al. (22 additional authors not shown)

Abstract: Specifically selected to leverage the unique ultraviolet capabilities of the Hubble Space Telescope, the Hubble Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) is a Director's Discretionary program of approximately 1000 orbits - the largest ever executed - that produced a UV spectroscopic library of O and B stars in nearby low metallicity galaxies and accreting low mass… ▽ More Specifically selected to leverage the unique ultraviolet capabilities of the Hubble Space Telescope, the Hubble Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) is a Director's Discretionary program of approximately 1000 orbits - the largest ever executed - that produced a UV spectroscopic library of O and B stars in nearby low metallicity galaxies and accreting low mass stars in the Milky Way. Observations from ULLYSES combined with archival spectra uniformly sample the fundamental astrophysical parameter space for each mass regime, including spectral type, luminosity class, and metallicity for massive stars, and the mass, age, and disk accretion rate for low-mass stars. The ULLYSES spectral library of massive stars will be critical to characterize how massive stars evolve at different metallicities; to advance our understanding of the production of ionizing photons, and thus of galaxy evolution and the re-ionization of the Universe; and to provide the templates necessary for the synthesis of integrated stellar populations. The massive star spectra are also transforming our understanding of the interstellar and circumgalactic media of low metallicity galaxies. On the low-mass end, UV spectra of T Tauri stars contain a plethora of diagnostics of accretion, winds, and the warm disk surface. These diagnostics are crucial for evaluating disk evolution and provide important input to assess atmospheric escape of planets and to interpret powerful probes of disk chemistry, as observed with ALMA and JWST. In this paper we motivate the design of the program, describe the observing strategy and target selection, and present initial results. △ Less

Submitted 7 April, 2025; originally announced April 2025.

Comments: 58 pages, 20 figures

An Investigation of Disk Thickness in M51 from H-alpha, Pa-alpha, and Mid-Infrared Power Spectra

Authors: Bruce G. Elmegreen, Daniela Calzetti, Angela Adamo, Karin Sandstrom, Daniel Dale, Varun Bajaj, Martha L. Boyer, Ana Duarte-Cabral, Ryan Chown, Matteo Correnti, Julianne J. Dalcanton, Bruce T. Draine, Brandt Gaches, John S. Gallagher III, Kathryn Grasha, Benjamin Gregg, Leslie K. Hunt, Kelsey E. Johnson, Robert Kennicutt, Jr., Ralf S. Klessen, Adam K. Leroy, Sean Linden, Anna F. McLeod, Matteo Messa, Goran Ostlin , et al. (5 additional authors not shown)

Abstract: Power spectra (PS) of high-resolution images of M51 (NGC 5194) taken with the Hubble Space Telescope and the James Webb Space Telescope have been examined for evidence of disk thickness in the form of a change in slope between large scales, which map two-dimensional correlated structures, and small scales, which map three-dimensional correlated structures. Such a slope change is observed here in H… ▽ More Power spectra (PS) of high-resolution images of M51 (NGC 5194) taken with the Hubble Space Telescope and the James Webb Space Telescope have been examined for evidence of disk thickness in the form of a change in slope between large scales, which map two-dimensional correlated structures, and small scales, which map three-dimensional correlated structures. Such a slope change is observed here in H-alpha, and possibly Pa-alpha, using average PS of azimuthal intensity scans that avoid bright peaks. The physical scale of the slope change occurs at ~120 pc and ~170 pc for these two transitions, respectively. A radial dependence in the shape of the H-alpha PS also suggests that the length scale drops from ~180 pc at 5 kpc, to ~90 pc at 2 kpc, to ~25 pc in the central ~kpc. We interpret these lengths as comparable to the thicknesses of the star-forming disk traced by HII regions. The corresponding emission measure is ~100 times larger than what is expected from the diffuse ionized gas. PS of JWST Mid-IR Instrument (MIRI) images in 8 passbands have more gradual changes in slope, making it difficult to determine a specific value of the thickness for this emission. △ Less

Submitted 9 April, 2025; v1 submitted 7 April, 2025; originally announced April 2025.

Comments: 20 pages, 7 figures, accepted by Astrophysical Journal April 5, 2025

Detection of Deuterated Hydrocarbon Nanoparticles in the Whirlpool Galaxy, M51

Authors: B. T. Draine, Karin Sandstrom, Daniel A. Dale, J. -D. T. Smith, Ryan Chown, Grant P. Donnelly, Sara E. Duval, Cory M. Whitcomb, Angela Adamo, L. Armus, Danielle A. Berg, Torsten Böker, Alberto D. Bolatto, Martha L. Boyer, Daniela Calzetti, B. G. Elmegreen, Brandt A. L. Gaches, Karl D. Gordon, L. K. Hunt, R. C. Kennicutt, Ralf S. Klessen, Thomas S. -Y. Lai, Adam K. Leroy, Sean T. Linden, Alex Pedrini , et al. (7 additional authors not shown)

Abstract: Deuteration of hydrocarbon material, including polycyclic aromatic hydrocarbons (PAHs), has been proposed to account for the low gas-phase abundances of D in the interstellar medium. JWST spectra of four star-forming regions in M51 show an emission feature, with central wavelength $\sim$4.647$μ$m and FWHM 0.0265$μ$m, corresponding to the C-D stretching mode in aliphatic hydrocarbons. The emitting… ▽ More Deuteration of hydrocarbon material, including polycyclic aromatic hydrocarbons (PAHs), has been proposed to account for the low gas-phase abundances of D in the interstellar medium. JWST spectra of four star-forming regions in M51 show an emission feature, with central wavelength $\sim$4.647$μ$m and FWHM 0.0265$μ$m, corresponding to the C-D stretching mode in aliphatic hydrocarbons. The emitting aliphatic material is estimated to have (D/H)$_{\rm aliph}\approx 0.17\pm0.02$ -- a factor $\sim$$10^4$ enrichment relative to the overall interstellar medium (ISM). On $\sim$$50\,$pc scales, deuteration levels toward four H$\,$II regions in M51 are 2-3 times higher than in the Orion Bar photodissociation region (PDR), with implications for the processes responsible for the formation and evolution of hydrocarbon nanoparticles, including PAHs. The deuteration of the aliphatic material is found to anticorrelate with helium ionization in the associated H$\,$II, suggesting that harsh FUV radiation may act to lower the deuteration of aliphatics in PDRs near massive stars. No evidence is found for deuteration of aromatic material, with (D/H)$_{\rm arom} \lesssim 0.016$: deuteration of the aliphatic material exceeds that of the aromatic material by at least a factor 10. The observed levels of deuteration may account for the depletion of D observed in the Galactic interstellar medium. If so, the $4.65μ$m feature may be detectable in absorption. △ Less

Submitted 13 May, 2025; v1 submitted 3 April, 2025; originally announced April 2025.

Comments: Revised. corresponding to published paper

Journal ref: Astrophys. J. Letters, 948:L42 (2025)

The Gaseous Blowout of the 30 Doradus Starburst Region in the LMC

Authors: Suraj Poudel, April Horton, Jo Vazquez, Kathleen A. Barger, Frances H. Cashman, Andrew J. Fox, Nicolas Lehner, Scott Lucchini, Dhanesh Krishnarao, N. M. McClure-Griffiths, Elena D'Onghia, Jason Tumlinson, Ananya Goon Tuli, Lauren Sdun, Stone Gebhart, Katherine Anthony, Bryce Cole, Jacco Th. van Loon, Julia Roman-Duval, Yik Ki Ma, Callum Lynn, Min-Young Lee, Denis Leahy

Abstract: Widespread galactic winds emanate from the Large Magellanic Cloud (LMC), with the 30 Doradus starburst region generating the fastest and most concentrated gas flows. We report on the gas distribution, kinematics, and ionization conditions of the near-side outflow along 8 down-the-barrel sightlines using UV absorption-line observations from the HST's ULLYSES program for this region along with H I 2… ▽ More Widespread galactic winds emanate from the Large Magellanic Cloud (LMC), with the 30 Doradus starburst region generating the fastest and most concentrated gas flows. We report on the gas distribution, kinematics, and ionization conditions of the near-side outflow along 8 down-the-barrel sightlines using UV absorption-line observations from the HST's ULLYSES program for this region along with H I 21-cm observations from the GASS and GASKAP surveys. We find that within 1.7 degrees from the center of 30 Doradus, the wind reaches maximum speeds of $100-150\,\text{km}\,\text{s}^{-1}$ from the LMC's disk. The total integrated column densities of low-ions (O I, Si II, and Fe II) in the blueshifted wind, up to $v_{\rm LSR}=150\,\text{km}\,\text{s}^{-1}$, are highest near the center and decline radially outward. We estimate an outflow mass of $M_{\rm outflow,\,Si II}\approx(5.7-8.6)\,\times 10^{5} M_{\odot}$, outflow rate of $\dot{M}_{\rm outflow}\gtrsim0.02 M_{\odot}\,\text{yr}^{-1}$, and mass loading factor of $η\gtrsim0.10$ within 0.52 degrees from the center of 30 Doradus. The observed ion ratios$-$together with photoionization modeling$-$reveal that this wind is roughly $40-97\%$ photoionized. The metallicities and dust depletion patterns of the high-velocity absorbers at $v_{\rm LSR}\approx+120\,\text{km}\,\text{s}^{-1}$ can be explained by either a foreground Milky Way (MW) halo cloud or an outflow from the LMC. For the high-ions, Si IV and C IV are broader and kinematically offset from the low-ions, suggesting turbulent mixing layers (TMLs) existing in the wind. Finally, our hydrodynamical simulations of the Magellanic Clouds (MCs) and MW system suggest that the Magellanic Corona can protect the LMC winds from the ram-pressure forces exerted by the MW's halo. △ Less

Submitted 7 March, 2025; originally announced March 2025.

Comments: 40 pages, 23 figures, accepted for publication in ApJ

A Catalog of Stellar and Dust Properties for 500,000 Stars in the Southwest Bar of the Small Magellanic Cloud

Authors: Petia Yanchulova Merica-Jones, Karl Gordon, Karin Sandstrom, Claire E. Murray, L. Clifton Johnson, Julianne J. Dalcanton, Julia Roman-Duval, Jeremy Chastenet, Benjamin F. Williams, Daniel R. Weisz, Andrew E. Dolphin

Abstract: We present a catalog of individual stellar and dust extinction properties along close to 500,000 sight lines in the southwest bar of the Small Magellanic Cloud (SMC). The catalog is based on multiband Hubble Space Telescope photometric data spanning near-ultraviolet to near-infrared wavelengths from the Small Magellanic Cloud Investigation of Dust and Gas Evolution survey (SMIDGE) covering a 100 x… ▽ More We present a catalog of individual stellar and dust extinction properties along close to 500,000 sight lines in the southwest bar of the Small Magellanic Cloud (SMC). The catalog is based on multiband Hubble Space Telescope photometric data spanning near-ultraviolet to near-infrared wavelengths from the Small Magellanic Cloud Investigation of Dust and Gas Evolution survey (SMIDGE) covering a 100 x 200 pc area. We use the probabilistic technique of the Bayesian Extinction And Stellar Tool (BEAST) to model the spectral energy distributions of individual stars in SMIDGE and include the effects of observational uncertainties in the data. We compare BEAST-derived dust extinction properties with tracers of the interstellar medium, such as the emission from the 12CO (2-1) transition (I(CO)), the dust mass surface density (Σdust) from far-IR emission, the H I column density (N(HI)) from the 21cm transition, and the mass fraction of polycyclic aromatic hydrocarbons (PAHs; qPAH, derived from IR emission). We find that the dust extinction (A(V )) in the SMIDGE field is strongly correlated with Σdust and I(CO), and less so with N(HI) and qPAH, and suggest potential explanations. Our extinction measurements are also sensitive to the presence of the 2175 Å bump in the extinction curve toward UV bright stars. While most do not show evidence for the bump, we identify ~200 lines of sight that are 2175 Å bump candidates. Furthermore, we find distinct structures in the dust extinction-distance distributions that provide insights into the 3D geometry of the SMC. △ Less

Submitted 10 January, 2025; originally announced January 2025.

Comments: 20 pages, 13 figures, 4 tables

Journal ref: ApJ 978 144 (2025)

A first taste of MEAD (Measuring Extinction and Abundances of Dust) -- I. Diffuse Milky Way interstellar dust extinction features in JWST infrared spectra

Authors: Marjorie Decleir, Karl D. Gordon, Karl A. Misselt, Burcu Günay, Julia Roman-Duval, Sascha T. Zeegers

Abstract: We present the initial results of MEAD (Measuring Extinction and Abundances of Dust), with a focus on the dust extinction features observed in our JWST near- and mid-infrared spectra of nine diffuse Milky Way sightlines ($1.2 \leq A(V) \leq 2.5$). For the first time, we find strong correlations between the 10 $μ$m silicate feature strength and the column densities of Mg, Fe and O in dust. This is… ▽ More We present the initial results of MEAD (Measuring Extinction and Abundances of Dust), with a focus on the dust extinction features observed in our JWST near- and mid-infrared spectra of nine diffuse Milky Way sightlines ($1.2 \leq A(V) \leq 2.5$). For the first time, we find strong correlations between the 10 $μ$m silicate feature strength and the column densities of Mg, Fe and O in dust. This is consistent with the well-established theory that Mg- and Fe-rich silicates are responsible for this feature. We obtained an average stoichiometry of the silicate grains in our sample of Mg:Fe:O = 1.1:1:11.2, constraining the grain composition. We find variations in the feature properties, indicating that different sightlines contain different types of silicates. In the average spectrum of our sample, we tentatively detect features around 3.4 and 6.2 $μ$m, which are likely caused by aliphatic and aromatic/olefinic hydrocarbons, respectively. If real, to our knowledge, this is the first detection of hydrocarbons in purely diffuse sightlines with $A(V) \leq 2.5$, confirming the presence of these grains in diffuse environments. We detected a 3 $μ$m feature toward HD073882, and tentatively in the sample average, likely caused by water ice (or solid-state water trapped on silicate grains). If confirmed, to our knowledge, this is the first detection of ice in purely diffuse sightlines with $A(V) \leq 2.5$, supporting previous findings that these molecules can exist in the diffuse ISM. △ Less

Submitted 18 December, 2024; originally announced December 2024.

Comments: 20 pages, 9 figures, 3 tables, accepted for publication in AJ, for the data used in this work, see https://doi.org/10.5281/zenodo.14286122 , for the code developed for this work, see https://doi.org/10.5281/zenodo.14291651

Scylla IV: Intrinsic Stellar Properties and Line-of-Sight Dust Extinction Measurements Towards 1.5 Million Stars in the SMC and LMC

Authors: Christina W. Lindberg, Claire E. Murray, Petia Yanchulova Merica-Jones, Caroline Bot, Clare Burhenne, Yumi Choi, Christopher J. R. Clark, Roger E. Cohen, Karoline M. Gilbert, Steven R. Goldman, Karl D. Gordon, Alec S. Hirschauer, Kristen B. W. McQuinn, Julia C. Roman-Duval, Karin M. Sandstrom, Elizabeth Tarantino, Benjamin F. Williams

Abstract: By analyzing the spectral energy distributions (SEDs) of resolved stars in nearby galaxies, we can constrain their stellar properties and line-of-sight dust extinction. From the Scylla survey, we obtain ultraviolet to near-infrared photometry from Wide Field Camera 3 onboard the {\it Hubble Space Telescope} for more than 1.5 million stars in the SMC and LMC. We use the Bayesian Extinction and Stel… ▽ More By analyzing the spectral energy distributions (SEDs) of resolved stars in nearby galaxies, we can constrain their stellar properties and line-of-sight dust extinction. From the Scylla survey, we obtain ultraviolet to near-infrared photometry from Wide Field Camera 3 onboard the {\it Hubble Space Telescope} for more than 1.5 million stars in the SMC and LMC. We use the Bayesian Extinction and Stellar Tool (BEAST) to analyze the multi-band SEDs of these sources and characterize their initial masses, ages, metallicities, distances, and line-of-sight extinction properties (e.g.~$A_V$, $R_V$). We apply quality cuts and perform validation simulations to construct a catalog of over 550,000 stars with high-reliability SED fits, which we use to analyze the stellar content and extinction properties of the SMC and LMC. We detect stars with masses as low as 0.6 $M_{\odot}$. BEAST stellar age distributions show a jump in observed stars around 6 Gyrs ago, which agrees with star-formation histories. Extinctions ($A_V$) in both galaxies follow a log-normal distribution. We compare $A_V$ with ancillary gas and dust tracers like $HI$, $H_α$, and far infrared (FIR) dust emission and find positive correlations on a field-by-field basis. We convert observed $A_V$ to predicted dust surface densities using the Draine et. al. (2014) model and find $A_V$-based dust surface densities are a factor of $\sim$2.5 lower than observed FIR-based dust surface densities, a correction factor similar to other studies. △ Less

Submitted 11 February, 2025; v1 submitted 25 October, 2024; originally announced October 2024.

Comments: Accepted to ApJ, updated to 38 pages, 20 figures

Scylla I: A pure-parallel, multi-wavelength imaging survey of the ULLYSES fields in the LMC and SMC

Authors: Claire E. Murray, Christina W. Lindberg, Petia Yanchulova Merica-Jones, Benjamin F. Williams, Roger E. Cohen, Karl D. Gordon, Kristen B. W. McQuinn, Yumi Choi, Clare Burhenne, Karin M. Sandstrom, Caroline Bot, L. Clifton Johnson, Steven R. Goldman, Christopher J. R. Clark, Julia C. Roman-Duval, Karoline M. Gilbert, J. E. G. Peek, Alec S. Hirschauer, Martha L. Boyer, Andrew E. Dolphin

Abstract: Scylla is a deep Hubble Space Telescope survey of the stellar populations, interstellar medium and star formation in the LMC and SMC. As a pure-parallel complement to the Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) survey, Scylla obtained 342 orbits of ultraviolet (UV) through near-infrared (IR) imaging of the LMC and SMC with Wide Field Camera 3. In this paper, we d… ▽ More Scylla is a deep Hubble Space Telescope survey of the stellar populations, interstellar medium and star formation in the LMC and SMC. As a pure-parallel complement to the Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) survey, Scylla obtained 342 orbits of ultraviolet (UV) through near-infrared (IR) imaging of the LMC and SMC with Wide Field Camera 3. In this paper, we describe the science objectives, observing strategy, data reduction procedure, and initial results from our photometric analysis of 96 observed fields. Although our observations were constrained by ULYSSES primary exposures, we imaged all fields in at least two filters (F475W and F814W), and 64% of fields in at least three and as many as seven WFC3 filters spanning the UV to IR. Overall, we reach average 50% completeness of $m_{\rm F225W}=26.0$, $m_{\rm F275W}=26.2$, $m_{\rm F336W}=26.9$, $m_{\rm F475W}=27.8$, $m_{\rm F814W}=25.5$, $m_{\rm F110W}=24.7$, and $m_{\rm F160W}=24.0$ Vega magnitudes in our photometric catalogs, which is faintward of the ancient main sequence turnoff in all filters. The primary science goals of Scylla include characterizing the structure and properties of dust in the MCs, as well as their spatially-resolved star formation and chemical enrichment histories. Our images and photometric catalogs, which represent the widest-area coverage of MCs with HST photometry to date, are available as a high-level science product at the Barbara A. Mikulski Archive for Space Telescopes. △ Less

Submitted 15 October, 2024; originally announced October 2024.

Comments: 25 pages, 16 figures, 8 tables. Accepted for publication in ApJS

DUNE: Dust depletion UNified method across cosmic time and Environments

Authors: Christina Konstantopoulou, Annalisa De Cia, Jens-Kristian Krogager, Cédric Ledoux, Julia Roman-Duval, Edward B. Jenkins, Tanita Ramburuth-Hurt, Anna Velichko

Abstract: We present a novel method to characterize dust depletion, namely, the depletion of metals into dust grains. We used observed correlations among relative abundances combining a total of 17 metals in diverse galactic environments, including the Milky Way (MW), Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and damped Lyman-$α$ absorbers (DLAs) towards quasars and gamma-ray bursts (GRBs)… ▽ More We present a novel method to characterize dust depletion, namely, the depletion of metals into dust grains. We used observed correlations among relative abundances combining a total of 17 metals in diverse galactic environments, including the Milky Way (MW), Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and damped Lyman-$α$ absorbers (DLAs) towards quasars and gamma-ray bursts (GRBs). We only considered the relative abundances of metals that qualify as tracers of dust and we used all available dust tracers. We find linear correlations among all studied dust tracers in a multidimensional space, where each dimension corresponds to an individual dust tracer. The fit to the linear correlations among the dust tracers describes the tendencies of different elements when depleting into dust grains. We determined the overall strength of dust depletion, $Δ$, along individual lines of sight, based on the correlations among different dust tracers. We avoided any preference for specific dust tracers or any other assumptions by including all available dust tracers in this multidimensional space. We also determined the dust depletion of Kr, C, O, Cl, P, Zn, Ge, Mg, Cu, Si, Fe, Ni, and Ti. Finally, we offer simple guidelines for the application of the method to the study of the observed patterns of abundances and relative abundances. This has allowed for a straightforward determination of the overall strength of depletion and the dust depletion of individual elements. We also obtained an estimate for the gas-phase metallicity and identified any additional deviations due to the nucleosynthesis of specific stellar populations. Thus, we have established a unified methodology for characterizing dust depletion across cosmic time and diverse galactic environments, offering a valuable new approach to the study of dust depletion in studies of the chemical evolution of galaxies. △ Less

Submitted 8 October, 2024; originally announced October 2024.

Comments: Accepted for publication in A&A

A radial decrease in kinetic temperature measured with formaldehyde in 30 Doradus

Authors: Remy Indebetouw, Tony Wong, Suzanne Madden, Marta Sewilo, Julia Roman-Duval, Melanie Chevance, Monica Rubio

Abstract: Feedback from star formation is a critical component of the evolution of galaxies and their interstellar medium. At parsec scales internal to molecular clouds, however, the observed signatures of that feedback on the physical properties of CO-emitting gas have often been weak or inconclusive. We present sub-parsec observations of H2CO in the 30 Doradus region, which contains the massive star clust… ▽ More Feedback from star formation is a critical component of the evolution of galaxies and their interstellar medium. At parsec scales internal to molecular clouds, however, the observed signatures of that feedback on the physical properties of CO-emitting gas have often been weak or inconclusive. We present sub-parsec observations of H2CO in the 30 Doradus region, which contains the massive star cluster R136 that is clearly exerting feedback on its neighboring gas. H2CO provides a direct measure of gas kinetic temperature, and we find a trend of decreasing temperature with projected distance from R136 that may be indicative of gas heating by the stars. While it has been suggested that mechanical heating affects H2CO-measured temperature, we do not observe any correlation between TK and line width. The lack of an enhancement in mechanical feedback close to R136 is consistent with the absence of a radial trend in gravitational boundedness seen the ALMA CO observations. Estimates of cosmic ray flux in the region are quite uncertain but can plausibly explain the observed temperatures, if R136 itself is the dominant local source of energetic protons. The observations presented here are also consistent with the H2CO-emitting gas near R136 being dominated by direct radiation from R136 and photoelectric heating in the photodissociation regions. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: accepted to ApJ

PDRs4All VIII: Mid-IR emission line inventory of the Orion Bar

Authors: Dries Van De Putte, Raphael Meshaka, Boris Trahin, Emilie Habart, Els Peeters, Olivier Berné, Felipe Alarcón, Amélie Canin, Ryan Chown, Ilane Schroetter, Ameek Sidhu, Christiaan Boersma, Emeric Bron, Emmanuel Dartois, Javier R. Goicoechea, Karl D. Gordon, Takashi Onaka, Alexander G. G. M. Tielens, Laurent Verstraete, Mark G. Wolfire, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Jan Cami, Sara Cuadrado , et al. (113 additional authors not shown)

Abstract: Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observat… ▽ More Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Journal ref: A&A 687, A86 (2024)

A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk

Authors: Olivier Berné, Emilie Habart, Els Peeters, Ilane Schroetter, Amélie Canin, Ameek Sidhu, Ryan Chown, Emeric Bron, Thomas J. Haworth, Pamela Klaassen, Boris Trahin, Dries Van De Putte, Felipe Alarcón, Marion Zannese, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Jan Cami, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea , et al. (121 additional authors not shown)

Abstract: Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of… ▽ More Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Journal ref: Science, 383, 6686, 2024

METAL-Z: Measuring dust depletion in low metalicity dwarf galaxies

Authors: Aleksandra Hamanowicz, Kirill Tchernyshyov, Julia Roman-Duval, Edward B. Jenkins, Marc Rafelski, Karl D. Gordon, Yong Zheng, Miriam Garcia, Jessica Werk

Abstract: The cycling of metals between interstellar gas and dust is a critical aspect of the baryon cycle of galaxies, yet our understanding of this process is limited. This study focuses on understanding dust depletion effects in the low metallicity regime (< 20% Zo) typical of cosmic noon. Using medium-resolution UV spectroscopy from the COS onboard the Hubble Space Telescope, gas-phase abundances and de… ▽ More The cycling of metals between interstellar gas and dust is a critical aspect of the baryon cycle of galaxies, yet our understanding of this process is limited. This study focuses on understanding dust depletion effects in the low metallicity regime (< 20% Zo) typical of cosmic noon. Using medium-resolution UV spectroscopy from the COS onboard the Hubble Space Telescope, gas-phase abundances and depletions of iron and sulfur were derived toward 18 sightlines in local dwarf galaxies IC 1613 and Sextans A. The results show that the depletion of Fe and S is consistent with that found in the Milky Way, LMC and SMC. The depletion level of Fe increases with gas column density, indicating dust growth in the interstellar medium (ISM). The level of Fe depletion decreases with decreasing metallicity, resulting in the fraction of iron in gas ranging from 3% in the MW to 9% in IC 1613 and ~19% in Sextans A. The dust-to-gas and dust-to-metal ratios (D/G, D/M) for these dwarf galaxies were estimated based on the MW relations between the depletion of Fe and other elements. The study finds that D/G decreases only slightly sub-linearly with metallicity, with D/M decreasing from 0.41 +/- 0.05 in the MW to 0.11 +/- 0.11 at 0.10 Zo (at log N(H) = 21 cm-2). The trend of D/G vs. metallicity using depletion in local systems is similar to that inferred in Damped Ly-alpha systems from abundance ratios but lies higher than the trend inferred from FIR measurements in nearby galaxies. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Braving the Storm: Quantifying Disk-wide Ionized Outflows in the Large Magellanic Cloud with ULLYSES

Authors: Yong Zheng, Kirill Tchernyshyov, Knut Olsen, Yumi Choi, Chad Bustard, Julia Roman-Duval, Robert Zhu, Enrico M. Di Teodoro, Jessica Werk, Mary Putman, Anna F. McLeod, Yakov Faerman, Raymond C. Simons, Joshua Peek

Abstract: The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regre… ▽ More The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($Σ_{\rm SFR}$), and the outflows with higher $Σ_{\rm SFR}$ tend to be more ionized. Considering outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $η\gtrsim 0.15$. We compare the LMC's outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto Σ_{\rm SFR}^{0.23}$ over a wide range of star-forming conditions ($Σ_{\rm SFR}\sim10^{-4.5}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way. △ Less

Submitted 4 August, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ. Key findings can be found in Figures 9-12. Main updates from the previous version include new estimates on star formation rate surface densities. Normalized SII, SiIV, and CIV line spectra derived for this work are published as a High Level Science Product called LMC-FLOWS (doi: 10.17909/hz0m-np43), available on website: https://archive.stsci.edu/hlsp/lmc-flows

α-element enhancements in the Magellanic Interstellar Medium: evidence for recent star formation

Authors: Annalisa De Cia, Julia Roman-Duval, Christina Konstantopoulou, Pasquier Noterdaeme, Tanita Ramburuth-Hurt, Anna Velichko, Andrew J. Fox, Cédric Ledoux, Patrick Petitjean, Iris Jermann, Jens-Kristian Krogager

Abstract: Important questions on the chemical composition of the neutral ISM in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) are still open. It is usually assumed that their metallicity is uniform and equal to that measured in hot stars and HII regions, but direct measurements on the neutral ISM have not been performed until now. Deriving the metallicity from the observed metal abundanc… ▽ More Important questions on the chemical composition of the neutral ISM in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) are still open. It is usually assumed that their metallicity is uniform and equal to that measured in hot stars and HII regions, but direct measurements on the neutral ISM have not been performed until now. Deriving the metallicity from the observed metal abundances is not straightforward because they also depend on the depletion of metals into dust as well as nucleosynthesis effects such as $α$-element enhancement. We collect literature column densities of TiII, NiII, CrII, FeII, MnII, SiII, CuII, MgII, SII, PII, ZnII, and OI in the neutral ISM towards 32 and 22 hot stars in the LMC and SMC. We measure the metallicity, dust depletion, and $α$-element enhancements in the neutral ISM in the LMC and SMC. We find $α$-element enhancements in the neutral ISM in most systems, on average 0.26 dex (0.35 dex) for the LMC (SMC), and Mn under-abundance in the SMC (on average $-0.35$ dex). These are higher than for stars at similar metallicities. The observed $α$-element enhancements and Mn under-abundance are likely due to bursts of star formation, more recently than ~1 Gyr ago, that enriched the ISM from core-collapse supernovae. We find total neutral ISM metallicities that are mostly consistent with hot stars metallicity, on average [M/H]$_{\rm tot} = -0.33$ ($-0.83$), in the LMC (SMC). In six systems, however, we find significantly lower metallicities, two out of 32 in the LMC (with ~16% solar) and four out of 22 in the SMC (3 and 10% solar), two of which are in the outskirts of the SMC near the Magellanic Bridge, a region known for having a lower metallicity. With the exception of lines of sight towards the Magellanic Bridge, the neutral gas in the LMC and SMC appears fairly well mixed in terms of metallicity. [abridged] △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication on A&A. 25 pages, 16 figures, 4 tables

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

Where do stars explode in the ISM? -- The distribution of dense gas around evolved massive stars in M33

Authors: Sumit K. Sarbadhicary, Jordan Wagner, Eric W. Koch, Ness Mayker Chen, Adam K. Leroy, Natalia Lahén, Erik Rosolowsky, Kathryn F. Neugent, Chang-Goo Kim, Laura Chomiuk, Julianne J. Dalcanton, Laura A. Lopez, Nickolas M. Pingel, Remy Indebetouw, Thomas G. Williams, Elizabeth Tarantino, Jennifer Donovan Meyer, Evan D. Skillman, Adam Smercina, Amanda A. Kepley, Eric J. Murphy, Jay Strader, Tony Wong, Snežana Stanimirović, Vicente Villanueva , et al. (5 additional authors not shown)

Abstract: The effect of supernovae (SNe) on star-formation in the interstellar medium (ISM) depends sensitively on where SNe explode with respect to ISM clouds. Observationally, SN ISM environments characterized by spatially-resolved gas maps can empirically guide the placement of SNe in subgrid models, but unfortunately such measurements remain scarce, as SNe are rare and often distant. Here we demonstrate… ▽ More The effect of supernovae (SNe) on star-formation in the interstellar medium (ISM) depends sensitively on where SNe explode with respect to ISM clouds. Observationally, SN ISM environments characterized by spatially-resolved gas maps can empirically guide the placement of SNe in subgrid models, but unfortunately such measurements remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach -- mapping the ISM around evolved massive stars that are soon to explode. These provide a substantially larger sample of `explosion sites' (than just historical SNe) in nearby galaxies that have high-resolution atomic and molecular ISM maps from Jansky VLA and ALMA. We demonstrate this technique in the well-resolved Local Group spiral M33 by analyzing the 50 pc-scale projected ISM densities around red supergiants (RSGs, 8-30 M$_{\odot}$ stars) Wolf-Rayet stars (WRs, $>$30 M$_{\odot}$ stars), and supernova remnants (SNRs). We find a \emph{mass-dependent} correlation between stars and gas clouds, with atleast 45\% of WRs and upto 77\% of RSGs having no detectable H$_2$ at their pixel locations. In the sample with H$_2$ detections, we find that more massive younger progenitors are coincident with denser gas. We show that the density distributions for stars $>$15 M$_{\odot}$ are statistically distinct from random alignment of stars and gas in M33. Our work provides the first observationally-derived estimate of the fraction of the SN-producing stellar population correlated with ISM density peaks. We demonstrate how this can be compared with galaxy simulations, and advocate similar comparisons to the community for constraining sub-grid models. △ Less

Submitted 26 November, 2025; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 34 pages, 14 figures. Accepted to ApJ. The density distributions are available as data behind figures in the journal publication. Please feel free to contact us in the meantime if you would like to use them

PDRs4All III: JWST's NIR spectroscopic view of the Orion Bar

Authors: Els Peeters, Emilie Habart, Olivier Berne, Ameek Sidhu, Ryan Chown, Dries Van De Putte, Boris Trahin, Ilane Schroetter, Amelie Canin, Felipe Alarcon, Bethany Schefter, Baria Khan, Sofia Pasquini, Alexander G. G. M. Tielens, Mark G. Wolfire, Emmanuel Dartois, Javier R. Goicoechea, Alexandros Maragkoudakis, Takashi Onaka, Marc W. Pound, Silvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma , et al. (113 additional authors not shown)

Abstract: (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion… ▽ More (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 52 pages, 30 figures, submitted to A&A

Journal ref: A&A 685, A74 (2024)

PDRs4All IV. An embarrassment of riches: Aromatic infrared bands in the Orion Bar

Authors: Ryan Chown, Ameek Sidhu, Els Peeters, Alexander G. G. M. Tielens, Jan Cami, Olivier Berné, Emilie Habart, Felipe Alarcón, Amélie Canin, Ilane Schroetter, Boris Trahin, Dries Van De Putte, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem El-Yajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (114 additional authors not shown)

Abstract: (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory o… ▽ More (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $μ$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $μ$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'. △ Less

Submitted 5 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 25 pages, 10 figures, to appear in A&A

Journal ref: A&A 685, A75 (2024)

PDRs4All II: JWST's NIR and MIR imaging view of the Orion Nebula

Authors: Emilie Habart, Els Peeters, Olivier Berné, Boris Trahin, Amélie Canin, Ryan Chown, Ameek Sidhu, Dries Van De Putte, Felipe Alarcón, Ilane Schroetter, Emmanuel Dartois, Sílvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Jan Cami, Sara Cuadrado, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (117 additional authors not shown)

Abstract: The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation fron… ▽ More The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate. △ Less

Submitted 2 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: A&A 685, A73 (2024)

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- V: First results on the magnetic field orientation of galaxies

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Mehrnoosh Tahani, Daniel A. Dale, Ignacio del Moral Castro, Julia Roman-Duval, Pamela M. Marcum, John E. Beckman, Kandaswamy Subramanian, Sarah Eftekharzadeh, Leslie Proudfit

Abstract: We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from… ▽ More We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($Ψ_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($ζ$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $Ψ_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($ζ_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($ζ_{\rm{154μm}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields. △ Less

Submitted 13 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

The Quest for the Missing Dust: II -- Two Orders of Magnitude of Evolution in the Dust-to-Gas Ratio Resolved Within Local Group Galaxies

Authors: Christopher J. R. Clark, Julia C. Roman-Duval, Karl D. Gordon, Caroline Bot, Matthew W. L. Smith, Lea M. Z. Hagen

Abstract: We explore evolution in the dust-to-gas ratio with density within four well-resolved Local Group galaxies - the LMC, SMC, M31, and M33. We do this using new ${\it Herschel}$ maps, which restore extended emission that was missed by previous ${\it Herschel}$ reductions. This improved data allows us to probe the dust-to-gas ratio across 2.5 orders of magnitude in ISM surface density. We find signific… ▽ More We explore evolution in the dust-to-gas ratio with density within four well-resolved Local Group galaxies - the LMC, SMC, M31, and M33. We do this using new ${\it Herschel}$ maps, which restore extended emission that was missed by previous ${\it Herschel}$ reductions. This improved data allows us to probe the dust-to-gas ratio across 2.5 orders of magnitude in ISM surface density. We find significant evolution in the dust-to-gas ratio, with dust-to-gas varying with density within each galaxy by up to a factor 22.4. We explore several possible reasons for this, and our favored explanation is dust grain growth in denser regions of ISM. We find that the evolution of the dust-to-gas ratio with ISM surface density is very similar between M31 and M33, despite their large differences in mass, metallicity, and star formation rate; conversely, we find M33 and the LMC to have very different dust-to-gas evolution profiles, despite their close similarity in those properties. Our dust-to-gas ratios address previous disagreement between UV- and FIR-based dust-to-gas estimates for the Magellanic Clouds, removing the disagreement for the LMC, and considerably reducing it for the SMC - with our new dust-to-gas measurements being factors of 2.4 and 2.0 greater than the previous far-infrared estimates, respectively. We also observe that the dust-to-gas ratio appears to fall at the highest densities for the LMC, M31, and M33; this is unlikely to be an actual physical phenomenon, and we posit that it may be due to a combined effect of dark gas, and changing dust mass opacity. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: Accepted for publication in the Astrophysical Journal

Far-ultraviolet Dust Extinction and Molecular Hydrogen in the Diffuse Milky Way Interstellar Medium

Authors: Dries Van De Putte, Stefan I. B. Cartledge, Karl D. Gordon, Geoffrey C. Clayton, Julia Roman-Duval

Abstract: We aim to compare variations in the full-UV dust extinction curve (912-3000 Angstrom), with the HI/H$_2$/total H content along diffuse Milky Way sightlines, to investigate possible connections between ISM conditions and dust properties. We combine an existing sample of 75 UV extinction curves based on IUE and FUSE data, with atomic and molecular column densities measured through UV absorption. The… ▽ More We aim to compare variations in the full-UV dust extinction curve (912-3000 Angstrom), with the HI/H$_2$/total H content along diffuse Milky Way sightlines, to investigate possible connections between ISM conditions and dust properties. We combine an existing sample of 75 UV extinction curves based on IUE and FUSE data, with atomic and molecular column densities measured through UV absorption. The H$_2$ column density data are based on existing Lyman-Werner absorption band models from earlier work on the extinction curves. Literature values for the HI column density were compiled, and improved for 23 stars by fitting a Ly$α$ profile to archived spectra. We discover a strong correlation between the H$_2$ column and the far-UV extinction, and the underlying cause is a linear relationship between H$_2$ and the strength of the far-UV rise feature. This extinction does not scale with HI, and the total H column scales best with $A(V)$ instead. The carrier of the far-UV rise therefore coincides with molecular gas, and further connections are shown by comparing the UV extinction features to the molecular fraction. Variations in the gas-to-extinction ratio $N(\rm{H})/A(V)$ correlate with the UV-to-optical extinction ratio, and we speculate this could be due to coagulation or shattering effects. Based on the H$_2$ temperature, the strongest far-UV rise strengths are found to appear in colder and denser sightlines. △ Less

Submitted 12 October, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 24 pages, 9 figures, accepted for publication in ApJ; fix missing text in generated pdf due to broken tex command

The 30 Doradus Molecular Cloud at 0.4 pc Resolution with the Atacama Large Millimeter/submillimeter Array: Physical Properties and the Boundedness of CO-emitting Structures

Authors: Tony Wong, Luuk Oudshoorn, Eliyahu Sofovich, Alex Green, Charmi Shah, Rémy Indebetouw, Margaret Meixner, Alvaro Hacar, Omnarayani Nayak, Kazuki Tokuda, Alberto D. Bolatto, Mélanie Chevance, Guido De Marchi, Yasuo Fukui, Alec S. Hirschauer, K. E. Jameson, Venu Kalari, Vianney Lebouteiller, Leslie W. Looney, Suzanne C. Madden, Toshikazu Onishi, Julia Roman-Duval, Mónica Rubio, A. G. G. M. Tielens

Abstract: We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with pr… ▽ More We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with previous studies, we find that the central region of the cloud has higher line widths at fixed size relative to the rest of the molecular cloud and to other LMC clouds, indicating an enhanced level of turbulent motions. However, there is no clear trend in gravitational boundedness (as measured by the virial parameter) with distance from R136. Structures observed in $^{13}$CO are spatially coincident with filaments and are close to a state of virial equilibrium. In contrast, CO structures vary greatly in virialization, with low CO surface brightness structures outside of the main filamentary network being predominantly unbound. The low surface brightness structures constitute ~10% of the measured CO luminosity; they may be shredded remnants of previously star-forming gas clumps, or alternatively the CO-emitting parts of more massive, CO-dark structures. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 24 pages; published in ApJ on 15 Jun 2022; associated data can be found at https://doi.org/10.13012/B2IDB-1671495_V1

Journal ref: ApJ, 932, 47 (2022)

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. IV. Calibration of Dust Depletions vs Abundance Ratios in the Milky Way and Magellanic Clouds and Application to Damped Lyman-alpha Systems

Authors: Julia Roman-Duval, Edward B. Jenkins, Kirill Tchernyshyov, Christopher J. R. Clark, Annalisa De Cia, Karl D. Gordon, Aleksandra Hamanowicz, Vianney Lebouteiller, Marc Rafelski, Karin Sandstrom, Jessica Werk, Petia Yanchulova Merica-Jones

Abstract: The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase a… ▽ More The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLA systems. We find that the fraction of metals in dust, the dust-to-gas-ratio, and total abundances are 2-5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance vs. metallicity remains only slightly sub-linear for all existing depletion calibrations, contrary to what is inferred from FIR, 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20\% solar will be needed to resolve this tension. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 37 pages; 22 figures

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

Authors: Enrique Lopez-Rodriguez, Sui Ann Mao, Rainer Beck, Alejandro S. Borlaff, Evangelia Ntormousi, Konstantinos Tassis, Daniel A. Dale, Julia Roman-Duval, Kandaswamy Subramanian, Sergio Martin-Alvarez, Pamela M. Marcum, Susan E. Clark, William T. Reach, Doyal A. Harper, Ellen G. Zweibel

Abstract: We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation an… ▽ More We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy. △ Less

Submitted 7 July, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 43 pages, 28 figures. Accepted for publication to ApJ

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

Authors: Enrique Lopez-Rodriguez, Melanie Clarke, Sachin Shenoy, William Vacca, Simon Coude, Ryan Arneson, Peter Ashton, Sarah Eftekharzadeh, Rainer Beck, John E. Beckman, Alejandro S. Borlaff, Susan E. Clark, Daniel A. Dale, Sergio Martin-Alvarez, Evangelia Ntormousi, William T. Reach, Julia Roman-Duval, Konstantinos Tassis, Doyal A. Harper, Pamela M. Marcum

Abstract: We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipe… ▽ More We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. III. Interstellar Depletions, Dust-to-Metal, and Dust-to-Gas Ratios Versus Metallicity

Authors: J. Roman-Duval, E. B. Jenkins, K. Tchernyshyov, C. J. R. Clark, A. De Cia, K. D. Gordon, A. Hamanowicz, V. Lebouteiller, M. Rafelski, K. Sandstrom, J. Werk, P. Yanchulova Merica-Jones

Abstract: The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of re… ▽ More The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 26 pages, 12 figures

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Authors: Olivier Berné, Émilie Habart, Els Peeters, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Emeric Bron, Jan Cami, Stéphanie Cazaux, Emmanuel Dartois, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Yoko Okada, Takashi Onaka, Massimo Robberto, Markus Röllig, Alexander G. G. M. Tielens, Silvia Vicente, Mark G. Wolfire, Felipe Alarcon, C. Boersma, Ameélie Canin, Ryan Chown, Daniel Dicken , et al. (112 additional authors not shown)

Abstract: Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation… ▽ More Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: Submitted to PASP