Search

We present $\sim300$ stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI ($M_V = -7.5$) and Andromeda XXVIII ($M_V = -8.8$) derived using metallicity-sensitive Calcium H & K narrow-band Hubble Space Telescope imaging. These are the first individual stellar metallicities in And~XVI (95 stars). Our And~XXVIII sample (191 stars) is a factor of $\sim15$ increase over literature metallicities. For And~XVI, we measure $\langle \mbox{[Fe/H]}\rangle = -2.17^{+0.05}_{-0.05}$, $\sigma_{\mbox{[Fe/H]}}=0.33^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.23\pm0.15$ dex $R_e^{-1}$. We find that And XVI is more metal-rich than MW UFDs of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure $\langle \mbox{[Fe/H]}\rangle = -1.95^{+0.04}_{-0.04}$, $\sigma_{\mbox{[Fe/H]}}=0.34^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.46 \pm 0.10$~dex~$R_e^{-1}$, placing it on the dwarf galaxy mass-metallicity relation. Neither galaxy has a metallicity distribution function with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint ($L \lesssim 10^6~L_{\odot}$) galaxies outside the Milky Way halo. Both galaxies' gradients are consistent with predictions from the FIRE simulations, where an age-gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal poor ($\mbox{[Fe/H]} < -3.0$) star candidates, which make up 7% of And~XVI's MDF and 6% of And~XXVIII's., Comment: 17 pages, 5 figures, 4 tables, ApJ submitted; comments welcome!

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 $\Gamma$, 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{\Gamma} = 1.49\pm0.18$, with an intrinsic scatter of $\sigma^{2}_{\Gamma} = 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{\Gamma}$ 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 ($\Gamma = 1.30$) IMF, as well as the Salpeter ($\Gamma = 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., Comment: Accepted for publication in ApJ. 9 Figures, 1 Table

The tip of the red giant (TRGB) is a standardizable candle and is identifiable as the discontinuity at the bright extreme of the red giant branch (RGB) stars in color-magnitude diagram (CMD) space. The TRGB-based distance method has been calibrated and used to measured distances to galaxies out to $D\leq20$ Mpc with the $I$-band equivalent Hubble Space Telescope ($HST$) $F814W$ filter, and as an important rung in the distance ladder to measure the Hubble constant, $H_0$. In the infrared (IR), the TRGB apparent magnitude ranges from $1-2$ magnitudes brighter than in the optical, and now with the IR James Webb Space Telescope ($JWST$) observatory the feasible distance range of the TRGB method can be extended to $\sim50$ Mpc. However, in the IR the TRGB luminosity depends to varying degrees on stellar metallicity and age. In this study we standardize the TRGB luminosity using stellar colors as a proxy for metallicity/age to derive color-based corrections for the $JWST$ Near-Infrared Camera (NIRCam) short wavelength (SW) filters $F090W$, $F115W$, $F150W$ and the long wavelength (LW) filters $F277W$, $F356W,$ and $F444W$. We provide recommended filter combinations for distance measurements depending on the requisite precision. For science requiring high precision ($\leq1\%$ in distance) and robustness we recommend measuring the TRGB in $F090W$ vs $F090W-F150W$ or $F115W$ vs. $F115W-F277W$ with the caveat that even with $JWST$ long integration times will be necessary at further distances. If lower precision ($>1.5\%$ in distance) can be tolerated, or if shorter integration times are desirable, we recommend measuring the TRGB in either $F115W$ or $F150W$ paired with $F356W$. We do not recommend $F444W$ for precision TRGB measurements due to its lower angular resolution., Comment: 25 pages, 11 figures

The Tip of the Red Giant Branch (TRGB)-based distance method in the I band is one of the most efficient and precise techniques for measuring distances to nearby galaxies (D <= 15 Mpc). The TRGB in the near infrared (NIR) is 1 to 2 magnitudes brighter relative to the I band, and has the potential to expand the range over which distance measurements to nearby galaxies are feasible. Using Hubble Space Telescope (HST) imaging of 12 fields in 8 nearby galaxies, we determine color-based corrections and zero points of the TRGB in the Wide Field Camera 3 IR (WFC3/IR) F110W and F160W filters. First, we measure TRGB distances in the I band equivalent Advanced Camera System (ACS) F814W filter from resolved stellar populations with the HST. The TRGB in the ACS F814W filter is used for our distance anchor and to place the WFC3/IR magnitudes on an absolute scale. We then determine the color dependence (a proxy for metallicity/age) and zero point of the NIR TRGB from photometry of WFC3/IR fields which overlap with the ACS fields. The new calibration is accurate to ~1% in distance, relative to the F814W TRGB. Validating the accuracy of the calibrations, we find that the distance modulus for each field using the NIR TRGB calibration agrees with the distance modulus of the same fields as determined from the F814W TRGB. This is a JWST preparatory program and the work done here will directly inform our approach to calibrating the TRGB in JWST NIRCam and NIRISS photometric filters., Comment: 23 pages, 8 figures, accepted by the Astrophysical Journal (ApJ)

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., Comment: 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

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., Comment: 15 page, 12 figures, submitted to AAS Journals

The tip of the red giant branch (TRGB) allows for the measurement of precise and accurate distances to nearby galaxies, based on the brightest ascent of low-mass red giant branch stars before they undergo the helium flash. With the advent of JWST, there is great promise to utilize the technique to measure galaxy distances out to at least 50 Mpc, significantly further than HST's reach of 20 Mpc. However, with any standard candle, it is first necessary to provide an absolute reference. Here we use Cycle 1 data to provide an absolute calibration in the F090W filter. F090W is most similar to the F814W filter commonly used for TRGB measurements with HST, which had been adopted by the community due to minimal dependence from the underlying metallicities and ages of stars. The imaging we use was taken in the outskirts of NGC 4258, which has a direct geometrical distance measurement from the Keplerian motion of its water megamaser. Utilizing several measurement techniques, we find $M_{TRGB}^{F090W}$ = -4.362 $\pm$ 0.033 (stat) $\pm$ 0.045 (sys) mag (Vega) for the metal-poor TRGB. We also perform measurements of the TRGB in two Type Ia supernova hosts, NGC 1559, and NGC 5584. We find good agreement between our TRGB distances and previous distance determinations to these galaxies from Cepheids ($\Delta$ = 0.01 $\pm$ 0.06 mag), with these differences being too small to explain the Hubble tension ($\sim$0.17 mag). As a final bonus, we showcase the serendipitous discovery of a faint dwarf galaxy near NGC 5584., Comment: 20 pages, 12 figures, accepted to ApJ

We present high-definition observations with the James Webb Space Telescope of >1000 Cepheids in a geometric anchor of the distance ladder, NGC4258, and in 5 hosts of 8 SNe~Ia, a far greater sample than previous studies with JWST. These galaxies individually contain the largest samples of Cepheids, an average of >150 each, producing the strongest statistical comparison to those previously measured with the Hubble Space Telescope in the NIR. They also span the distance range of those used to determine the Hubble constant with HST, allowing us to search for a distance-dependent bias in HST measurements. The superior resolution of JWST negates crowding noise, the largest source of variance in the NIR Cepheid Period-Luminosity relations (Leavitt laws) measured with HST. Together with the use of two-epochs to constrain Cepheid phases and three filters to remove reddening, we reduce the dispersion in the Cepheid PL relations by a factor of 2.5. We find no significant difference in the mean distance measurements determined from HST and JWST, with a formal difference of -0.01+/-0.03 mag. This result is independent of zeropoints and analysis variants including metallicity dependence, local crowding, choice of filters, and relation slope. We can reject the hypothesis of unrecognized crowding of Cepheid photometry from HST that grows with distance as the cause of the ``Hubble Tension'' at 8.2 sigma, i.e., greater confidence than that of the Hubble Tension itself. We conclude that errors in photometric measurements of Cepheids across the distance ladder do not significantly contribute to the Tension., Comment: ApJ Letters, accepted

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., Comment: 12 pages, 10 figures, 1 tables, accepted for publication at MNRAS

We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope (HST) narrow-band (F395N) Calcium H & K (CaHK) imaging. Our sample is a factor of $\sim7$ larger than what is published. Our main findings are: (i) A global metallicity distribution function (MDF) with $\langle \mbox{[Fe/H]} \rangle = -1.55 \pm 0.04$ and $\sigma_{\mbox{[Fe/H]}}=0.54\pm0.03$; (ii) A metallicity gradient of $-0.54 \pm 0.07$ dex $R_e^{-1}$ ($-2.1 \pm 0.3$ dex kpc$^{-1}$) over the extent of our imaging ($\sim 2.5 R_e$), which is steeper than literature measurements. Our finding is consistent with predicted gradients from the publicly-available FIRE-2 simulations, in which bursty star formation creates stellar population gradients and dark matter cores; (iii) Tucana's bifurcated RGB has distinct metallicities: a blue RGB with $\langle \mbox{[Fe/H]} \rangle = -1.78 \pm 0.06$ and $\sigma_{\mbox{[Fe/H]}}=0.44^{+0.07}_{-0.06}$, and a red RGB with $\langle \mbox{[Fe/H]} \rangle = -1.08 \pm 0.07$ and $\sigma_{\mbox{[Fe/H]}}=0.42 \pm 0.06$; (iv) At fixed stellar mass, Tucana is more MR than MW satellites by $\sim 0.4$ dex, but its blue RGB is chemically comparable to the satellites. Tucana's MDF appears consistent with star-forming isolated dwarfs, though MDFs of the latter are not as well-populated; (v) $\sim2$% of Tucana's stars have $\mbox{[Fe/H]} < -3$ and 20% $\mbox{[Fe/H]} > -1$. We provide a catalog for community spectroscopic follow-up., Comment: Replaced with ApJ published version; 23 pages, 18 figures

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., Comment: 27 pages, 10 figures, 1 table

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., Comment: Accepted to AJ

Understanding the interplay of stellar feedback and turbulence in the interstellar medium (ISM) is essential to modeling the evolution of galaxies. To determine the timescales over which stellar feedback drives turbulence in the ISM, we performed a spatially resolved, multi-wavelength study of the nearby star-forming dwarf galaxy UGC 4305 (aka Holmberg II). As indicators of turbulence on local scales (400 pc), we utilized ionized gas velocity dispersion derived from IFU H$\alpha$ observations and atomic gas velocity dispersion and energy surface densities derived from HI synthesis observations with the Very Large Array. These indicators of turbulence were tested against star formation histories over the past 560 Myr derived from Color-Magnitude Diagrams (CMD) using Spearman's rank correlation coefficient. The strongest correlation identified at the 400 pc scale is between measures of HI turbulence and star formation 70-140 Myr ago. We repeated our analysis of UGC 4305's current turbulence and past star formation activity on multiple physical scales ($\sim$560, and 800 pc) to determine if there are indications of changes in the correlation timescale with changes to the physical scale. No notable correlations were found at larger physical scales emphasizing the importance of analyzing star formation driven turbulence as a local phenomenon., Comment: 17 pages, 9 figure, accepted to AJ

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., Comment: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

High-resolution JWST observations can test confusion-limited HST observations for a photometric bias that could affect extragalactic Cepheids and the determination of the Hubble constant. We present JWST NIRCAM observations in two epochs and three filters of >330 Cepheids in NGC4258 (which has a 1.5% maser-based geometric distance) and in NGC5584 (host of SNIa 2007af), near the median distance of the SH0ES HST SNIa host sample and with the best leverage among them to detect such a bias. JWST provides far superior source separation from line-of-sight companions than HST in the NIR to largely negate confusion or crowding noise at these wavelengths, where extinction is minimal. The result is a remarkable >2.5x reduction in the dispersion of the Cepheid P-L relations, from 0.45 to 0.17 mag, improving individual Cepheid precision from 20% to 7%. Two-epoch photometry confirmed identifications, tested JWST photometric stability, and constrained Cepheid phases. The P-L relation intercepts are in very good agreement, with differences (JWST-HST) of 0.00+/-0.03 and 0.02+/-0.03 mag for NGC4258 and NGC5584, respectively. The difference in the determination of H_0 between HST and JWST from these intercepts is 0.02+/-0.04 mag, insensitive to JWST zeropoints or count-rate non-linearity thanks to error cancellation between rungs. We explore a broad range of analysis variants (including passband combinations, phase corrections, measured detector offsets, and crowding levels) indicating robust baseline results. These observations provide the strongest evidence yet that systematic errors in HST Cepheid photometry do not play a significant role in the present Hubble Tension. Upcoming JWST observations of >12 SNIa hosts should further refine the local measurement of the Hubble constant., Comment: submitted to ApJ, comments welcome

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., Comment: 28 pages, 4 tables, 18 figures, accepted for publication in The Astrophysical Journal Supplement Series

We report the discovery of two ultra-faint dwarf galaxies, Leo M and Leo K, that lie outside the halo of the Milky Way. Using Hubble Space Telescope imaging of the resolved stars, we create color-magnitude diagrams reaching the old main sequence turn-off of each system and (i) fit for structural parameters of the galaxies; (ii) measure their distances using the luminosity of the Horizontal Branch stars; (iii) estimate integrated magnitudes and stellar masses; and (iv) reconstruct the star formation histories. Based on their location in the Local Group, neither galaxy is currently a satellite of the Milky Way, although Leo K is located ~26 kpc from the low-mass galaxy Leo T and these two systems may have had a past interaction. Leo M and Leo K have stellar masses of 1.8 (+0.3/-0.2) x 10^4 Msun and 1.2+/-0.2 x 10^4 Msun, and were quenched 10.6 (+2.2/-1.1) Gyr and 12.8 (+0.1/-4.2) Gyr ago, respectively. Given that the galaxies are not satellites of the MW, it is unlikely that they were quenched by environmental processing. Instead, given their low stellar masses, their early quenching timescales are consistent with the scenario that a combination of reionization and stellar feedback shut-down star formation at early cosmic times., Comment: 12 pages, 9 figures, 1 table

The nearby face-on star forming spiral galaxy NGC 6946 is known as the Fireworks Galaxy due to its hosting an unusually large number of supernova. We analyze its resolved near-ultraviolet (NUV) stellar photometry measured from images taken with the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) with F275W and F336W filters. We model the color-magnitude diagrams (CMD) of the UV photometry to derive the spatially-resolved star formation history (SFH) of NGC 6946 over the last 25 Myr. From this analysis, we produce maps of the spatial distribution of young stellar populations and measure the total recent star formation rate (SFR) of nearly the entire young stellar disk. We find the global SFR(age$\leq$25 Myr)=$13.17 \substack{+0.91 \\-0.79} M_\odot/\rm yr$. Over this period, the SFR is initially very high ($23.39\substack{+2.43\\-2.11} M_\odot/\rm yr$ between 16-25 Myr ago), then monotonically decreases to a recent SFR of $5.31\substack{+0.19\\-0.17} M_\odot/\rm yr$ in the last 10 Myr. This decrease in global star formation rate over the last 25 Myr is consistent with measurements made with other SFR indicators. We discuss in detail two of the most active regions of the galaxy, which we find are responsible for 3% and 5% of the total star formation over the past 6.3 Myr., Comment: 19 pages, 11 figures, accepted for publication in ApJ

We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the number of metallicities in all known MW UFDs. We provide the first well-populated MDFs for all galaxies in this sample, with [Fe/H] ranging from -3.0 dex to -2.0 dex, and $\sigma_{\mbox{[Fe/H]}}$ ranging from 0.3 dex to 0.7 dex. We find a nearly constant [Fe/H] ~ -2.6 over 3 decades in luminosity (~$10^2 - 10^5 L_{\odot}$), suggesting that the mass-metallicity relationship does not hold for such faint systems. We find a larger fraction (24%) of EMP ([Fe/H]<-3.0) stars across our sample compared to the literature (14%), but note that uncertainties in our most metal-poor measurements make this an upper limit. We find 19% of stars in our UFD sample to be metal-rich ([Fe/H]>-2.0), consistent with the sum of literature spectroscopic studies. MW UFDs are known to be predominantly >13 Gyr old, meaning that all stars in our sample are truly ancient, unlike metal-poor stars in the MW, which have a range of possible ages. Our UFD metallicities are not well-matched to known streams in the MW, providing further evidence that known MW substructures are not related to UFDs. We include a catalog of our stars to encourage community follow-up studies, including priority targets for ELT-era observations., Comment: 37 pages, 20 figures, ApJ accepted and published

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., Comment: 22 pages. Accepted to the Astrophysical Journal

Using color-magnitude diagrams from deep archival Hubble Space Telescope imaging, we self-consistently measure the star formation history of Eridanus II (Eri II), the lowest-mass galaxy ($M_{\star}(z=0) \sim 10^5 M_{\odot}$) known to host a globular cluster (GC), and the age, mass, and metallicity of its GC. The GC ($\sim13.2\pm0.4$ Gyr, $\langle$[Fe/H]$\rangle = -2.75\pm0.2$ dex) and field (mean age $\sim13.5\pm0.3$ Gyr, $\langle$[Fe/H]$\rangle = -2.6\pm0.15$ dex) have similar ages and metallicities. Both are reionization-era relics that formed before the peak of cosmic star and GC formation ($z\sim2-4$). The ancient star formation properties of Eri II are not extreme and appear similar to $z=0$ dwarf galaxies. We find that the GC was $\lesssim4$ times more massive at birth than today and was $\sim$10% of the galaxy's stellar mass at birth. At formation, we estimate that the progenitor of Eri II and its GC had $M_{\rm UV} \sim -7$ to $-12$, making it one of the most common type of galaxy in the early Universe, though it is fainter than direct detection limits, absent gravitational lensing. Archaeological studies of GCs in nearby low-mass galaxies may be the only way to constrain GC formation in such low-mass systems. We discuss the strengths and limitations in comparing archaeological and high redshift studies of cluster formation, including challenges stemming from the Hubble Tension, which introduces uncertainties into the mapping between age and redshift., Comment: 11 pages, 3 Figures, 1 Table. Accepted for publication in ApJ

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., Comment: 20 pages, 7 figures, 2 tables, Accepted at ApJ

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., Comment: 32 pages, 18 figures, 4 tables

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., Comment: 5 pages, 1 figure

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., Comment: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

We report the discovery of an ultrafaint dwarf (UFD) galaxy, Pegasus W, located on the far side of the Milky Way-M31 system and outside the virial radius of M31. The distance to the galaxy is 915 (+60/-91) kpc, measured using the luminosity of horizontal branch (HB) stars identified in Hubble Space Telescope optical imaging. The galaxy has a half-light radius (r_h) of 100 (+11/-13) pc, M_V = -7.20 (+0.17/-0.16) mag, and a present-day stellar mass of 6.5 (+1.1/-1.4) x 10^4 Msun. We identify sources in the color-magnitude diagram (CMD) that may be younger than ~500 Myr suggesting late-time star formation in the UFD galaxy, although further study is needed to confirm these are bona fide young stars in the galaxy. Based on fitting the CMD with stellar evolution libraries, Pegasus W shows an extended star formation history (SFH). Using the tau_90 metric (defined as the timescale by which the galaxy formed 90% of its stellar mass), the galaxy was quenched only 7.4 (+2.2/-2.6) Gyr ago, which is similar to the quenching timescale of a number of UFD satellites of M31 but significantly more recent than the UFD satellites of the Milky Way. Such late-time quenching is inconsistent with the more rapid timescale expected by reionization and suggests that, while not currently a satellite of M31, Pegasus W was nonetheless slowly quenched by environmental processes., Comment: 15 pages, 10 figures, 2 tables

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 $\tau_{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., Comment: 8 pages, 4 figures, accepted for publication in MNRAS, Published by Oxford University Press on behalf of the Royal Astronomical Society

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., Comment: Accepted for publication in RNAAS

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., Comment: accepted for publication in ApJ, 37 pages, 18 figures, 6 tables

We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system., Comment: 32 pages, 15 figures, accepted for publication on ApJ

Stellar feedback is fundamental to the modeling of galaxy evolution as it drives turbulence and outflows in galaxies. Understanding the timescales involved are critical for constraining the impact of stellar feedback on the interstellar medium (ISM). We analyzed the resolved star formation histories along with the spatial distribution and kinematics of the atomic and ionized gas of four nearby star-forming dwarf galaxies (NGC 4068, NGC 4163, NGC 6789, UGC 9128) to determine the timescales over which stellar feedback drives turbulence. The four galaxies are within 5 Mpc and have a range of properties including current star formation rates of 0.0005 to 0.01 M$_{\odot}$ yr$^{-1}$, log(M$_*$/M$_{\odot}$) between 7.2 and 8.2, and log(M$_{HI}$/M$_\odot$) between 7.2 and 8.3. Their Color-Magnitude Diagram (CMD) derived star formation histories over the past 500 Myrs were compared to their atomic and ionized gas velocity dispersion and HI energy surface densities as indicators of turbulence. The Spearman's rank correlation coefficient was used to identify any correlations between their current turbulence and their past star formation activity on local scales ($\sim$400 pc). The strongest correlation found was between the HI turbulence measures and the star formation rate 100-200 Myrs ago. This suggests a coupling between the star formation activity and atomic gas on this timescale. No strong correlation between the ionized gas velocity dispersion and the star formation activity between 5-500 Myrs ago was found. The sample and analysis are the foundation of a larger program aimed at understanding the timescales over which stellar feedback drives turbulence., Comment: 21 pages, 19 figures. Accepted for publication in AJ

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 $\alpha = -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., Comment: 18 pages, 15 figures, 1 table, Accepted to ApJ (February 2, 2022)

We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more stars than current spectroscopic MDF determinations. We find that Eri II has a mean metallicity of [Fe/H]$=$-2.50$^{+0.07}_{-0.07}$ and a dispersion of $\sigma_{\mbox{[Fe/H]}}=0.42^{+0.06}_{-0.06}$, which are consistent with spectroscopic MDFs, though more precisely constrained owing to a larger sample. We identify a handful of extremely metal-poor star candidates (EMP; [Fe/H] $< -3$) that are marginally bright enough for spectroscopic follow up. Eri II's MDF appears well-described by a leaky box chemical evolution model. We also compute an updated orbital history for Eri II using Gaia eDR3 proper motions, and find that it is likely on first infall into the Milky Way. Our findings suggest that Eri II underwent an evolutionary history similar to that of an isolated galaxy. Compared to MDFs for select cosmological simulations of similar mass galaxies, we find that Eri II has a lower fraction of stars with [Fe/H] $< -3$, though such comparisons should currently be treated with caution due to a paucity of simulations, selection effects, and known limitations of Ca H&K for EMPs. This study demonstrates the power of deep HST CaHK imaging for measuring the MDFs of UFDs., Comment: Resubmitted to ApJ with revisions following a positive referee report of the initial draft; 29 pages, 19 figures, and 3 tables, including appendix. Comments welcome!

The Survey of HI in Extremely Low-mass Dwarfs (SHIELD) includes a volumetrically complete sample of 82 gas-rich dwarfs with M_HI~<10^7.2 Msun selected from the ALFALFA survey. We are obtaining extensive follow-up observations of the SHIELD galaxies to study their gas, stellar, and chemical content, and to better understand galaxy evolution at the faint end of the HI mass function. Here, we investigate the properties of 30 SHIELD galaxies using Hubble Space Telescope imaging of their resolved stars and Westerbork Synthesis Radio Telescope observations of their neutral hydrogen. We measure tip of the red giant branch (TRGB) distances, star formation activity, and gas properties. The TRGB distances are up to 4x greater than estimates from flow models, highlighting the importance of velocity-independent distance indicators in the nearby universe. The SHIELD galaxies are in under-dense regions, with 23% located in voids; one galaxy appears paired with a more massive dwarf. We quantify galaxy properties at low masses including stellar and HI masses, SFRs, sSFRs, SFEs, birthrate parameters, and gas fractions. The lowest mass systems lie below the mass thresholds where stellar mass assembly is predicted to be impacted by reionization. Even so, we find the star formation properties follow the same trends as higher mass gas-rich systems, albeit with a different normalization. The HI disks are small (<0.7 kpc) making it difficult to measure the HI rotation using standard techniques; we develop a new methodology and report the velocity extent, and its associated spatial extent, with robust uncertainties., Comment: 43 pages, 32 figures, 6 tables

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$\alpha$, [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., Comment: 28 pages, 13 figures, 4 tables

The large and diffuse galaxies NGC1052-DF2 and NGC1052-DF4 have been found to have very low dark matter content and a population of luminous globular clusters. Accurate distance measurements are key to interpreting these observations. Recently, the distance to NGC1052-DF4 was found to be $20.0\pm 1.6$ Mpc by identifying the tip of the red giant branch (TRGB) in 12 orbits of Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) imaging. Here we present 40 orbits of HST ACS data for NGC1052-DF2 and use these data to measure its TRGB. The TRGB is readily apparent in the color-magnitude diagram. Using a forward model that incorporates photometric uncertainties, we find a TRGB magnitude of $m_{\rm F814W, TRGB} = 27.67 \pm 0.10$ mag. The inferred distance is $D_{\rm TRGB} = 22.1 \pm 1.2$ Mpc, consistent with the previous surface brightness fluctuation distances to the bright elliptical galaxy NGC1052. The new HST distance rules out the idea that some of NGC1052-DF2's unusual properties can be explained if it were at $\sim 13$ Mpc; instead, it implies that the galaxy's globular clusters are even more luminous than had been derived using the previous distance of 20 Mpc. The distance from NGC1052-DF2 to NGC1052-DF4 is well-determined at $2.1\pm 0.5$ Mpc, significantly larger than the virial diameter of NGC1052. We discuss the implications for formation scenarios of the galaxies and for the external field effect, which has been invoked to explain the intrinsic dynamics of these objects in the context of modified Newtonian dynamics., Comment: Accepted for publication in ApJ Letters

The Extragalactic Distance Database (EDD) was created as a repository for high quality, redshift-independent distances. A key component of EDD is the Color Magnitude Diagrams/Tip of the Red Giant Branch (CMDs/TRGB) catalog, which provides information on the stellar content of nearby galaxies observed with the Hubble Space Telescope (HST). Here we provide a decadal update to this catalog, which has now doubled in size to over 500 galaxies. We highlight the additions to our data reduction and analysis techniques, and provide examples of the science that has been made possible with this large data set. We find the TRGB to be a reliable measure for distance, and we aim to extend its distance coverage with HST to every galaxy within 10 Mpc. In the near-future, the combination of the James Webb Space Telescope and the Nancy Grace Roman Space Telescope will dramatically increase the number of targets within our grasp., Comment: 8 figures, 20 pages. Accepted to AJ. Materials available at edd.ifa.hawaii.edu

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 $\alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($\alpha = -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.

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., Comment: 38 pages, 6 tables, 25 figures, accepted for publication in ApJS

The Hubble Space Telescope (HST) survey Measuring Young Stars in Space and Time (MYSST) entails some of the deepest photometric observations of extragalactic star formation, capturing even the lowest mass stars of the active star-forming complex N44 in the Large Magellanic Cloud. We employ the new MYSST stellar catalog to identify and characterize the content of young pre-main-sequence (PMS) stars across N44 and analyze the PMS clustering structure. To distinguish PMS stars from more evolved line of sight contaminants, a non-trivial task due to several effects that alter photometry, we utilize a machine learning classification approach. This consists of training a support vector machine (SVM) and a random forest (RF) on a carefully selected subset of the MYSST data and categorize all observed stars as PMS or non-PMS. Combining SVM and RF predictions to retrieve the most robust set of PMS sources, we find $\sim26,700$ candidates with a PMS probability above 95% across N44. Employing a clustering approach based on a nearest neighbor surface density estimate, we identify 18 prominent PMS structures at $1$ $\sigma$ significance above the mean density with sub-clusters persisting up to and beyond $3$ $\sigma$ significance. The most active star-forming center, located at the western edge of N44's bubble, is a subcluster with an effective radius of $\sim 5.6$ pc entailing more than 1,100 PMS candidates. Furthermore, we confirm that almost all identified clusters coincide with known H II regions and are close to or harbor massive young O stars or YSOs previously discovered by MUSE and Spitzer observations., Comment: 29 pages, 21 figures, accepted for publication in AJ

In order to better understand the role of high-mass stellar feedback in regulating star formation in giant molecular clouds, we carried out a Hubble Space Telescope (HST) Treasury Program "Measuring Young Stars in Space and Time" (MYSST) targeting the star-forming complex N44 in the Large Magellanic Cloud (LMC). Using the F555W and F814W broadband filters of both the ACS and WFC3/UVIS, we built a photometric catalog of 461,684 stars down to $m_\mathrm{F555W} \simeq 29$ mag and $m_\mathrm{F814W} \simeq 28$ mag, corresponding to the magnitude of an unreddened 1 Myr pre-main-sequence star of $\approx0.09$ $M_\odot$ at the LMC distance. In this first paper we describe the observing strategy of MYSST, the data reduction procedure, and present the photometric catalog. We identify multiple young stellar populations tracing the gaseous rim of N44's super bubble, together with various contaminants belonging to the LMC field population. We also determine the reddening properties from the slope of the elongated red clump feature by applying the machine learning algorithm RANSAC, and we select a set of Upper Main Sequence (UMS) stars as primary probes to build an extinction map, deriving a relatively modest median extinction $A_{\mathrm{F555W}}\simeq0.77$ mag. The same procedure applied to the red clump provides $A_{\mathrm{F555W}}\simeq 0.68$ mag., Comment: 29 pages, 15 figures, accepted for publication in AJ

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 $\mu$ = 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., Comment: 16 pages, 10 figures. Accepted for publication in ApJ

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., Comment: To appear in ApJ. Data will be soon on MAST

Extremely metal-poor (XMP) galaxies are low-mass, star-forming galaxies with gas-phase oxygen abundances below 12+log(O/H) = 7.35 (~1/20 Zsun). Galaxy evolution scenarios suggest three pathways to form an XMP: (1) secular evolution at low galaxy masses, (2) slow evolution in voids, or (3) dilution of measured abundances from infall of pristine gas. The recently discovered XMP galaxy Leoncino, with an oxygen abundance below 3% Zsun, provides an opportunity to explore these different scenarios. Using Hubble Space Telescope imaging of the resolved stellar populations of Leoncino, we measure the distance to the galaxy to be D=12.1 (+1.7/-3.4) Mpc and find that Leoncino is located in an under-dense environment. Leoncino has a compact morphology, hosts a population of young, massive stars, has a high gas-to-star mass ratio, and shows signs of interaction with a galaxy nearby on the sky, UGC 5186. Similar to nearly all XMP galaxies known in the nearby universe, Leoncino is offset from the Luminosity-Metallicity (LZ) relation. Yet, Leoncino is consistent with the stellar Mass-Metallicity (MZ) relation defined by Local Volume galaxies. Thus, our results suggest that the offset from the LZ relation is due to higher recent star formation, likely triggered by a minor interaction, while the low oxygen abundance is consistent with the expectation that low-mass galaxies will undergo secular evolution marked by inefficient star formation and metal-loss via galactic winds. This is in contrast to XMP galaxies that are outliers in both the LZ and MZ relations; in such cases, the low oxygen abundances are best explained by dilution due to the infall of pristine gas. We also discuss why quiescent XMP galaxies are underrepresented in current surveys., Comment: 20 pages, 8 figures, 1 table

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$., Comment: Accepted to ApJ. Main text is 25 pages, including 15 figures

Previous studies have shown that the large, diffuse galaxies NGC1052-DF2 and NGC1052-DF4 both have populations of unusually luminous globular clusters as well as a very low dark matter content. Here we present newly-obtained deep Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) imaging of one of these galaxies, NGC1052-DF4. We use these data to measure the distance of the galaxy from the location of the tip of the red giant branch (TRGB). We find a rapid increase in the number of detected stars fainter than $m_{F814W} \sim 27.3$, which we identify as the onset of the red giant branch. Using a forward modeling approach that takes the photometric uncertainties into account, we find a TRGB magnitude of $m_{F814W,\rm TRGB}=27.47 \pm 0.16$. The inferred distance, including the uncertainty in the absolute calibration, is $D_{\rm TRGB}=20.0 \pm 1.6$ Mpc. The TRGB distance of NGC1052-DF4 is consistent with the previously-determined surface brightness fluctuation distance of $D_{\rm SBF}=18.7\pm 1.7$ Mpc to NGC1052-DF2 and is consistent with the distance of the bright elliptical galaxy NGC1052. We conclude that the unusual properties of these galaxies cannot be explained by distance errors., Comment: Accepted for publication in ApJ Letters, in press. Figure 1 shows the color image of the galaxy. The main result is shown in Figure 4

We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($\tau \sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites., Comment: 7 pages, 3 figures, 1 table. Submitted to ApJL and posted after addressing the referee's report. Comments welcome

We present the star formation history (SFH) of the isolated (D~970 kpc) Local Group dwarf galaxy WLM measured from color-magnitude diagrams constructed from deep Hubble Space Telescope imaging. Our observations include a central (0.5 $r_h$) and outer field (0.7 $ r_h$) that reach below the oldest main sequence turnoff. WLM has no early dominant episode of star formation: 20% of its stellar mass formed by ~12.5 Gyr ago (z~5). It also has an SFR that rises to the present with 50% of the stellar mass within the most recent 5 Gyr (z<0.7). There is evidence of a strong age gradient: the mean age of the outer field is 5 Gyr older than the inner field despite being only 0.4 kpc apart. Some models suggest such steep gradients are associated with strong stellar feedback and dark matter core creation. The SFHs of real isolated dwarf galaxies and those from the the Feedback In Realistic Environment suite are in good agreement for $M_{\star}(z=0) \sim 10^7-10^9 M_{\odot}$, but in worse agreement at lower masses ($M_{\star}(z=0) \sim 10^5-10^7 M_{\odot}$). These differences may be explainable by systematics in the models (e.g., reionization model) and/or observations (HST field placement). We suggest that a coordinated effort to get deep CMDs between HST/JWST (crowded central fields) and WFIRST (wide-area halo coverage) is the optimal path for measuring global SFHs of isolated dwarf galaxies., Comment: 13 pages, 13 Figures, 4 Tables. Re-submitted to MNRAS after addressing the referee's comments