V1031 Ori: Hierarchical Triple System

• V1031 Ori is a hierarchical triple system comprising a circular Algol-type eclipsing binary and a distant companion in a highly eccentric orbit.
• Precision TESS photometry and spectroscopy yield mass and radius measurements of the inner A-type binary with accuracies near 0.5%–0.7%.
• Analysis of eclipse timing variations and light-curve residuals reveals tidally perturbed δ Scuti p-modes, indicating orbital splitting effects.

V1031 Ori (HD 38735, HR 2001, HIP 27341, MCA 22) is a hierarchical triple system comprising an Algol-type eclipsing binary and a distant visual companion. The inner pair is a circular eclipsing binary of two A-type stars with a period of about $3.4$ days, while the outer subsystem follows a long-period, highly eccentric orbit. Analyses combining TESS photometry, published spectroscopy, interferometry, and eclipse-timing variations have established precise fundamental parameters for the eclipsing pair and identified short-period oscillations interpreted as tidally perturbed $p$-modes from the secondary component (Lee, 2021, Zasche et al., 2014).

1. System configuration and identification

V1031 Ori is described as a hierarchical triple star system with a circular eclipsing inner binary $(\mathrm{AB})$ and a more distant companion $(\mathrm{C})$. In the earlier combined-orbit study it is also characterized as an Algol-type eclipsing binary and a visual binary discovered by McAlister et al. (1983), with a magnitude difference of about $1.5$ mag between the two components (Zasche et al., 2014).

The inner binary period is given approximately as $P \approx 3.4$ days in the triple-system study, while the TESS-based analysis reports $P_{\rm orb} = 3.405592 \pm 0.000018$ d for the eclipsing pair. Previous distance estimates quoted for the system were $215\pm25$ pc and $205\pm36$ pc; the later TESS-based estimate of $210 \pm 26$ pc is stated to be consistent with previous measurements (Lee, 2021, Zasche et al., 2014).

Within the 2014 sample of triple systems, V1031 Ori is singled out by the outer orbit’s combination of long period and very high eccentricity. The outer period is reported as $p$0 yr, and the eccentricity as $p$1, making it the longest and most eccentric outer orbit among the three systems examined in that study (Zasche et al., 2014).

2. Fundamental parameters of the eclipsing pair

The TESS-plus-spectroscopy analysis determined the masses and radii of the eclipsing pair to about $p$2 and $p$3 precision, respectively. The resulting component parameters are as follows (Lee, 2021).

Parameter	Primary	Secondary
Mass ($p$4)	$p$5	$p$6
Radius ($p$7)	$p$8	$p$9
$(\mathrm{AB})$0 (cgs)	$(\mathrm{AB})$1	$(\mathrm{AB})$2
$(\mathrm{AB})$3 (K)	$(\mathrm{AB})$4	$(\mathrm{AB})$5

The mass ratio is reported as $(\mathrm{AB})$6, and the orbital inclination as $(\mathrm{AB})$7 deg. The secondary is therefore both more massive and substantially larger than the primary (Lee, 2021).

The rotational velocities are explicitly described as sub-synchronous. The observed projected rotational velocities are $(\mathrm{AB})$8 km/s and $(\mathrm{AB})$9 km/s, compared with synchronous values $(\mathrm{C})$0 km/s and $(\mathrm{C})$1 km/s. This identifies a measurable departure from synchronous rotation in both stars (Lee, 2021).

3. TESS photometry and light-curve modelling

The photometric analysis used high-precision TESS Sector 6 data with 2-min cadence, comprising 14,871 time-series data points. The observed light curve shows deep eclipses, ellipsoidal modulation, and clear short-period oscillations superimposed on the eclipsing light curve (Lee, 2021).

Binary light-curve modelling was carried out with the Wilson-Devinney (W-D) code. The binary parameters, including the mass ratio and effective temperatures, were fixed; the bolometric albedo and gravity darkening were fixed to $(\mathrm{C})$2 for a radiative envelope; and the system was confirmed to be in a circular orbit. The third-light contribution from the tertiary is reported as

$(\mathrm{C})$3

at quadrature. After subtraction of the binary model, the residuals exhibit clear periodicities (Lee, 2021).

The same study gives the eclipse ephemeris

$(\mathrm{C})$4

where $(\mathrm{C})$5 is the epoch. This provides the timing framework used in the residual analysis and phase-dependent interpretation of the pulsations (Lee, 2021).

4. Oscillation content and frequency extraction

Frequency extraction from the eclipse-subtracted residuals was performed with PERIOD04 using prewhitening and a significance threshold of $(\mathrm{C})$6. The analysis was carried out up to the Nyquist limit of $(\mathrm{C})$7, and only out-of-eclipse portions were used to avoid eclipse-distortion (Lee, 2021).

A total of 23 significant frequencies were identified in two ranges: 16 frequencies in the low-frequency region $(\mathrm{C})$8 and 7 frequencies in the high-frequency region $(\mathrm{C})$9. The low frequencies are described as mostly artefacts or aliases, likely due to imperfect de-trending and binary residual effects, and are considered unlikely to be true $1.5$0-modes for these stars. By contrast, the high frequencies are interpreted as genuine $1.5$1-modes of $1.5$2 Scuti-type pulsation, with pulsation periods $1.5$3 d and pulsation constants $1.5$4 d (Lee, 2021).

The seven quoted $1.5$5-mode-region frequencies are:

Mode	Frequency ($1.5$6)	Amplitude (mmag)
$1.5$7	11.49342	1.81
$1.5$8	12.08130	1.35
$1.5$9	11.19949	0.18
$P \approx 3.4$0	12.66458	0.24
$P \approx 3.4$1	11.78507	0.26
$P \approx 3.4$2	12.37294	0.16
$P \approx 3.4$3	10.91244	0.09

An important interpretive point is that not all detected frequencies are treated symmetrically. The TESS analysis explicitly warns that sixteen frequencies in the $P \approx 3.4$4-mode region may be aliases and artefacts due to imperfect removal of the systematic trends and the binary effects from the TESS data, whereas only the seven frequencies in the $P \approx 3.4$5-mode region are treated as secure oscillation signatures (Lee, 2021).

5. Tidal perturbation, orbital splitting, and oblique pulsation

The high-frequency $P \approx 3.4$6-modes are systematically separated by the orbital frequency, with $P \approx 3.4$7. They are also split from the nearest exact orbital multiple by a nearly constant offset,

$P \approx 3.4$8

with the individual offsets listed for the seven $P \approx 3.4$9-mode frequencies. The paper emphasizes that these modes are not tidally induced in the classical sense of lying at exact orbital harmonics; rather, they are interpreted as tidally perturbed pulsations, i.e.

$P_{\rm orb} = 3.405592 \pm 0.000018$0

or equivalently

$P_{\rm orb} = 3.405592 \pm 0.000018$1

This distinction is central to the physical interpretation (Lee, 2021).

The pulsation amplitudes are clearly modulated with the binary star orbit. The amplitudes are largest during primary eclipse, when the pulsating secondary is in view, and they diminish during secondary eclipse and at quadratures because of partial occultation and changing geometry. The analysis therefore identifies the secondary as the pulsating component and attributes the observed orbital dependence to variation in the visible oscillating surface area (Lee, 2021).

The same phase-dependent behaviour is used to argue that the pulsation axis could be aligned with the tidal axis, defined by the line joining the two stellar centres. By analogy with HD 74423 and CO Cam, the paper interprets the secondary as an oblique pulsator in the context of tidal distortion. A plausible implication is that V1031 Ori belongs to the small set of eclipsing binaries in which tidal geometry directly structures the observable pulsation spectrum rather than merely perturbing stellar structure in a more diffuse way (Lee, 2021).

6. Outer-orbit solution from interferometry and eclipse-timing variation

The 2014 study analyzed V1031 Ori with a combined approach based on interferometry and apparent period variation. The method joins positional measurements of the visual double with an $P_{\rm orb} = 3.405592 \pm 0.000018$2 analysis of eclipse minima timings, interpreting the timing offsets through the Light-Time Effect (LITE). The astrometric orbit was fitted using least-squares and the simplex algorithm, and a simultaneous fit was performed for both data sets with shared orbital parameters. A stated advantage of the technique is that no spectroscopic monitoring of the visual orbit is required, despite the long outer period (Zasche et al., 2014).

The final outer-orbit parameters reported for V1031 Ori are:

Parameter	Value
Outer period $P_{\rm orb} = 3.405592 \pm 0.000018$3 (yr)	$P_{\rm orb} = 3.405592 \pm 0.000018$4
Eccentricity $P_{\rm orb} = 3.405592 \pm 0.000018$5	$P_{\rm orb} = 3.405592 \pm 0.000018$6
Inclination $P_{\rm orb} = 3.405592 \pm 0.000018$7 (deg)	$P_{\rm orb} = 3.405592 \pm 0.000018$8
$P_{\rm orb} = 3.405592 \pm 0.000018$9 (deg)	$215\pm25$0
$215\pm25$1 (deg)	$215\pm25$2
$215\pm25$3 (mas)	$215\pm25$4
$215\pm25$5 ($215\pm25$6)	$215\pm25$7
$215\pm25$8 ($215\pm25$9)	$205\pm36$0
$205\pm36$1 (AU)	$205\pm36$2
$205\pm36$3 (AU)	$205\pm36$4

The quoted LITE expression is

$205\pm36$5

where $205\pm36$6 is the true anomaly, $205\pm36$7 is eccentricity, $205\pm36$8 is the argument of periastron, and $205\pm36$9 is the LITE amplitude. The simultaneous fit minimizes the total chi-square, defined as the sum of the chi-squares from the $210 \pm 26$0 curve and the visual-orbit fit (Zasche et al., 2014).

This solution is described as much improved over earlier, highly uncertain ones: the outer period changed from a previous $210 \pm 26$1 yr estimate to $210 \pm 26$2 yr once observations near periastron became available. The visual orbit is stated to be much better defined because one revolution since discovery is now covered by observations. At the same time, the minima variation remains only mildly constrained because most eclipse timings are far from periastron, and the distance was not independently determined in this analysis because the amplitudes in both the LITE and the visual orbit were not sufficiently well constrained. The paper notes that additional timings, especially as the system approaches future periastron in 2035, together with new radial velocities, would improve the amplitude determination and could potentially enable an independent distance measurement (Zasche et al., 2014).