Showing 1–50 of 210 results for author: Williams, B F

Search and Study of the Brightest Stars in the Galaxy IC 342

Authors: O. N. Sholukhova, N. A. Tikhonov, Y. N. Solovyeva, A. N. Sarkisian, A. S. Vinokurov, A. T. Valcheva, P. L. Nedialkov, D. V. Bizyaev, B. F. Williams, V. D. Ivanov

Abstract: We have selected candidate massive stars in the galaxy IC 342 based on archival images from the Hubble Space Telescope and images from the 2 m telescope at the Natioal Astronomical Observatory Rozhen, Bulgaria. Spectral observations of 24 out of 27 selected stars are carried out with the 6 m BTA telescope at the SAO RAS and with the 3.5 m Apache Point Observatory telescope (USA) as part of the pro… ▽ More We have selected candidate massive stars in the galaxy IC 342 based on archival images from the Hubble Space Telescope and images from the 2 m telescope at the Natioal Astronomical Observatory Rozhen, Bulgaria. Spectral observations of 24 out of 27 selected stars are carried out with the 6 m BTA telescope at the SAO RAS and with the 3.5 m Apache Point Observatory telescope (USA) as part of the program for searching bright massive stars in galaxies outside the Local Group. Our analysis reveals that 12 objects have spectra lacking prominent features, except for the emission lines of the surrounding nebulae and are identified as single supergiants of classes O9 to F5 or spatially unresolved young compact clusters. One source with an absorption spectrum probably belongs to our Galaxy. The spectra of seven other objects show features typical of Wolf-Rayet stars or compact clusters containing Wolf-Rayet stars. Another source is a compact supernova remnant. Two other objects are tentatively classified as cold LBV candidates, and one object is classified as a B[e]-supergiant candidate. △ Less

Submitted 21 June, 2024; originally announced June 2024.

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). VI. The High-Mass Stellar Initial Mass Function of M33

Authors: Tobin M. Wainer, Benjamin F. Williams, L. Clifton Johnson, Daniel R. Weisz, Julianne J. Dalcanton, Anil C. Seth, Andrew Dolphin, Meredith J. Durbin, Eric F. Bell, Zhuo Chen, Puragra Guhathakurta, Eric W. Koch, Christina W. Lindberg, Erik Rosolowsky, Karin M. Sandstrom, Evan D. Skillman, Adam Smercina, Estephani E. TorresVillanueva

Abstract: We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). W… ▽ More We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overlineΓ = 1.49\pm0.18$, with an intrinsic scatter of $σ^{2}_Γ = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overlineΓ$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($Γ= 1.30$) IMF, as well as the Salpeter ($Γ= 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 9 Figures, 1 Table

The Vela Pulsar Progenitor Was Most Likely a Binary Merger

Authors: Jeremiah W. Murphy, Andres F. Barrientos, Rene Andrae, Joseph Guzman, Benjamin F. Williams, Julianne J. Dalcanton, Brad Koplitz

Abstract: Stellar evolution theory restricted to single stars predicts a minimum mass for core-collapse supernovae (CCSNe) of around eight solar masses; this minimum mass corresponds to a maximum age of around 45 million years for the progenitor and the coeval population of stars. Binary evolution complicates this prediction. For example, an older stellar population around 100 million years could contain st… ▽ More Stellar evolution theory restricted to single stars predicts a minimum mass for core-collapse supernovae (CCSNe) of around eight solar masses; this minimum mass corresponds to a maximum age of around 45 million years for the progenitor and the coeval population of stars. Binary evolution complicates this prediction. For example, an older stellar population around 100 million years could contain stellar mergers that reach the minimum mass for core collapse. Despite this clear prediction by binary evolution, there are few, if any CCSNe associated with a distinctly older stellar population...until now. The stellar population within 150 pc of the Vela Pulsar is inconsistent with single-star evolution only; instead, the most likely solution is that the stellar population is $\ge$80 Myr old, and the brightest stars are mass gainers and/or mergers, the result of binary evolution. The evidence is as follows. Even though the main sequence is clearly dominated by a $\ge$80-Myr-old population, a large fraction of the corresponding red giants is missing. The best-fitting single-star model expects 51.5 red giants, yet there are only 22; the Poisson probability of this is $1.7 \times 10^{-6}$. In addition, there is an overabundance of bright, young-looking stars (25-30 Myrs old), yet there is not a corresponding young main sequence (MS). Upon closer inspection, the vast majority of the young-looking stars show either past or current signs of binary evolution. These new results are possible due to exquisite Gaia parallaxes and a new age-dating software called {\it Stellar Ages}. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 13 pages, 7 figures, 1 table. submitted to ApJ

The JWST Resolved Stellar Populations Early Release Science Program VII. Stress Testing the NIRCam Exposure Time Calculator

Authors: A. Savino, M. Gennaro, A. E. Dolphin, D. R. Weisz, M. Correnti, J. Anderson, R. Beaton, M. L. Boyer, R. E. Cohen, A. A. Cole, M. J. Durbin, C. T. Garling, M. C. Geha, K. M. Gilbert, J. Kalirai, N. Kallivayalil, K. B. W. McQuinn, M. J. B. Newman, H. Richstein, E. D. Skillman, J. T. Warfield, B. F. Williams

Abstract: We empirically assess estimates from v3.0 of the JWST NIRCam Exposure Time Calculator (ETC) using observations of resolved stars in Local Group targets taken as part of the Resolved Stellar Populations Early Release Science (ERS) Program. For bright stars, we find that: (i) purely Poissonian estimates of the signal-to-noise ratio (SNR) are in good agreement between the ETC and observations, but no… ▽ More We empirically assess estimates from v3.0 of the JWST NIRCam Exposure Time Calculator (ETC) using observations of resolved stars in Local Group targets taken as part of the Resolved Stellar Populations Early Release Science (ERS) Program. For bright stars, we find that: (i) purely Poissonian estimates of the signal-to-noise ratio (SNR) are in good agreement between the ETC and observations, but non-ideal effects (e.g., flat field uncertainties) are the current limiting factor in the photometric precision that can be achieved; (ii) source position offsets, relative to the detector pixels, have a large impact on the ETC saturation predictions and introducing sub-pixel dithers in the observation design can improve the saturation limits by up to ~1 mag. For faint stars, for which the sky dominates the error budget, we find that the choice in ETC extraction strategy (e.g., aperture size relative to point spread function size) can affect the exposure time estimates by up to a factor of 5. We provide guidelines for configuring the ETC aperture photometry to produce SNR predictions in line with the ERS data. Finally, we quantify the effects of crowding on the SNRs over a large dynamic range in stellar density and provide guidelines for approximating the effects of crowding on SNRs predicted by the ETC. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ. 18 pages, 9 Figures, 2 Tables

A Panchromatic Study of the X-ray Binary Population in NGC 300 on Sub-Galactic Scales

Authors: Breanna A. Binder, Rosalie Williams, Jacob Payne, Michael Eracleous, Alexander Belles, Benjamin F. Williams

Abstract: The population-wide properties and demographics of extragalactic X-ray binaries (XRBs) correlate with the star formation rates (SFRs), stellar masses ($M_{\star}$), and environmental factors (such as metallicity, $Z$) of their host galaxy. Although there is evidence that XRB scaling relations ($L_X$/SFR for high mass XRBs [HMXBs] and $L_X$/$M_{\star}$ for low mass XRBs [LMXBs]) may depend on metal… ▽ More The population-wide properties and demographics of extragalactic X-ray binaries (XRBs) correlate with the star formation rates (SFRs), stellar masses ($M_{\star}$), and environmental factors (such as metallicity, $Z$) of their host galaxy. Although there is evidence that XRB scaling relations ($L_X$/SFR for high mass XRBs [HMXBs] and $L_X$/$M_{\star}$ for low mass XRBs [LMXBs]) may depend on metallicity and stellar age across large samples of XRB-hosting galaxies, disentangling the effects of metallicity and stellar age from stochastic effects, particularly on subgalactic scales, remains a challenge. We use archival X-ray through IR observations of the nearby galaxy NGC 300 to self-consistently model the broadband spectral energy distribution and examine radial trends in its XRB population. We measure a current ($<$100 Myr) SFR of 0.18$\pm$0.08 $M_{\odot}$ yr$^{-1}$ and $M_{\star}$= $(2.15^{+0.26}_{-0.14})\times10^9$ $M_{\odot}$. Although we measure a metallicity gradient and radially resolved star formation histories that are consistent with the literature, there is a clear excess in the number of X-ray sources below $\sim10^{37}$ erg s$^{-1}$ that are likely a mix of variable XRBs and additional background AGN. When we compare the subgalactic $L_X$/SFR ratios as a function of $Z$ to the galaxy-integrated $L_X$-SFR-$Z$ relationships from the literature, we find that only the regions hosting the youngest ($\lesssim$30 Myr) HMXBs agree with predictions, hinting at time evolution of the $L_X$-SFR-$Z$ relationship. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted to AAS Journals; 22 pages, 12 figures, 6 tables

The JWST Resolved Stellar Populations Early Release Science Program V. DOLPHOT Stellar Photometry for NIRCam and NIRISS

Authors: Daniel R. Weisz, Andrew E. Dolphin, Alessandro Savino, Kristen B. W. McQuinn, Max J. B. Newman, Benjamin F. Williams, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Karin M. Sandstrom, Andrew A. Cole, Jack T. Warfield, Evan D. Skillman, Roger E. Cohen, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, Michael C. Cooper, Julianne J. Dalcanton , et al. (16 additional authors not shown)

Abstract: We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy),… ▽ More We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 30 pages, 17 figures. Accepted to ApJS. Data products to be hosted on MAST. For DOLPHOT/JWST tutorials, see https://dolphot-jwst.readthedocs.io/en/latest/ . For more program and DOLPHOT info, see https://ers-stars.github.io

The JWST Resolved Stellar Populations Early Release Science Program VI. Identifying Evolved Stars in Nearby Galaxies

Authors: Martha L. Boyer, Giada Pastorelli, Léo Girardi, Paola Marigo, Andrew E. Dolphin, Kristen B. W. McQuinn, Max J. B. Newman, Alessandro Savino, Daniel R. Weisz, Benjamin F. Williams, Jay Anderson, Roger E. Cohen, Matteo Correnti, Andrew A. Cole, Marla C. Geha, Mario Gennaro, Nitya Kallivayalil, Evan N. Kirby, Karin M. Sandstrom, Evan D. Skillman, Christopher T. Garling, Hannah Richstein, Jack T. Warfield

Abstract: We present an investigation of evolved stars in the nearby star-forming galaxy WLM, using NIRCam imaging from the JWST resolved stellar populations early-release science (ERS) program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants (RSGs) and thermally-pulsing asymptotic giant branch (TP-AGB) stars from one another, while als… ▽ More We present an investigation of evolved stars in the nearby star-forming galaxy WLM, using NIRCam imaging from the JWST resolved stellar populations early-release science (ERS) program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants (RSGs) and thermally-pulsing asymptotic giant branch (TP-AGB) stars from one another, while also providing a reasonable separation of the primary TP-AGB subtypes: carbon-rich C-type stars and oxygen-rich M-type stars. The classification scheme we present here agrees very well with the well-established Hubble Space Telescope (HST) medium-band filter technique. The ratio of C to M-type stars (C/M) is 0.8$\pm$0.1 for both the new JWST and the HST classifications, which is within one sigma of empirical predictions from optical narrow-band CN and TiO filters. The evolved star colors show good agreement with the predictions from the PARSEC$+$COLIBRI stellar evolutionary models, and the models indicate a strong metallicity dependence that makes stellar identification even more effective at higher metallicity. However, the models also indicate that evolved star identification with NIRCam may be more difficult at lower metallicies. We test every combination of NIRCam filters using the models and present additional filters that are also useful for evolved star studies. We also find that $\approx$90\% of the dusty evolved stars are carbon-rich, suggesting that carbonaceous dust dominates the present-day dust production in WLM, similar to the findings in the Magellanic Clouds. These results demonstrate the usefulness of NIRCam in identifying and classifying dust-producing stars without the need for mid-infrared data. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 15 page, 12 figures, submitted to AAS Journals

JWST Reveals Star Formation Across a Spiral Arm in M33

Authors: Joshua Peltonen, Erik Rosolowsky, Thomas G. Williams, Eric W. Koch, Andrew Dolphin, Jeremy Chastenet, Julianne J. Dalcanton, Adam Ginsburg, L. Clifton Johnson, Adam K. Leroy, Theo Richardson, Karin M. Sandstrom, Sumit K. Sarbadhicary, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI… ▽ More Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI imaging mosaics at 5.6 and 21 microns that cover a significant portion of one of M33's spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep ALMA CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new DOLPHOT MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further shows YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass. △ Less

Submitted 18 December, 2023; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 12 pages, 10 figures, 1 tables, accepted for publication at MNRAS

Monitoring the X-ray Variability of Bright X-ray Sources in M33

Authors: Rebecca Kyer, Shelby Albrecht, Benjamin F. Williams, Kyros Hinton, Breanna Binder, Margaret Lazzarini, Kristen Garofali, Bret Lehmer, Michael Eracleous, Paul P. Plucinsky, Vallia Antoniou

Abstract: We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected i… ▽ More We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected in previous X-ray catalogs of M33 and discuss their possible nature. The final catalog includes 26 known HMXB candidates identified in the literature. We extend the baseline of the X-ray light curves up to 21 years by including archival X-ray observations of these sources. We compare the detection and non-detection epochs of the sources to suites of simulated source duty cycles and infer that most of our detected sources have duty cycles > 30%. We find only four sources whose detection patterns are consistent with having duty cycles below 30%. This large fraction of sources with high duty cycles is unexpected for a population of HMXBs, thus more frequent X-ray monitoring will likely reveal many more low duty cycle HMXBs in M33. △ Less

Submitted 5 December, 2023; originally announced December 2023.

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

The JWST Resolved Stellar Populations Early Release Science Program IV: The Star Formation History of the Local Group Galaxy WLM

Authors: Kristen. B. W. McQuinn, Max J. B. Newman, Alessandro Savino, Andrew E. Dolphin, Daniel R. Weisz, Benjamin F. Williams, Martha L. Boyer, Roger E. Cohen, Matteo Correnti, Andrew A. Cole, Marla C. Geha, Mario Gennaro, Nitya Kallivayalil, Karin M. Sandstrom, Evan D. Skillman, Jay Anderson, Alberto Bolatto, Michael Boylan-Kolchin, Christopher T. Garling, Karoline M. Gilbert, Leo Girardi, Jason S. Kalirai, Alessandro Mazzi, Giada Pastorelli, Hannah Richstein , et al. (1 additional authors not shown)

Abstract: We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galax… ▽ More We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering a SFH and AMR with the NIRCam F090W-F150W filter combination and provide validation of the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation re-ignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo~A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ~5 Gyr; the stellar mass then rapidly doubled in ~2.5 Gyr, followed by constant star formation over the last ~5 Gyr. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 27 pages, 10 figures, 1 table

TREX: Kinematic Characterisation of a High-Dispersion Intermediate-Age Stellar Component in M33

Authors: L. R. Cullinane, Karoline M. Gilbert, Puragra Guhathakurta, A. C. N. Quirk, Ivanna Escala, Adam Smercina, Benjamin F. Williams, Erik Tollerud, Jessamine Qu, Kaela McConnell

Abstract: The dwarf galaxy Triangulum (M33) presents an interesting testbed for studying stellar halo formation: it is sufficiently massive so as to have likely accreted smaller satellites, but also lies within the regime where feedback and other "in-situ" formation mechanisms are expected to play a role. In this work, we analyse the line-of-sight kinematics of stars across M33 from the TREX survey with a v… ▽ More The dwarf galaxy Triangulum (M33) presents an interesting testbed for studying stellar halo formation: it is sufficiently massive so as to have likely accreted smaller satellites, but also lies within the regime where feedback and other "in-situ" formation mechanisms are expected to play a role. In this work, we analyse the line-of-sight kinematics of stars across M33 from the TREX survey with a view to understanding the origin of its halo. We split our sample into two broad populations of varying age, comprising 2032 "old" red giant branch (RGB) stars, and 671 "intermediate-age" asymptotic giant branch (AGB) and carbon stars. We find decisive evidence for two distinct kinematic components in both old and intermediate-age populations: a low-dispersion (~22 km/s) disk-like component co-rotating with M33's HI gas, and a significantly higher-dispersion component (~50-60 km/s) which does not rotate in the same plane as the gas and is thus interpreted as M33's stellar halo. While kinematically similar, the fraction of stars associated with the halo component differs significantly between the two populations: this is consistently ~10% for the intermediate age population, but decreases from ~34% to ~10% as a function of radius for the old population. We additionally find evidence that the intermediate-age halo population is systematically offset from the systemic velocity of M33 by ~25 km/s, with a preferred central LOS velocity of ~-155 km/s. This is the first detection and characterisation of an intermediate-age halo in M33, and suggests in-situ formation mechanisms, as well as potentially tidal interactions, have helped shaped it. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: 27 pages, 8 figures. Accepted by ApJ

Considering the Single and Binary Origins of the Type IIP SN 2017eaw

Authors: K. Azalee Bostroem, Emmanouil Zapartas, Brad Koplitz, Benjamin F. Williams, Debby Tran, Andrew Dolphin

Abstract: Current population synthesis modeling suggests that 30-50% of Type II supernovae originate from binary progenitors, however, the identification of a binary progenitor is challenging. One indicator of a binary progenitor is that the surrounding stellar population is too old to contain a massive single star.Measurements of the progenitor mass of SN 2017eaw are starkly divided between observations ma… ▽ More Current population synthesis modeling suggests that 30-50% of Type II supernovae originate from binary progenitors, however, the identification of a binary progenitor is challenging. One indicator of a binary progenitor is that the surrounding stellar population is too old to contain a massive single star.Measurements of the progenitor mass of SN 2017eaw are starkly divided between observations made temporally close to core-collapse which show a progenitor mass of 13-15 solar masses (final helium core mass of 4.4 to 6.0 solar masses - which is a more informative property than initial mass) and those from the stellar population surrounding the SN which find M<10.8 solar masses (helium core mass <3.4 solar masses). In this paper, we reanalyze the surrounding stellar population with improved astrometry and photometry, finding a median age of 16.8 (+3.2, -1.0) Myr for all stars younger than 50 Myr (helium core mass of 4.7 solar masses) and 85.9 (+3.2, -6.5) Myr for stars younger than 150 Myr. 16.8 Myr is now consistent with the helium core mass range derived from the temporally near explosion observations for single stars. Applying the combined constraints to population synthesis models, we determine that the probability of the progenitor of SN 2017eaw being an initially single-star is 65% compared to 35% for prior binary interaction. 85.9 Myr is inconsistent with any formation scenarios. We demonstrate that combining progenitor age constraints with helium core mass estimates from red supergiant SED modeling, late-time spectra, and indirectly from light curve modeling can help to differentiate single and binary progenitor scenarios and provide a framework for the application of this technique to future observations. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: Accepted to AJ

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). V. The Structure of M33 in Resolved Stellar Populations

Authors: Adam Smercina, Julianne J. Dalcanton, Benjamin F. Williams, Meredith J. Durbin, Margaret Lazzarini, Eric F. Bell, Yumi Choi, Andrew Dolphin, Karoline Gilbert, Puragra Guhathakurta, Eric W. Koch, Hans-Walter Rix, Erik Rosolowsky, Anil Seth, Evan D. Skillman, Daniel R. Weisz

Abstract: We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$… ▽ More We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$ completely distinct from what is seen in broadband optical imaging, and the first-ever confirmation of a bar in M33. We estimate a bar extent of $\sim$1 kpc. The two spiral arms are asymmetric in orientation and strength, and likely represent the innermost impact of the recent tidal interaction responsible for M33's warp at larger scales. The flocculent multi-armed morphology for which M33 is known is only visible in the young upper main sequence population, which closely tracks the morphology of the ISM. We investigate the stability of M33's disk, finding $Q{\sim}1$ over the majority of the disk. We fit multiple components to the old stellar density distribution and find that, when considering recent stellar kinematics, M33's bulk structure favors the inclusion of an accreted halo component, modeled as a broken power-law. The best-fit halo model has an outer power-law index of $-$3 and accurately describes observational evidence of M33's stellar halo from both resolved stellar spectroscopy in the disk and its stellar populations at large radius. Integrating this profile yields a total halo stellar mass of ${\sim}5{\times}10^8\ M_{\odot}$, giving a total stellar halo mass fraction of 16%, most of which resides in the innermost 2.5 kpc. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

The Panchromatic Hubble Andromeda Treasury XXI. The Legacy Resolved Stellar Photometry Catalog

Authors: Benjamin F. Williams, Meredith Durbin, Dustin Lang, Julianne J. Dalcanton, Andrew E. Dolphin, Adam Smercina, Petia Yanchulova Merica-Jones, Daniel R. Weisz, Eric F. Bell, Karoline M. Gilbert, Leo Girardi, Karl Gordon, Puragra Guhathakurta, L. Clifton Johnson, Tod R. Lauer, Anil Seth, Evan Skillman

Abstract: We present the final legacy version of stellar photometry for the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We have reprocessed all of the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) imaging from the PHAT survey using an improved method that optimi… ▽ More We present the final legacy version of stellar photometry for the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We have reprocessed all of the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) imaging from the PHAT survey using an improved method that optimized the survey depth and chip gap coverage by including all overlapping exposures in all bands in the photometry. An additional improvement was gained through the use of charge transfer efficiency (CTE) corrected input images, which provide more complete star finding as well as more reliable photometry for the NUV bands, which had no CTE correction in the previous version of the PHAT photometry. While this method requires significantly more computing resources and time than earlier versions where the photometry was performed on individual pointings, it results in smaller systematic instrumental completeness variations as demonstrated by cleaner maps in stellar density, and it results in optimal constraints on stellar fluxes in all bands from the survey data. Our resulting catalog has 138 million stars, 18% more than the previous catalog, with lower density regions gaining as much as 40% more stars. The new catalog produces nearly seamless population maps which show relatively well-mixed distributions for populations associated with ages older than 1-2 Gyr, and highly structured distributions for the younger populations. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 28 pages, 4 tables, 18 figures, accepted for publication in The Astrophysical Journal Supplement Series

Extragalactic Star Cluster Science with the Nancy Grace Roman Space Telescope's High Latitude Wide Area Survey and the Vera C. Rubin Observatory

Authors: Kristen C. Dage, Christopher Usher, Jennifer Sobeck, Ana L. Chies Santos, Róbert Szabó, Marta Reina-Campos, Léo Girardi, Vincenzo Ripepi, Marcella Di Criscienzo, Ata Sarajedini, Will Clarkson, Peregrine McGehee, John Gizis, Katherine Rhode, John Blakeslee, Michele Cantiello, Christopher A. Theissen, Annalisa Calamida, Ana Ennis, Nushkia Chamba, Roman Gerasimov, R. Michael Rich, Pauline Barmby, Annette M. N. Ferguson, Benjamin F. Williams

Abstract: The Nancy Grace Roman Telescope's High Latitude Wide Area Survey will have a number of synergies with the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST), particularly for extragalactic star clusters. Understanding the nature of star clusters and star cluster systems are key topics in many areas of astronomy, chief among them stellar evolution, high energy astrophysics, galaxy asse… ▽ More The Nancy Grace Roman Telescope's High Latitude Wide Area Survey will have a number of synergies with the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST), particularly for extragalactic star clusters. Understanding the nature of star clusters and star cluster systems are key topics in many areas of astronomy, chief among them stellar evolution, high energy astrophysics, galaxy assembly/dark matter, the extragalactic distance scale, and cosmology. One of the challenges will be disentangling the age/metallicity degeneracy because young ($\sim$Myr) metal-rich clusters have similar SEDs to old ($\sim$Gyr) metal-poor clusters. Rubin will provide homogeneous, $ugrizy$ photometric coverage, and measurements in the red Roman filters will help break the age-metallicity and age-extinction degeneracies, providing the first globular cluster samples that cover wide areas while essentially free of contamination from Milky Way stars. Roman's excellent spatial resolution will also allow measurements of cluster sizes. We advocate for observations of a large sample of galaxies with a range of properties and morphologies in the Rubin/LSST footprint matching the depth of the LSST Wide-Fast-Deep field $i$ band limit (26.3 mag), and recommend adding the F213 filter to the survey. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: white paper submitted for Roman CCS input

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

Multiwavelength Characterization of the High Mass X-ray Binary Population of M33

Authors: Margaret Lazzarini, Kyros Hinton, Cheyanne Shariat, Benjamin F. Williams, Kristen Garofali, Julianne J. Dalcanton, Meredith Durbin, Vallia Antoniou, Breanna Binder, Michael Eracleous, Neven Vulic, Jun Yang, Daniel Wik, Aria Gasca, Quetzalcoatl Kuauhtzin

Abstract: We present multi-wavelength characterization of 65 high mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalog, which covers the inner disk of M33 a… ▽ More We present multi-wavelength characterization of 65 high mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalog, which covers the inner disk of M33 at near infrared, optical, and near ultraviolet wavelengths. We perform spectral energy distribution (SED) fitting on multi-band photometry for each point source optical counterpart to measure its physical properties including mass, temperature, luminosity, and radius. We find that the majority of the HMXB companion star candidates are likely B-type main sequence stars, suggesting that the HMXB population of M33 is dominated by Be-XRBs, as is seen in other Local Group galaxies. We use spatially-resolved recent star formation history (SFH) maps of M33 to measure the age distribution of the HMXB candidate sample and the HMXB production rate for M33. We find a bimodal distribution for the HMXB production rate over the last 80 Myr, with a peak at $\sim$10 Myr and $\sim$40 Myr, which match theoretical formation timescales for the most massive HMXBs and Be X-ray binaries (Be-XRBs), respectively. We measure an HMXB production rate of 107$-$136 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 50 Myr and 150$-$199 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 80 Myr. For sources with compact object classifications from overlapping NuSTAR observations, we find a preference for giant/supergiant companion stars in BH-HMXBs and main sequence companion stars in neutron star HMXBs (NS-HMXBs). △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ, 37 pages, 6 figures, 8 tables

Clusters, Clouds, and Correlations: Relating Young Clusters to Giant Molecular Clouds in M33 and M31

Authors: Joshua Peltonen, Erik Rosolowsky, L. Clifton Johnson, Anil C. Seth, Julianne Dalcanton, Eric F. Bell, Jonathan Braine, Eric W. Koch, Margaret Lazzarini, Adam K. Leroy, Evan D. Skillman, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249… ▽ More We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

The Spatial Correlation of High Mass X-ray Binaries and Young Star Clusters in Nearby Star-Forming Galaxies

Authors: Breanna A. Binder, Ashley K. Anderson, Kristen Garofali, Margaret Lazzarini, Benjamin F. Williams

Abstract: We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we… ▽ More We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we find that the peak value of the spatial correlation function correlates strongly with the specific star formation rate of the host galaxy. We additionally perform an X-ray stacking analysis of 211 non-X-ray detected YSCs in M31 and 463 YSCs in M33. We do not detect excess X-ray emission at the stacked cluster locations down to 3$σ$ upper limits of $\sim10^{33}$ erg s$^{-1}$ (0.35-8 keV) in both galaxies, which strongly suggests that dynamical formation within YSCs is not a major HMXB formation channel. We interpret our results in the context of (1) the recent star formation histories of the galaxies, which may produce differences in the demographics of compact objects powering the HMXBs, and (2) the differences in natal kicks experienced by compact objects during formation, which can eject newly-formed HMXB from their birth clusters. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Accepted to MNRAS; 11 pages, 7 figures

The Panchromatic Hubble Andromeda Treasury XX: The Disk of M31 is Thick

Authors: Julianne J. Dalcanton, Eric F. Bell, Yumi Choi, Andrew E. Dolphin, Morgan Fouesneau, Léo Girardi, David W. Hogg, Anil C. Seth, Benjamin F. Williams

Abstract: We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geomet… ▽ More We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geometrical projection leads to differences in the numbers of stars on the near and far sides of the thin dust layer. The fraction of reddened stars will thus differ from the 50% prediction, with a deviation that becomes larger for puffier disks. We map the fraction of reddened red giant branch (RGB) stars across M31, which shows prominent dust lanes on only one side of the major axis. The fraction of reddened stars varies systematically from 20% to 80%, which requires that these stars have an exponential scale height h_z that is 0.14+/-0.015 times the exponential scale length (h_r~5.5kpc). M31's RGB stars must therefore have h_z=770+/-80pc, which is far thicker than the Milky Way's thin disk, but comparable to its thick disk. The lack of a significant thin disk in M31 is unexpected, but consistent with its interaction history and high disk velocity dispersion. We suggest that asymmetric reddening be used as a generic criteria for identifying ``thick disk'' dominated systems, and discuss prospects for future 3-dimensional tomographic mapping of the gas and stars in M31. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 22 pages. Accepted to the Astrophysical Journal

A 1.3% distance to M33 from HST Cepheid photometry

Authors: Louise Breuval, Adam G. Riess, Lucas M. Macri, Siyang Li, Wenlong Yuan, Stefano Casertano, Tarini Konchady, Boris Trahin, Meredith J. Durbin, Benjamin F. Williams

Abstract: We present a low-dispersion period-luminosity relation (PL) based on 154 Cepheids in Messier 33 (M33) with Hubble Space Telescope (HST) photometry from the PHATTER survey. Using high-quality ground-based light curves, we recover Cepheid phases and amplitudes for multi-epoch HST data and we perform template fitting to derive intensity-averaged mean magnitudes. HST observations in the SH0ES near-inf… ▽ More We present a low-dispersion period-luminosity relation (PL) based on 154 Cepheids in Messier 33 (M33) with Hubble Space Telescope (HST) photometry from the PHATTER survey. Using high-quality ground-based light curves, we recover Cepheid phases and amplitudes for multi-epoch HST data and we perform template fitting to derive intensity-averaged mean magnitudes. HST observations in the SH0ES near-infrared Wesenheit system significantly reduce the effect of crowding relative to ground-based data, as seen in the final PL scatter of $σ$ = 0.11 mag. We adopt the absolute calibration of the PL based on HST observations in the Large Magellanic Cloud (LMC) and a distance derived using late-type detached eclipsing binaries to obtain a distance modulus for M33 of $μ$ = 24.622 $\pm$ 0.030 mag (d = 840 $\pm$ 11 kpc), a best-to-date precision of 1.3%. We find very good agreement with past Cepheid-based measurements. Several TRGB estimates bracket our result while disagreeing with each other. Finally, we show that the flux contribution from star clusters hosting Cepheids in M33 does not impact the distance measurement and we find only 3.7% of the sample is located in (or nearby) young clusters. M33 offers one of the best sites for the cross-calibration of many primary distance indicators. Thus, a precise independent geometric determination of its distance would provide a valuable new anchor to measure the Hubble constant. △ Less

Submitted 31 March, 2023; originally announced April 2023.

Comments: Submitted to ApJ, comments welcome

Journal ref: ApJ 951 118 2023

The Masses of Supernova Remnant Progenitors in M33

Authors: Brad Koplitz, Jared Johnson, Benjamin F. Williams, Mariangelly Diaz-Rodriguez, Jeremiah W. Murphy, Margaret Lazzarini, Joseph Guzman, Julianne J. Dalcanton, Andrew Dolphin, Meredith Durbin

Abstract: Using resolved optical stellar photometry from the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey, we measured the star formation history (SFH) near the position of 85 supernova remnants (SNRs) in M33. We constrained the progenitor masses for 60 of these SNRs, finding the remaining 25 remnants had no local SF in the last 56 Myr consistent with core-collapse SNe… ▽ More Using resolved optical stellar photometry from the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey, we measured the star formation history (SFH) near the position of 85 supernova remnants (SNRs) in M33. We constrained the progenitor masses for 60 of these SNRs, finding the remaining 25 remnants had no local SF in the last 56 Myr consistent with core-collapse SNe (CCSNe), making them potential Type Ia candidates. We then infer a progenitor mass distribution from the age distribution, assuming single star evolution. We find that the progenitor mass distribution is consistent with being drawn from a power-law with an index of $-2.9^{+1.2}_{-1.0}$. Additionally, we infer a minimum progenitor mass of $7.1^{+0.1}_{-0.2}\ M_{\odot}$ from this sample, consistent with several previous studies, providing further evidence that stars with ages older than the lifetimes of single 8 $M_{\odot}$ stars are producing supernovae. △ Less

Submitted 1 April, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 20 pages, 7 figures, 2 tables, Accepted at ApJ

A Multiwavelength Classification and Study of Red Supergiant Candidates in NGC 6946

Authors: Jared R. Johnson, Brad Koplitz, Benjamin F. Williams, Julianne J. Dalcanton, Andrew Dolphin, Leo Girardi

Abstract: We have combined resolved stellar photometry from Hubble Space Telescope (\emph{HST}), \emph{Spitzer}, and \emph{Gaia} to identify red supergiant (RSG) candidates in NGC~6946, based on their colors, proper motions, visual morphologies, and spectral energy distributions. We start with a large sample of 17,865 RSG candidates based solely on \emph{HST} near-infrared photometry. We then chose a small… ▽ More We have combined resolved stellar photometry from Hubble Space Telescope (\emph{HST}), \emph{Spitzer}, and \emph{Gaia} to identify red supergiant (RSG) candidates in NGC~6946, based on their colors, proper motions, visual morphologies, and spectral energy distributions. We start with a large sample of 17,865 RSG candidates based solely on \emph{HST} near-infrared photometry. We then chose a small sample of 385 of these candidates with Spitzer matches for more detailed study. Using evolutionary models and isochrones, we isolate a space where RSGs would be found in our photometry catalogs. We then visually inspect each candidate and compare to Gaia catalogs to identify and remove foreground stars. As a result, we classify 95 potential RSGs, with 40 of these being in our highest-quality sample. We fit the photometry of the populations of stars in the regions surrounding the RSGs to infer their ages. Placing our best candidate RSG stars into three age bins between 1 and 30 Myr, we find 27.5\% of the candidates falling between 1-10 Myr, 37.5\% between 10-20 Myr, and 35\% 20-30 Myr. A comparison of our results to the models of massive star evolution shows some agreement between model luminosities and the luminosities of our candidates for each age. Three of our candidates appear significantly more consistent with binary models than single-star evolution models. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 32 pages, 18 figures, 4 tables

The JWST Resolved Stellar Populations Early Release Science Program III: Photometric Star-Galaxy Separations for NIRCam

Authors: Jack T. Warfield, Hannah Richstein, Nitya Kallivayalil, Roger E. Cohen, Alessandro Savino, Martha L. Boyer, Christopher T. Garling, Mario Gennaro, Kristen B. W. McQuinn, Max J. B. Newman, Jay Anderson, Andrew A. Cole, Matteo Correnti, Andrew E. Dolphin, Marla C. Geha, Karin M. Sandstrom, Daniel R. Weisz, Benjamin F. Williams

Abstract: We present criteria for separately classifying stars and unresolved background galaxies in photometric catalogs generated with the point spread function (PSF) fitting photometry software DOLPHOT from images taken of Draco II, WLM, and M92 with the Near Infrared Camera (NIRCam) on JWST. Photometric quality metrics from DOLPHOT in one or two filters can recover a pure sample of stars. Conversely, co… ▽ More We present criteria for separately classifying stars and unresolved background galaxies in photometric catalogs generated with the point spread function (PSF) fitting photometry software DOLPHOT from images taken of Draco II, WLM, and M92 with the Near Infrared Camera (NIRCam) on JWST. Photometric quality metrics from DOLPHOT in one or two filters can recover a pure sample of stars. Conversely, colors formed between short-wavelength (SW) and long-wavelength (LW) filters can be used to effectively identify pure samples of galaxies. Our results highlight that the existing DOLPHOT output parameters can be used to reliably classify stars in our NIRCam data without the need to resort to external tools or more complex heuristics. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 5 pages, 1 figure

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

Authors: Daniel R. Weisz, Kristen B. W. McQuinn, Alessandro Savino, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Andrew E. Dolphin, Karin M. Sandstrom, Andrew A. Cole, Benjamin F. Williams, Evan D. Skillman, Roger E. Cohen, Max J. B. Newman, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, M. C. Cooper, Julianne J. Dalcanton, Aaron L. Dotter , et al. (17 additional authors not shown)

Abstract: We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We descr… ▽ More We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

A Panchromatic Study of Massive Stars in the Extremely Metal-Poor Local Group Dwarf Galaxy Leo A

Authors: Maude Gull, Daniel R. Weisz, Peter Senchyna, Nathan R. Sandford, Yumi Choi, Anna F. McLeod, Kareem El-Badry, Ylva Götberg, Karoline M. Gilbert, Martha Boyer, Julianne J. Dalcanton, Puragra GuhaThakurta, Steven Goldman, Paola Marigo, Kristen B. W. McQuinn, Giada Pastorelli, Daniel P. Stark, Evan Skillman, Yuan-sen Ting, Benjamin F. Williams

Abstract: We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region… ▽ More We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry. △ Less

Submitted 28 December, 2022; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 35 pages, 18 figures

A detailed star formation history for the extremely diffuse Andromeda XIX dwarf galaxy

Authors: Michelle L. M. Collins, Benjamin F. Williams, Erik J. Tollerud, Eduardo Balbinot, Karoline M. Gilbert, Andrew Dolphin

Abstract: We present deep imaging of the ultra-diffuse Andromeda XIX dwarf galaxy from the Advance Camera for Surveys on the Hubble Space Telescope which resolves its stellar populations to below the oldest main sequence turn-off. We derive a full star formation history for the galaxy using MATCH, and find no evidence of star formation in the past 8 Gyr. We calculate a quenching time of $τ_{90}=9.7\pm0.2$~G… ▽ More We present deep imaging of the ultra-diffuse Andromeda XIX dwarf galaxy from the Advance Camera for Surveys on the Hubble Space Telescope which resolves its stellar populations to below the oldest main sequence turn-off. We derive a full star formation history for the galaxy using MATCH, and find no evidence of star formation in the past 8 Gyr. We calculate a quenching time of $τ_{90}=9.7\pm0.2$~Gyr, suggesting Andromeda~XIX ceased forming stars very early on. This early quenching, combined with its extremely large half-light radius, low density dark matter halo and lower than expected metallicity make it a unique galaxy within the Local Group and raises questions about how it formed. The early quenching time allows us to rule out feedback from bursty star formation as a means to explain its diffuse stellar population and low density dark matter halo. We find that the extended stellar population, low density halo and star formation could be explained by either tidal interactions (such as tidal shocking) or by late dry mergers, with the latter also explaining its low metallicity. Proper motions and detailed abundances would allow us to distinguish between these two scenarios. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS, Published by Oxford University Press on behalf of the Royal Astronomical Society

Resolved SPLASH Chemodynamics in Andromeda's PHAT Stellar Halo and Disk: On the Nature of the Inner Halo Along the Major Axis

Authors: Ivanna Escala, Amanda C. N. Quirk, Puragra Guhathakurta, Karoline M. Gilbert, J. Leigh Wojno, Lara Cullinane, Benjamin F. Williams, Julianne Dalcanton

Abstract: Stellar kinematics and metallicity are key to exploring formation scenarios for galactic disks and halos. In this work, we characterized the relationship between kinematics and photometric metallicity along the line-of-sight to M31's disk. We combined optical HST/ACS photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey with Keck/DEIMOS spectra from the Spectroscopic and Photome… ▽ More Stellar kinematics and metallicity are key to exploring formation scenarios for galactic disks and halos. In this work, we characterized the relationship between kinematics and photometric metallicity along the line-of-sight to M31's disk. We combined optical HST/ACS photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey with Keck/DEIMOS spectra from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. The resulting sample of 3512 individual red giant branch stars spans 4-19 projected kpc, making it a useful probe of both the disk and inner halo. We separated these stars into disk and halo populations by modeling the line-of-sight velocity distributions as a function of position across the disk region, where $\sim$73% stars have a high likelihood of belonging to the disk and $\sim$14% to the halo. Although stellar halos are typically thought to be metal-poor, the kinematically identified halo contains a significant population of stars ($\sim$29%) with disk-like metallicity ([Fe/H]$_{\rm phot}$ $\sim$ $-0.10$). This metal-rich halo population lags the gaseous disk to a similar extent as the rest of the halo, indicating that it does not correspond to a canonical thick disk. Its properties are inconsistent with those of tidal debris originating from the Giant Stellar Stream merger event. Moreover, the halo is chemically distinct from the phase-mixed component previously identified along the minor axis (i.e., away from the disk), implying contributions from different formation channels. These metal-rich halo stars provide direct chemodynamical evidence in favor of the previously suggested "kicked-up" disk population in M31's inner stellar halo. △ Less

Submitted 12 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Accepted by AJ. Conclusions on page 21. 19 figures, 2 tables, 5 appendices

The JWST Resolved Stellar Populations Early Release Science Program I.: NIRCam Flux Calibration

Authors: Martha L. Boyer, Jay Anderson, Mario Gennaro, Marla Geha, Kristen B. Wingfield McQuinn, Erik Tollerud, Matteo Correnti, Max J. Brenner Newman, Roger E. Cohen, Nitya Kallivayalil, Rachel Beaton, Andrew A. Cole, Andrew Dolphin, Jason S. Kalirai, Karin M. Sandstrom, Alessandro Savino, Evan D. Skillman, Daniel R. Weisz, Benjamin F. Williams

Abstract: We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improve… ▽ More We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improved zeropoints for the ERS filters and show that alternate zeropoints derived by the community also improve the calibration significantly. We also find that the detector offsets appear to be time variable by up to at least 0.1 mag. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: Accepted for publication in RNAAS

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). IV. Star Cluster Catalog

Authors: L. Clifton Johnson, Tobin M. Wainer, Estephani E. TorresVillanueva, Anil C. Seth, Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Daniel R. Weisz, Eric F. Bell, Puragra Guhathakurta, Evan Skillman, Adam Smercina

Abstract: We construct a catalog of star clusters from Hubble Space Telescope images of the inner disk of the Triangulum Galaxy (M33) using image classifications collected by the Local Group Cluster Search, a citizen science project hosted on the Zooniverse platform. We identify 1214 star clusters within the Hubble Space Telescope imaging footprint of the Panchromatic Hubble Andromeda Treasury: Triangulum E… ▽ More We construct a catalog of star clusters from Hubble Space Telescope images of the inner disk of the Triangulum Galaxy (M33) using image classifications collected by the Local Group Cluster Search, a citizen science project hosted on the Zooniverse platform. We identify 1214 star clusters within the Hubble Space Telescope imaging footprint of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. Comparing this catalog to existing compilations in the literature, 68% of the clusters are newly identified. The final catalog includes multi-band aperture photometry and fits for cluster properties via integrated light SED fitting. The cluster catalog's 50% completeness limit is ~1500 solar masses at an age of 100 Myr, as derived from comprehensive synthetic cluster tests. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 24 pages, 15 figures, 9 tables, Accepted to ApJ

walter: A Tool for Predicting Resolved Stellar Population Observations with Applications to the Roman Space Telescope

Authors: Lachlan Lancaster, Sarah Pearson, Benjamin F. Williams, Kathryn V. Johnston, Tjitske K. Starkenburg, Erin Kado-Fong, Anil C. Seth, Eric F. Bell

Abstract: Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the… ▽ More Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the role of accretion in galactic formation. Upcoming wide-field surveys promise to extend this ability to all galaxies within 10~Mpc, drastically increasing our capability to decipher galaxy evolution and enabling statistical studies of galaxies' stellar populations. To facilitate the optimum use of these upcoming capabilities we develop a simple formalism to predict the density of resolved stars for an observation of a stellar population at fixed surface brightness and population parameters. We provide an interface to calculate all quantities of interest to this formalism via a public release of the code: \texttt{walter}. This code enables calculation of (i) the expected number density of detected stars, (ii) the exposure time needed to reach certain population features, such as the horizontal branch, and (iii) an estimate of the crowding limit, among other features. These calculations will be very useful for planning surveys with NASA's upcoming Nancy Grace Roman Space Telescope (Roman, formerly WFIRST), which we use for example calculations throughout this work. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 16 pages, 6 figures, submitted to AJ. Comments welcome

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

Authors: Margaret Lazzarini, Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Adam Smercina, Eric F. Bell, Yumi Choi, Andrew Dolphin, Karoline Gilbert, Puragra Guhathakurta, Erik Rosolowsky, Evan Skillman, O. Grace Telford, Daniel Weisz

Abstract: We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center… ▽ More We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center of M33 along the major and semi-major axes, respectively. We divide the PHATTER optical survey into 2005 regions that measure 24 arcsec, $\sim$100 pc, on a side and fit color magnitude diagrams for each region individually to measure the spatially resolved SFH of M33 within the PHATTER footprint. There are significant fluctuations in the SFH on small spatial scales and also galaxy-wide scales that we measure back to about 630 Myr ago. We observe a more flocculent spiral structure in stellar populations younger than about 80 Myr, while the structure of the older stellar populations is dominated by two spiral arms. We also observe a bar in the center of M33, which dominates at ages older than about 80 Myr. Finally, we find that the mean star formation rate (SFR) over the last 100 Myr within the PHATTER footprint is 0.32$\pm$0.02 M$_{\odot}$ yr$^{-1}$. We measure a current SFR (over the last 10 Myr) of 0.20$\pm$0.03 M$_{\odot}$ yr$^{-1}$. This SFR is slightly higher than previous measurements from broadband estimates, when scaled to account for the fraction of the D25 area covered by the PHATTER survey footprint. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: accepted for publication in ApJ, 37 pages, 18 figures, 6 tables

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

Authors: Sukanya Chakrabarti, Alex Drlica-Wagner, Ting S. Li, Neelima Sehgal, Joshua D. Simon, Simon Birrer, Duncan A. Brown, Rebecca Bernstein, Alberto D. Bolatto, Philip Chang, Kyle Dawson, Paul Demorest, Daniel Grin, David L. Kaplan, Joseph Lazio, Jennifer Marshall, Eric J. Murphy, Scott Ransom, Brant E. Robertson, Rajeev Singh, Anže Slosar, Tommaso Treu, Yu-Dai Tsai, Benjamin F. Williams

Abstract: We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measu… ▽ More We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 31 pages, 6 figures, 2 tables

The Candidate Progenitor Companion Star of the Type Ib/c SN 2013ge

Authors: Ori D. Fox, Schuyler D. Van Dyk, Benjamin F. Williams, Maria Drout, Emmanouil Zapartas, Nathan Smith, Dan Milisavljevic, Jennifer E. Andrews, K. Azalee Bostroem, Alexei V. Filippenko, Sebastian Gomez, Patrick L. Kelly, S. E. de Mink, Justin Pierel, Armin Rest, Stuart Ryder, Niharika Sravan, Lou Strolger, Qinan Wang, Kathryn E. Weil

Abstract: This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source's spectral energy distribution is… ▽ More This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source's spectral energy distribution is most consistent with that of a B5 I supergiant. While binary models tend to predict OB-type stars for stripped-envelope companions, the location of the source on a color-magnitude diagram (CMD) places it redward of its more likely position on the main sequence (MS). The source may be temporarily out of thermal equilibrium, or a cool and inflated non-MS companion, which is similar to the suggested companion of Type Ib SN 2019yvr that was constrained from pre-SN imaging. We also consider other possible physical scenarios for the source, including a fading SN, circumstellar shock interaction, line of site coincidence, and an unresolved host star cluster, all of which will require future observations to more definitively rule out. Ultimately, the fraction of surviving companions ("binary fraction") will provide necessary constraints on binary evolution models and the underlying physics. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted to ApJL. 9 pages, 4 figures, 1 table

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). III. The Mass Function of Young Star Clusters in M33

Authors: Tobin M. Wainer, L. Clifton Johnson, Anil C. Seth, Estephani E. TorresVillanueva, Julianne J. Dalcanton, Meredith J. Durbin, Andrew Dolphin, Daniel R. Weisz, Benjamin F. Williams

Abstract: We measure the star cluster mass function for the Local Group galaxy M33. We use the catalog of stellar clusters selected from the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. We analyze 711 clusters in M33 with $\rm 7.0 < log(Age/yr) < 8.5$, and log($M/M_{\odot}$) $>$ 3.0 as determined from color-magnitude diagram fits to individual stars. The M33 cluster m… ▽ More We measure the star cluster mass function for the Local Group galaxy M33. We use the catalog of stellar clusters selected from the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. We analyze 711 clusters in M33 with $\rm 7.0 < log(Age/yr) < 8.5$, and log($M/M_{\odot}$) $>$ 3.0 as determined from color-magnitude diagram fits to individual stars. The M33 cluster mass function is best described by a Schechter function with power law slope $α= -2.06^{+0.14}_{-0.13}$, and truncation mass log($M_c/M_{\odot}$) $= 4.24^{+0.16}_{-0.13}$. The data show strong evidence for a high-mass truncation, thus strongly favoring a Schechter function fit over a pure power law. M33's truncation mass is consistent with the previously identified linear trend between $M_c$, and star formation rate surface density, \SigSFR. We also explore the effect that individual cluster mass uncertainties have on derived mass function parameters, and find evidence to suggest that large cluster mass uncertainties have the potential to bias the truncation mass of fitted mass functions on the one sigma level. △ Less

Submitted 2 February, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: 18 pages, 15 figures, 1 table, Accepted to ApJ (February 2, 2022)

A Census of Thermally-Pulsing AGB stars in the Andromeda Galaxy and a First Estimate of their Contribution to the Global Dust Budget

Authors: Steven R. Goldman, Martha L. Boyer, Julianne Dalcanton, Iain McDonald, Leo Girardi, Benjamin F. Williams, Sundar Srinivasan, Karl Gordon

Abstract: We present a near-complete catalog of the metal-rich population of Thermally-Pulsing Asymptotic Giant Branch stars in the northwest quadrant of M31. This metal-rich sample complements the equally complete metal-poor Magellanic Cloud AGB catalogs produced by the SAGE program. Our catalog includes HST wide-band photometry from the Panchromatic Hubble Andromeda Treasury survey, HST medium-band photom… ▽ More We present a near-complete catalog of the metal-rich population of Thermally-Pulsing Asymptotic Giant Branch stars in the northwest quadrant of M31. This metal-rich sample complements the equally complete metal-poor Magellanic Cloud AGB catalogs produced by the SAGE program. Our catalog includes HST wide-band photometry from the Panchromatic Hubble Andromeda Treasury survey, HST medium-band photometry used to chemically classify a subset of the sample, and Spitzer mid- and far-IR photometry that we have used to isolate dust-producing AGB stars. We have detected 346,623 AGB stars; these include 4,802 AGB candidates producing considerable dust, and 1,356 AGB candidates that lie within clusters with measured ages, and in some cases metallicities. Using the Spitzer data and chemical classifications made with the medium-band data, we have identified both carbon- and oxygen-rich AGB candidates producing significant dust. We have applied color--mass-loss relations based on dusty AGB stars from the LMC to estimate the dust injection by AGB stars in the PHAT footprint. Applying our color relations to a subset of the chemically-classified stars producing the bulk of the dust, we find that ~97.8% of the dust is oxygen-rich. Using several scenarios for the dust lifetime, we have estimated the contribution of AGB stars to the global dust budget of M31 to be 0.9-35.5%, which is in line with previous estimates in the Magellanic Clouds. Follow-up observations of the M31 AGB candidates with the JWST will allow us to further constrain stellar and chemical evolutionary models, and the feedback and dust production of metal-rich evolved stars. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJS. 35 pages, 25 Figures

The Stellar Age Dependence of X-ray Emission from Normal Star-Forming Galaxies in the GOODS Fields

Authors: Woodrow Gilbertson, Bret Lehmer, Keith Doore, Rafael Eufrasio, Antara Basu-Zych, William Brandt, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Bin Luo, Paolo Tozzi, Fabio Vito, Benjamin F. Williams, Yongquan Xue

Abstract: The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve wit… ▽ More The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve with redshift across z = 0-2 following L_HMXB/SFR ~ 1 + z and L_LMXB/M_star ~ (1 + z)^{2-3}. However, these measurements alone do not directly reveal the physical mechanisms behind the redshift evolution of X-ray binaries (XRBs). We derive star-formation histories for a sample of 344 normal galaxies in the CDFs, using spectral energy distribution (SED) fitting of FUV-to-FIR photometric data, and construct a self-consistent, age-dependent model of the X-ray emission from the galaxies. Our model quantifies how X-ray emission from hot gas and XRB populations vary as functions of host stellar-population age. We find that (1) the ratio L_X/M_star declines by a factor of ~1000 from 0-10 Gyr and (2) the X-ray SED becomes harder with increasing age, consistent with a scenario in which the hot gas contribution to the X-ray SED declines quickly for ages above 10 Myr. When dividing our sample into subsets based on metallicity, we find some indication that L_X/M_star is elevated for low-metallicity galaxies, consistent with recent studies of X-ray scaling relations. However, additional statistical constraints are required to quantify both the age and metallicity dependence of X-ray emission from star-forming galaxies. △ Less

Submitted 6 December, 2021; originally announced December 2021.

The TREX Survey: Kinematical Complexity Throughout M33's Stellar Disk and Evidence for a Stellar Halo

Authors: Karoline M. Gilbert, Amanda C. N. Quirk, Puragra Guhathakurta, Erik Tollerud, Jennifer Wojno, Julianne J. Dalcanton, Meredith J. Durbin, Anil Seth, Benjamin F. Williams, Justin T. Fung, Pujita Tangirala, Ibrahim Yusufali

Abstract: We present initial results from a large spectroscopic survey of stars throughout M33's stellar disk. We analyze a sample of 1667 red giant branch (RGB) stars extending to projected distances of $\sim 11$ kpc from M33's center ($\sim 18$ kpc, or $\sim 10$ scale lengths, in the plane of the disk). The line-of-sight velocities of RGB stars show the presence of two kinematical components. One componen… ▽ More We present initial results from a large spectroscopic survey of stars throughout M33's stellar disk. We analyze a sample of 1667 red giant branch (RGB) stars extending to projected distances of $\sim 11$ kpc from M33's center ($\sim 18$ kpc, or $\sim 10$ scale lengths, in the plane of the disk). The line-of-sight velocities of RGB stars show the presence of two kinematical components. One component is consistent with rotation in the plane of M33's HI disk and has a velocity dispersion ($\sim 19$ km s$^{-1}$) consistent with that observed in a comparison sample of younger stars, while the second component has a significantly higher velocity dispersion. A two-component fit to the RGB velocity distribution finds that the high dispersion component has a velocity dispersion of $59.3^{+2.6}_{-2.5}$ km s$^{-1}$ and rotates very slowly in the plane of the disk (consistent with no rotation at the $<1.5σ$ level), which favors interpreting it as a stellar halo rather than a thick disk population. A spatial analysis indicates that the fraction of RGB stars in the high-velocity-dispersion component decreases with increasing radius over the range covered by the spectroscopic sample. Our spectroscopic sample establishes that a significant high-velocity-dispersion component is present in M33's RGB population from near M33's center to at least the radius where M33's HI disk begins to warp at 30$'$ ($\sim 7.5$ kpc) in the plane of the disk. This is the first detection and spatial characterization of a kinematically hot stellar component throughout M33's inner regions. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: 33 pages, 19 figures. Accepted for publication in ApJ

The Hough Stream Spotter: A New Method for Detecting Linear Structure in Resolved Stars and Application to the Stellar Halo of M31

Authors: Sarah Pearson, Susan E. Clark, Alexis J. Demirjian, Kathryn V. Johnston, Melissa K. Ness, Tjitske K. Starkenburg, Benjamin F. Williams, Rodrigo A. Ibata

Abstract: Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS),… ▽ More Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS), and apply it to the Pan-Andromeda Archaeological Survey (PAndAS) photometric data of resolved stars in M31's stellar halo. We first demonstrate that our code can re-discover known dwarf streams in M31. We then use the HSS to blindly identify 27 linear GC stream-like structures in the PAndAS data. For each HSS GC stream candidate, we investigate the morphologies of the streams and the colors and magnitudes of all stars in the candidate streams. We find that the five most significant detections show a stronger signal along the red giant branch in color-magnitude diagrams (CMDs) than spurious non-stream detections. Lastly, we demonstrate that the HSS will easily detect globular cluster streams in future Nancy Grace Roman Space Telescope data of nearby galaxies. This has the potential to open up a new discovery space for GC stream studies, GC stream gap searches, and for GC stream-based constraints on the nature of dark matter. △ Less

Submitted 3 January, 2022; v1 submitted 30 June, 2021; originally announced July 2021.

Comments: Accepted to ApJ, 34 pages, 18 figures

The Masses of Supernova Remnant Progenitors in NGC 6946

Authors: Brad Koplitz, Jared Johnson, Benjamin F. Williams, Knox S. Long, William P. Blair, Jeremiah W. Murphy, Andrew Dolphin, Tristan Hillis

Abstract: We constrained the progenitor masses for 169 supernova remnants, 8 historically observed supernovae, and the black hole formation candidate in NGC 6946, finding that they are consistent with originating from a standard initial mass function. Additionally, there were 16 remnants that showed no sign of nearby star formation consistent with a core-collapse supernova, making them good Type Ia candidat… ▽ More We constrained the progenitor masses for 169 supernova remnants, 8 historically observed supernovae, and the black hole formation candidate in NGC 6946, finding that they are consistent with originating from a standard initial mass function. Additionally, there were 16 remnants that showed no sign of nearby star formation consistent with a core-collapse supernova, making them good Type Ia candidates. Using $Hubble$ $Space$ $Telescope$ broadband imaging, we measured stellar photometry of ACS/WFC fields in F435W, F555W, F606W, and F814W filters as well as WFC3/UVIS fields in F438W, F606W, and F814W. We then fitted this photometry with stellar evolutionary models to determine the ages of the young populations present at the positions of the SNRs and SNe. We then infer a progenitor mass probability distribution from the fitted age distribution. For 37 SNRs we tested how different filter combinations affected the inferred masses. We find that filters sensitive to H$α$, [N II], and [S II] gas emission can bias mass estimates for remnants that rely on our technique. Using a KS-test analysis on our most reliable measurements, we find the progenitor mass distribution is well-matched by a power-law index of $-2.6^{+0.5}_{-0.6}$, which is consistent with a standard initial mass function. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 28 pages, 13 figures, 4 tables

The Wolf-Rayet + Black Hole Binary NGC 300 X-1: What is the Mass of the Black Hole?

Authors: Breanna A. Binder, Janelle M. Sy, Michael Eracleous, Dimitris M. Christodoulou, Sayantan Bhattacharya, Rigel Cappallo, Silas Laycock, Paul P. Plucinsky, Benjamin F. Williams

Abstract: We present new X-ray and UV observations of the Wolf-Rayet + black hole binary system NGC 300 X-1 with the Chandra X-ray Observatory and the Hubble Space Telescope Cosmic Origins Spectrograph. When combined with archival X-ray observations, our X-ray and UV observations sample the entire binary orbit, providing clues to the system geometry and interaction between the black hole accretion disk and… ▽ More We present new X-ray and UV observations of the Wolf-Rayet + black hole binary system NGC 300 X-1 with the Chandra X-ray Observatory and the Hubble Space Telescope Cosmic Origins Spectrograph. When combined with archival X-ray observations, our X-ray and UV observations sample the entire binary orbit, providing clues to the system geometry and interaction between the black hole accretion disk and the donor star wind. We measure a binary orbital period of 32.7921$\pm$0.0003 hr, in agreement with previous studies, and perform phase-resolved spectroscopy using the X-ray data. The X-ray light curve reveals a deep eclipse, consistent with inclination angles of $i=60-75^{\circ}$, and a pre-eclipse excess consistent with an accretion stream impacting the disk edge. We further measure radial velocity variations for several prominent FUV spectral lines, most notably He II $λ$1640 and C IV $λ$1550. We find that the He II emission lines systematically lag the expected Wolf-Rayet star orbital motion by a phase difference $Δφ\sim0.3$, while C IV $λ$1550 matches the phase of the anticipated radial velocity curve of the Wolf-Rayet donor. We assume the C IV $λ$1550 emission line follows a sinusoidal radial velocity curve (semi-amplitude = 250 km s$^{-1}$) and infer a BH mass of 17$\pm$4 M$_{\odot}$. Our observations are consistent with the presence of a wind-Roche lobe overflow accretion disk, where an accretion stream forms from gravitationally focused wind material and impacts the edge of the black hole accretion disk. △ Less

Submitted 13 February, 2021; originally announced February 2021.

Comments: Accepted to ApJ. 22 pages, 16 figures, 7 tables

Progenitor Mass Distribution for 22 Historic Core-Collapse Supernovae

Authors: Mariangelly Díaz-Rodríguez, Jeremiah W. Murphy, Benjamin F. Williams, Julianne J. Dalcanton, Andrew E. Dolphin

Abstract: We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are… ▽ More We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_\textrm{min}}~=~8.60^{+0.37}_{-0.41} \ M_\odot$ and the slope of the progenitor mass distribution is $α= -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($α= -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses. △ Less

Submitted 22 June, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) I. Ultraviolet to Infrared Photometry of 22 Million Stars in M33

Authors: Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Dustin Lang, Leo Girardi, Adam Smercina, Andrew Dolphin, Daniel R. Weisz, Yumi Choi, Eric F. Bell, Erik Rosolowsky, Evan Skillman, Eric W. Koch, Christine W. Lindberg, Lea Hagen, Karl D. Gordon, Anil Seth, Karoline Gilbert, Puragra Guhathakurta, Tod Lauer, Luciana Bianchi

Abstract: We present panchromatic resolved stellar photometry for 22 million stars in the Local Group dwarf spiral Triangulum (M33), derived from Hubble Space Telescope (HST) observations with the Advanced Camera for Surveys (ACS) in the optical (F475W, F814W), and the Wide Field Camera 3 (WFC3) in the near ultraviolet (F275W, F336W) and near-infrared (F110W, F160W) bands. The large, contiguous survey area… ▽ More We present panchromatic resolved stellar photometry for 22 million stars in the Local Group dwarf spiral Triangulum (M33), derived from Hubble Space Telescope (HST) observations with the Advanced Camera for Surveys (ACS) in the optical (F475W, F814W), and the Wide Field Camera 3 (WFC3) in the near ultraviolet (F275W, F336W) and near-infrared (F110W, F160W) bands. The large, contiguous survey area covers $\sim$14 square kpc and extends to 3.5 kpc (14 arcmin, or 1.5-2 scale lengths) from the center of M33. The PHATTER observing strategy and photometry technique closely mimic those of the Panchromatic Hubble Andromeda Treasury (PHAT), but with updated photometry techniques that take full advantage of all overlapping pointings (aligned to within $<$5-10 milliarcseconds) and improved treatment of spatially-varying point spread functions. The photometry reaches a completeness-limited depth of F475W$\sim$28.5 in the lowest surface density regions observed in M33 and F475W$\sim$26.5 in the most crowded regions found near the center of M33. We find the young populations trace several relatively tight arms, while the old populations show a clear, looser two-armed structure. We present extensive analysis of the data quality including artificial star tests to quantify completeness, photometric uncertainties, and flux biases. This stellar catalog is the largest ever produced for M33, and is publicly available for download by the community. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 38 pages, 6 tables, 25 figures, accepted for publication in ApJS

Three-Dimensional Structure and Dust Extinction in the Small Magellanic Cloud

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

Abstract: We examine the three-dimensional structure and dust extinction properties in a ~ 200 pc $\times$ 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). We model a deep Hubble Space Telescope optical color-magnitude diagram (CMD) of red clump and red giant branch stars to infer the dust extinction and galactic structure. We model the distance distribution of the stellar component w… ▽ More We examine the three-dimensional structure and dust extinction properties in a ~ 200 pc $\times$ 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). We model a deep Hubble Space Telescope optical color-magnitude diagram (CMD) of red clump and red giant branch stars to infer the dust extinction and galactic structure. We model the distance distribution of the stellar component with a Gaussian and find a centroid distance of 65.2 kpc (distance modulus $μ$ = 19.07 mag) with a FWHM $\approx$ 11.3 kpc. This large extent along the line of sight reproduces results from previous studies using variable stars and red clump stars. Additionally, we find an offset between the stellar and dust distributions, with the dust on the near side relative to the stars by 3.22 $^{+1.69}_{-1.44}$ kpc, resulting in a 73% reddened fraction of stars. Modeling the dust layer with a log-normal $A_V$ distribution indicates a mean extinction $\langle A_V \rangle$ = 0.41 $\pm$ 0.09 mag. We also calculate $A_V/N_H$ = 3.2 - 4.2 $\times10^{-23}$ mag cm$^2$ H$^{-1}$ which is significantly lower than the Milky Way value but is comparable to previous SMC dust-to-gas ratio measurements. Our results yield the first joint dust extinction and 3D geometry properties in a key region in the SMC. This study demonstrates that CMD modeling can be a powerful tool to simultaneously constrain dust extinction and geometry properties in nearby galaxies. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: 16 pages, 10 figures. Accepted for publication in ApJ

Multiwavelength Characterization of the High Mass X-ray Binary Population of M31

Authors: Margaret Lazzarini, Benjamin F. Williams, Meredith Durbin, Julianne Dalcanton, Vallia Antoniou, Breanna A. Binder, Michael Eracleous, Paul P. Plucinsky, Manami Sasaki, Neven Vulic

Abstract: We present our analysis of high quality high mass X-ray binary (HMXB) candidates in M31 selected from point-source optical-counterpart candidates from the Chandra-PHAT survey catalog. We fit the spectral energy distributions (SEDs) of optical counterpart candidates using the Bayesian Extinction and Stellar Tool (BEAST). We used the best-fit luminosity, effective temperature, radius and dust redden… ▽ More We present our analysis of high quality high mass X-ray binary (HMXB) candidates in M31 selected from point-source optical-counterpart candidates from the Chandra-PHAT survey catalog. We fit the spectral energy distributions (SEDs) of optical counterpart candidates using the Bayesian Extinction and Stellar Tool (BEAST). We used the best-fit luminosity, effective temperature, radius and dust reddening for the companion stars in combination with the local star formation history, dust maps of M31, published X-ray spectral fits from XMM-Newton observations, IR colors, and Chandra X-ray hardness ratios to determine our best sample of HMXB candidates. The age distribution of the HMXB sample appears peaked between 10 and 50 Myr, consistent with findings in other nearby galaxies. Using the age distribution and mean SFR, we find that 80$-$136 HMXBs were produced per unit of star formation rate over the last 50 Myr and 89$-$163 HMXBs were produced per unit of star formation rate over the last 80 Myr, if we expand the assumed age limit beyond the lifetimes of single massive stars. We also calculate the HMXB production rate (HMXBs/M$_{\odot}$) over time, which ranges from $7 \times 10^{-7}$ to $4 \times 10^{-6}$ HMXBs/M$_{\odot}$ over the last 80 Myr, in agreement with both theoretical predictions and measured production rates in other galaxies. △ Less

Submitted 19 November, 2020; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: accepted for publication in ApJ, 24 pages, 5 figures, 4 tables

Mapping the Escape Fraction of Ionizing Photons Using Resolved Stars: A Much Higher Escape Fraction for NGC 4214

Authors: Yumi Choi, Julianne J. Dalcanton, Benjamin F. Williams, Evan D. Skillman, Morgan Fouesneau, Karl D. Gordon, Karin M. Sandstrom, Daniel R. Weisz, Karoline M. Gilbert

Abstract: We demonstrate a new method for measuring the escape fraction of ionizing photons using Hubble Space Telescope imaging of resolved stars in NGC 4214, a local analog of high-redshift starburst galaxies that are thought to be responsible for cosmic reionization. Specifically, we forward model the UV through near-IR spectral energy distributions of $\sim$83,000 resolved stars to infer their individua… ▽ More We demonstrate a new method for measuring the escape fraction of ionizing photons using Hubble Space Telescope imaging of resolved stars in NGC 4214, a local analog of high-redshift starburst galaxies that are thought to be responsible for cosmic reionization. Specifically, we forward model the UV through near-IR spectral energy distributions of $\sim$83,000 resolved stars to infer their individual ionizing flux outputs. We constrain the local escape fraction by comparing the number of ionizing photons produced by stars to the number that are either absorbed by dust or consumed by ionizing the surrounding neutral hydrogen in individual star-forming regions. We find substantial spatial variation in the escape fraction (0-40%). Integrating over the entire galaxy yields a global escape fraction of 25% (+16%/-15%). This value is much higher than previous escape fractions of zero reported for this galaxy. We discuss sources of this apparent tension, and demonstrate that the viewing angle and the 3D ISM geometric effects are the cause. If we assume the NGC 4214 has no internal dust, like many high-redshift galaxies, we find an escape fraction of 59% (an upper-limit for NGC 4214). This is the first non-zero escape fraction measurement for UV-faint (M$_{\rm FUV}$ = -15.9) galaxies at any redshift, and supports the idea that starburst UV-faint dwarf galaxies can provide a sufficient amount of ionizing photons to the intergalactic medium. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 29 pages, 18 figures, 3 tables, accepted for publication in ApJ

PHAT XX. AGB stars and other cool giants in M31 star clusters

Authors: Leo Girardi, Martha L. Boyer, L. Clifton Johnson, Julianne J. Dalcanton, Philip Rosenfield, Anil C. Seth, Evan D. Skillman, Daniel R. Weisz, Benjamin F. Williams, Antara Raaghavi Bhattacharya, Alessandro Bressan, Nelson Caldwell, Yang Chen, Andrew E. Dolphin, Morgan Fouesneau, Steven Goldman, Puragra Guhathakurta, Paola Marigo, Sagnick Mukherjee, Giada Pastorelli, Amanda Quirk, Monika Soraisam, Michele Trabucchi

Abstract: The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of… ▽ More The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of the thermally-pulsing AGB stars as well as early-AGB stars and other luminous cool giants expected in young stellar populations (e.g. massive red supergiants, and intermediate-mass red helium-burning stars). The AGB stars can be differentiated, a posteriori, using the ages already estimated for our cluster sample. 937 candidates are found within the cluster aperture radii, half (450) of which are very likely cluster members. Cross-matching with additional databases reveals two carbon stars and ten secure variables among them. The field-corrected age distribution reveals the presence of young supergiants peaking at ages smaller than 100 Myr, followed by a long tail of AGB stars extending up to the oldest possible ages. This long tail reveals the general decrease in the numbers of AGB stars from initial values of 50e-6/Msun at 100 Myr down to 5e-6/Msun at 10 Gyr. Theoretical models of near-solar metallicity reproduce this general trend, although with localized discrepancies over some age intervals, whose origin is not yet identified. The entire catalogue is released together with finding charts to facilitate follow-up studies. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: To appear in ApJ. Data will be soon on MAST

MCR-TRGB: A Multiwavelength-Covariant, Robust Tip of the Red Giant Branch Measurement Method

Authors: M. J. Durbin, R. L. Beaton, J. J. Dalcanton, B. F. Williams, M. L. Boyer

Abstract: We present a new method to measure colors and magnitudes of the tip of the red giant branch in multiple bandpasses simultaneously by fitting an n-dimensional Gaussian to photometry of candidate tip stars. We demonstrate that this method has several advantages over traditional edge detection, particularly in regimes where the TRGB magnitude is strongly color-dependent, as is the case in the near-in… ▽ More We present a new method to measure colors and magnitudes of the tip of the red giant branch in multiple bandpasses simultaneously by fitting an n-dimensional Gaussian to photometry of candidate tip stars. We demonstrate that this method has several advantages over traditional edge detection, particularly in regimes where the TRGB magnitude is strongly color-dependent, as is the case in the near-infrared. We apply this method to a re-reduction of a set of optical and near-IR HST data originally presented in Dalcanton et al. (2012). The re-reduction takes advantage of the increased depth and accuracy in the NIR photometry enabled by simultaneous reduction with higher resolution optical data in crowded fields (Williams et al. 2014). We compare three possible absolute calibrations of the resulting apparent TRGB measurements, one adopting the same distance moduli as in Dalcanton et al. (2012), and two based on predicted TRGB absolute magnitudes from two widely-used, modern sets of model isochrones. We find systematic offsets among the model absolute calibrations at the ~0.1 mag level, in line with previous investigations. The models also have difficulty reproducing the optical-NIR color-magnitude behavior of our measurements, making these observations a useful benchmark for future improvements. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 32 pages, 21 figures, accepted to ApJ

Mass-to-Light Ratios of Spatially Resolved Stellar Populations in M31

Authors: O. Grace Telford, Julianne J. Dalcanton, Benjamin F. Williams, Eric F. Bell, Andrew E. Dolphin, Meredith J. Durbin, Yumi Choi

Abstract: A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial… ▽ More A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial variations in $M_\star/L$ and their dependence on color, star formation history, and dust across the disk of M31, using a map of $M^\mathrm{CMD}_\star$ derived from color-magnitude diagrams of resolved stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. First, we find comparable scatter in $M_\star/L$ for the optical and mid-IR, contrary to the common idea that $M_\star/L$ is less variable in the IR. Second, we confirm that $M_\star/L$ is correlated with color for both the optical and mid-IR and report color vs. $M_\star/L$ relations (CMLRs) in M31 for filters used in the Sloan Digital Sky Survey (SDSS) and Widefield Infrared Survey Explorer (WISE). Third, we show that the CMLR residuals correlate with recent SFH, such that quiescent regions are offset to higher $M_\star/L$ than star-forming regions at a fixed color. The mid-IR CMLR, however, is not linear due to the high scatter of $M_\star/L$ in star-forming regions. Finally, we find a flatter optical CMLR than any SPS-based CMLRs in the literature. We show this is an effect of dust geometry, which is typically neglected but should be accounted for when using optical data to map $M_\star/L$. △ Less

Submitted 24 January, 2020; originally announced January 2020.

Comments: Accepted to ApJ. Main text is 25 pages, including 15 figures

The Hubble Space Telescope Advanced Camera for Surveys Emission Line Survey of Andromeda. I: Classical Be Stars

Authors: M. Peters, J. P. Wisniewski, B. F. Williams, J. R. Lomax, Y. Choi, M. Durbin, L. C. Johnson, A. R. Lewis, J. Lutz, T. A. A. Sigut, A. Wallach, J. J. Dalcanton

Abstract: We present results from a 2-epoch HST H$α$ emission line survey of the Andromeda Galaxy that overlaps the footprint of the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We find 552 (542) classical Be stars and 8429 (8556) normal B-type stars in epoch # 1 (epoch # 2), yielding an overall fractional Be content of 6.15% $\pm$0.26% (5.96% $\pm$0.25%). The fractional Be content decreased with s… ▽ More We present results from a 2-epoch HST H$α$ emission line survey of the Andromeda Galaxy that overlaps the footprint of the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We find 552 (542) classical Be stars and 8429 (8556) normal B-type stars in epoch # 1 (epoch # 2), yielding an overall fractional Be content of 6.15% $\pm$0.26% (5.96% $\pm$0.25%). The fractional Be content decreased with spectral sub-type from $\sim$23.6% $\pm$2.0% ($\sim$23.9% $\pm$2.0%) for B0-type stars to $\sim$3.1% $\pm$0.34% ($\sim$3.4% $\pm$0.35%) for B8-type stars in epoch # 1 (epoch # 2). We observe a clear population of cluster Be stars at early fractional main sequence lifetimes, indicating that a subset of Be stars emerge onto the ZAMS as rapid rotators. Be stars are 2.8x rarer in M31 for the earliest sub-types compared to the SMC, confirming that the fractional Be content decreases in significantly more metal rich environments (like the Milky Way and M31). However, M31 does not follow a clear trend of Be fraction decreasing with metallicity compared to the Milky Way, which may reflect that the Be phenomenon is enhanced with evolutionary age. The rate of disk-loss or disk-regeneration episodes we observe, 22% $\pm$ 2% yr$^{-1}$, is similar to that observed for seven other Galactic clusters reported in the literature, assuming these latter transient fractions scale by a linear rate. The similar number of disk-loss events (57) as disk-renewal events (43) was unexpected since disk dissipation time-scales can be $\sim$2x the typical time-scales for disk build-up phases. △ Less

Submitted 9 January, 2020; v1 submitted 11 December, 2019; originally announced December 2019.

Comments: Accepted to AJ