Showing 1–50 of 217 results for author: Dalcanton, J J

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

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

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

Submitted 17 June, 2024; originally announced June 2024.

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

The Vela Pulsar Progenitor Was Most Likely a Binary Merger

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

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

Submitted 6 June, 2024; originally announced June 2024.

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

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine learning and Simulations

Authors: Jonah C. Rose, Paul Torrey, Francisco Villaescusa-Navarro, Mariangela Lisanti, Tri Nguyen, Sandip Roy, Kassidy E. Kollmann, Mark Vogelsberger, Francis-Yan Cyr-Racine, Mikhail V. Medvedev, Shy Genel, Daniel Anglés-Alcázar, Nitya Kallivayalil, Bonny Y. Wang, Belén Costanza, Stephanie O'Neil, Cian Roche, Soumyodipta Karmakar, Alex M. Garcia, Ryan Low, Shurui Lin, Olivia Mostow, Akaxia Cruz, Andrea Caputo, Arya Farahi , et al. (5 additional authors not shown)

Abstract: We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- f… ▽ More We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of Warm Dark Matter, each comprised of 1024 simulations generated using the Arepo code. One suite consists of uniform-box simulations covering a $(25~h^{-1}~{\rm M}_\odot)^3$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the Warm Dark Matter particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and Convolutional Neural Networks, to demonstrate how such simulation suites can be used to disentangle the effects of dark matter and baryonic physics on galactic properties. The DREAMS project can be extended further to include different dark matter models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of dark matter on these scales. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 28 pages, 8 figures, DREAMS website: https://www.dreams-project.org

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

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

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

Submitted 5 February, 2024; originally announced February 2024.

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

JWST Reveals Star Formation Across a Spiral Arm in M33

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

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

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

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

Orbital Torus Imaging: Acceleration, density, and dark matter in the Galactic disk measured with element abundance gradients

Authors: Danny Horta, Adrian M. Price-Whelan, David W. Hogg, Kathryn V. Johnston, Lawrence Widrow, Julianne J. Dalcanton, Melissa K. Ness, Jason A. S. Hunt

Abstract: Under the assumption of a simple and time-invariant gravitational potential, many Galactic dynamics techniques infer the Milky Way's mass and dark matter distribution from stellar kinematic observations. These methods typically rely on parameterized potential models of the Galaxy and must take into account non-trivial survey selection effects, because they make use of the density of stars in phase… ▽ More Under the assumption of a simple and time-invariant gravitational potential, many Galactic dynamics techniques infer the Milky Way's mass and dark matter distribution from stellar kinematic observations. These methods typically rely on parameterized potential models of the Galaxy and must take into account non-trivial survey selection effects, because they make use of the density of stars in phase space. Large-scale spectroscopic surveys now supply information beyond kinematics in the form of precise stellar label measurements (especially element abundances). These element abundances are known to correlate with orbital actions or other dynamical invariants. Here, we use the Orbital Torus Imaging (OTI) framework that uses abundance gradients in phase space to map orbits. In many cases these gradients can be measured without detailed knowledge of the selection function. We use stellar surface abundances from the APOGEE survey combined with kinematic data from the Gaia mission. Our method reveals the vertical ($z$-direction) orbit structure in the Galaxy and enables empirical measurements of the vertical acceleration field and orbital frequencies in the disk. From these measurements, we infer the total surface mass density, $Σ$, and midplane volume density, $ρ_0$, as a function of Galactocentric radius and height. Around the Sun, we find $Σ_{\odot}(z=1.1$ kpc)$=72^{+6}_{-9}$M$_{\odot}$pc$^{-2}$ and $ρ_{\odot}(z=0)=0.081^{+0.015}_{-0.009}$ M$_{\odot}$pc$^{-3}$ using the most constraining abundance ratio, [Mg/Fe]. This corresponds to a dark matter contribution in surface density of $Σ_{\odot,\mathrm{DM}}(z=1.1$ kpc)$=24\pm4$ M$_{\odot}$pc$^{-2}$, and in total volume mass density of $ρ_{\odot,\mathrm{DM}}(z=0)=0.011\pm0.002$ M$_{\odot}$pc$^{-3}$. Moreover, using these mass density values we estimate the scale length of the low-$α$ disc to be $h_R=2.24\pm0.06$kpc. △ Less

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

Comments: Accepted for publication in ApJ. 19 pages, 11 figures, 3 Tables

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

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

Abstract: Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around th… ▽ More Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use them

Dissecting the Gaia HR diagram II. The vertical structure of the star formation history across the Solar Cylinder

Authors: Alessandro Mazzi, Léo Girardi, Michele Trabucchi, Julianne J. Dalcanton, Rodrigo Luger, Paola Marigo, Andrea Miglio, Guglielmo Costa, Yang Chen, Giada Pastorelli, Morgan Fouesneau, Simone Zaggia, Alessandro Bressan, Piero Dal Tio

Abstract: Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Ga… ▽ More Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Galactic Disc stellar component, which is well-described by power-laws with indices ranging from $1/2$ to $2/3$. The vertically-integrated star formation rate falls from $(1.147 \pm 0.039)\times10^{-8}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at earlier times down to $(6.2 \pm 3.0) \times10^{-9}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at present times, but we find a significant peak of star formation in the 2 to 3 Gyr age bin. The total mass of stars formed per unit area over time is $118.7 \pm 6.2\, \text{M}_{\odot} \text{pc}^{-2}$, which is nearly twice the present stellar mass derived from kinematics within 1~kpc from the Galactic Plane, implying a high degree of matter recycling in successive generations of stars. The method is then modified by adopting an age-dependent correlation between the SFH across the different slices, which results in less noisy and more symmetrical results without significantly changing the previously mentioned quantities. This appears to be a promising way to improve SFH recovery in external galaxies. △ Less

Submitted 23 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Empirical 2MASS-WFC3/IR filter transformations from synthetic photometry

Authors: M. J. Durbin, R. L. Beaton, A. J. Monson, B. Swidler, J. J. Dalcanton

Abstract: Near-infrared bandpasses on spaceborne observatories diverge from their ground-based counterparts as they are free of atmospheric telluric absorption. Available transformations between respective filter systems in the literature rely on theoretical stellar atmospheres, which are known to have difficulties reproducing observed spectral energy distributions of cool giants. We present new transformat… ▽ More Near-infrared bandpasses on spaceborne observatories diverge from their ground-based counterparts as they are free of atmospheric telluric absorption. Available transformations between respective filter systems in the literature rely on theoretical stellar atmospheres, which are known to have difficulties reproducing observed spectral energy distributions of cool giants. We present new transformations between the 2MASS $JHK_S$ and HST WFC3/IR F110W, F125W, & F160W photometric systems based on synthetic photometry of empirical stellar spectra from four spectral libraries. This sample comprises over 1000 individual stars, which together span nearly the full HR diagram and sample stellar populations from the solar neighborhood out to the Magellanic Clouds, covering a broad range of ages, metallicities, and other relevant stellar properties. In addition to global color-dependent transformations, we examine band-to-band differences for cool, luminous giant stars in particular, including multiple types of primary distance indicators. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 24 pages, 15 figures, accepted to AJ

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

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

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

Submitted 22 August, 2023; originally announced August 2023.

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

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

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

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

Submitted 18 July, 2023; originally announced July 2023.

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

KPM: A Flexible and Data-Driven K-Process Model for Nucleosynthesis

Authors: Emily J. Griffith, David W. Hogg, Julianne J. Dalcanton, Sten Hasselquist, Bridget Ratcliffe, Melissa Ness, David H. Weinberg

Abstract: The element abundance pattern found in Milky Way disk stars is close to two-dimensional, dominated by production from one prompt process and one delayed process. This simplicity is remarkable, since the elements are produced by a multitude of nucleosynthesis mechanisms operating in stars with a wide range of progenitor masses. We fit the abundances of 14 elements for 48,659 red-giant stars from AP… ▽ More The element abundance pattern found in Milky Way disk stars is close to two-dimensional, dominated by production from one prompt process and one delayed process. This simplicity is remarkable, since the elements are produced by a multitude of nucleosynthesis mechanisms operating in stars with a wide range of progenitor masses. We fit the abundances of 14 elements for 48,659 red-giant stars from APOGEE DR17 using a flexible, data-driven K-process model -- dubbed KPM. In our fiducial model, with $K=2$, each abundance in each star is described as the sum of a prompt and a delayed process contribution. We find that KPM with $K=2$ is able to explain the abundances well, recover the observed abundance bimodality, and detect the bimodality over a greater range in metallicity than previously has been possible. We compare to prior work by Weinberg et al. (2022), finding that KPM produces similar results, but that KPM better predicts stellar abundances, especially for elements C+N and Mn and for stars at super-solar metallicities. The model fixes the relative contribution of the prompt and delayed process to two elements to break degeneracies and improve interpretability; we find that some of the nucleosynthetic implications are dependent upon these detailed choices. We find that moving to four processes adds flexibility and improves the model's ability to predict the stellar abundances, but doesn't qualitatively change the story. The results of KPM will help us to interpret and constrain the formation of the Galaxy disk, the relationship between abundances and ages, and the physics of nucleosynthesis. △ Less

Submitted 8 December, 2023; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: 33 pages, 13 figures, 3 tables

RomAndromeda: The Roman Survey of the Andromeda Halo

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

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

Submitted 21 June, 2023; originally announced June 2023.

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

Origins of the Evil Eye: M64's Stellar Halo Reveals the Recent Accretion of an SMC-mass Satellite

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Jeremy Bailin, Julianne J. Dalcanton, Roelof S. de Jong, Richard D'Souza, Katya Gozman, In Sung Jang, Antonela Monachesi, David Nidever, Colin T. Slater

Abstract: M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's out… ▽ More M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's outer disk have shown no evidence of such an event, leading to other interpretations, such as a "flyby" interaction with the distant diffuse satellite Coma P. We present Subaru Hyper Suprime-Cam observations of M64's stellar halo, which resolve its stellar populations and reveal a spectacular radial shell feature, oriented $\sim$30$^{\circ}$ relative to the major axis and along the rotation axis of the outer gas disk. The shell is $\sim$45 kpc southeast of M64, while a similar but more diffuse plume to the northwest extends to $>$100 kpc. We estimate a stellar mass and metallicity for the southern shell of $M_{\star} {=} 1.80~{\pm}~0.54{\times}10^8~M_{\odot}$ and [M/H] $=$ $-$1.0, respectively, and a similar mass of $1.42~{\pm}~0.71{\times}10^8 M_{\odot}$ for the northern plume. Taking into account the accreted material in M64's inner disk, we estimate a total stellar mass for the progenitor satellite of $M_{\rm \star,prog}~{\simeq}~5{\times}10^8~M_{\odot}$. These results suggest that M64 is in the final stages of a minor merger with a gas-rich satellite strikingly similar to the SMC, in which M64's accreted counter-rotating gas originated, and which is responsible for the formation of its dusty inner star-forming disk. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 5 figures, accepted for publication in ApJ Letters

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

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

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

Submitted 25 May, 2023; originally announced May 2023.

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

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

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

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

Submitted 17 April, 2023; originally announced April 2023.

Comments: 22 pages. Accepted to the Astrophysical Journal

The Masses of Supernova Remnant Progenitors in M33

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

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

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

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

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

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

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

Submitted 6 March, 2023; originally announced March 2023.

Comments: 32 pages, 18 figures, 4 tables

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

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

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

Submitted 11 January, 2023; originally announced January 2023.

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

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

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

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

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

Comments: 35 pages, 18 figures

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

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

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

Submitted 24 August, 2022; originally announced August 2022.

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

The Spectroscopy and H-band Imaging of Virgo cluster galaxies (SHIVir) Survey: Data Catalogue and Kinematic Profiles

Authors: Nathalie N. -Q. Ouellette, Stéphane Courteau, Jon A. Holtzman, Michael McDonald, Michele Cappellari, Joel C. Roediger, Patrick Côté, Julianne J. Dalcanton, Elena Dalla Bontà, Laura Ferrarese, R. Brent Tully, Connor Stone, Eric W. Peng

Abstract: The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reductio… ▽ More The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reduction, and analysis for this SHIVir sample. We have used and optimised the \texttt{pPXF} routine to extract stellar kinematics from our data. Ultimately, resolved kinematic profiles (rotation curves and velocity dispersion profiles) are available for 133 SHIVir galaxies. A comprehensive database of photometric and kinematic parameters for the SHIVir sample is presented with: grizH magnitudes, effective surface brightnesses, effective and isophotal radii, rotational velocities, velocity dispersions, and stellar and dynamical masses. Parameter distributions highlight some bimodal distributions and possible sample biases. A qualitative study of resolved extended velocity dispersion profiles suggests a link between the so-called ``sigma-drop'' kinematic profile and the presence of rings in lenticular S0 galaxies. Rising dispersion profiles are linked to early-type spirals or dwarf ellipticals for which a rotational component is significant, whereas peaked profiles are tied to featureless giant ellipticals. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 21 pages, 16 figures, accepted for publication in MNRAS (doi: 10.1093/mnras/stac1347). Supplemental material available at https://www.queensu.ca/academia/courteau/links

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 2356-2375

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

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

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

Submitted 22 June, 2022; originally announced June 2022.

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

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

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

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

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

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

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

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

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

Submitted 29 October, 2021; originally announced October 2021.

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

Progenitor Mass Distribution for 22 Historic Core-Collapse Supernovae

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

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

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

Old and New Major Mergers in the SOSIMPLE galaxy, NGC 7135

Authors: Thomas A. Davison, Harald Kuntschner, Bernd Husemann, Mark A. Norris, Julianne J. Dalcanton, Alessandra De Rosa, Pierre-Alain Duc, Stefano Bianchi, Pedro R. Capelo, Cristian Vignali

Abstract: The simultaneous advancement of high resolution integral field unit spectroscopy and robust full-spectral fitting codes now make it possible to examine spatially-resolved kinematic, chemical composition, and star-formation history from nearby galaxies. We take new MUSE data from the Snapshot Optical Spectroscopic Imaging of Mergers and Pairs for Legacy Exploration (SOSIMPLE) survey to examine NGC… ▽ More The simultaneous advancement of high resolution integral field unit spectroscopy and robust full-spectral fitting codes now make it possible to examine spatially-resolved kinematic, chemical composition, and star-formation history from nearby galaxies. We take new MUSE data from the Snapshot Optical Spectroscopic Imaging of Mergers and Pairs for Legacy Exploration (SOSIMPLE) survey to examine NGC 7135. With counter-rotation of gas, disrupted kinematics and asymmetric chemical distribution, NGC 7135 is consistent with an ongoing merger. Though well hidden by the current merger, we are able to distinguish stars originating from an older merger, occurring 6-10 Gyr ago. We further find a gradient in ex-situ material with galactocentric radius, with the accreted fraction rising from 0% in the galaxy centre, to ~7% within 0.6 effective radii. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Comments: 13 pages, 9 figures, Accepted for publication in MNRAS

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

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

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

Submitted 4 January, 2021; originally announced January 2021.

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

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

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

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

Submitted 21 October, 2020; originally announced October 2020.

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

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

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

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

Submitted 3 September, 2020; originally announced September 2020.

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

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Large Magellanic Cloud

Authors: Giada Pastorelli, Paola Marigo, Léo Girardi, Bernhard Aringer, Yang Chen, Stefano Rubele, Michele Trabucchi, Sara Bladh, Martha L. Boyer, Alessandro Bressan, Julianne J. Dalcanton, Martin A. T. Groenewegen, Thomas Lebzelter, Nami Mowlavi, Katy L. Chubb, Maria-Rosa L. Cioni, Richard de Grijs, Valentin D. Ivanov, Ambra Nanni, Jacco Th. van Loon, Simone Zaggia

Abstract: Reliable models of the thermally pulsing asymptotic giant branch (TP-AGB) phase are of critical importance across astrophysics, including our interpretation of the spectral energy distribution of galaxies, cosmic dust production, and enrichment of the interstellar medium. With the aim of improving sets of stellar isochrones that include a detailed description of the TP-AGB phase, we extend our rec… ▽ More Reliable models of the thermally pulsing asymptotic giant branch (TP-AGB) phase are of critical importance across astrophysics, including our interpretation of the spectral energy distribution of galaxies, cosmic dust production, and enrichment of the interstellar medium. With the aim of improving sets of stellar isochrones that include a detailed description of the TP-AGB phase, we extend our recent calibration of the AGB population in the Small Magellanic Cloud (SMC) to the more metal rich Large Magellanic Cloud (LMC). We model the LMC stellar populations with the TRILEGAL code, using the spatially-resolved star formation history derived from the VISTA survey. We characterize the efficiency of the third dredge-up by matching the star counts and the $K_{\rm s}$-band luminosity functions of the AGB stars identified in the LMC. In line with previous findings, we confirm that, compared to the SMC, the third dredge-up in AGB stars of the LMC is somewhat less efficient, as a consequence of the higher metallicity. The predicted range of initial mass of C-rich stars is between $M_{\rm i} \approx 1.7 - 3~\mathrm{M}_{\odot}$ at $Z_{\rm i} = 0.008$. We show how the inclusion of new opacity data in the carbon star spectra will improve the performance of our models. We discuss the predicted lifetimes, integrated luminosities and mass-loss rate distributions of the calibrated models. The results of our calibration are included in updated stellar isochrones publicly available. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS

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

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

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

Submitted 17 August, 2020; originally announced August 2020.

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

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

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

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

Submitted 15 June, 2020; originally announced June 2020.

Comments: 32 pages, 21 figures, accepted to ApJ

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

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

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

Submitted 24 January, 2020; originally announced January 2020.

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

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

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

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

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

Comments: Accepted to AJ

Stellar Feedback and Resolved Stellar IFU Spectroscopy in the nearby Spiral Galaxy NGC 300

Authors: Anna F. McLeod, J. M. Diederik Kruijssen, Daniel R. Weisz, Peter Zeidler, Andreas Schruba, Julianne J. Dalcanton, Steven N. Longmore, Mélanie Chevance, Christopher M. Faesi, Nell Byler

Abstract: We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from… ▽ More We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from ground-based data at the distance of 2 Mpc. For the two star-forming complexes, in which no O-type stars had previously been identified, we find a total of 13 newly identified O-type stars in the mass range 15-50 M$_{\odot}$, as well as 4 Wolf-Rayet stars. We use the derived massive stellar content to analyze the impact of stellar feedback on the HII regions. As already found for HII regions in the Magellanic Clouds, the dynamics of the analyzed NGC 300 HII regions are dominated by a combination of the pressure of the ionized gas and stellar winds. By comparing the derived ionized gas mass loading factors to the total gas mass loading factor across the NGC 300 disk, we find that the latter is an order of magnitude higher, either indicating very early evolutionary stages for these HII regions, or being a direct result of the multi-phase nature of feedback-driven bubbles. Moreover, we analyze the relation between the star formation rate and the pressure of the ionized gas as derived from small (<100 pc) scales, as both quantities are systematically overestimated when derived on galactic scales. With the wealth of ongoing and upcoming IFU instruments and programs, this study serves as a pathfinder for the systematic investigation of resolved stellar feedback in nearby galaxies, and it delivers the necessary analysis tools to enable massive stellar content and feedback studies sampling an unprecedented range of HII region properties across entire galaxies in the nearby Universe. △ Less

Submitted 16 January, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: Accepted for publication in ApJ

Star formation at the edge of the Local Group: a rising star formation history in the isolated galaxy WLM

Authors: Saundra M. Albers, Daniel R. Weisz, Andrew A. Cole, Andrew E. Dolphin, Evan D. Skillman, Benjamin F. Williams, Michael Boylan-Kolchin, James S. Bullock, Julianne J. Dalcanton, Philip F. Hopkins, Ryan Leaman, Alan W. McConnachie, Mark Vogelsberger, Andrew Wetzel

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 13 pages, 13 Figures, 4 Tables. Re-submitted to MNRAS after addressing the referee's comments

Combined Effects of Rotation and Age Spreads on Extended Main Sequence Turn Offs

Authors: Seth Gossage, Charlie Conroy, Aaron Dotter, Ivan Cabrera-Ziri, Andrew E. Dolphin, Nate Bastian, Julianne J. Dalcanton, Paul Goudfrooij, L. Clifton Johnson, Benjamin F. Williams, Philip Rosenfield, Jason Kalirai, Morgan Fouesneau

Abstract: The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar… ▽ More The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar models to simulate observed Hubble Space Telescope and Gaia star cluster data. We model the effect of rotation as a non-parametric distribution, allowing for maximum flexibility. In our models the slow rotators comprise the blueward, and fast rotators the redward portion of the eMSTO. We simulate data under three scenarios: non-rotating eSF, a range of rotation rates with a single age, and a combination of age and rotation effects. We find that two of the five clusters (the youngest and oldest) favor an age spread, but these also achieve the overall worst fits of all clusters. The other three clusters show comparable statistical evidence between rotation and an age spread. In all five cases, a rotation rate distribution alone is capable of qualitatively matching the observed eMSTO structure. In future work, we aim to compare our predicted Vsin(i) with observations in order to better constrain the physics related to stellar rotation. △ Less

Submitted 18 November, 2019; v1 submitted 25 July, 2019; originally announced July 2019.

Comments: 17 pages, 11 figures

Fast and inefficient star formation due to short-lived molecular clouds and rapid feedback

Authors: J. M. Diederik Kruijssen, Andreas Schruba, Mélanie Chevance, Steven N. Longmore, Alexander P. S. Hygate, Daniel T. Haydon, Anna F. McLeod, Julianne J. Dalcanton, Linda J. Tacconi, Ewine F. van Dishoeck

Abstract: The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in… ▽ More The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in modern galaxy formation simulations. The key question is why the timescale for depleting molecular gas through star formation in galaxies ($t_{\rm dep}\approx 2$ Gyr) exceeds the dynamical timescale of GMCs by two orders of magnitude. Either most of a GMC's mass is converted into stars over many dynamical times, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale. Here we report our observation that molecular gas and star formation are spatially de-correlated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2-3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution. △ Less

Submitted 23 May, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: Published in the May 23 issue of Nature; this is the authors' version before final edits; low-resolution versions of the Supplementary Videos are available as ancillary files in the arXiv source (see the link in the top right of this page), whereas the full-resolution versions are available with the Nature publication

Journal ref: Nature 569, 519 (2019)

Self-consistent predictions for LIER-like emission lines from post-AGB stars

Authors: Nell Byler, Julianne J. Dalcanton, Charlie Conroy, Benjamin D. Johnson, Jieun Choi, Aaron Dotter, Philip Rosenfield

Abstract: Early type galaxies (ETGs) frequently show emission from warm ionized gas. These Low Ionization Emission Regions (LIERs) were originally attributed to a central, low-luminosity active galactic nuclei. However, the recent discovery of spatially-extended LIER emission suggests ionization by both a central source and an extended component that follows a stellar-like radial distribution. For passively… ▽ More Early type galaxies (ETGs) frequently show emission from warm ionized gas. These Low Ionization Emission Regions (LIERs) were originally attributed to a central, low-luminosity active galactic nuclei. However, the recent discovery of spatially-extended LIER emission suggests ionization by both a central source and an extended component that follows a stellar-like radial distribution. For passively-evolving galaxies with old stellar populations, hot post-Asymptotic Giant Branch (AGB) stars are the only viable extended source of ionizing photons. In this work, we present the first prediction of LIER-like emission from post-AGB stars that is based on fully self-consistent stellar evolution and photoionization models. We show that models where post-AGB stars are the dominant source of ionizing photons reproduce the nebular emission signatures observed in ETGs, including LIER-like emission line ratios in standard optical diagnostic diagrams and H$α$ equivalent widths of order 0.1-3 angstroms. We test the sensitivity of LIER-like emission to the details of post-AGB models, including the mass loss efficiency and convective mixing efficiency, and show that line strengths are relatively insensitive to post-AGB timescale variations. Finally, we examine the UV-optical colors of the models and the stellar populations responsible for the UV-excess observed in some ETGs. We find that allowing as little as 3% of the HB population to be uniformly distributed to very hot temperatures (30,000 K) produces realistic UV colors for old, quiescent ETGs. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: ApJ accepted. 20 pages, 8 figures

The Masses of Supernova Remnant Progenitors in M83

Authors: Benjamin F. Williams, Tristan J. Hillis, William P. Blair, Knox S. Long, Jeremiah W. Murphy, Andrew Dolphin, Rubab Khan, Julianne J. Dalcanton

Abstract: We determine the ages of the young, resolved stellar populations at the locations of 237 optically-identified supernova remnants in M83. These age distributions put constraints on the progenitor masses of the supernovae that produced 199 of the remnants. The other 38 show no evidence for having a young progenitor and are therefore good Type Ia SNR candidates. Starting from Hubble Space Telescope b… ▽ More We determine the ages of the young, resolved stellar populations at the locations of 237 optically-identified supernova remnants in M83. These age distributions put constraints on the progenitor masses of the supernovae that produced 199 of the remnants. The other 38 show no evidence for having a young progenitor and are therefore good Type Ia SNR candidates. Starting from Hubble Space Telescope broadband imaging, we measured resolved stellar photometry of seven archival WFC3/UVIS fields in F336W, F438W, and F814W. We generate color-magnitude diagrams of the stars within 50 pc of each SNR and fit them with stellar evolution models to obtain the population ages. From these ages we infer the progenitor mass that corresponds to the lifetime of the most prominent age that is $<$50 Myr. In this sample, there are 47 SNRs with best-fit progenitor masses $>$15 M$_{\odot}$, and 5 of these are $>$15 M$_{\odot}$ at 84% confidence. This is the largest collection of high-mass progenitors to date, including our highest-mass progenitor inference found so far, with a constraint of $<$8 Myr. Overall, the distribution of progenitor masses has a power-law index of $-3.0^{+0.2}_{-0.7}$, steeper than Salpeter initial mass function ($-2.35$). It remains unclear whether the reason for the low number of high-mass progenitors is due to the difficulty of finding and measuring such objects or because only a fraction of very massive stars produce supernovae. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: Accepted for publication in ApJ, 33 pages, 7 figures, 3 tables

A Dramatic Decrease in Carbon Star Formation in M31

Authors: M. L. Boyer, B. F. Williams, B. Aringer, Y. Chen, J. J. Dalcanton, L. Girardi, P. Guhathakurta, P. Marigo, K. A. G. Olsen, P. Rosenfield, D. R. Weisz

Abstract: We analyze resolved stellar near-infrared photometry of 21 HST fields in M31 to constrain the impact of metallicity on the formation of carbon stars. Observations of nearby galaxies show that the carbon stars are increasingly rare at higher metallicity. Models indicate that carbon star formation efficiency drops due to the decrease in dredge-up efficiency in metal-rich thermally-pulsing Asymptotic… ▽ More We analyze resolved stellar near-infrared photometry of 21 HST fields in M31 to constrain the impact of metallicity on the formation of carbon stars. Observations of nearby galaxies show that the carbon stars are increasingly rare at higher metallicity. Models indicate that carbon star formation efficiency drops due to the decrease in dredge-up efficiency in metal-rich thermally-pulsing Asymptotic Giant Branch (TP-AGB) stars, coupled to a higher initial abundance of oxygen. However, while models predict a metallicity ceiling above which carbon stars cannot form, previous observations have not yet pinpointed this limit. Our new observations reliably separate carbon stars from M-type TP-AGB stars across 2.6-13.7 kpc of M31's metal-rich disk using HST WFC3/IR medium-band filters. We find that the ratio of C to M stars (C/M) decreases more rapidly than extrapolations of observations in more metal-poor galaxies, resulting in a C/M that is too low by more than a factor of 10 in the innermost fields and indicating a dramatic decline in C star formation efficiency at metallicities higher than [M/H] $\approx$ -0.1 dex. The metallicity ceiling remains undetected, but must occur at metallicities higher than what is measured in M31's inner disk ([M/H] $\gtrsim$ +0.06 dex). △ Less

Submitted 24 May, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 16 pages, 13 Figures; text clarifications in response to the referee. Results are unchanged; accepted for publication in ApJ

Constraining the thermally-pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud

Authors: Giada Pastorelli, Paola Marigo, Léo Girardi, Yang Chen, Stefano Rubele, Michele Trabucchi, Bernhard Aringer, Sara Bladh, Alessandro Bressan, Josefina Montalbán, Martha L. Boyer, Julianne J. Dalcanton, Kjell Eriksson, Martin A. T. Groenewegen, Susanne Höfner, Thomas Lebzelter, Ambra Nanni, Philip Rosenfield, Peter R. Wood, Maria-Rosa L. Cioni

Abstract: The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed wi… ▽ More The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed with the TRILEGAL code. The strength of our approach relies on the detailed spatially-resolved star formation history of the SMC, derived from the deep near-infrared photometry of the VISTA survey of the Magellanic Clouds, as well as on the capability to quickly and accurately explore a wide variety of parameters and effects with the COLIBRI code for the TP-AGB evolution. Adopting a well-characterized set of observations -- star counts and luminosity functions -- we set up a calibration cycle along which we iteratively change a few key parameters of the TP-AGB models until we eventually reach a good fit to the observations. Our work leads to identify two best-fitting models that mainly differ in the efficiencies of the third dredge-up and mass loss in TP-AGB stars with initial masses larger than about 3 M$_{\odot}$. On the basis of these calibrated models we provide a full characterization of the TP-AGB stellar population in the SMC in terms of stellar parameters (initial masses, C/O ratios, carbon excess, mass-loss rates). Extensive tables of isochrones including these improved models are publicly available. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: accepted for publication in MNRAS. 29 pages, 24 figures, 2 appendices. One online appendix available as ancillary file on this page

Asymmetric Drift in the Andromeda Galaxy (M31) as a Function of Stellar Age

Authors: Amanda C. N. Quirk, Puragra Guhathakurta, Laurent Chemin, Claire E. Dorman, Karoline M. Gilbert, Anil C. Seth, Benjamin F. Williams, Julianne J. Dalcanton

Abstract: We analyze the kinematics of Andromeda's disk as a function of stellar age by using photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and spectroscopy from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We use HI 21-cm and CO ($\rm J=1 \rightarrow 0$) data to examine the difference between the deprojected rotation velocity of the gas… ▽ More We analyze the kinematics of Andromeda's disk as a function of stellar age by using photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and spectroscopy from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We use HI 21-cm and CO ($\rm J=1 \rightarrow 0$) data to examine the difference between the deprojected rotation velocity of the gas and that of the stars. We divide the stars into four stellar age bins, from shortest lived to longest lived: massive main sequence stars (0.03 Gyr), more luminous intermediate mass asymptotic giant branch (AGB) stars (0.4 Gyr), less luminous intermediate mass AGB stars (2 Gyr), and low mass red giant branch stars (4 Gyr). There is a clear correlation between the offset of the stellar and the gas rotation velocity, or the asymmetric drift: the longer lived populations lag farther behind the gas than short lived populations. We also examine possible causes of the substructure in the rotation curves and find that the most significant cause of scatter in the rotation curves comes from the tilted ring model being an imperfect way to account for the multiple warps in Andromeda's disk. △ Less

Submitted 20 November, 2018; v1 submitted 16 November, 2018; originally announced November 2018.

Comments: 30 pages, 24 figures, 2 tables; accepted to ApJ November 16, 2018

Journal ref: The Astrophysical Journal, 871:11 (20pp), 2019 January 20

Spatially Resolved Metal Loss from M31

Authors: O. Grace Telford, Jessica K. Werk, Julianne J. Dalcanton, Benjamin F. Williams

Abstract: As galaxies evolve, they must enrich and exchange gas with the surrounding medium, but the timing of these processes and how much gas is involved remain poorly understood. In this work, we leverage metals as tracers of past gas flows to constrain the history of metal ejection and redistribution in M31. This roughly $L*$ galaxy is a unique case where spatially resolved measurements of the gas-phase… ▽ More As galaxies evolve, they must enrich and exchange gas with the surrounding medium, but the timing of these processes and how much gas is involved remain poorly understood. In this work, we leverage metals as tracers of past gas flows to constrain the history of metal ejection and redistribution in M31. This roughly $L*$ galaxy is a unique case where spatially resolved measurements of the gas-phase and stellar metallicity, dust extinction, and neutral ISM gas content are all available, enabling a census of the present metal mass. We combine spatially resolved star formation histories from the Panchromatic Hubble Andromeda Treasury survey with a metal production model to calculate the history of metal production in M31. We find that $1.8\times10^9 \, M_\odot$ of metals, or 62\% of the metal mass formed within $r < 19 \,\mathrm{kpc}$, is missing from the disk in our fiducial model, implying that the M31 disk has experienced significant gaseous outflows over its lifetime. Under a conservative range of model assumptions, we find that between 3\% and 88\% of metals have been lost ($1.9\times10^7 - 6.4\times10^9 \, M_\odot$), so metals are missing even when all model parameters are chosen to favor metal retention. We show that the missing metal mass could be harbored in M31's CGM if the majority of the metals reside in a hot gas phase. Finally, we find that some metal mass produced in the past 1.5 Gyr in the central $\sim5 \,\mathrm{kpc}$ has likely been redistributed to larger radii within the disk. △ Less

Submitted 20 May, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: In press at ApJ; updated to match published version. Figure 7 shows key result

Comparing Chandra and Hubble in the Northern Disk of M31

Authors: Benjamin F. Williams, Margaret Lazzarini, Paul Plucinsky, Manami Sasaki, Vallia Antoniou, Neven Vulic, Michael Eracleous, Knox S. Long, Breanna Binder, Julianne J. Dalcanton, Alexia R. Lewis, Daniel R. Weisz

Abstract: The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square… ▽ More The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: 24 pages, 12 figures, 7 tables, accepted for publication in ApJS

Constraints for the Progenitor Masses of Historic Core-Collapse Supernovae

Authors: Benjamin F. Williams, Tristan J. Hillis, Jeremiah W. Murphy, Karoline Gilbert, Julianne J. Dalcanton, Andrew E. Dolphin

Abstract: We age-date the stellar populations associated with 12 historic nearby core-collapse supernovae (CCSNe) and 2 supernova impostors, and from these ages, we infer their initial masses and associated uncertainties. To do this, we have obtained new HST imaging covering these CCSNe. Using these images, we measure resolved stellar photometry for the stars surrounding the locations of the SNe. We then fi… ▽ More We age-date the stellar populations associated with 12 historic nearby core-collapse supernovae (CCSNe) and 2 supernova impostors, and from these ages, we infer their initial masses and associated uncertainties. To do this, we have obtained new HST imaging covering these CCSNe. Using these images, we measure resolved stellar photometry for the stars surrounding the locations of the SNe. We then fit the color-magnitude distributions of this photometry with stellar evolution models to determine the ages of any young existing populations present. From these age distributions, we infer the most likely progenitor mass for all of the SNe in our sample. We find ages between 4 and 50 Myr, corresponding to masses from 7.5 to 59 solar masses. There were no SNe that lacked a young population within 50~pc. Our sample contains 4 type Ib/c SNe; their masses have a wide range of values, suggesting that the progenitors of stripped-envelope SNe are binary systems. Both impostors have masses constrained to be $\lesssim$7.5 solar masses. In cases with precursor imaging measurements, we find that age-dating and precursor imaging give consistent progenitor masses. This consistency implies that, although the uncertainties for each technique are significantly different, the results of both are reliable to the measured uncertainties. We combine these new measurements with those from our previous work and find that the distribution of 25 core-collapse SNe progenitor masses is consistent with a standard Salpeter power-law mass function, no upper mass cutoff, and an assumed minimum mass for core-collapse of 7.5~M$_{\odot}$. △ Less

Submitted 21 March, 2018; originally announced March 2018.

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

Progenitor Mass Distribution for Core-Collapse Supernova Remnants in M31 and M33

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

Abstract: Using the star formation histories (SFHs) near 94 supernova remnants (SNRs), we infer the progenitor mass distribution for core-collapse supernovae. We use Bayesian inference and model each SFH with multiple bursts of star formation (SF), one of which is assumed to be associated with the SNR. Assuming single-star evolution, the minimum mass of CCSNe is $7.33^{+0.02}_{-0.16}$ $\text{M}_\odot$, the… ▽ More Using the star formation histories (SFHs) near 94 supernova remnants (SNRs), we infer the progenitor mass distribution for core-collapse supernovae. We use Bayesian inference and model each SFH with multiple bursts of star formation (SF), one of which is assumed to be associated with the SNR. Assuming single-star evolution, the minimum mass of CCSNe is $7.33^{+0.02}_{-0.16}$ $\text{M}_\odot$, the slope of the progenitor mass distribution is $α= -2.96^{+0.45}_{-0.25}$, and the maximum mass is greater than $\text{M}_\textrm{max} > 59$ $\text{M}_\odot$ with a 68% confidence. While these results are consistent with previous inferences, they also provide tighter constraints. The progenitor distribution is somewhat steeper than a Salpeter initial mass function ($α$ = -2.35). This suggests that either SNR catalogs are biased against the youngest SF regions, or the most massive stars do not explode as easily as lower mass stars. If SNR catalogs are biased, it will most likely affect the slope but not the minimum mass. The uncertainties are dominated by three primary sources of uncertainty, the SFH resolution, the number of SF bursts, and the uncertainty on SF rate in each age bin. We address the first two of these uncertainties, with an emphasis on multiple bursts. The third will be addressed in future work. △ Less

Submitted 29 August, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 15 pages, 8 figures

Panchromatic Hubble Andromeda Treasury. XIX. The Ancient Star Formation History of the M31 Disk

Authors: Benjamin F. Williams, Andrew E. Dolphin, Julianne J. Dalcanton, Daniel R. Weisz, Eric F. Bell, Alexia R. Lewis, Philip Rosenfield, Yumi Choi, Evan Skillman, Antonela Monachesi

Abstract: We map the star formation history across M31 by fitting stellar evolution models to color-magnitude diagrams of each 83$''{\times}$83$''$ (0.3$\times$1.4 kpc, deprojected) region of the PHAT survey outside of the innermost 6$'{\times}$12$'$ portion. We find that most of the star formation occurred prior to $\sim$8 Gyr ago, followed by a relatively quiescent period until $\sim$4 Gyr ago, a subseque… ▽ More We map the star formation history across M31 by fitting stellar evolution models to color-magnitude diagrams of each 83$''{\times}$83$''$ (0.3$\times$1.4 kpc, deprojected) region of the PHAT survey outside of the innermost 6$'{\times}$12$'$ portion. We find that most of the star formation occurred prior to $\sim$8 Gyr ago, followed by a relatively quiescent period until $\sim$4 Gyr ago, a subsequent star formation episode about 2 Gyr ago and a return to relative quiescence. There appears to be little, if any, structure visible for populations with ages older than 2 Gyr, suggesting significant mixing since that epoch. Finally, assuming a Kroupa IMF from 0.1$-$100 M$_{\odot}$, we find that the total amount of star formation over the past 14 Gyr in the area over which we have fit models is 5${\times}$10$^{10}$ M$_{\odot}$. Fitting the radial distribution of this star formation and assuming azimuthal symmetry, (1.5$\pm$0.2)${\times}$10$^{11}$ M$_{\odot}$ of stars have formed in the M31 disk as a whole, (9$\pm$2)${\times}$10$^{10}$ M$_{\odot}$ of which has likely survived to the present after accounting for evolutionary effects. This mass is about one fifth of the total dynamical mass of M31. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: 49 pages, 12 tables, 30 figures, accepted for publication in ApJ

The Spectroscopy and H-band Imaging of Virgo cluster galaxies (SHIVir) Survey: Scaling Relations and the Stellar-to-Total Mass Relation

Authors: Nathalie N. -Q. Ouellette, Stéphane Courteau, Jon A. Holtzman, Aaron A. Dutton, Michele Cappellari, Julianne J. Dalcanton, Michael McDonald, Joel C. Roediger, James E. Taylor, R. Brent Tully, Patrick Côté, Laura Ferrarese, Eric W. Peng

Abstract: We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tu… ▽ More We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Accepted to ApJ. 28 pages, 16 figures