sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Richer, Harvey B.

Normalized to: Richer, H.

87 article(s) in total. 179 co-authors, from 1 to 43 common article(s). Median position in authors list is 2,0.

[1]  [pdf] - 2103456
Carbon stars as standard candles: I. The luminosity function of carbon stars in the Magellanic Clouds
Comments: 10 pages, 12 figures. Accepted for publication in MNRAS. Results unchanged. Correction of some uncertainty values in Eq. (24) and Eq. (25). Added some acknowledgements. Minor changes to the text. Codes and catalogues available on GitLab, see
Submitted: 2020-05-12, last modified: 2020-05-27
Our goal in this paper is to derive a carbon-star luminosity function that will eventually be used to determine distances to galaxies at $50$-$60$ Mpc and hence yield a value of the Hubble constant. Cool N-type carbon stars exhibit redder near-infrared colours than oxygen-rich stars. Using Two Micron All Sky Survey near-infrared photometry and the Gaia Data Release 2, we identify carbon stars in the Magellanic Clouds (MC) and the Milky Way (MW). Carbon stars in the MC appear as a distinct horizontal feature in the near-infrared ($(J-K_s)_0$, $M_J$) colour-magnitude diagram. We build a colour selection ($1.4 < (J-K_s)_0 < 2$) and derive the luminosity function of the colour-selected carbon stars. We find the median absolute magnitude and the dispersion, in the J band, for the Large Magellanic Cloud and the Small Magellanic Cloud to be, respectively, ($\bar{M_J} = -6.284 \pm 0.004$, $\sigma = 0.352 \pm 0.005$) and ($\bar{M_J} = -6.160 \pm 0.015$, $\sigma = 0.365 \pm 0.014$). The difference between the MC may be explained by the lower metallicity of the Small Magellanic Cloud, but in any case it provides limits on the type of galaxy whose distance can be determined with this technique. To account for metallicity effects, we developed a composite magnitude, named C, for which the error-weighted mean C magnitude of both the MC are equal. Thanks to the next generation of telescopes (JWST, ELT, TMT), carbon stars could be detected in MC-type galaxies at distances out to $50$-$60$ Mpc. The final goal is to eventually try and improve the measurement of the Hubble constant while exploring the current tensions related to its value.
[2]  [pdf] - 2097209
Small and Moderate Aperture Telescopes for Research and Education
Comments: White paper submitted to the Canadian Long Range Plan LRP2020
Submitted: 2019-10-01
In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade.
[3]  [pdf] - 1925099
A Massive Magnetic Helium Atmosphere White Dwarf Binary in a Young Star Cluster
Comments: 10 pages, 7 figures, accepted for publication in Astrophysical Journal
Submitted: 2019-06-11
We have searched the Gaia DR2 catalogue for previously unknown hot white dwarfs in the direction of young open star clusters. The aim of this experiment was to try and extend the initial-final mass relation (IFMR) to somewhat higher masses, potentially providing a tension with the Chandrasekhar limit currently thought to be around 1.38 M$_{\odot}$. We discovered a particularly interesting white dwarf in the direction of the young $\sim$150 Myr old cluster Messier 47 (NGC 2422). All Gaia indicators (proper motion, parallax, location in the Gaia colour-magnitude diagram) suggest that it is a cluster member. Its spectrum, obtained from Gemini South, yields a number of anomalies: it is a DB (helium-rich atmosphere) white dwarf, it has a large magnetic field (2.5 MG), is of high mass ($\sim$1.06 M$_\odot$) and its colours are very peculiar --- particularly the redder ones ($r$, $i$, $z$ and $y$), which suggest that it has a late-type companion. This is the only magnetized, detached binary white dwarf with a non-degenerate companion of any spectral type known in or out of a star cluster. If the white dwarf is a cluster member, as all indicators suggest, its progenitor had a mass just over 6 M$_\odot$. It may, however, be telling an even more interesting story than the one related to the IFMR, one about the origin of stellar magnetic fields, Type I supernovae and gravitational waves from low mass stellar systems.
[4]  [pdf] - 1877314
A Multi-Mass Velocity Dispersion Model of 47 Tucanae Indicates No Evidence for an Intermediate Mass Black Hole
Comments: Updated version since referee report includes: -- More thorough discussion of procedure throughout body of text --New and updated figures for clarity of process and results --New tables included with model parameters --Extended treatment of stellar remnant retention fraction and its impact on results. (13 pages, 9 figures)
Submitted: 2018-07-09, last modified: 2019-05-03
In this paper, we analyze stellar proper motions in the core of the globular cluster 47 Tucanae to explore the possibility of an intermediate-mass black hole (IMBH) influence on the stellar dynamics. Our use of short-wavelength photometry affords us an exceedingly clear view of stellar motions into the very center of the crowded core, yielding proper motions for $>$50,000 stars in the central 2'. We model the velocity dispersion profile of the cluster using an isotropic Jeans model. The density distribution is taken as a central IMBH point mass added to a combination of King templates. We individually model the general low-mass cluster objects (main sequence/giant stars), as well as the concentrated populations of heavy binary systems and dark stellar remnants. Using unbinned likelihood model fitting, we find that the inclusion of the concentrated populations in our model plays a crucial role in fitting for an IMBH mass. The concentrated binaries and stellar-mass black holes (BHs) produce a sufficient velocity dispersion signal in the core so as to make an IMBH unnecessary to fit the observations. We additionally determine that a stellar-mass BH retention fraction of $\gtrsim 8.5\%$ becomes incompatible with our observed velocities in the core.
[5]  [pdf] - 1868205
Hunting for Brown Dwarfs in the Globular Cluster M4: 2nd epoch HST NIR observations
Comments: 13 pages, including 7 figures. Accepted for publication in MNRAS
Submitted: 2019-04-05
We present an analysis of the second epoch HST WFC3 F110W near-Infrared (NIR) imaging data of the globular cluster M4. The new dataset suggests that one of the previously suggested four brown dwarf candidates in this cluster is indeed a high-probability cluster member. The position of this object in the NIR colour magnitude diagrams (CMDs) is in the white dwarf/brown dwarf area. The source is too faint to be a low-mass main sequence star, but, according to theoretical considerations, also most likely somewhat too bright to be a bona-fide brown dwarf. Since we know that the source is a cluster member, we determined a new optical magnitude estimate at the position the source should have in the optical image. This new estimate places the source closer to the white dwarf sequence in the optical-NIR CMD and suggests that it might be a very cool (T_eff < 4500 K) white dwarf at the bottom of the white dwarf cooling sequence in M4, or a white dwarf/brown dwarf binary. We cannot entirely exclude the possibility that the source is a very massive, bright brown dwarf, or a very low-mass main sequence star, however, we conclude that we still have not convincingly detected a brown dwarf in a globular cluster, but we expect to be very close to the start of the brown dwarf cooling sequence in this cluster. We also note that the main sequence ends at F110W approx.22.5 mag in the proper-motion cleaned CMDs, where completeness is still high.
[6]  [pdf] - 1850851
Hunting for ancient brown dwarfs: the developing field of brown dwarfs in globular clusters
Comments: White paper submitted for Astro2020 Decadal Survey. 8 pages, 2 figures
Submitted: 2019-03-15
The detection of brown dwarfs in globular star clusters will allow us to break the degeneracies in age, mass and composition that affect our current models, and therefore to constrain the physics of their atmospheres and interiors. Furthermore, detecting brown dwarfs will help us constrain the properties of the clusters themselves, as they carry information about the clusters' age and dynamics. Great advancements in this field are to be expected in the next ten years, thanks to the extraordinary sensitivity in the infrared of upcoming telescopes like JWST and the ELTs.
[7]  [pdf] - 1825567
The HST Large Programme on NGC 6752. I. Serendipitous discovery of a dwarf Galaxy in background
Comments: 5 pages, 1 table, 3 figures (2 in low resolution). Published on MNRAS Letters and freely available Related press releases: and , (and only in italian: )
Submitted: 2019-02-01
As part of a large Hubble Space Telescope investigation aiming at reaching the faintest stars in the Galactic globular cluster NGC 6752, an ACS/WFC field was the subject of deep optical observations reaching magnitudes as faint as V~30. In this field we report the discovery of Bedin I a dwarf spheroidal galaxy too faint and too close to the core of NGC 6752 for detection in earlier surveys. As it is of broad interest to complete the census of galaxies in the local Universe, in this Letter we provide the position of this new object along with preliminary assessments of its main param eters. Assuming the same reddening as for NGC 6752, we estimate a distance modulus of (m-M)_0 = 29.70+/-0.13 from the observed red giant branch, i.e., 8.7 (+0.5 -0.7) Mpc, and size of ~840x340pc, about 1/5 the size of the LMC. A comparison of the observed colour-magnitude diagram with synthetic counterparts that account for the galaxy distance modulus, reddening and photometric errors, suggests the presence of an old (~13Gyr) and metal poor ([Fe/H]~-1.3) population. This object is most likely a relatively isolated satellite dwarf spheroidal galaxy of the nearby great spiral NGC 6744, or potentially the most distant isolated dwarf spheroidal known with a secure distance.
[8]  [pdf] - 1823111
The HST Large Programme on NGC6752. II. Multiple populations at the bottom of the main sequence probed in NIR
Comments: 7 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2019-01-22
Historically, multiple populations in Globular Clusters (GCs) have been mostly studied from ultraviolet and optical filters down to stars that are more massive than ~0.6 solar masses. Here we exploit deep near-infrared (NIR) photometry from the Hubble Space Telescope to investigate multiple populations among M-dwarfs in the GC NGC6752. We discovered that the three main populations (A, B and C), previously observed in the brightest part of the color-magnitude diagram, define three distinct sequences that run from the main-sequence (MS) knee towards the bottom of the MS (~0.15 solar masses). These results, together with similar findings on NGC2808, M4, and omega Centauri, demonstrate that multiple sequences of M-dwarfs are common features of the color-magnitude diagrams of GCs. The three sequences of low-mass stars in NGC6752 are consistent with stellar populations with different oxygen abundances. The range of [O/Fe] needed to reproduce the NIR CMD of NGC6752 is similar to the oxygen spread inferred from high-resolution spectroscopy of red-giant branch (RGB) stars. The relative numbers of stars in the three populations of M-dwarfs are similar to those derived among RGB and MS stars more massive than ~0.6 solar masses. As a consequence, the evidence that the properties of multiple populations do not depend on stellar mass is a constraint for the formation scenarios.
[9]  [pdf] - 1783826
Distances to the globular clusters 47 Tucanae and NGC 362 using Gaia DR2 parallaxes
Comments: 15 pages, 9 figures, accepted for publication in ApJ
Submitted: 2018-07-18, last modified: 2018-09-10
Using parallaxes from Gaia DR2, we estimate the distance to the globular clusters 47 Tuc and NGC 362, taking advantage of the background stars in the Small Magellanic Cloud and quasars to account for various parallax systematics. We found the parallax to be dependent on the Gaia DR2 G apparent magnitude for stars with $13<\textrm{G}<18$ , where brighter stars have a lower parallax zero point than fainter stars. The distance to 47 Tuc was found to be $4.45\pm0.01\pm0.12$ kpc, and for NGC 362 $8.54\pm0.20\pm0.44$ kpc with random and systematic errors listed respectively. This is the first time a precise distance measurement directly using parallaxes has been done for either of these two globular clusters.
[10]  [pdf] - 1602746
Deep HST Imaging in 47 Tucanae: A Global Dynamical Model
Comments: 16 pages, 13 figures, accepted by Astrophysical Journal
Submitted: 2017-10-29
Multi-epoch observations with ACS and WFC3 on HST provide a unique and comprehensive probe of stellar dynamics within 47 Tucanae. We confront analytic models of the globular cluster with the observed stellar proper motions that probe along the main sequence from just above 0.8 to 0.1M$_\odot$ as well as white dwarfs younger than one gigayear. One field lies just beyond the half-light radius where dynamical models (\eg lowered Maxwellian distributions) make robust predictions for the stellar proper motions. The observed proper motions in this outer field show evidence for anisotropy in the velocity distribution as well as skewness; the latter is evidence of rotation. The measured velocity dispersions and surface brightness distributions agree in detail with a rotating, anisotropic model of the stellar distribution function with mild dependence of the proper-motion dispersion on mass. However, the best fitting models under-predict the rotation and skewness of the stellar velocities. In the second field, centered on the core of the cluster, the mass segregation in proper motion is much stronger. Nevertheless the model developed in the outer field can be extended inward by taking this mass segregation into account in a heuristic fashion. The proper motions of the main-sequence stars yield a mass estimate of the cluster of $1.31 \pm 0.02 \times 10^6 \mathrm{M}_\odot$ at a distance of 4.7 kpc. By comparing the proper motions of a sample of giant and sub-giant stars with the observed radial velocities we estimate the distance to the cluster kinematically to be $4.29 \pm 0.47$ kpc.
[11]  [pdf] - 1588782
The Onset of Convective Coupling and Freezing in the White Dwarfs of 47 Tucanae
Comments: 6 pages, 5 figures. Submitted to MNRAS
Submitted: 2017-09-23, last modified: 2017-10-13
Using images from the Hubble Space Telescope Advanced Camera for Surveys, we measure the rate of cooling of white dwarfs in the globular cluster 47 Tucanae and compare it to modelled cooling curves. We examine the effects of the outer convective envelope reaching the nearly isothermal degenerate core and the release of latent heat during core crystallisation on the white dwarf cooling rates. For white dwarfs typical of 47 Tuc, the onset of these effects occur at similar times. The latent heat released during crystallisation is a small heat source. In contrast, the heat reservoir of the degenerate core is substantially larger. When the convective envelope reaches the nearly isothermal interior of the white dwarf, the star becomes brighter than it would be in the absence of this effect. Our modelled cooling curves that include this convective coupling closely match the observed luminosity function of the white dwarfs in 47 Tuc.
[12]  [pdf] - 1587324
Chemical Mapping of the Milky Way With The Canada-France Imaging Survey: A Non-parametric Metallicity-Distance Decomposition of the Galaxy
Comments: 21 pages, 17 figures, accepted for publication in ApJ
Submitted: 2017-08-21
We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.
[13]  [pdf] - 1587321
The Canada-France Imaging Survey: First results from the u-band component
Comments: 12 pages, 11 figures, accepted for publication in ApJ
Submitted: 2017-08-21
The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.
[14]  [pdf] - 1535263
Globular cluster absolute ages from cooling brown dwarfs
Submitted: 2017-01-31
Globular clusters are the oldest conglomerates of stars in our Galaxy and can be useful laboratories to test theories from stellar evolution to cosmology. In this paper, we present a new method to estimate the absolute age of a globular cluster from observations of its brown dwarfs. The transition region between the end of the main sequence and the brown dwarf regime is characterized by a dearth of objects as function of magnitude. The brightest of the cooling brown dwarfs is easily identified by an increase in density in the color-magnitude diagram as you go fainter in magnitudes, and these brightest brown dwarfs get fainter with age. By identifying the brightest brown dwarfs, it is thus possible to determine the age of a globular cluster within a 1 Gyr precision with four-sigma confidence. This new method, which is independent of current methods of age estimation and which does not rely on the knowledge of the cluster's distance from Earth, will become feasible thanks to the high spatial resolution and incredible infrared sensitivity of the James Webb Space Telescope.
[15]  [pdf] - 1502111
Formation and Evolution of Blue Stragglers in 47 Tucanae
Comments: 11 pages, 4 figures, 6 tables, accepted in ApJ
Submitted: 2016-09-07
Blue stragglers (BSS) are stars whose position in the Color-Magnitude Diagram (CMD) places them above the main sequence turn-off (TO) point of a star cluster. Using data from the core of 47 Tuc in the ultraviolet (UV), we have identified various stellar populations in the CMD, and used their radial distributions to study the evolution and origin of BSS, and obtain a dynamical estimate of the mass of BSS systems. When we separate the BSS into two samples by their magnitude, we find that the bright BSS show a much more centrally concentrated radial distribution and thus higher mass estimate (over twice the TO mass for these BSS systems), suggesting an origin involving triple or multiple stellar systems. In contrast, the faint BSS are less concentrated, with a radial distribution similar to the main sequence (MS) binaries, pointing to the MS binaries as the likely progenitors of these BSS. Putting our data together with available photometric data in the visible and using MESA evolutionary models, we calculate the expected number of stars in each evolutionary stage for the normal evolution of stars and the number of stars coming from the evolution of BSS. The results indicate that BSS have a post-MS evolution comparable to that of a normal star of the same mass and a MS BSS lifetime of about 200-300 Myr. We also find that the excess population of asymptotic giant branch (AGB) stars in 47 Tuc is due to evolved BSS.
[16]  [pdf] - 1447765
Dynamical estimate of post main sequence stellar masses in 47 Tucanae
Comments: 8 pages, 7 figures, accepted for publication in The Astrophysical Journal
Submitted: 2016-05-18
We use the effects of mass segregation on the radial distribution of different stellar populations in the core of 47 Tucanae to find estimates for the masses of stars at different post main sequence evolutionary stages. We take samples of main sequence (MS) stars from the core of 47 Tucanae, at different magnitudes (i.e. different masses), and use the effects of this dynamical process to develop a relation between the radial distance (RD) at which the cumulative distribution reaches the 20th and 50th percentile, and stellar mass. From these relations we estimate the masses of different post MS populations. We find that mass remains constant for stars going through the evolutionary stages between the upper MS up to the horizontal branch (HB). By comparing RDs of the HB stars with stars of lower masses, we can exclude a mass loss greater than 0.09M during the red giant branch (RGB) stage at nearly the 3{\sigma} level. The slightly higher mass estimates for the asymptotic giant branch (AGB) are consistent with the AGB having evolved from somewhat more massive stars. The AGB also exhibits evidence of contamination by more massive stars, possibly blue stragglers (BSS), going through the RGB phase. We do not include the BSS in this paper due to the complexity of these objects, instead, the complete analysis of this population is left for a companion paper. The process to estimate the masses described in this paper are exclusive to the core of 47 Tuc.
[17]  [pdf] - 1392860
Constraining white dwarf structure and neutrino physics in 47 Tucanae
Comments: 12 pages, 15 figures
Submitted: 2016-02-19
We present a robust statistical analysis of the white dwarf cooling sequence in 47 Tucanae. We combine HST UV and optical data in the core of the cluster, Modules for Experiments in Stellar Evolution (MESA) white dwarf cooling models, white dwarf atmosphere models, artificial star tests, and a Markov Chain Monte Carlo (MCMC) sampling method to fit white dwarf cooling models to our data directly. We use a technique known as the unbinned maximum likelihood to fit these models to our data without binning. We use these data to constrain neutrino production and the thickness of the hydrogen layer in these white dwarfs. The data prefer thicker hydrogen layers $(q_\mathrm{H}=3.2\e{-5})$ and we can strongly rule out thin layers $(q_\mathrm{H}=10^{-6})$. The neutrino rates currently in the models are consistent with the data. This analysis does not provide a constraint on the number of neutrino species.
[18]  [pdf] - 1347463
Deep near-IR observations of the Globular Cluster M4: Hunting for Brown Dwarfs
Comments: ApJ accepted, 28 pages including 16 figures
Submitted: 2015-06-03, last modified: 2015-11-16
We present an analysis of deep HST/WFC3 near-IR (NIR) imaging data of the globular cluster M4. The best-photometry NIR colour-magnitude diagram (CMD) clearly shows the main sequence extending towards the expected end of the Hydrogen-burning limit and going beyond this point towards fainter sources. The white dwarf sequence can be identified. As such, this is the deepest NIR CMD of a globular cluster to date. Archival HST optical data were used for proper-motion cleaning of the CMD and for distinguishing the white dwarfs (WDs) from brown dwarf (BD) candidates. Detection limits in the NIR are around F110W approx 26.5 mag and F160W approx27 mag, and in the optical around F775W approx 28 mag. Comparing our observed CMDs with theoretical models, we conclude that we have reached beyond the H-burning limit in our NIR CMD and are probably just above or around this limit in our optical-NIR CMDs. Thus, any faint NIR sources that have no optical counterpart are potential BD candidates, since the optical data are not deep enough to detect them. We visually inspected the positions of NIR sources which are fainter than the H-burning limit in F110W and for which the optical photometry did not return a counterpart. We found in total five sources for which we did not get an optical measurement. For four of these five sources, a faint optical counterpart could be visually identified, and an upper optical magnitude was estimated. Based on these upper optical magnitude limits, we conclude that one source is likely a WD, one source could either be a WD or BD candidate, and the remaining two sources agree with being BD candidates. For only one source no optical counterpart could be detected, which makes this source a good BD candidate. We conclude that we found in total four good BD candidates.
[19]  [pdf] - 1296119
When do stars in 47 Tucanae lose their mass?
Comments: 8 pages, 4 figures, version accepted by the Astrophysical Journal
Submitted: 2015-02-25, last modified: 2015-08-05
By examining the diffusion of young white dwarfs through the core of the globular cluster 47 Tucanae, we estimate the time when the progenitor star lost the bulk of its mass to become a white dwarf. According to stellar evolution models of the white-dwarf progenitors in 47 Tucanae, we find this epoch to coincide approximately with the star ascending the asymptotic giant branch ($3.0 \pm 8.1$ Myr before the tip of the AGB) and more than ninety million years after the helium flash (with ninety-percent confidence). From the diffusion of the young white dwarfs we can exclude the hypothesis that the bulk of the mass loss occurs on the red-giant branch at the four-sigma level. Furthermore, we find that the radial distribution of horizontal branch stars is consistent with that of the red-giant stars and upper-main-sequence stars and inconsistent with the loss of more than 0.2 solar masses on the red-giant branch at the six-sigma level.
[20]  [pdf] - 1265787
Constraining Neutrino Cooling using the Hot White Dwarf Luminosity Function in the Globular Cluster 47 Tucanae
Comments: 13 pages, 12 figures, to appear in the Astrophysical Journal
Submitted: 2015-07-20
We present Hubble Space Telescope observations of the upper part (T_eff> 10 000 K) of the white dwarf cooling sequence in the globular cluster 47 Tucanae and measure a luminosity function of hot white dwarfs. Comparison with previous determinations from large scale field surveys indicates that the previously determined plateau at high effective temperatures is likely a selection effect, as no such feature is seen in this sample. Comparison with theoretical models suggests that the current estimates of white dwarf neutrino emission (primarily by the plasmon channel) are accurate, and variations are restricted to no more than a factor of two globally, at 95% confidence. We use these constraints to place limits on various proposed exotic emission mechanisms, including a non-zero neutrino magnetic moment, formation of axions, and emission of Kaluza-Klein modes into extra dimensions.
[21]  [pdf] - 1224346
A Measurement of Diffusion in 47 Tucanae
Comments: 12 pages, 10 figures, accept for publication in the Astrophysical Journal (added a short passage on blue stragglers descendants and corrected a minor typographical error in Eq. 15)
Submitted: 2015-02-06, last modified: 2015-02-20
Using images from the Hubble Space Telescope Wide-Field Camera 3, we measure the rate of diffusion of stars through the core of the globular cluster 47 Tucanae using a sample of young white dwarfs identified in these observations. This is the first direct measurement of diffusion due to gravitational relaxation. We find that the diffusion rate $\kappa\approx 10-13$ arcsecond$^2$ Myr$^{-1}$ is consistent with theoretical estimates of the relaxation time in the core of 47 Tucanae of about 70 Myr.
[22]  [pdf] - 1179607
Comparing the White Dwarf Cooling Sequences in 47 Tuc and NGC 6397
Comments: 7 pages, 4 figures, accepted for publication in the Astrophysical Journal
Submitted: 2013-09-30
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.
[23]  [pdf] - 1173448
Quantifying mass segregation and new core radii for 54 milky way globular clusters
Comments: 7 pages, 5 figures, Accepted for publication in ApJ - Aug. 15, 2013
Submitted: 2013-08-16
We present core radii for 54 Milky Way globular clusters determined by fitting King-Michie models to cumulative projected star count distributions. We find that fitting star counts rather than surface brightness profiles produces results that differ significantly due to the presence of mass segregation. The sample in each cluster is further broken down into various mass groups, each of which is fit independently, allowing us to determine how the concentration of each cluster varies with mass. The majority of the clusters in our sample show general agreement with the standard picture that more massive stars will be more centrally concentrated. We find that core radius vs. stellar mass can be fit with a two parameter power-law. The slope of this power-law is a value that describes the amount of mass segregation present in the cluster, and is measured independently of our distance from the cluster. This value correlates strongly with the core relaxation time and physical size of each cluster. Supplementary figures are also included showing the best fits and likelihood contours of fit parameters for all 54 clusters.
[24]  [pdf] - 1173123
An Age Difference of 2 Gyr between a Metal-Rich and a Metal-Poor Globular Cluster
Comments: Main Article: 10 pages, 4 figures; Supplementary Info 15 pages, 5 figures. Nature, Aug 1, 2013
Submitted: 2013-07-31
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.
[25]  [pdf] - 1171837
A Dynamical Signature of Multiple Stellar Populations in 47 Tucanae
Comments: Accepted for publication in the Astrophysical Journal Letters, 5 pages, 4 figures
Submitted: 2013-06-05
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.
[26]  [pdf] - 1158298
Ultra-Deep Hubble Space Telescope Imaging of the Small Magellanic Cloud: The Initial Mass Function of Stars with M <~ 1 Msun
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
Submitted: 2012-12-05, last modified: 2012-12-31
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 \alpha = -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^{\alpha}). This is shallower than the Salpeter slope of \alpha = -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.
[27]  [pdf] - 1151803
Deep HST Imaging in NGC 6397: Stellar Dynamics
Comments: 25 pages, 20 figures, accepted for publication in Astrophysical Journal
Submitted: 2012-10-02, last modified: 2012-10-04
Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr.
[28]  [pdf] - 1151567
An Empirical Measure of the Rate of White Dwarf Cooling in 47 Tucanae
Comments: 10 pages, 16 figures, accepted for publication in ApJ
Submitted: 2012-09-21
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.
[29]  [pdf] - 508801
Proper Motions and Internal Dynamics in the Core of the Globular Cluster M71
Comments: 5 pages, 4 figures, accepted for publication in Astrophysical Journal Letters
Submitted: 2012-03-23
We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of anisotropy in the motions and no radial variation in the proper motions with respect to distance from the cluster center. We also set an upper limit on any central black hole to be ~150 Msun at 90% confidence level.
[30]  [pdf] - 1092201
A Deep, Wide-Field, and Panchromatic View of 47 Tuc and the SMC with HST: Observations and Data Analysis Methods
Comments: Accepted for Publication in the Astronomical Journal. 23 pages, 13 figures, and 4 tables. Companion paper by K. Woodley is astro-ph/0374778
Submitted: 2011-12-06
In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.
[31]  [pdf] - 1092200
The Spectral Energy Distributions of White Dwarfs in 47 Tucanae: The Distance to the Cluster
Comments: Accepted to The Astronomical Journal, 13 Figures, 2 Tables. Figures 3 and 6 are figure sets, each composed of 59 subfigures (to appear in the electronic journal). This is a Companion paper to the article ID: submit/0375614
Submitted: 2011-12-06
We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample.
[32]  [pdf] - 1052163
Density Variations in the NW Star Stream of M31
Comments: ApJ revised version submitted
Submitted: 2011-02-17
The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the density of the stream varies between zero and about three times the mean along its length on scales of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is less than 10^-5. As a control sample we search for density variations at precisely the same location in stars with metallicity higher than the stream, [Fe/H]=[0, -0.5] and find no variations above the expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, acting on an approximately 10Gyr old star stream. The effects of a single set of halo substructure realizations are shown for illustration, reserving a statistical comparison for another study.
[33]  [pdf] - 1034313
The ACS Survey of Galactic Globular Clusters. X. New Determinations of Centers for 65 Clusters
Comments: 17 pages, 8 figures, Accepted for publication in AJ, supplementary material will be available upon publication
Submitted: 2010-08-16, last modified: 2010-08-27
We present new measurements of the centers for 65 Milky Way globular clusters. Centers were determined by fitting ellipses to the density distribution within the inner $2\arcmin$ of the cluster center, and averaging the centers of these ellipses. The symmetry of clusters was also analyzed by comparing cumulative radial distributions on opposite sides of the cluster across a grid of trial centers. All of the determinations were done with stellar positions derived from a combination of two single-orbit ACS images of the core of the cluster in $F606W$ and $F814W$. We find that the ellipse-fitting method provides remarkable accuracy over a wide range of core sizes and density distributions, while the symmetry method is difficult to use on clusters with very large cores, or low density. The symmetry method requires a larger field, or a very sharply peaked density distribution.
[34]  [pdf] - 1018218
Star clusters as laboratories for stellar and dynamical evolution
Comments: 27 pages, 8 figures. To appear as invited review article in a special issue of the Phil. Trans. Royal Soc. A: Ch. 4 "Star clusters as tracers of galactic star-formation histories" (ed. R. de Grijs). Fully peer reviewed. LaTeX, requires rspublic.cls style file
Submitted: 2009-11-04
Open and globular star clusters have served as benchmarks for the study of stellar evolution due to their supposed nature as simple stellar populations of the same age and metallicity. After a brief review of some of the pioneering work that established the importance of imaging stars in these systems, we focus on several recent studies that have challenged our fundamental picture of star clusters. These new studies indicate that star clusters can very well harbour multiple stellar populations, possibly formed through self-enrichment processes from the first-generation stars that evolved through post-main-sequence evolutionary phases. Correctly interpreting stellar evolution in such systems is tied to our understanding of both chemical-enrichment mechanisms, including stellar mass loss along the giant branches, and the dynamical state of the cluster. We illustrate recent imaging, spectroscopic and theoretical studies that have begun to shed new light on the evolutionary processes that occur within star clusters.
[35]  [pdf] - 901695
The Masses of Population II White Dwarfs
Comments: Accepted for Publication in Astrophys. J. on Aug. 05th, 2009. 19 pages including 9 figures and 2 tables (journal format)
Submitted: 2009-09-11
Globular star clusters are among the first stellar populations to have formed in the Milky Way, and thus only a small sliver of their initial spectrum of stellar types are still burning hydrogen on the main-sequence today. Almost all of the stars born with more mass than 0.8 M_sun have evolved to form the white dwarf cooling sequence of these systems, and the distribution and properties of these remnants uniquely holds clues related to the nature of the now evolved progenitor stars. With ultra-deep HST imaging observations, rich white dwarf populations of four nearby Milky Way globular clusters have recently been uncovered, and are found to extend an impressive 5 - 8 magnitudes in the faint-blue region of the H-R diagram. In this paper, we characterize the properties of these population II remnants by presenting the first direct mass measurements of individual white dwarfs near the tip of the cooling sequence in the nearest of the Milky Way globulars, M4. Based on Gemini/GMOS and Keck/LRIS multiobject spectroscopic observations, our results indicate that 0.8 M_sun population II main-sequence stars evolving today form 0.53 +/- 0.01 M_sun white dwarfs. We discuss the implications of this result as it relates to our understanding of stellar structure and evolution of population II stars and for the age of the Galactic halo, as measured with white dwarf cooling theory.
[36]  [pdf] - 901696
The Spectral Types of White Dwarfs in Messier 4
Comments: Accepted for Publication in Astrophys. J. 11 pages including 4 figures and 2 tables (journal format)
Submitted: 2009-09-11
We present the spectra of 24 white dwarfs in the direction of the globular cluster Messier 4 obtained with the Keck/LRIS and Gemini/GMOS spectrographs. Determining the spectral types of the stars in this sample, we find 24 type DA and 0 type DB (i.e., atmospheres dominated by hydrogen and helium respectively). Assuming the ratio of DA/DB observed in the field with effective temperature between 15,000 - 25,000 K, i.e., 4.2:1, holds for the cluster environment, the chance of finding no DBs in our sample due simply to statistical fluctuations is only 6 X 10^(-3). The spectral types of the ~100 white dwarfs previously identified in open clusters indicate that DB formation is strongly suppressed in that environment. Furthermore, all the ~10 white dwarfs previously identified in other globular clusters are exclusively type DA. In the context of these two facts, this finding suggests that DB formation is suppressed in the cluster environment in general. Though no satisfactory explanation for this phenomenon exists, we discuss several possibilities.
[37]  [pdf] - 901583
The remnants of galaxy formation from a panoramic survey of the region around M31
Comments: Published in Nature. Supplementary movie available at
Submitted: 2009-09-02
In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.
[38]  [pdf] - 21119
The Dynamical Effects of White Dwarf Birth Kicks in Globular Star Clusters
Comments: 5 pages, accepted to ApJL; minor changes to reflect accepted version
Submitted: 2009-02-06, last modified: 2009-03-02
Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 suggest that WDs receive a kick of a few km/s shortly before they are born. Using our Monte Carlo cluster evolution code, which includes accurate treatments of all relevant physical processes operating in globular clusters, we study the effects of the kicks on their host cluster and on the WD population itself. We find that in clusters whose velocity dispersion is comparable to the kick speed, WD kicks are a significant energy source for the cluster, prolonging the initial cluster core contraction phase significantly so that at late times the cluster core to half-mass radius ratio is a factor of up to ~ 10 larger than in the no-kick case. WD kicks thus represent a possible resolution of the large discrepancy between observed and theoretically predicted values of this key structural parameter. Our modeling also reproduces the observed trend for younger WDs to be more extended in their radial distribution in the cluster than older WDs.
[39]  [pdf] - 14048
A Glimpse into the Past: The Recent Evolution of Globular Clusters
Comments: 5 pages, 2 figures. Accepted for publication in ApJ Letters. Higher resolution Figure 1 available at
Submitted: 2008-06-30
We present the serendipitous discovery of 195 extragalactic globular clusters (GCs) in one of the deepest optical images ever obtained, a 126 orbit HST/ACS imaging study of the nearby Galactic GC NGC 6397. The distant GCs are all found surrounding a bright elliptical galaxy in the field, and are among the faintest objects detected in the image, with magnitudes 26 < F814W < 30. We measure the redshift of the parent elliptical galaxy, using GMOS on Gemini South, to be z = 0.089 (375 Mpc). This galaxy, and its associated clusters, therefore ranks as one of the most distant such systems discovered to date. The measured light from these clusters was emitted 1.2 Gyr ago (the lookback time) and therefore the optical properties hold clues for understanding the evolution of GCs over the past Gyr. We measure the color function of the bright GCs and find that both a blue and red population exist, and that the colors of each sub-population are redder than GCs in local elliptical galaxies of comparable luminosity. For the blue clusters, the observed color difference from z = 0.089 to today is only slightly larger than predictions from stellar evolution (e.g., changes in the luminosity and color of the main-sequence turnoff and the morphology of the horizontal branch). A larger color difference is found in the red clusters, possibly suggesting that they are very metal-rich and/or significantly younger than 12 Gyr.
[40]  [pdf] - 13856
A trio of new Local Group galaxies with extreme properties
Comments: 32 pages, 7 figures (ApJ preprint format). Accepted for publication in ApJ
Submitted: 2008-06-25
We report on the discovery of three new dwarf galaxies in the Local Group. These galaxies are found in new CFHT/MegaPrime g,i imaging of the south-western quadrant of M31, extending our extant survey area to include the majority of the southern hemisphere of M31's halo out to 150 kpc. All these galaxies have stellar populations which appear typical of dwarf spheroidal (dSph) systems. The first of these galaxies, Andromeda XVIII, is the most distant Local Group dwarf discovered in recent years, at ~1.4 Mpc from the Milky Way (~ 600 kpc from M31). The second galaxy, Andromeda XIX, a satellite of M31, is the most extended dwarf galaxy known in the Local Group, with a half-light radius of r_h ~ 1.7 kpc. This is approximately an order of magnitude larger than the typical half-light radius of many Milky Way dSphs, and reinforces the difference in scale sizes seen between the Milky Way and M31 dSphs (such that the M31 dwarfs are generally more extended than their Milky Way counterparts). The third galaxy, Andromeda XX, is one of the faintest galaxies so far discovered in the vicinity of M31, with an absolute magnitude of order M_V ~ -6.3. Andromeda XVIII, XIX and XX highlight different aspects of, and raise important questions regarding, the formation and evolution of galaxies at the extreme faint-end of the luminosity function. These findings indicate that we have not yet sampled the full parameter space occupied by dwarf galaxies, although this is an essential pre-requisite for successfully and consistently linking these systems to the predicted cosmological dark matter sub-structure.
[41]  [pdf] - 1000661
A Deep HST Study of the Globular Cluster NGC 6397: Reduction Methods
Comments: 35 pages, preprint style
Submitted: 2008-03-05
We describe here the reduction methods that we developed to study the faintest red dwarfs and white dwarfs in an outer field of NGC6397, which was observed by \hst for 126 orbits in 2005. The particular challenge of this data set is that the faintest stars are not readily visible in individual exposures, so special care must be taken to combine the information in all the exposures in order to identify and measure them. Unfortunately, it is hard to find the faintest stars without also finding a large number of faint galaxies, so we developed specialized tools to distinguish between the point-like stars and the barely resolved galaxies. We found that artificial-star tests, while obviously necessary for completeness determination, can also play an important role in helping us optimize our finding and measuring algorithms. Although this paper focuses on this data set specifically, many of the techniques are new and might find application in other work, particularly when a large number of images is available for a single field.
[42]  [pdf] - 1000654
The binary fraction in the globular cluster NGC 6397
Comments: 19 pages, 7 figures, accepted in AJ
Submitted: 2008-02-29, last modified: 2008-03-04
Using Hubble Space Telescope (HST) observations of the globular cluster NGC6397, we constrain the cluster's binary fraction. The observations consist of two fields: the primary science field, a single ACS pointing centered approximately 5' from the cluster center; and the parallel field, a single WFPC2 field centered on the cluster center. Using the exquisite photometric precision of these observations, we determine the binary fraction in these regions of the cluster by examining stars lying off the main sequence. The binary fraction is constrained to be 0.012 +/- 0.004 in the ACS field, and to be 0.051 +/- 0.010 in the WFPC field. N-body simulations by Hurley et al. (2007) suggest that the binary fraction remains nearly constant beyond the half-mass radius for the lifetime of the cluster. In the context of these simulations, our results suggest that NGC 6397 had a primordial binary fraction of only ~1%.
[43]  [pdf] - 10503
Deep ACS Imaging in the Globular Cluster NGC6397: Dynamical Models
Comments: 34 pages, 11 figures, submitted to AJ, companion paper to 0708.4030
Submitted: 2008-02-28
We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.
[44]  [pdf] - 2586
The Initial-Final Mass Relation: Direct Constraints at the Low Mass End
Comments: 18 pages, 12 figures, 3 tables. Submitted to Astrophys. J. Revised from first version
Submitted: 2007-06-27, last modified: 2007-10-16
The initial-final mass relation represents a mapping between the mass of a white dwarf remnant and the mass that the hydrogen burning main-sequence star that created it once had. The relation thus far has been constrained using a sample of ~40 stars in young open clusters, ranging in initial mass from ~2.75 -- 7 Msun, and shows a general trend that connects higher mass main-sequence stars with higher mass white dwarfs. In this paper, we present CFHT/CFH12K photometric and Keck/LRIS multiobject spectroscopic observations of a sample of 22 white dwarfs in two older open clusters, NGC 7789 (t = 1.4 Gyr) and NGC 6819 (t = 2.5 Gyr). We measure masses for the highest S/N spectra by fitting the Balmer lines to atmosphere models and place the first direct constraints on the low mass end of the initial-final mass relation. Our results indicate that the observed general trend at higher masses continues down to low masses, with M_initial = 1.16 Msun main-sequence stars forming M_final = 0.53 Msun white dwarfs (including our data from the very old open cluster, NGC 6791). This extention of the relation represents a four fold increase in the total number of hydrogen burning stars for which the integrated mass loss can now be calculated, assuming a Salpeter initial mass function. The new leverage at the low mass end is used to derive a purely empirical initial-final mass relation without the need for any indirectly derived anchor points. The sample of white dwarfs in these clusters also shows several very interesting systems that we discuss further: a DB (helium atmosphere) white dwarf, a magnetic white dwarf, a DAB (mixed hydrogen/helium atmosphere or a double degenerate DA+DB) white dwarf(s), and two possible equal mass DA double degenerate binary systems.
[45]  [pdf] - 122232
Stellar Evolution in NGC 6791: Mass Loss on the Red Giant Branch and the Formation of Low Mass White Dwarfs
Comments: 15 pages, 9 figures, 2 tables. Accepted for publication in Astrophys. J. Very minor changes from first version
Submitted: 2007-05-07, last modified: 2007-10-04
We present the first detailed study of the properties (temperatures, gravities, and masses) of the NGC 6791 white dwarf population. This unique stellar system is both one of the oldest (8 Gyr) and most metal-rich ([Fe/H] ~ 0.4) open clusters in our Galaxy, and has a color-magnitude diagram (CMD) that exhibits both a red giant clump and a much hotter extreme horizontal branch. Fitting the Balmer lines of the white dwarfs in the cluster, using Keck/LRIS spectra, suggests that most of these stars are undermassive, <M> = 0.43 +/- 0.06 Msun, and therefore could not have formed from canonical stellar evolution involving the helium flash at the tip of the red giant branch. We show that at least 40% of NGC 6791's evolved stars must have lost enough mass on the red giant branch to avoid the flash, and therefore did not convert helium into carbon-oxygen in their core. Such increased mass loss in the evolution of the progenitors of these stars is consistent with the presence of the extreme horizontal branch in the CMD. This unique stellar evolutionary channel also naturally explains the recent finding of a very young age (2.4 Gyr) for NGC 6791 from white dwarf cooling theory; helium core white dwarfs in this cluster will cool ~3 times slower than carbon-oxygen core stars and therefore the corrected white dwarf cooling age is in fact ~7 Gyr, consistent with the well measured main-sequence turnoff age. These results provide direct empirical evidence that mass loss is much more efficient in high metallicity environments and therefore may be critical in interpreting the ultraviolet upturn in elliptical galaxies.
[46]  [pdf] - 5398
On the Radial Distribution of White Dwarfs in the Globular Cluster NGC 6397
Comments: 13 pages including 1 table and 3 figures. Accepted for publication in the MNRAS Letters
Submitted: 2007-09-26
We have examined the radial distribution of white dwarfs over a single HST/ACS field in the nearby globular cluster NGC 6397. In relaxed populations, such as in a globular cluster, stellar velocity dispersion, and hence radial distribution, is directly dependent on stellar masses. The progenitors of very young cluster white dwarfs had a mass of ~0.8 solar masses, while the white dwarfs themselves have a mass of ~0.5 solar masses. We thus expect young white dwarfs to have a concentrated radial distribution (like that of their progenitors) that becomes more extended over several relaxation times to mimic that of ~0.5 solar mass main-sequence stars. However, we observe young white dwarfs to have a significantly extended radial distribution compared to both the most massive main sequence stars in the cluster and also to old white dwarfs.
[47]  [pdf] - 4376
Deep ACS Imaging in the Globular Cluster NGC 6397: The Cluster Color Magnitude Diagram and Luminosity Function
Comments: 43 pages including 4 tables and 12 diagrams. Figures 2 and 3 have been bitmapped. Accepted for publication in the Astronomical Journal
Submitted: 2007-08-29
We present the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.
[48]  [pdf] - 88867
The White Dwarf Cooling Sequence of NGC6397
Comments: 56 pages, 30 figures
Submitted: 2007-01-25, last modified: 2007-02-20
We present the results of a deep Hubble Space Telescope (HST) exposure of the nearby globular cluster NGC6397, focussing attention on the cluster's white dwarf cooling sequence. This sequence is shown to extend over 5 magnitudes in depth, with an apparent cutoff at magnitude F814W=27.6. We demonstrate, using both artificial star tests and the detectability of background galaxies at fainter magnitudes, that the cutoff is real and represents the truncation of the white dwarf luminosity function in this cluster. We perform a detailed comparison between cooling models and the observed distribution of white dwarfs in colour and magnitude, taking into account uncertainties in distance, extinction, white dwarf mass, progenitor lifetimes, binarity and cooling model uncertainties. After marginalising over these variables, we obtain values for the cluster distance modulus and age of \mu_0 = 12.02 \pm 0.06 and T_c = 11.47 \pm 0.47Gyr (95% confidence limits). Our inferred distance and white dwarf initial-final mass relations are in good agreement with other independent determinations, and the cluster age is consistent with, but more precise than, prior determinations made using the main sequence turnoff method. In particular, within the context of the currently accepted \Lambda CDM cosmological model, this age places the formation of NGC6397 at a redshift z=3, at a time when the cosmological star formation rate was approaching its peak.
[49]  [pdf] - 89263
Probing the Faintest Stars in a Globular Star Cluster
Comments: 12 pages, 4 figures. Full Resolution Figures in Science, 2006, 313, 936
Submitted: 2007-02-07
NGC 6397 is the second closest globular star cluster to the Sun. Using 5 days of time on the Hubble Space Telescope, we have constructed the deepest ever color-magnitude diagram for this cluster. We see a clear truncation in each of its two major stellar sequences. Faint red main sequence stars run out well above our observational limit and near to the theoretical prediction for the lowest mass stars capable of stable hydrogen-burning in their cores. We also see a truncation in the number counts of faint blue stars, namely white dwarfs. This reflects the limit to which the bulk of the white dwarfs can cool over the lifetime of the cluster. There is also a turn towards bluer colors in the least luminous of these objects. This was predicted for the very coolest white dwarfs with hydrogen-rich atmospheres as the formation of H2 causes their atmospheres to become largely opaque to infrared radiation due to collision-induced absorption.
[50]  [pdf] - 88910
The Space Motion of the Globular Cluster NGC 6397
Comments: 5 pages including 3 figures, accepted for publication in the Astrophysical Journal Letters. Very minor changes in V2. typos fixed
Submitted: 2007-01-26, last modified: 2007-02-06
As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.
[51]  [pdf] - 67252
The Initial-Final Mass Relationship: Spectroscopy of White Dwarfs in NGC 2099 (M37)
Comments: 5 pages including 3 figures and 1 table, accepted in ApJ Letters. Minor additions to content and typos corrected
Submitted: 2004-09-07, last modified: 2004-12-06
We present new observations of very faint white dwarfs (WDs) in the rich open star cluster NGC 2099 (M37). Following deep, wide field imaging of the cluster using CFHT, we have now obtained spectroscopic observations of candidate WDs using both GMOS on Gemini and LRIS on Keck. Of our 24 WD candidates (all fainter than V = 22.4), 21 are spectroscopically confirmed to be bona fide WDs, 4-5 of which are most likely field objects. Fitting 18 of the 21 WD spectra with model atmospheres, we find that most WDs in this cluster are quite massive (0.7-0.9 Msun), as expected given the cluster's young age (650 Myr) and, hence, high turnoff mass (~2.4 Msun). We determine a new initial-final mass relationship and almost double the number of existing data points from previous studies. The results indicate that stars with initial masses between 2.8 and 3.4 Msun lose 70-75% of their mass through stellar evolution. For the first time, we find some evidence of a metallicity dependence on the initial-final mass relationship.
[52]  [pdf] - 67253
The Dearth of Massive, Helium-Rich White Dwarfs in Young Open Star Clusters
Comments: 4 pages including 1 figure and 1 table, accepted in ApJ Letters. References added and discussion expanded in Section 5
Submitted: 2004-09-07, last modified: 2004-12-06
Spectra have been obtained of 21 white dwarfs (WDs) in the direction of the young, rich open star cluster NGC 2099. This represents an appreciable fraction (>30%) of the cluster's total WD population. The mean derived mass of the sample is 0.8 Msun - about 0.2 Msun larger than the mean seen among field WDs. A surprising result is that all of the NGC 2099 WDs have hydrogen-rich atmospheres (DAs) and none exhibit helium-rich ones (DBs), or any other spectral class. The number ratio in the field at the temperatures of the NGC 2099 WDs is DA/DB ~ 3.5. While the probability of seeing no DB WDs in NGC 2099 solely by chance is ~2%, if we include WDs in other open clusters of similar age it then becomes highly unlikely that the dearth of DB WDs in young open clusters is just a statistical fluctuation. We explore possible reasons for the lack of DBs in these clusters and conclude that the most promising scenario for the DA/DB number ratio discrepancy in young clusters is that hot, high-mass WDs do not develop large enough helium convection zones to allow helium to be brought to the surface and turn a hydrogen-rich WD into a helium-rich one.
[53]  [pdf] - 68694
Carbon Stars and other Luminous Stellar Populations in M33
Comments: 20 pages, 20 figures. Accepted for publication in The Astronomical Journal
Submitted: 2004-11-03
The M33 galaxy is a nearby, relatively metal-poor, late-type spiral. Its proximity and almost face-on inclination means that it projects over a large area on the sky, making it an ideal candidate for wide-field CCD mosaic imaging. Photometry was obtained for more than 10^6 stars covering a 74' x 56' field centered on M33. Main sequence (MS), supergiant branch (SGB), red giant branch (RGB) and asymptotic giant branch (AGB) populations are identified and classified based on broad-band V and I photometry. Narrow-band filters are used to measure spectral features allowing the AGB population to be further divided into C and M-star types. The galactic structure of M33 is examined using star counts, colour-colour and colour-magnitude selected stellar populations. We use the C to M-star ratio to investigate the metallicity gradient in the disk of M33. The C/M-star ratio is found to increase and then flatten with increasing galactocentric radius in agreement with viscous disk formation models. The C-star luminosity function is found to be similar to M31 and the SMC, suggesting that C-stars should be useful distance indicators. The ``spectacular arcs of carbon stars'' in M33 postulated recently by Block et al. (2004) are found in our work to be simply an extension of M33's disk.
[54]  [pdf] - 62310
HST Observations of the Main Sequence of M4
Comments: 37 pages, 14 figures, submitted to AJ
Submitted: 2004-01-21
This paper presents new results from a photometric study of the main-sequence stars in M4, the closest globular cluster to the sun. Multi-field, multi-epoch observations at approximately 1, 2, and 6 core-radii were obtained with the WFPC2 camera on the HST. The multi-epoch observations allowed us to clean the data on the basis of proper motion. In all the fields the cluster main sequence can be traced to at least V = 27.0 but there remains a trail of stars to the limit of the data near V = 30 in the deepest outer field. There is no evidence that we have reached the end of the hydrogen burning main sequence in any of our fields, however, there is some indication that very few stars remain to be detected in the deepest data. The scatter about the cluster main sequence yields a surprisingly small and variable binary fraction; f_b = 2% in the inner parts of the cluster falling to the 1% range outside. For the currently visible main sequence stars, the cluster mass function is very flat (alpha =0.1) in the outer field and flattens further in the inner fields suggesting well developed mass segregation. The observed variation in the mass function is broadly consistent with isotropic, multi-mass Michie-King models. Because we have a large sample of white dwarfs in the outer field we are able to show that the cluster IMF above 0.8Msun was considerably steeper than the present day mass function for low mass stars. Two appendicies are included in this contribution; a detailed discussion of the techniques used to reduce the data and a direct comparison between the cluster stars and those belonging to the inner spheroid of the Galaxy yielding a Galactocentric distance of 7.5+/-0.6 kpc.
[55]  [pdf] - 62313
HST Observations of the White Dwarf Cooling Sequence of M4
Comments: many pages, 43 postscript figures, submitted to ApJ
Submitted: 2004-01-21
We investigate in detail the white dwarf cooling sequence of the globular cluster Messier 4. In particular we study the influence of various systematic uncertainties, both observational and theoretical, on the determination of the cluster age from the white dwarf cooling sequence. These include uncertainties in the distance to the cluster and the extinction along the line of sight, as well as the white dwarf mass, envelope and core compositions and the white dwarf --main sequence mass relation. We find that fitting to the full two-dimensional colour-magnitude diagram offers a more robust method for age determination than the traditional method of fitting the one-dimensional white dwarf luminosity function. After taking into account the various uncertainties, we find a best fit age of 12.1 Gyr, with a 95% lower limit of 10.3 Gyr. We also perform fits using two other sets of cooling models from the literature. The models of Chabrier et al (2000) yield an encouragingly similar result, although the models of Salaris et al (2000) do not provide as good a fit. Our results support our previous determination of a delay between the formation of the Galactic halo and the onset of star formation in the Galactic disk.
[56]  [pdf] - 62316
Concerning the White Dwarf Cooling Age of M4: A Response to the Paper by De Marchi et al. on "A Different Interpretation of Recent HST Observations"
Comments: 14 pages, 5 figures, AJ in press
Submitted: 2004-01-21
We respond to the paper by De Marchi, Paresce, Straniero and Moroni (2003) on the white dwarf cooling age of M4. The authors question the data analysis and interpretation that led to the conclusions in Hansen et al. (2002). In their paper, De Marchi et al. are unable to obtain photometry as deep as ours from the same data set and therefore assert that only a lower limit to the white dwarf cooling age for this cluster of approximately 9 Gyr can be obtained. In this short contribution we show that shortcomings in the data analysis and reduction techniques of De Marchi et al. are responsible for their inability to reach the photometry limits that our study reports. In a forthcoming paper in which the complete techniques for age determination with white dwarfs are laid out, we demonstrate that their method of fitting the luminosity function gives a spuriously low white dwarf cooling age.
[57]  [pdf] - 55944
The Galactic Inner Halo: Searching for White Dwarfs and Measuring the Fundamental Galactic Constant, Vo/Ro
Comments: 31 pages, including 6 diagrams and 2 tables. Accepted for publication in ApJ
Submitted: 2003-04-01, last modified: 2003-10-03
We establish an extragalactic, zero-motion frame of reference within the deepest optical image of a globular star cluster, an HST 123-orbit exposure of M4 (GO 8679, cycle 9). The line of sight beyond M4 (l,b (deg) = 351,16) intersects the inner halo (spheroid) of our Galaxy at a tangent-point distance of 7.6 kpc (for Ro = 8 kpc). We isolate these spheroid stars from the cluster based on their proper motions over the 6-year baseline between these and previous epoch HST data (GO 5461, cycle 4). Distant background galaxies are also found on the same sight line using image-morphology techniques. This fixed reference frame allows us to independently determine the fundamental Galactic constant, Vo/Ro = 25.3 +/- 2.6 km/s/kpc, thus providing a velocity of the Local Standard of Rest, v = 202.7 +/- 24.7 km/s for Ro = 8.0 +/- 0.5 kpc. Secondly, the galaxies allow a direct measurement of M4's absolute proper motion, mu_total = 22.57 +/- 0.76 mas/yr, in excellent agreement with recent studies. The clear separation of galaxies from stars in these deep data also allow us to search for inner-halo white dwarfs. We model the conventional Galactic contributions of white dwarfs along our line of sight and predict 7.9 (thin disk), 6.3 (thick disk) and 2.2 (spheroid) objects to the limiting magnitude at which we can clearly delineate stars from galaxies (V = 29). An additional 2.5 objects are expected from a 20% white dwarf dark halo consisting of 0.5 Mo objects, 70% of which are of the DA type. After considering the kinematics and morphology of the objects in our data set, we find the number of white dwarfs to be consistent with the predictions for each of the conventional populations. However, we do not find any evidence for dark halo white dwarfs.
[58]  [pdf] - 57322
The CFHT Open Star Cluster Survey. IV. Two Rich, Young Open Star Clusters: NGC 2168 (M35) and NGC 2323 (M50)
Comments: 34 pages, including 12 diagrams and 5 tables. Accepted for publication in AJ. Minor typos corrected
Submitted: 2003-06-11, last modified: 2003-10-03
We continue our study of rich Galactic clusters by presenting deep CCD observations of both NGC 2168 (M35) and NGC 2323 (M50). Both clusters are found to be rich (NGC 2168 contains at least 1000 stars brighter than V = 22 and NGC 2323 contains approximately 2100 stars brighter than our photometric limit of V = 23) and young (age of NGC 2168 = 180 Myrs, age of NGC 2323 = 130 Myrs). The color-magnitude diagrams for the clusters exhibit clear main sequences stretching over 14 magnitudes in the V, B-V plane. Comparing these long main sequences with those of earlier clusters in the survey, as well as with the Hyades, has allowed for accurate distances to be established for each cluster (dist. of NGC 2168 = 912 +/- 70/65 pc, dist. of NGC 2323 = 1000 +/- 81/75 pc). Analysis of the luminosity and mass functions suggest that despite their young ages, both clusters are somewhat dynamically relaxed exhibiting signs of mass-segregation. This is especially interesting in the case of NGC 2323, which has an age of only 1.3 times the dynamical relaxation time. The present photometry is also deep enough to detect all of the white dwarfs in both clusters. We discuss some interesting candidates which may be the remnants of quite massive (M > 5 Mo) progenitor stars. The white dwarf cooling age of NGC 2168 is found to be in good agreement with the main-sequence turn-off age. These objects are potentially very important for setting constraints on the white dwarf initial-final mass relationship and upper mass limit for white dwarf production.
[59]  [pdf] - 56378
Searching for Variability in the Globular Cluster Messier 4
Comments: 36 pages including 12 figures, accepted to the Astronomical Journal. References added; typos fixed; paper made more concise at referee's recommendation
Submitted: 2003-04-25, last modified: 2003-10-03
Time-series data taken with the Hubble Space Telescope of a field six core radii (~5') from the center of the globular cluster Messier 4, and covering a period of about 10 weeks in early 2001, have been analyzed in search of variable objects. Various criteria were employed to select candidate variable stars. Period searches were performed on the selected candidates using phase dispersion minimization (PDM). The reliability of the PDM search results was tested using synthetic light curves of eclipsing binary stars and sinusoidal light curves of different periods. Results from this analysis showed that there are probably no eclipsing binary stars or periodic variables in our field with periods on the order of a few hours to a few days, down to limiting magnitudes of V~25 and I~24, which is consistent with the absence of contact binaries such as W Ursae Majoris systems. However, one candidate variable star does show an increase in brightness of ~0.1 magnitudes in both bandpasses, which seems to last for a few days. Possible explanations concerning the nature of this object include a binary system with a white dwarf primary and a low-mass main sequence secondary, or a BY Draconis variable star. We are able to set an upper limit to the observed fraction of photometric variability in this dataset of 0.05%.
[60]  [pdf] - 59329
The Pulsar-White Dwarf-Planet System in Messier 4: Improved Astrometry
Comments: 9 pages, 2 figures, accepted for publication in the Astrophysical Journal (Letters)
Submitted: 2003-09-17
A young and undermassive white dwarf has been identified as the possible companion to the millisecond pulsar PSR B1620-26 in Messier 4. This association is important as it then helps constrain the mass of the third body in the system to be of order a few times that of Jupiter. The presence of this planet in M4 has critical implications for planetary formation mechanisms in metal-poor environments such as globular clusters and the early Universe. The identification of the white dwarf is purely via the agreement in position between it and the pulsar and was limited by the accuracy of the pointing of HST which is +/-0.7 arcsec. We have redetermined the position of the white dwarf using ground-based data tied to USNOB-1.0 and find that the pulsar and white dwarf are now coincident to within 0.12 +/- 0.13 arcsec further strengthening the case for association between the two. We have also attempted to improve the proper motion measurement of the white dwarf by a maximum likelihood analysis of the stellar positions measured over a baseline of 5 years. While the errors are reduced by almost a factor of 6 from our previous work, we still have not resolved the cluster's intrinsic dispersion in proper motion. Thus the proper motion of the white dwarf with respect to the cluster itself is still not known although it is very small and is within 2-sigma of that of the cluster internal dispersion.
[61]  [pdf] - 58052
A Young White Dwarf Companion to Pulsar 1620-26: Evidence for Early Planet Formation
Comments: 4 pages, published in Science
Submitted: 2003-07-17
The pulsar B1620-26 has two companions, one of stellar mass and one of planetary mass. We detected the stellar companion using Hubble Space Telescope observations. The color and magnitude of the stellar companion indicate that it is an undermassive white dwarf ($0.34 \pm 0.04 M_{\odot}$) of age $480\pm 140$ Myr. This places a constraint on the recent history of this triple system and supports a scenario in which the current configuration arose through a dynamical exchange interaction in the cluster core. This implies that planets may be relatively common in low-metallicity globular clusters and that planet formation is more widespread and happened earlier than previously believed.
[62]  [pdf] - 49160
The Lower Main Sequence and Mass Function of the Globular Cluster Messier 4
Comments: 12 pages, 3 diagrams. To be published in the Astrophysical Journal Letters August 2002. Several minor edits and a different figure 2 from original submission
Submitted: 2002-05-06, last modified: 2002-07-10
The deepest optical image ever in a globular star cluster, a Hubble Space Telescope (HST) 123 orbit exposure in a single field of Messier 4, was obtained in 2 filters (F606W, F814W) over a 10 week period in early 2001. A somewhat shallower image obtained in 1995 allowed us to select out cluster and field objects via their proper motion displacement resulting in remarkably clean color-magnitude diagrams that reach to V = 30, I = 28. The cluster main sequence luminosity function contains very few stars fainter than MV = 15.0, MI = 11.8 which, in both filters, is more than 2 magnitudes brighter than our limit. This is about the faintest luminosity seen among field Population II subdwarfs of the same metallicity. However, there remains a sprinkling of potential cluster stars to lower luminosity all the way down to our limiting magnitudes. These latter objects are significantly redder than any known metal-poor field subdwarf. Comparison with the current generation of theoretical stellar models implies that the masses of the lowest luminosity cluster stars observed are near 0.09 Msun. We derive the mass function of the cluster in our field and find that it is very slowly rising towards the lowest masses with no convincing evidence of a turnover even below 0.1 Msun. The formal slope between 0.65 and 0.09 Msun is alpha = 0.75 (Salpeter of 2.35) with a 99% confidence interval 0.55 - 1.05. A consistency check between these slopes and the number of observed cluster white dwarfs yields a range of possible conclusions, one of which is that we have indeed seen the termination of the white dwarf cooling sequence in M4.
[63]  [pdf] - 49817
The white dwarf cooling sequence in the old open cluster NGC 188
Comments: 6 pages, 5 figures, accepted on Astronomy and Astrophysics
Submitted: 2002-06-12
We develop the white dwarf luminosity function (LF) of the old open cluster NGC 188 in order to determine a lower limit to the age of the cluster by using the faint end of the cooling sequence. To produce an extensive sequence of the cooling white dwarfs we imaged four contiguous HST-WFPC2 fields in the center of the cluster in the F555W and F814W filters. After imposing selection criteria on the detected objects we found a white dwarf cooling sequence (down to V ~26.5) including 28 candidate white dwarfs in the cluster. The exposures are not deep enough to reach the end of this sequence, but the results of our analysis allow us to establish a lower limit to the age of the cluster independently of the isochrone fit to the cluster turnoff. The most ancient white dwarfs found are ~4 Gyr old, an age that is set solely by the photometric limit of our data. Classical methods provide an estimate of \~7 Gyr (Sarajedini et al., 1999).
[64]  [pdf] - 316284
The White Dwarf Cooling Sequence of the Globular Cluster Messier 4
Comments: 14 pages, 4 diagrams. Accepted for publication in the Astrophysical Journal Letters
Submitted: 2002-05-06
We present the white dwarf sequence of the globular cluster M4, based on a 123 orbit Hubble Space Telescope exposure, with limiting magnitude V = 30, I = 28. The white dwarf luminosity function rises sharply for I >25.5, consistent with the behaviour expected for a burst population. The white dwarfs of M4 extend to approximately 2.5 magnitudes fainter than the peak of the local Galactic disk white dwarf luminosity function. This demonstrates a clear and significant age difference between the Galactic disk and the halo globular cluster M4. Using the same standard white dwarf models (Hansen 1999) to fit each luminosity function yields ages of 7.3 +/- 1.5 Gyr for the disk and 12.7 +/- 0.7 Gyr for M4 (2-sigma statistical errors).
[65]  [pdf] - 44867
The CFHT Open Star Cluster Survey. III. The White Dwarf Cooling Age of the Rich Open Star Cluster NGC 2099 (M37)
Comments: 15 Figures and 6 Tables -- 48 pages. Accepted for publication in AJ. Higher resolution images available at
Submitted: 2001-09-20
We present deep CCD photometry of the very rich, intermediate aged (similar to the Hyades) open star cluster NGC 2099 (M37). The V, B-V color-magnitude diagram (CMD) for the cluster shows an extremely well populated and very tightly constrained main-sequence extending over 12 magnitudes from the turn-off. The CFH12K photometry for this cluster is faint enough (V ~ 24.5) to detect the remnants of the most massive progenitor cluster stars under the Type I SNe limit. Therefore, the CMD of the cluster also exhibits a well defined white dwarf `clump' caused by the decreased rate of cooling of these stars as they age, and a subsequent gap with very few objects. The termination point (Mv = 11.95 +/- 0.30) of the white dwarf luminosity function gives a white dwarf cooling age of 566 +/- 154/176 Myrs which is in excellent agreement with the main-sequence turn-off isochrone age (520 Myrs). By carefully accounting for errors, we show that the cooling age confirms that models including convective core overshooting are preferred for young-intermediate aged clusters. We also derive the reddening (E(B-V) = 0.21 +/- 0.03) and distance ((m-M)v = 11.55 +/- 0.13) to NGC 2099 by matching main-sequence features in the cluster to a new fiducial main-sequence for the Hyades. As a continuing part of the goals of the CFHT Open Star Cluster Survey to better understand dynamical processes of open clusters, we also fit a King model to the cluster density distribution and investigate the cluster main-sequence luminosity and mass functions in increasing concentric annuli. We find some evidence for mass segregation within the boundary of NGC 2099 as expected given the cluster's age relative to the dynamical age. The present global mass function for the cluster is found to be shallower than a Salpeter IMF.
[66]  [pdf] - 41902
The CFHT Open Star Cluster Survey I -- Cluster Selection and Data Reduction
Comments: 4 Figures and 3 Tables available in higher resolution at ADS. Accepted for Publication in AJ -- Kalirai et al. 2001a, 122, 257. Updated
Submitted: 2001-04-10, last modified: 2001-07-19
We present this paper in conjuction with the following as the first results in the CFHT Open Star Cluster Survey. This survey is a large BVR imaging data set of 19 open star clusters in our Galaxy. This data set was taken with the CFH12K mosaic CCD (42' X 28') and the majority of the clusters were imaged under excellent photometric, sub-arcsecond seeing conditions. The combination of multiple exposures extending to deep (V ~ 25) magnitudes with short (< 10 second) frames allows for many different studies. The primary aim of this survey is to catalogue the white dwarf stars in these clusters and establish observational constraints on the initial-final mass relationship for these stars and the upper mass limit to white dwarf production. Additionally, we hope to better determine the properties of the clusters, such as age and distance, and also test evolution and dynamical theories by analyzing luminosity and mass functions. In order to more easily incorporate this data in further studies, we have produced a catalogue of positions, magnitudes, colors, and stellarity confidence for all stars in each cluster of the survey. This paper is intended both as a source for the astronomical community to obtain information on the clusters in the survey and as a detailed reference of reduction procedures for further publications of individual clusters. We discuss the methods employed to reduce the data and compute the photometric catalogue. The scientific results for each individual cluster and global results from the study of the entire survey will be presented in further publications (results for the rich old open cluster NGC 6819 immediately follow this paper).
[67]  [pdf] - 41903
The CFHT Open Star Cluster Survey II -- Deep CCD Photometry of the Old Open Star Cluster NGC 6819
Comments: 15 Figures and 6 Tables available in higher resolution from ADS. Accepted for Publication in AJ -- Kalirai et al. 2001b, 122, 266. Updated
Submitted: 2001-04-10, last modified: 2001-07-19
We present analysis of deep CCD photometry for the very rich, old open star cluster NGC 6819. These CFH12K data results represent the first of nineteen open star clusters which were imaged as a part of the CFHT Open Star Cluster Survey. We find a tight, very rich, main-sequence and turn-off consisting of over 2900 cluster stars in the V, B-V color-magnitude diagram (CMD). Main-sequence fitting of the un-evolved cluster stars with the Hyades star cluster yields a distance modulus of (m-M)v = 12.30 +/- 0.12, for a reddening of E(B-V) = 0.10. These values are consistent with a newly calculated theoretical stellar isochrone of age 2.5 Gyrs, which we take to be the age of the cluster. Detailed star counts indicate a much larger cluster extent (R = 9.5' +/- 1.0'), by a factor of ~2 over some previous estimates. Incompleteness tests confirm a slightly negatively sloped luminosity function extending to faint (V ~ 23) magnitudes which is indicative of a dynamically evolved cluster. Further luminosity function and mass segregation tests indicate that low mass objects (M < 0.65Mo) predominate in the outer regions of the cluster, 3.5 < R < 9.5. The estimation of the number of white dwarfs in NGC 6819 are in good agreement with the observed number. For those white dwarf candidates which pass both a statistical and image classification tests, we show comparisons to white dwarf isochrones and cooling models which suggest the need for spectroscopy to confirm the white dwarf nature of the brighter objects.
[68]  [pdf] - 43450
Searching for the Galactic Dark Matter
Comments: To appear in proceedings of ``The Dark Universe: Matter, Energy and Gravity" held at the Space Telescope Science Institute April 2 - 5, 2001
Submitted: 2001-07-04
A straightforward interpretation of the MACHO microlensing results in the direction of the Magellanic Clouds suggests that an important fraction of the baryonic dark matter component of our Galaxy is in the form of old white dwarfs. If correct, this has serious implications for the early generations of stars that formed in the Universe and also on the manner in which galaxies formed and enriched themselves in heavy elements. I examine this scenario in some detail and in particular explore whether the searches currently being carried out to locate local examples of these MACHOs can shed any light at all on this scenario.
[69]  [pdf] - 108009
Isochrones and Luminosity Functions for Old White Dwarfs
Comments: 7 pages, 8 tables, accepted for publication in the Astrophysical Journal
Submitted: 1999-08-30
Using a new grid of models of cooling white dwarfs, we calculate isochrones and luminosity functions in the Johnson-Kron/Cousins and HST filter sets for systems containing old white dwarfs. These new models incorporate a non-grey atmosphere which is necessary to properly describe the effects of molecular opacity at the cool temperatures of old white dwarfs. The various functions calculated and extensively tabulated and plotted are meant to be as utilitarian as possible for observers so all results are listed in quantities that observers will obtain. The tables and plots developed should eventually prove critical in interpreting the results of HST's Advanced Camera observations of the oldest white dwarfs in nearby globular clusters, in understanding the results of searches for old white dwarfs in the Galactic halo, and in determining ages for star clusters of all ages using white dwarfs. As a practical application we demonstrate the use of these results by deriving the white dwarf cooling age of the old Galactic cluster M67.
[70]  [pdf] - 107942
Faint, moving objects in the Hubble Deep Field: components of the dark halo?
Comments: 6 pages, using emulateapj, including 2 colour figures, accepted for publication in ApJ Letters
Submitted: 1999-08-24
The deepest optical image of the sky, the Hubble Deep Field (HDF), obtained with the Hubble Space Telescope (HST) in December 1995, has been compared to a similar image taken in December 1997. Two very faint, blue, isolated and unresolved objects are found to display a substantial apparent proper motion, 23+/-5 mas/yr and 26+/-5 mas/yr; a further three objects at the detection limit of the second epoch observations may also be moving. Galactic structure models predict a general absence of stars in the color-magnitude range in which these objects are found. However, these observations are consistent with recently-developed models of old white dwarfs with hydrogen atmospheres, whose color, contrary to previous expectations, has been shown to be blue. If these apparently moving objects are indeed old white dwarfs with hydrogen atmospheres and masses near 0.5 M_Sun, they have ages of approximately 12 Gyr, and a local mass density that is sufficient, within the large uncertainties arising from the small size of the sample, to account for the entire missing Galactic dynamical mass.
[71]  [pdf] - 107145
If the MACHOs are Old White Dwarfs
Comments: 6 pages, 1 figure, Paper presented at the Frascati Workshop 1999, Multifrequency Behaviour of High Energy Cosmic Sources, May 24 - 29, 1999, Vulcano, Italy
Submitted: 1999-06-26
The microlensing experiments in the direction of the LMC seem to be indicating that about 60% of the dark matter in the Galactic halo is tied up in objects whose masses are about half the mass of the Sun. This mass is a natural one for old white dwarfs although other possibilities do exist. Using a grid of newly constructed models of cooling white dwarfs which incorporate for the first time the effects of molecular opacity in the stellar atmosphere, and assuming that such white dwarfs make up the entire Galactic dark matter, I predict the numbers of old white dwarfs expected in various surveys currently being conducted. In particular, I note the number to be expected in the Hubble Deep Field (HDF), the deepest image of the sky yet obtained.
[72]  [pdf] - 102880
Ages for Globular Clusters in the Outer Galactic Halo: The Second Parameter Clusters Palomar~3, Palomar~4 and Eridanus
Comments: 20 figures, 6 tables, 59 postscript pages
Submitted: 1998-09-14
We have used the WFPC2 camera on the Hubble Space Telescope to obtain photometry of the outer-halo globular clusters Palomar 3, Palomar 4, and Eridanus. These three are classic examples of the "second parameter" anomaly because of their red horizontal-branch morphologies in combination with their low-to-intermediate metallicities. Our color-magnitude diagrams in (V, V-I) reach $V_{lim} \simeq 27.0$, clearly delineating the subgiant and turnoff regions and about three magnitudes of the unevolved main sequences. The slopes and dereddened colors of the giant branches are consistent with published [Fe/H] estimates that rank the clusters (Pal 3, Eridanus, Pal 4) in order of increasing metallicity, with all three falling near or between the abundance values of the classic nearby halo clusters M3 and M5. Differential fits of their color-magnitude diagrams are made to each other and to M3 and M5 for relative age determinations. We find that the three outer-halo cluster CMDs differ from the nearby clusters in a way that is consistent with their being younger by $\sim 1.5 - 2 $Gyr, if we have correctly estimated the clusters' chemical-abundance ratios. Conversely, the inferred age difference could be smaller ($\ltsim 1 $Gyr) if either [Fe/H] or [$\alpha$/Fe] for the outer-halo clusters is significantly lower than we have assumed. Possible age spreads of order 1$ $Gyr among both the nearby and outer-halo clusters may also be present.
[73]  [pdf] - 101765
The White Dwarf Cooling Age of M67
Comments: 15 pages plus 3 diagrams, minor corrections, Accepted for publication in the Astrophysical Journal Letters, to be published September 10, 1998
Submitted: 1998-06-11, last modified: 1998-07-22
A deep imaging survey covering the entire $23\arcmin$ diameter of the old open cluster M67 to $V = 25$ has been carried out using the mosaic imager (UHCam) on the Canada-France-Hawaii Telescope. The cluster color-magnitude diagram (CMD) can be traced from stars on its giant branch at $M_{V} = +1$ down through main sequence stars at least as faint as $M_{V} = 13.5$. Stars this low in luminosity have masses below $0.15 M_{\odot}$. A modest white dwarf (WD) cooling sequence is also observed commencing slightly fainter than $M_V = 10$ and, after correction for background galaxy and stellar field contamination, terminating near $M_V = 14.6$. The observed WDs follow quite closely a theoretical cooling sequence for $0.7 M_{\odot}$ pure carbon core WDs with hydrogen-rich atmospheres (DA WDs). The cooling time to an $M_V$ of 14.6 for such WDs is 4.3 Gyr which we take as the WD cooling age of the cluster. A fit of a set of isochrones to the cluster CMD indicates a turnoff age of 4.0 Gyr. The excellent agreement between these results suggests that ages derived from white dwarf cooling should be considered as reliable as those from other dating techniques. The WDs currently contribute about 9% of the total cluster mass but the number seen appears to be somewhat low when compared with the number of giants observed in the cluster.
[74]  [pdf] - 101404
HST Photometry of the Globular Cluster M4
Comments: 14 pages, 7 figures, accepted for publication in ApJS
Submitted: 1998-05-14
This paper presents a detailed description of the acquisition and processing of a large body of imaging data for three fields in the globular cluster M4 taken with the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope. Analysis with the ALLFRAME package yielded the deepest photometry yet obtained for this cluster. The resulting data-set for 4708 stars (positions and calibrated photometry in V, I, and, in two fields, U) spanning approximately six cluster core radii is available on the AAS CD-ROM (or email a request to RAI). The scientific analysis is deferred to three companion papers, which investigate the significant white dwarf population discovered and the main sequence population.
[75]  [pdf] - 98410
Deep HST V- and I-Band Observations of Two Globular Clusters in the Halo of M31
Comments: 29 pages, 21 Postscript figures, uses aaspp4.sty, to be published in the October 1997 AJ
Submitted: 1997-08-25
We present deep (V ~= 27) V- and I-band stellar photometry of G302 and G312, two globular star clusters in the halo of M31. These data were obtained using the Hubble Space Telescope's Wide Field/Planetary Camera 2. We find iron abundances of [Fe/H] = -1.85 +/- 0.12 for G302 and [Fe/H] = -0.56 +/- 0.03 for G312, consistent with spectroscopic measurements. The color-magnitude diagrams for each cluster show no evidence for an intermediate-aged population of stars, or a second parameter effect in the morphology of the horizontal branch. G302 shows no evidence for a color gradient but the inner regions of G312 are bluer than the outer regions. G312 shows no evidence of ellipticity or an extended halo of unbound stars. G302 has a projected ellipticity of 0.195 +/- 0.012 with the projected major axis oriented towards the center of M31. G302 also shows evidence of an extended asymmetric stellar halo extending to at least twice the fitted Michie-King tidal radius. The amount of mass beyond the tidal radius of G302 is consistent with the stellar escape rates which have been predicted by N-body simulations of globular clusters in the Galactic tidal field.
[76]  [pdf] - 97872
NGC 2419, M92, and the Age Gradient in the Galactic Halo
Comments: 32 pages, aastex, 9 postscript figures; accepted for publication in AJ, September 1997. Also available by e-mail from
Submitted: 1997-07-03
The WFPC2 camera on HST has been used to obtain deep main sequence photometry of the low-metallicity ([Fe/H]=-2.14), outer-halo globular cluster NGC 2419. A differential fit of the NGC 2419 CMD to that of the similarly metal-poor \ standard cluster M92 shows that they have virtually identical principal sequences and thus the same age to well within 1 Gyr. Since other low-metallicity clusters throughout the Milky Way halo have this same age to within the 1-Gyr precision of the differential age technique, we conclude that the earliest star (or globular cluster) formation began at essentially the same time everywhere in the Galactic halo throughout a region now almost 200 kpc in diameter. Thus for the metal-poorest clusters in the halo there is no detectable age gradient with Galactocentric distance. To estimate the absolute age of NGC 2419 and M92, we fit newly computed isochrones transformed through model-atmosphere calculations to the (M_V,V-I) plane, with assumed distance scales that represent the range currently debated in the literature. Unconstrained isochrone fits give M_V(RR) = 0.55 \pm 0.06 and a resulting age of 14 to 15 Gyr. Incorporating the full effects of helium diffusion would further reduce this estimate by about 1 Gyr. A distance scale as bright as M_V(RR) = 0.15 for [Fe/H] = -2, as has recently been reported, would leave several serious problems which have no obvious solution in the context of current stellar models.
[77]  [pdf] - 96809
Formation of the Milky Way
Comments: 8 pages. To appear in Journal of the Korean Astronomical Society. Available also at
Submitted: 1997-03-07
We review observational evidence bearing on the formation of a prototypical large spiral galaxy, the Milky Way. New ground- and space-based studies of globular star clusters and dwarf spheroidal galaxies provide a wealth of information to constrain theories of galaxy formation. It appears likely that the Milky Way formed by a combination of rapid, dissipative collapse and mergers, but the relative contributions of these two mechanisms remain controversial. New evidence, however, indicates that initial star and star cluster formation occurred simultaneously over a volume that presently extends to twice the distance of the Magellanic Clouds.
[78]  [pdf] - 317000
White Dwarfs in Globular Clusters: HST Observations of M4
Comments: 46 pages, latex, no figures included, figures available at size 2.7Mb. To be published in the Astrophysical Journal
Submitted: 1997-02-19
Using WFPC2 on the Hubble Space Telescope, we have isolated a sample of 258 white dwarfs (WDs) in the Galactic globular cluster M4. Fields at three radial distances from the cluster center were observed and sizeable WD populations were found in all three. The location of these WDs in the color-magnitude diagram, their mean mass of 0.51($ \pm 0.03$)M$_{\odot}$, and their luminosity function confirm basic tenets of stellar evolution theory and support the results from current WD cooling theory. The WDs are used to extend the cluster main-sequence mass function upward to stars that have already completed their nuclear evolution. The WD/red dwarf binary frequency in M4 is investigated and found to be at most a few percent of all the main-sequence stars. The most ancient WDs found are about 9 Gyr old, a level which is set solely by the photometric limits of our data. Even though this is less than the age of M4, we discuss how these cooling WDs can eventually be used to check the turnoff ages of globular clusters and hence constrain the age of the Universe.
[79]  [pdf] - 95987
The Initial Mass Function: Now and Then
Comments: 10 pages, postscript. To appear in "Proceedings 7th Annual Astrophysics Conference in Maryland-STAR FORMATION NEAR AND FAR"
Submitted: 1996-11-22
We examine whether existing data in clusters, both old and young, and in the field of the Galactic disk and halo is consistent with a universal slope for the initial mass function (IMF). The most reasonable statement that can be made at the current time is that there is no strong evidence to support a claim of any real variations in this slope. If the IMF slope is universal then this in itself is remarkable implying that variations in metallicity, gas density or other environmental factors in the star formation process play no part in determining the slope of the mass function.
[80]  [pdf] - 94838
Deep HST V- and I-Band Observations of the Halo of M31: Evidence for Multiple Stellar Populations
Comments: 31 pages AASTeX v4.0, 10 Figures avaliable at
Submitted: 1996-06-12
We present deep V & I WFPC2 photometry in two fields in the M31 halo 32' and 50' from the center of M31 along the SE minor axis at the locations of the M31 GCs G302 and G312 respectively. The M31 halo LFs are not consistent with a single high-metallicity population but are consistent with a mix of 50% to 75% metal-rich stars and 25% to 50% metal-poor stars. This agrees with the RGB morphology, the HB luminosity, and the RR Lyrae stars in the M31 halo. The RGB morphology indicates a metallicity spread of -2 < [m/H] < -0.2 with most stars having [m/H] = -0.6, making the M31 halo significantly more metal- rich than either the Galactic halo or the M31 GC system. The HB is dominated by a red clump similar to the 47 Tuc HB but a small number of blue HB stars are visible, supporting the conclusion that there is a metal-poor component to the M31 halo. The number of HB stars is smaller than would be expected from the observed metallicity distribution but it is not clear if this is due to the photometric limits of our data or a 2nd parameter effect. We find a He abundance of Y >~ 0.20 to 0.27, comparable with the Galaxy. LFs show weak evidence that the R = 50' field contains a higher fraction of metal-poor stars than the R = 32' field but the metallicity distributions of the RGB stars strongly suggest that both fields have the same mix of stellar populations.
[81]  [pdf] - 94583
Globular Clusters in NGC 1275
Comments: 25 pages, latex, postscript figures, tarred, Unix compressed, postscript version of paper and figures available at
Submitted: 1996-05-07
We present the results of a deep photometric study of the outer halo of NGC 1275, the highly active cD galaxy at the center of the Perseus cluster. We find a modest excess of faint ($R > 22.5$) starlike objects in its halo, indicating a population of old-halo globular clusters. However, the total estimated cluster population corresponds to a specific frequency of $S_N = 4.4 \pm 1.2$, no larger than that of normal giant ellipticals and three times lower than that of other central cD galaxies such as M87. We discuss several ideas for the origin of this galaxy. Our results reinforce the view that high $S_N$ (ie: highly efficient globular cluster formation) is not associated with cooling flows, or with recent starburst or merger phenomena.
[82]  [pdf] - 93927
Age Gradient and the Second Parameter Problem in the Galactic Halo
Comments: 26 pages, AASTeX v3.0 with 3 Postscript figures
Submitted: 1996-01-12
We establish a framework for determing absolute ages of Galactic globular clusters and then use these ages to investigate the age-metallicity and age-Galactocentric distance relations for the 36 clusters with the most reliable age data. The clusters span Galactocentric distances from 4 through 100 kpc and cover a metallicity range from $[Fe/H] = -0.6$ to $-2.3$. Adopting currently plausible choices for the relation between cluster metallicity and horizontal-branch luminosity, and alpha-enhancement ratios, we find that the majority of the globular clusters form an age distribution with a dispersion $\sigma(t)$ about $10^9$ years, and a total age spread smaller than 4 Gyr. Clusters in the lowest metallicity group ($[Fe/H] < -1.8$) appear to be the same age to well within 1 Gyr at all locations in the Milky Way halo, suggesting that star formation began throughout the halo nearly simultaneously in its earliest stages. We find no statistically significant correlation between mean cluster age and Galactocentric distance (no age gradient) from 4 to 100 kpc. The correlation between cluster age and horizontal-branch type suggests that causes in addition to metallicity and age are required to understand the distribution of stars along the horizontal branches in globular cluster color-magnitude diagrams.
[83]  [pdf] - 93047
Hubble Space Telescope Observations of White Dwarfs in the Globular Cluster M4
Comments: To appear in the Astrophysical Journal Letters, 451, September 20, 1995. The text is also available as a postscript file from with postscript files of the figures available in
Submitted: 1995-07-26
With the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope, we have discovered in M4 (NGC 6121, C 1620-264) the first extensive sequence of cooling white dwarfs seen in a globular cluster. Adopting a distance modulus of (m-M)_V = 12.65 and a reddening of E(B-V) = 0.37, we show that the sequence, which extends over 9 < M_U < 13, is comprised of white dwarfs of mass \sim 0.5 M_{\odot}. The total mass loss from the present turnoff to the white dwarf sequence is 0.31 M_{\odot} and the intrinsic dispersion in the mean mass appears to be < 0.05 M_{\odot}. Both the location of the white dwarf cooling sequence in the cluster color-magnitude diagram and the cumulative luminosity function attest to the basic correctness and completeness of the physics in theoretical models for the upper three magnitudes of the observed white dwarf cooling sequence. To test the theory in globular clusters at cooling ages beyond \sim 3 \times 10^8 years will require deeper and more complete data.
[84]  [pdf] - 92843
Counting Pairs of Faint Galaxies
Comments: 471 kb uuencoded, postscript file with figures included, 32 pages, To appear in the Nov., 1995 issue of Astrophysical Journal
Submitted: 1995-06-07
The number of close pairs of galaxies observed to faint magnitude limits, when compared to nearby samples, determines the interaction or merger rate as a function of redshift. The prevalence of mergers at intermediate redshifts is fundamental to understanding how galaxies evolve and the relative population of galaxy types. Mergers have been used to explain the excess of galaxies in faint blue counts above the numbers expected from no-evolution models. Using deep CFHT ($I\leq24$) imaging of a ``blank'' field we find a pair fraction which is consistent with the galaxies in our sample being randomly distributed with no significant excess of ``physical'' close pairs. This is contrary to the pair fraction of $34\%\pm9\%$ found by Burkey {\it et al.} for similar magnitude limits and using an identical approach to the pair analysis. Various reasons for this discrepancy are discussed. Colors and morphologies of our close pairs are consistent with the bulk of them being random superpositions although, as indicators of interaction, these criteria are uncertain due to contamination of field galaxies and difficulty in seeing structure at intermediate redshifts, respectively. This observed lack of close pairs is probably linked with the decreasing amplitude of the angular correlation function at faint limits. If our faint samples are comprised of galaxies which have a higher average redshift than brighter samples studied by other workers then either the merger rate has been overestimated or there is a change in its behaviour from what is observed at brighter magnitude limits.
[85]  [pdf] - 1469242
Late-Type Stars in M31. I.: A Photometric Study of AGB Stars and Metallicity Gradients.
Comments: UUencoded compressed postscript, 3 Figs. available on request. (Contact
Submitted: 1995-03-10
We have imaged five 7\arcmin \x 7\arcmin\ fields in M31 spanning galactocentric radii from 4 to 32 kpc along the SW-major axis. The fields were observed through two broad-band (\V\ and \I) and two narrow-band (\CN\ and \TiO) filters. The broad-band data were used to construct \IvsVI\ color-magnitude diagrams (CMDs) and, in some of our fields, we found significant numbers of stars in the Cepheid instability strip. A distance modulus for the Cepheids in the middle field was found that agreed well with other values in the literature values. The width of the giant branch (GB) in the \IvsVI\ CMD of all 5 fields was investigated, and we show that in four of the fields a likely explanation for the GB width is a combination of {\it both} metallicity and mass variations. Using the broad-band data, the asymptotic giant branch (AGB) luminosity functions (LFs) were measured in the five fields, and we show that differences exist between these LFs. We speculate on how the different star forming histories in the fields may lead to the observed AGB LFs and GB widths. Using the narrow-band data along with the broad-band data we separated the AGB stars into carbon-rich (C) and oxygen-rich (M) types. The carbon stars LFs were used to obtain an estimate for the distance modulus of M31 which agrees with the value derived from Cepheids. The ratio of C- to M-stars (C/M) is believed to be an indicator of gaseous chemical abundance at the time of formation of these stars. We show that the C/M ratio increases smoothly with galactocentric distance, suggesting an inverse correlation with metallicity. This is the first demonstration of this effect within a single extragalactic system. We find that differences in the width of the GB and the AGB LFs do not significantly affect the C/M ratio. We consider the effect of the increasing C/M ratio on the ISM in M31, and cite evidence in favor of a model where the grain composition in M31 is a function of galactocentric distance.
[86]  [pdf] - 92370
The Double Cluster G185 in M31
Comments: 22 pages, ~1Mb postscript file
Submitted: 1995-02-10
We have identified a small globular cluster in M31 located approximately 4 arcseconds northwest of the M31 globular cluster G185. While several multiple globular clusters have been observed in the Magellanic Clouds none have been found in the Galaxy or in M31. We estimate the probability of such a chance line-of-sight alignment occuring near the nucleus of M31 to be 0.09 +/- 0.03 and find no obvious indication of any tidal deformation in either cluster, as would be expected if the clusters were interacting. Two-dimensional modelling suggests G185 has a King (1966) [AJ, 71, 64] concentration of c = 1.11 +/- 0.08 while the companion has c = 0.67 +/- 0.17 and is physically smaller than G185. Both objects have integrated dereddened colors similar to those of Galactic globular clusters.
[87]  [pdf] - 90918
Star counts in NGC 6397
Comments: 16 pages + 7 pages of tables, LaTeX using AASTeX macros, 11 figures available by request, IoA preprint
Submitted: 1993-09-01
I-band CCD images of a large area of the nearby globular cluster NGC~6397 have been used to construct a surface density profile and two luminosity and mass functions. The surface density profile extends out to 14\arcm from the cluster center and shows no sign of a tidal cutoff. The inner profile is a power-law with slope -0.8 steepening to -1.7 outside of 1\arcm. The mass functions are for fields at 4\arcm\ and 11\arcm from the cluster center and confirm the upturn in the mass function for stars less massive than about 0.4 M\solar. There appears to be an excess of low-mass stars over higher-mass stars in the outer field with respect to the inner, in qualitative agreement with expectations for mass segregation.