🔭 astrophysics

Rediscussion of eclipsing binaries. Paper XXIX. The F-type twin system BS Draconis

This paper presents a high-precision analysis of the F-type twin eclipsing binary BS Draconis using 40 sectors of TESS data and spectroscopy, yielding accurate masses, radii, and an age of 1.6 Gyr while identifying the primary eclipse and establishing the system's potential as a celestial clock with an eclipse timing scatter of only 0.37 seconds.

Original authors: John Southworth

Published 2026-03-13

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: John Southworth

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine two stars dancing a perfect, eternal waltz around each other. They are so close, so identical in their steps, and so perfectly matched in size and weight that for decades, astronomers couldn't tell which one was leading and which was following. They were like identical twins wearing the same outfit, spinning in a dark room.

This paper is the story of how John Southworth finally solved the mystery of this cosmic dance, a system called BS Draconis.

Here is the breakdown of the discovery, translated into everyday language:

1. The "Twin" Mystery

For a long time, scientists knew BS Draconis was a pair of stars (a binary system) that passed in front of each other from our point of view, causing them to dim slightly. This is called an eclipsing binary.

Think of it like two identical twins walking past a streetlamp. Every time one walks in front of the other, the light dims. Because the twins were so similar, previous astronomers couldn't tell which twin was slightly bigger or hotter. They were stuck guessing.

2. The Super-Sharp Camera (TESS)

The key to solving this was a space telescope called TESS (Transiting Exoplanet Survey Satellite). Imagine TESS as a camera with a super-sharp lens that never blinks. It watched this star system for a long time, taking 40 separate "movies" (called sectors) of the stars dancing.

Because the data was so clean and high-quality, Southworth could see something tiny that no one else had noticed before: The two eclipses weren't exactly the same depth.

One time the stars crossed, the light dipped just a tiny bit more (by 0.007 magnitudes) than the other time.

The Analogy: Imagine two twins walking past a light. One is wearing a slightly thicker coat. When the "thicker coat" twin blocks the light, the room gets a tiny bit darker than when the other twin blocks it.
The Result: This tiny difference proved that one star (Star A) is slightly hotter, larger, and heavier than the other (Star B). The mystery of "who is who" was finally solved!

3. The Cosmic Scale

With the mystery solved, Southworth could measure the stars with incredible precision.

Mass: Star A weighs about 1.3 times what our Sun weighs. Star B is just a hair lighter.
Size: They are both about 1.4 times the size of our Sun.
Precision: The measurements are so accurate that the error margin is less than the width of a human hair compared to the size of a planet. This is like weighing a car and knowing the weight to within a few grams.

4. The Perfect Clock

Because these stars dance in such a perfect circle and their eclipses are so sharp and deep (like a triangle shape rather than a soft curve), they make an excellent celestial clock.

Southworth calculated exactly when the stars would eclipse each other for years to come. The predictions are so precise that if you set a timer based on these stars, it would only be off by 0.37 seconds over a long period.

Why it matters: Astronomers can use BS Draconis to check if their own clocks (and even the clocks on satellites like TESS) are running fast or slow. It's a "timekeeper" for the universe.

5. The Age and Personality

By comparing these stars to computer models of how stars grow up, the author figured out:

Age: They are about 1.6 billion years old. (For context, our Sun is about 4.6 billion years old, so these stars are in their "young adult" phase).
Personality: They are a bit "metal-poor" (in astronomy, "metals" are elements heavier than hydrogen and helium). They are like a slightly less rich version of our Sun's neighborhood.
Activity: The stars have some "sunspots" and magnetic activity (like solar flares), but they are surprisingly calm. They aren't throwing tantrums; they are just doing their steady dance.

Summary

In short, this paper is about using a high-tech space camera to finally tell apart two identical cosmic twins. By spotting a tiny difference in how they block light, the author measured their exact weight and size, proved one is the "leader" of the pair, and turned the whole system into a super-accurate clock that helps astronomers keep time across the galaxy.

It's a reminder that even when things look identical, there's always a tiny detail waiting to be found if you look closely enough.

1. Problem Statement

BS Draconis (BS Dra) is a detached eclipsing binary (dEB) system consisting of two nearly identical F3 V stars in a 3.36-day circular orbit. Historically, the system has been difficult to analyze definitively because the two components are so similar that previous studies could not determine which star was the "primary" (the hotter, more massive component).

Previous Limitations: Earlier analyses (e.g., Popper & Etzel, Ibanoğlu et al., Milone et al.) struggled with indistinguishable eclipse depths, leading to ambiguous identifications of the primary star and imprecise measurements of physical parameters (masses and radii) due to limited data quality and scattered light curves.
Goal: To resolve the ambiguity regarding the primary component, determine the system's physical properties with high precision, and assess its utility as a celestial clock.

2. Methodology

The study combines high-precision space-based photometry with a re-analysis of historical and modern spectroscopic data.

Photometric Data:
- Utilized 40 sectors of TESS (Transiting Exoplanet Survey Satellite) observations, covering the system extensively within its continuous viewing zone.
- Data processing involved converting flux to differential magnitudes, removing median offsets, and filtering out low-quality data points.
- Modeling: The jktebop code (version 44) was used to model the light curves. The analysis fitted each TESS sector separately to account for variations in third light (background contamination) and instrumental effects.
- Parameters Fitted: Orbital period ( $P$ ), reference time of primary minimum ( $T_0$ ), sum and ratio of fractional radii ( $r_A+r_B$ , $k=r_B/r_A$ ), central surface brightness ratio ( $J$ ), third light ( $L_3$ ), and orbital inclination ( $i$ ). Limb darkening was modeled using the power-2 law.
- Cadence Handling: Special numerical integration was applied to Sector 41 (600s cadence) to correct for lower sampling rates, ensuring consistency with 120s cadence data.
Spectroscopic Data:
- Re-analyzed radial velocity (RV) data from three sources: Fitzgerald (19 photographic spectra), Popper (17 photographic spectra), and Milone et al. (M05, 27 echelle spectra).
- RVs were fitted using jktebop assuming a circular orbit. The study tested models with a single systemic velocity versus separate systemic velocities for each star.
- The M05 dataset was selected as the primary spectroscopic basis due to its higher number of data points and precision compared to the photographic datasets.
Physical Property Derivation:
- Combined photometric radii, orbital inclination, period, and spectroscopic velocity amplitudes ( $K_A, K_B$ ) using the jktabsdim code to calculate masses, densities, and surface gravities.
- Estimated distance using surface brightness calibrations against $BV$ and $JHK_s$ magnitudes, cross-referenced with Gaia DR3 parallax.
- Determined age and metallicity by comparing observed properties against the PARSEC 1.2 evolutionary models.

3. Key Contributions

Definitive Primary Identification: For the first time, the study conclusively identifies the primary eclipse. The primary eclipse is found to be 0.007 mag deeper than the secondary, proving that Star A is the hotter, larger, and more massive component.
High-Precision Parameters: Achieved mass and radius precisions of 1.3% and 0.4%, respectively, significantly improving upon previous literature (which were often limited to ~2% or worse).
Celestial Clock Utility: Demonstrated that BS Dra's deep, V-shaped eclipses allow for timing measurements with an extremely low scatter of 0.37 seconds around a linear ephemeris, making it a potential standard for checking timing accuracy in future missions (e.g., PLATO).
Stellar Activity Analysis: Detected chromospheric activity via Ca II H & K line infilling but confirmed the absence of significant starspot activity in the TESS light curves.

4. Key Results

Physical Properties:

Orbit: Circular, Period $P = 3.364012273$ days.
Star A (Primary):
- Mass: $1.305 \pm 0.015 M_\odot$
- Radius: $1.409 \pm 0.006 R_\odot$
- Effective Temperature: $6630 \pm 150$ K
Star B (Secondary):
- Mass: $1.284 \pm 0.017 M_\odot$
- Radius: $1.400 \pm 0.006 R_\odot$
- Effective Temperature: $6620 \pm 150$ K
Systemic Properties:
- Mass Ratio ( $M_B/M_A$ ): $0.9837 \pm 0.0077$
- Distance: $191.3 \pm 2.9$ pc (consistent with Gaia DR3 parallax of $195.34 \pm 0.50$ pc).
- Age: $1600 \pm 300$ Myr.
- Metallicity: Slightly sub-solar ( $Z \approx 0.014$ ), consistent with previous suggestions of mild metal deficiency.

Ephemeris and Timing:

Ephemeris: $Min I = BJDTDB\ 2459728.011509(2) + 3.364012273(4)E$ .
Timing Precision: The scatter of eclipse midpoints around the linear fit is only 0.37 seconds.
Period Change: No significant period change detected; a $3\sigma$ upper limit of $9.6 \times 10^{-11}$ s s $^{-1}$ was set.

Stellar Activity:

Ca II H & K lines show infilling indicative of chromospheric activity.
No evidence of starspots or pulsations (specifically $\delta$ Scuti or $\gamma$ Doradus) was found in the TESS residuals (limit of $10^{-5}$ mag).

5. Significance

Benchmark for Stellar Models: The high precision of the mass and radius measurements (better than 2%) provides a rigorous test for stellar evolution models, particularly for F-type twin systems. The results support a slightly sub-solar metallicity and an age of ~1.6 Gyr.
Astrometric and Timing Standard: The system's stability and deep eclipses make it an ideal "celestial clock." Its 0.37s timing precision allows it to be used for cross-calibrating timing data from space missions like TESS and the upcoming PLATO mission, helping to identify systematic errors in satellite timing.
Resolution of Historical Ambiguity: The paper resolves decades of uncertainty regarding the identity of the primary component in BS Dra, correcting previous assumptions and providing a definitive physical characterization of the "twin" system.