A&A 404, 317-323 (2003)
DOI: 10.1051/0004-6361:20030421
R. L. Smart - M. G. Lattanzi - B. Bucciarelli - G. Massone - R. Casalegno - G. Chiumiento - R. Drimmel - L. Lanteri - F. Marocco - A. Spagna
Istituto Nazionale di Astrofisica (INAF), Osservatorio Astronomico di Torino, Strada Osservatorio 20, 10025 Pino Torinese, Italy
Received 21 January 2003 / Accepted 17 March 2003
Abstract
We present parallax determinations for six white dwarf
candidates in the Torino Observatory Parallax Program. The absolute
parallaxes are found with precisions at the 2-3 milliarcsec level.
For WD 1126+185 we find a distance incompatible with being a white
dwarf, implying an incorrect classification. For WD 2216+484 we find
our distance is consistent with a simple DA white dwarf rather than a
composite system as previously proposed in the literature. In general
it is found that the published photometric distance is an overestimate
of the distance found here.
Key words: white dwarfs - stars: distances - astrometry
In this article we discuss the instrumentation, observational procedures, calibrations, reduction techniques and first results for six white dwarf candidates of the TOPP list.
This telescope, filter and CCD combination was fixed, giving us a stable system which is required for astrometric studies of this nature. The TOPP was provided with 40% of the observational time on this telescope during the period 10 November 1994 to 2 July 2001. In Fig. 1 we plot the change in temperature, humidity, seeing and number of frames observed over the course of this period. In July 2001 the CCD was changed and the program for these targets was halted. The telescope is presently undergoing automation refitment and when that is stable we will begin another program with the new, higher resolution, CCD.
![]() |
Figure 1: The nightly number of frames taken, seeing, humidity and temperature conditions. Day number is from the 11 November 1994, i.e. Julian date 2 449 667.0. |
Open with DEXTER |
During the observational period the only major maintenance was a periodic realuminisation of the mirrors on 14/04/1997, 24/05/1999 and 04/12/2000. This required a dismounting of the mirrors and transportation to the Asiago Observatory. This procedure should not have an impact on the parallax program but we search for correlations with residuals in the solutions and these dates as a matter of routine.
For each target we make one 30 s exposure, then measure both the position and maximum pixel count for the target. The telescope is then repointed so that the target is in the center of the CCD and the exposure time is calculated in a simple manner to optimize the counts for the target and reference stars. In this way we insure that spatial variations on the CCD have minimal effect and that nightly variations of the sky are to a reasonable extent accounted for. If the target is not too bright or faint and the reference stars not too sparse, then we choose a maximum pixel count that on a bad night will not require more than a five minute exposure and on a good night defaults to this optimal exposure time. We attempt to restrict all observations to within 30 min of the meridian and to make two exposures for each target per night.
This object list and cleaned image is then processed to find the most
important measured quantity in this project - the centroid. In Smart
et al. (1999) we discuss the experiments carried out
to determine the image centroids and conclude there that the best
procedure is an in-house two-dimensional Gaussian profile fitter. We
estimate by comparison of consecutive frames that the centroiding
precision is better than 0.02 pixels (10 mas) for stars within 5 magnitudes of saturation and decreases rapidly for fainter stars.
The largest source of systematic error calibrated in MON92 was the effect of Differential Color Refraction (DCR). This occurs because the frame defined by the reference stars has a different average color than the target star. The position of the target star will be refracted differently than the mean field by the atmosphere and the size of this difference will be a function of the hour angle. This is especially important when observations were carried out far from the meridian to maximize the parallax factor.
In MON92 this effect was discussed at length and we are currently using data taken over the course of the program to calibrate the effect in a similar way. However we note that the combination of our CCD and longer wavelength filter means that our DCR will be much lower than theirs. In Fig. 2 we plot the combined CCD/filter/telescope efficency. In Smart et al. (1999) we show that the DCR of the Torino program is negligible in the case of QSOs which have similar colors to these white dwarfs and are close to the mean color of the reference stars. As a result of our low DCR, the restriction that observations be made close to the meridian, and a low color difference between the target and the mean reference star color, we believe we are justified in ignoring the DCR correction for these particular objects. Future work will include the DCR.
![]() |
Figure 2: The combined CCD quantum efficiency, optical transmission and filter transmission for the USNO Flagstaff 1.5 m reflector and the OATo 1.05 m reflector. |
Open with DEXTER |
For the relative parallaxes found in this work we have used a simple iterative approach. We will eventually use a more rigorous procedure that will find directly the absolute parallax using a global central overlap procedure similar to that described in Eichhorn (1997). Since the stars under discussion here have parallaxes of tens of milliarcseconds we expect refinements such as the use of the central overlap to have minimal effect.
The selection of frames and reference stars is completely automatic. We have the ability to manually select the reference stars but testing on simulations has shown that an intelligent choice of selection criteria is sufficient and has the advantage of being free of observer bias. Initially all frames with an hour angle greater than one hour and all frames used to set up the target are dropped. A frame from the middle of the sequence with a high number of stars is chosen as the base-frame.
From the measured coordinates (x,y) of objects on the base frame we
calculate standard coordinates ()
by a linear transformation
using stars common to a reduction of the Digitized Sky Survey
(DSS). This sets the plate scale and orientation which allows us
to use parallax factors in standard coordinates. The precision of this
transformation is limited by the quality of the DSS solution; for
the orientation any errors less than 1 degree do not affect our use of
the standard coordinate parallax factors and we assume any errors in
the scale are less than 1%.
Using pattern recognition software we find 6 stars in common with the base frame and the other frames; any frame without at least six common stars is dropped. We then orient each frame to the standard coordinate system of the base frame using these 6 stars. Then we find all stars in common with the base frame and pick the 12 most frequent stars including the target star, all frames without these stars are dropped. These 12 stars define the reference system for the target star. Obviously if all the frames have a larger subset of stars in common then we will use that, but usually to have at least 12 stars we have to drop some of the frames.
We then find the frame parameters A, B, C using the common stars by
solving the transformation:
Then for all stars on the base frame we find a preliminary parallax
and proper motion using:
From this fit we have predicted positions of the reference stars which we can then use to iterate again through Eqs. (1) and (2). Before doing this iteration we attempt to identify bad frames and stars by finding the root mean square of the residuals per frame and per star. Every frame or star with a root mean square greater than three sigma from the mean we delete from the next iteration.
We iterate for at least 4 iterations and then stop when the change in the parallax is less than one mas. We attempted various forms of reweighting as we iterated both in the calculation of the plate parameters and in the calculation of the final parallax - however we found the best results by only applying weights in the calculation of the plate parameters in the final iteration. We have also reduced the sequences using Gaussfit (Jefferys et al. 1988), constraining the parallaxes and proper motions of the reference stars to zero, and this provide similar results, increasing our confidence in the solutions.
Using this reduction procedure we have implicitly assumed that the proper motions and parallaxes of the reference stars are uncorrelated within the frame. This will not be true for the proper motions, which in a given direction will be dominated by a combination of differential galactic rotation and solar motion, but this will only add a systematic error to the proper motions and we can ignore it as our main goal here is the parallax.
The quantity calculated above is the parallax of the target star
relative to the anonymous reference stars. For each reference star we
have also calculated a parallax and proper motion using the same
reference system as that of the parallax star. Inevitably the
reference star parameters (
,
)
are
correlated with the definition of the reference system; however with at
least 12 stars in each system we consider this correlation to be
negligible.
From the target star's relative parallax we should subtract the mean parallax of the reference stars weighted by the weights used in the plate solutions and the dependences (Schlesinger 1911). This enforces the assumption that the reference frame is at a mean parallax of zero and the use of weights in this step allows us to consider each reference star as having unit weight in the final step for the correction from relative to absolute parallax. However, the dependences are all nearly equal and we expect this mean to be zero. We have found that usually the mean is much smaller than the variance of the various parallaxes and we have decided not to apply this correction as it is not sensibly different from zero.
As part of the calibration we observed standard magnitudes for all stars in each field. Using the Mendez & van Altena (1996) galaxy model we can calculate the most probable distance of each reference star based on their magnitude. The mean of these distances is an estimate of the correction to absolute parallax that we must add to the calculated targets relative parallax. In Table 1 we list the range and mean of the reference star magnitudes and the corresponding relative to absolute correction (RAC) for each field.
Table 1: Magnitude ranges and means in the I filter, and relative to absolute corrections in milliarcseconds as found from the Mendez & van Altena galaxy model.
After the final iteration the formal error of the star parameters come directly from the final covariance matrix. As we used a nominal input for the centroiding error we multiply the terms of the covariance matrix by the error of unit weight.
Based on various considerations (see for example Smart et al. 1997) we estimate the error in the RAC found using a galaxy model to be about 30%. Therefore to the formal relative parallax error above we add in quadrature 30% of the RAC to provide us with a final error.
Table 2: Parallaxes and proper motions for the six white dwarf candidates in the Torino Program. See text for explanation of columns.
We note that the use of the DSS to transfer the frame to the sky is basically a setting of the scale and orientation of the CCD. From an examination of the plate parameter errors and tests against external catalogs we have found that the scale of the DSS is correct to better than 1%. This error is invisible to the calculation of the relative parallax and cannot be easily estimated on a case by case basis. We consider it a systematic error that for all of the targets here will be below the value of the formal error (being at most 1% of the targets parallax) and therefore if the user requires a maximum error estimate we suggest that they add 1% of the targets parallax to the quoted error.
The error on the proper motions is more difficult to access as it is affected directly by the orientation, scale, and the fact that they are relative as we have made no attempt to account for systematic motions of the reference stars. Any use of the proper motions quoted here that requires high precision should be done with care.
![]() |
Figure 3: Observations of 0423+120 along with the fitted path. The X marks the position of the base plate. |
Open with DEXTER |
In Fig. 3 we present as a typical example the observed
and fitted path for the target 0423+120. The average residual is 12
and 10 milliarcsec in
and
respectively.
Two targets, 1716+020 and 1840-111, are part of wide binary systems with M-dwarf stars. In both cases we find the difference in proper motions measured in these pairs are less than the errors. As can be seen from Table 3 both these systems have parallaxes in good agreement. The systematic direction of the existing difference, where the red star has larger parallax, is consistent with what is expected by not applying the DCR. Nevertheless the difference is of the order of the errors justifying our neglect of the DCR for these objects.
Table 3: Parallaxes of two binary systems.
Table 4: Observations from the Fourth Yale Parallax Catalogue for the determination of the parallax of 1716+020.
0322-019 (LHS1547) is listed in the Astronomisches Rechen-Institut
database for nearby stars (http://www.ari.uni-heidelberg.de/aricns/index.htm, hereafter ARICNS)
as a DZ9 with a photometric parallax of 62 mas, consistent
with our measurement of
.
0423+120 (EJ 169) is listed in the ARICNS as a DC8 with a photometric
parallax of mas, significantly different from our value of 58 mas. This is derived from the absolute magnitude of McCook & Sion
(1999, hereafter MCC99) and the apparent magnitude in the
ARICNS.
1126+185 (GJ 3667) is listed as a DC9 with a photometric
parallax of mas, inconsistent with our zero value.
We have used the finding chart provided in
Green et al. (1986) where it was first identified
as a composite DC white dwarf based
on at least two spectra which were not published. Subsequent work with
new spectroscopy by Putney (1997, Fig. 2l) also
confirmed this result with a spectra that indicated "a composite of a
DC and a cool star''. We have discussed these observations with
various spectroscopists in light of this object's near zero proper
motion and parallax, but have not been able to come up with other
suggestions. As no other star in the TOPP field has a significant
proper motion we therefore assume this is a case of
misclassification and/or misidentification.
1716+020 (YPC 3926.01A, LHS3278) is listed in the Fourth Yale
Parallax Catalog (van Altena et al. 2001) with a mean parallax of mas and a spectral type of DA6 which is inconsistent with
our value of
mas. However, as we noted before,
this object has a companion red dwarf (LHS3279) and to calculate the
Yale mean parallax the various observations of both stars were
averaged. In Table 4 we list the individual observations.
Considering 1716+020 alone, our measurement is consistent with its
mean parallax, whereas inconsistency arises when the parallax of the
red dwarf companion is considered. This difference may be due to an
incorrect treatment of the refraction effect in the other programs
where it is much more important than ours. We also note that the
projected distance in all cases is less than 0.016 parsecs while only
our values can produce distances that are consistent with a line of
sight difference of this order.
1840-111 (GJ 2139) is listed in the ARICNS as a DA5 with a
photometric parallax of mas, inconsistent with our value of 35 mas. This is a derived value from the MCC99 work and is discussed
below. The higher error of this target may be a result of it's lower
declination combined with the high field density, which could adversely
affect the centroiding by the crowding when the seeing was not
optimal.
2216+484 (GD 402) has no parallax estimates in the literature. In Bergeron et al. (1990), based on anomalies in the spectra, they conclude that this object is probably a composite system of a DA and a DC white dwarf. From the apparent magnitude and our parallax we derive an absolute magnitude of 13.76 in V which is consistent with a single DA star and contrary to the thesis of a double system which is estimated in Bergeron et al. to be of absolute magnitude 12.94. An examination of the point spread function during good seeing does not show any indication of duplicity. The apparent magnitude and the astrometric residuals of the system are stable over the program. These observations are not conclusive proof of singularity but they do put restrictions on orbit orientation in the case of a composite system.
Table 5: Absolute magnitudes calculated from TOPP parallaxes and apparent magnitudes along with those derived in MCC99 from spectrophotometry, the code indicates the method used: 2 = multichannel spectrophotometric colors, 3 = uvby colors.
In Table 5 we list the absolute magnitudes derived from the parallaxes and apparent magnitudes quoted here along with those derived in MCC99 from spectrophotometry. We have not included 1126+185 and we note that the 2216+484 MCC99 value is based on the assumption that it is a single DA7 star. It is beyond the scope of this paper to discuss the relative merits of absolute magnitudes derived from spectrophotometry versus those from parallaxes and the sample here is too small for such a comparison. However, there is a suggestion that the absolute magnitudes derived from uvby colors are underestimates (i.e. numerically larger) than the true magnitudes - assuming that there is not a systematic error of 30% in our results.
Apart from 1716+020 we do not have any direct comparisons to other trigonometric parallaxes determinations. We can however make an indirect comparison. In Fig. 4 we plot the white dwarf tract in the H-R diagram as derived from both HIPPARCOS stars and those published in MON92 with the TOPP results over-plotted. We believe this adequately shows that our results are consistent with other high precision parallax determinations.
![]() |
Figure 4: The white dwarf sequence in the H-R diagram. The asterisks are from MON92, triangles from HIPPARCOS and pluses with error bars this paper. |
Open with DEXTER |
As we have stated the reduction technique used here is not optimal. The central overlap method developed by Eichhorn (1997) and extended by Rapaport (1998) is preferred. This method will allow the simultaneous adjustment of frame parameters, star parameters and the correction to the absolute frame using all the information available. It will mean that we will not have to use formal errors as weights, assume a mean parallax for the reference frame or neglect the dependencies. Significant external information is required for this procedure and it will probably only be important for stars which have parallaxes much closer to our precision limit where the relative error is larger. We do plan to implement this procedure in the future development of this program.
We have shown that the Torino program is able to produce useful parallaxes with precisions and in a magnitude range that complement the HIPPARCOS results and the USNO program. We plan to improve the reduction procedure which will only improve the precision and produce parallaxes for the other 100+ targets observed over the last 6 years. Once the automatization of the REOSC telescope is complete we expect to be able to continue the program looking at the many exotic objects being discovered by the large all-sky surveys in progress.
We have measured the first parallax for five of six white dwarfs
candidates. One object, 1126+185, is obviously misclassified,
another, 2216+484, is inconsistent with previous observations
that suggest it is a composite system. The overall photometric
magnitude calibration of McCook & Sion (1999) appears
to overestimate the objects magnitude. These discrepancies with
published work highlight the need for more direct parallax
determinations which are fundamental tests of current theories.
Note added in proof: Since going to press the ARICNS parallax for WD 0423+120 has been updated with a new photometric parallax of 52 mas, consistent with the value derived here (Jahreiss personal communication).
Acknowledgements
First thanks go to the Consiglio Nazionale delle Ricerche for initially providing the CCD, and to the former Consiglio per la Ricerca Astronomica current INAF and the Directors of OATo for their longterm support of this program. A six year observational program could not have been carried out without the support of a number of the OATo staff. In particular we wish to thank Massimo Villata, Francesco Porcu, Fillipo Racioppi, Leonardo Corcione, Marco Del Bò, Francesco Salvati and the QSO monitoring group for observations made in service mode. Thanks to Bill Jefferys and Michel Rapaport for many useful discussions, also Dave Monet and Conrad Dahn for providing us with an initial kick start for this project. We would also like to thank the referee, Dr. van Leeuwen, whose comments helped improve the clarity of this paper. Finally, RLS is grateful to Heinrich Eichhorn who before his death was a constant source of guidance and through his work continues to inspire today. RLS acknowledges the past support of the Royal Society and Torino University.