Osculating Stars in Binary Asteroseismology

• Osculating stars are binary systems in which tidal forces couple stellar pulsations with orbital perturbations, enabling unified asteroseismic analyses.
• Researchers use high-precision light-curve modeling, radial velocities, and Fourier spectrum analysis to identify pulsation modes and decode orbital dynamics.
• Integrating binary modeling with seismic inferences yields actionable insights into stellar interiors, mass transfer dynamics, and orbital evolution.

“Osculating stars” — Editor’s term — may be used for stellar systems in which stellar oscillations, tides, and orbital perturbations must be analyzed together rather than as separable phenomena. In current arXiv literature, the clearest instances are tidally tilted pulsators in very close binaries, heartbeat stars with tidally excited oscillations, eccentric binaries containing solar-like oscillators, and Roche-lobe-overflowing binaries whose evolution is computed with osculating orbital theory. Across these cases, the common feature is that mode visibility, mode excitation, and orbital evolution are all conditioned by the companion’s gravity, so that asteroseismology and binary dynamics become a single inference problem (Handler et al., 2022).

1. Scope and phenomenological classes

The most sharply defined subclass is the tidally tilted pulsator: a pulsating component of a very close binary whose pulsation axis has been pulled out of the stellar rotation axis and into the orbital plane by the companion’s tidal field. Handler et al. described this group on the basis of the first three representatives, HD 74423, CO Cam, and TIC 63328020. In these systems the pulsation axis lies in, or close to, the orbital plane, so the star is effectively viewed under all possible pulsation–inclination angles over one orbit (Handler et al., 2022).

A second major subclass is the heartbeat star. Heartbeat stars are eccentric binary stars in short period orbits whose light curves are shaped by tidal distortion, reflection, and Doppler beaming. Some heartbeat stars exhibit tidally excited oscillations, making them laboratories for the forced response of stellar modes to orbital harmonics. KOI-54 is the canonical example in which the light curve contains strong ellipsoidal variability during periastron passage together with pulsations at perfect harmonics of the orbital frequency and additional nonharmonic pulsations (Fuller, 2017).

A third class comprises eccentric binaries with solar-like oscillating red giants. In these systems, stochastic p-mode excitation by convection coexists with binary light-curve signatures such as ellipsoidal modulation and periastron brightening. Kepler data identified 18 eccentric red-giant binaries with $20\,{\rm d}<P<440\,{\rm d}$ and $0.2Beck et al., 2014).

A fourth, dynamically distinct class comprises interacting binaries modeled with the theory of osculating orbits. Here the focus is not on pulsation visibility but on the secular evolution of the orbit under perturbing forces arising from Roche lobe overflow, stream gravity, and linear-momentum exchange. In conservative Algol evolution, the osculating prescription predicts substantially tighter post-mass-transfer binaries than the classical $J_{\rm orb}={\rm const}$ prescription (Davis et al., 2014).

Class	Defining feature	Representative systems
Tidally tilted pulsators	Pulsation axis in the orbital plane	HD 74423, CO Cam, TIC 63328020
Heartbeat / TEO systems	Eccentric binaries with tidally excited oscillations	KOI-54, KIC 5006817
Oscillating red-giant binaries	Solar-like oscillations plus eccentric-binary modulation	Kepler heartbeat giants
Osculating-orbit interacting binaries	RLOF evolution with perturbing stream forces	Conservative Algols

2. Geometry, forcing, and the tidal control of oscillations

The defining geometry of tidally tilted pulsators is set by the orbital inclination $i$, the obliquity $\beta$ of the pulsation axis relative to the orbital normal, and the orbital phase $\phi_{\rm orb}(t)$, with $\phi_{\rm orb}=0$ at superior conjunction. Since the pulsation axis lies in the orbital plane, $\beta \approx 90^\circ$, and the instantaneous angle between the pulsation axis and the line of sight satisfies

$$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$

which for $\beta=90^\circ$ simplifies to

$0.2

As $0.2Handler et al., 2022).

The tidal forcing itself is naturally expressed in spherical harmonics. For a close companion, the time-dependent tidal potential in the rotating frame of the primary is

$0.2

with the $0.2$J_{\rm orb}={\rm const}$0 acoustic modes. When the forcing frequency $J_{\rm orb}={\rm const}$1 approaches a free acoustic or gravity-mode frequency, the mode is resonantly excited (Handler et al., 2022).

In eccentric binaries the forcing becomes strongly phase concentrated. For KOI-54, the companion’s tidal potential is sharply peaked near periastron, producing a brief “periastron kick” of the stellar envelope. This periastron pumping excites the dynamical tide, that is, standing g-mode oscillations whose natural frequencies happen to lie near integer harmonics of the orbital frequency. In that framework the spatial coupling is measured by the overlap integral $J_{\rm orb}={\rm const}$2, the temporal coupling by the Hansen coefficient $J_{\rm orb}={\rm const}$3, and the response is controlled by the Lorentzian factor

$J_{\rm orb}={\rm const}$4

so that amplitudes and phases encode detuning and damping rather than the orbital frequencies themselves (Burkart et al., 2011).

3. Mode visibility, multiplets, and identification

A central observational advantage of tidally tilted pulsators is that mode identification can be carried out geometrically. Because $J_{\rm orb}={\rm const}$5 is known as a function of orbital phase, the observed amplitude of an $J_{\rm orb}={\rm const}$6 mode may be written

$J_{\rm orb}={\rm const}$7

where

$J_{\rm orb}={\rm const}$8

is the intrinsic visibility integral and $J_{\rm orb}={\rm const}$9 is the limb-darkening law. Over one orbit, both the amplitude and pulsation phase trace signatures that depend only on $i$0, so a simultaneous fit directly yields mode identification. The paper characterizes this as analogous to eclipse-mapping methods and to oblique pulsator analysis in roAp stars; it further states that mode identification is virtually guaranteed by the complete aspect sampling, with no traditional $i$1 ambiguity (Handler et al., 2022).

In the Fourier spectrum of a tidally tilted pulsator, a pulsation frequency $i$2 appears as a multiplet spaced by the orbital frequency $i$3: $i$4 For HD 74423, the dominant axisymmetric $i$5 Scuti mode at $i$6 is accompanied by sidebands at $i$7, with $i$8, and at $i$9. Extraction proceeds by iterative pre-whitening: removal of the binary-orbit term, identification of pulsation peaks in the residuals, and grouping into multiplets with spacing $\beta$0. The centroid and sideband ratios encode $\beta$1, while the orbital-phase dependence of amplitudes and phases yields $\beta$2 (Handler et al., 2022).

For heartbeat stars, the key diagnostic is similar in spirit but different in implementation. In tidal asteroseismology the observed frequencies are exactly harmonics $\beta$3, and the stellar interior is constrained primarily by amplitudes and phases rather than by the frequencies themselves. In KOI-54, the phase of a harmonic relative to periastron is sensitive to detuning and damping, and the nonharmonic pulsations can be produced by nonlinear three-mode coupling. The observed nonharmonic peaks satisfy $\beta$4 for a large harmonic $\beta$5, which is the signature of parametric resonance in the formalism of parent and daughter modes (Burkart et al., 2011).

A related, more general treatment factors out the equilibrium tide, incorporates rotation using the traditional approximation, includes non-adiabatic effects in the surface luminosity perturbation, allows for spin–orbit misalignment, and correctly sums over contributions from many modes. In that framework the total photometric signal at a given harmonic is a coherent sum over all contributing $\beta$6 modes, and the theory can be used to test whether an observed oscillation is consistent with a chance resonance or requires resonance locking (Fuller, 2017).

4. Binary–seismic inference

The principal astrophysical gain of these systems is the combination of binary modeling with asteroseismology. For tidally tilted pulsators, binary modeling from radial velocities, eclipse shapes, and the ellipsoidal light curve supplies precise masses, radii, orbital inclination $\beta$7, and Roche geometry. Combined with straightforward mode identification, this yields tight constraints on stellar interior structure, including the sound-speed profile and convective-core size, in a regime of strong tidal distortion not probed by single stars. The same framework tests which $\beta$8 modes are tidally trapped, constrains tidal realignment timescales, and probes the internal viscosity of A/F stars (Handler et al., 2022).

In red-giant eccentric binaries, the seismic observables are the frequency of maximum power $\beta$9 and the large frequency separation $\phi_{\rm orb}(t)$0. Kepler long-cadence photometry is typically used, since it covers $\phi_{\rm orb}(t)$1 without aliasing. One determines $\phi_{\rm orb}(t)$2 by fitting a Gaussian profile to the oscillation envelope above the background and estimates $\phi_{\rm orb}(t)$3 from the autocorrelation of the power spectrum or from identified radial modes. Standard scaling relations then provide first estimates of mass and radius: $\phi_{\rm orb}(t)$4 Light-curve fitting codes such as PHOEBE and ELC then model period, eccentricity, argument of periastron, inclination, scaled radii, temperature ratio, tidal distortions, and Doppler beaming. The data summary states that asteroseismic priors on $\phi_{\rm orb}(t)$5 and $\phi_{\rm orb}(t)$6 drastically reduce degeneracies in the light-curve fit, yielding precise $\phi_{\rm orb}(t)$7, $\phi_{\rm orb}(t)$8, and $\phi_{\rm orb}(t)$9 (Beck et al., 2014).

The wider asteroseismic context is that $\phi_{\rm orb}=0$0 scales to first order as the square root of the mean density and $\phi_{\rm orb}=0$1 scales with the acoustic cutoff frequency, hence with surface gravity and effective temperature. Extraction of $\phi_{\rm orb}=0$2 and $\phi_{\rm orb}=0$3 from high-precision space photometry, together with $\phi_{\rm orb}=0$4 from spectroscopy or photometry, enables either direct scaling-relation estimates or grid-based Bayesian modeling. Reported typical precisions for red giants are $\phi_{\rm orb}=0$5 in radius, $\phi_{\rm orb}=0$6 in mass, and $\phi_{\rm orb}=0$7–$\phi_{\rm orb}=0$8 in age, though the same source emphasizes calibration issues in $\phi_{\rm orb}=0$9, surface terms, mass loss, uncertain input physics, and selection biases (Miglio et al., 2014).

5. Osculating-orbit theory in interacting binaries

In the strict dynamical sense, the “osculating” component of the subject is the treatment of interacting binaries through osculating orbital elements. In the Davis et al. formalism, the relative motion of donor and accretor is written

$\beta \approx 90^\circ$0

where $\beta \approx 90^\circ$1 and $\beta \approx 90^\circ$2 are perturbing accelerations arising from the change of linear momentum due to ejection and accretion, the gravitational attraction of the stream on each star, and center-of-mass acceleration under asymmetric mass change. The stream is split into material that reaches the companion and material that falls back onto the donor, and the stream-gravity term is obtained by discretizing the ballistic stream into particles and integrating their gravitational pull (Davis et al., 2014).

The osculating elements then evolve according to perturbation equations. For semi-major axis $\beta \approx 90^\circ$3 and eccentricity $\beta \approx 90^\circ$4,

$\beta \approx 90^\circ$5

$\beta \approx 90^\circ$6

with $\beta \approx 90^\circ$7 the true anomaly. In circular orbits, only $\beta \approx 90^\circ$8 contributes to $\beta \approx 90^\circ$9, while $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$0 when averaged over one orbit (Davis et al., 2014).

Spin evolution is likewise coupled to tides and mass transfer: $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$1 For the mass-transfer contribution,

$$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$2

For the donor’s Roche-lobe overflow mass loss,

$$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$3

The principal quantitative result is that post-mass-transfer periods are typically shorter by a factor $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$4 than in the classical point-mass $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$5 prescription. For a $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$6 system with initial $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$7 d, the final period is $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$8 versus $$\cos \theta_{\rm obs}(t)=\sin i\,\cos\beta\,\sin\phi_{\rm orb}(t)+\cos i\,\sin\beta ,$$9; for a $\beta=90^\circ$0 system with initial $\beta=90^\circ$1 d, the final period is $\beta=90^\circ$2 versus $\beta=90^\circ$3. During the rapid RLOF phase, the donor spins down faster than tides can re-synchronize, because $\beta=90^\circ$4–$\beta=90^\circ$5 yr while $\beta=90^\circ$6–$\beta=90^\circ$7 yr. The resulting sub-synchronous rotation causes self-accretion of $\beta=90^\circ$8–0.20 of the ejected stream, further enhancing orbital shrinkage (Davis et al., 2014).

This result has direct bearing on any encyclopedia treatment of “osculating stars,” because it supplies the rigorous dynamical framework for binaries in which the companion’s perturbation is not a small correction to a fixed Kepler problem but part of the primary stellar evolution.

6. Scientific significance, extensions, and unresolved issues

The scientific significance of these systems lies in the fact that they convert orbital geometry into seismic leverage. Tidally tilted pulsators offer nearly complete aspect sampling over a single orbital cycle; heartbeat stars encode mode physics in forced amplitudes and phases; eccentric red-giant binaries allow seismic priors to regularize light-curve and radial-velocity solutions; and osculating-orbit binaries expose the dynamical back-reaction of mass transfer on orbital elements. This suggests that “osculating stars,” in the present editorial sense, are best regarded as a boundary domain between asteroseismology, binary-star modeling, and tidal dynamics rather than as a single pulsator class (Handler et al., 2022).

Several unresolved issues are explicitly identified in the cited literature. In tidal asteroseismology, one must decide whether an observed tidally excited oscillation is explained by a chance resonance with a stellar mode or by a resonance locking process; the formalism of Fuller derives both the probability theory for chance resonances and the locked-amplitude solution. The same work states that tidally excited oscillations are more visible in hot stars with $\beta=90^\circ$9, because g-modes propagate closer to the photosphere and have larger surface flux perturbations there (Fuller, 2017).

At even stronger interaction, asynchronous stellar coalescence produces another extension of the subject. Hydrodynamic simulations show that resonance crossings with high-azimuthal-order fundamental modes, $0.2wave with $0.2Macleod et al., 2018).

A different frontier treats stars themselves as detectors of external perturbations. Sun-like and red-giant oscillators possess low-degree quadrupole modes that can be driven by gravitational waves, potentially probing the $0.2Lopes et al., 2015).

Taken together, these developments define a research program in which orbital perturbations are not merely sources of noise for stellar pulsation studies. They are the organizing principle of the observable signal, the route to mode identification, and, in interacting binaries, part of the secular evolution itself.