Search

The Local Ultraviolet to Infrared Treasury (LUVIT) is a Hubble Space Telescope program that combines newly acquired data in the near ultraviolet (NUV), optical, and near infrared (NIR) with archival optical and NIR imaging to produce multiband panchromatic resolved stellar catalogs for 23 pointings in 22 low-mass, star-forming galaxies ranging in distance from the outskirts of the Local Group to ~3.8 Mpc. We describe the survey design, detail the LUVIT broadband filter observations and the archival datasets included in the LUVIT reductions, and summarize the simultaneous multiband data reduction steps. The spatial distributions and color-magnitude diagrams (CMDs) from the resulting stellar catalogs are presented for each target, from the NUV to the NIR. We demonstrate in which regions of the CMDs stars with NUV and optical, optical and NIR, and NUV through NIR detections reside. For each target, we use the results from artificial star tests to measure representative completeness, bias, and total photometric uncertainty as a function of magnitude in each broadband filter. We also assess which LUVIT targets have significant spatial variation in the fraction of stars recovered at a given magnitude. The panchromatic LUVIT stellar catalogs will provide a rich legacy dataset for a host of resolved stellar population studies., Comment: 48 pages, 14 figures, 8 tables, accepted for publication in ApJS

We present archival HST and JWST ultraviolet through near infrared time series photometric observations of a massive minimal-contact binary candidate in the metal-poor galaxy WLM ($Z = 0.14 Z_{\odot}$). This discovery marks the lowest metallicity contact binary candidate observed to date. We determine the nature of the two stars in the binary by using the eclipsing binary modeling software (PHysics Of Eclipsing BinariEs; PHOEBE) to train a neural network to fit our observed panchromatic multi-epoch photometry. The best fit model consists of two hot MS stars ($T_1=29800^{+2300}_{-1700}$ K, $M_1=16^{+2}_{-3}~M_{\odot}$, and $T_2=18000^{+5000}_{-5000}$ K, $M_2=7^{+5}_{-3}~M_{\odot}$). We discuss plausible evolutionary paths for the system, and suggest the system is likely to be currently in a contact phase before ultimately ending in a merger. Future spectroscopy will help to further narrow down evolutionary pathways. This work showcases a novel use of data of JWST and HST imaging originally taken to characterize RR Lyrae. We expect time series imaging from LSST, BlackGEM, etc. to uncover similar types of objects in nearby galaxies., Comment: comments welcome

We demonstrate how near infrared (NIR) imaging of resolved luminous asymptotic giant branch (AGB) stars can be used to measure well-constrained star formation histories (SFHs) across cosmic time. Using UKIRT J and K-band imaging of M31, we first show excellent agreement over the past $\sim8$ Gyr between the PHAT SFH of M31's outer disk derived from a deep optical color-magnitude diagram (CMD; $\sim3.3\times10^{7}$ stars with $M_{F814W} \lesssim +2$), and our spatially-matched SFH based only on modeling AGB stars on a NIR CMD ($\sim2.3\times10^{4}$ stars with $M_{J} \lesssim -5$). We find that only hundreds of AGB stars are needed for reliable SFH recovery, owing to their excellent age sensitivity in the NIR. We then measure the spatially resolved SFH of M31's inner stellar halo ($D_{M31, projected} \sim20-30$ kpc) using $\sim10^4$ AGB stars. We find: (i) a dominant burst of star formation across M31's inner stellar halo from $4-5$ Gyr ago and lower level, spatially distributed star formation $\sim1-2$ Gyr ago; (ii) a younger 'quenching time' in the vicinity of NGC 205 ($\sim1$ Gyr ago) than near M32 ($\sim1.6$ Gyr ago); (iii) $M_{\star}\sim4\pm0.5\times10^9 M_{\odot}$ formed over the past $\sim8$ Gyr. We discuss some caveats and the enormous potential of resolved AGB stars in the NIR for measuring SFHs back to ancient epochs ($\sim13$ Gyr ago) in galaxies to large distances ($D\gtrsim20$ Mpc) with JWST, Roman, and Euclid., Comment: 25 pages, 19 figures, submitted to AAS Journals, typos in abstract were fixed in v2

We examine star-formation and dust properties for a sample of 660 galaxies at $1.37\leq z\leq 2.61$ in the MOSDEF survey by dividing them into groups with similarly-shaped spectral energy distributions (SEDs). For each group, we combine the galaxy photometry into a finely-sampled composite SED, and stack their spectra. This method enables the study of more complete galaxy samples, including galaxies with very faint emission lines. We fit these composite SEDs with Prospector to measure the stellar attenuation and SED-based star-formation rates (SFRs). We also derive emission-line properties from the spectral stacks, including Balmer decrements, dust-corrected SFRs, and metallicities. We find that stellar attenuation correlates most strongly with mass, while nebular attenuation correlates strongly with both mass and SFR. Furthermore, the excess of nebular compared to stellar attenuation correlates most strongly with SFR. The highest SFR group has 2 mag of excess nebular attenuation. Our results are consistent with a model in which star-forming regions become more dusty as galaxy mass increases. To explain the increasing excess nebular attenuation, we require a progressively larger fraction of star formation to occur in highly-obscured regions with increasing SFR. This highly-obscured star formation could occur in dusty clumps or central starbursts. Additionally, as each galaxy group represents a different evolutionary stage, we study their locations on the UVJ and SFR-mass diagrams. As mass increases, metallicity and dust attenuation increase, while sSFR decreases. However, the most massive group moves towards the quiescent region of the UVJ diagram, while showing less obscuration, potentially indicating removal of dust., Comment: 21 pages, 7 figures, accepted for publication in ApJ

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 present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log $M/M_\odot=10.2-11.2$) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.26$\pm0.04$ dex in [C/H], 0.16$\pm0.03$ dex in [Fe/H], and 0.07$\pm0.04$ dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies forming at higher redshift consistently show higher [Mg/Fe] and lower [Fe/H] and [Mg/H], regardless of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily driven by lower-redshift samples naturally including galaxies with longer star-formation timescales. In contrast, the lower [Mg/H] likely reflects earlier-forming galaxies expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is somewhat lower at z=1-3 compared to the lower redshift relation. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies., Comment: Accepted to ApJ; 22 pages, 6 figures, 3 tables

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 Astro2020 Decadal Survey recognizes time-domain astronomy as a key science area over the next decade and beyond. With over 30 years of HST data and the potential for 20 years of JWST operations, these flagship observatories offer an unparalleled prospect for a half-century of space-based observations in the time domain. To take best advantage of this opportunity, STScI charged a working group to solicit community input and formulate strategies to maximize the science return in time-domain astronomy from these two platforms. Here, the HST/JWST Long-Term Monitoring Working Group reports on the input we received and presents our recommendations to enhance the scientific return for time-domain astronomy from HST and JWST. We suggest changes in policies to enable and prioritize long-term science programs of high scientific value. As charged, we also develop recommendations based on community input for a JWST Director's Discretionary Time program to observe high-redshift transients., Comment: 17 pages, 3 figures, minor typos corrected, working group information at https://outerspace.stsci.edu/display/HPR/Long-term+variability+monitoring+strategies+for+HST+and+JWST

We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ($z=1-3$, $H_{AB}\le23$), as well as 53 star-forming galaxies at $z=1-4$. With 16~hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log$M_*/M_{\odot}\sim10.2-11.5$) and detailed star-formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages ($0.8-3.3$~Gyr) and star formation timescales ($\sim0.5-4$~Gyr), and four out of 20 galaxies were already quiescent by $z=3$. On average, the $z>1.75$ $[z<1.75]$ galaxies formed 50\% of their stellar mass before $z=4$ $[z=3]$. Furthermore, the typical SFHs of galaxies in these two redshift bins ($z_{\text{mean}}=2.2~[1.3]$) indicate that galaxies at higher redshift formed earlier and over shorter star-formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, number density calculations imply that mergers and/or late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances, resolved stellar populations and kinematics., Comment: Accepted in ApJ; 25 pages, 14 figures, 2 tables (excluding appendices)

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

We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of the galaxies in our sample are on first infall with $>$90\% certainty. The observed star formation histories (SFHs) of these first-infall dwarfs are generally consistent with infalling dwarfs in simulations. The remaining four galaxies have crossed the virial radius of either the MW or M31. When we compare their star formation (SF) and orbital histories we find tentative agreement between the inferred pattern of SF with the timing of dynamical events in the orbital histories. For Leo~I, SF activity rises as the dwarf crosses the MW's virial radius, culminating in a burst of SF shortly before pericenter ($\approx1.7$~Gyr ago). The SF then declines after pericenter, but with some smaller bursts before its recent quenching ($\approx0.3$~Gyr ago). This shows that even small dwarfs like Leo~I can hold on to gas reservoirs and avoid quenching for several Gyrs after falling into their host, which is longer than generally found in simulations. Leo~II, NGC~6822, and IC~10 are also qualitatively consistent with this SF pattern in relation to their orbit, but more tentatively due to larger uncertainties., Comment: 27 pages, 9 Figures, 8 Tables, Submitted to ApJ

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

The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful `near-far' technique for studying the reionisation-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass (Mstar < 10^9 Msun) galaxies across 11 Milky Way- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterise their progenitors at z ~ 6 - 9, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with Mstar ~ 10^5 Msun (10^9 Msun) had ~1 (~30) progenitors at z ~ 7, and its main progenitor comprised ~100% (~50%) of the total stellar mass of all its progenitors at z ~ 7. We show that although only ~ 15% of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at z ~ 6 - 9. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just the low-mass galaxies at z = 0 and not the MW/M31 hosts correctly recovers the slope of the SMF down to Mstar ~ 10^4.5 Msun at z > 6. Thus, we validate the `near-far' approach as an unbiased tool for probing low-mass reionisation-era galaxies., Comment: 16 pages, 6 figures, 3 appendix figures, and 1 table. Submitted to MNRAS

Omega Centauri ($\omega$ Cen) is the most massive globular cluster of the Milky Way and has been the focus of many studies that reveal the complexity of its stellar populations and kinematics. However, most previous studies have used photometric and spectroscopic datasets with limited spatial or magnitude coverage, while we aim to investigate it having full spatial coverage out to its half-light radius and stars ranging from the main sequence to the tip of the red giant branch. This is the first paper in a new survey of $\omega$ Cen that combines uniform imaging and spectroscopic data out to its half-light radius to study its stellar populations, kinematics, and formation history. In this paper, we present an unprecedented MUSE spectroscopic dataset combining 87 new MUSE pointings with previous observations collected from guaranteed time observations. We extract spectra of more than 300,000 stars reaching more than two magnitudes below the main sequence turn-off. We use these spectra to derive metallicity and line-of-sight velocity measurements and determine robust uncertainties on these quantities using repeat measurements. Applying quality cuts we achieve signal-to-noise ratios of 16.47/73.51 and mean metallicity errors of 0.174/0.031 dex for the main sequence stars (18 mag $\rm < mag_{F625W}<$22 mag) and red giant branch stars (16 mag $<\rm mag_{F625W}<$10 mag), respectively. We correct the metallicities for atomic diffusion and identify foreground stars. This massive spectroscopic dataset will enable future studies that will transform our understanding of $\omega$ Cen, allowing us to investigate the stellar populations, ages, and kinematics in great detail., Comment: 27 pages, 18 figures, 3 tables, published in ApJ, the catalog is available in the online material of the published article; typos corrected in this version

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

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 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 investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at $1.3\leq z\leq2.6$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We divide galaxies with a detected H$\alpha$ emission line and coverage of H$\beta$ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), then stack their spectra. From each stack, we measure Balmer decrement and gas-phase metallicity, then we compute median \AV and UV continuum spectral slope ($\beta$). First, we find that none of the dust properties (Balmer decrement, \AV, $\beta$) vary with axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model where attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases SFR and lowers metallicity, but leaves dust column density mostly unchanged. We quantify this idea using the Kennicutt-Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation., Comment: 19 pages, 13 figures, accepted for publication in ApJ

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

Large-scale surveys will provide spectroscopy for $\sim$50 million resolved stars in the Milky Way and Local Group. However, these data will have a high degree of heterogeneity and most will be low-resolution ($R<10000$), posing challenges to measuring consistent and reliable stellar labels. Here, we introduce a framework for identifying and remedying these issues. By simultaneously fitting the full spectrum and Gaia photometry with the Payne, we measure $\sim$40 abundances for 8 red giants in M15. From degraded quality Keck/HIRES spectra, we evaluate trends with resolution and S/N and find that (i) $\sim$20 abundances are recovered consistently within $\lesssim$0.1 dex agreement and with $\lesssim$0.05-0.15~dex systematic uncertainties from $10000\lesssim R\lesssim80000$; (ii) for 9 elements (C, Mg, Ca, Sc, Ti, Fe, Ni, Y, Nd), this systematic precision and accuracy extends down to $R\sim2500$; and (iii) while most elements do not exhibit strong S/N-dependent systematics, there are non-negligible biases for 4 elements (C, Mg, Ca, and Dy) below $\text{S/N}\sim10$ pixel$^{-1}$. We compare statistical uncertainties from MCMC sampling to the easier-to-compute Cram\'er-Rao bounds and find that they agree for $\sim$75% of elements, indicating the latter to be a reliable and faster way to estimate uncertainties. Our analysis illustrates the great promise of low-resolution spectroscopy for stellar chemical abundance work, and ongoing improvements to stellar models (e.g., 3D-NLTE physics) will only further extend its viability to more elements and to higher precision and accuracy., Comment: 46 pages, 26 figures, submitted to ApJS. Comments welcome!

We present elemental abundance patterns (C, N, Mg, Si, Ca, Ti, V, Cr, Fe, Co, and Ni) for a population of 135 massive quiescent galaxies at $z\sim0.7$ with ultra-deep rest-frame optical spectroscopy drawn from the LEGA-C survey. We derive average ages and elemental abundances in four bins of stellar velocity dispersion ($\sigma_v$) ranging from 150$~$km$\,$s$^{-1}$ to 250$~$km$\,$s$^{-1}$ using a full-spectrum hierarchical Bayesian model. The resulting elemental abundance measurements are precise to 0.05$\,$dex. The majority of elements, as well as the total metallicity and stellar age, show a positive correlation with $\sigma_v$. Thus, the highest dispersion galaxies formed the earliest and are the most metal-rich. We find only mild or non-significant trends between [X/Fe] and $\sigma_v$, suggesting that the average star-formation timescale does not strongly depend on velocity dispersion. To first order, the abundance patterns of the $z\sim0.7$ quiescent galaxies are strikingly similar to those at $z\sim0$. However, at the lowest velocity dispersions the $z\sim0.7$ galaxies have slightly enhanced N, Mg, Ti, and Ni abundance ratios and earlier formation redshifts than their $z\sim0$ counterparts. Thus, while the higher-mass quiescent galaxy population shows little evolution, the low-mass quiescent galaxies population has grown significantly over the past six billion years. Finally, the abundance patterns of both $z\sim0$ and $z\sim0.7$ quiescent galaxies differ considerably from theoretical prediction based on a chemical evolution model, indicating that our understanding of the enrichment histories of these galaxies is still very limited., Comment: 24 pages, 12 figures, accepted for publication in ApJ

We build a statistical framework to infer the global properties of the satellite system of the Andromeda galaxy (M31) from the properties of individual dwarf galaxies located in the Pan-Andromeda Archaelogical Survey (PAndAS) and the previously determined completeness of the survey. Using forward modeling, we infer the slope of the luminosity function of the satellite system, the slope of its spatial density distribution, and the size-luminosity relation followed by the dwarf galaxies. We find that the slope of the luminosity function is $\beta=-1.5\pm0.1$. Combined with the spatial density profile, it implies that, when accounting for survey incompleteness, M31 hosts $92_{-26}^{+19}$ dwarf galaxies with $M_\textrm{V}<-5.5$ and a sky-projected distance from M31 between 30 and 300kpc. We conclude that many faint or distant dwarf galaxies remain to be discovered around Andromeda, especially outside the PAndAS footprint. Finally, we use our model to test if the higher number of satellites situated in the hemisphere facing the Milky Way could be explained simply by the detection limits of dwarf galaxy searches. We rule this out at $>99.9\%$ confidence and conclude that this anisotropy is an intrinsic feature of the M31 satellite system. The statistical framework we present here is a powerful tool to robustly constrain the properties of a satellite system and compare those across hosts, especially considering the upcoming start of the Euclid or Rubin large photometric surveys that are expected to uncover a large number of dwarf galaxies in the Local Volume., Comment: Submitted to ApJ - 12 pages, 6 figures, 2 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 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., Comment: 35 pages, 18 figures

Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The number and frequency of microlensing events therefore opens a window into the population of stars and compact objects, as well as high-redshift stars. To assemble the first statistical sample of stars in order to constrain the initial mass function (IMF) of massive stars at redshift z=0.7-1.5, the abundance of primordial black holes in galaxy-cluster dark matter, and the IMF of the stars making up the intracluster light, we are carrying out a 192-orbit program with the Hubble Space Telescope called "Flashlights," which is now two-thirds complete owing to scheduling challenges. We use the ultrawide F200LP and F350LP long-pass WFC3 UVIS filters and conduct two 16-orbit visits separated by one year. Having an identical roll angle during both visits, while difficult to schedule, yields extremely clean subtraction. Here we report the discovery of more than a dozen bright microlensing events, including multiple examples in the famous "Dragon Arc" discovered in the 1980s, as well as the "Spocks" and "Warhol" arcs that have hosted already known supergiants. The ultradeep observer-frame ultraviolet-through-optical imaging is sensitive to hot stars, which will complement deep James Webb Space Telescope infrared imaging. We are also acquiring Large Binocular Telescope LUCI and Keck-I MOSFIRE near-infrared spectra of the highly magnified arcs to constrain their recent star-formation histories.

We present novel constraints on the underlying galaxy formation physics (e.g., mass loading factor, star formation history, metal retention) at $z\gtrsim7$ for the low-mass ($M_*\sim10^5$ M$_\odot$) Local Group ultra-faint dwarf galaxy (UFD) Eridanus {\sc II} (Eri II). Using a hierarchical Bayesian framework, we apply a one-zone chemical evolution model to Eri II's CaHK-based photometric metallicity distribution function (MDF; [Fe/H]) and find that the evolution of Eri II is well-characterized by a short, exponentially declining star-formation history ($\tau_\text{SFH}=0.39\pm_{0.13}^{0.18}$ Gyr), a low star-formation efficiency ($\tau_\text{SFE}=27.56\pm_{12.92}^{25.14}$ Gyr), and a large mass-loading factor ($\eta=194.53\pm_{42.67}^{33.37}$). Our results are consistent with Eri II forming the majority of its stars before the end of reionization. The large mass-loading factor implies strong outflows in the early history of Eri II and is in good agreement with theoretical predictions for the mass-scaling of galactic winds. It also results in the ejection of $>$90\% of the metals produced in Eri II. We make predictions for the distribution of [Mg/Fe]-[Fe/H] in Eri II as well as the prevalence of ultra metal-poor stars, both of which can be tested by future chemical abundance measurements. Spectroscopic follow-up of the highest metallicity stars in Eri II ($\text{[Fe/H]} > -2$) will greatly improve model constraints. Our new framework can readily be applied to all UFDs throughout the Local Group, providing new insights into the underlying physics governing the evolution of the faintest galaxies in the reionization era., Comment: 20 pages; 12 figures, submitted to MNRAS

Andromeda (And) XXV has previously been reported as a dwarf spheroidal galaxy (dSph) with little-to-no dark matter. However, the uncertainties on this result were significant. In this study, we double the number of member stars and re-derive the kinematics and mass of And XXV. We find that And XXV has a systemic velocity of $\nu_\mathrm{r}=-107.7\pm1.0 \mathrm{~km s}^{-1}$ and a velocity dispersion of $\sigma_\nu=4.5\pm1.0\mathrm{~km s}^{-1}$. With this better constrained velocity dispersion, we derive a mass contained within the half-light radius of $M(r< r_\mathrm{h})=6.9^{+3.2}_{-2.8}\times10^6\mathrm{~M}_\odot$. This mass corresponds to a mass-to-light ratio of $\mathrm{[M/L]}_\mathrm{r_\mathrm{h}}=37^{+17}_{-15}\mathrm{~M}_\odot/\mathrm{L}_\odot$, demonstrating, for the first time, that And XXV has an unambiguous dark matter component. We also measure the metallicity of And XXV to be $\mathrm{[Fe/H]}=-1.9\pm0.1$$\mathrm{~}$dex, which is in agreement with previous results. Finally, we extend the analysis of And XXV to include mass modelling using GravSphere. We find that And XXV has a low central dark matter density, $\rho_\mathrm{DM}(150\mathrm{pc})= 2.7^{+1.8}_{-1.6}\times10^7\mathrm{~M}_\odot\mathrm{kpc}^{-3}$, making And XXV a clear outlier when compared to other Local Group (LG) dSphs of the similar stellar mass. In a companion paper, we will explore whether some combination of dark matter cusp-core transformations and/or tides can explain And XXV's low density., Comment: 13 pages, 8 figures (7 main, 1 appendix). Submitted to MNRAS

We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and astrometry constrains the companion mass to $M_2 = 9.62\pm 0.18\,M_{\odot}$. The spectroscopic orbit alone sets a minimum companion mass of $M_2>5\,M_{\odot}$; if the companion were a $5\,M_{\odot}$ star, it would be $500$ times more luminous than the entire system. These constraints are insensitive to the mass of the luminous star, which appears as a slowly-rotating G dwarf ($T_{\rm eff}=5850\,\rm K$, $\log g = 4.5$, $M=0.93\,M_{\odot}$), with near-solar metallicity ($\rm [Fe/H] = -0.2$) and an unremarkable abundance pattern. We find no plausible astrophysical scenario that can explain the orbit and does not involve a black hole. The orbital period, $P_{\rm orb}=185.6$ days, is longer than that of any known stellar-mass black hole binary. The system's modest eccentricity ($e=0.45$), high metallicity, and thin-disk Galactic orbit suggest that it was born in the Milky Way disk with at most a weak natal kick. How the system formed is uncertain. Common envelope evolution can only produce the system's wide orbit under extreme and likely unphysical assumptions. Formation models involving triples or dynamical assembly in an open cluster may be more promising. This is the nearest known black hole by a factor of 3, and its discovery suggests the existence of a sizable population of dormant black holes in binaries. Future Gaia releases will likely facilitate the discovery of dozens more., Comment: 30 pages, 15 figures, 1-2 black holes. Accepted to MNRAS

We present the method behind HubPUG, a software tool built for recovering systemic proper motions (PMs) of Hubble Space Telescope (HST) fields with two epochs of observations by utilizing stars observed by Gaia as a foreground frame of reference. HST PM experiments have typically relied on the use of distant background galaxies or quasi-stellar objects (QSOs) as stationary sources against which to measure PMs. Without consistent profiles, background galaxies are more difficult to centroid, but benefit on-aggregate from their large numbers. QSOs, though they can be fit with stellar point-spread functions, are sparse, with most fields containing none. Historically, the use of stars as references against which to measure PMs would have been difficult because they have individual PMs of their own. However, Gaia has now provided positions and PMs for over 1.4 billion stars, which are much more likely to be well-imaged in the fields around targets versus background sources, have predictable stellar profiles, and require less observing time per-image for good signal-to-noise. This technique allows us to utilize the power of Gaia to measure the PM of targets too faint for Gaia to observe itself. We have recovered PMs for the Milky Way satellites Sculptor and Draco with comparable uncertainties over HST-only and Gaia-only measurements, limited primarily by the current capabilities of the Gaia data. We also show the promise of this method for satellites of M31 with a new PM measurement for Andromeda VII.

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 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., Comment: 24 pages, 15 figures, 9 tables, Accepted to ApJ

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

We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$), self-bound stellar clusters form at $11\gtrsim z \gtrsim 5$ in progenitors with $M_{\rm vir} \approx 10^9\,{\rm M}_\odot$. Cluster formation is triggered when at least $10^7\,{\rm M}_\odot$ of dense, turbulent gas reaches $\Sigma_{\rm gas} \approx 10^4\, {\rm M}_\odot\, {\rm pc}^{-2}$ as a result of the compressive effects of supernova feedback or from cloud-cloud collisions. The clusters can survive for $2-3\,{\rm Gyr}$; absent numerical effects, they would likely survive substantially longer, perhaps to $z=0$. The longest-lived clusters are those that form at significant distance -- several hundreds of pc -- from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter halos as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in \textit{Hubble Space Telescope} lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with the \textit{James Webb Space Telescope}. By contrast, the star clusters' host galaxies will remain undetectable., Comment: 14 pages, 5 figures, submitted to MNRAS

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)

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program., Comment: 69 pages, 43 figures

We provide uniform RR Lyrae-based distances to 39 dwarf galaxies in and around the Local Group. We determine distances based on a Bayesian hierarchical model that uses periods and magnitudes of published RR Lyrae in dwarf galaxies and is anchored to well-studied Milky Way (MW) RR Lyrae with spectroscopic metallicities and Gaia eDR3 parallaxes. Gaia eDR3 parallaxes for the anchor sample are a factor of 2, on average, more precise than DR2 parallaxes, and allow for a much better constrained period-luminosity-metallicity relation. While $\sim75$% of our distances are within 1-$\sigma$ of recent literature RR Lyrae distances, our distances are also $\sim2$-$3$ times more precise than literature distances, on average. On average, our distances are $\sim0.05$ mag closer than literature distances, as well as $\sim0.06$ mag closer than distances derived using a theoretical period-luminosity-metallicity relation. These discrepancies are largely due to our eDR3 parallax anchor. We show that the Hipparcos-anchored RR Lyrae distance scale of Carretta et al. (2000) over-predicts distances to MW RR Lyrae by $\sim0.05$ mag. The largest uncertainties in our distances are (i) the lack of direct metallicity measurements for RR Lyrae and (ii) the heterogeneity of published RR Lyrae photometry. We provide simple formulae to place new dwarf galaxies with RR Lyrae on a common distance scale with this work. We provide public code that can easily incorporate additional galaxies and data from future surveys, providing a versatile framework for cartography of the local Universe with RR Lyrae., Comment: 34 pages, 15 figures, Published in ApJ

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!

Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$ arcmin$^{2}$ (${\sim}12$ kpc$^{2}$) of the nearby galaxy NGC 300 to quantify the effect of stellar feedback as a function of the local galactic environment. We extract spectra from emission line regions identified within dendrograms, combine emission line ratios and line widths to distinguish between HII regions, planetary nebulae, and supernova remnants, and compute their ionised gas properties, gas-phase oxygen abundances, and feedback-related pressure terms. For the HII regions, we find that the direct radiation pressure ($P_\mathrm{dir}$) and the pressure of the ionised gas ($P_{HII}$) weakly increase towards larger galactocentric radii, i.e. along the galaxy's (negative) abundance and (positive) extinction gradients. While the increase of $P_{HII}$ with galactocentric radius is likely due to higher photon fluxes from lower-metallicity stellar populations, we find that the increase of $P_\mathrm{dir}$ is likely driven by the combination of higher photon fluxes and enhanced dust content at larger galactocentric radii. In light of the above, we investigate the effect of increased pre-supernova feedback at larger galactocentric distances (lower metallicities and increased dust mass surface density) on the ISM, finding that supernovae at lower metallicities expand into lower-density environments, thereby enhancing the impact of supernova feedback., Comment: Accepted for publication in MNRAS

We present LAMOST J0140355+392651 (hereafter J0140), a close ($P_{\rm orb} = 3.81$ hours) binary containing a bloated, low-mass ($M \approx 0.15 M_{\odot}$) proto-white dwarf (WD) and a massive ($M\approx 0.95\,M_{\odot}$) WD companion. The system's optical light curve is dominated by large-amplitude ellipsoidal variability but also exhibits additional scatter, likely driven by pulsations. The proto-WD is cooler ($T_{\rm eff} = 6800\pm 100$ K) and more puffy ($\log\left[g/\left({\rm cm\,s^{-2}}\right)\right]=4.74\pm0.07$) than any known extremely low mass (ELM) WD, but hotter than any known cataclysmic variable (CV) donor. It either completely or very nearly fills its Roche lobe ($R/R_{{\rm Roche\,lobe}}=0.99\pm0.01$), suggesting ongoing or recently terminated mass transfer. No dwarf nova-like outbursts have been observed. The spectrum is dominated by the proto-WD but shows tentative hints of H$\alpha$ emission, perhaps due to accretion onto the massive WD. The properties of the system are well-matched by MESA binary evolution models of CVs with donors that underwent significant nuclear evolution before the onset of mass transfer. In these models, the bloated proto-WD is either still losing mass via stable Roche lobe overflow or was doing so until very recently. In either case, it is evolving toward higher temperatures at near-constant luminosity to become an ELM WD. If the system is detached, mass transfer likely ended when the donor became too hot for magnetic braking to remain efficient. Evolutionary models predict that the binary will shrink to $P_{\rm orb}\lesssim 10$ minutes within a few Gyr, when it will either merge or become an AM CVn binary. J0140 provides an observational link between the formation channels of CVs, ELM WDs, detached ultracompact WD binaries, and AM CVn systems., Comment: 19 pages, 14 figures. Accepted to MNRAS

Planck data provide precise constraints on cosmological parameters when assuming the base $\Lambda$CDM model, including a $0.17\%$ measurement of the age of the Universe, $t_0=13.797 \pm 0.023\,{\rm Gyr}$. However, the persistence of the "Hubble tension" calls the base $\Lambda$CDM model's completeness into question and has spurred interest in models such as Early Dark Energy (EDE) that modify the assumed expansion history of the Universe. We investigate the effect of EDE on the redshift-time relation $z \leftrightarrow t$ and find that it differs from the base $\Lambda$CDM model by at least ${\approx} 4\%$ at all $t$ and $z$. As long as EDE remains observationally viable, any inferred $t \leftarrow z$ or $z \leftarrow t$ quoted to a higher level of precision do not reflect the current status of our understanding of cosmology. This uncertainty has important astrophysical implications: the reionization epoch - $10>z>6$ - corresponds to disjoint lookback time periods in the base $\Lambda$CDM and EDE models, and the EDE value of $t_0=13.25 \pm 0.17~{\rm Gyr}$ is in tension with published ages of some stars, star clusters, and ultra-faint dwarf galaxies. However, most published stellar ages do not include an uncertainty in accuracy (due to, e.g., uncertain distances and stellar physics) that is estimated to be $\sim7-10\%$, potentially reconciling stellar ages with $t_{0,\rm EDE}$. We discuss how the big data era for stars is providing extremely precise ages ($<1\%$) and how improved distances and treatment of stellar physics such as convection could result in ages accurate to $4-5\%$, comparable to the current accuracy of $t \leftrightarrow z$. Such precise and accurate stellar ages can provide detailed insight into the high-redshift Universe independent of a cosmological model., Comment: 20 pages, 9 figures; MNRAS (accepted). V2: minor revisions and clarifications; conclusions unchanged

Andromeda XXI (And XXI) has been proposed as a dwarf spheroidal galaxy with a central dark matter density that is lower than expected in the Standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmology. In this work, we present dynamical observations for 77 member stars in this system, more than doubling previous studies to determine whether this galaxy is truly a low density outlier. We measure a systemic velocity of $v_r=-363.4\pm1.0\,{\rm kms}^{-1}$ and a velocity dispersion of $\sigma_v=6.1^{+1.0}_{-0.9}\,{\rm kms}^{-1}$, consistent with previous work and within $1\sigma$ of predictions made within the modified Newtonian dynamics framework. We also measure the metallicity of our member stars from their spectra, finding a mean value of ${\rm [Fe/H]}=-1.7\pm0.1$~dex. We model the dark matter density profile of And~XXI using an improved version of \GravSphere, finding a central density of $\rho_{\rm DM}({\rm 150 pc})=2.7_{-1.7}^{+2.7} \times 10^7 \,{\rm M_\odot\,kpc^{-3}}$ at 68\% confidence, and a density at two half light radii of $\rho_{\rm DM}({\rm 1.75 kpc})=0.9_{-0.2}^{+0.3} \times 10^5 \,{\rm M_\odot\,kpc^{-3}}$ at 68\% confidence. These are both a factor ${\sim}3-5$ lower than the densities expected from abundance matching in $\Lambda$CDM. We show that this cannot be explained by `dark matter heating' since And~XXI had too little star formation to significantly lower its inner dark matter density, while dark matter heating only acts on the profile inside the half light radius. However, And~XXI's low density can be accommodated within $\Lambda$CDM if it experienced extreme tidal stripping (losing $>95\%$ of its mass), or if it inhabits a low concentration halo on a plunging orbit that experienced repeated tidal shocks., Comment: Accepted to MNRAS, 17 pages, 12 figures

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

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

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