Tri-axial Pulsators in Binary Stars

• Tri-axial pulsators are tidally distorted stars in close binaries whose oscillation modes align with three mutual perpendicular axes (x, y, and z), distinguishing them from isolated pulsators.
• They exhibit orbital-phase-dependent amplitude and phase modulations, enabling clear identification of dipole and quadrupole pulsation modes aligned with tidal and rotational axes.
• Perturbative analyses show that tidal forces and rotational flattening in synchronized binaries create standing-wave geometries, advancing three-dimensional asteroseismology in such systems.

Searching arXiv for papers on tri-axial pulsators and related tidal asteroseismology. Tri-axial pulsators are pulsating stars in close binaries whose nonradial oscillation modes are organized about three mutually orthogonal axes of a tidally distorted star rather than about a single symmetry axis. In the current formulation, these axes are the tidal axis $x$, the in-plane axis $y$ perpendicular to both the tidal axis and the orbital angular momentum axis, and the orbital/rotation axis $z$. The class emerged observationally with TIC 184743498 and TIC 435850195, and was placed on a broader theoretical footing by perturbative calculations showing that synchronized tidally distorted stars can support standing-wave pulsation geometries aligned with the three principal axes of a triaxial ellipsoid (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024). Subsequent analysis of EL CMi extended the picture beyond dipole modes by reporting the first detected quadrupole Tidally Tilted Standing mode (Handler et al., 28 Jul 2025).

1. Definition and relation to tidally tilted pulsation

The defining property of a tri-axial pulsator is that its pulsation spectrum cannot be explained by modes aligned with only one axis, or even two axes, but instead requires three distinct pulsation axes in the same star. In the observational papers, these are identified as dipole-mode families $Y_{10x}$, $Y_{10y}$, and $Y_{10z}$, where $x$ points toward the companion, $y$ lies in the orbital plane perpendicular to $x$, and $z$ is the orbital/rotation axis (Zhang et al., 2023, Jayaraman et al., 2024).

This class is distinct from ordinary isolated pulsators, whose mode geometry is usually fixed relative to the stellar spin axis and therefore does not generate orbital-phase-dependent multiplet structure. It is also more general than the earlier tidally tilted pulsator picture. Tidally tilted pulsators in close binaries show pulsation axes tilted into the orbital plane by the tidal bulge, but tri-axial pulsators exhibit modes aligned along multiple orthogonal axes, including the $y$0-axis that is perpendicular to both the tidal axis and the orbital angular momentum axis (Jayaraman et al., 2024). In that sense, tri-axial pulsation is not merely an obliquity effect; it is a mode-geometry consequence of stellar triaxiality.

The 2024 perturbative analysis generalized this further by arguing that tidally distorted stars are, in the relevant p-mode regime, intrinsically triaxial pulsators. In that treatment, the equilibrium figure is set by the combined action of tidal elongation and rotational flattening, so the oscillation eigenfunctions inherit the star’s three principal ellipsoidal axes rather than a single rotational symmetry axis (Fuller et al., 2024).

2. Dynamical origin in tidally distorted synchronized binaries

The theoretical framework starts from linear perturbation theory for oscillations in a distorted star, written as a generalized eigenvalue problem involving the kinetic-energy operator $y$1, the potential-energy operator $y$2, and the Coriolis operator $y$3 (Fuller et al., 2024):

$y$4

In the spherical limit, each mode is labeled by a single spherical harmonic $y$5. Rotation and tides perturb that description in qualitatively different ways. The Coriolis force is axisymmetric about the spin axis and only couples modes with the same $y$6, recovering the familiar rotational splitting at first order. By contrast, tidal and centrifugal distortions break spherical symmetry and introduce couplings between different $y$7 values within a multiplet. Because the dominant tidal distortion is quadrupolar, the crucial couplings are between $y$8 values differing by $y$9 or $z$0, with the $z$1 components producing the off-diagonal coupling between $z$2 and $z$3 that is central to dipole triaxiality (Fuller et al., 2024).

For $z$4 pressure modes in synchronized binaries, the tidal coupling can dominate over Coriolis splitting. The $z$5 component remains effectively uncoupled and becomes the $z$6 mode aligned with the spin/orbital axis. The $z$7 pair mixes into two standing-wave combinations aligned with the two horizontal axes. In the observationally motivated notation used across the literature, these are the $z$8 and $z$9 modes (Zhang et al., 2023, Fuller et al., 2024). The key geometric outcome is therefore a dipole triplet reorganized into three orthogonal axis families rather than the usual $Y_{10x}$0 rotationally split set.

A closely related formulation appears in the stellar-discovery papers, where the coupled $Y_{10x}$1 problem is written in terms of the amplitudes $Y_{10x}$2 and $Y_{10x}$3, with tidal terms $Y_{10x}$4 and $Y_{10x}$5 coupling the $Y_{10x}$6 pair. In the strong-coupling limit, the mixed $Y_{10x}$7 components become standing waves proportional to the Cartesian directions $Y_{10x}$8 and $Y_{10x}$9, while the $Y_{10y}$0 component remains proportional to $Y_{10y}$1 (Zhang et al., 2023). The later EL CMi work described the same physical objects as Tidally Tilted Standing (TTS) modes, emphasizing that the spatial pattern is nearly fixed relative to the companion rather than propagating around the star (Handler et al., 28 Jul 2025).

A major theoretical conclusion is the separation between p-mode and g-mode behavior. The perturbative analysis argues that pressure modes in sufficiently close synchronized binaries should often be strongly tidally tilted and therefore triaxial, whereas gravity modes should generally remain aligned with the stellar spin axis because Coriolis effects dominate their geometry (Fuller et al., 2024). This suggests that tri-axial pulsation is not expected to be a universal property of all mode classes in all close binaries.

3. Observational diagnostics and mode identification

The primary observational signature of tri-axial pulsation is orbital-phase-dependent amplitude and phase modulation. Because the pulsation pattern is fixed in the binary frame, the observer samples a changing aspect angle over the orbit. For the $Y_{10y}$2- and $Y_{10y}$3-aligned dipole modes, the observed amplitude passes through minima where the visible hemisphere changes sign, and the phase correspondingly jumps by about $Y_{10y}$4 cycles or $Y_{10y}$5 radians. In Fourier space, this modulation produces pairs of peaks separated by exactly $Y_{10y}$6; by contrast, the $Y_{10y}$7-aligned dipole mode has no orbital modulation and appears as a singlet (Fuller et al., 2024).

The discovery paper for TIC 184743498 supplied the canonical empirical version of this signature. The star is a $Y_{10y}$8 Scuti pulsator in a tight eclipsing binary with $Y_{10y}$9, and each of the nine principal pulsation modes shows amplitude modulation and $Y_{10z}$0-rad phase shifts twice per orbital cycle (Zhang et al., 2023). Five modes have amplitude maxima at eclipses and phase jumps near the ellipsoidal light variation peaks, matching dipole modes aligned with the tidal axis $Y_{10z}$1. Four others are shifted by $Y_{10z}$2 in orbital phase, with amplitude maxima at the ellipsoidal-variation peaks and phase jumps near eclipse, identifying them with the orthogonal in-plane $Y_{10z}$3-axis (Zhang et al., 2023).

The echelle diagram is a practical tool in this context. In TIC 435850195, the echelle phase is defined as

$Y_{10z}$4

so peaks separated by integer multiples of $Y_{10z}$5 align in the diagram (Jayaraman et al., 2024). A dipole doublet is then operationally a pair of peaks sharing the same echelle phase and separated by $Y_{10z}$6. That structure is the principal observational hallmark of the standing-wave $Y_{10z}$7 and $Y_{10z}$8 modes.

An important diagnostic use of the amplitude and phase behavior is the exclusion of alternative mode identifications. In TIC 184743498, the four modes with maxima at ellipsoidal-variation peaks might superficially resemble $Y_{10z}$9-type modes, but the paper argues that such an interpretation fails at the measured inclination $x$0, because a tidally tilted $x$1 mode would produce a strong central peak that is not observed (Zhang et al., 2023). In TIC 435850195, the analogous $x$2 explanation is excluded because the observed doublets lack the detectable central peak expected at the measured inclination, while rotational splitting is rejected because the splittings are exactly $x$3 to very high precision and would require an implausibly small Ledoux constant; eclipse mapping is also found insufficient to explain the strong doublets (Jayaraman et al., 2024). These system-specific exclusions established that tri-axial classification rests on geometry-sensitive diagnostics rather than on peak counting alone.

4. Empirical systems and the emergence of the class

The first two confirmed systems established the empirical definition of the class and showed that the $x$4-axis family is not an isolated peculiarity but a repeatable phenomenon.

System	Core observational features	Mode interpretation
TIC 184743498	Tight binary with $x$5; nine principal modes between 38 and 56 d$x$6; $x$7-rad phase shifts twice per orbit	Five $x$8, four $x$9, two singlet-like $y$0 (Zhang et al., 2023)
TIC 435850195	Sixteen robustly detected multiplets; fourteen dipole doublets split by $y$1; systematic orbital modulation	Eight $y$2, six $y$3, additional likely $y$4 modes (Jayaraman et al., 2024)
EL CMi	Orbital sideband multiplets in TESS photometry; stable orbital phase modulation; binary characterized with RVs and PHOEBE2	$y$5, $y$6, and quadrupole $y$7 TTS mode (Handler et al., 28 Jul 2025)

TIC 184743498 was identified as the first tri-axial stellar pulsator because its oscillation spectrum required three pulsation axes in the same star. The analysis used TESS Sectors 61 and 62 after removing the first 30 orbital harmonics. The star exhibits nine dominant pulsation modes between about 38 and 56 d$y$8, organized into 11 multiplets or doublets in the echelle diagram, plus two singlet modes. The two singlets near 52.5418 and 55.6671 d$y$9 show little or no significant orbital modulation in amplitude or phase and are interpreted as $x$0-like modes rather than radial modes, because the inferred radial-mode spectrum is not dense enough to place both as radial overtones without requiring a much more evolved star than implied by the SED and binary constraints (Zhang et al., 2023).

TIC 435850195 provided the second confirmed case and sharpened the classification criteria. The light curve from TESS Sector 56 shows primary eclipses, strong ellipsoidal variations, and orbital-phase-dependent pulsation amplitudes and phases. The joint fit to the TESS light curve and the SED from Gaia, Pan-STARRS1, WISE, 2MASS, Tycho-2, SDSS, and GALEX yielded a primary with $x$1, $x$2, $x$3 K, and $x$4, and a secondary with $x$5, $x$6, $x$7 K, and $x$8, at an age of $x$9 Myr and inclination $z$0 (Jayaraman et al., 2024). The primary is described as slightly evolved off the zero-age main sequence and fills a bit over 50% of its Roche lobe, placing it in the regime where $z$1 Scuti pulsation and strong tidal distortion coexist.

A notable empirical advance in TIC 435850195 is the explicit identification of six $z$2 modes as the novel tri-axial class. Their amplitude minima occur at primary and secondary eclipse, their maxima at the ellipsoidal-variation maxima, and they show the expected $z$3 phase jump at minimum amplitude (Jayaraman et al., 2024). This behavior is the $z$4-shifted analogue of the $z$5 pattern and therefore the cleanest direct signature of pulsation about the orthogonal in-plane axis.

5. Extension to Tidally Tilted Standing modes and EL CMi

EL CMi extended tri-axial pulsation from a dipole-only observational class to one that includes higher-degree TTS modes. The system is an eclipsing mass-transferring binary containing a $z$6 Scuti pulsator as the hotter component. TESS observations from Sectors 7, 34, 61, and 88 revealed dominant mode groups $z$7, $z$8, and $z$9, plus a weaker $y$00, with orbital sideband multiplets such as $y$01, $y$02, and $y$03 (Handler et al., 28 Jul 2025).

The amplitude and phase variations over orbital phase were reconstructed using the method of Jayaraman et al. (2022), and the crucial result was that although amplitudes vary somewhat from sector to sector, the phase modulation over the orbit is essentially the same in all sectors. This enabled robust mode identification: $y$04 was identified as a dipole $y$05 mode, $y$06 as a dipole $y$07 mode, and $y$08 as a quadrupole TTS mode of type $y$09 (Handler et al., 28 Jul 2025). The paper describes this as the first observational detection of a quadrupole TTS mode.

In the EL CMi analysis, the quadrupole label

$y$10

denotes a standing pattern with maxima and minima offset from the tidal axis by an angle of $y$11, rather than a traveling $y$12 wave around the rotation axis (Handler et al., 28 Jul 2025). This is consistent with the general theoretical statement that higher-degree modes in tidally distorted stars are also standing modes, but with more complex observable modulation patterns than the dipole case (Fuller et al., 2024).

The binary characterization reinforced the astrophysical significance of the detection. New spectroscopy with NOT/ALFOSC gave $y$13 and $y$14, while PHOEBE2 modeling of the TESS light curve and radial velocities yielded approximately $y$15, $y$16, $y$17, $y$18, $y$19 K, $y$20 K, $y$21, and $y$22 (Handler et al., 28 Jul 2025). A separate SED plus light-curve analysis with Gaia DR3 parallax, extinction estimates, broad-band photometry, and MIST tracks gave a broadly consistent picture and supported the conclusion that the secondary is a stripped donor star undergoing ongoing low-rate mass transfer and likely evolving into a low-mass helium white dwarf (Handler et al., 28 Jul 2025). This showed that tri-axial pulsation is relevant not only to detached synchronized binaries but also to post-mass-transfer or weakly mass-transferring systems.

6. Asteroseismic significance, ambiguities, and scope of the term

The central asteroseismic consequence of tri-axial pulsation is that orbital amplitude and phase modulation provide an additional mode-identification channel. In the dipole case, exact $y$23 doublets with characteristic phase reversals distinguish the $y$24- and $y$25-aligned modes, while singlets identify the $y$26-aligned member of the triplet (Fuller et al., 2024). This can make mode identification easier than in isolated pulsators, especially for $y$27 modes, because the binary orbit supplies a geometrical modulation that labels the mode family.

At the same time, the theory paper emphasizes that the new geometry can introduce fresh ambiguities. The $y$28 and $y$29 modes can look very similar observationally, as can $y$30 and $y$31, while $y$32 may masquerade as a radial mode or as $y$33 because all appear as singlets without orbital modulation (Fuller et al., 2024). A plausible implication is that future tri-axial asteroseismology will require simultaneous modeling of tidal coupling, centrifugal distortion, rotation, and eclipse geometry rather than relying on conventional single-star mode taxonomy.

The theory also situates tri-axial pulsation within a wider continuum of tidal mode geometry. In the moderate-distortion regime, the natural outcome is triaxial standing modes aligned with the star’s three ellipsoidal axes. In more strongly distorted systems, higher-order tidal terms can produce “single-sided” or tidally trapped pulsations that are asymmetric across the $y$34-$y$35 plane (Fuller et al., 2024). This suggests that tri-axial pulsation is not an isolated anomaly but part of a broader family of non-axisymmetric binary-star oscillation phenomena.

A recurring misconception is to treat “tri-axial pulsator” as a generic label for any system with three prominent pulsation frequencies. In the current literature, the term is much narrower: it denotes close-binary pulsators whose modes are demonstrably organized by tidal geometry about three perpendicular axes and whose orbital-phase-dependent amplitude and phase behavior reveals that geometry (Zhang et al., 2023, Jayaraman et al., 2024). Another possible source of confusion is terminological rather than astrophysical. Outside stellar pulsation, “tri-axial” has been used for three-axis magnetoresistance in $y$36 nanodevices and for tri-axial time-dependent magnetic-field calibration using adiabatically evolving atomic spins, where it refers to three-dimensional response or control rather than stellar oscillation geometry (Zhang et al., 2018, Bevilacqua et al., 2024). In astrophysics, by contrast, the term specifically denotes a tidally organized pulsation eigenbasis in a close binary.

Taken together, the observational discoveries and the perturbative theory establish tri-axial pulsators as a distinct close-binary asteroseismic class. The empirical sequence from TIC 184743498 to TIC 435850195 and EL CMi indicates that three-axis tidal mode organization is reproducible, diagnostically accessible in TESS-quality photometry, and extensible from dipole modes to at least one quadrupole TTS mode (Zhang et al., 2023, Jayaraman et al., 2024, Handler et al., 28 Jul 2025). This suggests a route toward genuinely three-dimensional asteroseismology of tidally distorted stars, provided that the mode physics is modeled in the binary frame rather than imported unchanged from isolated-star pulsation theory.