A&A 421, 763-770 (2004)
DOI: 10.1051/0004-6361:20047009

The Liverpool-Edinburgh high proper motion catalogue[*]

R. S. Pokorny1 - H. R. A. Jones1 - N. C. Hambly2 - D. J. Pinfield1


1 - Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead CH41 1LD, UK
2 - Wide Field Astronomy Unit, Institute for Astronomy, School of Physics, University of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK

Received 5 January 2004 / Accepted 31 March 2004

Abstract
We present a machine selected catalogue of 11 289 objects with proper motions exceeding 0.18 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$ and an R-band faint magnitude limit of 19.5 mag. The catalogue was produced using SuperCOSMOS digitized R-Band ESO and UK Schmidt Plates in 287 Schmidt fields covering almost 7000 square degrees ($\sim$17$\%$ of the whole sky) at the South Galactic Cap. The catalogue includes UK Schmidt BJ and I magnitudes for all of the stars as well as 2MASS magnitudes for 10 447 of the catalogue stars. We also show that the NLTT is $\sim$95$\%$ complete for Dec > -32.5$^{\circ }$.

Key words: stars: kinematics - stars: statistics - stars: subdwarfs - Galaxy: solar neighbourhood - stars: late-type - stars: low-mass, brown dwarfs

1 Introduction

A number of recent studies have shown the power of proper motion surveys, particularly when combined with infra-red photometry, to identify various disk and halo populations, accurately characterize the nearby stellar census, identify nearby brown dwarfs, find common proper motion binaries and pick out members of kinematic moving groups. We briefly summarise below some recent studies in these areas.

Reid & Cruz (2002), Reid et al. (2002), Cruz & Reid (2002) and Reid et al. (2003) used NLTT data cross-referenced with the 2MASS Point Source Catalogue to search for nearby stars, both for follow up observation of specific stellar types, and for determination of the luminosity function, the mass-luminosity relationship, the stellar contribution to the local mass density, the velocity distribution and the stellar multiplicity statistics.

Rojo & Ruiz (2003) used 112 stars taken from three areas of the Calan-ESO survey (Ruiz et al. 2001), identified with medium resolution spectroscopy, determined their Galactic velocities and constructed luminosity functions for the disk and halo populations using the 1/ $V_{{\rm max}}$ method.

Reyle et al. (2002) took the sample of 693 stars from Scholz et al. (2000) and identified 301 of these stars in the DENIS catalogue. They used the I-J colour to estimate the distances of these stars and found that at least 107 of them were within 50 pc, including 15 within 25 pc and one star, LHS124, which could be at 7.2 pc. Lepine et al. (2002) report the discovery of an M 8.5 star with proper motion $\mu $ = 2.38 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$ at a distance of 13.9 $\pm$ 3.5 pc. Lepine et al. (2003) also reports the discovery of over 30 more stars within 25 pc. These studies enhance the local census of stars and enable more accurate determinations of the luminosity functions of both the disk and halo populations.

Salim & Gould (2002, 2003; and Gould & Salim 2003) have obtained improved astrometry and photometry for stars in the NLTT with Dec > -17$^{\circ }$ by matching that catalogue with the Hipparcos, Tycho-2 & Starnet catalogues for bright stars, and the USNO-A and 2MASS catalogues for fainter stars. This revision to the NLTT has allowed them to improve the characterization and statistics for the NLTT.

Silvestri et al. (2002) took a sample of 13 white dwarf - M dwarf high proper motion binaries with what they describe as "halo-like'' velocities. 12 of these 13 M dwarfs have near-solar metallicity, from which they conclude that these systems are part of the kinematic tail of the thick disk population. They further state that this shows that the white dwarf components are not, therefore, part of the Galactic dark matter halo.

Ribas (2003) has suggested that the high proper motion field brown dwarf LP 944-20 is probably a member of the Castor moving group. Assuming this to be the case allows an age and a metallicity to be applied, which means that the properties of LP 944-20 may be compared to models.

An object which is present in our catalogue is a third component to the M dwarf white dwarf binary LHS4039/LHS4040 which Scholz et al. (2004) recently identified as an active M 8.5 dwarf.

In this paper we expand on the previously published first release of the Liverpool-Edinburgh High Proper Motion Catalogue (Pokorny et al. 2003, hereafter PJH03). The aim of the project was to produce an internally homogeneous large area survey with well defined completeness limits. We chose to use only the two R-band datasets in order to prevent losing faint objects with extreme colours which would not show up in a search using the BJ and I band plates. PJH03 discussed the importance of proper motion surveys and described the survey data, search technique and the statistical tests conducted on the survey. We now present the fully completed version of the Liverpool-Edinburgh High Proper Motion Catalogue, update our analysis of the completeness levels for the survey and present 2MASS photometry for the catalogue sources.


  \begin{figure}
\par\includegraphics[angle=180,width=12.2cm,clip]{0009fg1.eps}\end{figure} Figure 1: Plot of the fields used in this survey. Fields shaded in dark were used, fields shaded in light were available, but could not be used either because of a short epoch difference between the ESO and UK plates or because they contained too few bright stars to calculate the full relative astrometry. Unshaded fields below Dec = $-20^{\circ }$ could not be used because they are too close to the Galactic plane or cover the Magellanic clouds.
Open with DEXTER

Section 2 summarises the proper motion measurements and the scanned photographic plate magnitudes. Section 3 details how we cross matched our proper motion sources with the 2MASS database. Section 4 discusses ongoing and planned future work with the catalogue.

2 The catalogue

2.1 Proper motion data

We have now completed the survey and catalogue, using a further 156 R-band Schmidt plates digitised by the SuperCOSMOS machine at the Wide Field Astronomy Unit of the Royal Observatory Edinburgh (Hambly et al. 1998, 2001a,b,c). The complete catalogue contains 11 899 objects (approximately twice that of the first catalogue release) with proper motions exceeding 0.18 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$ and R-band magnitude $\leq $ 19.5 mag in both plates, of which 610 are secondary or subsequent detections giving 11 289 separate objects. Multiple detections are left out of all but one of the plots in this paper. The full survey used 287 Schmidt fields with plate centres at and below -20$^{\circ }$ declination covering almost 7000 square degrees ($\sim$17$\%$ of the whole sky). This comprises all of the available sky out of the Galactic plane south of -20$^{\circ }$. Figure 1 shows the distribution of the fields used for the survey. Figure 2 shows the distribution of epoch differences for the R-band plates that were used for the survey, whilst Fig. 3 shows the distribution of high proper motion objects on an equal area plot of the sky. No fields with epoch differences of less than 3 years were used. The mean epoch difference is $\sim$8.5. Almost all of the stars have BJ and I magnitudes, and 2MASS J, H and  $K_{{\rm S}}$ magnitudes (see Sect. 3). In the magnitude columns in Table 1 a value of -1.00 indicates that no plate was scanned in for that area and a value of 0.00 indicates that no object match was found in that passband.

In addition, the second release has updated $\sigma _\mu $ values, where we note that in the previous release we used an incorrect formula that generally overestimated these values. It should also be noted that the $R_{{\rm ESO}}$, $R_{{\rm UK}}$ and I band magnitudes will differ slightly from those given by the SuperCOSMOS Sky Survey (SSS) (http://www-wfau.roe.ac.uk/sss/index.html). This is because the SSS magnitudes are updated as improved calibrations become available.

2.2 Completeness

As discussed in PJH03 a plot of log(Nc) against log($\mu $) should have a slope of -3 for a complete sample of a population with constant space density. Figure 4 is such a plot for all non-multiple detections in the catalogue. The least squares fit for this data has a slope of -3.057. This gives us confidence that our search procedures, as described in PJH03, are robust.

Figure 5 is a histogram of the fractional proper motion errors ( $\sigma _\mu $/$\mu $) for all objects in the catalogue. It is truncated at 0.33 because we removed all objects with fractional errors greater than this value in order to minimize any false detections, although the binning means that it appears to truncate at 0.34. We also visually inspected all objects with proper motions $\mu $ $\ge$ 1.0 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$. The simulated star completeness test detailed in PJH03 applies generally to the Schmidt plate data and the search process used for the whole of the survey. The results of that test are therefore applicable to the results for this second section of the survey. A version of the completeness test was conducted on a sample strip of fields running across the Galactic plane, to assess the minimum Galactic latitude it would be feasible to conduct the survey down to. The results indicate that although the completeness started to drop at $\mid $b$\mid $ $\sim$ 35$^{\circ }$, the search process would still run with reasonable efficiency to $\mid $b$\mid $ $\sim$ 20$^{\circ }$. The survey was therefore conducted on all available fields down to this limit. Table 2 shows how the relative completeness changes with low Galactic latitude, assuming that the number of proper motion stars is invariant with Galactic latitude, and the completeness for $\mid $b$\mid $ > 35$^{\circ }$ is 95$\%$. As discussed in PJH03 the NLTT (Luyten 1979b) and LHS (Luyten 1979a) catalogues can be split into two areas; the northern area comprising all fields with Dec $\geq $ -30$^{\circ }$ which are covered by the Palomar Observatory Sky Survey (POSS, Cabanela et al. 2003), and the southern area which comprises all fields with Dec $\leq $ -35$^{\circ }$ which are covered by the Bruce Proper Motion Survey (BPM, Williams 2000). The drop in completeness shown in Table 2 can be seen in the plot comparing the number of stars in our catalogue to the number of stars in the NLTT for the northern POSS area of the NLTT (Fig. 6). The section of our catalogue which corresponds to the BPM survey area is still significantly more complete than the BPM (Fig. 7). For comparison Fig. 8 shows the number of stars in this catalogue and the NLTT for -32.5$^{\circ }$ $\leq $ Dec $\leq $ -17.5$^{\circ }$ and $\mid $b$\mid $ $\geq $ 40$^{\circ }$. This shows that for R $\ga$ 13 this catalogue (for $\mid $b$\mid $ $\geq $ 40$^{\circ }$) and the NLTT (northern POSS section) contain almost exactly the same numbers of stars, and can therefore be assumed to have the same completeness in these areas. For this catalogue the completeness has been shown to be $\sim$95$\%$ for $\mid $b$\mid $ $\geq $ 35$^{\circ }$, which suggests that the northern section of the NLTT (Dec $\geq $ -32.5$^{\circ }$) is also $\sim$95$\%$ complete. Gould & Salim (2003) use comparison with the Hipparcos and Tycho-2 catalogues to estimate the completenes of the NLTT as a function of Galactic latitude and find that it is $\geq $95$\%$ complete for $\mid $b$\mid $ $\geq $ 20$^{\circ }$ and V $\la$ 11.5, which is in good agreement with our findings for fainter NLTT stars at mid and high Galactic latitudes.


  \begin{figure}
\par\includegraphics[angle=270,width=7.7cm,clip]{0009fg2.eps}\end{figure} Figure 2: Histogram of the epoch differences between the ESO and UK R-band plates for all of the fields used in the survey.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=7cm,clip]{0009fg3.eps}\end{figure} Figure 3: Equal area plot of all objects in the catalogue.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=7.7cm,clip]{0009fg4.eps}\end{figure} Figure 4: Cumulative number of stars in this survey as a function of their proper motions. The offset straight line has a slope of -3. Error bars are $\surd$N $_{{\rm C}}$.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=7.9cm,clip]{0009fg5.eps}\end{figure} Figure 5: Histogram of the fractional error on proper motion ( $\sigma _\mu $/$\mu $) for all entries in the catalogue.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=7.7cm,clip]{0009fg6.eps}\end{figure} Figure 6: Number of high proper motion stars per 1.5 magnitudes in R for the Northern POSS areas of the NLTT: Solid line - this catalogue, dashed line - NLTT.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=7.7cm,clip]{0009fg7.eps}\end{figure} Figure 7: Number of high proper motion stars per 1.5 magnitudes in R for the Southern BPM areas of the NLTT: Solid line - this catalogue, dashed line - NLTT.
Open with DEXTER

3 Cross matching with 2MASS

We have used the All-sky-data-release of the Two Micron All Sky Survey (2MASS) Point Source Catalogue (Cutri et al. 2003), to provide near infrared photometry for our catalogue sources. 2MASS provides broadband J, H, and $K_{{\rm S}}$ photometry for point-like and extended sources over the full celestial sphere. The primary distinction between 2MASS and many other photometric systems is that the Ks-filter transmission truncates at long wavelengths to avoid terrestrial H2O absorption (see Persson et al. 1998). The effective wavelength of the $K_{{\rm S}}$ filter is 2.15 microns (cf. 2.19 microns for standard K filters). A secondary difference is in the J-band filter. The short and long wavelength cutoffs of this filter extend into the atmospheric H2O absorption features at 1.1 and 1.4 microns, which means that 2MASS J magnitudes can be affected by atmospheric water vapour content (by up to $\sim$0.1 mag). Carpenter (2001) derives transforms between 2MASS magnitudes and those in other photometric systems, and finds only small differences between 2MASS $K_{{\rm S}}$ magnitudes and standard K magnitudes. For a source to appear in the 2MASS catalogue, it must either have a signal-to-noise ratio (SNR) > 5 ( $\sigma _m <$ 0.22) in each band, or a SNR > 7 ( $\sigma _m <$ 0.16) in at least one band with fainter detections in the other(s). The 2MASS survey completeness limits are J=15.8, H=15.1 and $K_{{\rm S}}=$ 14.3. The SNR = 5 limits are J=17.2, H=16.2 and $K_{{\rm S}}=$ 15.8. And the commonly quoted SNR = 10 ( $\sigma _m <$ 0.11) limits are J=16.1, H=15.2 and  $K_{{\rm S}}=$ 14.9.

Table 1: Sample table from the full online catalogue, columns are; catalogue ID, RA, Dec (J2000), multiple detection flag, ESO/SRC field number, ESO R-band plate epoch, Epoch difference of observations, total $\mu $ (arcsec per year), error on $\mu $ (arcsec per year), position angle, magnitudes - BJ, $R_{{\rm ESO}}$, $R_{{\rm UK}}$, I, J, H, $K_{{\rm S}}$. For clarity four columns have been omitted, these are the components of the proper motion in RA and Dec and their errors. In the magnitude columns -1.00 indicates that no plate was scanned in for that area and 0.00 indicates that no match was found in that passband.

Table 2: Relative average completeness at different Galactic latitudes assuming 95$\%$ completeness for $\mid $b$\mid $ $\geq $ 35$^{\circ }$.


  \begin{figure}
\par\includegraphics[angle=270,width=7.7cm,clip]{0009fg8.eps}\end{figure} Figure 8: Number of high proper motion stars per 1.5 magnitudes in R for -32.5$^{\circ }$ $\leq $ Dec $\leq $ -17.5$^{\circ }$ and $\mid $b$\mid $ $\geq $ 40$^{\circ }$: Solid line - this catalogue, dashed line - NLTT.
Open with DEXTER

In order to allow a positional match with the 2MASS catalogue, we therefore applied proper-motion corrections to our catalogue coordinates to transform them to epoch 1999.5 (An approximate epoch for 2MASS (north or south) is 1999.5 $\pm$ 2). In order to search for 2MASS counterparts, we defined an error ellipse for each source that would cover its likely location on the 2MASS sky. The axes of these error ellipses were calculated as follows. We began with a base level representing the intrinsic positional uncertainty in the coordinates derived for both our survey and 2MASS. The main source of uncertainty in this base level comes from the possibility that some sources could be partly resolved multiple systems. Although 2MASS will initially flag a multiple system (e.g. a binary) with a $\sim$1 arcsec separation as an extended source (using shape analysis; see 2MASS web-site http://www.ipac.caltech.edu/2mass), it may then subsequently remove it from the extended source catalogue (and put it in the point source catalogue) using additional shape analysis (e.g. symmetry metrics). Therefore, the 2MASS PSC will have some multiple sources with wider angular separation. The worst seeing accepted for 2MASS results in a FWHM of $\sim$4 arcsec, and so it is possible that a binary with an angular separation of this order could appear in the PSC. The relative colours and brightnesses of the multiple components in such a system could result in optical and NIR source positions that differ by up to this level. This positional uncertainty dwarfs the astrometric uncertainties (up to 0.2 and 0.26 arcsec for our survey and 2MASS respectively), and we consequently set our base level uncertainty at this 4 arcsec value. In addition to this base level, there are the uncertainties associated with our proper motion measurements. We duly multiplied these proper motion uncertainties (in RA and Dec) by the difference between the epoch of their observation and 1999.5 $\pm$ 2 (to cover the full 2MASS base-line), and combined the results with the base level value to give final error ellipse axis in both coordinates for each source. We then searched for 2MASS sources within the error ellipses of our transformed catalogue coordinates using the GATOR tool provided by the Infrared Science Archive (IRSA). We added the stipulations that a 2MASS source must be at least 6 arcsec away from its nearest neighbour (to avoid spatial confusion from line-of-site proximity), and that it has not been flagged as a contamination/confusion affected source. Where multiple matches were found, the 2MASS counterpart nearest to the centre of the ellipse was selected.


  \begin{figure}
\par\includegraphics[angle=90,width=7.6cm,clip]{0009fg9.eps}\end{figure} Figure 9: Plot of $K_{{\rm S}}$ (2MASS) against $R_{{\rm ESO}}-K_{{\rm S}}$ for all stars with $K_{{\rm S}}$ magnitudes. Vertical solid lines indicate the $R_{{\rm ESO}}-K_{{\rm S}}$ limits for the stellar types indicated. The horizontal dashed line indicates the SNR = 5 limit for the 2MASS  $K_{{\rm S}}$ magnitude, the diagonal dashed line indicates the R= 19.5 mag cutoff of our catalogue. Large points are stars with $\mu $ $\geq $ 1.0 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$. The square indicates the position of the benchmark brown dwarf LP 944-20 (Tinney 1998).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=14.1cm,clip]{0009fg10.eps}\end{figure} Figure 10: Plot of reduced proper motion, HB, against $B-K_{{\rm S}}$ for all stars in the survey with B and $K_{{\rm S}}$ magnitudes. The lines showing the different regions of the plot are for guidance only.
Open with DEXTER

In order to check for possible mis-matches amongst the highest proper motion sources, we carried out the following test. We selected the 38 highest proper motion sources (>1 arcsec/yr) from our catalogue, and cross-matched them with 2MASS as before, but this time with the added stipulation that a 2MASS source must be redder than $J-K_{{\rm S}}=$ 1.0 to provide a match. This tested the possibility that some high proper motion red sources (most likely nearby late M/L dwarfs) were being missed by chance alignment with bluer background sources. However, no additional 2MASS matches were found that had not already been identified without the colour stipulation. Our 2MASS matching procedure obtained 2MASS photometry for 10 870 of the 11 899 sources in our catalogue, or $\sim$92$\%$.

Figure 9 is the $K_{{\rm S}}$ versus $R_{{\rm UK}}-K_{{\rm S}}$ colour-magnitude diagram for all stars in the catalogue with 2MASS $K_{{\rm S}}$ magnitudes, and demonstrates the ability of this catalogue to identify late M and early L class brown dwarfs. Figure 10 is the HBJ versus $B_{J}-K_{{\rm S}}$ reduced proper motion diagram for the catalogue (HBJ= B + 5 + 5log $\mu $ which can be shown to be the same as HBJ= BJ + 5log  $V_{\rm T}-3.379$ where $V_{\rm T}$ is the tangential velocity of the star; PJH03).

4 Discussion

The main aims of high proper motion surveys are to preferentially select nearby stars and Halo stars and to differentiate between the Halo and Disk populations. As shown by Salim & Gould (2002) and PJH03, optical-infrared reduced proper motion diagrams give much better differentiation of the disk, halo and white dwarf populations than optical-optical reduced proper motion diagrams. However, as passbands further into the infrared are used white dwarfs are lost due to their bluer colours (see Fig. 10), and the optimum reduced proper motion - passband combination for white dwarf selection seems to be HBJ versus B-I. Selection of objects near the Hydrogen burning limit (i.e. late M dwarfs and brown dwarfs) is best done with an optical-infrared colour magnitude diagram (see Fig. 9). Reliable differentiation of the various populations of stars is of great importance to studies of Galactic structure and dynamics, and the Halo luminosity function.

In PJH03 we used the HBJ, B-I and colour-colour plots to differentiate the disk and halo populations. Compared to Table 6 of Reid (1984) our numbers fell between those of models A and D and excluded models B and C. Although we have more than doubled the number of fields in the survey since that analysis the added fields were almost all at intermediate and low Galactic latitudes, and the predictions from Reid are for polar proper motion surveys. This means that we have added very few fields at high Galactic latitude, that the relative numbers of halo and disk stars at high latitude has not changed and that the new fields have therefore not made any difference to that result.

We have recently produced an extended version of the catalogue with a lower proper motion limit of $\mu $ $\ge$ 0.1 $\hbox {$^{\prime \prime }$ }$  ${\rm yr}^{-1}$ and an R-band faint limit of R $\le$ 21 for the specific purpose of identifying proper motion binaries and potential moving group members. This extended version of the catalogue will not be made generally available because of the high relative uncertainties on the proper motions and the difficulty in characterizing the completeness that these limits produce. The fainter magnitude limit of the extended catalogue means that we will be able to use it to identify potential field brown dwarfs which would be too faint to appear in the basic catalogue. This is well demonstrated by the recovery of Kelu-1 in the extended catalogue with an R-band magnitude of  $R_{\rm {UK}}=$ 19.75.

Acknowledgements
The authors would like to thank all of the staff of the SuperCOSMOS project at the ROE and Mike Read for producing Fig. 1. R.S.P. acknowledges support of a Liverpool John Moores University studentship. D.J.P. acknowledges support of a PPARC fellowship. The authors also would like to thank the referee, Dr. Ralf Scholz, for his extremely helpful comments on this paper.

References



Copyright ESO 2004