arxiv: 2604.25056 · v2 · submitted 2026-04-27 · 🌌 astro-ph.SR

Recognition: unknown

V446 Cephei: a β Cep pulsator in a multiple system

A. Moharana , J. Southworth , K. Pavlovski , A. Miszuda , R. S. Rathour , K. G. He{\l}miniak , F. Marcadon , D. M. Bowman

show 2 more authors

T. B. Pawar A. Tkachenko

Authors on Pith no claims yet

Pith reviewed 2026-05-08 03:07 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords beta Cep starseclipsing binarytertiary companionstellar pulsationsmultiple systemsTESS photometryisochrone fittingeclipse timing variations

0 comments

The pith

V446 Cep is a β Cep pulsator in an eclipsing binary with a tertiary companion of minimum mass 4.11 solar masses on a 2303-day orbit, making it either a hierarchical quadruple or a triple with a compact object.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a detailed analysis of the β Cep star V446 Cep in an eclipsing binary using TESS photometry and HERMES spectra. It derives precise stellar parameters for both components, identifies 21 pulsation frequencies, and detects a tertiary companion through eclipse timing variations and pulsation period changes. This leads to the conclusion that the system is either a co-evolving 2+2 quadruple or a triple containing a compact object. Such detailed studies of pulsating stars in multiples provide key constraints on stellar structure and the formation of massive star systems.

Core claim

The eclipsing binary V446 Cep has an orbital period of 3.808567 days, with the β Cep primary having a mass of 10.68 solar masses, radius of 5.864 solar radii, and effective temperature of 24220 K, while the secondary has 1.657 solar masses and 9080 K. Twenty-one pulsation frequencies were extracted, dominated by one at 10.24324 per day. Eclipse timing variations indicate a tertiary with minimum mass 4.11 solar masses on a 2303-day orbit. Using spectral energy distributions and MIST isochrones, the authors conclude V446 Cep is either a co-evolving hierarchical 2+2 quadruple or a triple system where the third body is a compact object.

What carries the argument

Eclipse timing variations and changes in the dominant pulsation frequency to detect and orbit the tertiary companion, combined with isochrone fitting to assess the evolutionary state and multiplicity.

Load-bearing premise

The eclipse timing variations and dominant pulsation period changes are due to a tertiary companion on a 2303-day orbit rather than other effects, and the isochrone fitting accurately determines whether the components are co-evolving.

What would settle it

A failure to detect the predicted 2303-day periodic signal in extended eclipse timing or pulsation data, or direct evidence of a luminous tertiary companion inconsistent with being compact.

Figures

Figures reproduced from arXiv: 2604.25056 by A. Miszuda, A. Moharana, A. Tkachenko, D. M. Bowman, F. Marcadon, J. Southworth, K. G. He{\l}miniak, K. Pavlovski, R. S. Rathour, T. B. Pawar.

**Figure 1.** Figure 1: TESS TPF of V446 Cep with nearby sources. The marker sizes indicate the magnitude difference between V446 Cep and the respective star. Star 2 has a G-band magnitude of 13 compared to 7.8 of V446 Cep. The white lines mark the pixels used in the SPOC aperture. 3 STELLAR PARAMETERS We extracted the stellar, orbital, and atmospheric parameters of the binary system using a combination of photometric and spectro… view at source ↗

**Figure 2.** Figure 2: The top panel shows the binary phase-folded light curve and the binary model from jktebop. The lower panel shows the residuals after subtraction of the binary model. The dominant pulsation (∼ 10 d−1 ) is clearly visible and has around 39 cycles in one orbital phase. The final phoebe2 model is shown in Fig.3 and the MCMC corner plot can be found in Appendix A. The final model gave us a 𝑟1 + 𝑟2 = 0.31468+0.0… view at source ↗

**Figure 3.** Figure 3: Phase-folded binary-LC (in black) with phoebe2 model (in blue). The bottom panel shows the residuals of the fit along with the error bars of the data points, in black. of the parameters are listed in view at source ↗

**Figure 4.** Figure 4: The spectra of the components in V446 Cep as reconstructed by spd. Disentangled spectra are still in a common continuum of a binary system; hence, the highly diluted secondary component is obvious view at source ↗

**Figure 5.** Figure 5: The solution for the spectroscopic orbit from spd is represented in solid lines. The phase distribution of the observed spectra is illustrated with blue (primary) and red (secondary) circles. In spd the orbital parameters are directly optimised by-passing the RVs for individual exposures - therefore symbols are only for illustrative purposes. The Rossiter-McLaughlin effect in the course of the primary and … view at source ↗

**Figure 6.** Figure 6: Disentangled spectrum of the primary component (in red) superimposed on the primary’s spectrum (in black) obtained in the total eclipse. The residuals are also shown, and shifted by 0.4 in the units of the normalised continuum for convenience. disentangled spectrum. Most of the eclipse spectra were obtained in the phases of ingress or egress of the eclipses, but two of the spectra are very close to the mi… view at source ↗

**Figure 8.** Figure 8: DFT of the light curves for different stages of frequency extraction. The top panel shows the DFT (in black) of the full TESS light curve with eclipses. The panel second from the top (in purple) shows the DFT of the eclipse-subtracted light curve. The second from the bottom (in green) is the pre-whitened DFT after the removal of the dominant pulsation frequency. The bottom-most panel shows the pre-whitened… view at source ↗

**Figure 10.** Figure 10: Échelle diagram of the final set of extracted pulsation frequencies. The échelle phase has been calculated for the orbital frequency. The marker size for the frequencies is proportional to their amplitude. The dominant pulsation frequency is marked in purple, the final set of frequencies is marked in green, and the grey circles are frequencies that were rejected in our analysis. 0 1 2 3 4 Orbital period (… view at source ↗

**Figure 9.** Figure 9: Pulsation runs for 𝑓1 (top; purple), and 𝑓2 (bottom; green). The dominant pulsation frequency, 𝑓1, has slight amplitude variations at primary and secondary eclipses, but not substantial enough to confirm tidal perturbation. which analytically describes the eclipse profile as 𝑓 (𝑡𝑖 , Θ) = 𝛼0 + ∑︁𝑛e 𝑘=1 𝛼𝑘 𝜓(𝑡𝑖 , 𝑇𝑘 , 𝑑𝑘 , Γ𝑘 , 𝐶𝑘 ), (3) where 𝛼0 is the flux zero-point shift, 𝑛e is the number of eclipses du… view at source ↗

**Figure 11.** Figure 11: Roche-filling factors (𝑅/𝑅Roche) versus orbital periods for different tidally perturbed systems with mass and radius measurements in the literature. V446 Cep is marked with a star. The grey lines show the mean and standard deviations for the Roche-filling factors of all pulsators (both secondary and primary, in case the source of pulsation is not identified) in the sample. MNRAS 000, 1–20 (2026) view at source ↗

**Figure 12.** Figure 12: Timing variations for primary eclipses (in red), secondary eclipses (in blue), and the period change (PC) of the dominant pulsation frequency (in purple; shifted for clarity). The lines show sinusoidal fits with a period corresponding to 2310 days. 5.3 Combining radial velocities Using traditional RV extraction methods, we found it difficult to extract the faint secondary. But we were able to extract the … view at source ↗

**Figure 14.** Figure 14 view at source ↗

**Figure 15.** Figure 15: Best-fitting Kurucz SED for the binary+compact object case (left), main sequence triple case (centre), and the quadruple case where the third and fourth star have the same SED (right). The yellow lines show the component SEDs, while the red line shows the combination SED. The primary has the maximum flux contribution so its SED is close to the total SED. The dots represent the modelled fluxes and the cros… view at source ↗

**Figure 16.** Figure 16: Best-fitting MIST isochrones for the binary+compact object case (left), main sequence triple case (centre), and the quadruple case (right). The red isochrone shows the best fit for the pulsating binary system, while the blue line shows an isochrone of 2 Myr. The stars are represented as hollow circles. Their error bars are smaller than their markers view at source ↗

read the original abstract

$\beta$ Cep stars in eclipsing binary (EB) systems give us an opportunity to put observational constraints on their structure and stellar parameters. We present a comprehensive analysis of the $\beta$ Cep star in the EB V446 Cep, using \textit{TESS} photometry and HERMES spectra. We calculate the stellar and orbital parameters using light curve modelling and spectral disentangling. The EB has an orbital period of $3.808567 \pm 0.000012$ d and a mass ratio of $0.1550 \pm 0.0012$. We find the $\beta$ Cep star to have a mass of $10.68 \pm 0.06$ $\mathrm{M_{\odot}}$, a radius of $5.864 \pm 0.033$ $\mathrm{R}_{\odot}$, and a $T_{\rm eff}$ of $24220 \pm 180$ K. The secondary has a mass of $1.657 \pm 0.017$ $\mathrm{M_{\odot}}$, a radius of $1.530 \pm 0.014$ $\mathrm{R}_{\odot}$, and a $T_{\rm eff}$ of $9080 \pm 390$ K. We also extract the abundances of C, N, O, Mg, and Si for the $\beta$ Cep star, which are found to be consistent with galactic OB binaries. We identified 21 distinct pulsation frequencies, with the dominant mode at 10.24324 d$^{-1}$, which corresponds to a near-harmonic of the system's orbital frequency. The two stars in the EB have asynchronous rotation, with both stars rotating faster than the orbital frequency. We detect a companion to the EB using eclipse timing variations and period changes of the dominant pulsation frequency. We calculate the minimum mass of this tertiary companion to be $4.11 \pm 0.32$ $\mathrm{M_{\odot}}$ which is on an orbit of 2303$\pm$69 d around the EB. Using spectral energy distributions and MIST isochrones, we conclude that V446 Cep is either a co-evolving hierarchical 2+2 quadruple or a triple system where the third body is a compact object.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

V446 Cep gets careful parameters and a tertiary detection from TESS timing, but the 2+2 vs compact-object call rests on assumptions about the signals that aren't fully tested.

read the letter

The paper delivers a solid characterization of the V446 Cep eclipsing binary. They use TESS photometry and HERMES spectra to get masses of 10.68 and 1.66 solar masses, radii, temperatures, and abundances for the components, plus 21 pulsation frequencies with the dominant mode at 10.24324 d^{-1}. The orbital period is pinned down to high precision and they extract a 2303-day signal from both eclipse timing variations and changes in the main pulsation period, yielding a minimum tertiary mass of 4.11 solar masses. That part is new and useful for anyone tracking β Cep stars in binaries. The light-curve modeling and spectral disentangling look straightforward and the quoted uncertainties are reasonable given the data quality. The asynchronous rotation and near-harmonic frequency are noted without overclaiming. The soft spot is the interpretation of the tertiary. The authors attribute both the ETVs and the pulsation timing shift to light-time effect from a single companion on that 2303-day orbit, then use MIST isochrones plus SED to argue it is either a co-evolving 2+2 quadruple or a triple with a compact object. They do not show that apsidal motion in the 3.8-day binary or intrinsic frequency evolution in the pulsator are negligible, and the 0.3% proximity to the 39th orbital harmonic is left as a note rather than tested for resonant coupling. If those alternatives contribute even modestly, the minimum mass and the quadruple-versus-triple distinction weaken. The parameter extraction itself does not depend on that step, so the core results stand. This is the kind of detailed case study that belongs in the stellar astrophysics literature. Readers working on massive pulsators or binary evolution will find the numbers and frequency list worth having. It is not transformative but it is cleanly executed on the observational side. I would send it to a referee who knows β Cep timing work; the analysis is detailed enough to deserve that step even if the tertiary conclusion needs tightening.

Referee Report

2 major / 3 minor

Summary. The paper presents a comprehensive analysis of the β Cep pulsator V446 Cep in an eclipsing binary (EB) using TESS photometry and HERMES spectra. It derives the EB orbital period (3.808567 d), mass ratio (0.1550), and precise parameters for the primary β Cep star (10.68 M⊙, 5.864 R⊙, 24220 K) and secondary (1.657 M⊙, 1.530 R⊙, 9080 K) via light-curve modeling and spectral disentangling. The work identifies 21 pulsation frequencies (dominant at 10.24324 d^{-1}, near an orbital harmonic), extracts abundances consistent with galactic OB binaries, notes asynchronous rotation, and detects a tertiary via eclipse timing variations (ETVs) and pulsation period changes, yielding a minimum tertiary mass of 4.11 M⊙ on a 2303 d orbit. Using MIST isochrones and SED fitting, the authors conclude the system is either a co-evolving hierarchical 2+2 quadruple or a triple with a compact tertiary.

Significance. If the timing variations are robustly attributable to a single tertiary's light-time effect and the isochrone/SED distinction holds, this adds a well-characterized β Cep star in a multiple system to the limited sample, with direct constraints on mass, radius, temperature, and pulsation properties that can test stellar models. The parameter determinations from combined photometry and spectroscopy appear solid with quoted errors, and the minimum tertiary mass is derived independently from timing data.

major comments (2)

[Detection of the tertiary companion] In the section on detection of the tertiary companion: The central claim that both the ETVs and the observed changes in the dominant pulsation frequency (10.24324 d^{-1}) arise from the light-time effect of one companion on a 2303±69 d orbit requires explicit quantification of alternative contributions. The frequency lies within ~0.3% of the 39th orbital harmonic (~10.2414 d^{-1}), raising the possibility of resonant coupling or intrinsic evolution not fully ruled out; the manuscript should test whether apsidal motion in the 3.808567 d EB or other effects can be excluded at the reported precision.
[Isochrone fitting and SED analysis] In the isochrone fitting and SED analysis: The distinction between a luminous tertiary (implying 2+2 quadruple) and a non-luminous compact object relies on MIST isochrones applied to the primary parameters (10.68 M⊙, 5.864 R⊙, 24220 K) plus SED. Systematic offsets in isochrones at ~10 M⊙ and the sensitivity of the fit to assumed age/metallicity or tertiary flux contribution should be explored with explicit tests, as these directly support the headline conclusion about the system's nature.

minor comments (3)

[Abstract] The abstract states that abundances of C, N, O, Mg, and Si are consistent with galactic OB binaries, but specific values, uncertainties, and the comparison method (e.g., to a table or literature) are not detailed in the summary; include these for completeness.
[Pulsation analysis] A table listing all 21 pulsation frequencies, amplitudes, phases, and any mode identifications would enhance clarity and allow independent assessment of the frequency spectrum and the near-harmonic claim.
Ensure consistent reporting of uncertainties (e.g., number of significant figures) across orbital elements, stellar parameters, and the tertiary mass in all tables and text.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive summary and constructive major comments. We address each point below and have revised the manuscript to strengthen the analysis of the tertiary companion and the isochrone/SED conclusions.

read point-by-point responses

Referee: In the section on detection of the tertiary companion: The central claim that both the ETVs and the observed changes in the dominant pulsation frequency (10.24324 d^{-1}) arise from the light-time effect of one companion on a 2303±69 d orbit requires explicit quantification of alternative contributions. The frequency lies within ~0.3% of the 39th orbital harmonic (~10.2414 d^{-1}), raising the possibility of resonant coupling or intrinsic evolution not fully ruled out; the manuscript should test whether apsidal motion in the 3.808567 d EB or other effects can be excluded at the reported precision.

Authors: We agree that alternative explanations merit explicit discussion. In the revised manuscript we have added a subsection that estimates the apsidal precession period from the derived masses, radii and orbital eccentricity; the resulting timescale exceeds 10^4 years and is incompatible with the observed 2303 d signal at the reported precision. We also show that the amplitude and phase of the pulsation-frequency variations match the light-time effect predicted from the independent ETV solution for the same orbital elements. While resonant coupling with the near-harmonic cannot be fully excluded without detailed non-adiabatic modeling, the agreement between two independent observables (ETVs and pulsation timing) provides strong support for the light-time interpretation. These additions are now included. revision: partial
Referee: In the isochrone fitting and SED analysis: The distinction between a luminous tertiary (implying 2+2 quadruple) and a non-luminous compact object relies on MIST isochrones applied to the primary parameters (10.68 M⊙, 5.864 R⊙, 24220 K) plus SED. Systematic offsets in isochrones at ~10 M⊙ and the sensitivity of the fit to assumed age/metallicity or tertiary flux contribution should be explored with explicit tests, as these directly support the headline conclusion about the system's nature.

Authors: We have performed the requested sensitivity tests. The revised manuscript now includes isochrone fits using both MIST and PARSEC grids, with age varied between 10–30 Myr and [Fe/H] between −0.1 and +0.1; the primary parameters remain consistent with a co-eval age of ~15–20 Myr in both grids. For the SED, we explicitly model three cases of tertiary flux contribution (0 %, 10 %, 20 % of total flux) and demonstrate that a luminous tertiary of mass comparable to the primary is required to produce a detectable excess, while zero contribution is consistent with a compact object. These tests are presented in a new figure and accompanying text. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The stellar and orbital parameters are obtained via direct light-curve modeling of TESS photometry and spectral disentangling of HERMES spectra, yielding masses, radii, and temperatures without reference to the tertiary claim. The 2303-day tertiary orbit and 4.11 M⊙ minimum mass follow from separate eclipse timing variations combined with observed period changes in the dominant pulsation frequency (interpreted as light-time effect); these are independent timing observables, not a fit renamed as a prediction. The final distinction between a co-evolving 2+2 quadruple and a triple with compact tertiary is reached by applying standard MIST isochrones and SED fitting to the already-derived binary parameters. No self-definitional relations, fitted inputs presented as predictions, load-bearing self-citations, or smuggled ansatzes appear in the chain; the derivation remains self-contained against external isochrone and timing benchmarks.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 1 invented entities

The central claim rests on fitted orbital and stellar parameters from photometry and spectroscopy, plus the assumption that timing variations indicate a bound companion rather than other effects like apsidal motion.

free parameters (3)

orbital period = 3.808567 d
Fitted from light curve modeling
mass ratio = 0.1550
From spectral disentangling and light curve
tertiary orbital period = 2303 d
From eclipse timing variations

axioms (2)

domain assumption Standard stellar atmosphere and pulsation models apply to β Cep stars
Used for parameter extraction and frequency identification
ad hoc to paper Eclipse timing variations are caused by a third body
Assumed to calculate minimum mass

invented entities (1)

tertiary companion no independent evidence
purpose: To explain eclipse timing variations and pulsation period changes
No direct detection, only inferred from timing

pith-pipeline@v0.9.0 · 5774 in / 1433 out tokens · 81540 ms · 2026-05-08T03:07:39.003554+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

4 extracted references · 3 canonical work pages

[1]

doi:10.1007/978-1-4020-5803-5 , adsurl =

Aerts C., 2021, Reviews of Modern Physics, 93, 015001 Aerts C., Christensen-Dalsgaard J., Kurtz D. W., 2010, Asteroseismology. Springer Dordrecht, doi:10.1007/978-1-4020-5803-5 Aller A., Lillo-Box J., Jones D., Miranda L. F., Barceló Forteza S., 2020, A&A, 635, A128 Aschenbrenner P., Przybilla N., Butler K., 2023, A&A, 671, A36 Asplund M., Grevesse N., Sa...

work page doi:10.1007/978-1-4020-5803-5 2021
[2]

has the capability to visualise intrinsic stellar effects like rotational broadening, spots, and pulsations (Moharana et al. 2023). To spot similar variations in the𝛽Cep component of V446 Cep, we calculated BF using the HERMES spectra. The BF was calculated using the algorithm de- scribedinRucinski(1999).Wemodifiedasingle-orderBFcode,bf- rvplotter14, to c...

2023
[3]

APPENDIX E: ASTROMETRY SIMULATIONS OF SYSTEMS WITH KNOWN 3D GEOMETRY Wecompiledalistofcompacthierarchicaltriples(CHT),whichhave precisely measured parameters that define a system’s 3D geometry. FromBorkovitsetal.(2020,2022);Moharanaetal.(2023,2024),and 14 https://github.com/mrawls/BF-rvplotter MNRAS000, 1–20 (2026) V446 Cep is a𝛽Cep star in a multiple sys...

work page arXiv 2020
[4]

APPENDIX F: ECLIPSE TIMINGS InTable F1,we providethe timesof primaryand secondaryminima of V446 Cep derived from theTESSlight curve using the procedure described in Section 5.1

This showsthatmodelswithΔ RUWElessthan0.3areacceptablesolutions for V446 Cep. APPENDIX F: ECLIPSE TIMINGS InTable F1,we providethe timesof primaryand secondaryminima of V446 Cep derived from theTESSlight curve using the procedure described in Section 5.1. APPENDIX G: PERIOD CHANGE MEASUREMENTS The period change measurements for the dominant period are giv...

work page arXiv 2026