Source author record

Morgan Fouesneau

Morgan Fouesneau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

A stellar stream remnant of a globular cluster below the metallicity floor

Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo.

preprint2022arXiv

An empirical model of the Gaia DR3 selection function

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky. On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project: $\texttt{https://github.com/gaia-unlimited/gaiaunlimited}$

preprint2022arXiv

Estimating $\left[ α/ \text{Fe} \right]$ from Gaia low-resolution BP/RP spectra using the ExtraTrees algorithm

Gaia Data Release 3 will contain more than a billion sources with positions, parallaxes, and proper motions. In addition, for hundreds of millions of stars, it will include low-resolution blue photometer (BP) and red photometer (RP) spectra. Obtained by dispersing light with prisms, these spectra have resolutions that are too low to allow us to measure individual spectral lines and bands. However, the combined BP/RP spectra can be used to estimate some stellar properties such as $T_{\text{eff}}$, $\log{g}$, and $\left[ \text{M} / \text{H} \right]$. We investigate the feasibility of using the ExtraTrees algorithm to estimate the alpha element to iron abundance ratio $\left[ α/ \text{H} \right]$ from low-resolution BP/RP spectra. To infer $\left[ α/ \text{H} \right]$ from the spectra, we created regression models trained on two samples: a set of synthetic spectra and a set of observed spectra from stars that have known $\left[ α/ \text{H} \right]$ since they have been observed using HERMES. We applied each model to the other sample and to a larger observed sample to assess the performance of the models. In addition, we used our models to analyse stars from the Gaia-Enceladus structure. We find that models using the ExtraTrees algorithm can be used to estimate $\left[ α/ \text{H} \right]$ from low-resolution BP/RP spectra of cool stars. However, they do this by exploiting correlations between $\left[ α/ \text{H} \right]$ and other parameters, rather than the causal effect of $\left[ α/ \text{H} \right]$ on the spectrum. Hence, they are unlikely to be useful in studies that attempt to distinguish stars that only differ in $\left[ α/ \text{H} \right]$.

preprint2022arXiv

Improving white dwarfs as chronometers with Gaia parallaxes and spectroscopic metallicities

White dwarfs (WDs) offer unrealized potential in solving two problems in astrophysics: stellar age accuracy and precision. WD cooling ages can be inferred from surface temperatures and radii, which can be constrained with precision by high quality photometry and parallaxes. Accurate and precise Gaia parallaxes along with photometric surveys provide information to derive cooling and total ages for vast numbers of WDs. Here we analyse 1372 WDs found in wide binaries with MS companions, and report on the cooling and total age precision attainable in these WD+MS systems. The total age of a WD can be further constrained if its original metallicity is known, because the main-sequence (MS) progenitor lifetime depends on metallicity at fixed mass, yet metallicity is unavailable via spectroscopy of the WD. We show that incorporating spectroscopic metallicity constraints from 38 wide binary MS companions substantially decreases internal uncertainties in WD total ages compared to a uniform constraint. Averaged over the 38 stars in our sample, the total (internal) age uncertainty improves from 21.04% to 16.77% when incorporating the spectroscopic constraint. Higher mass WDs yield better total age precision; for 8 WDs with zero age main sequence masses $\geq$ 2.0 solar masses, the mean uncertainty in total ages improves from 8.61% to 4.54% when incorporating spectroscopic metallicities. We find that it is often possible to achieve 5% total age precision for WDs with progenitor masses above 2.0 solar masses if parallaxes with $\leq$ 1% precision and Pan-STARRS $g, r$, and $i$ photometry with $\leq$ 0.01 mag precision are available.

preprint2022arXiv

The Pristine survey -- XVII. The C-19 stream is dynamically hot and more extended than previously thought

The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($σ_{\rm [Fe/H]}$ $<$ 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream's properties. One of these new member stars is located 30 $°$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $σ_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$°\pm0.08°$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.

preprint2020arXiv

A Gaia early DR3 mock stellar catalog: Galactic prior and selection function

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besançon Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

preprint2020arXiv

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

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

preprint2019arXiv

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

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

preprint2016arXiv

A radial age gradient in the geometrically thick disk of the Milky Way

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

preprint2016arXiv

Panchromatic Hubble Andromeda Treasury XVI. Star Cluster Formation Efficiency and the Clustered Fraction of Young Stars

We use the Panchromatic Hubble Andromeda Treasury (PHAT) survey dataset to perform spatially resolved measurements of star cluster formation efficiency ($Γ$), the fraction of stellar mass formed in long-lived star clusters. We use robust star formation history and cluster parameter constraints, obtained through color-magnitude diagram analysis of resolved stellar populations, to study Andromeda's cluster and field populations over the last $\sim$300 Myr. We measure $Γ$ of 4-8% for young, 10-100 Myr old populations in M31. We find that cluster formation efficiency varies systematically across the M31 disk, consistent with variations in mid-plane pressure. These $Γ$ measurements expand the range of well-studied galactic environments, providing precise constraints in an HI-dominated, low intensity star formation environment. Spatially resolved results from M31 are broadly consistent with previous trends observed on galaxy-integrated scales, where $Γ$ increases with increasing star formation rate surface density ($Σ_{\mathrm{SFR}}$). However, we can explain observed scatter in the relation and attain better agreement between observations and theoretical models if we account for environmental variations in gas depletion time ($τ_{\mathrm{dep}}$) when modeling $Γ$, accounting for the qualitative shift in star formation behavior when transitioning from a H$_2$-dominated to a HI-dominated interstellar medium. We also demonstrate that $Γ$ measurements in high $Σ_{\mathrm{SFR}}$ starburst systems are well-explained by $τ_{\mathrm{dep}}$-dependent fiducial $Γ$ models.

preprint2016arXiv

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

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.

preprint2015arXiv

A Global Star Forming Episode in M31 2-4 Gyr Ago

We have identified a major global enhancement of star formation in the inner M31 disk that occurred between 2-4 Gyr ago, producing $\sim$60% of the stellar mass formed in the past 5 Gyr. The presence of this episode in the inner disk was discovered by modeling the optical resolved star color-magnitude diagrams of low extinction regions in the main disk of M31 (3$<$R$<$20 kpc) as part of the Panchromatic Hubble Andromeda Treasury. This measurement confirms and extends recent measurements of a widespread star formation enhancement of similar age in the outer disk, suggesting that this burst was both massive and global. Following the galaxy-wide burst, the star formation rate of M31 has significantly declined. We briefly discuss possible causes for these features of the M31 evolutionary history, including interactions with M32, M33 and/or a merger.

preprint2015arXiv

Panchromatic Hubble Andromeda Treasury XII. Mapping Stellar Metallicity Distributions in M31

We present a study of spatial variations in the metallicity of old red giant branch stars in the Andromeda galaxy. Photometric metallicity estimates are derived by interpolating isochrones for over seven million stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. This is the first systematic study of stellar metallicities over the inner 20 kpc of Andromeda's galactic disk. We see a clear metallicity gradient of $-0.020\pm0.004$ dex/kpc from $\sim4-20$ kpc assuming a constant RGB age. This metallicity gradient is derived after correcting for the effects of photometric bias and completeness and dust extinction and is quite insensitive to these effects. The unknown age gradient in M31's disk creates the dominant systematic uncertainty in our derived metallicity gradient. However, spectroscopic analyses of galaxies similar to M31 show that they typically have small age gradients that make this systematic error comparable to the 1$σ$ error on our metallicity gradient measurement. In addition to the metallicity gradient, we observe an asymmetric local enhancement in metallicity at radii of 3-6 kpc that appears to be associated with Andromeda's elongated bar. This same region also appears to have an enhanced stellar density and velocity dispersion.

preprint2015arXiv

Panchromatic Hubble Andromeda Treasury. XIV. The Period-Age Relationship of Cepheid Variables in M31 Star Clusters

We present a sample of 11 M31 Cepheids in stellar clusters, derived from the overlap of the Panchromatic Hubble Andromeda Treasury (PHAT) cluster catalog and the Pan-STARRS1 (PS1) disk Cepheid catalog. After identifying the PS1 Cepheids in the HST catalog, we calibrate the PS1 mean magnitudes using the higher resolution HST photometry, revealing up to 1 magnitude offsets due to crowding effects in the ground-based catalog. We measure ages of the clusters by performing single stellar population fits to their color-magnitude diagrams (CMDs) excluding their Cepheids. From these cluster age measurements, we derive an empirical period-age relation which agrees well with the existing literature values. By confirming this relation for M31 Cepheids, we justify its application in high-precision pointwise age estimation across M31.

preprint2015arXiv

PHAT Stellar Cluster Survey. II. Andromeda Project Cluster Catalog

We construct a stellar cluster catalog for the Panchromatic Hubble Andromeda Treasury (PHAT) survey using image classifications collected from the Andromeda Project citizen science website. We identify 2,753 clusters and 2,270 background galaxies within ~0.5 deg$^2$ of PHAT imaging searched, or ~400 kpc$^2$ in deprojected area at the distance of the Andromeda galaxy (M31). These identifications result from 1.82 million classifications of ~20,000 individual images (totaling ~7 gigapixels) by tens of thousands of volunteers. We show that our crowd-sourced approach, which collects >80 classifications per image, provides a robust, repeatable method of cluster identification. The high spatial resolution Hubble Space Telescope images resolve individual stars in each cluster and are instrumental in the factor of ~6 increase in the number of clusters known within the survey footprint. We measure integrated photometry in six filter passbands, ranging from the near-UV to the near-IR. PHAT clusters span a range of ~8 magnitudes in F475W (g-band) luminosity, equivalent to ~4 decades in cluster mass. We perform catalog completeness analysis using >3000 synthetic cluster simulations to determine robust detection limits and demonstrate that the catalog is 50% complete down to ~500 solar masses for ages <100 Myr. We include catalogs of clusters, background galaxies, remaining unselected candidates, and synthetic cluster simulations, making all information publicly available to the community. The catalog published here serves as the definitive base data product for PHAT cluster science, providing a census of star clusters in an L$^*$ spiral galaxy with unmatched sensitivity and quality.

preprint2015arXiv

Red giant masses and ages derived from carbon and nitrogen abundances

We show that the masses of red giant stars can be well predicted from their photospheric carbon and nitrogen abundances, in conjunction with their spectroscopic stellar labels log g, Teff, and [Fe/H]. This is qualitatively expected from mass-dependent post main sequence evolution. We here establish an empirical relation between these quantities by drawing on 1,475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from APOGEE DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional r.m.s. errors of about 14% (typically 0.2 Msun). From these masses, we derive ages with r.m.s errors of 40%. This empirical model allows us for the first time to make age determinations (in the range 1-13 Gyr) for vast numbers of giant stars across the Galaxy. We apply our model to 52,000 stars in APOGEE DR12, for which no direct mass and age information was previously available. We find that these estimates highlight the vertical age structure of the Milky Way disk, and that the relation of age with [alpha/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

preprint2015arXiv

Testing Density Wave Theory with Resolved Stellar Populations around Spiral Arms in M81

Stationary density waves rotating at a constant pattern speed $Ω_{\rm P}$ would produce age gradients across spiral arms. We test whether such age gradients are present in M81 by deriving the recent star formation histories (SFHs) of 20 regions around one of M81's grand-design spiral arms. For each region, we use resolved stellar populations to determine the SFH by modeling the observed color-magnitude diagram (CMD) constructed from archival Hubble Space Telescope (HST) F435W and F606W imaging. Although we should be able to detect systematic time delays in our spatially-resolved SFHs, we find no evidence of star formation propagation across the spiral arm. Our data therefore provide no convincing evidence for a stationary density wave with a single pattern speed in M81, and instead favor the scenario of kinematic spiral patterns that are likely driven by tidal interactions with the companion galaxies M82 and NGC 3077.

preprint2015arXiv

The High-Mass Stellar Initial Mass Function in M31 Clusters

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]

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

preprint2015arXiv

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

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]

preprint2014arXiv

Studying Gender in Conference Talks -- data from the 223rd meeting of the American Astronomical Society

We present a study on the gender balance, in speakers and attendees, at the recent major astronomical conference, the American Astronomical Society meeting 223, in Washington, DC. We conducted an informal survey, yielding over 300 responses by volunteers at the meeting. Each response included gender data about a single talk given at the meeting, recording the gender of the speaker and all question-askers. In total, 225 individual AAS talks were sampled. We analyze basic statistical properties of this sample. We find that the gender ratio of the speakers closely matched the gender ratio of the conference attendees. The audience asked an average of 2.8 questions per talk. Talks given by women had a slightly higher number of questions asked (3.2$\pm$0.2) than talks given by men (2.6$\pm$0.1). The most significant result from this study is that while the gender ratio of speakers very closely mirrors that of conference attendees, women are under-represented in the question-asker category. We interpret this to be an age-effect, as senior scientists may be more likely to ask questions, and are more commonly men. A strong dependence on the gender of session chairs is found, whereby women ask disproportionately fewer questions in sessions chaired by men. While our results point to laudable progress in gender-balanced speaker selection, we believe future surveys of this kind would help ensure that collaboration at such meetings is as inclusive as possible.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury V: Ages and Masses of the Year 1 Stellar Clusters

We present ages and masses for 601 star clusters in M31 from the analysis of the six filter integrated light measurements from near ultraviolet to near infrared wavelengths, made as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We derive the ages and masses using a probabilistic technique, which accounts for the effects of stochastic sampling of the stellar initial mass function. Tests on synthetic data show that this method, in conjunction with the exquisite sensitivity of the PHAT observations and their broad wavelength baseline, provides robust age and mass recovery for clusters ranging from $\sim 10^2 - 2 \times 10^6 M_\odot$. We find that the cluster age distribution is consistent with being uniform over the past $100$ Myr, which suggests a weak effect of cluster disruption within M31. The age distribution of older ($>100$ Myr) clusters fall towards old ages, consistent with a power-law decline of index $-1$, likely from a combination of fading and disruption of the clusters. We find that the mass distribution of the whole sample can be well-described by a single power-law with a spectral index of $-1.9 \pm 0.1$ over the range of $10^3-3 \times 10^5 M_\odot$. However, if we subdivide the sample by galactocentric radius, we find that the age distributions remain unchanged. However, the mass spectral index varies significantly, showing best fit values between $-2.2$ and $-1.8$, with the shallower slope in the highest star formation intensity regions. We explore the robustness of our study to potential systematics and conclude that the cluster mass function may vary with respect to environment.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury. VII. The Steep Mid-Ultraviolet to Near-Infrared Extinction Curve in the Central 200 pc of the M31 Bulge

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5μm), using Swift/UVOT and HST WFC3/ACS observations in thirteen bands. Taking advantage of the high angular resolution of the HST WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (RV =2.4-2.5) are steeper than the average Galactic one (RV =3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic Bulge (RV~2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size <2 pc) exhibits a strong UV bump (extinction at 2175A), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central super-massive black hole, resulting in the observed steepened extinction curve.

preprint2014arXiv

The Supernova Progenitor Mass Distributions of M31 and M33: Further Evidence for an Upper Mass Limit

Using Hubble Space Telescope (HST) photometry to measure star formation histories, we age-date the stellar populations surrounding supernova remnants (SNRs) in M31 and M33. We then apply stellar evolution models to the ages to infer the corresponding masses for their supernova progenitor stars. We analyze 33 M33 SNR progenitors and 29 M31 SNR progenitors in this work. We then combine these measurements with 53 previously published M31 SNR progenitor measurements to bring our total number of progenitor mass estimates to 115. To quantify the mass distributions, we fit power laws of the form $dN/dM \propto M^{-α}$. Our new, larger sample of M31 progenitors follows a distribution with $α= 4.4\pm 0.4$, and the M33 sample follows a distribution with $α= 3.8^{+0.4}_{-0.5}$. Thus both samples are consistent within the uncertainties, and the full sample across both galaxies gives $α= 4.2\pm 0.3$. Both the individual and full distributions display a paucity of massive stars when compared to a Salpeter initial mass function (IMF), which we would expect to observe if all massive stars exploded as SN that leave behind observable SNR. If we instead fix $α= 2.35$ and treat the maximum mass as a free parameter, we find $M_{max} \sim 35-45M_{sun}$, indicative of a potential maximum cutoff mass for SN production. Our results suggest that either SNR surveys are biased against finding objects in the youngest (<10 Myr old) regions, or the highest mass stars do not produce SNe.

preprint2012arXiv

Analysing star cluster populations with stochastic models: the HST/WFC3 sample of clusters in M83

The majority of clusters in the Universe have masses well below 10^5 Msun. Hence their integrated fluxes and colors can be affected by the random presence of a few bright stars introduced by stochastic sampling of the stellar mass function. Specific methods are being developed to extend the analysis of cluster SEDs into the low-mass regime. In this paper, we apply such a method to observations of star clusters, in the nearby spiral galaxy M83. We reassess ages and masses of a sample of 1242 objects for which UBVIHalpha fluxes were obtained with the HST/WFC3 images. Synthetic clusters with known properties are used to characterize the limitations of the method. The ensemble of color predictions of the discrete cluster models are in good agreement with the distribution of observed colors. We emphasize the important role of the Halpha data in the assessment of the fraction of young objects, particularly in breaking the age-extinction degeneracy that hampers an analysis based on UBVI only. We find the mass distribution of the cluster sample to follow a power-law of index -2.1 +/-0.2, and the distribution of ages a power-law of index -1.0 +/-0.2 for M > 10^3.5 Msun and ages between 10^7 and 10^9 yr. An extension of our main method, that makes full use of the probability distributions of age and mass of the individual clusters, is explored. It produces similar power-law slopes and will deserve further investigation. Although the properties derived for individual clusters significantly differ from those obtained with traditional, non-stochastic models in ~30% of the objects, the first order aspect of the age and mass distributions are similar to those obtained previously for this M83 sample in the range of overlap of the studies. We extend the power-law description to lower masses with better mass and age resolution and without most of the artifacts produced by the classical method.

preprint2012arXiv

PHAT Stellar Cluster Survey I. Year 1 Catalog and Integrated Photometry

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an on-going Hubble Space Telescope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, representing more than a factor of four increase over previous catalogs within the current survey area (390 arcmin^2). This work presents results derived from the first \sim25% of the survey data; we estimate that the final sample will include \sim2500 clusters. For the Year 1 objects, we present a catalog with positions, radii, and six-band integrated photometry. Along with a general characterization of the cluster luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and highlight a number of individually interesting clusters found in the Year 1 search.

preprint2012arXiv

Supernova Remnant Progenitor Masses in M31

Using HST photometry, we age-date 59 supernova remnants (SNRs) in the spiral galaxy M31 and use these ages to estimate zero-age main sequence masses (MZAMS) for their progenitors. To accomplish this, we create color-magnitude diagrams (CMDs) and use CMD fitting to measure the recent star formation history (SFH) of the regions surrounding cataloged SNR sites. We identify any young coeval population that likely produced the progenitor star and assign an age and uncertainty to that population. Application of stellar evolution models allows us to infer the MZAMS from this age. Because our technique is not contingent on precise location of the progenitor star, it can be applied to the location of any known SNR. We identify significant young SF around 53 of the 59 SNRs and assign progenitor masses to these, representing a factor of 2 increase over currently measured progenitor masses. We consider the remaining 6 SNRs as either probable Type Ia candidates or the result of core-collapse progenitors that have escaped their birth sites. The distribution of recovered progenitor masses is bottom heavy, showing a paucity of the most massive stars. If we assume a single power law distribution, dN/dM proportional to M^alpha, we find a distribution that is steeper than a Salpeter IMF (alpha=-2.35). In particular, we find values of alpha outside the range -2.7 to -4.4 inconsistent with our measured distribution at 95% confidence. If instead we assume a distribution that follows a Salpeter IMF up to some maximum mass, we find that values of M_max greater than 26 Msun are inconsistent with the measured distribution at 95% confidence. In either scenario, the data suggest that some fraction of massive stars may not explode. The result is preliminary and requires more SNRs and further analysis. In addition, we use our distribution to estimate a minimum mass for core collapse between 7.0 and 7.8 Msun.

preprint2012arXiv

The ACS Nearby Galaxy Survey Treasury. X. Quantifying the Star Cluster Formation Efficiency of Nearby Dwarf Galaxies

We study the relationship between the field star formation and cluster formation properties in a large sample of nearby dwarf galaxies. We use optical data from the Hubble Space Telescope and from ground-based telescopes to derive the ages and masses of the young (t_age < 100Myr) cluster sample. Our data provides the first constraints on two proposed relationships between the star formation rate of galaxies and the properties of their cluster systems in the low star formation rate regime. The data show broad agreement with these relationships, but significant galaxy-to-galaxy scatter exists. In part, this scatter can be accounted for by simulating the small number of clusters detected from stochastically sampling the cluster mass function. However, this stochasticity does not fully account for the observed scatter in our data suggesting there may be true variations in the fraction of stars formed in clusters in dwarf galaxies. Comparison of the cluster formation and the brightest cluster in our sample galaxies also provide constraints on cluster destruction models.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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)

preprint2010arXiv

Accounting for Stochastic Fluctuations when Analysing Integrated Light of Star Clusters. I: First Systematics

Star clusters are studied widely both as benchmarks for stellar evolution models and in their own right. Cluster age and mass distributions within galaxies are probes of star formation histories, and of cluster formation and disruption processes. The vast majority of clusters in the Universe is small, and it is well known that the integrated fluxes and colors have broad probability distributions, due to small numbers of bright stars. This paper goes beyond the description of predicted probability distributions, and presents results of the analysis of cluster energy distributions in an explicitly stochastic context. The method developed is Bayesian. It provides posterior probability distributions in the age-mass-extinction space, using multi-wavelength photometric observations and a large collection of Monte-Carlo simulations of clusters of finite stellar masses. Both UBVI and UBVIK datasets are considered, and the study conducted in this paper is restricted to the solar metallicity. We first reassess and explain errors arising from the use of standard analysis methods, which are based on continuous population synthesis models: systematic errors on ages and random errors on masses are large, while systematic errors on masses tend to be smaller. The age-mass distributions obtained after analysis of a synthetic sample are very similar to those found for real galaxies in the literature. The Bayesian approach on the other hand, is very successful in recovering the input ages and masses. Taking stochastic effects into account is important, more important for instance than the choice of adding or removing near-IR data in many cases. We found no immediately obvious reason to reject priors inspired by previous (standard) analyses of cluster populations in galaxies, i.e. cluster distributions that scale with mass as M^-2 and are uniform on a logarithmic age scale.