skip to main content

gaia data release 3 documentation

7.5 Orbital solutions from spectroscopy: SB2

7.5.3 Underlying orbital models

Four kinds of orbital solutions are considered here and they are labelled SB2, SB2C, TrendSB2 and StochasticSB2. We describe them below.

Orbital models of type SB2

The general eccentric Keplerian model is, for the spectroscopic channel, expressed in terms of the Campbell coefficients (see the appendices in Halbwachs et al. 2023). The model radial velocity is given by, for one of the component labelled 1,

RV1(t)=γ+K1[cos(v(t)+ω)+ecos(ω)] (7.36)

where v(t) is the true anomaly which is deduced from the eccentric anomaly E by

tanv2=1+e1-etanE2 (7.37)

which in turn is expressed as a function of the mean anomaly M

M=E-esin(E)=2πP(t-T0) (7.38)

where T0 is the time of passage at periastron. The ω is the longitude of periastron (for component 1) measured in the sense of the orbit with its origin at the ascending node. In this model, the eccentricity e enters as a very non-linear parameter. Concerning the other component (2), the curve can be modelled by

RV2(t)=γ-K2[cos(v(t)+ω)+ecos(ω)]. (7.39)

The parameters to determine are P, γ, K1, K2, e, ω, T0 (respectively, the period, the centre-of-mass velocity, the semi-amplitude of component 1, the semi-amplitude of component 2, the eccentricity, the longitude of periastron and the time of passage at periastron). The first step of the computation is devoted to the determination of the period. Some significance tests are made on the very existence of the period. If the period is not significant, the model is fitting noise and is considered as being invalid. This is a statistical decision. It is interesting to note that in general

RV1-γRV2-γ=-K1K2. (7.40)

From this property, one deduces that in a RV2 versus RV1 diagram, the data should occupy a straight line whose slope is function of the mass ratio (noise effects excepted).

Orbital models of type SB2C

Another particular kind of orbit is the circular one, i.e. e=0. In this case, the model is much more simple and is expressed by the following equation

RV1(t)=γ+K1[cos(2πP(t-T0))] (7.41)
RV2(t)=γ-K2[cos(2πP(t-T0))] (7.42)

where here T0 is the time of maximum velocity of the first component and the ω is fixed at zero corresponding to the extremum velocities. The free parameters to determine are restricted to P, γ, K1, K2, T0. The model is linear except, of course, for the determination of the period and for the constraint that RV1 and RV2 should be in antiphase. The significance of the period is here also determined.

Models of type TrendSB2

The equivalent of the TrendSB1 (see Section 7.4.3) formalism can also be extended to the difference RV2-RV1. However, no SB2 object entered this category for Gaia DR3 and this model is not further detailed.

Models of type StochasticSB2

When the fit of a TrendSB2 solution is not possible and the fit of an orbital solution corresponding to a significant period neither, the solution is classified as StochasticSB2. The RV is considered as constant but with an extra dispersion compared to the RV epoch uncertainties. This extra dispersion is estimated for the difference RV2-RV1, if possible.