← Latest papers

🔭 astrophysics

Rediscussion of Eclipsing Binaries. Paper XXX. The Slightly Evolved F-type System BK Pegasi

This paper presents high-precision measurements of the masses and radii of the slightly evolved F-type detached eclipsing binary BK Pegasi using TESS light curves and spectroscopic data, confirming that these properties align with theoretical stellar evolutionary models for an age of 2.65 Gyr while also providing an updated ephemeris for the system.

Original authors: A. C. Kutluay, J. Southworth

Published 2026-03-24

📖 5 min read🧠 Deep dive

Original authors: A. C. Kutluay, J. 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, BK Pegasi, locked in a cosmic waltz around a common center of gravity. They are so close that, from our perspective on Earth, they periodically pass in front of one another, dimming the total light we see. This makes them an eclipsing binary, a natural cosmic laboratory where we can weigh and measure stars with incredible precision.

This paper is like a high-tech "re-check" of this star system. The authors, Ahmet Cem Kutluay and John Southworth, decided to take a fresh look using the latest data from NASA's TESS satellite (a space telescope designed to hunt for planets, but great for watching stars too) combined with old data from previous decades.

Here is the story of their findings, broken down into simple concepts:

1. The Cosmic Scale Model

Usually, measuring how heavy a star is or how big it is requires guessing and using complex math that can go wrong. But with an eclipsing binary, nature does the math for us.

The Analogy: Imagine you are trying to guess the size and weight of two dancers spinning in the dark. If you can't see them, it's hard. But if they are holding hands and spinning in front of a bright spotlight, you can see exactly how much light they block when they pass each other. By timing exactly when they block the light and how much light is lost, you can calculate their exact size and weight.
The Result: Using this method, the authors measured the stars with extreme precision.
- Star A (The Big One): It's about 1.4 times the mass of our Sun and nearly 2 times the Sun's radius.
- Star B (The Small One): It's about 1.25 times the mass of our Sun and 1.46 times the Sun's radius.
- Why it matters: These measurements are accurate to within 0.3%. That's like measuring the distance from London to New York and being off by less than the length of a football field. This is the "gold standard" for measuring stars.

2. The "Older" Star is the Bigger One

There is a twist in this story. Usually, the bigger star is also the hotter, bluer, and more energetic one. But in BK Pegasi, the bigger star (Star A) is actually slightly cooler than its smaller partner.

The Analogy: Think of a balloon. As you blow more air into it, it gets bigger. But if you stop blowing and let it sit, it might start to sag or change shape. Star A is "sagging" because it is getting older. It has used up most of the hydrogen fuel in its core and is starting to puff up into a "sub-giant" (a star that is about to leave its main life stage). Even though it's bigger, it's cooling down as it expands, making it look slightly dimmer than its younger, tighter partner.

3. The Cosmic Clock

The stars orbit each other every 5.49 days. The authors acted like detectives, collecting "timestamps" of when the eclipses happened. They used:

New data from the TESS satellite (taking photos in 2022 and 2024).
Old data from the 1990s and 2000s.
The Analogy: Imagine you are trying to predict when a train will arrive at a station. You have a schedule from 1990, a schedule from 2010, and a live camera feed from today. By comparing all of them, you can see if the train is running exactly on time, or if it's slowly speeding up or slowing down.
The Finding: The train is running very close to schedule, but there is a tiny bit of "jitter" in the timing. The authors suspect this might be due to "star spots" (like sunspots, but on these stars) or a very slow wobble in the orbit (apsidal motion), but they need more data to be sure.

4. The Age and Makeup

The authors compared their precise measurements against computer models of how stars evolve.

The Analogy: It's like comparing a real car to a blueprint. If the real car has 100,000 miles on it, the blueprint should show a car that looks a bit worn.
The Match: The stars fit a model for a system that is about 2.65 billion years old (roughly half the age of our Solar System).
The Discrepancy: There is a small conflict. The chemical makeup (metallicity) measured by the authors' models suggests the stars have a bit more "heavy elements" than previous studies claimed. It's like the blueprint says the car is made of a specific alloy, but a previous mechanic said it was made of a slightly different one. The authors suggest their new, more precise measurements might correct the older guess.

5. The Distance

Finally, they calculated how far away BK Pegasi is.

The Result: They found it is about 301 light-years away.
The Check: They compared this to data from the Gaia satellite (which measures star positions with laser-like precision). The Gaia data said it was about 306 light-years away. The authors' result is slightly closer, but the two numbers are close enough to agree that they are both very good estimates.

The Bottom Line

This paper is a "quality control" report for the universe. By using modern space telescopes to re-examine old star systems, the authors have sharpened our understanding of how stars grow, age, and behave. They confirmed that BK Pegasi is a pair of slightly evolved stars, measured their sizes and weights with record-breaking accuracy, and provided a better clock for when their eclipses will happen in the future.

It's a reminder that even with advanced technology, the best way to understand the universe is often to look at the same old things with new, sharper eyes.

1. Problem Statement

The study focuses on BK Pegasi (BK Peg), a double-lined detached eclipsing binary (dEB) system consisting of two late-F type stars. While previous studies (Popper, Demircan et al., and Clausen et al. [CL10]) had established basic orbital and physical parameters, discrepancies remained regarding:

Precise Physical Dimensions: Previous mass and radius measurements showed significant variations (e.g., masses ranging from 1.43/1.28 $M_\odot$ to 1.41/1.25 $M_\odot$ ).
Evolutionary Status: Conflicting interpretations existed regarding whether the stars were still on the main sequence or evolving toward the sub-giant stage, and the precise age of the system.
Orbital Dynamics: The system has a slight eccentricity, raising questions about apsidal motion and the stability of the orbital ephemeris.
Data Limitations: Older studies relied on ground-based photometry and photographic spectra, which often lack the precision and continuous coverage required for high-precision astrophysical modeling.

The primary objective was to re-analyze BK Peg using high-precision space-based photometry (TESS) combined with existing spectroscopic data to refine the system's physical properties to a level of precision rarely achieved for normal stars.

2. Methodology

The authors employed a multi-faceted approach combining new photometric data with a synthesis of historical spectroscopic data:

Photometric Data:
- Utilized TESS (Transiting Exoplanet Survey Satellite) data from Sector 57 (October 2022) with 120-second cadence.
- Data was processed using the lightkurve package, applying Simple Aperture Photometry (SAP) and filtering out low-quality points.
- Light curves were manually trimmed to remove gaps and poorly sampled eclipse regions.
- Modeling: The jktebop code was used to model the light curves. Parameters fitted included fractional radii, orbital inclination, eccentricity ( $e$ ), argument of periastron ( $\omega$ ), and limb darkening coefficients (using the power-2 law).
- Uncertainties were estimated using 10,000 Monte Carlo simulations.
Spectroscopic Data:
- Re-analyzed radial velocity (RV) data from two previous studies: Popper (12) (photographic spectra) and CL10 (FIES spectrograph).
- Four distinct RV analysis scenarios were tested: Popper only, CL10 only, combined data with separate systemic velocities ( $V_\gamma$ ), and combined data with a common $V_\gamma$ .
- The jktebop code was used to fit the RV curves, fixing the orbital period and eccentricity parameters derived from the photometry.
Physical Property Derivation:
- The jktabsdim code was used to derive absolute physical properties (mass, radius, density, surface gravity) by combining the photometric fractional radii with the spectroscopic semi-amplitudes ( $K_A, K_B$ ) and effective temperatures ( $T_{\text{eff}}$ ) from CL10.
- Distance was calculated using surface brightness calibrations (Kervella et al.) applied to BV and 2MASS JHKs magnitudes.
Theoretical Comparison:
- Results were compared against PARSEC 1.2 stellar evolutionary models to determine the system's age and metallicity ( $Z$ ).

3. Key Contributions

High-Precision Measurements: The study provides the most precise measurements of mass and radius for BK Peg to date, achieving a precision of 0.3% for both quantities.
Updated Ephemeris: A new, highly precise linear ephemeris was established by combining TESS mid-eclipse times with a comprehensive collection of historical mid-eclipse times from literature and the VarAstro database.
Resolution of Evolutionary Status: The authors definitively place the system at an age of 2.65 Gyr with a metallicity of $Z = 0.017$ (slightly super-solar), resolving previous ambiguities about the evolutionary stage of the primary star.
Data Integration: The paper successfully integrates TESS photometry with older spectroscopic datasets, demonstrating a robust methodology for refining dEB parameters.

4. Key Results

Physical Parameters:

Masses: $M_A = 1.411 \pm 0.004 \, M_\odot$ , $M_B = 1.254 \pm 0.004 \, M_\odot$ .
Radii: $R_A = 1.990 \pm 0.004 \, R_\odot$ , $R_B = 1.460 \pm 0.004 \, R_\odot$ .
Orbital Period: $P = 5.48991130 \pm 0.00000012$ days.
Eccentricity: $e = 0.00368 \pm 0.00004$ (very low, nearly circular).
Effective Temperatures: $T_{\text{eff},A} = 6265 \pm 85$ $T_{eff, A} = 6265 \pm 85$ K, $T_{\text{eff},B} = 6320 \pm 90$ $T_{eff, B} = 6320 \pm 90$ K.
- Note: The primary (Star A) is more massive and larger but slightly cooler than the secondary, a signature of its evolution toward the sub-giant branch.

Distance and Reddening:

Distance: $301.2 \pm 3.6$ pc (derived from surface brightness).
Comparison: This is slightly lower than the Gaia DR3 parallax-derived distance of $306.3 \pm 1.7$ pc.
Reddening: $E(B-V) = 0.04 \pm 0.02$ mag.

Evolutionary Modeling:

Models with solar metallicity ( $Z=0.014$ ) fit the masses/radii but failed to match the temperatures.
Models with $Z = 0.017$ provided the best fit for mass, radius, and temperature simultaneously, suggesting an age of 2650 $\pm$ 50 Myr.
This age and metallicity contradict the lower metallicity ($[Fe/H] = -0.12$) previously reported by CL10, suggesting the need for re-evaluation of spectroscopic abundances.

Orbital Dynamics:

The residuals of the eclipse times show an excess scatter ( $\chi^2_\nu = 4.32$ ), but no clear evidence of apsidal motion or pulsations was found in the TESS data (limit of 0.2 mmag).
The authors note that the argument of periastron ( $\omega$ ) is highly degenerate due to the low eccentricity.

5. Significance

Benchmark for Stellar Models: The extreme precision of the mass and radius measurements (0.3%) makes BK Peg a critical benchmark for testing stellar evolutionary models, particularly for F-type stars in the late main-sequence/sub-giant transition.
Metallicity Tension: The discrepancy between the metallicity required to fit the evolutionary models ( $Z=0.017$ ) and the spectroscopic measurement ( $[Fe/H] \approx -0.12$ ) highlights potential systematic issues in stellar abundance determinations or model physics (e.g., convective core overshooting).
Methodological Validation: The study demonstrates the efficacy of combining TESS photometry with legacy spectroscopic data to refine binary parameters, a strategy applicable to the broader catalog of detached eclipsing binaries (DEBCat).
Future Work: The paper identifies the need for continued monitoring of eclipse times to detect apsidal motion and refine the ephemeris, as well as further spectroscopic analysis to resolve the metallicity conflict.