long - lived low - mass stars are a local relic of the earliest stages of the galaxy s formation . dubbed the population ii in contrast to the disk population i , these stars make up the galactic halo ( sometimes called spheroid ) .
although comprising only a very small faction of the stars in the galaxy , population ii stars provide one of the most important records of galactic history .
the population ii luminosity function is of particular interest since its derivation is a requirement in most studies of the stellar mass function .
determination of the stellar mass function as function of galactic epoch and metallicity provides significant constraints upon the theory of star formation .
evolution of the mass function with time or metallicity might be expected since conditions in star - forming clouds were likely much different in the protogalaxy ( c.f .
@xcite ) . moreover , since halo brown dwarfs have faded far below the threshold for even near - infrared detection , the form of the subdwarf mass function near the hydrogen - burning limit represents the most effective technique ( in the absence of gravitational lensing detections ) for estimating the population ii dark matter contribution .    despite the intuitive simplicity of the process of counting stars to determine the luminosity function , there have substantial disagreements between different studies in both the normalization and slope of the luminosity function ( see the review of @xcite ) .
two recent surveys have measured the space density of very - low - mass population ii stars , the m subdwarfs .
counts of faint , distant m subdwarfs using hst ( @xcite ) predict @xmath0 times fewer stars at @xmath1 than are found locally by dahn et al .
dahn et al .
( 1995 ) base their study upon a sample of the local high proper motion stars in the lhs catalog ( @xcite ) for which they have measured accurate trigonometric parallaxes .    in this paper , we determine the luminosity function of the field population ii m subdwarfs using a new proper motion survey based upon plates taken as part of the first and second palomar sky surveys . in section  [ general : halolf ] , we discuss some general problems with the selection of halo stars . in section  [ data :
halolf ] , we describe the data used in our survey and the methods used to select a halo sample . in section  [ lf : halolf ] , we use our data to measure the m subdwarf luminosity function . in section  [ discussion : halolf ] , we compare our luminosity function to other published luminosity functions and discuss its significance . in section  [ conclusion : halolf ] , we summarize our results . in the appendix
, we describe some stars that may be of particular interest .
measurement of the population ii luminosity function is considerable more difficult than the measurement of the population i luminosity function , for the simple reason that the population ii represents a tiny fraction of stars locally .
one must either count very faint stars at large distances from the galactic plane ( @xcite ) , in which case one sacrifices the ability to obtain spectroscopic data , or else determine some way to identify the local population ii stars by excluding the local @xmath2 of stars that belong to the disk .
the latter course requires an efficient method of eliminating the disk stars before time - consuming followup observations are obtained . as members of a high velocity - dispersion , low rotation population , subdwarfs tend to have high heliocentric velocities , and choosing an initial sample of proper motion stars increases the halo contribution from @xmath3 to @xmath4 ( schmidt 1975 ) .
such selection criteria introduce some kinematic bias , but lead to a much more manageable sample .
nevertheless , such a sample is still dominated by old disk ( population i ) and `` thick disk '' ( intermediate population ii , ipii ) stars .    as recognized by schmidt ( 1975 ) and further discussed by bahcall & casertano ( 1986 , hereafter bc ) , even a small proportion of contamination by high velocity disk stars can lead a gross overestimate of the space density of the halo . in order to ensure a pure halo sample for his study , schmidt imposed a restrictive velocity criterion ( @xmath5 km / s ) and then applied a correction to account for halo stars with smaller velocities .
bc used monte carlo simulations to estimate a correction factor of 3.03 for a cutoff velocity of @xmath6 km s@xmath7 ( this correction is dependent on the kinematics assumed for the underlying population ) .
most recently , this technique has been used to derive a halo luminosity function from the lhs catalog stars with @xmath8 and @xmath9 ( @xcite ) .
if trigonometric parallaxes are available for the target stars , then the distances and tangential velocities are estimated easily and the appropriate kinematic criteria can be applied .
this is the case for the schmidt ( 1975 ) sample and the dahn et al .
( 1995 ) samples .
measurement of a useful trigonometric parallax requires years of observing  and the the target must be relatively close ( @xmath10 pc ) .
expansion of the sample to cover a larger space volume requires a less direct method of distance estimation .
photometric parallaxes provide that method , with the calibrating color - absolute magnitude relations provided by subdwarfs with high quality trigonometric parallaxes .
the efficiency of the initial selection of halo candidates can be increased by use of the `` reduced proper motion '' , h , first used by luyten ( 1939 ) .
it is defined as @xmath11 where m is the apparent magnitude and m is the absolute magnitude .
note that @xmath12 is distance independent and determined entirely by observables .
if we assume a color - magnitude relationship , we can plot constant tangential velocity contours on the data .
figure  [ figure - rpm - demo ] illustrates the concept , where data are plotted for 2111 m dwarfs within @xmath1330 parsecs of the sun(reid et al .
1995 ; hawley et al . 1996 ) .
( note that the increasing proper motion bias for redder stars in the gliese & jahreiss 1991 nearby star catalog is clearly evident . )
the lines illustrate where a disk population with halo - like velocities will lie , and how halo stars will be even more easily distinguished in the diagram .    as subluminous stars in the ( m@xmath14 , ( r - i ) ) and ( m@xmath14 , ( v - i ) ) planes , halo stars have systematically larger h than disk dwarfs at a given color , mimicking the effects of increasing @xmath15 . in contrast ,
the coolest metal - poor subdwarfs are superluminous in the ( m@xmath16 , ( b - v ) ) plane ( gizis 1997 ) , and a cutoff in the ( h@xmath16 , ( b - v ) ) diagram tends to discriminate against lower - abundance stars .
we have therefore defined our sample of candidate halo subdwarfs as those stars with reduced proper motion , h@xmath14 , which exceed the value predicted for _ disk _ dwarfs with @xmath17 km s@xmath7 and of the appropriate ( v - i ) or ( r - i ) color . since the halo stars are subluminous , this selection corresponds to a cutoff velocity that is a function of metallicity - an effect which is taken into account in the subsequent analysis .
we utilized the techniques described above to define our sample .
proper motion stars were identified using the first and second epoch 48-inch schmidt photographic plates , as described in ( section  [ mumeasure ] ) .
an initial cut of the sample was made to select large tangential velocity stars ( section  [ rpm ] ) .
we then obtained spectra of each candidate halo star to estimate its metallicity and radial velocity ( section  [ spectroscopy : halolf ] ) .
table  [ table - plates ] lists the fields analyzed in the current survey . in each , the effective area is @xmath18 square degrees , with most having plate material spanning a baseline of at least 40 years . in each case , the first epoch plates were taken with the 48-in .
oschin telescope as part of the first palomar observatory sky survey ( hereafter possi ) .
for the most part , the second epoch plates were taken in the 1980s by either the 48-in .
oschin telescope ( in the course of the palomar observatory second sky survey , possii ) or by the u.k .
schmidt telescope ( ukst ) .
these second epoch plates were used by tinney , reid , & mould ( 1993 , hereafter trm ) in their determination of the disk luminosity function .
all of the possi plates were scanned by the apm facility at the institute of astronomy , cambridge , while the possii / ukst plates were scanned by cosmos at the royal observatory , edinburgh .
one additional field , centered near the north galactic pole , was included . as described by reid ( 1990 ) ,
the second epoch plates for this field were taken at palomar in 1976 , and all plate material was scanned by cosmos .
each of these scans yields positions ( both pixel x , y positions on the plate and right ascension , declination positions on the sky ) , magnitudes , and morphological classification as stars , galaxies , `` merged stars '' or noise .
all objects classified as noise were excluded from further analysis .    as a preliminary to identifying proper motion stars , all objects ( stars and galaxies ) identified in the scans of the plates taken at the same epoch must be matched . in the case of the first epoch
possi data , this is straightforward , since the o and e plates ( corresponding roughly to photometric b and r observations ) were taken on the same night .
thus there is no real motion between the time of the observations and the plate scans can be matched by demanding positional coincidence to a high level of precision . for our second epoch
trm data , the iiiaf ( photometric r ) and iv - n ( photometric i ) plates were not observed on the same night .
trm paired the objects identified on their scans using a 3 arcsec box .
since the trm plates were taken @xmath19 years apart , this procedure loses the highest ( @xmath20 arcsec yr@xmath7 ) proper motion stars . in our survey , this is not a concern because the 40 year separation between the first and second epoch , with no plates at intermediate epochs , prevents our identifying reliably stars with high proper motion due to the large number of possible pairings .
the effective upper limit for our survey is @xmath21=0.375 arcsec yr@xmath7 , and that limit is taken into account in our luminosity function analysis .
we should note that the trm detections are restricted to the region within 3 degrees of the plate center in order to avoid the vignetted regions near the edge of the plates , and that the cosmos scans did not cover the entire plate .
the appropriate area of the survey is given in table  [ table - plates ] .
moreover , the the possi field centers do not correspond to the possii / ukst field centers .
thus , between one and four possi plates are paired with the trm fields , as listed in table  [ table - plates ] .
the identification of the proper motion stars was made as follows .
first , an initial pairing of objects within a 3 arcsecond box ( @xmath22 ) was made . in each field
, the surface density of celestial objects is sufficiently low that matching within this radius results in an unambiguous pairing .
figure  [ figure - close513 ] illustrates the agreement between the @xmath23 and e magnitudes for one representative pairing ( possii field 513 and possi e102 ) .
most of the outliers evident in this plot are correct pairings of extended objects ( galaxies or merged stars )  they are outliers due to differences in the definitions of magnitudes used by the two machines , as well as the different plate characteristics . for this study , we are interested only in stars , so such outliers are unimportant .
having paired the unambiguous matches , we obtain a list of unmatched objects  these objects have no counterparts within 3 arcseconds at the second epoch and are good candidates for proper motion stars . eliminating paired objects
, we find all possible matches within 15 arcsec of unmatched objects .
we exclude pairings which have magnitudes inconsistent with the e @xmath23 relation delineated by the close matches  inconsistent being defined as a discrepancy of from @xmath241 magnitude @xmath25 to @xmath26 magnitudes at @xmath27 .
these limits are `` conservative '' in the sense that some false pairings are included in order to avoid excluding true proper motion stars
.    the vast majority of objects on the plate should show insignificant proper motion between the two epochs , but directly matching the right ascension and declination catalog positions leads to many ( false ) proper motion objects .
those spurious motions arise from systematic errors in the relative positions of the two catalogs ( caused by differing plate centers , inadequacy of the plate model , plate variations , telescope differences , etc . )
these systematic positional errors can be corrected by determining the transformation between the original pixel positions for the two epochs ( @xmath28 and @xmath29 ) . both our own experience with the trm fields and the reid ( 1990 ) analysis showed that a polynomial solution over the entire plate works poorly , producing many false proper motion stars .
analysis over smaller regions , however , is effective .
each plate was divided into approximately one by one degree regions and the positional transformation determined using second - order two - dimensional polynomials .
only stars with @xmath30 matches were used for the initial co - ordinate transformation , iterating twice to exclude objects with large ( @xmath31 ) residuals .
we then used these polynomial solutions to determine the proper motions of _ _ a__ll stars in the field . the stars that prove to have @xmath32 arcsec yr@xmath7
are included in our sample of proper motion stars .
the selection of the candidate halo stars is described in section  [ rpm ] .
these proper motions are measured relative to faint stars , which have typical distances of @xmath33 pc .
thus we expect the reference frame to have a small , but non - zero , mean proper motion .
we derive the conversion from this relative frame to an absolute ( extragalactic ) frame using measurements of objects classified as galaxies .
the latter are identified using the @xmath34-parameter defined by picard ( 1991 ) .
we take @xmath35 for these purposes , which provides a sample dominated by galaxies .
the upper cutoff is taken to exclude the `` fuzziest '' galaxies  we found that these `` fuzzy '' galaxies have large spurious proper motions , which we attribute to a combination of the different plate material used at the different epochs , and the different scanning and centroiding techniques used at the apm and cosmos . in figure
[ figure - absmu ] , we plot the distribution of galaxy proper motions and stellar proper motions for the fields .
it is clear that the galaxies do indeed show a systematic shift with respect to the stars .
table  [ table - absmu ] lists the median proper motions ( @xmath36 ) derived for the galaxies , which at most amount to @xmath37 mas  yr@xmath7 .
we have checked these values by comparison with the projection of the solar motion and rotational lag for the ipii and halo populations .
the direction of motion derived by the galaxies ( i.e. , the _ sign _ of the proper motions ) agrees with the expectations .
these values given in table  [ table - absmu ] are used to correct the relative proper motions to an absolute system : @xmath38 we also give the observed standard deviation of the galaxian motions , with values in the range @xmath39 mas  yr@xmath7 . these provide an estimate of the uncertainties of our proper motions , but we expect that the extended galaxies should have poorer positions , and hence less accurate proper motions , than is the case for stars .
it is evident in table  [ table - absmu ] that field 831 shows the largest random errors and has absolute proper motions correction that are discrepant at the @xmath40 mas level from the adjoining fields .
r. mendez ( private communication ) has kindly computed model proper motions for this direction based upon the galactic structure model of mendez & van altena ( 1996 ) .
the predicted mean stellar proper motions of ( @xmath41 ) = ( @xmath42 ) mas  yr@xmath7 are in excellent agreement with our observations in fields 829 and 832 .
we take this as confirmation that our simple correction technique using galaxies is adequate for our purposes .
the data for field 831 are evidently of lower astrometric quality , but the accuracy remains adequate for identifying the high proper motion stars analyzed in this study .
the nltt catalog ( @xcite ) offers an opportunity to check both the completeness and accuracy of our measurements . as noted by many authors ,
identifying the nltt stars can be difficult since some of the positions are quite poor ( for example , we find a star that matches the magnitudes and proper motion of lp 621 - 64 , but is 10 _ arcminutes _ east of luyten s position ) . excluding stars which are too bright , too faint , or too blue ( i.e. , too faint to be seen on the possii iv - n plates )
, we find that 122 of 137 nltt stars are recovered , indicating a completeness of at least @xmath43 .
we find that the standard deviations of @xmath44 and @xmath45 are in the range 10 to 25 mas yr@xmath7 .
luyten estimated his motions to be accurate to @xmath46 mas yr@xmath7 .
the combination of our uncertainties estimated from the galaxies and the nltt uncertainties account for the observed scatter .
selection of halo candidates using reduced proper motion requires photometry in addition to the proper motions .
we have used the photometry from the second epoch plates which provides us with r and i measurements for each star ( except in the ngp field , for which we have v and i ) .
we adopt the photometric calibration determined by trm using ccd measurements for each plate . because their ( `` palomar '' ) calibration is on a unique system
, we describe it briefly here .
a complete description is given by trm .
the best defined photometric system in r and i for late - type stars is the cousins system ( @xcite ) is cousins r - i , @xmath47 is trm s palomar system , and @xmath48 is gunn r - i . ] however , m dwarfs have strong color terms in this system because the r@xmath49 band filter extends far into the i band  as a result , the effective wavelength of the r@xmath49 filter moves redward for cooler stars , as discussed by bessell ( 1986 ) .
the possii iiiaf emulsion plus filter combination has a sharper long - wavelength cutoff at @xmath50 , so the color term differs significantly from the standard cousins system .
the possii response is close to the gunn r filter defined by thuan & gunn ( 1976 ) .
therefore , trm s ccd calibrations of the fields used gunn r@xmath51 and i@xmath51 filters but observed _ cousins _ late - type standards .
these ccd magnitudes were matched to the cosmos photographic magnitudes .
they found that no significant i color term was necessary to calibrate the observations but , as expected , a strong r color term was necessary .
thus , their photometric observations consist of an r magnitude , denoted @xmath23 , an i magnitude , denoted @xmath52 , and a color , denoted @xmath47 , for each star .
the relation between the standard cousins color and the palomar color is @xmath53 @xmath54 note that trm use the relation @xmath55 , but examination of their figure  2 indicates that there is no offset between @xmath52 and @xmath56 for @xmath57 , which is the color of the halo stars in our survey .
the v and i photometric calibration for the ngp field is described in reid ( 1990 ) .
the reduced proper motion diagram for field 513 is shown in figure  [ figure - rpm513 ] .
the solid line marks the disk main sequence defined by trm , @xmath58 offset to match @xmath59 km s@xmath7 .
stars falling below this line ( larger h@xmath14 ) are identified as halo candidates for spectroscopic followup .
recently , we have pointed out that the disk main sequence shows evidence for a `` kink '' at spectral type m4.0 v ( @xcite ; @xcite ) . while this feature is not included in the trm calibration
, it occurs at @xmath60 , so it is not relevant to the current selection of candidate halo subdwarfs .
while our tangential velocity limit is set at @xmath6 km s@xmath7 for solar abundance stars , m subdwarfs are selected at smaller @xmath15 .
the esdm enter the candidate pool for @xmath61 km s@xmath7 and the brighter sdm at @xmath62 km s@xmath7 .
the most likely contaminants , modestly metal - poor ( @xmath63 \approx -0.6 $ ] ) intermediate population ii m dwarfs , lie @xmath64 magnitudes below the disk main sequence ( @xcite , figure 7 ) , and therefore enter at @xmath65 km s@xmath7 .
we obtained optical spectra of the candidate halo stars in order to determine radial velocities and metallicities .
the hale 200-in . and the las campanas du pont 100-in .
telescopes were used to obtain data in the wavelength range @xmath66 with @xmath67 resolution .
these parameters were chosen to correspond to the observations of nearby m dwarfs ( @xcite ; @xcite ) and m subdwarfs ( @xcite ) that we have already published .
for the 200-in .
observations , we used the double spectrograph ( oke and gunn 1982 ) . in august 1995 ,
the blue camera was set to observe @xmath68 and the red camera was set to @xmath69 using @xmath70 l / mm gratings blazed at @xmath71 and @xmath72 respectively . in october 1995 , a new red camera was installed in the double spectrograph and was used in all subsequent runs to observe the region @xmath73 at @xmath74 with the @xmath70 l / mm grating blazed at @xmath72 . with this setup , the blue camera was used to cover @xmath75 with a @xmath76 l / mm grating blazed at @xmath77 . in practice ,
the latter data were useful only for the bluest ( @xmath78 ) stars .
the du pont observations used the modular spectrograph with a 1200 line grating blazed at @xmath79 .    from these spectra , we determined radial velocities and spectral types . during each observing run , we observed k and m radial velocity standard stars drawn from the list of marcy & benitz ( 1989 ) .
a velocity for each candidate halo star was determined by cross - correlating with the best match standard star .
comparison of subdwarfs with known radial velocities observed during the same runs shows that the accuracy is @xmath80 km s@xmath7 ( @xcite ) . using these radial velocities , we measured tio and cah bandstrength indices , as defined originally by reid et al .
each index measures the ratio of the flux ( @xmath81 ) within the absorption feature to nearby pseudo - continuum points .
table  [ table - index ] lists the wavelength definitions of the features ( w ) and the pseudo - continuum ( s1,s2 ) points .
we have defined a spectral classification system using these indices ( @xcite ) .
stars are classified as disk stars ( m v ) , m subdwarfs ( sdm ) , or extreme m subdwarfs ( esdm ) .
we have shown on the basis of model atmosphere calculations ( @xcite ) that these empirical classifications correspond to @xmath63 \approx -0.5 $ ] , @xmath63 \approx -1.2\pm 0.3 $ ] , and @xmath63 \approx -2.0 \pm 0.5 $ ] ( where @xmath63\approx[fe / h]$ ] ; @xcite ) . each halo candidate has been classified using this system .
in this section , we calculate the i - band luminosity function from the data described in section  [ data : halolf ] .
the luminosity function , @xmath82 , is defined as the number of stars per cubic parsec per magnitude in a bin centered at @xmath83 .
we use schmidt s ( 1968 ) @xmath84 technique to estimate @xmath82 @xmath85 @xmath86 where @xmath87 is the solid angle on the sky of the survey , @xmath88 is the maximum distance that a given star could have been detected in this survey , and @xmath89 is the corresponding volume . a minimum distance , @xmath90 , also appears to account for the distance at which the star would have evaded detection due to the upper proper motion limit of @xmath91 yr@xmath7 .
the maximum distance may be limited by either the lower proper motion limit ( @xmath92 yr@xmath7 ) or the limiting magnitudes ( @xmath93 ) , hence for a star , distance d , magnitudes ( r , i ) , @xmath94 ; \rm{dex}[0.2(i_0-i)]\right)\ ] ] we estimate the associated error by assuming that each star contributes an uncertainty of @xmath84 ( @xcite ) .
then @xmath95    as noted in section  [ general : halolf ] , we apply a tangential velocity cutoff in order to ensure a `` pure '' halo sample .
we adopt a value of @xmath96 km s@xmath7 in our main analysis . as a result
, a correction must be applied to the apparent luminosity function in order to obtain the space of density of halo stars .
thus , @xmath97 the discovery fraction @xmath98 is a function of position on the sky and the adopted @xmath99 . in our case , we compute @xmath98 for each field in the survey .
this procedure is described in section  [ kinematic : halolf ] .
a check on the completeness of the survey can be made by considering @xmath100 which considers the ratio of the volume @xmath101 corresponding to the star s actual distance @xmath102 to @xmath89 . in the case of a spatially uniform sample , such as the halo over the volume probed by this survey
, we expect @xmath103 .
we find that eight of our nine fields are within @xmath104 of 0.5 for sdm / esdm with @xmath105 km s@xmath7 , and conclude that , as a whole , the survey is complete statistically .
we use the photometry and parallax data compiled by gizis ( 1997 ) to determine color - magnitude relations appropriate for each spectral class of star ( m v ; sdm ; esdm ) .
there are many sdm and esdm with accurate distance determinations available , mainly due to the efforts of the usno ccd parallax program ( @xcite ) , which allows us to determine the following relations : @xmath106 @xmath107 these relations are valid for @xmath108 .
our survey is aimed at the fainter , redder m subdwarfs , but for comparisons to other surveys of higher mass stars it is useful to have an estimate of the main sequence for the bluer subdwarfs . for these stars , the above linear relations are no longer appropriate , due to an inflection in the main sequence at @xmath109 ( c.f .
, @xcite ; @xcite )
. unfortunately , there are few parallax k subdwarfs and fewer with red ( ri ) photometry .
the hipparcos mission , however , has determined parallaxes for many g subdwarfs .
we therefore fit a linear relation to the hipparcos g subdwarfs supplemented by the available k subdwarfs from the fourth edition of the yale parallax catalog ( van altena et al .
1996 ) . in this region
, models predict that for @xmath63 \lesssim -1 $ ] the absolute magnitudes are relatively insensitive to metallicity ( @xcite ) .
we select only stars with @xmath110 km s@xmath7 to ensure a halo sample .
the data are shown are in figure  [ figure - hipp ] .
there are 29 subdwarfs with @xmath111 , @xmath112 , and @xmath113 .
we find @xmath114 this slope is consistent with the slope measured in globular clusters by santiago et al .
( 1996 ) .
the above relations can only be used for the ngp field , where we have v and i photometry .
unfortunately most sdm and esdm lack r photometry , while we have only r and i data for the remaining fields in our survey .
hence we have determined the relationship between @xmath115 and @xmath47 using equation  [ ri - ri ] and bessell s ( 1990 ) vri@xmath49 photometry of nearby stars .
we find : @xmath116 applying this transformation to the @xmath115 color of each parallax star , we find that @xmath117 @xmath118 relations  [ ms - ri - sdm ] and  [ ms - ri - esdm ] are valid for @xmath119 .
it is possible that the r - i , v - i transformation for the extreme m subdwarfs is not the same as that for disk m dwarfs . while the currently available photometry ( summarized in @xcite ) is unable to resolve this issue , we have obtained spectrophotometry ( @xcite ) which indicates that there is no significant color term for the sdm stars or the earliest ( esdk7 ; esdm0.0 ) extreme subdwarfs .    for the g subdwarfs , we transform equation  [ ms - vi - sdg ] and obtain : @xmath120 equation  [ ms - ri - sdg ] applies for @xmath121 .
the luminosity function computed using the techniques described above , although correct in shape , represents only the fraction of the halo that has @xmath122 km  s@xmath7 . given a kinematic model of the halo , we can calculate the discovery fraction ( @xmath98 ) for each field using monte carlo simulations .
bc have shown that variations of @xmath123 km  s@xmath7 in the kinematic parameters ( @xmath124,@xmath125 , @xmath126,@xmath127 ) imply variations of less than @xmath128 in @xmath98 ; sadly , published data for halo stars span a much larger range of kinematics .
we therefore compare the observed kinematics of our survey stars to the predictions of representative halo models .
analyses of local ( d@xmath129 1 kpc ) halo stars generally derive mild prograde rotation ( e.g. , @xmath130 km s@xmath7 , norris 1986 ) , although recent studies based on more distant ( @xmath245 kpc ) halo stars have found evidence for a retrograde rotating halo ( see the review of majewski 1993 ) .
we compare our data to the models of bahcall & casertano ( 1986 , bc ) , norris ( 1986 , n ) , casertano et al .
( 1990 , crb ) , layden et al .
( 1996 , l ) , and beers & sommer - larson ( 1995 , bsl ) .
the last is of considerable interest as it features a retrograde rotating halo .    for each of these models
, we generated a monte carlo simulation of our survey with @xmath131 detected stars .
the halo velocities are represented by gaussian distributions in u , v , and w.,@xmath132 ) , v represents motion in the direction of galactic rotation , and w represents motion perpendicular to the galactic plane . ] in table  [ table - compareuvw ] , we compare the predictions of the u , v , and w velocity distributions of the models to the actual data for the 35 metal - poor stars with @xmath133 .
we emphasize that these are the predicted velocity characteristics of the kinematically - biased observed sample , not the true halo kinematics .
we also give the probability that the observed data and the model simulations are drawn from the same distributions , as determined by kolmogorov - smirnov tests .
the comparison is most sensitive to the observed v velocities , which are most important in determining @xmath98 .
it clear that all the published models agree fairly well with the observed u and w velocity distribution , although the n model is significantly poorer for u , due to its smaller velocity dispersion ( @xmath134 km s@xmath7 ) .
the v velocities favor prograde halos , with the bc model ( @xmath135 km s@xmath7 ) yielding too small a mean v velocity and the n model ( @xmath136 km s@xmath137 ) yielding too large a mean velocity .
the bsl model s retrograde halo is much less likely than the prograde models .
of the published models , the layden et al .
model gives the highest probability when all three velocity components are considered .
guided by the comparison above , we have constructed a composite kinematic model , taking the parameters derived by norris as a starting point and increasing both @xmath125 to match other studies , and the halo rotation to a value midway between bc and n. the resulting parameters are @xmath138 in km s@xmath7 and we refer to this as our model a. these changes are within @xmath139 of norris s estimated uncertainties .
the values of @xmath140 for the bc and n models lie within @xmath141 of the composite values , which are listed in table  [ table - corr ] .
values of @xmath140 based on the layden et al . halo model
are on average @xmath142 smaller .
we have also calculated models including the ipii ( bsl , crb ) . only one in @xmath143 ipii stars will enter our @xmath110 km s@xmath7 sample ( although we caution that this number depends upon the assumption of gaussian kinematics ) .
since the ipii consists of approximately @xmath144 of stars locally , it will contribute less than @xmath145 of our observed halo sample and thus will not bias significantly our measurements .
this can also been seen in table 1 of bc , where it is shown that the kinematic properties of stars in eggen s ( 1979a , 1980 ) proper motion catalog with @xmath146 , @xmath147 and @xmath148 km s@xmath7 are indistinguishable .
we expect that the ipii contribution is much less than this , since direct spectroscopic observation shows that most ipii m dwarfs both locally ( nearby stars with @xmath149 km s@xmath7 , gizis 1997 ) and at height of a few kiloparsecs ( @xcite ) are not classified as sdm or esdm .
we reject all spectroscopic m  v stars from our sample , and therefore we expect negligible ipii contamination .    for each field ,
we calculate the luminosity function using the equations  [ equation - firstlf ] to  [ equation - lastlf ] using only stars with @xmath150 km s@xmath7 .
we combine data for all fields to derive our best estimate luminosity function , calculated for 0.5 magnitude bins , which is listed in table  [ table - lf ] .
these conservative selection criteria ensure a pure population ii halo sample but suffers from the uncertainty of the kinematic correction .
since we have obtained spectroscopy for all of our reduced proper motion selected stars , we have information on the actual space density of metal - poor very - low - mass stars with lower velocities ( @xmath151 km s@xmath7 for esdm ; @xmath62 km s@xmath7 for sdm ) .
we have therefore computed the space density of sdm and esdm separately , using the appropriate velocity cutoff for each metallicity class .
the completeness corrections , @xmath152 , for these velocity cutoffs are only @xmath153 , since virtually all halo stars have @xmath62 km / s .
the resulting luminosity functions for the sdm and esdm are listed in table  [ table - lf - lowvtan ] .
this luminosity function is remarkably consistent with that estimated from the corrected `` pure '' halo sample .
the normalizations differ by no more that @xmath142 from that predicted by our composite model .
we discuss these low-@xmath15 samples in more detail in section  [ discussion : halolf ] .
our measurement of the population ii luminosity function is aimed at only the limited range in luminosity corresponding to the fainter m subdwarfs .
nonetheless , these results offer a means of resolving the difference , noted in the section  [ halolf - intro ] , between gould et al.s hst starcount analysis at @xmath154 kpc and dahn et al.s local trigonometric parallax sample . before discussing that issue
, we consider whether our results can be combined with previous measurements of the bright end of the halo luminosity function .
since previous studies of the halo luminosity function have presented @xmath155 for all metallicities combined , we have transformed our data using the relation @xmath156 and the color - magnitude relations given in section  [ section : colormag ] .
the results for the combined sdm and esdm luminosity function is plotted in figure  [ figure - field - mv ] .
the other luminosity functions shown are the local trigonometric parallax samples of schmidt ( 1975 ) and dahn et al . , the local photometric parallax sample of bahcall & casertano ( 1986 ) , and the gould et al .
hst starcount analysis .
we have not plotted dawson s ( 1986 ) statistical analysis of the lhs stars , since it is superceded by dahn et al.s calculations , which have the advantage of improved data for individual stars .
schmidt s results have been updated using literature v - band photometry , but there are only 0 - 2 stars per bin , except for the last two bins which have five stars apiece .
it should also be noted that , except for the hst data , these estimates of @xmath157(m@xmath16 ) are based on kinematically - selected samples .
dahn et al s @xmath15 correction factors , @xmath98 , agree to within 5% with the values predicted by our model a , and we have also used that model to adjust schmidt s data .
the bahcall & casertano @xmath157(m@xmath16 ) plotted is derived using their published correction factors .
adopting @xmath98 appropriate to model a reduces the inferred space densities by @xmath158 , increasing the discrepancy with respect to the other studies ( including our own results ) .
it is obvious from figure  [ figure - field - mv ] that there are potential problems involved in combining data from the several investigations to derive a halo luminosity function spanning the full range in absolute magnitude . of the surveys which include brighter stars near the turnoff , schmidt s
luminosity function ( as noted originally by schmidt ! )
has too few stars to measure reliably the shape of the luminosity function .
the bc dataset has sufficient stars , but there is a clear , systematic offset of almost a factor of two between their derived space densities and dahn et al s results .
bahcall & casertano s analysis is based on eggen s ( 1979a , 1980 ) survey of southern ( @xmath159 ) stars with @xmath160 yr@xmath7 and v@xmath161 , but a comparison of both the relative numbers and @xmath162 distributions of such stars shows no evidence for any significant increase in incompleteness at southern declinations . on the other hand ,
the color - magnitude relation ( from eggen , 1979b ) which bc use to estimate distances is a poor description of the lower pop .
ii main sequence ( m@xmath163 ) , and therefore that portion of @xmath157(m@xmath16 ) should probably be excluded .
the calibration is likely to be more reliable for the more luminous stars , since the main sequence is consistent with the hipparcos data in figure  [ figure - hipp ] for the g subdwarfs .
one notable aspect of the eggen / bc sample is that it may be incomplete in the galactic plane , since it is difficult to identify proper motion stars against the dense stellar background .
we have simulated halo samples using the bc , a , and l models for both `` all - sky '' ( @xmath159 ) and `` non - plane '' ( @xmath164 and @xmath165 ) selection . the kinematic correction factors ( @xmath140 for @xmath110 km s@xmath7 ) for the `` non - plane '' samples are greater than the `` all - sky '' samples by between 8% and 12% depending on the model .
however , our analysis of the eggen ( 1979a , 1980 ) data shows that the derived halo space density for stars with @xmath166 is actually @xmath167 higher than that derived for the entire eggen sample .
this suggests that the bc luminosity function is underestimated by @xmath4 due to incompleteness in the galactic plane .
this would account for about a third of the discrepancy evident in figure  [ figure - field - mv ] .    naively combining the higher - luminosity stars from bc with dahn et al . and
our own results for fainter stars suggests a rather steep luminosity function .
globular clusters provide a potential means of testing this result .
figure  [ figure - gci ] compares @xmath157(m@xmath14 ) from the field against results for two of the better - studied clusters , using data from piotto et al .
( 1997 ) adjusted to be consistent with the new distances derived by reid ( 1997 ) on the basis of main - sequence fitting .
that comparison suggests that the bc results may underestimate the space densities of stars with @xmath168 , consistent with the offset between the bc and dahn et al results plotted in figure  [ figure - field - mv ] .
it should of course be noted that the field and cluster luminosity functions may differ  the present - day luminosity ( mass ) function of globular clusters could be affected by dynamical evolution , or alternatively the initial mass function of field and cluster stars could be different .
extension of the dahn et al .
sample to brighter absolute magnitudes ( e.g. , expanding the schmidt sample ) should provide the most effective means of confirming the shape of the field luminosity function .
the dahn et al . luminosity function , transformed via equation  [ ms - vi - esdm ] is also plotted , but the transformation is uncertain since we have assumed all stars to be esdm , which is not the case  and in any case , this preliminary dahn et al .
luminosity function represents a range of metallicities , so it may not be fair to compare it to the globular clusters .
nevertheless it is clear that the dahn et al .
luminosity function peaks at a fainter absolute magnitude than the globular clusters .
as noted above , figure  [ figure - field - mv ] suggests that the gould et al ( 1998 ) hst analysis underestimates the local density of late - type halo subdwarfs by a factor of four . to underscore that point , figure
[ figure - field - mi ] plots our luminosity function against a @xmath82 representation of the hst data that was computed by gould ( private comm . ) using the same data and techniques used to derive the @xmath169 luminosity function of gould et al .
while our study is limited to a restricted range of luminosity , the space densities , derived from a completely independent sample of nearby subdwarfs , are clearly strongly supportive of the lhs - star analysis .
a possible resolution of the discrepancy between the distant - halo and local - halo analyses may lie in the composite nature of the galactic halo .
sommer - larsen & zhen ( 1990 ) estimate that at least 40 percent of the local subdwarfs reside in a highly - flattened component , which would essentially be absent from the hst dataset . taken in isolation ,
this additional component can account only partially for the factor of four offset in the derived densities .
however , its presence underlines the substantial extrapolation involved in transforming the observed number density at large distances from the plane to a local halo subdwarf density
. an additional contribution could be due to the choice of local calibrating subdwarfs for the gould et al .
analysis , which may be biased towards higher metallicities .
this would lead to overestimated m@xmath16 and distance , and hence an underestimate of the space density . for a discussion of possible systematic errors , see gould et al.s discussion .
the kinematics of the very - low - mass halo stars are consistent with those found for higher mass metal - poor stars .
the strongest conclusion is that prograde rotating halos are favored .
if the expanded halo sample that includes lower tangential velocity stars is included , then kinematics like those of model a is favored over layden et al.s rr lyrae kinematics , since the expanded luminosity function is more consistent with the model a corrected pure halo sample ( tables  [ table - lf ] and  [ table - lf - lowvtan ] ) . there
_ may _ be an excess of sdm at low velocities over the halo model a predictions . since metallicities of @xmath63 \approx -1.2 \pm 0.3 $ ]
are expected to have some ipii contribution , this is not too surprising . if the two best determined bins are used , then the excess is only @xmath170 .
chiba & yoshii ( 1998 ) find an excess of @xmath171 in the range @xmath172 \le -1.0 $ ] based upon hipparcos data for red giant and rr lyrae stars , in excellent agreement with our results .
the mass function can be estimated using the observed luminosity function and a mass - luminosity relation as follows : @xmath173 this is often parameterized as a power - law with the form @xmath174 since there are no direct determinations of mass for any sdm or esdm stars , we must rely upon theoretical mass - luminosity relations .
these must be viewed with caution  for example , reid & gizis ( 1997 ) have shown that , in the case of the baraffe et al solar - abundance models , the slope deduced from nearby - star data for the disk mass function depends on the passband used in the analysis .
our data do not provide sufficient range in @xmath83 to estimate the mass function reliably , even given a good mass - luminosity relationship , but it does allow us to estimate the relative numbers of sdm and esdm stars at a mass of @xmath175 .
we transform the observed pure halo @xmath83 luminosity functions into mass functions using the baraffe et al .
( 1997 ) models .
we adopt the @xmath63= -1.0 $ ] and @xmath63= -1.5 $ ] interior models for the sdm and esdm respectively , and fit a power - law to the mass function in order to obtain the best estimate of the numbers .
the ratio of sdm to esdm at @xmath176 is 0.38 . in comparison ,
we obtain a ratio of 0.36 if we adopt the alexander et al .
( 1997 ) models , and 0.60 for the dantona & mazzitelli ( 1996 ) tracks  but the latter match the observed @xmath177 main sequence poorly .
these results are consistent with the carney et al .
( 1994 ) distribution of metallicities for higher mass stars ( mainly g subdwarfs )  @xmath178 of their stars in a pure ( @xmath179 km s@xmath7 ) halo sample have @xmath63>-1.5 $ ] .
thus , the metallicity distribution for the very - low - mass stars of the halo ( @xmath180 ) is apparently the same as that of stars on the upper main sequence ( @xmath181 ) . barring a cosmic conspiracy of offsetting parameters ,
the most reasonable explanation is that the mass function in the halo is independent of metallicity .
although our data is inadequate to constrain the mass function well , we note that the best - fit parameters are @xmath182 for the esdm and @xmath183 for the sdm .
these compare with a power - law index of @xmath184 for disk stars in the mass range @xmath185 ( @xcite ) .
the consistency of the metallicity distribution and kinematics of the low - mass stars ( sdm and esdm ) is not surprising , but nevertheless important .
studies of the halo , whether using evolved stars such as rr lyrae or red giants , or main sequence stars ( f , g subdwarfs ) , assume that the properties of stars in the narrow mass range from just below the turnoff to the tip of the asymptotic giant branch ( @xmath186 , are representative of the halo population as a whole .
our results verify one aspect of this assumption for stars near the bottom of the halo main sequence .
we have derived a new estimate of the space density of very low mass metal - poor stars .
comparisons with other surveys shows that we are in agreement with the dahn et al .
( 1995 ) survey based upon trigonometric parallaxes of lhs catalog stars .
this implies that the local space density of metal - poor stars is higher than that predicted by hst observations of more distant stars .
our work provides measurements of the kinematics and metallicity distribution of very - low - mass stars .
these measurements are consistent with those from studies of higher - mass stars in the halo .
we find that the esdm are 2.5 times as common as sdm at @xmath187 .
this is consistent with the relative numbers of metal - poor ( @xmath63 \le -1.5 $ ] ) and metal - rich ( @xmath188 > -1.5 $ ] ) g subdwarfs found by carney et al . ( 1994 ) . note that the two metallicity scales are comparable because gizis & reid ( 1997 ) have found agreement between the carney et al .
metallicities of fgk subdwarf primaries with metallicities based upon g97 .
the agreement of the metallicity distribution at both low and high mass implies that the mass function is not a strong function of metallicity for ( @xmath189 ) halo stars .
this result is supported by the general similarity of the mass functions derived for the two metallicity bins considered ; however , those mass functions depend upon uncertain model mass - luminosity transformations .
we thank chris tinney and mike irwin for assistance with the plate scans and the staff of palomar observatory for their capable support .
we also thank m. schmidt for useful conversations and a. gould for providing the i - band hst luminosity function .
j.e.g . is grateful for partial support through greenstein and kingsley fellowships as well as nasa grants
go-05353.01 - 93a , go-05913.01 - 94a , and go-06344.01 - 95a , this work is based partly on photographic plates obtained at the palomar observatory 48-inch oschin telescope for the second palomar observatory sky survey which was funded by the eastman kodak company , the national geographic society , the samuel oschin foundation , the alfred sloan foundation , the national science foundation grants ast84 - 08225 , ast87 - 19465 ,
ast90 - 23115 and ast93 - 18984 , and the national aeronautics and space administration grants ngl 05002140 and nagw 1710 .
this research has made use of the simbad database , operated at cds , strasbourg , france .
most of our subdwarfs lie at distances of greater than 200 parsecs , and therefore are more difficult to study than lhs catalog subdwarfs .
three of our stars , however , are close enough that additional observations may be profitable .
all were previously identified by luyten ( 1979 - 1980 ) . the positions ( b1950 , at the epoch given in table  [ table - plates ] ) and other data for these stars
are given in table  [ table - stars ] .
spectra of the subdwarfs lp 622 - 7 and lp 382 - 40 are plotted in figure  [ figure - stars1 ] .
lp 622 - 7 has a proper motion too large to be included in our complete sample , and was observed at the palomar 60-in . in june 1995 . the radial velocity may be unreliable due to a problem with the calibrating arc lamp exposure .
lp 589 - 7 is one of the coolest subdwarfs known .
its spectral indices are tio5@xmath190 , cah1=@xmath191 , cah2=@xmath192 , and cah3@xmath193 . in figure
[ figure - stars2 ] , we show that this star is extremely similar to the esdm5.0 ( @xcite ) star lhs 3061 ( @xmath194 , @xcite ) .
the conspicuous tio absorption at @xmath195 distinguishes both these stars and the slightly cooler lhs 1742a ( @xmath196 , esdm5.5 ) from the other two very cool extreme subdwarfs lhs 205a ( @xmath194 , esdm5.0 ) and lhs 1826 ( @xmath83 unknown , esdm6.0 ; @xcite ) .
lp 589 - 7 is just as close as lhs 1742a , so it is an excellent candidate for a trigonometric parallax .
there is a hint of @xmath197 emission in our spectrum of lp 589 - 7 .
the feature , which may not real but lies at the correct wavelength , has an equivalent width of @xmath198 .
two sdm systems ( lhs 482 and lhs 2497 ) are known to show emission due to close companions ( @xcite ) .
since lp 589 - 7 is near the hydrogen burning limit , any companion must be either a very faint , old white dwarf ( as for the lhs 482 system ) , a nearly equal luminosity m companion ( as for lhs 2497 ) , or else a brown dwarf ( if it proves to be too faint to contribute detectable light ) .
alternatively , lp 589 - 7 could a `` young '' metal - poor star with an age of up to a few billion years .
if an sb2 , its estimated distance and tangential velocity should be increased by up to @xmath199 .
all this is speculative , however , since the @xmath197 feature is weak and may not be real .
alexander , d.r . ,
brocato , e. , cassisi , s. , castellani , v. , ciacio , f. , & deglinnocenti , s. 1997 , , 317 , 90 allard , f. , & hauschildt , p.h .
1995 , , 445 , 433 bahcall , j.n . , & casertano , s. 1986 , , 308 , 347 ( bc ) baraffe , i. , chabrier , g. , allard , f. , & hauschildt , p.h .
1997 , , 327 , 1054 baraffe , i. , chabrier , g. , allard , f. , & hauschildt , p.h .
1998 , , 337 , 403 beers , t.c .
, & sommer - larsen , j. 1995 , , 96 , 175 bessell , m.s .
1986 , , 98 , 1303 carney , b.w . ,
latham , d.w . ,
laird , j.b . , and aguilar , l.a .
1994 , , 107 , 2240 casertano , s. , ratnatunga , k.u . , bahcall , j.n .
1990 , , 357 , 435 chiba , m. , & yoshii , y. 1998 , , 115 , 168 cool , a.m. , piotto , g. , & king , i.r .
1996 , , 468 , , 655 cousins , a.w.j . 1973 ,
mem.r.a.s . , 77 , 223 dahn , c.c . , liebert , j. , harris , h.c . , &
guetter , h.h .
1995 , the bottom of the main sequence  and beyond , edited by c.g .
tinney ( berlin , springer ) , p.239 dantona , f. , & mazzitelli , i. 1996 , , 456 , 329 dawson , p.c .
1986 , , 311 , 984 eggen o.j .
1979a , , 39 , 89 eggen o.j .
1979b , , 230 , 786 eggen o.j .
1980 , , 43 , 457 felton , j.e . 1976 , ,
207 , 700 gizis , j.e .
1997 , , 113 , 806 gizis , j.e .
1998a , ph.d .
thesis , california institute of technology gizis , j.e .
1998b , , 115 , 2053 gizis , j.e . , &
reid , i.n .
1996 , , 111 , 365 gizis , j.e . , & reid , i.n .
1997 , , 109 , 894 gliese , w. , & jahreiss , h. 1991 , preliminary version of the third catalog of nearby stars ( cns3 ) gould , a. , flynn , c. , bahcall , j.h .
1998 , , in press hawley , s.l . , gizis , j.e . , & reid , i.n .
1996 , , 112 , 2799 layden , a.c . ,
hanson , r.b . , hawley , s.l . ,
klemola , a.r . , & hanley , c.j .
1996 , , 112 , 2110 luyten , w.j .
1939 , publ .
minnesota , 2 , 121 luyten , w.j .
1979 , catalogue of stars with proper motions exceeding 0.5 annually ( lhs ) , ( university of minnesota , minneapolis , minnesota ) luyten w. j. , 1979 - 1980 , new luyten catalogue of stars with proper motions larger than two tenths of an arcsecond , ( minneapolis , university of minnesota ) , computer - readable version on adc selected astronomical catalogs vol.1 - cd - rom marcy , g.w . , & benitz , k.j .
1989 , , 344 , 441 majewski , s.r .
1993 , , 31 , 575 mendez , r.a .
, & van altena , w.f .
1996 , , 112 , 655 monet , d.g . , dahn , c.c .
, vrba , f.j . ,
harris , h.c . ,
pier , j.r . ,
luginbuhl , c.b . , &
ables , h.d .
1992 , , 103 , 638 mould , j.r .
1996 , , 108 , 35 murray , s.d .
, & lin , d.c .
1996 , , 467 , 728 norris , j. 1986 , , 61 , 667 oke , j.b . , & gunn , j.e . 1982 , , 94 , 586 paresce , f. , de marchi , g. , & romaniello , m. 1995 , , 440 , 216 picard , a. 1991 , , 102 , 445 piotto , g. , cool , a.m. , king , i.r .
1997 , , 113 , 1345 reid , n. 1990 , , 247 , 70 reid , i.n .
1997 , , 114 , 161 reid , i.n . , hawley , s.l . , & gizis , j.e .
1995 , , 110 , 1838 reid , i.n . , & gizis , j.e .
1997 , , 113 , 2246 reid , i.n . , gizis , j.e . , cohen , j.g . ,
pahre , m.a . ,
hogg , d.w . , cowie , l. , hu , e. , & songaila , a. 1997 , , 109 , 539 santiago , b.x . , gilmore , g. , elson , r.a.w .
1996 , , 281 , 871 schmidt , m. 1968 , , 151 , 393 schmidt , m. 1975 , , 202 , 22 sommer - larsen , j. & zhen , c. 1990 , , 242 , 10 stobie , r.s . ,
ishida , k. , & peacock , j.a .
1989 , , 238 , 709 tinney , c.g . ,
reid , i.n .
, & mould , j.r .
1993 , , 414 , 254 ( trm ) thuan , t.x . , & gunn , j.e .
1976 , , 88 , 543 van altena , w.f . , lee , j.t . , & hoffleit , e.d .
1995 , the general catalogue of trigonometric stellae parallaxes , fourth edition ( schenectady , l. davis press ) zinnecker , h. 1995 , the bottom of the main sequence  and beyond , edited by c.g .
tinney ( berlin , springer ) , p.257    rrlccclc [ table - plates ] 829 & 02 & 00 & + 00 & uk / pii & 1987.83 & eo852 & 1953.78 + 831 & 02 & 40 & + 00 & uk / pii & 1981.92 & eo1283 & 1954.90 + & & & & & & eo1453 & 1955.81 + 832 & 03 & 00 & + 00 & uk / pii & 1986.89 & eo363 & 1951.69 + & & & & & & eo1453 & 1955.81 + 262 & 10 & 15 & + 45 & possii & 1987.34 & eo672 & 1953.12 + ngp & 13 & 04 & + 29 & palomar & 1976.23 & eo1393 & 1955.29 + 868 & 15 & 00 & + 00 & uk / pii & 1987.30 & eo1402 & 1955.30 + & & & & & & eo1613 & 1957.32 + 513 & 15 & 00 & + 25 & possii & 1987.34 & eo102 & 1950.35
+ & & & & & & eo87 & 1950.30 + & & & & & & eo1390 & 1955.29 + & & & & & & eo1092 & 1954.49 + 889 & 22 & 00 & + 00 & possii & 1990.55 & eo1146 & 1954.57 + 890 & 22 & 20 & + 00 & possii & 1990.82 & eo364 & 1951.69 + & & & & & & eo1146 & 1954.57 +    crrrr [ table - absmu ] & mas & mas & mas & mas + 262 & 2.6 & 5.2 & 13.0 & 13.7 + 513 & 4.4 & 4.2 & 12.1 & 11.8 + 829 & -3.7 & 3.8 & 14.9 & 13.9 + 831 & -2.0 & 2.3 & 17.0 & 16.0 + 832 & -4.3 & 3.9 & 13.6 & 12.6 + 868 & 3.2 & 3.3 & 13.6 & 12.6 + 889 & -0.2 & 4.3 & 10.6 & 10.6 + 890 & -0.9 & 3.4 & 11.9 & 11.5 + ngp & 5.0 & 3.9 & & +      cccc [ table - corr ] 829 & 25.32 & 17.3 & 2.07 + 831 & 26.49 & 17.3 & 2.15 + 832 & 25.32 & 17.3 & 2.11 + 262 & 14.56 & 17.3 & 2.30 + ngp & 28.00 & 17.0 & 2.01 + 868 & 25.29 & 17.3 & 2.11 + 513 & 25.29 & 17.3 & 2.18 + 889 & 25.30 & 17.3 & 3.70 + 890 & 25.26 & 17.3 & 3.56 +    cccccccc [ table - compareuvw ] @xmath125 & 171 & 171 & 146 & 179 & 165 & 155 & 158 + probability ( u ) & & 0.43&0.15&0.55&0.29&0.19&0.31 + @xmath124 & -234 & -222&-244&-240&-252&-264&-233 + @xmath126 & 76 & 82 & 82 & 84 & 73 & 76 & 85 + probability ( v ) & & 0.37&0.49&0.75&0.22&0.05&0.87 + @xmath127 & 90 & 77 & 86 & 97 & 94 & 105 & 86 + probability ( w ) & & 0.38&0.49&0.32&0.38&0.28&0.49 +    rrlrlr [ table - lf ] & & 8.5 & 2.2 & 0.7 & 1.8 & 0.5 & 16 9.0 & 2.0 & 0.6 & 1.6 & 0.5 & 11 9.5 & 2.0 & 0.9 & 1.6 & 0.8 & 5 10.0 & 2.0 & 2.0 & 1.6 & 1.6 & 1 9.5 & 4.6 & 1.8 & 3.9 & 1.5 & 10 10.0 & 3.6 & 1.6 & 2.9 & 1.3 & 5 10.5 & 3.1 & 2.3 & 2.6 & 2.0 & 2    rrlrlr [ table - lf - lowvtan ] & & 8.5 & 2.6 & 0.8 & 2.5 & 0.7 & 20 9.0 & 2.4 & 0.8 & 2.2 & 0.7 & 14 9.5 & 4.4 & 1.3 & 4.1 & 1.2 & 12 10.0 & 3.1 & 1.6 & 2.9 & 1.5 & 4 10.5 & 3.6 & 3.6 & 3.4 & 3.4 & 1 9.5 & 6.9 & 1.8 & 6.6 & 1.7 & 20 10.0 & 2.5 & 1.0 & 2.4 & 1.0 & 6 10.5 & 4.1 & 2.1 & 3.9 & 2.0 & 4 11.0 & & & & & 0 11.5 & & & & & 0 12.0 & 15 . & 15 . & 15 . & 15 . &
1    lrrrrrrrrrrr [ table - stars ] lp 589 - 7 & 28.72333 & + 1.03647 & -0.12 & -0.33 & 15.90 & 1.38 & 11.77 & 64 & 106 & -65 & esdm5 + lp 382 - 40 & 223.53375 & + 23.83433 & -0.28 & -0.09 & 16.47 & 1.08 & 10.59 & 150 & 210 & -230 & esdm2.5/3 + lp 622 - 7 & 227.18958 & -1.02281 & + 0.20 & -0.43 & 13.87 & 0.87 & 9.69 & 69 & 155 & 300 & esdk7/m0 +