+ | + | + | + | + | + ... = ? | Eq. (0) |
ABSTRACT: We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch-Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral-energy distribution, and even complete light curves. We use a modified version of the Wilson-Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N-body integration.
If one has all kinds of observations at disposal, a joint A χ^{2} ($\chi^2$) metric and an optimisation algorithm (a simplex or simulated annealing) allows to search for a global minimum and construct very robust models of stellar systems. At the same time, our N-body model is free from artefacts which may arise if mutual gravitational interactions among all components are not self-consistently accounted for.
Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening).
The model description was published in Broz (2017), ApJS 230, 19. Moreover, an unoffcial Appendix C briefly describes dynamical effects that can be expected in N-body simulations. Yet another useful reference is Nemravova et al. (2016), A&A 594, A55.
DATA_20160413_GAMMA.tar.gz DATA_20160331_CLOSUREPHASE.tar.gz DATA_20160415_AMBER.tar.gz DATA_20160524_SYNTHETIC.tar.gz DATA_20160612_UBV.tar.gz DATA_20170802_TELLURIC.tar.gz xitau_simplex_chi2_GAMMA_201604131719.tar.gz xitau_simplex_chi2_CLOSUREPHASE_201604131705.tar.gz xitau_simplex_chi2_PERIOD_201704081157.tar.gz xitau_simplex_chi2_HEC88_201708041001.tar.gz swift_bs_fp_201606032142.tar.gz old_versions/
Figure 1: A χ^{2} ($\chi^2$) vs the number of iterations, showing a smooth convergence of simplex to a local minimum. We can distinguish individual contributions to χ^{2}, namely astrometry (SKY), radial velocities (RV), minima timings (TTV), eclipse durations (ECL), squared visibilities (VIS), closure phases (CLO), triple product amplitude (T3), light curve (LC), synthetic spectra (SYN), spectral-energy distribution (SED), and additional mass constraints (MASS; actually not visible). |