-
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.
-
Scylla at APOGEE: The Impact of Starbursts on the Chemical Evolution of the Magellanic Clouds
Authors:
Ivanna Escala,
Kristen B. W. McQuinn,
Sten Hasselquist,
Roger E. Cohen,
James W. Johnson,
Christian R. Hayes,
Clare Burhenne,
Alexander P. Ji,
Christina W. Lindberg,
Petia Yanchulova Merica-Jones,
Yumi Choi,
Andrew E. Dolphin,
Benjamin F. Williams,
Claire E. Murray
Abstract:
Owing to their proximity to the Milky Way, the Large and Small Magellanic Clouds (L/SMC) uniquely probe the evolution of low-mass galaxies undergoing mutual interactions. In this work, we investigate the connection between the star formation histories (SFHs) of the L/SMC measured from HST imaging in the Scylla survey and APOGEE chemical abundances. We model the chemical evolution of the L/SMC in t…
▽ More
Owing to their proximity to the Milky Way, the Large and Small Magellanic Clouds (L/SMC) uniquely probe the evolution of low-mass galaxies undergoing mutual interactions. In this work, we investigate the connection between the star formation histories (SFHs) of the L/SMC measured from HST imaging in the Scylla survey and APOGEE chemical abundances. We model the chemical evolution of the L/SMC in the [Mg/Fe]-[Fe/H] plane within a robust statistical framework to predict chemical abundance signatures resulting directly from starbursts in Scylla SFHs. Both the L/SMC rapidly enrich to high metallicity ([Fe/H] $\gtrsim$ $-1$) within 3 Gyr, followed by slower chemical evolution regulated by sequential starbursts, where the SMC may require higher Fe yields from Type Ia supernovae than the LMC. We also model the [Mg/Fe]-[Fe/H] plane to infer starburst properties across distinct spatial regions in the L/SMC. We identify dominant starbursts in the L/SMC with broadly similar timing, though the SMC may host an earlier burst, and larger burst strength in the LMC. The global starburst properties are nearly uniform across the LMC disk, whereas the dominant SMC population experiences a stronger and later-onset burst in its eastern wing compared to the main body. We also find evidence for a chemically distinct secondary population in the SMC that preferentially traces the foreground and may originate from the LMC. We discuss the implications of these results for the evolutionary history of the L/SMC and for starbursts in interacting low-mass galaxy pairs.
△ Less
Submitted 26 March, 2026;
originally announced March 2026.
-
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.
-
Crowning the Queen: Membership, Age, Rotation, and Activity for the Open Cluster Coma Berenices
Authors:
M. A. Agüeros,
J. L. Curtis,
A. Núñez,
C. Burhenne,
P. Rothstein,
B. J. Shaham,
K. Singh,
P. Bergeron,
M. Kilic,
K. R. Covey,
S. T. Douglas
Abstract:
Despite being only 85 pc away, the open cluster Coma Berenices (Coma Ber) has not been extensively studied. This is due in part to its sparseness and low proper motion, which together made Coma Ber's membership challenging to establish. Gaia data for $\approx$400 previously cataloged candidate cluster stars allowed us to identify $\approx$300 as members. With [Fe/H] measurements for nine members,…
▽ More
Despite being only 85 pc away, the open cluster Coma Berenices (Coma Ber) has not been extensively studied. This is due in part to its sparseness and low proper motion, which together made Coma Ber's membership challenging to establish. Gaia data for $\approx$400 previously cataloged candidate cluster stars allowed us to identify $\approx$300 as members. With [Fe/H] measurements for nine members, we found that Coma Ber has a solar metallicity, and then fit isochrones to its color--magnitude diagram to determine that it is 675$\pm$100 Myr old. With photometry obtained by the Transiting Exoplanet Survey Satellite (TESS) and Zwicky Transient Facility (ZTF), we measured rotation periods for 137 of Coma Ber's low-mass stars, increasing the sample of members with measured periods by a factor of six, and extending the rotational census for the cluster from its late F stars through to its fully convective M dwarfs. By measuring the equivalent width of the H$α$ line for $\approx$250 stars and collecting X-ray detections for $\approx$100 ($\approx$85% and $\approx$33% of the cluster's members, respectively), we characterized magnetic activity in Coma Ber and examined the dependence of chromospheric and of coronal activity on rotation in these stars. Despite having a metallicity that is 0.2 dex below that of their coeval cousins in Praesepe and the Hyades, low-mass stars in Coma Ber seem to follow a similar rotation--activity relation. In detail, however, there are differences that may provide further insight into the impact of metallicity on this still poorly understood relation.
△ Less
Submitted 3 September, 2025;
originally announced September 2025.
-
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.
-
Scylla III. The Outside-In Radial Age Gradient in the Small Magellanic Cloud and the Star Formation Histories of the Main Body, Wing and Outer Regions
Authors:
Roger E. Cohen,
Kristen B. W. McQuinn,
Claire E. Murray,
Benjamin F. Williams,
Yumi Choi,
Christina W. Lindberg,
Clare Burhenne,
Karl D. Gordon,
Petia Yanchulova Merica-Jones,
Caroline Bot,
Andrew E. Dolphin,
Karoline M. Gilbert,
Steven Goldman,
Alec S. Hirschauer,
Karin M. Sandstrom,
O. Grace Telford
Abstract:
The proximity of the Large and Small Magellanic Clouds (LMC and SMC) provides the opportunity to study the impact of dwarf-dwarf interactions on their mass assembly with a unique level of detail. To this end, we analyze two-filter broadband imaging of 83 Hubble Space Telescope (HST) pointings covering 0.203 deg$^2$ towards the SMC, extending out to $\sim$3.5 kpc in projection from its optical cent…
▽ More
The proximity of the Large and Small Magellanic Clouds (LMC and SMC) provides the opportunity to study the impact of dwarf-dwarf interactions on their mass assembly with a unique level of detail. To this end, we analyze two-filter broadband imaging of 83 Hubble Space Telescope (HST) pointings covering 0.203 deg$^2$ towards the SMC, extending out to $\sim$3.5 kpc in projection from its optical center. Lifetime star formation histories (SFHs) fit to each pointing independently reveal an outside-in age gradient such that fields in the SMC outskirts are older on average. We measure radial gradients of the lookback time to form 90%, 75% and 50% of the cumulative stellar mass for the first time, finding $δ$($τ_{90}$, $τ_{75}$, $τ_{50}$)/$δ$R = (0.61$^{+0.08}_{-0.07}$, 0.65$^{+0.09}_{-0.08}$, 0.82$^{+0.12}_{-0.16}$) Gyr/kpc assuming PARSEC evolutionary models and a commonly used elliptical geometry of the SMC, although our results are robust to these assumptions. The wing of the SMC deviates from this trend, forming 25\% of its cumulative mass over the most recent 3 Gyr due to a best-fit star formation rate that remains approximately constant. Our results are consistent with chemodynamical evidence of a tidally stripped SMC component in the foreground, and imply contributions to the observed SFH from multiple previous LMC-SMC interactions. We also compare our SMC SFH with results from a companion study of the LMC, finding that while the two galaxies present different internal, spatially resolved SFH trends, both the LMC and SMC have similar near-constant lifetime SFHs when viewed globally.
△ Less
Submitted 15 October, 2024;
originally announced October 2024.
-
Scylla II. The Spatially Resolved Star Formation History of the Large Magellanic Cloud Reveals an Inverted Radial Age Gradient
Authors:
Roger E. Cohen,
Kristen B. W. McQuinn,
Claire E. Murray,
Benjamin F. Williams,
Yumi Choi,
Christina W. Lindberg,
Clare Burhenne,
Karl D. Gordon,
Petia Yanchulova Merica-Jones,
Karoline M. Gilbert,
Martha L. Boyer,
Steven Goldman,
Andrew E. Dolphin,
O. Grace Telford
Abstract:
The proximity of the Magellanic Clouds provides the opportunity to study interacting dwarf galaxies near a massive host, and spatial trends in their stellar population properties in particular, with a unique level of detail. The Scylla pure parallel program has obtained deep (80% complete to >1 mag below the ancient main sequence turnoff), homogeneous two-filter Hubble Space Telescope (HST) imagin…
▽ More
The proximity of the Magellanic Clouds provides the opportunity to study interacting dwarf galaxies near a massive host, and spatial trends in their stellar population properties in particular, with a unique level of detail. The Scylla pure parallel program has obtained deep (80% complete to >1 mag below the ancient main sequence turnoff), homogeneous two-filter Hubble Space Telescope (HST) imaging sampling the inner star-forming disk of the Large Magellanic Cloud (LMC), the perfect complement to shallower, contiguous ground-based surveys. We harness this imaging together with extant archival data and fit lifetime star formation histories (SFHs) to resolved color-magnitude diagrams (CMDs) of 111 individual fields, using three different stellar evolutionary libraries. We validate per-field recovered distances and extinctions as well as the combined global LMC age-metallicity relation and SFH against independent estimates. We find that the present-day radial age gradient reverses from an inside-out gradient in the inner disk to an outside-in gradient beyond $\sim$2 disk scalelengths, supported by ground-based measurements. The gradients become relatively flatter at earlier lookback times, while the location of the inversion remains constant over an order of magnitude in lookback time, from $\sim$1$-$10 Gyr. This suggests at least one mechanism that predates the recent intense LMC-SMC interaction. We compare observed radial age trends to other late-type galaxies at fixed stellar mass and discuss similarities and differences in the context of potential drivers, implying strong radial migration in the LMC.
△ Less
Submitted 15 October, 2024;
originally announced October 2024.
-
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.