Online material
Appendix A: Seismic and orbital data
Frequencies for Star A.
Frequencies for Star B.
Mode linewidths, mode heights and mode amplitude for Star A.
Mode linewidths, mode heights and mode amplitude for Star B.
Speckle interferometric data of the orbital position of the binary.
Appendix B: Astrometric orbit
The coordinates of the orbit on the plane of the sky (x,y) are derived as follows: where A,B,F,G are the ThieleInnes elements, and (X,Y) are the position of the object in the plane of the orbit. The ThielesInnes elements are related to the semimajor axis a, to the inclination of the orbit i on the plane of sky and to the argument of the periastron (ω) and the longitude of the ascending nodes (Ω) as: The position of the orbit in the plane of the orbit is given by: where E is the eccentric anomaly and e is the orbit eccentricity. The eccentric anomaly can be found by solving the following equation: (B.9)where T_{0} is the orbit periodicity and T_{1} is the time of passage at periastron.
Appendix C: Derivation of the orbit
For the derivation of the orbit, we computed the log likelihood of the data D given the orbital parameters (=(Ω,ω,i,a,e,T_{1},T_{0})), as: (C.1)where D are the data, and are the output of the astrometric model of Appendix A at time t_{i}, and and are the observations at time t_{i}, and the σ_{i} are the estimated errors of the observations. There are two approaches used for deriving the orbital parameters:

Frequentist: minimisation of the log likelihood

Bayesian: Monte Carlo Markov Chain using the Metropolis Hasting algorithm.
The first approach is equivalent to the use of a nonlinear least square fit. The error bars are derived using the inverse of the Hessian matrix. This is the classical way of estimating parameters from a model and observation. The error bars were also verified using a MonteCarlo simulation of the orbit. When the estimated parameters are close to a minimum, the error bars from the Hessian and the MonteCarlo simulations are similar (See Appourchaux 2014).
Unfortunately, the error bars returned by the two methods gave inconsistent results which is why we implemented the second approach which is based upon a Bayesian framework. In that case, we derived the posterior probability () of the parameters using Bayes’ theorem as: (C.2)where is the a priori probability of the orbital parameters, and P(D) is the global normalisation likelihood. The derivation of the posterior probabilities can be done using the Metropolis Hasting algorithm (See as a starting point, Appourchaux 2014). We use a Markov Chain for exploring the space to go from a set to another set , assuming that either set have the same probabilities, i.e. . The MetropolisHasting algorithm then requires that we accept the new set using the following ratio: (C.3)This is simply the ratio of the likelihood given in Eq. (B.1). The new value is accepted if r<α with probability α (drawn from a uniform distribution) otherwise it is rejected.
We set 10 chains of 10 million points each, with the following starting points taken randomly: (C.4)The new set of parameters is computed from a random walk from the previous value as: (C.5)where is given by a multinomial normal distribution with independent parameters: (C.6)where α_{rate} is an adjustable parameter that is reduced by a factor 2 until the rate of acceptance of the new value t′ is above 25%. The proper convergence of the chains was verified using the GelmanRubin test (Gelman & Rubin 1992) as implemented by Ford (2006). After rejecting the initial burnin phase (10% of the chain), all values of the test of the 7 sets of parameters were below 1.1. Then the chains provide the posterior probability for each parameter. For all chains of each parameter, we computed the median and the credible intervals at 16% and 84%, corresponding to a 1σ interval for a normal distribution. The advantage of this
percentile definition over the mode (maximum of the posterior distribution) or the mean (average of the distribution) is that it is conservative with respect to any change of variable over these parameters.
Appendix D: Derivation of the mass of the binary system and the associated errors
Knowing the distance (via the parallax π), the semimajor axis a and the period of the system T_{0}, we can deduce from the Kepler’s third law the total mass of the binary system in units of the solar mass as (D.1)where π is in mas. The error bars can be computed assuming that the semimajor axis, the orbital period, and parallax are independent of each other as: (D.2)In order to derive the corresponding credible intervals for the mass of the system, we use a MonteCarlo simulation using the chains for the semimajor axis and the period from our Bayesian analysis, and use a randomised parallax as inferred from van Leeuwen (2007) for HIP 93511. The use of the MCMC chains explicitly includes the correlation between the orbital period and the semimajor axis for the final error propagation. The three values were then injected in Eq. (D.1) for getting the median and the credible intervals.
© ESO, 2015