Showing 1–50 of 105 results for author: Kalirai, J S

The structure of the stellar halo of the Andromeda galaxy explored with the NB515 for Subaru/HSC. I.: New Insights on the stellar halo up to 120 kpc

Authors: Itsuki Ogami, Mikito Tanaka, Yutaka Komiyama, Masashi Chiba, Puragra Guhathakurta, Evan N. Kirby, Rosemary F. G. Wyse, Carrie Filion, Karoline M. Gilbert, Ivanna Escala, Masao Mori, Takanobu Kirihara, Masayuki Tanaka, Miho N. Ishigaki, Kohei Hayashi, Myun Gyoon Lee, Sanjib Sharma, Jason S. Kalirai, Robert H. Lupton

Abstract: We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and CFHT/MegaCam g- & i-bands. We succeed in separating M31's halo stars from foreground contamination with $\sim$ 90 \% accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with… ▽ More We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and CFHT/MegaCam g- & i-bands. We succeed in separating M31's halo stars from foreground contamination with $\sim$ 90 \% accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with the Giant Southern Stream (GSS) based on their photometric metallicity estimates. We also produce the distance and photometric metallicity estimates for the known substructures. While these quantities for the GSS are reproduced in our study, we find that the North-Western stream shows a steeper distance gradient than found in an earlier study, suggesting that it is likely to have formed in an orbit closer to the Milky Way. For two streams in the eastern halo (Stream C and D), we identify distance gradients that had not been resolved. Finally, we investigate the global halo photometric metallicity distribution and surface brightness profile using the NB515-selected halo stars. We find that the surface brightness of the metal-poor and metal-rich halo populations, and the all population can be fitted to a power-law profile with an index of $α= -1.65 \pm 0.02$, $-2.82\pm0.01$, and $-2.44\pm0.01$, respectively. In contrast to the relative smoothness of the halo profile, its photometric metallicity distribution appears to be spatially non-uniform with nonmonotonic trends with radius, suggesting that the halo population had insufficient time to dynamically homogenize the accreted populations. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: 24 pages, 26 figures, 5 tables, submitted to MNRAS

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

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

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

Submitted 5 December, 2023; originally announced December 2023.

Comments: 27 pages, 10 figures, 1 table

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

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

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

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

Submitted 6 September, 2022; originally announced September 2022.

Comments: Accepted for publication in RNAAS

Kinematics and Metallicity of Red Giant Branch Stars in the Northeast Shelf of M31

Authors: Ivanna Escala, Karoline M. Gilbert, Mark Fardal, Puragra Guhathakurta, Robyn E. Sanderson, Jason S. Kalirai, Bahram Mobasher

Abstract: We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge… ▽ More We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge pattern agrees with expectations of a tidal shell formed in a radial merger and provides strong evidence in favor of predictions of Giant Stellar Stream (GSS) formation models in which the NE shelf originates from the second orbital wrap of the tidal debris. The observed concentric wedge patterns of the NE, West (W), and Southeast (SE) shelves corroborate this interpretation independently of the models. We do not detect a kinematical signature in the NE shelf region corresponding to an intact progenitor core, favoring GSS formation models in which the progenitor is completely disrupted. The shelf's photometric metallicity distribution implies that it is dominated by tidal material, as opposed to the phase-mixed stellar halo or the disk. The metallicity distribution ([Fe/H]$_{\rm phot}$ = $-0.42$ $\pm$ $0.01$) also matches the GSS, and consequently the W and SE shelves, further supporting a direct physical association between the tidal features. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 21 pages main text, 15 figures, 5 tables (including Appendix). Submitted to AJ

Relative Ages of Nine Inner Milky Way Globular Clusters from Proper Motion Cleaned Color-Magnitude Diagrams

Authors: Roger E. Cohen, Andrea Bellini, Luca Casagrande, Thomas M. Brown, Matteo Correnti, Jason S. Kalirai

Abstract: Our picture of the age-metallicity relation for Milky Way globular clusters (MWGCs) is still highly incomplete, and the majority of MWGCs lack self-consistent age measurements. Here, we exploit deep, homogenous multi-epoch Hubble Space Telescope (HST) imaging of nine MWGCs located towards the inner Milky Way to measure their relative ages, in most cases for the first time. Our relative age measure… ▽ More Our picture of the age-metallicity relation for Milky Way globular clusters (MWGCs) is still highly incomplete, and the majority of MWGCs lack self-consistent age measurements. Here, we exploit deep, homogenous multi-epoch Hubble Space Telescope (HST) imaging of nine MWGCs located towards the inner Milky Way to measure their relative ages, in most cases for the first time. Our relative age measurements are designed to be directly comparable to the large set of MWGC ages presented by VandenBerg et al. (2013, V13), using identical filters, evolutionary models, and bolometric corrections, extended to the higher extinction values relevant to our target clusters. Adopting the V13 MWGC age scale, our relative age measurements imply that our target clusters are consistently very old, with a mean age of 12.9$\pm$0.4 Gyr, with the exception of the young metal-rich MWGC NGC 6342. We perform two tests to validate the precision of our methodology, and discuss the implications of our target cluster loci in the MWGC age-metallicity plane. In addition, we use our fully self-consistent bolometric corrections to assess the systematic impact of variations in the total-to-selective extinction ratio $R_{V}$ on relative age measurements. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: AJ Accepted. 15 pages, 5 figures, 3 tables

Structure and Internal Kinematics of Nine Inner Milky Way Globular Clusters

Authors: Roger E. Cohen, Andrea Bellini, Mattia Libralato, Matteo Correnti, Thomas M. Brown, Jason S. Kalirai

Abstract: This study constitutes part of a larger effort aimed at better characterizing the Galactic globular clusters (GGCs) located towards the inner Milky Way bulge and disk. Here, we focus on internal kinematics of nine GGCs, obtained from space-based imaging over time baselines of $>$9 years. We exploit multiple avenues to assess the dynamical state of the target GGCs, constructing radial profiles of p… ▽ More This study constitutes part of a larger effort aimed at better characterizing the Galactic globular clusters (GGCs) located towards the inner Milky Way bulge and disk. Here, we focus on internal kinematics of nine GGCs, obtained from space-based imaging over time baselines of $>$9 years. We exploit multiple avenues to assess the dynamical state of the target GGCs, constructing radial profiles of projected stellar density, proper motion dispersion, and anisotropy. We posit that two-thirds (6/9) of our target GGCs are in an advanced dynamical state, and are close to (or have recently undergone) core collapse, supported by at least two lines of evidence: First, we find relatively steep proper motion dispersion profiles, in accord with literature values for core-collapsed GGCs. Second, we find that our sample is, in the mean, isotropic even out to their half-light radii, although one of our target clusters (NGC 6380) is tangentially anisotropic at $>$1$σ$ beyond its half-light radius, in accord with theoretical predictions for clusters evolving in strong tidal fields. Our proper motion dispersion and anisotropy profiles are made publicly available. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: AJ in press

Project AMIGA: The Circumgalactic Medium of Andromeda

Authors: Nicolas Lehner, Samantha C. Berek, J. Christopher Howk, Bart P. Wakker, Jason Tumlinson, Edward B. Jenkins, J. Xavier Prochaska, Ramona Augustin, Suoqing Ji, Claude-Andre Faucher-Giguere, Zachary Hafen, Molly S. Peeples, Kat A. Barger, Michelle A. Berg, Rongmon Bordoloi, Thomas M. Brown, Andrew J. Fox, Karoline M. Gilbert, Puragra Guhathakurta, Jason S. Kalirai, Felix J. Lockman, John M. O'Meara, D. J. Pisano, Joseph Ribaudo, Jessica K. Werk

Abstract: Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and… ▽ More Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: Submitted to the Astrophysical Journal. Comments welcome

Cool white dwarfs as standards for infrared observations

Authors: Nicola Pietro Gentile Fusillo, Pier-Emmanuel Tremblay, Ralph C. Bohlin, Susana E. Deustua, Jason S. Kalirai

Abstract: In the era of modern digital sky surveys, uncertainties in the flux of stellar standards are commonly the dominant systematic error in photometric calibration and can often affect the results of higher-level experiments. The Hubble Space Telescope (HST) spectrophotometry, which is based on computed model atmospheres for three hot (Teff > 30,000 K) pure-hydrogen (DA) white dwarfs, is currently cons… ▽ More In the era of modern digital sky surveys, uncertainties in the flux of stellar standards are commonly the dominant systematic error in photometric calibration and can often affect the results of higher-level experiments. The Hubble Space Telescope (HST) spectrophotometry, which is based on computed model atmospheres for three hot (Teff > 30,000 K) pure-hydrogen (DA) white dwarfs, is currently considered the most reliable and internally consistent flux calibration. However many next generation facilities (e.g. Harmoni on E-ELT, Euclid and JWST) will focus on IR observations, a regime in which white dwarf calibration has not yet been robustly tested. Cool DA white dwarfs have energy distributions that peak close to the optical or near-IR, do not have shortcomings from UV metal line blanketing, and have a reasonably large sky density ($\simeq$ 4 deg$^{-2}$ at G < 20), making them, potentially, excellent calibrators. Here we present a pilot study based on STIS+WFC3 observations of two bright DA white dwarfs to test whether targets cooler than current hot primary standards (Teff < 20,000 K) are consistent with the HST flux scale. We also test the robustness of white dwarf models in the IR regime from an X-shooter analysis of Paschen lines and by cross-matching our previously derived Gaia white dwarf catalogue with observations obtained with 2MASS, UKIDSS, VHS, and WISE. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 12 pages, 8 figures, 2 tables. Accepted for publication in MNRAS. The best fitting models to the two white dwarfs analyzed are available from \<http://deneb.astro.warwick.ac.uk/phrhcc/WD_models.tar>

Elemental Abundances in M31: First Alpha and Iron Abundance Measurements in M31's Giant Stellar Stream

Authors: Karoline M. Gilbert, Evan N. Kirby, Ivanna Escala, Jennifer Wojno, Jason S. Kalirai, Puragra Guhathakurta

Abstract: We present the first measurements of [Fe/H] and [$α$/Fe] abundances, obtained using spectral synthesis modeling, for red giant branch stars in M31's giant stellar stream. The spectroscopic observations, obtained at a projected distance of 17 kpc from M31's center, yielded 61 stars with [Fe/H] measurements, including 21 stars with [$α$/Fe] measurements, from 112 targets identified as M31 stars. The… ▽ More We present the first measurements of [Fe/H] and [$α$/Fe] abundances, obtained using spectral synthesis modeling, for red giant branch stars in M31's giant stellar stream. The spectroscopic observations, obtained at a projected distance of 17 kpc from M31's center, yielded 61 stars with [Fe/H] measurements, including 21 stars with [$α$/Fe] measurements, from 112 targets identified as M31 stars. The [Fe/H] measurements confirm the expectation from photometric metallicity estimates that stars in this region of M31's halo are relatively metal-rich compared to stars in the MW's inner halo: more than half the stars in the field, including those not associated with kinematically identified substructure, have [Fe/H] abundances $> -1.0$. The stars in this field are $α$-enhanced at lower metallicities, while [$α$/Fe] decreases with increasing [Fe/H] above metallicities of [Fe/H] $\gtrsim -0.9$. Three kinematical components have been previously identified in this field: the giant stellar stream, a second kinematically cold feature of unknown origin, and M31's kinematically hot halo. We compare probabilistic [Fe/H] and [$α$/Fe] distribution functions for each of the components. The giant stellar stream and the second kinematically cold feature have very similar abundance distributions, while the halo component is more metal-poor. Although the current sample sizes are small, a comparison of the abundances of stars in the giant stellar stream field with abundances of M31 halo and dSph stars from the literature indicate that the progenitor of the stream was likely more massive, and experienced a higher efficiency of star formation, than M31's existing dSphs or the dEs NGC147 and NGC185. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 23 pages, 12 figures, accepted for publication in ApJ

A Novel Approach to Constrain Rotational Mixing & Convective-Core Overshoot in Stars Using the Initial-Final Mass Relation

Authors: Jeffrey D. Cummings, Jason S. Kalirai, Jieun Choi, Cyril Georgy, Pier-Emmanuel Tremblay, Enrico Ramirez-Ruiz

Abstract: The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge… ▽ More The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge-up efficiencies, and convective-core overshoot may help explain the differences between models and observations. These parameters play an important role at the lowest masses (M$_{\rm initial}$ $<$ 3 M$_\odot$). At higher masses, only convective-core overshoot meaningfully affects white dwarf mass, but alone it likely cannot explain the observed white dwarf masses nor why the IFMR scatter is larger than observational errors predict. These higher masses, however, are also where rotational mixing in main sequence stars begins to create more massive cores, and hence more massive white dwarfs. This rotational mixing also extends a star's lifetime, making faster rotating progenitors appear like less massive stars in their semi-empirical age analysis. Applying the observed range of young B-dwarf rotations to the MIST or SYCLIST rotational models demonstrates a marked improvement in reproducing both the observed IFMR data and its scatter. The incorporation of both rotation and efficient convective-core overshoot significantly improves the match with observations. This work shows that the IFMR provides a valuable observational constraint on how rotation and convective-core overshoot affect the core evolution of a star. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 7 pages, 5 figures, accepted in ApJ Letters

Project AMIGA: Distance and Metallicity Gradients Along Andromeda's Giant Southern Stream from the Red Clump

Authors: Roger E. Cohen, Jason S. Kalirai, Karoline M. Gilbert, Puragra Guhathakurta, Molly S. Peeples, Nicolas Lehner, Thomas M. Brown, Luciana Bianchi, Kathleen A. Barger, John M. O'Meara

Abstract: The Giant Southern Stream (GSS) of M31, a keystone signature of a major accretion event, yields crucial constraints on M31 formation and evolution models. Currently, our understanding of the GSS, in terms of both its geometry and its chemistry, results from either wide-field imaging probing only a few magnitudes below the red giant branch tip, or deep imaging or spectroscopy of isolated regions. H… ▽ More The Giant Southern Stream (GSS) of M31, a keystone signature of a major accretion event, yields crucial constraints on M31 formation and evolution models. Currently, our understanding of the GSS, in terms of both its geometry and its chemistry, results from either wide-field imaging probing only a few magnitudes below the red giant branch tip, or deep imaging or spectroscopy of isolated regions. Here, we take an alternative approach, using Hubble Space Telescope (HST) imaging to characterize the horizontal branch red clump (RC) using unbinned maximum likelihood fits to luminosity functions (LFs) from observed color-magnitude diagrams (CMDs). Comparing the RC mean magnitude across three fields at projected distances of 21, 52 and 80 kpc from M31, we find a line of sight distance gradient identical to recent literature measurements in fields along the core. We also find tentative evidence that the line of sight distance dispersion increases with projected distance from M31. Meanwhile, the metallicity in the 52 kpc field westward of the GSS core is at least as high as that in the 21 kpc GSS core field, and the peak colors of the RC in these two fields imply identical metallicities to within 0.2 dex. We discuss implications for distance and metallicity gradients both along and perpendicular to the GSS in the context of recent ground-based photometric and spectroscopic results, including evidence for a dropoff in metallicity moving westward from the GSS, as well as prospects for further constraining stellar populations in the vicinity of the GSS. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 30 pages, 14 figures. AJ accepted. Figure quality has been reduced to meet arXiv size limits

The White Dwarf Initial-Final Mass Relation for Progenitor Stars From 0.85 to 7.5 M$_\odot$

Authors: Jeffrey D. Cummings, Jason S. Kalirai, P. -E. Tremblay, Enrico Ramirez-Ruiz, Jieun Choi

Abstract: We present the initial-final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These reobserved white dwarfs have masses ranging from 0.72 to 0.97 M$_\odot$, with initial masses… ▽ More We present the initial-final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These reobserved white dwarfs have masses ranging from 0.72 to 0.97 M$_\odot$, with initial masses from 3.0 to 3.65 M$_\odot$, where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter ($σ$ = 0.06 M$_\odot$) and spans from progenitors of 0.85 to 7.5 M$_\odot$. Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses ($>$5.5 M$_\odot$). The non-linearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 M$_{\rm final}$/M$_{\rm initial}$) and higher masses (0.11 M$_{\rm final}$/M$_{\rm initial}$) that are broken up by a steep slope (0.19 M$_{\rm final}$/M$_{\rm initial}$) between progenitors from 2.85 to 3.6 M$_\odot$. This IFMR shows total stellar mass loss ranges from 33\% of M$_{\rm initial}$ at 0.83 M$_\odot$ to 83\% of M$_{\rm initial}$ at 7.5 M$_\odot$. Testing this total mass loss for dependence on progenitor metallicity, however, finds no detectable sensitivity across the moderate range of -0.15 $<$ [Fe/H] $<$ +0.15. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

Improved Main Sequence Turnoff Ages of Young Open Clusters: Multicolor UBV Techniques & the Challenges of Rotation

Authors: Jeffrey D. Cummings, Jason S. Kalirai

Abstract: Main sequence turnoff ages in young open clusters are complicated by turnoffs that are sparse, have high binarity fractions, can be affected by differential reddening, and typically include a number of peculiar stars. Furthermore, stellar rotation can have a significant effect on a star's photometry and evolutionary timescale. In this paper we analyze in 12 nearby open clusters, ranging from ages… ▽ More Main sequence turnoff ages in young open clusters are complicated by turnoffs that are sparse, have high binarity fractions, can be affected by differential reddening, and typically include a number of peculiar stars. Furthermore, stellar rotation can have a significant effect on a star's photometry and evolutionary timescale. In this paper we analyze in 12 nearby open clusters, ranging from ages of 50 Myr to 350 Myr, how broadband UBV color-color relations can be used to identify turnoff stars that are Be stars, blue stragglers, certain types of binaries, or those affected by differential reddening. This UBV color-color analysis also directly measures a cluster's E(B-V) and estimates its [Fe/H]. The turnoff stars unaffected by these peculiarities create a narrower and more clearly defined cluster turnoff. Using four common isochronal models, two of which consider rotation, we fit cluster parameters using these selected turnoff stars and the main sequence. Comparisons of the photometrically fit cluster distances to those based on Gaia DR2 parallaxes find that they are consistent for all clusters. For older (>100 Myr) clusters, like the Pleiades and the Hyades, comparisons to ages based on the lithium depletion boundary method finds that these cleaned turnoff ages agree to within ~10% for all four isochronal models. For younger clusters, however, only the Geneva models that consider rotation fit turnoff ages consistent with lithium-based ages, while the ages based on non-rotating isochrones quickly diverge to become 30% to 80% younger. This illustrates the importance of rotation for deriving ages in the youngest (<100 Myr) clusters. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: Accepted for publication by AJ, 18 pages, 12 figures

The Age of the Old Metal-Poor Globular Cluster NGC6397 Using WFC3/IR Photometry

Authors: Matteo Correnti, Mario Gennaro, Jason S. Kalirai, Roger E. Cohen, Thomas M. Brown

Abstract: Globular Clusters (GCs) in the Milky Way represent the ideal laboratory to establish the age of the oldest stellar populations and to measure the color-magnitude relation of stars. Infrared (IR) photometry of these objects provides a new opportunity to accomplish this task. In particular, at low stellar masses, the stellar main sequence (MS) in an IR color-magnitude diagram (CMD) exhibits a sharp… ▽ More Globular Clusters (GCs) in the Milky Way represent the ideal laboratory to establish the age of the oldest stellar populations and to measure the color-magnitude relation of stars. Infrared (IR) photometry of these objects provides a new opportunity to accomplish this task. In particular, at low stellar masses, the stellar main sequence (MS) in an IR color-magnitude diagram (CMD) exhibits a sharp "kink" (due to opacity effects in M dwarfs), such that lower mass and cooler dwarfs become bluer in the F110W - F160W color baseline and not redder. This inversion of the color-magnitude relation offers the possibility to fit GC properties using IR imaging, and to reduce their uncertainties. Here, we used the IR channel of the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain new, deep high-resolution photometry of the old metal-poor GC NGC6397. From the analysis of the GC CMD, we revealed below the MS "kink" the presence of two MSs with different chemical composition. We derived the cluster fiducial line and we compared it with a grid of isochrones over a large range of parameter space, allowing age, metallicity, distance and reddening to vary freely within reasonable selected ranges. We derived an age of 12.6 Gyr with a random uncertainty sigma ~ 0.7 Gyr. These results confirm that the analysis of the IR color-magnitude of stars provide a valuable tool to measure the GC ages and offers a new venue to determine their absolute age to sub-Gyr accuracy with next generation IR telescopes. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: 9 pages in emulateapj LaTeX style, incl. 9 figures, 2 tables. Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:1603.05254

Deep Hubble Space Telescope Imaging of Globular Clusters Towards the Galactic Bulge: Observations, Data Reduction and Color-Magnitude Diagrams

Authors: Roger E. Cohen, Francesco Mauro, Javier Alonso-García, Maren Hempel, Ata Sarajedini, Antonio J. Ordoñez, Douglas Geisler, Jason S. Kalirai

Abstract: The Galactic globular clusters (GGCs) located towards the Galactic bulge have generally been excluded from large-scale photometric GGC surveys due to severe total and differential extinction. Here, we present an overview of a Hubble Space Telescope program designed to obtain deep, high spatial resolution multiband imaging of 16 poorly studied GGCs located towards the inner Galactic bulge and disk.… ▽ More The Galactic globular clusters (GGCs) located towards the Galactic bulge have generally been excluded from large-scale photometric GGC surveys due to severe total and differential extinction. Here, we present an overview of a Hubble Space Telescope program designed to obtain deep, high spatial resolution multiband imaging of 16 poorly studied GGCs located towards the inner Galactic bulge and disk. In this first paper of a series resulting from these observations, we give an overview of target cluster selection, observations and data reduction procedures for optimizing the resulting photometric catalogs. Artificial star tests are used to compare the respective advantages of different data reduction strategies in terms of photometric and astrometric precision and photometric incompleteness. We present the resulting color-magnitude diagrams (CMDs) of all target clusters in several color-magnitude planes, along with CMDs of comparison fields from parallel observations. For each target cluster, we summarize existing studies, and discuss their CMDs qualitatively in the context of these results. △ Less

Submitted 28 May, 2018; originally announced May 2018.

Comments: AJ accepted. Figure quality has been reduced to meet arXiv size limits

On the use of the main sequence knee (saddle) to measure globular cluster ages

Authors: S. Saracino, E. Dalessandro, F. R. Ferraro, B. Lanzoni, L. Origlia, M. Salaris, A. Pietrinferni, D. Geisler, J. S. Kalirai, M. Correnti, R. E. Cohen, F. Mauro, S. Villanova, C. Moni Bidin

Abstract: In this paper we review the operational definition of the so-called main sequence knee (MS-knee), a feature in the color-magnitude diagram (CMD) occurring at the low-mass end of the MS. The magnitude of this feature is predicted to be independent of age at fixed chemical composition. For this reason, its difference in magnitude with respect to the MS turn-off (MS-TO) point has been suggested as a… ▽ More In this paper we review the operational definition of the so-called main sequence knee (MS-knee), a feature in the color-magnitude diagram (CMD) occurring at the low-mass end of the MS. The magnitude of this feature is predicted to be independent of age at fixed chemical composition. For this reason, its difference in magnitude with respect to the MS turn-off (MS-TO) point has been suggested as a possible diagnostic to estimate absolute globular cluster (GC) ages. We first demonstrate that the operational definition of the MS-knee currently adopted in the literature refers to the inflection point of the MS (that we here more appropriately named MS-saddle), a feature that is well distinct from the knee and that cannot be used as its proxy. The MS-knee is only visible in near-infrared CMDs, while the MS-saddle can be also detected in optical-NIR CMDs. By using different sets of isochrones we then demonstrate that the absolute magni- tude of the MS-knee varies by a few tenths of a dex from one model to another, thus showing that at the moment stellar models may not capture the full systematic error in the method. We also demonstrate that while the absolute magnitude of the MS-saddle is almost coincident in different models, it has a systematic dependence on the adopted color combinations which is not predicted by stellar models. Hence, it cannot be used as a reliable reference for absolute age determination. Moreover, when statistical and systematic uncertainties are properly taken into ac- count, the difference in magnitude between the MS-TO and the MS-saddle does not provide absolute ages with better accuracy than other methods like the MS-fitting. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 31 pages, 12 figures, 4 tables. Accepted for publication in The Astrophysical Journal

Evidence of a non universal stellar Initial Mass Function. Insights from HST optical imaging of 6 Ultra Faint Dwarf Milky Way Satellites

Authors: Mario Gennaro, Kirill Tchernyshyov, Thomas M. Brown, Marla Geha, Roberto J. Avila, Puragra Guhathakurta, Jason S. Kalirai, Evan N. Kirby, Alvio Renzini, Joshua D. Simon, Jason Tumlinson, Luis C. Vargas

Abstract: Using deep HST/ACS observations, we demonstrate that the sub-solar stellar initial mass function (IMF) of 6 ultra-faint dwarf Milky Way Satellites (UFDs) is more bottom light than the IMF of the Milky Way disk. Our data have a lower mass limit of about 0.45 M$_{\odot}$, while the upper limit is $\sim 0.8$ M$_\odot$, set by the turn-off mass of these old, metal poor systems. If formulated as a sing… ▽ More Using deep HST/ACS observations, we demonstrate that the sub-solar stellar initial mass function (IMF) of 6 ultra-faint dwarf Milky Way Satellites (UFDs) is more bottom light than the IMF of the Milky Way disk. Our data have a lower mass limit of about 0.45 M$_{\odot}$, while the upper limit is $\sim 0.8$ M$_\odot$, set by the turn-off mass of these old, metal poor systems. If formulated as a single power law, we obtain a shallower IMF slope than the "Salpeter" value of $-2.3$, ranging from $-1.01$ for Leo IV, to $-1.87$ for Boötes I. The significance of such deviations depends on the galaxy and is typically 95\% or more. When modeled as a log-normal, the IMF fit results in a larger peak mass than in the Milky Way disk, however a Milky Way disk value for the characteristic system mass ($\sim0.22$ M$_{\odot}$) is excluded only at 68\% significance, and only for some UFDs in the sample. We find that the IMF slope correlates well with the galaxy mean metallicity and, to a lesser degree, with the velocity dispersion and the total mass. The strength of the observed correlations is limited by shot noise in the number of observed stars, but future space-based missions like JWST and WFIRST will both enhance the number of dwarf Milky Way Satellites that can be studied in such detail, and the observation depth for individual galaxies. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: Accepted for publication in ApJ on 18 Jan 2018

Stellar Stream and Halo Structure in the Andromeda Galaxy From a Subaru/Hyper Suprime-Cam Survey

Authors: Yutaka Komiyama, Masashi Chiba, Mikito Tanaka, Masayuki Tanaka, Takanobu Kirihara, Yohei Miki, Masao Mori, Robert H. Lupton, Puragra Guhathakurta, Jason S. Kalirai, Karoline Gilbert, Evan Kirby, Myun Gyoon Lee, In Sung Jang, Sanjib Sharma, Kohei Hayashi

Abstract: We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the N… ▽ More We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 24 pages, 18 figures, accepted for publication in ApJ

Global Properties of M31's Stellar Halo from the SPLASH Survey: III. Measuring the Stellar Velocity Dispersion Profile

Authors: Karoline M. Gilbert, Erik Tollerud, Rachael L. Beaton, Puragra Guhathakurta, James S. Bullock, Masashi Chiba, Jason S. Kalirai, Evan N. Kirby, Steven R. Majewski, Mikito Tanaka

Abstract: We present the velocity dispersion of red giant branch (RGB) stars in M31's halo, derived by modeling the line of sight velocity distribution of over 5000 stars in 50 fields spread throughout M31's stellar halo. The dataset was obtained as part of the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) Survey, and covers projected radii of 9 to 175 kpc from M31's center. A… ▽ More We present the velocity dispersion of red giant branch (RGB) stars in M31's halo, derived by modeling the line of sight velocity distribution of over 5000 stars in 50 fields spread throughout M31's stellar halo. The dataset was obtained as part of the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) Survey, and covers projected radii of 9 to 175 kpc from M31's center. All major structural components along the line of sight in both the Milky Way (MW) and M31 are incorporated in a Gaussian Mixture Model, including all previously identified M31 tidal debris features in the observed fields. The probability an individual star is a constituent of M31 or the MW, based on a set of empirical photometric and spectroscopic diagnostics, is included as a prior probability in the mixture model. The velocity dispersion of stars in M31's halo is found to decrease only mildly with projected radius, from 108 km/s in the innermost radial bin (8.2 to 14.1 kpc) to $\sim 80$ to 90 km/s at projected radii of $\sim 40$ to 130 kpc, and can be parameterized with a power-law of slope $-0.12\pm 0.05$. The quoted uncertainty on the power-law slope reflects only the precision of the method, although other sources of uncertainty we consider contribute negligibly to the overall error budget. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: Submitted to the Astrophysical Journal

Can magnetic fields suppress convection in the atmosphere of cool white dwarfs? A case study on WD2105-820

Authors: Nicola Pietro Gentile Fusillo, Pier-Emmanuel Tremblay, Stefan Jordan, Boris T. Gänsicke, Jason S. Kalirai, Jeffrey Cummings

Abstract: Around 10% of white dwarfs exhibit global magnetic structures with fields ranging from 1 kG to hundreds of MG. Recently, the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarfs showed that convection should be suppressed in their photospheres for magnetic fields with strengths B $\gtrsim$ 50 kG. These predictions are in agreement with our knowledge of stellar physics… ▽ More Around 10% of white dwarfs exhibit global magnetic structures with fields ranging from 1 kG to hundreds of MG. Recently, the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarfs showed that convection should be suppressed in their photospheres for magnetic fields with strengths B $\gtrsim$ 50 kG. These predictions are in agreement with our knowledge of stellar physics (e.g. energy transfer in strong magnetic field regions of the solar photosphere), but have yet to be directly confirmed from white dwarf observations. We obtained COS far-UV spectroscopy of the weakly magnetic, hydrogen-atmosphere, white dwarf WD2105-820 and of three additional non-magnetic, convective remnants (all in the $T_{\mathrm{eff}}$ range 9000-11,000 K). We fitted both the COS and the already available optical spectra with convective and radiative atmospheric models. As expected, we find that for two of the non-magnetic comparison stars only convective model fits predict consistent $T_{\mathrm{eff}}$ values from both the optical and the FUV spectra. In contrast, for WD2105-820 only the best fitting radiative model produced consistent results. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: 8 pages, 7 figures, 1 table, accepted for publication in MNRAS

Dissecting the Extended Main Sequence Turn-off of the Young Star Cluster NGC1850

Authors: Matteo Correnti, Paul Goudfrooij, Andrea Bellini, Jason S. Kalirai, Thomas H. Puzia

Abstract: We use the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain deep, high-resolution photometry of the young (~ 100 Myr) star cluster NGC1850 in the Large Magellanic Cloud. We analyze the cluster colour-magnitude diagram (CMD) and find that it hosts an extended main sequence turn-off (MSTO) and a double MS. We demonstrate that these features cannot be due to photometric errors, field… ▽ More We use the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain deep, high-resolution photometry of the young (~ 100 Myr) star cluster NGC1850 in the Large Magellanic Cloud. We analyze the cluster colour-magnitude diagram (CMD) and find that it hosts an extended main sequence turn-off (MSTO) and a double MS. We demonstrate that these features cannot be due to photometric errors, field star contamination, or differential reddening. From a comparison with theoretical models and Monte Carlo simulations, we show that a coeval stellar population featuring a distribution of stellar rotation rates can reproduce the MS split quite well. However, it cannot reproduce the observed MSTO region, which is significantly wider than the simulated ones. Exploiting narrow-band Halpha imaging, we find that the MSTO hosts a population of Halpha-emitting stars which are interpreted as rapidly rotating Be-type stars. We explore the possibility that the discrepancy between the observed MSTO morphology and that of the simulated simple stellar population (SSP) is caused by the fraction of these objects that are highly reddened, but we rule out this hypothesis. We demonstrate that the global CMD morphology is well-reproduced by a combination of SSPs that cover an age range of ~ 35 Myr as well as a wide variety of rotation rates. We derive the cluster mass and escape velocity and use dynamical evolution models to predict their evolution starting at an age of 10 Myr. We discuss these results and their implications in the context of the extended MSTO phenomenon. △ Less

Submitted 27 December, 2016; originally announced December 2016.

Comments: 15 pages, 15, figures, 1 table. Accepted for publication in MNRAS

The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

Authors: Karl D. Gordon, Morgan Fouesneau, Heddy Arab, Kirill Tchernyshyov, Daniel R. Weisz, Julianne J. Dalcanton, Benjamin F. Williams, Eric F. Bell, Luciana Bianchi, Martha Boyer, Yumi Choi, Andrew Dolphin, Leo Girardi, David W. Hogg, Jason S. Kalirai, Maria Kapala, Alexia R. Lewis, Hans-Walter Rix, Karin Sandstrom, Evan D. Skillman

Abstract: We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the s… ▽ More We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys. △ Less

Submitted 20 June, 2016; originally announced June 2016.

Comments: 20 pages, 19 figures, ApJ, in press

The Field White Dwarf Mass Distribution

Authors: P. -E. Tremblay, J. Cummings, J. S. Kalirai, B. T. Gaensicke, N. Gentile-Fusillo, R. Raddi

Abstract: We revisit the properties and astrophysical implications of the field white dwarf mass distribution in preparation of Gaia applications. Our study is based on the two samples with the best established completeness and most precise atmospheric parameters, the volume-complete survey within 20 pc and the Sloan Digital Sky Survey (SDSS) magnitude-limited sample. We explore the modelling of the observe… ▽ More We revisit the properties and astrophysical implications of the field white dwarf mass distribution in preparation of Gaia applications. Our study is based on the two samples with the best established completeness and most precise atmospheric parameters, the volume-complete survey within 20 pc and the Sloan Digital Sky Survey (SDSS) magnitude-limited sample. We explore the modelling of the observed mass distributions with Monte Carlo simulations, but find that it is difficult to constrain independently the initial mass function (IMF), the initial-to-final-mass relation (IFMR), the stellar formation history (SFH), the variation of the Galactic disk vertical scale height as a function of stellar age, and binary evolution. Each of these input ingredients has a moderate effect on the predicted mass distributions, and we must also take into account biases owing to unidentified faint objects (20 pc sample), as well as unknown masses for magnetic white dwarfs and spectroscopic calibration issues (SDSS sample). Nevertheless, we find that fixed standard assumptions for the above parameters result in predicted mean masses that are in good qualitative agreement with the observed values. It suggests that derived masses for both studied samples are consistent with our current knowledge of stellar and Galactic evolution. Our simulations overpredict by 40-50% the number of massive white dwarfs (M > 0.75 Msun) for both surveys, although we can not exclude a Salpeter IMF when we account for all biases. Furthermore, we find no evidence of a population of double white dwarf mergers in the observed mass distributions. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 15 pages, 16 figures, accepted for publication in MNRAS

Constraining Globular Cluster Age Uncertainties Using the IR Color-Magnitude Diagram

Authors: Matteo Correnti, Mario Gennaro, Jason S. Kalirai, Thomas M. Brown, Annalisa Calamida

Abstract: Globular Clusters (GCs) in the Milky Way are the primary laboratories for establishing the ages of the oldest stellar populations and for measuring the color-magnitude relation of stars. In infrared (IR) color-magnitude diagrams (CMDs), the stellar main sequence (MS) exhibits a "kink", due to opacity effects in M dwarfs, such that lower mass and cooler dwarfs become bluer in the IR color baseline.… ▽ More Globular Clusters (GCs) in the Milky Way are the primary laboratories for establishing the ages of the oldest stellar populations and for measuring the color-magnitude relation of stars. In infrared (IR) color-magnitude diagrams (CMDs), the stellar main sequence (MS) exhibits a "kink", due to opacity effects in M dwarfs, such that lower mass and cooler dwarfs become bluer in the IR color baseline. This diagnostic offers a new opportunity to model GC CMDs and to reduce uncertainties on cluster properties (e.g., their derived ages). In this context, we analyzed Hubble Space Telescope Wide Field Camera 3 IR archival observations of four GCs - 47Tuc, M4, NGC2808, and NGC6752 - for which the data are deep enough to fully sample the low-mass MS, reaching at least ~ 2 mag below the "kink". We derived the fiducial lines for each cluster and compared them with a grid of isochrones over a large range of parameter space, allowing age, metallicity, distance, and reddening to vary within reasonable selected ranges. The derived ages for the four clusters are respectively 11.6, 11.5, 11.2, and 12.1 Gyr and their random uncertainties are σ~ 0.7 - 1.1 Gyr. Our results suggest that the near-IR MS "kink", combined with the MS turn-off, provides a valuable tool to measure GC ages and offers a promising opportunity to push the absolute age of GCs to sub-Gyr accuracy with the next generation IR telescopes such as the James Webb Space Telescope and the Wide-Field Infrared Survey Telescope. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: 16 pages in emulateapj LaTeX style, incl. 10 figures (low res.), 4 tables. Accepted for publication in ApJ

An Ultramassive 1.28 M$_\odot$ White Dwarf in NGC 2099

Authors: Jeffrey D. Cummings, Jason S. Kalirai, P. -E. Tremblay, Enrico Ramirez-Ruiz, P. Bergeron

Abstract: With the Keck I Low-Resolution Imaging Spectrometer we have observed nine white dwarf candidates in the very rich open cluster NGC 2099 (M37). The spectroscopy shows seven to be DA white dwarfs, one to be a DB white dwarf, and one to be a DZ white dwarf. Three of these DA white dwarfs are consistent with singly evolved cluster membership: an ultramassive (1.28$^{+0.05}_{-0.08}$ M$_\odot$) and two… ▽ More With the Keck I Low-Resolution Imaging Spectrometer we have observed nine white dwarf candidates in the very rich open cluster NGC 2099 (M37). The spectroscopy shows seven to be DA white dwarfs, one to be a DB white dwarf, and one to be a DZ white dwarf. Three of these DA white dwarfs are consistent with singly evolved cluster membership: an ultramassive (1.28$^{+0.05}_{-0.08}$ M$_\odot$) and two intermediate-mass (0.70 and 0.75 M$_\odot$) white dwarfs. Analysis of their cooling ages allows us to calculate their progenitor masses and establish new constraints on the initial-final mass relation. The intermediate-mass white dwarfs are in strong agreement with previous work over this mass regime. The ultramassive white dwarf has $V$ = 24.5, $\sim$2 mag fainter than the other two remnants. The spectrum of this star has lower quality, so the derived stellar properties (e.g., T$_{\rm eff}$, log g) have uncertainties that are several times higher than the brighter counterparts. We measure these uncertainties and establish the star's final mass as the highest-mass white dwarf discovered thus far in a cluster, but we are unable to calculate its progenitor mass because at this high mass and cooler T$_{\rm eff}$ its inferred cooling age is highly sensitive to its mass. At the highest temperatures, however, this sensitivity of cooling age to an ultramassive white dwarf's mass is only moderate. This demonstrates that future investigations of the upper-mass end of the initial-final mass relation must identify massive, newly formed white dwarfs (i.e., in young clusters with ages 50-150 Myr). △ Less

Submitted 1 March, 2016; originally announced March 2016.

Comments: 6 pages, 4 figures

Two Massive White Dwarfs from NGC 2323 and the Initial-Final Mass Relation for Progenitors of 4 to 6.5 M$_\odot$

Authors: Jeffrey D. Cummings, Jason S. Kalirai, P. -E. Tremblay, Enrico Ramirez-Ruiz

Abstract: We have observed a sample of 10 white dwarf candidates in the rich open cluster NGC 2323 (M50) with the Keck Low-Resolution Imaging Spectrometer. The spectroscopy shows eight to be DA white dwarfs, with six of these having high S/N appropriate for our analysis. Two of these white dwarfs are consistent with singly evolved cluster membership, and both are high mass $\sim$1.07 M$_\odot$, and give equ… ▽ More We have observed a sample of 10 white dwarf candidates in the rich open cluster NGC 2323 (M50) with the Keck Low-Resolution Imaging Spectrometer. The spectroscopy shows eight to be DA white dwarfs, with six of these having high S/N appropriate for our analysis. Two of these white dwarfs are consistent with singly evolved cluster membership, and both are high mass $\sim$1.07 M$_\odot$, and give equivalent progenitor masses of 4.69 M$_\odot$. To supplement these new high-mass white dwarfs and analyze the initial-final mass relation (IFMR), we have also looked at 30 white dwarfs from publicly available data that are mostly all high-mass ($\gtrsim$0.9 M$_\odot$). These original published data exhibited significant scatter, and to test if this scatter is true or simply the result of systematics, we have uniformly analyzed the white dwarf spectra and have adopted thorough photometric techniques to derive uniform cluster parameters for their parent clusters. The resulting IFMR scatter is significantly reduced, arguing that mass-loss rates are not stochastic in nature and that within the ranges of metallicity and mass analyzed in this work mass loss is not highly sensitive to variations in metallicity. Lastly, when adopting cluster ages based on Y$^2$ isochrones, the slope of the high-mass IFMR remains steep and consistent with that found from intermediate-mass white dwarfs, giving a linear IFMR from progenitor masses between 3 to 6.5 M$_\odot$. In contrast, when adopting the slightly younger cluster ages based on PARSEC isochrones, the high-mass IFMR has a moderate turnover near an initial mass of 4 M$_\odot$. △ Less

Submitted 12 January, 2016; originally announced January 2016.

Comments: 14 pages, 10 figures

Initial-Final Mass Relation for 3 to 4 M$_\odot$ Progenitors of White Dwarfs from the Single Cluster NGC 2099

Authors: Jeffrey D. Cummings, Jason S. Kalirai, P. -E. Tremblay, Enrico Ramirez-Ruiz

Abstract: We have expanded the sample of observed white dwarfs in the rich open cluster NGC 2099 (M37) with the Keck Low-Resolution Imaging Spectrometer. Of 20 white dwarf candidates, the spectroscopy shows 19 to be true white dwarfs with 14 of these having high S/N. We find 11 of these 14 to be consistent with singly evolved cluster members. They span a mass range of $\sim$0.7 to 0.95 M$_\odot$, excluding… ▽ More We have expanded the sample of observed white dwarfs in the rich open cluster NGC 2099 (M37) with the Keck Low-Resolution Imaging Spectrometer. Of 20 white dwarf candidates, the spectroscopy shows 19 to be true white dwarfs with 14 of these having high S/N. We find 11 of these 14 to be consistent with singly evolved cluster members. They span a mass range of $\sim$0.7 to 0.95 M$_\odot$, excluding a low-mass outlier, corresponding to progenitor masses of $\sim$3 to 4 M$_\odot$. This region of the initial final mass relation (IFMR) has large scatter and a slope that remains to be precisely determined. With this large sample of white dwarfs that belong to a single age and metallicity population, we find an initial-final mass relation of (0.171$\pm$0.057)M$_{\rm initial}$+0.219$\pm$0.187 M$_\odot$, significantly steeper than the linear relation adopted over the full observed white dwarf mass range in many previous studies. Comparison of this new relation from the solar metallicity NGC 2099 to 18 white dwarfs in the metal-rich Hyades and Praesepe shows that their IFMR also has a consistently steep slope. This strong consistency also suggests that there is no significant metallicity dependence of the IFMR at this mass and metallicity range. As a result, the IFMR can be more reliably determined with this broad sample of 29 total white dwarfs giving M$_{\rm final}$=(0.163$\pm$0.022)M$_{\rm initial}$+0.238$\pm$0.071 M$_\odot$ from M$_{\rm initial}$ of 3 to 4 M$_\odot$. A steep IFMR in this mass range indicates that the full IFMR is nonlinear. △ Less

Submitted 25 May, 2015; originally announced May 2015.

Comments: 13 pages, 6 figures, 3 tables

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

Authors: Alexia R. Lewis, Andrew E. Dolphin, Julianne J. Dalcanton, Daniel R. Weisz, Benjamin F. Williams, Eric F. Bell, Anil C. Seth, Jacob E. Simones, Evan D. Skillman, Yumi Choi, Morgan Fouesneau, Puragra Guhathakurta, Lent C. Johnson, Jason S. Kalirai, Adam K. Leroy, Antonela Monachesi, Hans-Walter Rix, Andreas Schruba

Abstract: We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that th… ▽ More We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged] △ Less

Submitted 13 April, 2015; originally announced April 2015.

Comments: 23 pages, 17 figures, 2 tables, accepted for publication in ApJ

The High-Mass Stellar Initial Mass Function in M31 Clusters

Authors: Daniel R. Weisz, L. Clifton Johnson, Daniel Foreman-Mackey, Andrew E. Dolphin, Lori C. Beerman, Benjamin F. Williams, Julianne J. Dalcanton, Hans-Walter Rix, David W. Hogg, Morgan Fouesneau, Benjamin D. Johnson, Eric F. Bell, Martha L. Boyer, Dimitrios Gouliermis, Puragra Guhathakurta, Jason S. Kalirai, Alexia R. Lewis, Anil C. Seth, Evan D. Skillman

Abstract: We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope fo… ▽ More We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope for stars >2 Msun. For the ensemble of clusters, the distribution of stellar MF slopes is best described by $Γ=+1.45^{+0.03}_{-0.06}$ with a very small intrinsic scatter. The data also imply no significant dependencies of the MF slope on cluster age, mass, and size, providing direct observational evidence that the measured MF represents the IMF. This analysis implies that the high-mass IMF slope in M31 clusters is universal with a slope ($Γ=+1.45^{+0.03}_{-0.06}$) that is steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values. Using our inference model on select Milky Way (MW) and LMC high-mass IMF studies from the literature, we find $Γ_{\rm MW} \sim+1.15\pm0.1$ and $Γ_{\rm LMC} \sim+1.3\pm0.1$, both with intrinsic scatter of ~0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be universal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study represents the most robust measurement of the high-mass IMF slope to date. We have grafted the M31 high-mass IMF slope onto widely used sub-solar mass Kroupa and Chabrier IMFs and show that commonly used UV- and Halpha-based star formation rates should be increased by a factor of ~1.3-1.5 and the number of stars with masses >8 Msun are ~25% fewer than expected for a Salpeter/Kroupa IMF. [abridged] △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: 11 pages, 7 Figures, submitted to ApJ. Comments welcome

Tracing the Metal-Poor M31 Stellar Halo with Blue Horizontal Branch Stars

Authors: Benjamin F. Williams, Julianne J. Dalcanton, Eric F. BellKaroline M. Gilbert, Puragra Guhathakurta, Claire Dorman, Tod R. Lauer, Anil C. Seth, Jason S. Kalirai, Philip Rosenfield, Leo Girardi

Abstract: We have analyzed new HST/ACS and HST/WFC3 imaging in F475W and F814W of two previously-unobserved fields along the M31 minor axis to confirm our previous constraints on the shape of M31's inner stellar halo. Both of these new datasets reach a depth of at least F814W$<$27 and clearly detect the blue horizontal branch (BHB) of the field as a distinct feature of the color-magnitude diagram. We measur… ▽ More We have analyzed new HST/ACS and HST/WFC3 imaging in F475W and F814W of two previously-unobserved fields along the M31 minor axis to confirm our previous constraints on the shape of M31's inner stellar halo. Both of these new datasets reach a depth of at least F814W$<$27 and clearly detect the blue horizontal branch (BHB) of the field as a distinct feature of the color-magnitude diagram. We measure the density of BHB stars and the ratio of BHB to red giant branch stars in each field using identical techniques to our previous work. We find excellent agreement with our previous measurement of a power-law for the 2-D projected surface density with an index of 2.6$^{+0.3}_{-0.2}$ outside of 3 kpc, which flattens to $α<$1.2 inside of 3 kpc. Our findings confirm our previous suggestion that the field BHB stars in M31 are part of the halo population. However, the total halo profile is now known to differ from this BHB profile, which suggests that we have isolated the metal-poor component. This component appears to have an unbroken power-law profile from 3-150 kpc but accounts for only about half of the total halo stellar mass. Discrepancies between the BHB density profile and other measurements of the inner halo are therefore likely due to the different profile of the metal-rich halo component, which is not only steeper than the profile of the met al-poor component, but also has a larger core radius. These profile differences also help to explain the large ratio of BHB/RGB stars in our observations. △ Less

Submitted 26 January, 2015; originally announced January 2015.

Comments: 11 pages, 6 figures, accepted for publication in ApJ

Extended Main Sequence Turnoffs in Intermediate-Age Star Clusters: A Correlation Between Turnoff Width and Early Escape Velocity

Authors: Paul Goudfrooij, Leo Girardi, Vera Kozhurina-Platais, Jason S. Kalirai, Imants Platais, Thomas H. Puzia, Matteo Correnti, Alessandro Bressan, Rupali Chandar, Leandro Kerber, Paola Marigo, Stefano Rubele

Abstract: We present color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1 - 2 Gyr old) star clusters in the Magellanic Clouds, including 8 clusters for which new data was obtained. We find that ${\it all}$ star clusters in our sample feature extended main sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simpl… ▽ More We present color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1 - 2 Gyr old) star clusters in the Magellanic Clouds, including 8 clusters for which new data was obtained. We find that ${\it all}$ star clusters in our sample feature extended main sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simple stellar population (including unresolved binary stars). FWHM widths of the MSTOs indicate age spreads of 200-550 Myr. We evaluate dynamical evolution of clusters with and without initial mass segregation. Our main results are: (1) the fraction of red clump (RC) stars in secondary RCs in eMSTO clusters scales with the fraction of MSTO stars having pseudo-ages $\leq 1.35$ Gyr; (2) the width of the pseudo-age distributions of eMSTO clusters is correlated with their central escape velocity $v_{\rm esc}$, both currently and at an age of 10 Myr. We find that these two results are unlikely to be reproduced by the effects of interactive binary stars or a range of stellar rotation velocities. We therefore argue that the eMSTO phenomenon is mainly caused by extended star formation within the clusters; (3) we find that $v_{\rm esc} \geq 15$ km/s out to ages of at least 100 Myr for ${\it all}$ clusters featuring eMSTOs, while $v_{\rm esc} \leq 12$ km/s at all ages for two lower-mass clusters in the same age range that do ${\it not}$ show eMSTOs. We argue that eMSTOs only occur for clusters whose early escape velocities are higher than the wind velocities of stars that provide material from which second-generation stars can form. The threshold of 12-15 km/s is consistent with wind velocities of intermediate-mass AGB stars in the literature. △ Less

Submitted 24 November, 2014; v1 submitted 14 October, 2014; originally announced October 2014.

Comments: 21 pages in emulateapj LaTeX style, incl. 12 figures, 5 tables, published in ApJ

Journal ref: The Astrophysical Journal, 797, 35 (2014)

The Quenching of the Ultra-Faint Dwarf Galaxies in the Reionization Era

Authors: Thomas M. Brown, Jason Tumlinson, Marla Geha, Joshua D. Simon, Luis C. Vargas, Don A. VandenBerg, Evan N. Kirby, Jason S. Kalirai, Roberto J. Avila, Mario Gennaro, Henry C. Ferguson, Ricardo R. Munoz, Puragra Guhathakurta, Alvio Renzini

Abstract: We present new constraints on the star formation histories of six ultra-faint dwarf galaxies: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV, and Ursa Major I. Our analysis employs a combination of high-precision photometry obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, medium-resolution spectroscopy obtained with the DEep Imaging Multi-Object Spectrogr… ▽ More We present new constraints on the star formation histories of six ultra-faint dwarf galaxies: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV, and Ursa Major I. Our analysis employs a combination of high-precision photometry obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, medium-resolution spectroscopy obtained with the DEep Imaging Multi-Object Spectrograph on the W.M. Keck Observatory, and updated Victoria-Regina isochrones tailored to the abundance patterns appropriate for these galaxies. The data for five of these Milky Way satellites are best fit by a star formation history where at least 75% of the stars formed by z~10 (13.3 Gyr ago). All of the galaxies are consistent with 80% of the stars forming by z~6 (12.8 Gyr ago) and 100% of the stars forming by z~3 (11.6 Gyr ago). The similarly ancient populations of these galaxies support the hypothesis that star formation in the smallest dark matter sub-halos was suppressed by a global outside influence, such as the reionization of the universe. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: Accepted for publication in The Astrophysical Journal. 13 pages, 5 color and 3 grayscale figures, 2 tables

Global Properties of M31's Stellar Halo from the SPLASH Survey: II. Metallicity Profile

Authors: Karoline M. Gilbert, Jason S. Kalirai, Puragra Guhathakurta, Rachael L. Beaton, Marla C. Geha, Evan N. Kirby, Steven R. Majewski, Richard J. Patterson, Erik J. Tollerud, James S. Bullock, Mikito Tanaka, Masashi Chiba

Abstract: We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars bey… ▽ More We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo. △ Less

Submitted 12 September, 2014; originally announced September 2014.

Comments: Accepted for publication in the Astrophysical Journal. 21 pages, 16 figures

New Clues to the Cause of Extended Main Sequence Turn-Offs in Intermediate-Age Star Clusters in the Magellanic Clouds

Authors: Matteo Correnti, Paul Goudfrooij, Jason S. Kalirai, Leo Girardi, Thomas H. Puzia, Leandro Kerber

Abstract: We use the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) to obtain deep, high resolution images of two intermediate-age star clusters in the Large Magellanic Cloud of relatively low mass ($\approx$ $10^4$ $M_{\odot}$) and significantly different core radii, namely NGC2209 and NGC2249. For comparison purposes, we also re-analyzed archival HST images of NGC1795 and IC2146, two… ▽ More We use the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) to obtain deep, high resolution images of two intermediate-age star clusters in the Large Magellanic Cloud of relatively low mass ($\approx$ $10^4$ $M_{\odot}$) and significantly different core radii, namely NGC2209 and NGC2249. For comparison purposes, we also re-analyzed archival HST images of NGC1795 and IC2146, two other relatively low mass star clusters. From the comparison of the observed color-magnitude diagrams with Monte Carlo simulations, we find that the main sequence turnoff (MSTO) regions in NGC2209 and NGC2249 are significantly wider than that derived from simulations of simple stellar populations, while those in NGC1795 and IC2146 are not. We determine the evolution of the clusters' masses and escape velocities from an age of 10 Myr to the present age. We find that the differences among these clusters can be explained by dynamical evolution arguments if the currently extended clusters (NGC2209 and IC2146) experienced stronger levels of initial mass segregation than the currently compact ones (NGC2249 and NGC1795). Under this assumption, we find that NGC2209 and NGC2249 have estimated escape velocities $V_{\rm esc}$ $\geq$ 15 km s$^{-1}$ at an age of 10 Myr, large enough to retain material ejected by slow winds of first-generation stars, while the two clusters that do not feature extended MSTOs have $V_{\rm esc}$ $\leq$ 12 km s$^{-1}$ at that age. These results suggest that the extended MSTO phenomenon can be better explained by a range of stellar ages rather than a range of stellar rotation velocities or interacting binaries. △ Less

Submitted 25 July, 2014; originally announced July 2014.

Comments: 9 pages (in emulateapj format), 4 figures, accepted for publication in ApJ

White Dwarf Cosmochronology in the Solar Neighborhood

Authors: P. -E. Tremblay, J. S. Kalirai, D. R. Soderblom, M. Cignoni, J. Cummings

Abstract: The study of the stellar formation history in the solar neighborhood is a powerful technique to recover information about the early stages and evolution of the Milky Way. We present a new method which consists of directly probing the formation history from the nearby stellar remnants. We rely on the volume complete sample of white dwarfs within 20 pc, where accurate cooling ages and masses have be… ▽ More The study of the stellar formation history in the solar neighborhood is a powerful technique to recover information about the early stages and evolution of the Milky Way. We present a new method which consists of directly probing the formation history from the nearby stellar remnants. We rely on the volume complete sample of white dwarfs within 20 pc, where accurate cooling ages and masses have been determined. The well characterized initial-final mass relation is employed in order to recover the initial masses (1 < M/Msun < 8) and total ages for the local degenerate sample. We correct for moderate biases that are necessary to transform our results to a global stellar formation rate, which can be compared to similar studies based on the properties of main-sequence stars in the solar neighborhood. Our method provides precise formation rates for all ages except in very recent times, and the results suggest an enhanced formation rate for the solar neighborhood in the last 5 Gyr compared to the range 5 < Age (Gyr) < 10. Furthermore, the observed total age of ~10 Gyr for the oldest white dwarfs in the local sample is consistent with the early seminal studies that have determined the age of the Galactic disk from stellar remnants. The main shortcoming of our study is the small size of the local white dwarf sample. However, the presented technique can be applied to larger samples in the future. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: 25 pages, 10 figures, accepted for publication in the Astrophysical Journal

White Dwarfs In The UKIRT Infrared Deep Sky Survey Data Release 9

Authors: P. -E. Tremblay, S. K. Leggett, N. Lodieu, B. Freytag, P. Bergeron, J. S. Kalirai, H. -G. Ludwig

Abstract: We have identified eight to ten new cool white dwarfs from the Large Area Survey (LAS) Data Release 9 of the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS). The data set was paired with the Sloan Digital Sky Survey (SDSS) to obtain proper motions and a broad ugrizYJHK wavelength coverage. Optical spectroscopic observations were secured at Gemini Observatory and confirm… ▽ More We have identified eight to ten new cool white dwarfs from the Large Area Survey (LAS) Data Release 9 of the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS). The data set was paired with the Sloan Digital Sky Survey (SDSS) to obtain proper motions and a broad ugrizYJHK wavelength coverage. Optical spectroscopic observations were secured at Gemini Observatory and confirm the degenerate status for eight of our targets. The final sample includes two additional white dwarf candidates with no spectroscopic observations. We rely on improved 1D model atmospheres and new multi-dimensional simulations with CO5BOLD to review the stellar parameters of the published LAS white dwarf sample along with our additional discoveries. Most of the new objects possess very cool atmospheres with effective temperatures below 5000 K, including two pure-hydrogen remnants with a cooling age between 8.5 and 9.0 Gyr, and tangential velocities in the range 40 km/s < vtan < 60 km/s. They are likely thick disk 10-11 Gyr-old objects. In addition we find a resolved double degenerate system with vtan ~ 155 km/s and a cooling age between 3.0 and 5.0 Gyr. These white dwarfs could be disk remnants with a very high velocity or former halo G stars. We also compare the LAS sample with earlier studies of very cool degenerates and observe a similar deficit of helium-dominated atmospheres in the range 5000 < Teff (K) < 6000. We review the possible explanations for the spectral evolution from helium-dominated towards hydrogen-rich atmospheres at low temperatures. △ Less

Submitted 1 May, 2014; originally announced May 2014.

Comments: 40 pages, 14 figures, accepted for publication in the Astrophysical Journal

The Core Mass Growth and Stellar Lifetime of Thermally Pulsing Asymptotic Giant Branch Stars

Authors: Jason S. Kalirai, Paola Marigo, Pier-Emmanuel Tremblay

Abstract: We establish new constraints on the intermediate-mass range of the initial-final mass relation by studying white dwarfs in four young star clusters, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). We show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with $M_{\rm initial}$ = 1.6 to 2.0 $M_\odot$. At la… ▽ More We establish new constraints on the intermediate-mass range of the initial-final mass relation by studying white dwarfs in four young star clusters, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). We show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with $M_{\rm initial}$ = 1.6 to 2.0 $M_\odot$. At larger masses, the core-mass growth decreases steadily to $\sim$10% at $M_{\rm initial}$ = 3.4 $M_\odot$. These observations are in excellent agreement with predictions from the latest TP-AGB evolutionary models in Marigo et al. (2013). We also compare to models with varying efficiencies of the third dredge-up and mass loss, and demonstrate that the process governing the growth of the core is largely the stellar wind, while the third dredge-up plays a secondary, but non-negligible role. Based on the new white dwarf measurements, we perform an exploratory calibration of the most popular mass-loss prescriptions in the literature. Finally, we estimate the lifetime and the integrated luminosity of stars on the TP-AGB to peak at $t$ $\sim$ 3 Myr and $E$ = 1.2 $\times$ 10$^{10}$ $L_\odot$ yr for $M_{\rm initial}$ $\sim$ 2 $M_\odot$ ($t$ $\sim$ 2 Myr for luminosities brighter than the RGB tip at $\log(L/L_{\odot})$ $>$ 3.4), decreasing to $t$ = 0.4 Myr and $E$ = 6.1 $\times$ 10$^{9}$ $L_\odot$ yr for stars with $M_{\rm initial}$ $\sim$ 3.5 $M_\odot$. The implications of these results are discussed with respect to general population synthesis studies that require correct modeling of the TP-AGB phase of stellar evolution. △ Less

Submitted 16 December, 2013; originally announced December 2013.

Comments: 14 pages, 7 figures, 4 tables. Accepted for publication in ApJ

The Formation History of the Ultra-Faint Dwarf Galaxies

Authors: Thomas M. Brown, Jason Tumlinson, Marla Geha, Evan Kirby, Don A. VandenBerg, Jason S. Kalirai, Joshua D. Simon, Roberto J. Avila, Ricardo R. Munoz, Puragra Guhathakurta, Alvio Renzini, Henry C. Ferguson, Luis C. Vargas, Mario Gennaro

Abstract: We present early results from a Hubble Space Telescope survey of the ultra-faint dwarf galaxies. These Milky Way satellites were discovered in the Sloan Digital Sky Survey, and appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in the efforts to understand the missing satellite problem. Because they are the least luminous, most dark matter domina… ▽ More We present early results from a Hubble Space Telescope survey of the ultra-faint dwarf galaxies. These Milky Way satellites were discovered in the Sloan Digital Sky Survey, and appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in the efforts to understand the missing satellite problem. Because they are the least luminous, most dark matter dominated, and least chemically evolved galaxies known, the ultra-faint dwarfs are the best candidate fossils from the early universe. The primary goal of the survey is to measure the star-formation histories of these galaxies and discern any synchronization due to the reionization of the universe. We find that the six galaxies of our survey have very similar star-formation histories, and that each is dominated by stars older than 12 Gyr. △ Less

Submitted 2 October, 2013; originally announced October 2013.

Comments: 5 pages, 2 color figures. To appear in the proceedings of "EWASS 2013 Symposium 5: Local Group, Local Cosmology," held 8-9 July 2013 in Turku, Finland. Eds. S. Salvadori and M. Monelli. To be published in Memorie della Societa' Astronomica Italiana

Comparing the White Dwarf Cooling Sequences in 47 Tuc and NGC 6397

Authors: Harvey B. Richer, Ryan Goldsbury, Jeremy Heyl, Jarrod Hurley, Aaron Dotter, Jason S. Kalirai, Kristin A. Woodley, Gregory G. Fahlman, R. Michael Rich, Michael M. Shara

Abstract: Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most o… ▽ More Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included. △ Less

Submitted 30 September, 2013; originally announced October 2013.

Comments: 7 pages, 4 figures, accepted for publication in the Astrophysical Journal

An Age Difference of 2 Gyr between a Metal-Rich and a Metal-Poor Globular Cluster

Authors: B. M. S. Hansen, J. S. Kalirai, J. Anderson, A. Dotter, H. B. Richer, R. M. Rich, M. M. Shara, G. G. Fahlman, J. R. Hurley, I. R. King, D. Reitzel, P. B. Stetson

Abstract: Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the… ▽ More Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397. △ Less

Submitted 31 July, 2013; originally announced August 2013.

Comments: Main Article: 10 pages, 4 figures; Supplementary Info 15 pages, 5 figures. Nature, Aug 1, 2013

A Dynamical Signature of Multiple Stellar Populations in 47 Tucanae

Authors: Harvey B. Richer, Jeremy Heyl, Jay Anderson, Jason S. Kalirai, Michael M. Shara, Aaron Dotter, Gregory G. Fahlman, R. Michael Rich

Abstract: Based on the width of its main sequence, and an actual observed split when viewed through particular filters, it is widely accepted that 47 Tucanae contains multiple stellar populations. In this contribution, we divide the main-sequence of 47 Tuc into four color groups, which presumably represent stars of various chemical compositions. The kinematic properties of each of these groups is explored v… ▽ More Based on the width of its main sequence, and an actual observed split when viewed through particular filters, it is widely accepted that 47 Tucanae contains multiple stellar populations. In this contribution, we divide the main-sequence of 47 Tuc into four color groups, which presumably represent stars of various chemical compositions. The kinematic properties of each of these groups is explored via proper-motions, and a strong signal emerges of differing proper-motion anisotropies with differing main-sequence color; the bluest main-sequence stars exhibit the largest proper-motion anisotropy which becomes undetectable for the reddest stars. In addition, the bluest stars are also the most centrally concentrated. A simiilar analysis for SMC stars, which are located in the background of 47 Tuc on our frames, yields none of the anisotropy exhibited by the 47 Tuc stars. We discuss implications of these results for possible formation scenarios of the various populations. △ Less

Submitted 5 June, 2013; originally announced June 2013.

Comments: Accepted for publication in the Astrophysical Journal Letters, 5 pages, 4 figures

An extended main sequence turn-off in the Small Magellanic Cloud star cluster NGC411

Authors: Leo Girardi, Paul Goudfrooij, Jason S. Kalirai, Leandro Kerber, Vera Kozhurina-Platais, Stefano Rubele, Alessandro Bressan, Rupali Chandar, Paola Marigo, Imants Platais, Thomas H. Puzia

Abstract: Based on new observations with the Wide Field Camera 3 onboard the Hubble Space Telescope, we report the discovery of an extended main sequence turn-off (eMSTO) in the intermediate-age star cluster NGC411. This is the second case of an eMSTO being identified in a star cluster belonging to the Small Magellanic Cloud (SMC), after NGC419. Despite the present masses of these two SMC clusters differ by… ▽ More Based on new observations with the Wide Field Camera 3 onboard the Hubble Space Telescope, we report the discovery of an extended main sequence turn-off (eMSTO) in the intermediate-age star cluster NGC411. This is the second case of an eMSTO being identified in a star cluster belonging to the Small Magellanic Cloud (SMC), after NGC419. Despite the present masses of these two SMC clusters differ by a factor of 4, the comparison between their colour--magnitude diagrams (CMD) shows striking similarities, especially regarding the shape of their eMSTOs. The loci of main CMD features are so similar that they can be well described, in a first approximation, by the same mean metallicity, distance and extinction. NGC411, however, presents merely a trace of secondary red clump as opposed to its prominent manifestation in NGC419. This could be due either to the small number statistics in NGC411, or by the star formation in NGC419 having continued for 60 Myr longer than in NGC411. Under the assumption that the eMSTOs are caused by different generations of stars at increasing age, both clusters are nearly coeval in their first episodes of star formation. The initial period of star formation, however, is slightly more marked in NGC419 than in NGC411. We discuss these findings in the context of possible scenarios for the origin of eMSTOs. △ Less

Submitted 6 March, 2013; originally announced March 2013.

Comments: To appear in MNRAS, 10 pages

Ultra-Deep Hubble Space Telescope Imaging of the Small Magellanic Cloud: The Initial Mass Function of Stars with M <~ 1 Msun

Authors: Jason S. Kalirai, Jay Anderson, Aaron Dotter, Harvey B. Richer, Gregory G. Fahlman, Brad M. S. Hansen, Jarrod Hurley, I. Neill Reid, R. Michael Rich, Michael M. Shara

Abstract: We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits… ▽ More We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits. The stellar main sequence of the SMC is measured in the F606W, F814W color-magnitude diagram (CMD) down to ~30th magnitude, and is cleanly separated from the foreground star cluster population using proper motions. We simulate the SMC population by extracting stellar masses (single and unresolved binaries) from specific IMFs, and converting those masses to luminosities in our bandpasses. The corresponding photometry for these simulated stars is drawn directly from a rich cloud of 4 million artificial stars, thereby accounting for the real photometric scatter and completeness of the data. Over a continuous and well populated mass range of M = 0.37 - 0.93 Msun (i.e., down to a ~75% completeness limit at F606W = 28.7), we demonstrate that the IMF is well represented by a single power-law form with slope α= -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^α). This is shallower than the Salpeter slope of α= -2.35, which agrees with the observed stellar luminosity function at higher masses. Our results indicate that the IMF does {\it not} turn over to a more shallow power-law form within this mass range. We discuss implications of this result for the theory of star formation, the inferred masses of galaxies, and the (lack of a) variation of the IMF with metallicity. △ Less

Submitted 31 December, 2012; v1 submitted 5 December, 2012; originally announced December 2012.

Comments: 9 pages, 5 figures, 1 table. Accepted for publication in ApJ. Minor revisions to second version based on ApJ proofs. Error bar on IMF slope clarified as 3-sigma interval

The Panchromatic Hubble Andromeda Treasury IV. A Probabilistic Approach to Inferring the High Mass Stellar Initial Mass Function and Other Power-law Functions

Authors: Daniel R. Weisz, Morgan Fouesneau, David W. Hogg, Hans-Walter Rix, Andrew E. Dolphin, Julianne J. Dalcanton, Daniel T. Foreman-Mackey, Dustin Lang, L. Clifton Johnson, Lori C. Beerman, Eric F. Bell, Karl D. Gordon, Dimitrios Gouliermis, Jason S. Kalirai, Evan D. Skillman, Benjamin F. Williams

Abstract: We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associa… ▽ More We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, α, are unbiased and that the uncertainty, Δα, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on α and provide an analytic approximation for Δα as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <α>=2.46 with a 1-σ dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged) △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 21 pages, 15 figures, accepted to ApJ

Global Properties of M31's Stellar Halo from the SPLASH Survey. I. Surface Brightness Profile

Authors: Karoline M. Gilbert, Puragra Guhathakurta, Rachael L. Beaton, James Bullock, Marla C. Geha, Jason S. Kalirai, Evan N. Kirby, Steven R. Majewski, James C. Ostheimer, Richard J. Patterson, Erik J. Tollerud, Mikito Tanaka, Masashi Chiba

Abstract: We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This… ▽ More We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects. △ Less

Submitted 13 November, 2012; v1 submitted 11 October, 2012; originally announced October 2012.

Comments: Accepted for publication in ApJ; v2 matches published version. Only minor editorial changes made

Journal ref: The Astrophysical Journal, Volume 760, Issue 1, article id. 76, 21pp, (2012)

The Panchromatic Hubble Andromeda Treasury III. Measuring Ages and Masses of Partially Resolved Stellar Clusters

Authors: Lori C. Beerman, L. Clifton Johnson, Morgan Fouesneau, Julianne J. Dalcanton, Daniel R. Weisz, Anil C. Seth, Ben F. Williams, Eric F. Bell, Luciana C. Bianchi, Nelson Caldwell, Andrew E. Dolphin, Dimitrios A. Gouliermis, Jason S. Kalirai, Søren S. Larsen, Jason L. Melbourne, Hans-Walter Rix, Evan D. Skillman

Abstract: The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the stellar initial mass function, when inferred from the integrated color of low mass clusters (less than ~10^4 solar masses). We use simulated star clusters to show that these effects are minimized when the brightest, rapidly evolving stars in a cluster can be resolved, and the light of the fainter,… ▽ More The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the stellar initial mass function, when inferred from the integrated color of low mass clusters (less than ~10^4 solar masses). We use simulated star clusters to show that these effects are minimized when the brightest, rapidly evolving stars in a cluster can be resolved, and the light of the fainter, more numerous unresolved stars can be analyzed separately. When comparing the light from the less luminous cluster members to models of unresolved light, more accurate age estimates can be obtained than when analyzing the integrated light from the entire cluster under the assumption that the initial mass function is fully populated. We show the success of this technique first using simulated clusters, and then with a stellar cluster in M31. This method represents one way of accounting for the discrete, stochastic sampling of the stellar initial mass function in less massive clusters and can be leveraged in studies of clusters throughout the Local Group and other nearby galaxies. △ Less

Submitted 24 September, 2012; originally announced September 2012.

Comments: 11 pages, 8 figures, accepted for publication in ApJ

An Empirical Measure of the Rate of White Dwarf Cooling in 47 Tucanae

Authors: R. Goldsbury, J. Heyl, H. B. Richer, P. Bergeron, A. Dotter, J. S. Kalirai, J. MacDonald, R. M. Rich, P. B. Stetson, P. -E. Tremblay, K. A. Woodley

Abstract: We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models ar… ▽ More We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models are then used to determine the rate at which objects are leaving the main sequence, which must be the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result is an empirically derived relation between temperature ($T_{eff}$) and time ($t$) on the white dwarf cooling sequence. Comparing this result to theoretical cooling models, we find general agreement with the expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition to the Mestel cooling rate of $T_{eff} \propto t^{-0.4}$ is found to occur at hotter temperatures, and more abruptly than is predicted by any of these models. △ Less

Submitted 21 September, 2012; originally announced September 2012.

Comments: 10 pages, 16 figures, accepted for publication in ApJ

The Panchromatic Hubble Andromeda Treasury II. Tracing the Inner M31 Halo with Blue Horizontal Branch Stars

Authors: Benjamin F. Williams, Julianne J. Dalcanton, Eric F. Bell, Karoline M. Gilbert, Puragra Guhathakurta, Tod R. Lauer, Anil C. Seth, Jason S. Kalirai, Philip Rosenfield, Leo Girardi

Abstract: We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find… ▽ More We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galactocentric distances. If we assume a power-law profile r^(-α) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of α=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to α<1.2 inside of 3 kpc. This slope is consistent with previous measurements but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass of 2.1^{+1.7}_{-0.4} x 10^9 M_sun. These properties are comparable with both simulations of stellar halo formation formed by satellite disruption alone, and with simulations that include some in situ formation of halo stars. △ Less

Submitted 11 September, 2012; originally announced September 2012.

Comments: 15 pages, 1 table, 5 figures, accepted for publication in ApJ

Stellar Kinematics of the Andromeda II Dwarf Spheroidal Galaxy

Authors: Nhung Ho, Marla C. Geha, Ricardo R. Munoz, Puragra Guhathakurta, Jason S. Kalirai, Karoline M. Gilbert, Erik J. Tollerud, James S. Bullock, Rachael L. Beaton, Steve R. Majewski

Abstract: We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity disp… ▽ More We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs. △ Less

Submitted 10 September, 2012; originally announced September 2012.

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