# Physical parameters and ±0.2% parallax of the detached eclipsing binary V923 Scorpii

T. Pribulla, A. Mérand, P. Kervella, C. Cameron, C. Deen, P. Garcia, M.

Horrobin, J. Matthews, A. Moffat, O. Pfuhl, et al.

► To cite this version:

T. Pribulla, A. Mérand, P. Kervella, C. Cameron, C. Deen, et al. Physical parameters and ±0.2% parallax of the detached eclipsing binary V923 Scorpii. Astronomy & Astrophysics - A&A, 2018, 616, pp.A49. 10.1051/0004-6361/201730673. hal-02302455

HAL Id: hal-02302455

https://hal.science/hal-02302455v1

Submitted on 19 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Physical parameters and ±0.2% parallax of the detached eclipsing binary V923 Scorpii*

T. Pribulla<fnref n="1" />,**, A. Mérand<fnref n="2" />, P. Kervella<fnref n="3" />,<fnref n="4" />, C. Cameron<fnref n="5" />,<fnref n="6" />, C. Deen<fnref n="7" />,<fnref n="8" />, P. J. V. Garcia<fnref n="9" />,<fnref n="10" />, M. Horrobin<fnref n="11" />, J. M. Matthews<fnref n="12" />, A. F. J. Moffat<fnref n="13" />, O. Pfuhl<fnref n="8" />, S. M. Rucinski<fnref n="14" />, O. Straub<fnref n="4" />, and W. W. Weiss<fnref n="15" />

Received 22 February 2017 / Accepted 18 April 2018

## ABSTRACT

*Context.* V923 Sco is a bright (V = 5.91), nearby (π = 15.46 ± 0.40 mas) southern eclipsing binary. Because both components are slow rotators, the minimum masses of the components are known with 0.2% precision from spectroscopy. The system seems ideal for very precise mass, radius, and luminosity determinations and, owing to its proximity and long orbital period (~34.8 days), promises to be resolved with long-baseline interferometry.

*Aims.* The principal aim is very accurate determinations of absolute stellar parameters for both components of the eclipsing binary and a model-independent determination of the distance.

*Methods.* New high-precision photometry of both eclipses of V923 Sco with the MOST satellite was obtained. The system was spatially resolved with the VLTI AMBER, PIONIER, and GRAVITY instruments at nine epochs. Combining the projected size of the spectroscopic orbit (in km) and visual orbit (in mas) the distance to the system is derived. Simultaneous analysis of photometric, spectroscopic, and interferometric data was performed to obtain a robust determination of the absolute parameters.

*Results.* Very precise absolute parameters of the components were derived in spite of the parameter correlations. The primary component is found to be overluminous for its mass. Combining spectroscopic and interferometric observations enabled us to determine the distance to V923 Sco with better than 0.2% precision, which provides a stringent test of Gaia parallaxes.

*Conclusions.* It is shown that combining spectroscopic and interferometric observations of nearby eclipsing binaries can lead to extremely accurate parallaxes and stellar parameters.

**Key words.** binaries: spectroscopic – binaries: eclipsing – stars: fundamental parameters – parallaxes – stars: distances

## 1. Introduction

HR6327 (V = 5.91) was found to be a SB2 system by Bennet et al. (1963). The authors determined an orbital period

of P = 34.8189 days and estimated the spectral types of its components as F3IV-V and F3V. Using HR6327 as a comparison star, Bolton & Herbst (1976) serendipitously discovered a primary eclipse in the system of depth 0.35 mag. Kholopov et al. (1981) included the system in the GCVS as V923 Sco.

Fekel et al. (2011, hereafter F11) presented new high-dispersion CCD spectroscopic observations. Fifty-six radial-velocity measurements for both components led to the masses of the components; these measurements have about 0.2% precision and result in an orbital eccentricity of *e* = 0.472. The continuum flux ratio of the secondary to primary component was found to be 0.754 (Δmag = 0.31) at 6430 Å. The spectra showed that the

<footnote n="1">¹ Astronomical Institute, Slovak Academy of Sciences, 059 60 Tatranská Lomnica, Slovakia
e-mail: pribulla@ta3.sk</footnote>

<footnote n="2">² European Southern Observatory, Alonso de Córdova 3107, Casilla 19001, Santiago 19, Chile</footnote>

<footnote n="3">³ Unidad Mixta Internacional Franco-Chilena de Astronomía (CNRS UMI 3386), Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile</footnote>

<footnote n="4">⁴ LESIA (UMR 8109), Observatoire de Paris, PSL Research University, CNRS, UPMC, Univ. Paris-Diderot, 5 Place Jules Janssen, 92195 Meudon, France</footnote>

<footnote n="5">⁵ Department of Mathematics, Physics and Geology, Cape Breton University, 1250 Grand Lake Road, B1P 6L2 Sydney, Nova Scotia, Canada</footnote>

<footnote n="6">⁶ Canadian Coast Guard College, Department of Arts, Sciences, and Languages, Sydney, Nova Scotia B1R 2J6, Canada</footnote>

<footnote n="7">⁷ Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="8">⁸ Max Planck Institute for Extraterrestrial Physics, Giessenbachstr., 85748 Garching, Germany</footnote>

<footnote n="9">⁹ Faculdade de Engenharia, Universidade do Porto, rua Dr. Roberto Frias, 4200-465 Porto, Portugal</footnote>

<footnote n="10">¹⁰ CENTRA – Centro de Astrofísica e Gravitação, IST, Universidade de Lisboa, 1049-001 Lisboa, Portugal</footnote>

<footnote n="11">¹¹ Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany</footnote>

<footnote n="12">¹² Department of Physics and Astronomy, University of British Columbia, Vancouver, BC V6T1Z1, Canada</footnote>

<footnote n="13">¹³ Département de physique, Université de Montréal CP 6128, Succursale Centre-Ville, Montréal, QC H3C 3J7, Canada</footnote>

<footnote n="14">¹⁴ Department of Astronomy and Astrophysics, University of Toronto, 50 St George Street, Toronto, ON M5S 3H4, Canada</footnote>

<footnote n="15">¹⁵ University of Vienna, Institute for Astronomy, Türkenschanzstrasse 17, 1180 Vienna, Austria</footnote>

<footnote>* Based on observations made with ESO telescopes at the Paranal Observatory, under ESO program 091.D-0207 and data from the MOST satellite, a former Canadian Space Agency mission, jointly operated by Microsatellite Systems Canada Inc. (MSCI; formerly Dynacon Inc.), the University of Toronto Institute for Aerospace Studies, and the University of British Columbia; assistance was provided by the University of Vienna.</footnote>

<footnote>** ESO visiting scientist.</footnote>
components are slow rotators; the $v_1 \sin i = 5.2 \pm 1.0$ km s$^{-1}$ and $v_2 \sin i = 8.0 \pm 1.0$ km s$^{-1}$. A pseudosynchronous rotational velocity was estimated as $7$ km s$^{-1}$ for both components. The authors also indicated the system to be suitable for long-baseline interferometric observations.

The search for a possible secondary eclipse in the photometry of Bolton & Herbst (1976) at the predicted spectroscopic conjunction (at phase 0.3713 from the primary eclipse) was not conclusive owing to insufficient precision and quantity of the data. The HIPPARCOS photometry (HIP 83491) does not cover eclipses in the system and does not show any out-of-eclipse variability (95th and 5th brightness percentiles differ by only 0.02 mag). Unfortunately, the star is brighter than the magnitude range of ASAS<fnref n="1" />, $V = 8-14$. Searches in other all-sky databases yielded no useful photometry.

With very precise minimum component masses, the projected size of the orbit from spectroscopy, and the possibility of resolving the visual orbit with long-baseline interferometry the system promise not only precise absolute parameters of the components but also a model-independent determination of the distance.

The present paper is organized as follows: Sect. 2 describes high-precision satellite photometry and VLTI interferometry of the system, Sect. 3 gives a preliminary analysis of individual datasets, and Sect. 4 describes the simultaneous analysis of all datasets. The absolute parameters of the components and their evolutionary stage are discussed in Sects. 5 and 6, respectively.

## 2. New observations and data reduction

Having very precise minimum masses for the components from F11, $M_1 \sin^3 i = 1.4708(31) M_\odot$ and $M_2 \sin^3 i = 1.4178(23) M_\odot$, but only limited photometry (Bolton & Herbst 1976) the system was subject to new photometric and interferometric observations.

### 2.1. MOST photometry

The Microvariability and Oscillations of STars (MOST) microsatellite houses a 15-cm telescope, which feeds a CCD photometer through a single, broadband optical filter (350–700 nm). The initial post-launch performance was described by Matthews et al. (2004). Although the original mission goals were asteroseismology of bright ($V < 6$) solar-type stars, Wolf–Rayet and magnetic pulsating stars, MOST obtained photometry of all kinds of variable stars, transiting exoplanets but also of a few dozens of eclipsing binary stars (see Pribulla et al. 2010).

Unfortunately, V923 Sco is outside the so-called continuous viewing zone; thus, it was alternated with other objects to charge MOST batteries using built-in solar panels (see Walker et al. 2003). Hence, the photometry is noncontinuous and the data segments repeat with the orbital period of the satellite (101.4 min). The new observations were focused on the eclipses in the system using the spectroscopic prediction of F11.

The MOST frames were reduced in the following steps: (1) clipping two-sigma outliers in counts and telescope pointing position, which was repeated twice; (2) removing the power at all frequencies less than 3 c/d that have $S/N > 10$; (3) filtering of sky background and inter-pixel variations using polynomials; (4) filtering of orbital modulation of scattered light by removing patterns averaged across 35 bins phased at the MOST satellite

orbital period; (5) clipping three-sigma outliers in data averaged in 50-point bins; and (6) restoring power removed in step 2.

In 2012, the photometry was obtained from May 12 to May 18 (the secondary eclipse) and from May 31 to June 12 (the primary). Unfortunately, the 2012 data suffer from a saturation of the CCD at maximum light. Only part of the primary minimum is below the saturation level and could be used. The observations were therefore repeated in 2014 on May 1 and 2 (the primary eclipse) and July 23 and 24 (the secondary). The out-of-eclipse light curve shows about 1.8 mmag scatter (as determined from the best fit). The photometry covers both minima. Because of the eccentric orbit the primary minimum is about 0.38 mag deep and partial, while the secondary is a grazing eclipse with only $\Delta$mag = 0.013. The presence of the secondary eclipse is, however, crucial for reliable element determination.

### 2.2. Long-baseline interferometry

Using the projected size of the orbit from F11, $a \sin i = 0.2971(2)$ au, and the HIPPARCOS parallax $\pi = 15.46(40)$ mas (van Leeuwen 2007)<fnref n="2" />, the expected apparent size of the orbit is $a \sin i = 4.59$ mas. The angular sizes of the components assuming $R_1 = 2.0 R_\odot$ and $R_2 = 1.9 R_\odot$ (F11) are expected to be only 0.29 and 0.27 mas.

We observed V923 Sco over 16 epochs spread between April 2013 and June 2014, using the Very Large Telescope Interferometer (Berger et al. 2010; Mérand et al. 2014) equipped with the three-telescope beam combiner AMBER (Petrov et al. 2007), or the four-telescope beam combiner PIONIER (Le Bouquin et al. 2011). The 1.8 m auxiliary telescopes were positioned on the A1-G1-J3 and D0-H0-G1 baseline triplets for the AMBER observations, and the A1-G1-K0-J3 quadruplet for the single PIONIER observing epoch. These configurations provide ground baselines between 56 and 140 m, which are suitable for the resolution of the angular separation of the two components of V923 Sco. We also observed V923 Sco using the GRAVITY instrument (Gravity Collaboration et al. 2017) of the VLTI in June and July 2017. The operational wavelength band of GRAVITY is the K band. The science combiner includes a dispersive element that provides wavelength-dependent measurements of the interferometric quantities (e.g., fringe visibilities and phases) at spectral resolutions of $R = 20, 500, \text{and } 4500$. For our observations of V923 Sco, GRAVITY was used in medium spectral resolution ($R = 500$) and in single field mode, that is, its two beam combiners were recording fringes on V923 Sco.

The system was well resolved during nine epochs. The disks of the components were unresolved. Each observation consisted of three pointings: calibrator, V923 Sco, and calibrator. Two stars were used as calibrators: HD153368 ($\theta_H = 1.005 \pm 0.014$ mas, $\theta_K = 1.009 \pm 0.014$ mas<fnref n="3" />; Mérand et al. 2005, the AMBER and the GRAVITY observations), and HD159941 ($\theta_H = 1.081 \pm 0.015$ mas, $\theta_K = 1.089 \pm 0.015$ mas; Mérand et al. 2005, the AMBER and the PIONIER observations). Each pointing took 20 min, leading to a total of 1 hour per observation. The list of observations with the corresponding configurations is presented in Table 1.

The AMBER raw data have been reduced using `amdlib v3⁴`. The PIONIER raw data were reduced using `pndrs<fnref n="5" />`. The GRAVITY data were reduced using the standard

<footnote n="1">¹ All Sky Automated Survey, www.astrouw.edu.pl/</footnote>

<footnote n="2">² Original HIPPARCOS astrometry solution resulted in $\pi = 15.61(80)$ mas.</footnote>

<footnote n="3">³ $\theta_H$, $\theta_K$ - uniform disk angular diameter in the H and K passbands.</footnote>

<footnote n="4">⁴ http://www.jmmc.fr/data_processing_amber.htm</footnote>

<footnote n="5">⁵ http://www.jmmc.fr/data_processing_pionier.htm</footnote>
<table><caption>**Table 1.** Relative astrometric position, $\Delta X$ and $\Delta Y$, of the secondary component of V923 Sco determined from the VLTI observations.</caption><thead><tr><th>HJD</th><th>φ</th><th>ΔX<br>[mas]</th><th>ΔY<br>[mas]</th><th>a<br>[mas]</th><th>b<br>[mas]</th><th>PA<br>[deg]</th><th>F<sub>2</sub>/F<sub>1</sub></th><th>Bands</th><th>Stations</th><th>Combiner</th></tr></thead><tbody><tr><td>2 456 388.877</td><td>0.758</td><td>-4.251</td><td>+1.850</td><td>0.2633</td><td>0.1250</td><td>113.5</td><td>0.837(11)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 456 507.582</td><td>0.165</td><td>+2.635</td><td>-1.339</td><td>0.1936</td><td>0.1250</td><td>116.9</td><td>0.829(4)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 456 540.584</td><td>0.112</td><td>+3.206</td><td>-1.548</td><td>0.1132</td><td>0.0125</td><td>115.8</td><td>0.9(3)</td><td>H</td><td>A1-G1-K0-J3</td><td>PIONIER</td></tr><tr><td>2 456 727.837</td><td>0.487</td><td>-3.187</td><td>+1.114</td><td>0.4082</td><td>0.1250</td><td>109.3</td><td>0.789(7)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 456 728.852</td><td>0.516</td><td>-2.849</td><td>+0.894</td><td>0.1947</td><td>0.1250</td><td>107.4</td><td>0.899(20)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 456 740.864</td><td>0.861</td><td>-3.550</td><td>+1.587</td><td>0.2311</td><td>0.1250</td><td>114.1</td><td>0.840(11)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 456 817.776</td><td>0.069</td><td>+2.618</td><td>-1.193</td><td>0.1905</td><td>0.1250</td><td>114.5</td><td>0.811(14)</td><td>H + K</td><td>A1-G1-J3</td><td>AMBER</td></tr><tr><td>2 457 932.627</td><td>0.089</td><td>+2.853</td><td>-1.306</td><td>0.0300</td><td>0.0240</td><td>24.6</td><td>0.820(10)</td><td>K</td><td>A0-G1-J2-J3</td><td>GRAVITY</td></tr><tr><td>2 457 936.587</td><td>0.203</td><td>+2.572</td><td>-1.327</td><td>0.0128</td><td>0.0125</td><td>124.1</td><td>0.814(8)</td><td>K</td><td>A0-G1-J2-K0</td><td>GRAVITY</td></tr></tbody></table>

**Notes.** The orbital phases ($\phi$) correspond to the ephemeris HJD $2 456 779.83651 + 34.838646 \times E$ defined by the primary minimum from MOST photometry in 2014 and the orbital period from spectroscopy (F11). The table also gives the flux ratios $F_2/F_1$ of the components (the error of the last digit is given in parenthesis) and the corresponding photometric band(s). The flux ratio is the average for the photometric bands used. The uncertainty of the position is expressed by the error ellipse with major and minor axes $a, b$ and the position angle measured from north through east.

GRAVITY Data Reduction System version 1.0 (Lapeyrere et al. 2014). The binary separations were computed from the reduced files for each epoch using CANDID (Gallenne et al. 2015). Because the closure phase was always achieved (three baselines for the AMBER and six for the PIONIER instrument) giving the measure and direction of the object asymmetry, the secondary position is known without the ±180 degrees ambiguity affecting the case of two baselines or speckle interferometry.

As an illustration, Fig. 1 shows the result of the CANDID adjustment of a binary star model (red curves) to the interferometric observables produced by the GRAVITY instrument for the observation recorded on 2 July 2017.

The errors of the relative positions derived from interferometric observations strongly depend on the direction on the sky. This results from insufficient and nonuniform coverage of the aperture and causes the synthesized PSF to be elongated at a certain angle. Hence errors of all positional measurements were represented by ellipses. Their major axis and orientation were determined by bootstrapping of the data. For all epochs 1000 bootstrapping experiments were performed.

An additional source of position uncertainty of the interferometric observations is the precision of the wavelength calibration. The extra errors were estimated at 5%, 2.5%, and 0.1% for the AMBER, PIONIER, and GRAVITY beam combiners, respectively. The wavelength uncertainty was taken into account in the error ellipses.

The relative positions of the secondary component and the flux ratio ($F_2/F_1$) are listed in Table 1. The average flux ratio in the $H + K$ photometric band, $F_2/F_1 = 0.82 \pm 0.02$, is larger than the spectroscopic flux ratio. This means that the secondary component is redder (=cooler) than the primary.

## 2.3. Observed colors and surface temperatures

According to the 2MASS catalog (Skrutskie et al. 2006), $J = 5.160(29)$, $H = 4.981(24)$, and $K_s = 4.895(20)$ for V923 Sco. The observation was performed just after the primary eclipse at phase 0.00970(5) (ephemeris from Table 2). Using the transformation of Bessell & Brett (1988) $K_s$ transforms to $K = 4.934(20)$.

Feinstein & Forte (1974), studying the field of NGC 6281, listed the visual magnitudes and colors of V923 Sco as $V = 5.90$, $B-V = 0.40$, and $U-B = -0.02$ (three photoelectric observations)$^6$ and mentioned that it is a nonmember object. The authors

determined the distance to the cluster at $d = 560 \pm 30$ pc, and interstellar reddening of $E_{(B-V)} = 0.15 \pm 0.02$. Assuming a uniform distribution of the interstellar material toward NGC 6281, and the HIPPARCOS distance ($64.7^{+1.7}_{-1.6}$ pc), leads to $E_{(B-V)} = 0.0173(25)$ for V923 Sco. The error is dominated by the reducing uncertainty toward NGC 6281.

If we adopt $R_K = A_K/E_{B-V} = 0.346$ and $R_V = A_V/E_{B-V} = 3.09$ (see Rieke & Lebofsky 1985), and assume uncertainty of the observed $V$ magnitude as 0.01 mag, then the reduced apparent brightness of the system is $V = 5.846(12)$ and $K = 4.928(20)$ and the corresponding color index $(V - K)_0 = 0.918(23)$.

The authors of F11 found that the components essentially have solar iron abundances. Their determination of masses and radii of the components corresponds to $\log g_1 = 4.00$ (cgs) and $\log g_2 = 4.03$ (cgs). Using the calibration of Worthey & Lee (2011) for solar iron abundance and $\log g = 4.00$ and $(V - K)_0 = 0.918(23)$ we get $T_{\text{eff}} = 6820 \pm 40$ K. Using dereddened $(B-V)_0 = 0.383$ leads to 6750 K. The $(B-V)_0$ error can be estimated at about 0.02, which propagates to about an 80 K error of $T_{\text{eff}}$. The effective temperature corresponds to the combined color of the components.

# 3. Separate analysis of datasets

The published radial velocities, the MOST light curve and the interferometric visual orbit were first analyzed separately to obtain preliminary parameters and reasonable error estimates for the data independent of formal errors of individual observations.

## 3.1. Published spectroscopy

In multi-dataset modeling of the system, the old photographic spectroscopy of Bennet et al. (1963) was not included because of the very large scatter. Hence, only new CCD spectroscopy of F11 was taken into account.

The authors do not list the errors of individual radial-velocity measurements, but comparing variances of the separate solutions for the primary and secondary component, adopt relative weights $w_1 = 1.0$ and $w_2 = 0.3$, respectively. In our analysis these weights were transformed ($w = 1/\sigma^2$) to data uncertainties $\sigma_1 = 1.00$ km s$^{-1}$ and $\sigma_2 = 1.82$ km s$^{-1}$. The differential-correction optimization (i.e., the steepest-descent method in Press et al. 1986) of the radial velocities leads to practically the same spectroscopic elements as obtained by F11. The resulting $\chi^2 = 2.4$ for (112–7) degrees of freedom (d.o.f.) indicates

<footnote>$^6$ Times of observations not given but are consistent with out-of-eclipse magnitudes of Bolton & Herbst (1976).</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Closure phases (**left column**) and squared visibilities (**right column**) of V923 Sco measured with GRAVITY on 2 July 2017 (HJD=2457936.587). The best fit CANDID binary star model is represented with red curves in each panel. The recorded photometric spectrum is also shown in the *upper left panel* before and after the correction of the telluric absorption lines.</figcaption></figure>

that the true data uncertainties are 6.62 times smaller, that is, $\sigma_1 = 0.151$ km s$^{-1}$ and $\sigma_2 = 0.275$ km s$^{-1}$. The parameter errors, checked by the bootstrap and Monte Carlo techniques, are almost identical to those listed by F11. The minimum masses, $M_1 \sin^3 i = 1.4708(31) M_\odot$ and $M_2 \sin^3 i = 1.4178(23) M_\odot$ listed in their Table 5, obtained from the same spectroscopic elements, however, correspond to the solar-mass parameter, $GM_\odot$, smaller by about 0.23% than adopted by the International Astronomical Union (IAU). Using the IAU resolution B3 for solar and planetary properties (Prša et al. 2016) with $GM_\odot = 1.3271244(10) 10^{20}$ m$^3$ s$^{-2}$, we get the minimum masses as 1.4674(31) $M_\odot$ and 1.4146(23) $M_\odot$. The radial velocities define the orbital period as $P = 34.838648(95)$ days.

## 3.2. MOST light curve

Adopting radius estimates and the projected size of the major axis from F11 and assuming that we see the orbit edge on, the ratios of component radii and the instantaneous separation are about $R_1/a(1-e) = 0.059$ and $R_2/a(1-e) = 0.056$ at periastron, and only $R_1/a(1+e) = 0.021$ and $R_2/a(1+e) = 0.020$ at apastron. This means very small proximity effects (tidal deformation, mutual irradiation, and gravity darkening).

Because of the slow rotation, the ellipsoidal deformation of the components is also small. If $R_{\text{equ}}$ and $R_{\text{pol}}$ are the equato-

rial and polar radii and $\Omega/\Omega_{\text{crit}}$ is the ratio of the observed and critical angular velocity, the flattening of a rotating star can be approximated by (see Rozelot & Neiner 2009)

$$ \frac{R_{\text{equ}} - R_{\text{pol}}}{R_{\text{pol}}} = \frac{1}{2} \left( \frac{\Omega}{\Omega_{\text{crit}}} \right)^2 \qquad (1) $$

with

$$ \Omega_{\text{crit}} = \sqrt{\frac{GM}{R_{\text{equ}}^3}}, \qquad (2) $$

where $G$ is the gravitational constant and $M$ is the mass of the star. For the secondary component of V923 Sco, which rotates faster of the two, $v_2 \sin i = 8$ km s$^{-1}$, the flattening is only $2.27 10^{-4}$.

The reflection effect also has very low amplitude. The fraction of the reflected light can be simply estimated (see Sen 1948) as $(1/4)r_j^2$, where $r_j$ is the fractional radius of the reflecting component (the ratio of its radius and semimajor axis). If both components have the same luminosity, then the amplitude of the reflection effect is $\approx(1/8)r_j^2$. In the case of V923 Sco this gives an amplitude <0.4 mmag, which is below the MOST photometry precision.
<table><caption>**Table 2.** Photometric elements derived from the MOST observations of V923 Sco.</caption><thead><tr><th>Parameter</th><th></th><th>σ</th></tr></thead><tbody><tr><td>P [days]</td><td>34.838579</td><td>0.000008</td></tr><tr><td>T<sub>min</sub> [HJD]</td><td>2 456 779.83629</td><td>0.00012</td></tr><tr><td>i [deg]</td><td>87.687</td><td>0.005</td></tr><tr><td>r<sub>1</sub> + r<sub>2</sub></td><td>0.06182</td><td>0.00014</td></tr><tr><td>r<sub>2</sub>/r<sub>1</sub></td><td>0.8965</td><td>0.0028</td></tr><tr><td>e cos ω</td><td>-0.1828</td><td>0.0004</td></tr><tr><td>e sin ω</td><td>+0.4487</td><td>0.0017</td></tr><tr><td>I<sub>1</sub></td><td>0.22799</td><td>0.00013</td></tr><tr><td>I<sub>2</sub></td><td>0.2139</td><td>0.0013</td></tr><tr><td>χ<sup>2</sup>/d.o.f.</td><td>488/(493-9)</td><td>—</td></tr></tbody></table>

**Notes.** The last line gives $\chi^2$ and number of degrees of freedom (d.o.f.) for the mean error of each individual data point of 1.8 mmag.

Thus for the preliminary parameter estimate it is sufficient to model the system by two spherical, limb-darkened stars revolving in an eccentric orbit. A simple program integrating light from the visible surface on the eclipsed component was used.

Unfortunately, fitting a model light curve to the observations of a detached eclipsing binary showing partial eclipses suffers from strong correlations between the parameters (see Southworth et al. 2007). The phase shift of the secondary minimum defines $e \cos \omega$ very well, while $e \sin \omega$ is less well defined by the ratio of minima duration (see Binnendijk 1960, pages 328–329). One can easily determine the sum of fractional radii of the components, $r_1 + r_2$, but not their ratio (and hence individual values, which are strongly correlated). The correlation stems from the fact that a depth of a partial eclipse depends primarily on the impact factor (product of separation of components and $\cos i$) and the sum of the radii. The inclination angle, on the other hand, is well determined.

To arrive at a meaningful solution, the spectroscopic flux ratio $F_2/F_1 = 0.754$ (at 6430 Å) was used as an additional dataset to constrain the solutions. Its uncertainty was set (rather arbitrarily) to 0.001. The flux ratio of the components was determined using the following formula for the linear limb-darkening approximation (see Gray 2008):

$$ \frac{F_2}{F_1} = \left( \frac{r_2}{r_1} \right)^2 \frac{I_2 (1 - u_2/3)}{I_1 (1 - u_1/3)}, \quad (3) $$

where $r_1, r_2$ are the fractional radii of the components, $I_1, I_2$ are central intensities of the disks, and $u_1, u_2$ are the linear limb-darkening coefficients for the components. The coefficients were interpolated from dedicated tables of Claret et al. (2014) for the MOST satellite photometry as $u_1 = u_2 = 0.6194$ for $V_\xi = 2$ km s$^{-1}$, $T_{eff} = 6700$ K, $\log g = 4.0$, and [M/H] = 0.

In addition to the 2014 MOST data, only the nonsaturated part of the 2012 MOST light curve during the primary minimum has been used. Using the data from both years better defines the orbital period but also the shape of the primary minimum.

The free parameters were the time of the primary minimum $T_{min}$, the inclination angle $i$, the sum of fractional radii $r_1 + r_2$, the ratio of fractional radii $r_2/r_1 < 1$, $e \cos \omega$, $e \sin \omega$, and the central intensities of the disks, $I_1, I_2$. Because no additional component was detected in the spectroscopy, zero third light, $I_3$, was assumed. All data points were given the same weight.

The resulting parameters are listed in Table 2 and corresponding fits in Fig. 2. The separate fractional radii are

<table><caption>**Table 3.** Orbital elements of the visual orbit obtained from the VLTI, P - orbital period, e - eccentricity, i - inclination angle, $T_0$ - time of periastron passage, $\omega$ - longitude of periastron, $\Omega$ - longitude of the ascending node, and a - semimajor axis.</caption><thead><tr><th>Parameter</th><th></th><th>σ</th><th></th><th>σ</th></tr></thead><tbody><tr><td>P [days]</td><td>34.838646</td><td>—</td><td>34.838646</td><td>—</td></tr><tr><td>e</td><td>0.47204</td><td>—</td><td>0.439</td><td>0.014</td></tr><tr><td>i [deg]</td><td>87.51</td><td>0.23</td><td>87.70</td><td>0.26</td></tr><tr><td>T<sub>0</sub> [HJD]</td><td>2454272.1636</td><td>—</td><td>2454272.00</td><td>0.07</td></tr><tr><td>ω [deg]</td><td>292.853</td><td>—</td><td>293.5</td><td>1.3</td></tr><tr><td>Ω [deg]</td><td>114.41</td><td>0.27</td><td>114.43</td><td>0.29</td></tr><tr><td>a [mas]</td><td>4.750</td><td>0.013</td><td>4.59</td><td>0.05</td></tr><tr><td>χ<sup>2</sup>/d.o.f.</td><td>23.25/(18-3)</td><td>—</td><td>18.06/(18-6)</td><td>—</td></tr></tbody></table>

**Notes.** Two solutions are given: (i) for $P, e, T_0, \omega$ and (ii) for orbital period $P$ adopted from the spectroscopy.

$r_1 = 0.03260(9)$ and $r_2 = 0.02922(8)$. The eccentricity is $e=0.4845(15)$ and the longitude of the periastron passage $\omega=112.17(9)$ degrees. The computed flux ratio of the components, $F_2/F_1 = 0.75404$, is close to the input value.

### 3.3. VLTI visual orbit

The interferometric orbit of V923 Sco is the relative orbit of the fainter component around the brighter component of the eclipsing pair. Although the orientation toward the secondary component is determined without the 180-degree ambiguity, the longitude of the ascending angle $\Omega$ cannot be determined without spectroscopic observations. Because we only had the visual orbit, we would have two equally good solutions differing in $\omega$ and $\Omega$ by 180 degrees.

The longitude of periastron in the case of the spectroscopic orbit (112.853 ± 0.074 deg according to F11) is related to the primary component. Thus, we selected the solution with $\omega \sim 293$ degrees because the relative visual orbit is the orbit of the secondary around the primary.

Nine $\Delta X$ (toward east) and $\Delta Y$ (toward north) positions when the system was resolved with the VLTI were used. The uncertainty ellipses giving the positional uncertainty in the direction of the separation vector were used.

The parameters were adjusted using the differential-correction method. Because we only had nine separation vectors, we performed two solutions fixing some parameters from spectroscopy covering 6.5 yr. In the first solution only the orbital period was fixed and not adjusted. In the second solution the orbital period, eccentricity, time of the periastron passage, and longitude of the periastron were adopted.

The corresponding fit to the relative visual orbit is shown in Fig. 3. The value of reduced $\chi^2_r = 1.505$ indicates either slightly underestimated position errors or a deviating point(s). Combining the projected semimajor axis from spectroscopy, $a \sin i = 0.29705(19)$ au with apparent semimajor axis from interferometry (first solution), $a \sin i = 4.747(18)$ mas, we get the parallax of V923 Sco as 15.96(6) mas.

## 4. Simultaneous modeling of the datasets

The individual datasets, which were not obtained simultaneously, were light curve (HJD 2 456 083 and 2 456 863), radial velocities (2 453 092 – 2 455 465), and visual orbit (2 456 389 – 2 457 937). For a well-detached binary with ~34.8-day period we can, however, assume that the apsidal motion is too slow
<figure><img src="image_3.png" /><figcaption>Fig. 2. Best fits to the secondary (left) and primary (right) minimum of V923 Sco. The 2014 data are plotted with filled circles, the nonsaturated 2012 data with open triangles. For clarity the data errors are not shown during the primary minimum (smaller than the symbols). The phases correspond to the optimum ephemeris for the primary minimum in Table 2. The best fit assuming two limb-darkened spheres is plotted with a solid line. The dotted line corresponds to a simultaneous modeling of the light-curve, radial-velocity and visual-orbit data (Sect. 4).</figcaption></figure>

to cause significant changes of orbital elements throughout the whole time range of the observations (about 13.3 yr). In the case of the MOST observations, we used both 2012 (nonsaturated data in the primary minimum) and 2014 observations to define the orbital period better and to set constraints on the position of the secondary minimum.

Because the system is wide and proximity effects negligible (see Sect. 3.2), all observations were modeled assuming that components are two limb-darkened spheres. As discussed in Sect. 3.2, to constrain the ratio of radii, the flux ratio, $F_2/F_1 = 0.754$, from spectroscopy and, $F_2/F_1 = 0.82 \pm 0.02$, from interferometry were included as additional datasets. Since we did not know the uncertainty of the spectroscopic light ratio, we set this value rather arbitrarily to 0.001. The linear limb-darkening coefficient for the MOST passband was adopted from Tables 5 and 13 of Claret et al. (2014) after each iteration<fnref n="7" />. The minima last only about 0.013 (secondary) and 0.011 (primary) of an orbital period. Hence, while we used outside the minima a 1/360 step in phase, we used a phase step that was seven times finer during the minima to synthesize the light curve and photo-centric radial velocity.

The multi-dataset analysis requires proper definition of the merit function to be minimized. The most obvious is to use reduced $\chi^2$, which for the case of combining light-curve and radial-velocity data is written as

$$ \chi_r^2 = \frac{1}{N_{LC} + N_{RV}} \left[ \sum_{i=1}^{N_{LC}} \frac{[y_i(m) - y_i(o)]^2}{\sigma_i^2} + \sum_{i=1}^{N_{RV}} \frac{[z_i(m) - z_i(o)]^2}{\sigma_i^2} \right], \quad (4) $$

where $N_{LC}$, $N_{RV}$ are the numbers of data points for light-curve and radial-velocity data, $y$ are intensities, $z$ are radial velocities, while “m” denotes model and “o” observed. If the number of data points is markedly different between the datasets and there is a discrepancy in optimal parameters, the optimization gives higher weight to the dataset with the larger number of points. A solution to avoid this problem is to give all datasets the same

weight; this can be carried out by modifying Eq. (4) to

$$ \chi_r^2 = \frac{1}{N_{LC}} \sum_{i=1}^{N_{LC}} \frac{[y_i(m) - y_i(o)]^2}{\sigma_i^2} + \frac{1}{N_{RV}} \sum_{i=1}^{N_{RV}} \frac{[z_i(m) - z_i(o)]^2}{\sigma_i^2}. \quad (5) $$

The optimization was performed by the differential correction method. The modified sum of $\chi^2$ (Eq. (5)) was used as the merit function. The resulting parameters for all cases are listed in Table 4. The fit to the MOST light curve corresponding to the combined radial-velocity and light-curve solution is plotted in Fig. 2.

Most parameters are consistent within one or two $\sigma$ between the solutions. Small inconsistencies are visible, for example, in the time of periastron passage. The light and radial-velocity curves better define the orbit orientation ($\omega$) and eccentricity than the visual orbit. The position of the secondary eclipse observed by MOST perfectly defines $e \cos \omega$, while radial-velocity curves define $e \sin \omega$. On the other hand, the parallax of V923 Sco and the longitude of the ascending angle cannot be determined without the positional data. The temperature of the secondary component was determined with a very low formal error. The error estimate is affected by an unknown error of the spectroscopic flux ratio, $F_2/F_1$. Hence, it must be used with caution.

## 5. Absolute parameters and distance

The distance and temperature consistency of the component can be checked using the absolute parameters of the components and the bolometric correction. Because proximity effects can be neglected and the components are very close to a spherical shape (see Sect. 3.2) their luminosity can be satisfactorily derived from Stefan–Boltzmann’s law. The absolute parameters corresponding to the combined light-curve and radial-velocity solution are listed in Table 5. We used the IAU resolution B3 prescriptions for the solar and planetary properties (Prša et al. 2016)<fnref n="8" />.

<footnote n="7">⁷ The distance to the system was also optimized (defined primarily by the apparent size of the visual orbit and the size of the spectroscopic orbit).</footnote>

<footnote n="8">⁸ $R_{\odot} = 695,700$ km, $T_{\text{eff}\odot} = 5,772$ K.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Relative orbit of the secondary component of V923 Sco and its best fit assuming spectroscopic elements (*P*, *e*, *ω*, *T*₀) (top). The nodal line is plotted by dash dots. Character “A” shows the position of the ascending node. Close-up view of error ellipses (dot-dashed line) at 9 epochs are shown at *bottom*. The orbit is prograde (counterclockwise motion on the sky). North is up and east is left on both panels.</figcaption></figure>

Using bolometric corrections interpolated from Table 1 of Popper (1980), B.C. = -0.018 (for $T_1 = 6750$ K) and B.C. = -0.026 (for $T_2 = 6562$ K) and $M_{bol} = 4.74$, we get the absolute visual magnitudes as $M_{V1} = 2.547 \pm 0.026$ and $M_{V2} = 2.873 \pm 0.027$ and the combined absolute brightness $M_V = 1.945 \pm 0.037$. For the observed dereddened visual brightness $V = 5.846(12)$ this gives the distance $d = 60.28 \pm 1.10$ pc or $\pi = 16.59 \pm 0.30$ mas. The combined solution of radial-velocity and visual-orbit data resulted in a model-independent parallax $\pi = 15.975 \pm 0.026$ mas or distance $d = 62.59 \pm 0.09$ pc. To reach a model-independent distance, the components would have to be slightly cooler and thus less luminous. The model-independent distance deviates more than one sigma from the revised HIPPARCOS parallax of $15.46 \pm 0.40$ mas or $d = 64.7_{-1.6}^{+1.7}$ pc but is reasonably consistent with the original HIPPARCOS data reduction ($\pi = 15.61 \pm 0.80$ mas). The discrepancy can, partially,

<table><caption>**Table 4.** Multi-dataset modeling of V923 Sco observables.</caption><thead><tr><th>Parameter</th><th>LC+RV</th><th>LC+VO</th><th>VO+RV</th><th>LC+RV+VO</th></tr></thead><tbody><tr><td><i>P</i> [days]</td><td>34.838595(6)</td><td>34.838646</td><td>34.83870(7)</td><td>34.838593(6)</td></tr><tr><td><i>T</i><sub>0</sub> [HJD]</td><td>0.1716(7)</td><td>0.1736(10)</td><td>0.159(4)</td><td>0.1719(7)</td></tr><tr><td><i>i</i> [deg]</td><td>87.724(4)</td><td>87.727(4)</td><td>87.73(8)</td><td>87.725(4)</td></tr><tr><td><i>e</i></td><td>0.47243(15)</td><td>0.47161(19)</td><td>0.4721(4)</td><td>0.47242(14)</td></tr><tr><td><i>r</i><sub>1</sub></td><td>0.03182(4)</td><td>0.03180(5)</td><td>—</td><td>0.03182(4)</td></tr><tr><td><i>r</i><sub>2</sub></td><td>0.02897(3)</td><td>0.02897(3)</td><td>—</td><td>0.02897(3)</td></tr><tr><td><i>K</i><sub>1</sub> [km s<sup>-1</sup>]</td><td>51.673(20)</td><td>51.65</td><td>51.65(3)</td><td>51.673(20)</td></tr><tr><td><i>K</i><sub>2</sub> [km s<sup>-1</sup>]</td><td>53.597(20)</td><td>53.58</td><td>53.58(3)</td><td>53.597(20)</td></tr><tr><td><i>V</i><sub>0</sub> [km s<sup>-1</sup>]</td><td>-15.079(8)</td><td>—</td><td>-15.082(17)</td><td>-15.084(8)</td></tr><tr><td><i>ω</i> [deg]</td><td>112.964(13)</td><td>113.058(19)</td><td>112.83(7)</td><td>112.966(13)</td></tr><tr><td><i>Ω</i> [deg]</td><td>—</td><td>114.325(28)</td><td>114.34(10)</td><td>114.17(26)</td></tr><tr><td><i>π</i> [mas]</td><td>—</td><td>—</td><td>15.975(26)</td><td>16.018(9)</td></tr><tr><td><i>T</i><sub>eff2</sub> [K]</td><td>6562(3)</td><td>6565(4)</td><td>—</td><td>6565(3)</td></tr><tr><td>d.o.f.</td><td>606–12</td><td>510–9</td><td>128–10</td><td>622–12</td></tr><tr><td><i>χ</i><sub>r</sub><sup>2</sup> (LC)</td><td>1.009</td><td>1.146</td><td>—</td><td>1.005</td></tr><tr><td><i>χ</i><sub>r</sub><sup>2</sup> (RV)</td><td>1.194</td><td>—</td><td>1.064</td><td>1.198</td></tr><tr><td><i>χ</i><sub>r</sub><sup>2</sup> (FR)</td><td>0.042</td><td>0.042</td><td>—</td><td>0.047</td></tr><tr><td><i>χ</i><sub>r</sub><sup>2</sup> (VO)</td><td>—</td><td>2.001</td><td>1.971</td><td>1.990</td></tr></tbody></table>

**Notes.** Fixed parameters are listed without error; parameters irrelevant for a given dataset(s) are skipped. The number of data points is 493 for the MOST light curve (LC), 112 for the radial-velocity (RV) curves (56 epochs), and 18 (9 epochs) for the visual orbit (VO) and two flux ratios (FR). The parameters listed are as follows: *P* indicates the orbital period, *T*₀ the time of periastron passage (listed as HJD-2454272), *i* the inclination angle, *e* the orbital eccentricity, *T*₂ the secondary component temperature, *r*₁, *r*₂ the fractional radii of the components, *K*₁, *K*₂ the semiamplitudes of the radial-velocity changes, *V*₀ the mass-center velocity, *ω* the longitude of periastron passage, *Ω* the longitude of ascending node, *π* the parallax, *T*<sub>eff2</sub> the effective temperature of the secondary component, and d.o.f. the number of degrees of freedom and reduced *χ*<sub>r</sub><sup>2</sup> for individual datasets. The modeling of the MOST light curve included a normalization factor not listed in the table but very close to unity.

<table><caption>**Table 5.** Absolute parameters of the components of V923 Sco from the combined light-curve and radial-velocity curve solution.</caption><thead><tr><th>Parameter</th><th>σ</th><th></th></tr></thead><tbody><tr><td><i>R</i><sub>1</sub> [R<sub>☉</sub>]</td><td>2.0246</td><td>0.0026</td></tr><tr><td><i>M</i><sub>1</sub> [M<sub>☉</sub>]</td><td>1.4714</td><td>0.0014</td></tr><tr><td><i>L</i><sub>1</sub> [L<sub>☉</sub>]</td><td>7.67</td><td>0.18</td></tr><tr><td>log <i>g</i><sub>1</sub> [cgs]</td><td>3.9931</td><td>0.0012</td></tr><tr><td><i>R</i><sub>2</sub> [R<sub>☉</sub>]</td><td>1.8496</td><td>0.0019</td></tr><tr><td><i>M</i><sub>2</sub> [M<sub>☉</sub>]</td><td>1.4186</td><td>0.0013</td></tr><tr><td><i>L</i><sub>2</sub> [L<sub>☉</sub>]</td><td>5.72</td><td>0.14</td></tr><tr><td>log <i>g</i><sub>2</sub> [cgs]</td><td>4.0558</td><td>0.0010</td></tr><tr><td><i>a</i> [R<sub>☉</sub>]</td><td>63.941</td><td>0.014</td></tr></tbody></table>

**Notes.** We assume $T_{eff1} = 6750$ K and $T_{eff2} = 6562 \pm 3$ K. An additional temperature uncertainty resulting from the $V-K_s$ color $\sigma(T) = 40$ K is propagated to errors of the luminosities.

result from the photocenter motion in this high-eccentricity system not accounted for in the HIPPARCOS astrometry reductions. The photocenter semiamplitude can be found as

$$A = \alpha \frac{q - \frac{F_2}{F_1}}{(1 + q)(1 + \frac{F_2}{F_1})}, \qquad (6)$$

where $\alpha$ is the maximum separation, $q = M_2/M_1$ the mass ratio, and $F_2/F_1$ is the flux ratio. If $q = 0.964$, $\alpha = 5.36$ mas, and $F_2/F_1 = 0.754$ (assuming that the flux ratio in the HIPPARCOS
<figure><img src="image_5.png" /><figcaption>Fig. 4. Hertzsprung-Russell diagram for the components of V923 Sco. Left panel: isochrones and evolutionary tracks of Girardi et al. (2000); right panel: those of Demarque et al. (2004). The tracks assuming convective overshooting and solar metallicity were selected in both cases. The temperature uncertainties were set to 200 K. The absolute parameters of the components correspond to the combined solution of MOST light curve and published radial-velocity data.</figcaption></figure>

band was the same as the spectroscopic ratio at 6430 Å), we have A = 0.33 mas. This is not a negligible oscillation compared to the HIPPARCOS data precision.

The distance to V923 Sco can also be checked using the surface brightness-color relations derived from long-baseline interferometry. Relations giving empirical stellar angular diameters for a zero-magnitude star can be found in Boyajian et al. (2014). Using the polynomial coefficients from their Table 1 for $(V-K)_0 = 0.918(23)$ (Sect. 2.3) we get the zero magnitude angular radius as $\theta_{mv}=0 = 5.96 \pm 0.08$ mas. Combining the relation for the angular radius of a star (in mas),

$$ \theta = \frac{2000R}{214.94d'} \quad (7) $$

where R is the radius of a star in solar radii and d is the distance in pc, and using the obvious equation $V = -5 \log(\theta/\theta_0)$ we get for the magnitude of a star

$$ V = -4.8436 - 5 \log R + 5 \log d + 5 \log \theta_{mv=0}. \quad (8) $$

Assuming $d = 62.59$ pc (our combined radial-velocity and visual-orbit solution) we have $V_1 = 6.53(3)$ and $V_2 = 6.73(3)$ or $V = 5.87(3)$. This is consistent with the dereddened brightness of V923 Sco, $V = 5.846(12)$.

## 6. Evolutionary state

The evolutionary state of the components was studied by F11 using the Padova solar-abundance tracks (G2000, Girardi et al. 2000). Their Fig. 3 shows that the primary component of V923 Sco is overluminous for its mass. Its luminosity corresponds to about 5% larger theoretical mass. The secondary position was consistent with theory.

Because of the high inclination angle, $i \sim 87.7$ degrees, the true masses of the components are only 0.24% higher than the minimum masses of F11. Although the radii of the components of V923 Sco are now known with higher precision, the luminosity is more uncertain because of the use of color-temperature calibrations, which contain about 40 K error (Kervella et al. 2004).

The radii of the components also depend on the flux ratio. Having both the continuum flux ratio from the spectroscopy of F11 and the $H+K$ band flux ratio from new interferometry puts a strong constraint on the modeling.

While the luminosities and temperatures of the components are encumbered with the error related to color-temperature relations, their masses and radii are now known with 0.1-0.2% errors. This is enough to perform tests of stellar evolution models (see, e.g. Torres et al. 2010). V923 Sco was studied using (i) the G2000 and (ii) Yonsei-Yale ($Y^2$, Demarque et al. 2004) evolutionary tracks and isochrones. In both cases, models corresponding to solar chemical abundance were used. Both libraries use similar physics but the $Y^2$ models differ in a number of details such as the treatment of convective core overshooting, and helium- and heavy-element diffusion. The position of the components of V923 Sco in the $T_{eff} - \log g$ plane for both model libraries is shown in Fig. 4. The figure shows the primary component of V923 Sco to be overluminous for its mass for both model libraries, while the parameters of the secondary component are consistent with the theory.

The masses and radii of the components of V923 Sco are plotted together with data of other detached eclipsing binaries listed in the online version<fnref n="9" /> of DEBCat (Southworth 2015) in Fig. 5. The figure shows rather good agreement of the observed masses and radii for both components of V923 Sco with the theoretical isochrone for 2 Ga with about 0.1 Ga uncertainty.

## 7. Discussion and conclusions

Eclipsing binaries are crucial objects for astrophysics. Combining photometric and spectroscopic observations enables us to determine absolute parameters of stars that are necessary to test models of stellar evolution. Adding positional measurements and flux ratios (determined from spectroscopy and interferometry) gives an additional constraint and makes the results more robust.

We obtained high-precision photometry of both eclipses of V923 Sco with the MOST satellite and the spatially resolved

<footnote n="9">⁹ http://www.astro.keele.ac.uk/~jkt/debcat/</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Mass-radius diagram for V923 Sco (large symbols). The measurements are compared to theoretical isochrones of Girardi et al. (2000) for 100 Ma to 2 Ga for solar metallicity ($Y = 0.70, Z = 0.02$). The tracks assuming convective overshooting were selected. The vectors connect the positions of secondary to primary components in other well-studied detached binaries in DEBcat of Southworth (2015).</figcaption></figure>

orbit with the VLTI at nine epochs. Combining the projected size of the spectroscopic orbit, $a \sin i$, and that from the visual orbit we derive the distance to the system. Simultaneous analysis of photometric, spectroscopic, and interferometric data leads to the first reliable determination of the absolute parameters.

The main results of the combined analysis of published spectroscopic data, new satellite photometry, and VLTI interferometry of V923 Sco are as follows:

- First detection of the secondary eclipse of V923 Sco at phase $\phi_{II} = 0.370457(19)$, which defines the orientation of the orbit very well and sets a constraint on the photometric elements.

- Model-independent determination of the distance to the system, $d = 62.59 \pm 0.09$ pc, found by comparing $a \sin i$ from interferometric visual orbit and the projected major axis determined from spectroscopy. The error takes into account wavelength uncertainties of the interferometric combiners used.

- Model-dependent determination of the distance based on the physical parameters and apparent brightness of the system, $d = 60.28 \pm 1.10$ pc or $\pi = 16.59 \pm 0.30$ mas. This indicates that the components are slightly cooler than those found by spectroscopy.

- Determination of the fundamental parameters of the two components with superb accuracy (see Table 5).

- The mass-radius diagram for the system indicates that the age of the system is close to 2 Ga.

The formal precision of the model-independent parallax, 26 µas, is several times worse than can be expected for single objects in the bright-star regime of the Gaia spacecraft (see de Bruijne 2012). Still, V923 Sco is an important benchmark system to test the Gaia parallaxes (see Stassun & Torres 2016).

With reliable component radii the system can also be used to calibrate surface brightness-color relations (see Graczyk et al. 2017).

Although the combined solution resulted in very accurate component parameters and orbital elements, the solution

strongly depends on the spectroscopic flux ratio $F_2/F_1 = 0.754$ found by F11. The modeling, however, shows that the interferometric estimate in the near-infrared is consistent with the spectroscopic estimate. While the observed masses and radii show good accord with the 2 Ga isochrone, the primary component is overluminous for its mass.

**Acknowledgements.** The authors thank the anonymous referee for valuable comments that improved the manuscript significantly. This research took advantage of the SIMBAD and VIZIER databases at the CDS, Strasbourg (France), and NASA's Astrophysics Data System Bibliographic Services. TP acknowledges support of the Fizeau exchange program (funding from the European Union's FP7 research and innovation program under Grant Agreement 312430, OPTI-CON) and ESO visiting scientist program. This work has been supported by project VEGA 2/0031/18 and the Slovak Research and Development Agency under the contract No. APVV-15-0458. AFJM, JMM, and SMR are grateful for financial aid to NSERC (Canada). AFJM is also grateful for financial assistance to FQRNT (Quebec).

**Note added in proof.** The parallax of V923 Sco published in the second Gaia data release issued on April 25, 2018 is $\pi = 16.070 \pm 0.054$ mas (Gaia Collaboration et al. 2018). The remarkable agreement with our independent orbital parallax $\pi = 16.018 \pm 0.009$ mas ($\Delta\pi = 0.9$ sigma) is a strong indication of the reliability of Gaia DR2 parallaxes, even close to the brightness limit of the satellite (G(V923 Sco) = 5.792).

## References

Bennet, N. W. W., Evans, N. S., & Laing, J. D. 1963, R. Obs. Bull., 78

Berger, J.-P., Zins, G., Lazareff, B., et al. 2010, Proc. SPIE, 7734, 773435

Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134

Binnendijk, L.A. 1960, Properties of double stars (Philadelphia: University of Pennsylvania Press)

Bolton, C. T., & Herbst, W. 1976, AJ, 81, 339

Boyajian, T. S., van Belle, G., & von Braun, K. 2014, AJ, 147, 47

Claret, A., Dragomir, D., & Matthews, J. M. 2014, A&A, 567, A3

de Bruijne, J. H. J. 2012, Ap&SS, 341, 31

Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, 155, 667

Feinstein, A., & Forte, J. C. 1974, PASP, 86, 284

Fekel, F. C., Williamson, M. H., & Henry, G. W. 2011, AJ, 141, 145

Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

Gallenne, A., Mérand, A., Kervella, P., et al. 2015, A&A, 579, A68

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371

Graczyk, D., Konorski, P., Pietrzyński, G., et al. 2017, ApJ, 837, 7

Gravity Collaboration (Abuter, R., Accardo, M., et al.) 2017, A&A, 602, A94

Gray, D. F. 2008, The Observation and Analysis of Stellar Photospheres (Cambridge: Cambridge University Press)

Kervella, P., Thévenin, F., Di Folco, E., & Ségransan, D. 2004, A&A, 426, 297

Kervella, P., Bigot, L., Gallenne, A., & Thévenin, F. 2017, A&A, 597, A137

Kholopov, P. N., Samus, N. N., Kukarkina, N. P., Medvedeva, G. I., & Perova, N. B. 1981, Information Bulletin on Variable Stars, 1921

Lapeyrere, V., Kervella, P., Lacour, S., et al. 2014, Proc. SPIE, 9146, 91462D

Le Bouquin, J.-B., Berger, J.-P., Lazareff, B., et al. 2011, A&A, 535, A67

Matthews, J. M., Kuschnig, R., Guenther, D. B., et al. 2004, Nature, 430, 51

Mérand, A., Bordé, P., & Coudé du Foresto, V. 2005, A&A, 433, 1155

Mérand, A., Abuter, R., Aller-Carpentier, E., et al. 2014, Proc. SPIE, 9146, 91460J

Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1

Popper, D. M. 1980, ARA&A, 18, 115

Press, W. H., Flannery, B. P., & Teukolsky, S. A. 1986, Numerical recipes (Cambridge: Cambridge University Press)

Pribulla, T., Rucinski, S. M., Latham, D. W., et al. 2010, Astron. Nachr., 331, 397

Prša, A., Harmanec, P., Torres, G., et al. 2016, AJ, 152, 41

Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618

Rozelot, J.-P., & Neiner, C. 2009, Lect. Notes Phys., 765

Sen, H. K. 1948, Proc. Nat. Acad. Sci., 34, 311

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Southworth, J. 2015, ASP Conf. Ser., 496, 164

Southworth, J., Bruntt, H., & Buzasi, D. L. 2007, A&A, 467, 1215

Stassun, K. G., & Torres, G. 2016, ApJ, 831, L6

Torres, G., Andersen, J., & Giménez, A. 2010, A&ARv, 18, 67

van Leeuwen, F. 2007, A&A, 474, 653

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023

Worthey, G., & Lee, H.-c. 2011, ApJS, 193, 1