Applications of High-Resolution Astrometry to Galactic Studies
DISSERTATION
Presented in Partial Fulfillment of the Requirements for
the Degree Doctor of Philosophy in the
Graduate School of The Ohio State University
By
Samir Salim, M.S.
*****
The Ohio State University
2002
Dissertation Committee: Approved by
Professor Andrew Gould, Adviser
Professor Marc Pinsonneault Advisor Astronomy Graduate Program Professor David Weinberg
UMI Number: 3062651
______UMI Microform 3062651 Copyright 2002 by ProQuest Information and Learning Company. All rights reserved. This microform edition is protected against unauthorized copying under Title 17, United States Code. ______
ProQuest Information and Learning Company 300 North Zeeb Road PO Box 1346 Ann Arbor, MI 48106-1346 ABSTRACT
Astrometry is undergoing a revolution that started with Hipparcos satellite and is continuing with future missions like DIVA, FAME, GAIA or SIM;aswell as novel ground-based techniques. This dissertation investigates in detail several applications of high-resolution astrometry to studies of the Galaxy.
For many of the studies presented here a catalog of high proper motion stars is required. The largest such catalog, Luyten’s NLTT, constructed several decades ago, is of limited usefulness in its original form. We therefore construct its refined version, containing improved astrometry and photometry for the vast majority of the ∼ 59, 000 stars in NLTT. The bright end is constructed by matching NLTT to Hipparcos, Tycho-2, and Starnet; the faint end by matching to USNO-A and
2MASS. We improve Luyten’s 6 positions to better than 0 .1. Proper-motion errors in NLTT (∼ 25 mas yr−1) are improved to 5.5masyr−1. Further, with our improved optical/infrared photometry we construct a reduced proper motion (RPM) diagram that, in contrast to the original NLTT RPM diagram, allows for the first time the classification of NLTT stars into main-sequence (MS) stars, subdwarfs (SDs), and
ii white dwarfs (WDs). We use this diagram to analyze the properties of our revised, and of the original NLTT. We also produce a list of new candidate nearby WDs.
We use improved NLTT in selection of nearby microlensing candidates. SIM can measure a minute deflection in the source’s apparent position, and so provide a precise (1%) mass determination of the nearby lens star (in many cases a SD).
We search for lens-source encounters using Hipparcos, ACT and NLTT to select lenses, and USNO-A to select sources. Among 32 candidates from Hipparcos there are Proxima Centauri and Barnard’s star. For NLTT lenses, the distance and the impact parameter are more poorly known, leading to large uncertainties in the amount of SIM observing time. However, using single-epoch CCD observations of the candidates, and the information from the revised NLTT, we have considerably reduced the uncertainties and produced a reliable list of targets.
Three planned astrometry survey satellites, FAME, DIVA,andGAIA,all aim at observing magnitude-limited samples. We argue that substantial additional scientific opportunities are within a reach if a limited number of fainter targets is included. Thus, we can increase the number of late-M dwarfs, L dwarfs, and WDs with good parallaxes by an order of magnitude, and enable good determinations of local mass functions (MF). In most cases, the candidate dim dwarfs are not yet known, and we present various methods to identify them. The presented analysis applies to DIVA as well.
iii By combining SIM observations with ground-based photometry, one can completely solve microlensing events seen toward the Galactic bulge, allowing one to measure the mass, distance, and transverse velocity of ∼ 100 lenses to ∼ 5% precision. This would allow the MF of the bulge objects, both luminous and non-luminous (remnants), to be measured.
As an application of ground-based astrometry, we present a new geometrical method for measuring the distance to the Galactic center (R0) by solving for the orbit of individual stars bound to the Sgr A*. We identify stars to which the method may be applied, and show that two of them could produce 1-5% accuracy of
R0 after 15 years of positional and radial velocity measurements. Further, we show that combining the measurements of the three stars, while the common center and mass are constrained, produces considerable improvements in the R0 determination.
By modeling the probability distribution over allowed orbital parameter space, we
find that by 2010 the achievement of 3% precision is almost certain, while there is a 30% chance of obtaining 1% precision. These estimates would improve if the astrometry errors are reduced from the current 2 mas.
For studying the stellar halo we propose observing proper motions of faint horizontal branch stars with FAME. Using them as standard candles, halo rotation can be precisely (2 km s−1) mapped out to 25 kpc, and the clumps in the kinematic space (halo substructure) detected. Finally, we propose using SIM astrometric
iv microlensing (in a fashion similar to that applied to get bulge MF) to shed light on the nature of MACHO objects and their relation to the dark halo. By simultaneous observations of microlensing events towards the LMC from SIM and the Earth, we expect to measure distances and masses of these objects.
v To my parents Mirjana and Yahya
vi ACKNOWLEDGMENTS
I wish to express my gratitude to Andrew Gould for being a truly great adviser. I am also thankful to everyone at the Department of Astronomy for always being helpful, and for creating an excellent research environment.
This work was supported by JPL contract 1226901, grant AST 97-27520 from the NSF, grant NAG5-3111 from NASA, and Program for the Enhancement of
Graduate Studies (PEGS) from The Ohio State University.
I thank Staˇsa for her love, which is the greatest support one can have.
vii VITA
May 16, 1971 ...... Born– Belgrade,Yugoslavia
1996 ...... B.S.Astrophysics,UniversityofBelgrade
2000 ...... M.S.Astronomy,TheOhioStateUniversity
1997 – 1998 ...... UniversityFellow,TheOhioStateUniversity
1998 – 2002 ...... GraduateTeachingandResearch Associate,
The Ohio State University
PUBLICATIONS
Research Publications
1. S. Salim and A. Gould, “Sagittarius A* “Visual Binaries”: A Direct Measurement of the Galactocentric Distance.” Astroph. J., 523, 633, (1999).
2. A. Gould and S. Salim, “Photometric Microlens Parallaxes with the Space Interferometry Mission.” Astroph. J., 524, 794, (1999).
3. S. Salim and A. Gould, “Nearby Microlensing Events: Identification of the Candidates for the Space Interferometry Mission.” Astroph. J., 539, 241, (2000).
4. I. Ferrin et al.“Discovery of the Bright Trans-Neptunian Object 2000
EB173.” Astroph. J., 548, L243, (2001).
viii 5. C. Flynn, J. Sommer-Larsen, B. Fuchs, D. S. Graff and S. Salim, “A Search for Nearby Counterparts to the Moving Objects in the Hubble Deep Field.” Monthly Notices of the Royal Astronomical Society, 322, 533, (2001).
6. Z. Zheng, C. Flynn, A. Gould, J. N. Bahcall and S. Salim, “M Dwarfs from Hubble Space Telescope Counts. IV.” Astroph. J., 555, 393, (2001).
7. A. Gould and S. Salim, “Searching for Failed Supernovae With Astrometric Binaries.” Astroph. J., 572, 944, (2002).
8. S. Salim, A. Gould and R. P. Olling, “Astrometry Survey Missions Beyond the Magnitude Limit.” Astroph. J., 573, 631, (2002).
9. S. Salim and A. Gould, “Classifying Luyten Stars Using an Optical- Infrared Reduced Proper-Motion Diagram” Astroph. J., 575, L86, (2002).
FIELDS OF STUDY
Major Field: Astronomy
ix Table of Contents
Abstract...... ii
Dedication...... vi
Acknowledgments...... vii
Vita...... viii
ListofTables...... xvi
ListofFigures...... xxi
1 Introduction 1
1.1AstrometryandAstronomy...... 1
1.2TheHipparcosRevolution...... 7
1.3FutureAstrometryMissions...... 10
x 1.4 Future of Ground-based Astrometry? ...... 13
2 Refinement and Analysis of Luyten’s Catalog of High Proper
Motion Stars 16
2.1Introduction ...... 16
2.2BrightStarsinNLTT...... 23
2.2.1 Positional Errors of Bright NLTT Stars ...... 23
2.2.2 Strategy to Match Bright NLTT Stars ...... 28
2.2.3 Proper Motion Errors of Bright NLTT Stars ...... 33
2.3FaintStarsinNLTT...... 36
2.3.1 Strategy for Matching NLTT to USNO-A and 2MASS . . . . 36
2.3.2 AdditionalMatches ...... 47
2.3.3 Positional Errors of Faint NLTT Stars ...... 53
2.3.4 Proper Motion Errors of NLTT Stars ...... 55
2.3.5 CommonProperMotionBinaries...... 56
2.4RevisedNLTTCatalog ...... 58
xi 2.4.1 FormatoftheRevisedNLTT...... 58
2.4.2 Proper Motion Errors of the Revised NLTT ...... 63
2.5ClassificationofNLTTStars ...... 66
2.5.1 Classification with the Original NLTT ...... 66
2.5.2 ClassificationwiththeRevisedNLTT ...... 68
2.5.3 CandidateNearbyWhiteDwarfs...... 69
2.6CompletenessofNLTT ...... 70
2.6.1 Completeness at Bright Magnitudes ...... 70
2.6.2 Bright-end Completeness as a Function of Galactic
Coordinates...... 73
2.6.3 CompletenessatFaintMagnitudes...... 75
2.6.4 Completeness of Different Star Populations as a Function of
GalacticLatitude ...... 79
2.7DiscussionandConclusion...... 84
3 Precise Masses of Nearby Stars 111
xii 3.1Introduction ...... 111
3.2 Astrometric Microlensing with Space Interferometry Mission . . . . . 112
3.3SelectionofMicrolensingCandidates...... 114
3.3.1 PrinciplesofSelection...... 114
3.3.2 CatalogsUsedinSelection ...... 116
3.3.3 EstimatesofErrors ...... 123
3.3.4 SearchingfortheCandidateEvents ...... 134
3.3.5 CandidateEvents ...... 138
3.4ConfirmationofCandidateEvents ...... 143
3.4.1 RefiningEstimatesUsing2MASS...... 143
3.4.2 ObservingCampaign...... 144
3.5DiscussionandConclusion...... 146
4 Mass Function of Stellar-Mass Objects 183
4.1 Mass Function in the Solar Neighborhood ...... 183
4.1.1 Late-MandLdwarfs ...... 183
xiii 4.1.2 WhiteDwarfs ...... 187
4.2BulgeMassFunction...... 191
4.2.1 Bulge Mass Function and the Microlensing ...... 191
4.2.2 DegeneracyofMicrolensingEvents...... 192
4.2.3 Breaking the Degeneracy with Astrometric Microlensing . . . 196
4.2.4 Addition of Photometric Observations ...... 199
4.2.5 Simulations of SIM Observations ...... 203
4.2.6 Expected Precision of Mass Measurements ...... 213
5 Precise Distance to the Galactic Center 226
5.1Introduction ...... 226
5.1.1 Past Determinations of the Galactocentric Distance ...... 226
5.1.2 “VisualBinaries”AroundSgrA*...... 229
5.2PhysicalPrinciplesoftheMethod ...... 231
5.3MethodofErrorDetermination...... 234
xiv 5.4 Predictions of the Uncertainty of Galactocentric Distance
Determination ...... 238
5.4.1 IndividualStars ...... 238
5.4.2 CombinedSolution...... 242
5.4.3 Probability of Achieving a Given Precision ...... 243
5.5DiscussionandConclusion...... 247
6 Kinematics of Stellar and Nature of Dark Halo 258
6.1KinematicsoftheStellarHalo ...... 258
6.1.1 Introduction ...... 258
6.1.2 SelectingFaintA-typeStars ...... 259
6.1.3 HaloRotationMeasurement ...... 263
6.1.4 Substructure in the Galactic Halo ...... 267
6.2NatureoftheDarkMatterHalo ...... 269
6.2.1 Introduction ...... 269
xv 6.2.2 MACHO Masses and Distances from Astrometric
Microlensing ...... 270
6.2.3 Expected Precision of Distance and Mass Measurements . . . 272
Bibliography...... 284
xvi List of Tables
2.1BrightNLTTstarspositionprecisions ...... 86
2.2 Bright NLTT stars proper motion precisions ...... 87
2.3NumberofUSNO/2MASSmatches...... 88
2.4Newcandidatewhitedwarfscloserthan20pc...... 89
3.1 Hipparcos and ACT events - lens star properties (astrometry and
photometry) ...... 148
3.2 Hipparcos and ACT events - lens star properties (distance, kinematic
andphysicalproperties)...... 151
3.3 Hipparcos and ACT events - source star and event properties . . . . 154
3.4NLTTevents-lensstarproperties...... 157
3.5NLTT-sourcestarandeventproperties...... 167
xvii 4.1 Uncertainties of microlensing parameters, Bulge (I = 15, ψ = 225◦,
5hours)...... 219
6.1RotationofthehalofromBHBstars...... 276
6.2 Uncertainties of microlensing parameters, LMC (V = 20, arbitrary
ψ,20hours.)...... 277
xviii List of Figures
1.1 Astrometric accuracies of different missions ...... 15
2.1 Positional errors of bright stars in NLTT – a close-up ...... 91
2.2PositionalerrorsofbrightstarsinNLTT...... 92
2.3 Differences between photographic photometry from NLTT and USNO 93
2.4PositionalerrorsoffaintstarsinNLTT...... 94
2.5 Fraction of NLTT stars found in a rectangle as a function of V .... 95
2.6 Fraction of NLTT stars found in a rectangle as a function of declination 96
2.7NLTTpropermotionerrors...... 97
2.8 Difference between separation vectors of CPM binaries in NLTT and
theRevisedcatalog...... 98
xix 2.9 Difference in relative proper motions for NLTT CPM binaries . . . . 99
2.10 Differences between proper motions as given by Tycho-2 and our
measurements based on identifying the stars in USNO and 2MASS . . 100
2.11 Difference in the relative proper motion of the components of
subdwarfbinaries...... 101
2.12 Original reduced proper motion (RPM) diagram for NLTT stars . . . 102
2.13 Optical-infrared RPM diagram for NLTT stars ...... 103
2.14 Original NLTT RPM diagrams with stars classified using optical-
infraredRPMdiagram...... 104
2.15 Completeness of the combined Hipparcos and Tycho-2 catalogs with
respect to NLTT, as a function of RNLTT ...... 105
2.16 Completeness of bright-end NLTT with respect to Hipparcos/Tycho-
2,asafunctionofpropermotion...... 106
2.17 Completeness of bright-end NLTT as a function of Galactic latitude . 107
2.18 Completeness of bright-end NLTT as a function of Galactic longitude 108
2.19 Completeness of faint-end NLTT as a function of magnitude . . . . . 109
xx 2.20 Surface density of three different stellar populations as a function of
Galacticlatitude...... 110
3.1Mass-luminosityrelation...... 177
3.2 Reduced proper motion diagram for Hipparcos and NLTT stars . . . 178
3.3 Distance-modulus errors vs. distance modulus for Hipparcos stars . . 179
3.4 Hipparcos color-magnitude diagram with candidate events ...... 180
3.561CygA/Bmicrolensingevents...... 181
3.6 Apparent magnitude – proper motion diagram of candidate events . . 182
4.1 Histogram of late-type dwarfs with precise parallax estimates using
FAME ...... 224
4.2 Distribution over MV of white dwarfs with precise parallax estimates
using FAME ...... 225
5.1 Evolution of R0 fractional uncertainty from S0-1 orbits ...... 250
5.2 Evolution of σR0 from S0-1 orbit with different values of σvrad .....251
5.3 Contour plot of σR0 from S0-1 orbit (tobs =15yr)...... 252
xxi 5.4 Contour plot of σR0 from S0-1 orbit (tobs =30yr)...... 253
5.5 Contour plot of σR0 from S0-2 orbit (tobs =15yr)...... 254
5.6 Contour plot of σR0 from S0-2 orbit (tobs =30yr)...... 255
5.7 Evolution of σR0 inacombinedsolution...... 256
5.8 Probability of obtaining specific σR0 ...... 257
6.1PropermotiondistributionofA-typestars...... 282
6.2 Estimate of the measurement error of halo rotation ...... 283
xxii Chapter 1
Introduction
1.1. Astrometry and Astronomy
The Merriam-Webster Dictionary defines “astrometry” as: “a branch of astronomy that deals with measurements (as of positions and movements) of celestial bodies.” It also lists that the term first appears around 1859. From this definition it is clear that astrometry is just one of the branches of astronomy, but its relation to astronomy and the role it plays in it, and also the way this relation and role have changed throughout the history requires further elaboration.
We can say that to a great extent, from its earliest beginnings until the late
19th century astronomy was astrometry. The first time some ancient astronomer looked at the sky and recorded relative positions of stars in the form of a constellation, he or she was already making an astrometric measurement, albeit
1 a very crude one. And while there was not much to record in terms of stellar motions (except for the diurnal motion), following the apparently complex motions of planets prompted the first sound ideas about the physical underpinnings of the universe. Actually it was the lack of any observable relative motion of stars that imposed constraints on these early models of the universe – namely the lack of parallax was one of the key arguments that the Earth was motionless, and as such occupying the central position in the universe.
Despite the fact that the change of planetary positions was more interesting to observe, it was soon realized that in order to determine these changes accurately enough, there must also exist a good catalog of stellar positions. Most likely using some very simple apparatus for measuring relative stellar angles and the altitudes of stars passing the meridian, Hipparchus, in the 2nd century BC, was able to compile a catalog of stellar positions and apparent magnitudes of 850 naked-eye stars, with a typical accuracy of half a degree. The original catalog of Hipparchus was included in Ptolemey’s “Alamgest” three centuries later, and possibly expanded with some of his own measurements. This catalog became a standard for centuries to come.
The next major advancement in astronomy (or astrometry) that had profound consequences for the science as a whole (and actually for laying the foundations of modern science), came with Tycho Brahe in the 16th century, and his precise measurements of stellar and planetary positions. Using large instruments for
2 measuring angles, he was able to achieve arcminute precision – the limit of naked-eye observations. His observations of some 800 stars were compiled by
Kepler. The level of accuracy allowed Kepler to challenge the assumption of circular planetary orbits, and ultimately to formulate his laws of planetary motion, which paved the way for one of the greatest achievements in the history of science
– Newton’s theory of universal gravity.
With the introduction of telescopes, more and more of astronomical observations had aspects different from those of astrometry – telescopes allowed for the exploration of physical features of the celestial bodies, thus setting the stage for what we today call astrophysics, which practically became synonymous with astronomy. The introduction of telescopes meant further advances in astrometry as well, since more precise astrometric measurements become possible. Armed with a telescope and the Newton’s theory of gravity, a new field emerged – celestial mechanics – determining orbits of Solar System bodies and predicting their future positions. The successes of the celestial mechanics were demonstrated early on by the discovery of the periodicity of the Halley’s comet and ultimately by the discovery of Neptune through gravitational perturbations on Uranus, by Le Verrier and Adams in mid 19th century.
The aspects of astrometry related to stars also came to flourish. For the first time it was possible to measure a star’s own (proper) motion – in 1718 Halley found
3 that Arcturus, Procyon and Sirius could no longer be considered “fixed” stars. The
Earth’s motion was reflected in the abberation of light, discovered by Bradley, and more importantly, and what Bradley actually attempted to measure but it took a century before instruments became precise enough to do it, the parallax of another star was detected in 1838 by Bessell, thus obtaining the first distance to any star besides the Sun, and also ultimately proving that Earth orbits the Sun. This was also how astrometry touched upon the emerging field of astrophysics – measuring stellar distances enabled their physical brightness to be known and compared to the
Sun’s, leading to the concept of HR diagram, the main sequence, and giving clues to stellar evolution. Proper motions coupled with distances produce two components of the physical velocity, enabling studies of the kinematics and dynamics of the
Solar neighborhood, and later, of the Galaxy as a whole. Bessel is also credited with detecting wobbles in Sirius and Procyon, produced by companions, which at that time escaped optical detection .
Besides this aspect of high-precision (sub-arcsecond) astrometry, the charting of the skies that included more and more stars continued. Flamsteed produced a catalog of 3000 stars in the early 18th century (including Halley’s observations of stars in the southern hemisphere). This catalog reached a precision of 10 . Around the mid-19th century, Argelander, after working for 25 years, produced a massive
4 catalog and charts containing 300,000 northern stars to 10th magnitude. This work was later expanded to the south.
At approximately the same time the photographic plates came to be used in astronomy. Now positions of many stars could be measured on a single plate, and with yet greater precision. The realization of this capability of astrophotography led to the creation of a huge international program of producing a photographic atlas of the sky at a scale good enough for astrometric measurements, known as Carte du
Ciel. Although the work was divided among many observatories (mostly in Europe and some in South Africa, South America and Australia), the project proved to be overwhelming, not only for its technical requirements where the quality and the number of plates was concerned (each was 2◦ square), but even more so when it came to measuring the positions of the stars from the plates to produce the accompanying catalog (Astrographic Catalog). The atlas aimed at reaching 14th magnitude, and the catalog which was never completed, was supposed to include stars to 12th magnitude. It is widely believed that the commitment of European astronomers to this burdensome project led to their falling behind their American colleagues, who at that time concentrated on the studies of extragalactic “nebulae” that resulted in one of the most profound discoveries, that of the expanding
Universe and the Big Bang. It was a century after its conception that the project of Carte du Ciel reached its conclusion.
5 Later in the 20th century more efficient methods of photographing large swaths of the sky became possible with the introduction of Schmidt cameras and their large 6◦ field of view. The Palomar Observatory Sky Survey released in 1958 was a treasure trove for discoveries of every kind, from local asteroids to distant galaxy clusters. It was this kind of plates that William Luyten used in his half a century long effort to catalog all stars with high proper motion. He achieved this by blinking plates taken at different epochs, either manually or using simple plate measuring machines. This endeavor resulted in the catalog of nearly 60,000 stars moving faster than 0 .2yr−1, known as the New Luyten Two-Tenths (NLTT) catalog. In many aspects, such as depth, completeness and overall accuracy, this catalog is still unsurpassed. We found this catalog essential for many of the studies presented in this dissertation and that is why we pay special attention to it, and made considerable effort to bring this catalog up-to-date as to make it more useful for the variety of studies. Chapter 2 is thus completely devoted to this catalog.
Besides inability to produce deeper and more precise all-sky catalogs of stellar positions (and proper motions), difficulties also plagued further improvements in narrow-angle astrometry, which is responsible for yielding stellar parallaxes.
It is not only that the observations where technically difficult, but the blurring of images because of the Earth’s atmosphere set the limit to the accuracy of individual astrometric measurement and also to the ultimate achievable accuracy,
6 thus limiting the volume of space that could be explored. As late as in the 1980s, the most sophisticated astrometric observations augmented by CCDs could at best reach a parallax precision of a few milli-arcseconds (mas), but even this for a very limited number of stars.
1.2. The Hipparcos Revolution
Close to the end of the 20th century it appeared that astrometry, and its role in astronomy was all but dead. However, this condition was radically changed with the release of results from the Hipparcos astrometric satellite in 1997. Launched in
1989, the satellite was not able to reach the correct orbit and it appeared that the entire mission would fail. However, useful data where obtained in this faulty orbit as well. Without blurring from the Earth’s atmosphere it was possible to achieve very high precision even with small aperture of Hipparcos telescopes by repeatedly making hundreds of positional measurements of the same star. The satellite carried two instruments. The main instrument performed observations of some 120,000 preselected stars, achieving the average positional, parallax and proper motion precision of 1 mas. This sample (approximately one half of the total number) was magnitude-complete to V ∼ 8. Of the fainter stars, to the limiting magnitude of the instrument of V ∼ 11.5, only some where observed – mainly high proper motion stars drawn from NLTT, but also many astrometric standards for which
7 improved positions were desirable, and even stars that are being occulted by the
Moon (presumably in order to use occultation timings to refine the lunar orbit.)
This main catalog is known simply as the Hipparcos catalog (ESA 1997). The second instrument carried out measurements of all stars down to the magnitude limit, but with a 20 times lower precision. This is the Tycho catalog (ESA 1997), and it contains around 1 million stars.
Precise parallaxes and proper motions of tens of thousands of stars in the
Hipparcos catalog enabled many different studies. The HR diagram was calibrated with the unprecedented precision. The disputed distance to the Hyades was measured directly, leading to new estimates of the luminosity of stellar candles and redefining most of the cosmic distance ladder. Proper motions provided more secure estimates of Sun’s motion with respect to the local standard of rest, identified clumps of stars in the velocity space, and produced estimates of the local mass density. Many more important results stemming from Hipparcos catalog could be mentioned.
The larger, but less precise Tycho catalog also proved to be quite valuable.
Namely, it was realized that the density of Tycho stars is high enough to enable precise astrometric calibration of wide field photographic plates, like that of POSS, which are otherwise difficult to calibrate because of severe distortions. However, for first generation POSS plates, taken in the 1950s, the 20 mas Tycho positions
8 obtained in 1991 would degrade to several arcsec because of imprecise proper motions. This is where the material taken for Carte du Ciel Astrographic Catalog came to the rescue. With modern computer technology it was relatively easy to reduce these old plates. With the mean epoch of 1907, the Astrographic Catalog
(AC, Urban et al. 1998a) served as a perfect first epoch to compare with Tycho and derive proper motions with accuracy of 4 mas yr−1. This combination of AC and Tycho catalogs came to be known as the ACT catalog (Urban, Corbin, &
Wycoff 1998b). Later, the data from the Tycho experiment were rereduced to reach 0.5 mag fainter limit, doubling the number of stars, and again these positions were combined with AC and approximately 150 other smaller catalogs, to produce
Tycho-2 (Høg et al. 2000). With ACT stars serving as positional references, POSS-I and southern ESO/UKIRT plates were calibrated to produce the largest stellar catalog ever, USNO-A (Monet 1998), containing 5 × 108 stars. ACT also served as a basis for astrometric calibration of the two catalogs derived from truly digital sky surveys: optical SDSS (York et al. 2000), and infrared 2MASS (Skrutskie et al. 1997). We used USNO-A and 2MASS extensively to refine astrometry and photometry of NLTT stars.
9 1.3. Future Astrometry Missions
Building on the success of Hipparcos, new astrometric missions having greater accuracy and sensitivity were proposed. In the first group of these new missions we can include two missions aiming to achieve accuracies an order of magnitude better than that of Hipparcos – Full-Sky Astrometric Explorer (FAME) and DIVA.Both should be complete to R ∼ 15, totaling some 4 × 107 stars. The main expected scientific returns are that the primary distance calibrators such as Cepheids, RR
Lyrae and low-metallicity MS stars will be calibrated directly, and with a great precision. Astrometric monitoring over several years of mission lifetime will also enable detection of wobbles induced by planetary and substellar companions on the primaries. Unlike the radial velocity method, astrometric method produces masses without uncertainties due to orbital inclination. A wealth of information on the distribution of visible and dark mass will come from precise proper-motion measurements of tracer populations.
Neither FAME nor DIVA can transmit entire CCD frames they observe, therefore they require some pre-selection of stars that will be measured. FAME was supposed to have an input catalog with positions of stars brighter than R = 15, while DIVA has an on-board sensor that selects bright stars and transmits parts of the CCD array containing them. Our analysis has led us to conclude that
10 insisting on magnitude-limited samples, as it is planned, would lead to a loss of great scientific opportunities achievable by observing a limited number of fainter objects. These objects would have to be preselected based on scientific merit. In this dissertation (§4.1), I emphasize improvements achievable using this approach in obtaining the local mass function and luminosity calibration of late-M and L dwarfs, and of white dwarfs. Also, by selecting faint horizontal branch stars and measuring their proper motions, valuable information on the motion, and kinematic substructure of the Galactic halo is achievable (§6.1). Furthermore, observing proper motion of many faint quasars, whose motion is known a priori to be zero, enables multi-fold improvement in the definition of the reference frame. We have illustrated these principles using proposed characteristics of FAME. Recently, the funding for this satellite was terminated, however, the same principles can be applied to DIVA, scheduled for launch in 2004.
A more ambitious project, but using the same techniques as Hipparcos and
DIVA, is planned for early in the next decade. GAIA will have precision as high as
12 µas yr−1, and will measure 109 stars to V ∼ 20 (at which point precision falls to
160 µas yr−1). It will also measure radial velocities to V = 17. This will allow the full six-dimensional phase-space to be explored for a large part of the Galaxy. This will make feasible studies that give a full picture of the structure and evolution of the Galaxy. At this point it is too early to make any conclusions whether, as in the
11 case of FAME and DIVA, including still fainter samples of stars (beyond V = 20) is technically sound and scientifically justified.
Still greater accuracies will be achievable with Space Interferometry Mission
(SIM), scheduled for launch in 2011. It should be able to provide 4 µas yr−1 accuracy down to V = 20, regardless of magnitude. Unlike all the previous missions that are constantly scanning the sky and thus produce all-sky surveys, SIM will observe specific targets (some 40,000 during its 5 yr mission lifetime). Also, unlike the other missions, SIM uses optical interferometry to achieve this unprecedented accuracy. Scientific objectives, and therefore the targets to be observed will be defined by the winning proposals. Currently selected proposals (key projects) include a search for planets, obtaining direct distances to globular clusters, measuring proper motions of nearby galaxies, etc. One of the selected projects is astrometric measurement of microlensing events, headed by Andrew Gould.
Classical photometric microlensing is in itself a very powerful technique, especially in the statistical sense, but becomes even more powerful when astrometric effects are observed as well so that distances and masses of lenses can be determined. In this respect, we investigate SIM’s performance and principles required to obtain the bulge mass function of luminous and non-luminous objects (§4.2), and the distances and masses of MACHOs (§6.2). For rare cases of lensing by nearby
12 stars the only observable effect is astrometric and this allows their masses to be measured with high precision (Chapter 3).
Figure 1.1 compares the astrometric accuracies of different missions as a function of magnitude.
1.4. Future of Ground-based Astrometry?
It might seem that with the advent os space-based astrometry there is no place for astrometric studies made from the ground. However, this is not entirely true.
The techniques such as speckle interferometry and adaptive optics in many cases efficiently counter the deleterious effects of the atmospheric seeing. Even better results in achieving high-resolution imaging, and therefore enabling astrometry of high precision, are expected from ground-based optical and infrared interferometers that include telescopes of large apertures. Unlike space missions, these facilities will operate for longer periods of time, thus enabling studies that take more time
(for example, planets with longer periods). Also, they have greater collecting power than satellites and can operate in wavelengths not available to satellites.
Observing in the infrared is necessary in studying the Galactic center, which suffers too much extinction in the optical band. Astrometry of stars lying as close as 1000 AU from the radio source Sgr A* already give strong indications that these
13 stars move in the gravitational potential of a supermassive black hole associated with Sgr A*. In Chapter 5 we propose a method that would combine these astrometric observations with radial velocity measurements to achieve an accurate distance to the Galactic center, a fundamental parameter in Galactic astronomy.
In order for this method to be competitive, observations over a number of years are required, which is another reason why ground based observations are favored.
After decades of neglect, it is clear that once again astrometry, now of very high resolution, plays a crucial role in modern astronomy, and that no field of astronomy or astrophysics can claim to be complete without results provided by and made possible by astrometry.
14 0.1“ FK5
10mas
Hipparcos 1mas
.1mas FAME
Mission Accuracy
10uas GAIA SIM SIM (exceptional targets) 1uas
0.1uas 0 5 10 15 20 V magnitude
Fig. 1.1.— Comparison of astrometric accuracies of different missions as a function of magnitude. 15 Chapter 2
Refinement and Analysis of Luyten’s Catalog of High Proper Motion Stars
2.1. Introduction
More than two decades after its final compilation, the New Luyten Two-Tenths
Catalog (NLTT) of high proper-motion stars (µ>0 .18 yr) (Luyten 1979, 1980;
Luyten & Hughes 1980) and its better known subset, the Luyten Half-Second
Catalog (Luyten 1979, LHS) (µ>0 .5 yr), continue to be a vital source of astrometric data. They are still mined for nearby stars (Reid & Cruz 2002;
Jahreiss, Scholz, Meusinger, & Lehmann 2001; Scholz, Meusinger, & Jahreiss 2001;
Gizis & Reid 1997; Henry, Ianna, Kirkpatrick, & Jahreiss 1997), subdwarfs (Gizis
& Reid 1997; Ryan 1992, 1989), and white dwarfs (Reid, Liebert, & Schmidt
2001; Schmidt et al. 1999; Liebert, Dahn, Gresham, & Strittmatter 1979; Jones
1972; Luyten 1970, 1977). NLTT is at the center of the controversy over whether
16 halo white dwarfs can contribute significantly to the dark matter (Reid, Sahu, &
Hawley 2001; Oppenheimer et al. 2001; Flynn et al. 2001); and it is the primary source of candidates for astrometric microlensing events to be observed by the
Space Interferometry Mission (SIM) (Chapter 3). Despite the advent of many new proper-motion surveys, including Hipparcos (ESA 1997), Tycho-2 (Høg et al.
2000), Starnet (R¨oser 1996), UCAC1 (Zacharias et al. 2000), as well as deeper but more localized surveys: SuperCOSMOS Sky Survey in the south (Hambly et al. 2001; Hambly, Davenhall, Irwin, & MacGillivray 2001), Digital Sky Survey
(DSS)-based survey in the Galactic plane (L´epine, Shara, & Rich 2002), search for µ>0 .4yr stars in ∼ 1400 deg2 (Monet et al. 2000), EROS 2 proper motion search in ∼ 400 deg2 of high latitude fields (EROS Collaboration et al. 1999), and
MACHO search in 50 deg2 towards the bulge and the LMC (Alcock et al. 2001),
NLTT remains unchallenged as a deep, all-sky, proper motion catalog.
NLTT is an all-sky catalog with position and proper-motions (PPM) for stars
−1 above a proper-motion threshold of µlim = 180 mas yr . It extends to V ∼ 19 over
◦ much of the sky, although it is less deep (V ∼< 15) within about 10 of the Galactic plane and also in the celestial south (δ<−30◦). In addition to PPM, NLTT lists somewhat crude photographic photometry in two bands (BNLTT, RNLTT), rough
“spectral types” (usually based on photographic colors), as well as important notes on some individual stars, primarily common proper motion (CPM) binaries.
17 To be sure, the newer catalogs have superseded NLTT in certain domains. By observing a large fraction of the brighter (V ∼< 11) NLTT stars, Hipparcos obtained vastly superior astrometry and photometry for about 13% of NLTT, although in its magnitude limited survey (V<7.3–9.0) it did not find any significant number of new high proper-motion stars not already catalogued by Luyten. Tycho-2, which combined re-reduced Tycho observations of 2.5 million stars with 144 ground-based catalogs (most notably the early 20-century Astrographic Catalog) to derive proper motions, includes PPM and photometry for several thousand additional NLTT stars, and also contains several hundred previously unknown bright (V ∼< 11) high proper-motion stars. However, neither Hipparcos nor Tycho-2 probes anywhere near the faint (V ∼ 19) limit of NLTT. Moreover, their overlap with NLTT has never been systematically studied. UCAC (US Naval Observatory CCD Astrograph
Catalog) is in the process of delivering a new all-sky PPM catalog down to R ∼ 16 based on CCD observations. The first release (UCAC1) covers 80% of the Southern hemisphere. While its photometry is only for a single band (close to R), UCAC1 can easily be matched to 2MASS (Skrutskie et al. 1997) with its JHK photometry, effectively creating a multi-wavelength PPM catalog. Unfortunately, the current release, UCAC1, excludes most of the NLTT high proper-motion stars, due to absence of these stars in Southern parts of USNO-A2 catalog (Monet 1998), which was used as a first epoch. (USNO-A2, and its earlier version USNO-A1,
Monet 1996, were derived by measuring first generation Schmidt plates.) High
18 proper-motion stars pose an especially difficult problem for automated catalog construction: counterparts of slow-moving stars can be reliably identified at different epochs, but fast-moving ones are easily confused with pairs of unmatched, but unrelated stars, or even with spurious objects. These problems grow rapidly worse near the magnitude limit and towards the Galactic plane. Without either new and robust (but difficult to write) algorithms or a vast investment of manual labor, the only routes open in the face of these difficulties are to eliminate the potential high proper-motion stars from the new catalogs or to include them but to acknowledge that many may be spurious. These two approaches have been respectively adopted by UCAC1 and Starnet. (Starnet derives its proper motions for 4.3 million stars by combining the Nesterov, Kislyuk, & Potter 1990 reduction of Astrographic Catalog with GSC 1.0, Lasker et al. 1990.)
There are two other PPM catalogs whose release is promised in the near future, USNO-B (Monet 2000) and GSC-II (Lasker et al. 1998, also known as GSC
2.3 ). Both are based on photographic astrometry and photometry down to or beyond the NLTT limit of V ∼ 19. They use a combination of first and second generation Schmidt plates to determine proper motions. Neither project has stated explicitly how they will handle high proper-motion stars.
Ultimately, space missions such as FAME, DIVA, and/or GAIA will produce reliable catalogs of high proper-motion stars, since each star is observed hundreds
19 of times with great astrometric precision, thus eliminating the possibility of false matching. However, NLTT stars that are fainter than the FAME or DIVA survey limit of R ∼ 15 cannot be observed by them unless they are securely located prior to the mission. Hence, this is yet another reason for obtaining improved astrometry for these stars.
While NLTT has proven incredibly valuable, it also has significant shortcomings. As mentioned above, its two-band photometry is relatively crude, so that classification of stars using the NLTT reduced proper motion (RPM) diagram
(one of the main motivations for constructing the catalog) is rather uncertain: white dwarfs are easily confused for subdwarfs, as are subdwarfs for main-sequence stars. Very few of the NLTT stars have their optical magnitudes in standard bands available in the literature. This problem is not easily resolved by matching with other catalogs that have better photometry, such as USNO-A (photographic BUSNO and RUSNO: σ ∼ 0.3) or 2MASS (JHK: σ ∼ 0.03). This is because although a large fraction of NLTT stars have fairly precise (∼ 10 ) positions, a significant minority have much larger errors, making automated identification of counterparts in other catalogs quite difficult.
Some applications require much higher precision proper motions than NLTT’s characteristic 20–30 mas yr−1. For example, even when faint NLTT stars can be identified in USNO-A (thus establishing their POSS I 1950s positions to ∼ 250 mas)
20 the NLTT proper-motion errors propagate to create ∼ 1 .5 position errors in 2010, too large to reliably predict viable astrometric microlensing events to be observed by SIM (Chapter 3), and so measure precisely the mass of the lens (high proper motion star). These errors can be reduced to ∼ 100 mas, which is quite adequate for our purposes, by matching NLTT stars first to USNO-A and then to 2MASS.
In §2.2 we present the first step in the construction of a refined NLTT catalog: matching bright NLTT entries with those in more recent and higher-precision catalogs, primarily Hipparcos and Tycho-2, but also Starnet. By doing so, we characterize both NLTT’s PPM errors and its completeness at the bright end for the first time. These results are of interest in their own right, but they also serve to guide our search procedure for USNO-A/2MASS counterparts at fainter magnitudes.
That faint-end search is discussed in §2.3. The general approach is be to match (circa 1950) USNO-A stars with their (circa 2000) 2MASS counterparts using NLTT (1950 epoch) positions as a rough guide as to where to find USNO-A counterparts, and NLTT proper motions to predict the positions of their 2MASS counterparts. To correctly apply this approach and to understand its limitations, one must have a good grasp of the NLTT PPM error distributions, and for this it is best to compare NLTT with other PPM catalogs.
21 The bright-end and faint-end searches are complementary. On the one hand, Hipparcos/Tycho-2/Starnet become highly incomplete for V ∼> 12, and so cannot probe faint magnitudes. On the other hand, the POSS plates that were scanned to produce USNO-A saturate for V ∼< 11, leading to increasingly unreliable photometry, astrometry, and even identifications at bright magnitudes.
Nevertheless, these searches do have some overlap at intermediate magnitudes, and we will exploit this overlap to check each method against the other. Hence, we push the bright-end search as faint as we can, incorporating Starnet (with its spurious high-proper motion stars). And we push the faint-end search as bright as we can, eliminating only the Hipparcos/NLTT matches before beginning the search.
By combining the bright-end and faint-end searches we come to our revised catalog (Revised NLTT). The format and the properties of this catalog are described in §2.4.
The Revised NLTT has a wide range of potential applications. Most notably, the resulting V − J RPM diagram (where V is calibrated from BUSNO and RUSNO) permits a much more reliable separation of different classes of stars than did the original BNLTT − RNLTT RPM diagram. The construction and the analysis of RPM diagram is discussed in §2.5. Study of this diagram is useful in its own right, but can also guide the interpretation of RPM diagrams constructed of NLTT stars present in SDSS (York et al. 2000), making use of excellent SDSS photometry and
22 SDSS/USNO-A derived proper motions. A refined version of NLTT with improved astrometry and photometry could also be exploited to study binarity as a function of stellar population and to find new wide-binary companions of NLTT stars.
Finally, in §2.6 we present a thorough study of NLTT’s completeness. NLTT’s completeness (or lack thereof) is central to the debate over the possibility that a significant fraction of the dark matter can be in halo white dwarfs. We first demonstrate that NLTT is in fact complete at the bright end, and explore this completeness as a function of the position in the Galactic coordinate system. Then we argue that if NLTT is complete at the bright end (V ∼< 13), then one could use a bootstrap-type procedure to show that it is also close to complete at the faint end. Finally, we investigate the relative completeness of different types of stars in the Revised NLTT as a function of Galactic latitude.
2.2. Bright Stars in NLTT
2.2.1. Positional Errors of Bright NLTT Stars
To match NLTT stars with counterparts found in other catalogs, we have to start with NLTT’s positions and proper motions. To devise a search strategy, it is therefore crucial that we first evaluate the error distributions of these quantities.
23 Positions in NLTT are nominally given to two different levels of precision, indicated by a flag: in right ascension (RA), either 1 second or 6 seconds of time and in declination (DEC), either 0. 1or1 of arc. The inferior precisions are the rule in the south (δ<−45◦), and the superior precisions are the rule in the north, although for δ>80◦, the RA is given at nominally inferior precision because 6 seconds of time still corresponds to a relatively small arc.
What is the actual PPM error distribution for the ∼ 54, 000 NLTT entries with nominally good positional precision? To investigate this, we first select only
Tycho-2 stars with proper motions µ>180 mas yr−1, i.e., the same limit as of
NLTT. Using Tycho-2 proper motions (having typical errors of just 2 mas yr−1), we propagate their positions back to 1950, the epoch of NLTT. We then find all NLTT entries whose catalog coordinates lie within 3 of each Tycho-2 star – many times larger than the nominal positional accuracy. To ensure a high level of confidence that the counterpart is real, we restrict consideration to stars for which there is one and only one Tycho-2/NLTT match, and for which the unique NLTT match does not have a flag indicating an inferior-precision measurement. Finally, we eliminate
84 matches for which the Tycho-2 proper motion disagrees with the NLTT proper motion by more than 75 mas yr−1. Figures 2.1 and 2.2 show the position differences
(NLTT − Tycho-2) of the resulting 6660 unique matches, on small and large scales respectively. The central portion of this distribution (Fig. 2.1) exhibits an
24 obvious rectangle with dimensions (1 s × 6 ), i.e., exactly the precisions imposed by the rounding truncation of the catalog entries. Hence, one can deduce that for the majority of the entries within this rectangle, Luyten actually measured the positions to much better precision than they were recorded into NLTT. Also obvious in this figure is a halo of stars that have characteristic errors that are of order the discretization noise. We therefore model the distribution of Luyten’s position measurements as being composed of two populations, with intrinsic
Gaussian measurement errors (in arcsec) of σ1 and σ2. Given the discretization of the reported positions, the distribution of residuals in DEC (∆δ) should therefore be,
1 c + w/2 − ∆δ c − w/2 − ∆δ P (∆δ)= q erf −erf 2w σ σ 1 1 c + w/2 − ∆δ c − w/2 − ∆δ +(1 − q) erf −erf , (2.1) σ2 σ2 where c is the center and w is the width of the discretization box in DEC and q gives the relative normalizations of the two populations. We apply equation (2.1) to the N = 5495 stars lying in the strip 7 .1cosδ +0 .5 > ∆α>−7 .9cosδ − 0 .5, and fit only in the region shown in Figure 2.1, i.e., |∆δ| < 12 . (The reason for offsetting the center of this box from 0 will be made clear in Table 2.1.) We find that when w is left as a free parameter, the best fit value is consistent with w =6 , the value expected due to discretization. We then fix w and rederive the other
25 parameters. Next, we write a similar equation for the residuals in RA,
sec δ c + w/2 − ∆α c − w/2 − ∆α P (∆α)= q erf −erf 2w σ sec δ σ sec δ 1 1 c + w/2 − ∆α c − w/2 − ∆α +(1 − q) erf −erf . (2.2) σ2 sec δ σ2 sec δ
We apply equation (2.2) to the N = 5022 stars in the strip 3 .3 > ∆δ>−3 .7, and fit only in the region |∆α| < 12 sec δ. The width is again consistent with the expected value, w =15 , so we again hold w fixed. Table 2.1 shows our results.
Here N is the total number of stars in the subsample and N1 = Nq is the number in the good-precision (σ1) population. For both the α and δ directions, we find