Triaxial Pulsators in Close Binaries

• Triaxial pulsators are stars in close binaries whose oscillation modes align along three mutually perpendicular axes due to tidal distortions.
• They exhibit distinct orbital multiplets and phase modulation patterns that differentiate in-plane dipole modes from out-of-plane modes.
• This framework enables detailed three-dimensional asteroseismic analysis, refining mode identification and advancing our understanding of stellar structure.

Searching arXiv for papers on triaxial pulsators and closely related work. Triaxial pulsators are oscillating systems in which the observable pulsation geometry is organized around three distinct axes rather than a single symmetry axis. In contemporary stellar asteroseismology, the term refers most specifically to close-binary pulsators whose nonradial modes are aligned with three mutually perpendicular directions tied to the binary geometry, so that different dipole or higher-order modes behave as pulsations about the tidal axis, an in-plane perpendicular axis, and the orbital angular-momentum axis (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024). In a broader astrophysical usage, “triaxial pulsar” denotes a neutron star treated as a triaxial rigid body, with unequal principal moments of inertia, whose non-axisymmetry gives rise to continuous gravitational-wave emission and, if misaligned with the angular momentum, free or forced precession (Pathak et al., 2022, Gao et al., 2020). The phrase also appears in nuclear collective models for quadrupole motion about triaxial equilibrium shapes, where vibration and rotation coexist around $\gamma \approx 30^\circ$ (Yigitoglu et al., 2010). The modern stellar concept has developed rapidly through TESS discoveries and perturbative mode-coupling theory, culminating in the statement that tidally distorted stars in close synchronized binaries are generically triaxial pulsators (Fuller et al., 2024).

1. Stellar triaxial pulsators in close binaries

In the stellar context, a triaxial pulsator is a pulsating star in a close binary whose equilibrium figure is a triaxial ellipsoid because the rotation or centrifugal distortion is axisymmetric about the spin or orbital axis and the tidal distortion is axisymmetric about the line of centers to the companion (Fuller et al., 2024). The three principal axes are conventionally labeled $x$, $y$, and $z$, with the longest axis $x$ along the tidal line joining the two stars, the shortest axis $z$ along the spin or orbital axis, and the intermediate axis $y$ perpendicular to both (Jayaraman et al., 2024, Fuller et al., 2024). In this setting, low-degree pressure modes are no longer described by a single spherical harmonic $Y_{\ell m}$ about the rotation axis. Instead, each $\ell=1$ triplet is transformed into three eigenmodes aligned with the three principal axes of the triaxial figure, denoted $Y_{10x}$, $x$0, and $x$1 (Jayaraman et al., 2024, Fuller et al., 2024).

The essential observational consequence is orbital-phase-dependent amplitude and phase modulation. In tidally tilted pulsators, the tidal bulge defines a preferred axis in the orbital plane and the observer’s latitude relative to the pulsation axis cycles through $x$2 over the orbit, producing orbital multiplets in the Fourier spectrum (Jayaraman et al., 2024). Triaxial pulsators represent the more extreme case in which tilted modes coexist with modes aligned with the orbital axis, so that the star simultaneously supports oscillations about three perpendicular axes (Zhang et al., 2023, Jayaraman et al., 2024).

This stellar usage differs from the roAp literature. Rapidly oscillating Ap stars are classical laboratories for the oblique pulsator model because their pulsation axes are not aligned with their rotation axes, but the relevant geometry is usually discussed in terms of inclination $x$3 and magnetic obliquity $x$4, distorted dipole or quadrupole modes, and possible departures from axisymmetry rather than explicit tri-axial labeling (Zhong et al., 2024). The roAp case is therefore adjacent to, but not identical with, the close-binary triaxial pulsator framework.

2. Theoretical basis: tidal coupling and standing modes

The perturbative framework for stellar triaxial pulsators computes mode frequencies and geometries of tidally distorted stars while accounting for the Coriolis force and coupling between different azimuthal orders $x$5 induced by tidal distortion (Fuller et al., 2024). In the rotating star, the linear eigenproblem is written in operator form, and the tidal quadrupole $x$6 couples $x$7 components inside a given multiplet through off-diagonal matrix elements, whereas the Coriolis operator is diagonal in $x$8 and provides first-order rotational splitting (Fuller et al., 2024).

For synchronized close binaries, the dominant effect on dipole pressure modes is the coupling of the $x$9 and $y$0 components. In the strong tidal coupling limit, the uncoupled $y$1 mode remains aligned with the rotation axis and becomes the $y$2 mode, while the equal-weight combinations of $y$3 become standing modes aligned with the $y$4 and $y$5 axes (Fuller et al., 2024). The corresponding real-space flux perturbations are written as

$y$6

$y$7

with the $y$8 mode given by

$y$9

so that the three members of the dipole multiplet become standing modes aligned with the three principal axes of the triaxial ellipsoid (Fuller et al., 2024).

A closely related formulation was developed for tidally tilted pulsators and extended to tri-axial pulsators in the analysis of TIC 435850195 (Jayaraman et al., 2024). There, the tidal perturbation is approximated as

$z$0

which, when written in spherical harmonics about the $z$1-axis, couples the dipole modes with $z$2 and produces new eigenmodes

$z$3

$z$4

This formulation makes explicit that the traveling equatorial waves are converted into standing dipole modes whose axes lie in the orbital plane (Jayaraman et al., 2024).

For quadrupole modes, the behavior is more complex. In the general perturbative treatment, the $z$5 quintuplet separates into odd-$z$6 and even-$z$7 subsystems, yielding standing modes denoted $z$8, $z$9, and $x$0 (Fuller et al., 2024). In the EL CMi analysis, the first observationally detected quadrupole Tidally Tilted Standing mode is identified as $x$1, described as a superposition of $x$2 and $x$3 that behaves like an $x$4 pattern wrapped around the $x$5-axis, with maxima and minima offset from the tidal axis by an angle of $x$6, and as a standing mode that does not propagate around the $x$7-axis (Handler et al., 28 Jul 2025).

The theory also predicts a clear distinction between p-modes and g-modes. Pressure modes should exhibit triaxial behavior in stellar binaries close enough to be tidally synchronized, whereas gravity modes should remain aligned with the star’s spin axis because Coriolis effects dominate over tidal coupling for g-mode cavities (Fuller et al., 2024). This is a central selection rule for identifying genuine triaxial pulsation.

3. Spectral signatures and mode identification

The defining observational signatures of stellar triaxial pulsators are not merely nonradial modes in a close binary, but very specific orbital multiplets and orbital-phase modulation patterns. For dipole p-modes aligned with the orbital plane, the observed amplitudes and phases are modulated throughout the orbit, producing doublets in the power spectrum that are spaced by exactly twice the orbital frequency (Fuller et al., 2024). In the idealized dipole case, $x$8 and $x$9 each produce equal-amplitude doublets at $z$0, while $z$1 appears as a singlet with no orbital amplitude modulation (Jayaraman et al., 2024, Fuller et al., 2024).

The orbital-phase behavior distinguishes the $z$2 and $z$3 families. For $z$4 modes, amplitude maxima occur at eclipses and minima at quadratures, with $z$5 phase jumps at the minima; for $z$6 modes, the amplitude maxima occur at quadratures and minima at eclipses, again with $z$7 phase jumps at the minima (Zhang et al., 2023, Jayaraman et al., 2024). These two families therefore differ by a phase shift of $z$8 in orbital phase (Jayaraman et al., 2024).

A key diagnostic concerns exclusion of $z$9-type interpretations. In TIC 435850195, the observed doublets with no detectable central peak are inconsistent with simulated oblique $y$0 modes at the measured inclination, because those should produce a triplet whose central peak amplitude is roughly half that of the sidelobes at $y$1 (Jayaraman et al., 2024). The absence of such central peaks strongly supports the $y$2 interpretation.

Quadrupole modes introduce additional patterns. In the general theory, $y$3 behave like the in-plane dipoles and produce doublets at $y$4, while $y$5 produce doublets at $y$6 with four phase jumps per orbit, and $y$7 remains a singlet (Fuller et al., 2024). In more realistic models, slight mixing can produce triplets or asymmetric sidebands, but the harmonic content remains diagnostic (Fuller et al., 2024). EL CMi shows a broad multiplet around $y$8 with components at $y$9 and, in the full-eclipse data, weak components up to $Y_{\ell m}$0; the amplitude–phase behavior identifies it as the quadrupole $Y_{\ell m}$1 mode (Handler et al., 28 Jul 2025).

The methodology used across the current literature is consistent. TESS light curves are detrended or stripped of orbital harmonics, multifrequency fits are carried out, and échelle diagrams are constructed modulo the orbital frequency to reveal vertically aligned multiplets (Zhang et al., 2023, Jayaraman et al., 2024, Handler et al., 28 Jul 2025). Amplitude and phase are then reconstructed as functions of orbital phase for each multiplet and compared against theoretical templates (Zhang et al., 2023, Jayaraman et al., 2024, Handler et al., 28 Jul 2025). This combination of orbital-frequency spacing and orbital-phase modulation is the practical basis for triaxial mode identification.

4. Benchmark stellar systems

The current observational class is defined by a small number of close-binary $Y_{\ell m}$2 Sct systems observed with TESS.

System	Main claim	Key signatures
TIC 184743498	First tri-axial stellar pulsator	Five $Y_{\ell m}$3-like doublets, four $Y_{\ell m}$4-like doublets, two nearly unmodulated $Y_{\ell m}$5-like singlets (Zhang et al., 2023)
TIC 435850195	Second tri-axial pulsator	Fourteen dipole doublets separated by $Y_{\ell m}$6, eight identified as $Y_{\ell m}$7 and six as $Y_{\ell m}$8, plus two singlets consistent with $Y_{\ell m}$9 (Jayaraman et al., 2024)
EL CMi	Confirmation of triaxial pulsation theory	Two orthogonal dipole TTS modes and the first detected quadrupole TTS mode $\ell=1$0 (Handler et al., 28 Jul 2025)

TIC 184743498, an eclipsing binary with $\ell=1$1, was presented as the first clear example of a tri-axial pulsator (Zhang et al., 2023). Eleven of its prominent pulsation peaks form doublets split by $\ell=1$2 around inferred central frequencies, while two appear as singlets (Zhang et al., 2023). Five modes show amplitude maxima near the eclipses and minima near the ellipsoidal-light-variation maxima, with $\ell=1$3-rad phase jumps near quadrature, and were identified as $\ell=1$4 modes. Four modes show the complementary pattern, with amplitude minima at eclipses and maxima at quadrature, and were identified as $\ell=1$5 modes. Two singlets show no significant orbital-phase amplitude or phase modulation and were interpreted as $\ell=1$6 modes (Zhang et al., 2023). This tripartite structure established the phenomenological definition of the class.

TIC 435850195, an eclipsing binary with $\ell=1$7 d, was reported as the second-ever discovered tri-axial pulsator (Jayaraman et al., 2024). The TESS 200 s light curve revealed sixteen robustly detected pulsation multiplets, of which fourteen are dipole doublets whose two components are separated by $\ell=1$8 to within an rms fractional uncertainty $\ell=1$9 (Jayaraman et al., 2024). Eight of these were associated with $Y_{10x}$0 modes and six with $Y_{10x}$1 modes, while two strong singlets were interpreted as $Y_{10x}$2-like modes (Jayaraman et al., 2024). The system provided a particularly clean case for ruling out $Y_{10x}$3 alternatives and for showing the coexistence of all three dipole-axis families in one star.

EL CMi extended the class beyond dipole-only identification (Handler et al., 28 Jul 2025). It is an eclipsing close binary with orbital ephemeris

$Y_{10x}$4

containing a $Y_{10x}$5 Scuti primary and a donor star close to or filling its Roche lobe (Handler et al., 28 Jul 2025). The three dominant pulsation modes are centered near $Y_{10x}$6, $Y_{10x}$7, and $Y_{10x}$8, with orbital sidelobes that reveal the multiplet structure (Handler et al., 28 Jul 2025). The amplitude–phase behavior of $Y_{10x}$9 is consistent with a $x$00 dipole mode, that of $x$01 with a $x$02 dipole mode, and that of $x$03 with the quadrupole $x$04 mode, yielding the first detection of a quadrupole TTS oscillation mode in a triaxial pulsator (Handler et al., 28 Jul 2025).

A broader theoretical synthesis argues that previously discovered tidally tilted or “single-sided” pulsators are part of the same general phenomenon (Fuller et al., 2024, Handler et al., 28 Jul 2025). This suggests that the currently named tri-axial pulsators are the most explicit members of a wider population of tidally distorted close-binary pulsators.

5. Relation to roAp stars and the question of non-axisymmetry

Although roAp stars are not usually classified as tri-axial pulsators, they are relevant because they show how departures from a single symmetry axis manifest in multiplets, unequal sidelobes, and amplitude–phase modulation. Rapidly oscillating Ap stars are cool, chemically peculiar A-type stars with strong, organized magnetic fields and high-overtone, low-degree, non-radial p-modes with periods of 4.7–25.8 minutes (Zhong et al., 2024). In the oblique pulsator model, a central pulsation frequency $x$05 is split into components $x$06, where $x$07, and the geometry is characterized by the inclination $x$08 and magnetic obliquity $x$09 (Zhong et al., 2024).

In the TESS study of TIC 96315731, TIC 72392575, and TIC 318007796, quintuplets were analyzed by fitting rotational sidelobes and using amplitude ratios to infer $x$10 and $x$11 (Zhong et al., 2024). The paper uses the standard quadrupole relations

$x$12

together with spherical-harmonic decomposition into $x$13 components (Zhong et al., 2024). The resulting mode identifications were: distorted dipole for TIC 96315731, distorted quadrupole for TIC 72392575, and distorted dipole for TIC 318007796 (Zhong et al., 2024).

The relevance to triaxial pulsators lies in the departure from an ideal single-axis picture. The roAp paper explicitly notes that “triaxial pulsator” is not a standard term in the roAp literature, but proposes it as a heuristic for stars where the pulsation axis is not aligned with the magnetic axis, the mode is not axisymmetric in any simple frame, and the observed pattern cannot be described by a single symmetry axis (Zhong et al., 2024). However, the analysis remains conservative: it keeps the formal oblique pulsator model axisymmetric and does not fit non-axisymmetric $x$14 with $x$15 (Zhong et al., 2024).

This comparison is important because it delineates two related but distinct meanings of triaxiality in stellar pulsation. In tidally distorted close binaries, tri-axial pulsation is an explicit, coordinate-based classification of orthogonal pulsation axes tied to the stellar figure (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024). In roAp stars, “triaxial” is at most a heuristic label for increasingly complex multi-axis geometry and non-axisymmetric distortion (Zhong et al., 2024). The two literatures share the idea that multiplets alone do not uniquely identify pure $x$16 geometry and that mixed spherical-harmonic content is common (Zhong et al., 2024, Fuller et al., 2024).

6. Triaxial pulsars and gravitational-wave emission

Outside optical asteroseismology, the term “triaxial pulsar” has a distinct and older meaning: a rotating neutron star with unequal principal moments of inertia that emits continuous gravitational waves (Pathak et al., 2022). In this framework, the star is assumed to be a triaxial rigid rotator about a principal moment of inertia axis, with

$x$17

and an equatorial ellipticity

$x$18

so that the mass quadrupole varies periodically and radiates at

$x$19

for a star spinning at frequency $x$20 (Pathak et al., 2022).

The continuous-wave spin-down limit under this assumption is

$x$21

where $x$22 is the distance and $x$23 is the intrinsic spin-frequency derivative (Pathak et al., 2022). A recent study improved the inferred spin-down limits for 237 pulsars targeted by LVK by using the GalDynPsr package and a full Galactic potential model to correct for kinematic and dynamical effects in $x$24 and hence in $x$25 (Pathak et al., 2022). After exclusions, the clean sample contained 139 pulsars, and 93 of 139 showed an increase in $x$26 relative to Abbott et al. (2021), leading to larger $x$27 values (Pathak et al., 2022). The paper’s focus is not pulsation geometry in the asteroseismic sense, but the gravitational-wave implications of triaxial rotation.

A further extension concerns precession of triaxially deformed neutron stars (Gao et al., 2020). There the star is modeled as a rigid triaxial body with principal moments $x$28, $x$29, and $x$30, and free precession is described analytically in terms of Jacobi elliptic functions. The deformation parameters are

$x$31

and the angular velocities are written as

$x$32

with precession frequency parameter

$x$33

and elliptic modulus

$x$34

for wobble angle $x$35 (Gao et al., 2020). The gravitational-wave signal from such a precessing triaxial neutron star contains lines at

$x$36

showing that “triaxial pulsator” in this domain refers to a precessing gravitational-wave source rather than to an optical pulsating star (Gao et al., 2020).

This dual usage is a common source of confusion. In continuous-wave neutron-star astronomy, “triaxial pulsar” concerns body shape and mass quadrupole emission (Pathak et al., 2022, Gao et al., 2020). In close-binary asteroseismology, “triaxial pulsator” concerns the orientation of stellar oscillation eigenfunctions relative to the three principal axes of a tidally distorted star (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024).

7. Broader meanings and current significance

A third, more formal usage appears in nuclear structure theory. In the Bohr Hamiltonian with a Davidson potential in $x$37 and a steep harmonic oscillator in $x$38 centered at $x$39, a “triaxial pulsator” is a nucleus whose collective wave function is centered on a triaxial equilibrium shape but still has vibrational quanta in $x$40 and in small oscillations of $x$41, together with rotational excitations (Yigitoglu et al., 2010). The Hamiltonian is

$x$42

with triaxiality introduced by a potential minimum near $x$43 (Yigitoglu et al., 2010). In this usage, the term belongs to quadrupole collective dynamics rather than to stellar pulsation.

The coexistence of these three literatures explains why “triaxial pulsators” is not a universally standardized term. The most rapidly evolving and observationally active usage is the close-binary stellar one (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024, Handler et al., 28 Jul 2025). In that domain, several points are now established.

First, tidally distorted stars are generically triaxial pulsators for low-degree p-modes in sufficiently close synchronized binaries (Fuller et al., 2024). This is not a special-case anomaly but a geometric consequence of the triaxial equilibrium figure and of strong $x$44-coupling within a multiplet.

Second, the class is observationally identifiable. Equal-amplitude doublets at $x$45, specific orbital-phase amplitude and phase modulation, and the coexistence of in-plane and orbital-axis mode families provide an empirical diagnostic that is far more specific than ordinary rotational splitting (Zhang et al., 2023, Jayaraman et al., 2024, Fuller et al., 2024).

Third, the framework has moved beyond dipole-only phenomenology. EL CMi confirms the quadrupole TTS prediction through a detected $x$46 mode, showing that the theory extends naturally to higher multipoles (Handler et al., 28 Jul 2025).

Fourth, the class has direct asteroseismic promise. Because each radial order can, in principle, produce modes aligned with different principal axes, triaxial pulsators provide a new form of three-dimensional mode labeling. This suggests a route to detailed asteroseismic analyses of tidally tilted pulsators and, more broadly, to asteroseismic inferences of the structure of stars in close binaries before and after mass transfer and in three spatial dimensions (Fuller et al., 2024, Handler et al., 28 Jul 2025).

A plausible implication is that the present named examples represent only the easily recognized portion of a wider population of tidally distorted p-mode pulsators. The current literature already argues that several earlier tidally tilted or “single-sided” systems can be reinterpreted within the same general triaxial framework (Fuller et al., 2024, Handler et al., 28 Jul 2025). As TESS-like photometry, orbital-phase-resolved modeling, and forward seismic calculations improve, the term “triaxial pulsator” is likely to denote an increasingly broad and theoretically unified class of close-binary oscillators.