next previous
Up: Mass determination of

4. The specific processing

4.1. Synthetic presentation

We consider in this section the adaptation of the Hipparcos data processing to the case of short period astrometric binaries, in particular when the hippacentre and the photocentre are not alike. In this case, the reduction consists of producing the astrometry of the centre of mass and of determining the hippacentre path, yielding both tex2html_wrap_inline1601 and B. Otherwise, it yields only the scale tex2html_wrap_inline1605 of the photocentric orbit.

Let S be a program star of the Hipparcos mission. The main goal of the standard astrometric processing was to determine accurately the position of the star with respect to the centre of mass C of the solar system. When S is a single star, the position vector CS at epoch t may be modelled with the usual five astrometric parameters:
equation349
corresponding to the coordinates of position, the parallax, and the two components of the proper motion. Although our treatment makes use of the ecliptic reference frame, this choice is of no consequence in the following.

In the case of a short period binary, the instantaneous observation is tied to the hippacentre H and we have,
equation355
where G is the centre of mass of the binary and H the hippacentre, which is the same as the photocentre for very close binaries. The motion of G is well described by the five astrometric parameters tex2html_wrap_inline1625, while tex2html_wrap_inline1627 is given by Eq. (7). If everything is known on the double star, tex2html_wrap_inline1629 can be computed at any time and the raw observations duly corrected for the duplicity effect and subsequently solved for the astrometric parameters of the centre of mass. If only the relative orbit is known, one must adjust in a single process both the astrometric parameters of the centre of mass and the physical parameters tex2html_wrap_inline1631 and B of the binary. For separation less than tex2html_wrap_inline1635 0tex2html_wrap1641 3, the method will fail as the hippacentre and the photocentre merge and only tex2html_wrap_inline1639 could be determined.

4.2. Practical implementation

The equations relating the observations to the astrometric and physical parameters are all non-linear. An approximate value of the position, proper motion and parallax can be used to linearize the set of equations in the vicinity of this values. No particular problem arises in this step. As for the physical parameters tex2html_wrap_inline1643 and B one must be much more careful, and in most case the plain linearization is not sufficient and several tricks must be implemented to ensure the convergence during the iterations, and that the solutions we arrived at satisfies the non-linear system. With these safeguard we observed that the solution was not sensitive to the quality of the initial values of tex2html_wrap_inline1647 and B and that the non-linearity problem was overcome. The processing consisted primarily in the resolution of an over-determined system of N equations (the number of observations) with five or seven unknowns, by a weighted least squares. The various steps are outlined in Fig. 6 (click here).

As the magnitude difference is generally fairly well known, the input value of tex2html_wrap_inline1653 is not totally arbitrary. Consequently, a first guess of B can also be based on the mass luminosity relation for dwarf stars.

  figure366
Figure 6: Schematic principle of the specific processing of the short period binaries. In an ideal case a) where the intensity and mass fractions are perfectly known, we just use this information to shift the observations to the centre of mass and then solve for the astrometric parameters. Whenever it is not feasible b) one uses a non linear model with 7 unknowns (5 astrometric parameters for the centre of mass + 2 physical terms)


next previous
Up: Mass determination of

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr