A&A 448, 579-588 (2006)
DOI: 10.1051/0004-6361:20054113
D. Carollo1 - B. Bucciarelli1 - S. T. Hodgkin3 - M. G. Lattanzi1 - B. McLean2 - R. Morbidelli1 - R. L. Smart1 - A. Spagna1 - L. Terranegra4
1 - INAF - Osservatorio Astronomico di Torino,
10025 Pino Torinese, Italy
2 -
Space Telescope Science Institute, 3700 San Martin Drive
Baltimore, MD 21218, USA
3 -
Cambridge Astronomical Survey Unit, Institute of
Astronomy, Madingley Road,Cambridge, CB3 0HA, UK
4 -
INAF - Osservatorio Astronomico di Capodimonte, via
Moiariello 16, 80131 Napoli, Italy
Received 27 August 2005 / Accepted 2 November 2005
Abstract
The GSC-II white dwarf survey was designed
to identify faint and high proper motion objects, which we used to
define a new and independent sample of cool white dwarfs.
With this survey we aim to derive new constraints on the halo white dwarf space density.
Also, these data can provide information on the age of
thick disk and halo through the analysis of the luminosity
function. On the basis of astrometric and photometric
parameters,
we selected candidates with
and
in an area of 1150 square degrees.
Then, we separated white dwarfs
from late type dwarfs and subdwarfs by means of the reduced proper
motion diagram. Finally, spectroscopic
follow-up observations were carried out to confirm the
white dwarf nature of the selected candidates.
We found 41 white dwarfs of which 24 are new discoveries.
Here we present the full sample and for each object
provide positions, absolute proper motions, photometry, and
spectroscopy.
Key words: stars: white dwarfs - astrometry - techniques: spectroscopic - surveys
Understanding the nature of the Baryonic Dark Matter in the
galactic halo is a very fundamental subject of astrophysics. The
most favoured candidates are brown dwarfs, planets, ancient cool
white dwarfs (WDs), neutron stars and primordial black holes. All
these objects are known as MACHOs (Massive Compact Halo Objects)
and, in principle, could be detected by their gravitational
microlensing effect on field stars (Pacynski 1986).
The two major experiments for the detection of microlensing events
are MACHO (Alcock et al. 2000) and EROS (Afonso et al. 2003). They measured several events toward the
Galactic bulge and the Large (LMC) and Small Magellanic Clouds
(SMC). The excess of events towards the Magellanic Clouds is
directly related to the dark matter in the halo of our Galaxy, and
their interpretation is still under discussion. One potential
explanation for these events is that the lenses belong to the halo
and, in this case, the MACHO experiment seems to indicate that
of the dark halo is made of compact objects with mass
of
(Alcock et al. 2000). At the moment
this evaluation of the baryonic content of the dark halo is
controversial because, recently, the EROS-2 collaboration claimed
a new value of
,
using a refined set of data (Tisserand
& Milsztajn 2005). Moreover, new
photometric measures of the microlensing event MACHO-LMC-5
confirmed that the lens is an M 5 dwarf star of
and
could suggest another explanation for the microlensing events,
i.e. the lenses belong to a previously undetected component of the
disk of the Milky Way (Nguyen et al. 2004). If the
lenses belong to the halo, the natural candidates that explain the
MACHOs and EROS results are very old cool WDs with a mean mass
around
,
and several surveys have been carried out
in order to detect a significant number of Population II WDs (see
the review by Hansen & Liebert 2003 and Reid
2005).
Here we present a new survey based on the material used for the construction of the GSC-II (Guide Star Catalogue II) (McLean et al. 2000). This survey was developed in order to search for faint and high proper motion objects and to define a new sample of cool WDs. This paper presents the catalogue of the spectroscopically confirmed WDs. Section 2 describes the GSC-II material, astrometry, photometry, and proper motions. In Sect. 3 we present the selection criteria for the candidates, while in Sect. 4 we give details of the spectroscopic observations. Finally in Sect. 5 we describe the WD sample and discuss some interesting objects.
The Guide Star Catalog II is an all-sky database based on the
photographic surveys carried out with the Schmidt telescopes of
the Palomar and Anglo-Australian (AAO) observatories. All
plates were digitized at the STScI
(Space Telescope Science Institute) utilizing modified PDS type
scanning machines (15
m/pixel, i.e. 1 arcsec/pixel). This
GSC-II database, from which the public GSC2.2 catalogue was
exported, contains multi-epoch positions, multi-band magnitudes,
and classification for
1 billion objects to
mag derived from the analysis of
8000 plates (McLean
et al. 2000). Table 1 shows the main
characteristics of the material utilized for the construction of
the GSC-II. For each survey we report the observation epochs, the
photographic emulsion and filter that define the passband, and the
magnitude limit. For our survey we used 32 regions of the Northern
hemisphere
,
selected at high
galactic latitudes in order to avoid crowding, and concentrated
towards the North Galactic Pole (NGP). Figure 1
shows the distribution of the 32 regions. The total area covered
is about 1150 square degrees. For each region we used three second
epoch plates (POSS-II blue, red and infrared) and two first epoch
plates (POSS-I blue and red), and occasionally the Quick V and ER plates. To detect high proper motion objects, we chose regions
having overlapping POSS-II plates with epoch difference of
yr. Most of these plates were reprocessed
with the same pipeline used for the construction of the GSC-II, in
order to take advantage of the latest photometric calibration
available to the project.
In other cases, the data were extracted directly from the GSC-II
database (e.g. the POSS-I plates used to improve the precision of
the proper motions). Then, the object matching and proper motion
evaluation was performed using the procedure described in Spagna
et al. (1996).
Object positions were computed by means of astrometric
calibrations based on the reference stars extracted from the
Tycho-2 catalogue (Hog et al. 2000). A quadratic
polynomial was adopted to model the transformation from the plate
coordinates x, y, previously corrected for refraction, to the
standard coordinates from which object positions
were derived. Also, the typical distorsions which affect Schmidt
plates were corrected by means of astrometric masks which model
the 2D pattern of the position residuals with respect to the
reference catalog as derived from the average of about 100 plates.
The final accuracy of the absolute astrometry is better than
,
while relative astrometry has an accuracy of
.
The photometric calibrators used for GSC-II are the
GSPC-II, GSPC-I and Tycho catalogs, which sample, respectively,
the faint, bright, and very bright range of the GSC-II
photographic magnitudes (Bucciarelli et al. 2001). A photometric accuracy of 0.1-0.2 mag is
generally attained. Larger errors may be present for faint objects
close to the magnitude limit of the plate.
Table 1: Main characteristics of the plate material utilized for the construction of the GSC-II.
![]() |
Figure 1:
Distribution of the GSC-II fields (boxes) in equatorial
coordinates. Lines show the location of the NGP, the galactic
plane and the strips at
![]() ![]() |
Open with DEXTER |
Initially, relative proper motions were derived from POSS-II
plates only (epoch difference between 1-10 years) by applying the
procedure described in Spagna et al. (1996). In
practice, faint field stars (
)
have been selected
and used as reference stars to fit a third order polynomial and
transform the instrumental coordinates x,y from the secondary
plates to the reference plate, usually the one at intermediate
epoch. Then, relative proper motions
were computed
from the multi-epoch positions by means of a linear fit. Given a
position accuracy of
,
the POSS-II
plates provide proper motions with a typical accuracy,
,
of a few
hundredths of arcsec per years which is sufficient for the
identification of high proper motion objects and the selection of
our targets. For the candidates confirmed by the spectroscopic
follow-up (Table 2) final and more accurate
proper motions were computed by using multi-epoch positions
derived from a larger set of plates including the first epoch
POSS-I, the intermediate epoch Quick-V, the second epoch POSS-II
and ER (for the equatorial fields) (cf. Table 1). In order to
cross-match the objects on all these plates, the short baseline
proper motions previously computed were used to predict their
position at the different epochs. Thanks to the
40 years
between POSS-I and POSS-II plates, the typical precision of the
final proper motions is
mas yr-1 down to
.
Finally, absolute proper motions were derived by forcing the
extragalactic sources to have a null tangential motion.
As an external check we compared the final proper motions of the
targets in Table 2 with values from USNO-B. The
measurements of the two catalogues appear well correlated, a part
few outliers possibly due to wrong USNO-B solutions.
The rms on the total
proper motion difference is about
mas yr-1, after rejecting 7 objects with large proper motion
difference (
mas yr-1) and
large internal errors
mas yr-1).
![]() |
Figure 2: Final proper motion errors as a function of magnitude for one of the field towards the NGP included in this study. |
Open with DEXTER |
After selecting a reference plate, usually that of intermediate epoch, objects on different plates were matched with a variable search radius, which grew as a function of the epoch difference, so as to include objects moving up to 2.5 '' yr-1. Their relative proper motions were then derived. To minimize false candidates, we selected only objects which were matched on all three POSS-II plates.
White dwarf candidates were selected by means of various criteria.
We searched for faint (
)
and fast moving stars
(
yr-1)
.
Assuming that our survey is efficient down to
(i.e. 1.0-1.5 mag brighter than the plate limit),
we can detect ancient and faint WDs with
up to 30-40 pc. At this distance range, the imposed proper motion limits
select objects moving with tangential velocities in the interval
km s-1, which includes almost
all halo WDs and a large fraction of the thick disk WDs.
The previous criteria were implemented by means of an automatic
procedure which screens the 105 objects typically
available in each field and provides a list of a few hundred high
proper motion candidates. However, not all the selected targets
had a real proper motion due to mismatches, binaries, etc. Also,
the proper motion of the faintest objects were considered
suspicious, due to the sudden increase of the proper motion error
close to the magnitude limit of the plate. For these reasons the
proper motion of each preliminary candidates was confirmed by
visual inspection of the POSS-I and POSS-II plates. Finally, a
cross-correlation with other catalogues (2MASS, LHS, NLTT, SDSS)
was also performed.
A very useful parameter to further refine target selection
and separate candidates WDs from (sub)dwarfs is the reduced proper
motion (RPM), H, defined as (Luyten 1922)
![]() |
(1) |
The reduced proper motion, used in combination with colors, is
very powerful in separating the different population of the
Galaxy.
This is evident in Fig. 3 (left panel), which shows
the RPM diagram for one of the GSC-II regions located towards the
NGP. Here, the dots represent all the stars found in the
field (i.e. without any proper motion selection), and we notice
that they are clustered in two main populations. The group of
objects on the top right part of the diagram is composed of late
type disk dwarfs, while the population on the left with
is mostly subdwarfs.
White dwarfs are located on the bottom left part of the RPM
diagram. This classification is
confirmed by the position of the stellar isochrones and WD cooling
tracks shown in Fig. 3. The solid and dotted lines on the top
right part of the diagram show the loci of the disk dwarfs and of
the halo subdwarfs based on the 10 Gyr isochrones down to
of Baraffe et al. (1997,1998)
with [M/H]=0 and -1.5, respectively. Theoretical isochrones were
converted to photographic magnitudes BJ and
using the
color transformations adopted for the photometric calibration of
the GSC-II.
The dot-dashed, dot-dot-dashed, and dashed lines are based on the
cooling tracks for
WDs with hydrogen atmosphere
from Chabrier et al. (2000). They were obtained by
adopting tangential velocities of 38, 80, and 270 km s-1,
for the thin disk, thick disk and halo, respectively
.
The shape of the cooling tracks reflects the recent improvement of
the theory of atmosphere of cool WDs. In case of a pure hydrogen
atmosphere, the cooling sequence turns back toward the blue. This
effect is due to the collision-induced-absorption (CIA) in the
infrared produced by molecular hydrogen (H2) which gives a
redistribution of the flux toward shorter wavelength, and is
strongly present at temperatures
K (see, e.g.,
Chabrier et al. 2000).
![]() |
Figure 3: Left panel: RPM diagram for a GSC-II region located towards the NGP. Stellar isochrones and WD cooling tracks are plotted with different lines (see text). Right panel: RPM diagram for all the candidates selected in the 1150 square degrees surveyed for this study. Dots show the unobserved targets, while triangles and asterisks are the spectroscopically confirmed WDs and contaminating (sub)dwarfs, respectively. Crosses indicate stars unclassified due to the low signal to noise of their spectra. |
Open with DEXTER |
The right panel in Fig. 3 shows the reduced proper motion
diagram for all of the objects selected, in the 1150 square
degrees surveyed, after the application of the selection
procedures just described. This RPM diagram was utilized to refine
the sample of WDs candidates for spectroscopic follow-up. To
separate WDs from contaminating (sub)dwarfs, we used as boundary
line the cooling sequence for pure hydrogen atmosphere WDs (0.6 )
with the mean tangential velocity of the thin disk
fixed at the value in the direction of the NGP; this in practice,
corresponds to a kinematic selection of WDs faster than 38 km s-1. Besides a significant fraction of thin disk WDs, such
threshold includes almost all of the WDs with halo kinematics and
most of the thick disk WDs
.
We observed all the candidates lying below the boundary line; a
few objects above the threshold and some candidates with
redder than the cooling track turnoff (
)
were
also observed in order to check the efficiency of the selection
criteria and to take into account the photometric error.
In the figure, dots are used for the stars in the initial candidate sample which were later excluded from the spectroscopic follow-up list; confirmed WDs are plotted as triangles, while (sub)dwarfs are represented as asterisks. We note that all the confirmed WDs (except a couple of cases with slightly redder color) fall under the boundary line, while the contaminating (sub)dwarfs are distributed near the threshold. These results demonstrate the good quality of the data, the selection criteria and the accuracy of our final sample.
![]() |
Figure 4: Optical spectra. All the spectra are normalized at 5500 Å, shifted vertically from each other by arbitrary units, and ordered by increasing right ascension. Asterisks indicate spectra not fully calibrated, shown for classification only. |
Open with DEXTER |
![]() |
Figure 5: Optical spectra. All the spectra are normalized at 5500 Å, shifted vertically from each other by arbitrary units, and ordered by increasing right ascension. Asterisks indicate spectra not fully calibrated, shown for classification only. |
Open with DEXTER |
![]() |
Figure 6: Optical spectra. All the spectra are normalized at 5500 Å, shifted vertically from each other by arbitrary units, and ordered by increasing right ascension. Asterisks indicate spectra not fully calibrated, shown for classification only. |
Open with DEXTER |
![]() |
Figure 7: Optical spectra of the peculiar WDs. Left panel shows GSC2U J 131147.2+292348 which is the coolest DQ with extremely strong C2 bands. Right panel shows WD 1328+307 (G 165-7) which is a heavily blanketed DZ. |
Open with DEXTER |
Low dispersion spectra were measured to confirm the nature of the selected candidates and remove possible contaminants, like (sub)dwarfs.
The bulk of the spectroscopic follow-up was carried out at the 3.5 m Telescopio Nazionale Galileo (TNG, La Palma), with additional contribution of the service programmes of the 4.2 m William Herschel Telescope and the 2.4 m Isaac Newton Telescope (WHT and INT, La Palma).
We used the low resolution spectrograph DOLORES (Device Optimized
for LOw RESolution) installed at the Naysmith B focus of the TNG.
Half of our spectra were measured on the nights of April 16th-18th 2002, conditions were mostly good, the first half of the first
night suffered seeing of around 2 arcsec, but this settled
down over the rest of the run to be 1-1.5 arcsec. The
remaining objects were measured in a second run on October 2002,
in service time on February 2003 and in service time from August 2003 to January 2004. Most of our spectra were observed with the
low resolution, blue-optimized Grism 1 (or LR-B) with a long slit.
The transmission of the LR-B grism peaks at around 400-500 nm with
a central weavelength of 585.0 nm and useful coverage from 300 nm
to 880 nm and a dispersion of 2.8 Å per pixel. The camera is
equipped with a
Loral thinned and
back-illuminated CCD with 15 micron pixels. The field of view is
arcmin with a 0.275 arcsec/pix scale, hence a 1 arcsec slit projects to 3.6 pixels, or a resolution of
11 Å. Most of our targets were observed with a 1.0 arcsec
slit; on April 15th when seeing was poor we employed the 1.5 arcsec slit and the 2.0 arcsec slit. Each target was
pre-imaged to define its exact position (helpful with such high
proper motion stars). Typically, an exposure time of 3600 s and 300 s was sufficient to provides a suitable signal-to-noise
ratio,
,
for the faintest (
)
and the
brightest (
)
target respectively.
Spectrophotometric standards were observed a few times per night
to enable an approximate flux calibration, but most importantly to
help get the shape of the spectrum correct. The standards were
drawn from the list at http://www.ing.iac.es/Astronomy/observing/manuals/html_manuals/tech_notes/tn065-100/workflux.html
We measured the stars GD 140, BD+17 4708, Feige 67 and HZ 44. We obtained Helium arc lamp exposures to facilitate wavelength calibration. Halogen lamp flat fields were measured to enable us to measure the pixel-to-pixel sensitivity variations.
Spectra for 4 objects in Tables 2-4 were measured using the intermediate dispersion spectrographic and imaging system (ISIS) on the 4.2 m WHT. Observations were made during service nights on the 27th and 29th January 2001, 6th December 2003, and 3rd and 4th February 2004. ISIS is a two-arm spectrograph, employing the 5700 Å dichroic to split the blue and red light into two cameras. In the blue arm we used the R300B grating with the EEV12 CCD (4096 pixels in the wavelength direction) at a central wavelength of about 4400 Å. This setup achieves a dispersion 0.86 Å/pix and useful wavelength coverage of about 3000 Å. The slit width was matched to the seeing, typically 1 arcsec projecting to 4 detector pixels. In the red arm we used the R316R grating with the Marconi2 CCD (4610 pixels in the dispersion direction) at a central wavelength of about 6250 Å. The dispersion is 0.84 Å/pix covering 3858 Å, but around 1000 of these pixels suffer significant (50-80%) vignetting from the camera optics. Combining the red and blue arm spectra gave us spectra with useful coverage over the wavelength range 3400-8000 Å.
Spectra for 2 of the objects listed in the Tables 2-4 were measured in service
time on the 2.5 m INT. We employed the (currently decommissioned)
Intermediate Dispersion Spectrograph (IDS) mounted at the
Cassegrain focus. Observations were made during the nights of 8th April 2003 and 19th April 2003. The instrumental setup was the
same for both nights, using the 235 mm camera with the R150V grating and the EEV10 CCD. This grating and camera combination
gives a dispersion of 271.3 Å mm-1, which is equivalent to 3.66 Å/pixel for the EEV10 detector, and useful wavelength
coverage was from approximately 3300-8700 Å. The camera and
detector achieve a spatial scale of 0.4 arcsec per pixel. On
the first night, the slit was opened up to 2.5 arcsec (6 pixels). On the second night we observed with a 1.0 arcsec slit
(2.5 pixels). A total of 4 targets were measured on the first
night, with exposure times of between 1200 and 2700 s each,
as well as one standard star, SP1550+330. On the second night we
observed 2 targets, each with
s exposures, as well as
the flux standard SP1446+259. Calibration arc lamps and flat
fields were measured on both nights.
Data reduction for all spectra was performed within the IRAF and MIDAS data reduction package. The same procedure was followed for all objects wherever possible, regardless of the instrument and telescope used to acquire the data. All images were debiased using a mean zero second frame. Images were then flatfielded using lamp observations to remove the pixel-to-pixel sensitivity variations. 1-D spectra were extracted and background subtracted, then wavelength calibrated using arc lamp spectra.
Most of the spectra were flux calibrated using the standard star observations, although no correction was made for extinction. 16 spectra have not been flux calibrated at this time and are presented in Figs. 4-6 in their uncalibrated form. Thus the turn-over to the blue of these objects reflects the sensitivity of the instrument, rather than the intrinsic spectral shape of the object. This does not prevent us from identifying that these stars are WDs, but does caution against detailed spectral analysis at this time, which will be completed and presented in a future paper (Carollo et al., in preparation).
We secured spectra for 59 candidates: 40 were confirmed WDs, 16
were a mixture of late type dwarfs and subdwarfs, while the
spectra of the last 3 stars failed to reduce due to poor data
quality.
A couple of candidates below the boundary line (see Sect. 3) were
not observed for technical reasons. One of them (LHS 1093) is a
known WD (Bergeron et al. 1997). Thus, our sample
counts 41 WDs in total, 17 of which show no presence of the H line, and two are exotic examples of DQ and DZ.
Table 2 gives the astrometry of the WD sample; Col. 1 reports the GSC-II identification, while the remaining columns list: positions with accuracy in the range 0.2 ''-0.3'', absolute proper motions with a typical total error between 2 and 6 mas yr-1, and the epoch of the corresponding position.
Table 3 provides BJ and
photometry,
spectroscopic classification, and source telescope. Magnitudes are
between
and colors in the range
.
The last column of Table 3 lists
the cross-identification with other WD catalogues. Only 5 objects
are in common with the SDSS sample of cool WDs (Kilic et al. 2005), while 12 WDs were found in the Villanova
catalogue (McCook & Sion 1999) available at CDS
(http://cdsweb.u-strasbg.fr/). Therefore, our sample counts 24 newly discovered WDs.
Finally, Table 4 reports the cross-match with
Luyten's catalogues, the Data Release 4 (DR4) of the SDSS
(Adelman-McCarthy et al. 2005), and the
2MASS catalogue (Cutri et al. 2003); Col. 2 gives
the NLTT/LHS identification. In the remaining columns we provide
ugriz photometry for 19 WDs found in SDSS-DR4, and JHKmagnitudes for 19 objects found in 2MASS.
Figures 4-6 show the optical spectra of the WDs in the
GSC-II sample, except for the two peculiar WDs, plotted
in Fig. 7. Of the 41 WDs, 17 are featureless (DC) and 19 are DA. Two stars in the sample (LHS 5222 and LP
618-6) show carbon features and one of them (LP 618-6) is magnetic
(Schmidt et al. 2003).
Two stars (GSC2U J030339.2+140805 and LP 618-14) show the features of a composite spectrum which indicates the presence of a red companion. A third star (LP 618-1) is an unresolved degenerate double formed by a magnetic DA and a massive DC (Bergeron et al. 1993).
Also, we found 2 common proper motion systems composed of WDs only (LHS 5023 + LHS 5024 and LP 642-52 + LP 642-53). White dwarf GSC2U J071859+551406 is the faint component of another common proper motion pair with G 193-40 as the bright companion.
The left panel in Fig. 7 shows the coolest (5120 K) carbon rich WD known to date: GSC2U J131147.2+292348 described in Carollo et al. (2002,2003). Another similar DQ was recently discovered in the SDSS (Kleinman et al. 2004); its modelling by Dufour et al. (2005) yields a temperature of 5140 K, which increases to two the DQs with temperature below the expected cutoff of 6000 K.
The right panel shows the spectrum of GSC2U J133059.8+302955 (G 165-7) which is a peculiar DZ with an estimated temperature of 7500 K (Wehrse & Liebert 1980).
We have presented a new high proper motion survey to search for
ancient cool WDs. The survey is based on the GSC-II material and
covers an area of about 1150 square degrees. Using appropriate
selection criteria, the RPM diagram, and spectroscopic follow-up,
we identified 41 WDs, 17 of which show no H line, and two
of peculiar nature.
This new and independent WD sample will be used to infer new constraints on the halo WD space density (Carollo et al., in preparation) by means of a new method to separate thick disk from halo WDs based on the kinematic properties of the candidates (Spagna et al. 2004). Furthermore, this sample will be useful to increase the statistics of such rare objects and provide information on the age of thick disk and halo through the analysis of the luminosity function.
Acknowledgements
The GSC-II is a joint project of the Space Telescope Science Institute and the INAF-Osservatorio Astronomico di Torino. Space Telescope Science Institute is operated by AURA for NASA under contract NAS5-26555. Partial financial support to this research cames from the Italian CNAA and the Italian Ministry of Research (MIUR) through the COFIN-2001 program. This work is based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Centro Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias, and also based on observations made with the William Herschel Telescope and the Isaac Newton Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias.Special thanks to C. Loomis, A. Volpicelli and A. Zacchei for their valuable technical support to this project and to D. Koester for the many useful discussions. Finally, B.B. and M.G.L. acknowledge the support of the STScI through the Institute's Visitor Program for 2005.
Table 2: The sample of cool WDs from the GSC-II survey: astrometric parameters.
Table 3: The sample of cool WDs from the GSC-II survey: photometry, classification, source telescope, WD ID in other catalogues.
Table 4: The sample of cool WDs from the GSC-II survey: cross-match with Luyten's catalogues (LHS and NLTT). Photometry from SDSS and 2MASS.