3 Main Uses:
- Fundamental Properties
- Dynamic and Activity
- Internal Structure
Joel Ong
Hubble Fellow,
Univ. of Hawaiʻi at Mānoa
Center for Astrophysics | Harvard &
Smithsonian;
January 28, 2025
physics of stellar interiors
quantitative
astronomy & astrophysics
(RHD simulations courtesy of Joel D. Tanner)
Convection excites
pressure waves (p-modes).

\(\ell = 0\) MDI Doppler velocities
Power spectra of MDI dopplergrams
\[ \begin{aligned} {\Delta\nu_\odot} &\sim 135\ \mathrm{\mu Hz} \\ {\nu_{\text{max},\odot}} &\sim 3090\ \mathrm{\mu Hz} \end{aligned} \]
(roughly 5-minute oscillations)
p-mode frequencies satisfy \(\nu_{n\ell} \sim \Delta\nu\left(n + {\ell \over 2} + \epsilon_\ell(\nu)\right) + \mathcal{O}(1/\nu)\)
Stochastic,
broad-band
excitation
\[\vdots\]
Telescopes can only point at one star at time…


Required photometric stability not achievable from ground
\[ \begin{aligned} {\Delta\nu} & \sim 1/t_\text{cross} \sim \sqrt{M/R^3}\\ {\nu_{\text{max}}} &\sim{g/c_s} \sim {M/R^2\sqrt{T_\text{eff}}} \end{aligned} \]
\[V_\text{osc} \sim L / M\]
\[ \begin{aligned} {\Delta\nu} &\sim \sqrt{M/R^3}\\ {\nu_{\text{max}}} &\sim {M/R^2\sqrt{T_\text{eff}}} \end{aligned} \]
\[ \begin{aligned} {M \over M_\odot} &\sim \left(\nu \over \nu_{\text{max},\odot}\right)^{3}\left(\Delta\nu \over \Delta\nu_\odot\right)^{-4} \left(T_\text{eff} \over T_{\text{eff},\odot}\right)^{3/2} \\ {R \over R_\odot} &\sim \left(\nu \over \nu_{\text{max},\odot}\right)\left(\Delta\nu \over \Delta\nu_\odot\right)^{-2} \left(T_\text{eff} \over T_{\text{eff},\odot}\right)^{1/2} \end{aligned} \]
(better scaling relations: Ong & Basu 2019a,b)
Hare “Zebedee”, Cunha+ incl. Ong (2021);
HPC & pipeline: Ong+ 2021a

Precise measurements of field stars: \[ {\sigma_R \over R} \lesssim 2 \%;\ {\sigma_M \over M} \lesssim 5 \%;\ \sigma_\text{Age} \lesssim 0.4\ \mathrm{Gyr} \]
\[\begin{array}{c} {\scriptsize\text{Power}}\\ \big\uparrow \end{array}\]
\[\longrightarrow {\scriptsize\text{Frequency}}\]
Rotational inversions constrain differential rotation
(e.g. Backus & Gilbert 1968; Gough 1985;
Pijpers & Thompson 1992; Schunker 2016;
Ong 2024 etc.)
(e.g. Ong et al. 2024, Ong 2025)
Probes of the internal states of stars … now return constraints on stellar structure previously only theorized. —Astro 2020
Decadal Survey
(e.g. Bellinger+ 2017,2019;Ong & Basu
2019a,b;
Lindsay, Ong, Basu 2022†, 2023†,
2024†;
Vanlaer+ 2023; Buchele+ 2024, 2025)
†: mentoree paper

Bellinger+ 2019
Relative difference in isothermal sound speed
Asteroseismology is a low-investment
yet versatile conceptual instrument
that has revolutionised stellar metrology,
both global and local.
We are drowning in data.
big data deluge
physical interpretation
analytic theory
computational technique
statistical methodology
\[ {\large\textbf{Theory}} \longrightarrow \text{Single Stars} \longrightarrow \text{Populations} \longrightarrow \text{Astrophysics} \]
Technique
Development
Interpretation
Data Volume
Evolved stars dominate our asteroseismic sample.
(e.g. only \(\sim 100\) Kepler main-sequence stars)
\[\large V_\text{osc} \sim L/M\]
\(T_\text{eff}/\mathrm{K}\)
\(R/R_\odot\)
Pressure waves (p-modes)
propagate isotropically.
p-modes:
Characteristic overtone frequency spacing \(\Delta\nu\)
Buoyancy waves (g-modes)
propagate anisotropically.
g-modes:
Characteristic overtone period spacing \(\Delta\Pi_\ell\)
g-mode Period Spacing \(\Delta\Pi_1/\mathrm{s}\)
p-mode Frequency Spacing \(\Delta\nu/\mu\mathrm{Hz}\)
from Mosser+ (2014)
Single-star electron degeneracy sequence:
deviations → merger remnants?
(Rui and Fuller 2021, Deheuvels+ 2021)
clump stars
first-ascent RGs
\[\begin{array}{c} {\scriptsize\text{Power}}\\ \big\uparrow \end{array}\]
\[\longrightarrow {\scriptsize\text{Frequency}}\]
(proxy for age \(\to\))
Mixed modes exhibit
avoided crossings
between underlying p- and
g-modes.
\[\huge {\partial^2 \over \partial t^2}f(t, \mathbf{x}) = c_s^2 \nabla^2 f(t, \mathbf{x})\]
\[\huge -\omega^2 \hat{f}(\omega, \mathbf{k}) = -c_s^2 |\mathbf{k}|^2 \hat{f}(\omega, \mathbf{k})\]
Dispersion relation: \[\boxed{\huge \omega^2 = c_s^2 |\mathbf{k}|^2}\]
\[\huge {\partial^2 \over \partial t^2}f = \left(c_s^2 \nabla^2 - \omega^2_\text{ac}(\mathbf{x})\right) f\]
\[\huge \underbrace{-\nabla^2 f}_{\equiv k^2f} = -\left(\omega^2 - \omega^2_\text{ac}(\mathbf{x}) \over c_s^2\right) f\]
Wave propagates where \(k^2(r,
\omega) > 0\),
and decays where \(k^2(r, \omega)
< 0\).
\[{c_s^2 k_r^2 \sim \omega^2 \left(1 - {{\color{blue} S_\ell}^2 \over \omega^2}\right)\left(1 - {{\color{darkorange}N}^2 \over \omega^2}\right)}\]
\[\small N^2 = {- g}\left.{\partial \log \rho \over \partial s}\right|_P{\mathrm d s \over \mathrm d r}\] entropy gradient (\(=0\) in CZ)
\[\small S_\ell^2 = c_s^2 k_h^2 = {\ell(\ell+1) c_s^2 \over r^2}\] wave angular momentum
\[\Large \omega_-^2 \sim N^2 {k_h^2 \over |\mathbf{k}|^2}\]
\[{\color{red} \omega_g < N, S_\ell}\]
\[\Large \omega_+^2 \sim c_s^2 |\mathbf{k}|^2\]
\[{\color{gray} \omega_p > S_\ell, N}\]
Brute-force numerical solution (quantitative)
vs
JWKB approximation (qualitative)
\[\Large \hat H \left|n\right> = \left(\hat T + \hat V\right)\left|n\right> = E_n \left|n\right>\]
\[\large {\color{blue}{\hat H_1}} \left|n\right> = \left(\hat T + \hat V_1\right)\left|n\right> = {\color{blue}E_n \left|n\right>}\]
\[\large {\color{orange}{\hat H_2}} \left|n\right> = \left(\hat T + \hat V_2\right)\left|n\right> = {\color{orange}E_n \left|n\right>}\]
\[\huge \hat H = \hat T + \hat V_1 + \hat V_2 = {\color{blue}{\hat H_1}} + \hat V_2 = \hat V_1 + {\color{orange}{\hat H_2}}\]
\[\left|\psi_\text{mol}\right> = {\color{blue}\sum_i c_{1,i}\left|\psi_{1,i}\right>} + {\color{orange}\sum_j c_{2,j}\left|\psi_{2,j}\right>}.\]
Energy eigenvalues and mixing coefficients are solutions to the
Generalised Hermitian Eigenvalue Problem
\[ \mathbf{H} \begin{bmatrix} {\color{blue}\mathbf{c}_1} \\ {\color{orange}\mathbf{c}_2} \end{bmatrix} = E \cdot \mathbf{D} \begin{bmatrix} {\color{blue}\mathbf{c}_1} \\ {\color{orange}\mathbf{c}_2} \end{bmatrix} \] \[ \small H_{{\color{blue}i}{\color{orange}j}} = \langle {\color{blue}\psi_i}|\hat{H}|{\color{orange}\psi_j}\rangle;~~~~D_{{\color{blue}i}{\color{orange}j}} = \langle {\color{blue}\psi_i}|{\color{orange}\psi_j}\rangle \]
\[\Large \hat{\mathcal{L}} \xi_n = -\omega_n^2 \xi_n\]
\[\Large \hat{\mathcal{L}} = {\color{red}\hat{\mathcal{L}}_\gamma} + \hat{\mathcal{R}}_\gamma\]
\[\Large {\color{red}\hat{\mathcal{L}}_\gamma \xi_{\gamma,n} = -\omega^2_{\gamma,n}\xi_{\gamma,n}}\]
\[\Large \hat{\mathcal{L}} = \hat{\mathcal{R}}_\pi + {\color{grey}\hat{\mathcal{L}}_\pi}\]
\[\Large {\color{grey}\hat{\mathcal{L}}_\pi \xi_{\pi,n} = -\omega^2_{\pi,n}\xi_{\pi,n}}\]
\[\Large \hat{\mathcal{L}} = {\color{grey}\hat{\mathcal{L}}_\pi} + \hat{\mathcal{R}}_\pi = \hat{\mathcal{R}}_\gamma + {\color{red}\hat{\mathcal{L}}_\gamma}\]
\[\Large \xi_\text{mixed} \sim {\color{grey}c_\pi \xi_\pi} + {\color{red}c_\gamma \xi_\gamma}\]
\[\xi_\text{mixed} \sim {\color{grey} \sum_i c_{\pi, i} \xi_{\pi,i}} + {\color{red} \sum_j c_{\gamma, j} \xi_{\gamma,j}}\]
Mixed mode
frequencies and mixing coefficients can likewise be found by solving
the
Generalised Hermitian Eigenvalue Problem: \[
\scriptsize
\begin{bmatrix}
{\color{grey}\mathbf{L}_{\pi\pi}} & \mathbf{L}_{\pi\gamma} \\
\mathbf{L}_{\pi\gamma}^T & {\color{red}\mathbf{L}_{\gamma\gamma}}
\end{bmatrix}
\begin{bmatrix}
{\color{grey}\mathbf{c}_\pi} \\ {\color{red}\mathbf{c}_\gamma}
\end{bmatrix}
= -\omega^2 \begin{bmatrix}
\mathbb{1} & \mathbf{D} \\ \mathbf{D}^T & \mathbb{1}
\end{bmatrix}
\begin{bmatrix}
{\color{grey}\mathbf{c}_\pi} \\ {\color{red}\mathbf{c}_\gamma}
\end{bmatrix}.
\] \[
\small
L_{{\color{grey}i}{\color{red}j}} = \langle {\color{grey}\xi_i},
\hat{\mathcal{L}}{\color{red}\xi_j}\rangle;~~~~D_{{\color{grey}i}{\color{red}j}}
= \langle {\color{grey}\xi_i}, {\color{red}\xi_j}\rangle
\]
(Ong & Basu 2020; Ong & Gehan 2023; Ong 2025)
State of the art for determining (sub)giant
structure and properties
(Ong+ 2021a, b, c),
and internal rotation (Ong+ 2022;
Ong 2024, 2025).
We know more about giant cores
than about the core of our own Sun!
(e.g. Ong et al. 2024, Ong 2025)
\[\begin{array}{c} {\scriptsize\text{Power}}\\ \big\uparrow \end{array}\]
\[\longrightarrow {\scriptsize\text{Frequency}}\]
Li et al. Nat. 2022:
Asymmetric splittings probe
core magnetic fields
(Li+ 2023, Deheuvels+ 2023;
Rui, Ong, Mathis, 2024†)
Population studies
of rotation vs. magnetism
(Hatt, Ong, et al., 2024†)

\[\scriptsize \delta \nu_{\text{mag}, g, \ell=1} \sim {m^2 \over \nu^3}\]


†: mentoree paper
Saunders et
al. 2024
for \(\ell = 1,\)
\[\small\mathbf{J}_x \hat{=} {1 \over \sqrt{2}}\begin{bmatrix}0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0\end{bmatrix}; \mathbf{J}_y \hat{=} {1 \over \sqrt{2}}\begin{bmatrix}0 & i & 0 \\ -i & 0 & i \\ 0 & -i & 0\end{bmatrix}; \mathbf{J}_z \hat{=} \begin{bmatrix}-1 & 0 & 0 \\ 0 &0 &0 \\ 0& 0 &1\end{bmatrix}.\]
\[\left(-\mathbf{\Omega}_0^2 + 2 \omega m \mathbf{R} + \omega^2 \mathbf{\Delta}\right)\mathbf{c} = 0\]
For fixed \(n\) (to leading order):
\[(-\omega_0^2 \mathbb{1}_{2\ell+1} + 2 \omega \mathbf{J}_z R_{n,n} + \omega^2 \mathbb{1}_{2\ell+1})\mathbf{y} = 0\]
\[\implies \left(-\mathbf{\Omega_0}^2 \otimes \mathbb{1}_{2\ell+1} + 2 \omega \underbrace{\mathbf{R} \otimes \mathbf{J}_z}_{\tilde{\mathbf{R}}} + \omega^2 \mathbf{\Delta} \otimes \mathbb{1}_{2\ell+1}\right)\mathbf{x} = 0\]
\[\small \begin{aligned} \mathbf{R}_{n\ell, n\ell} &= b_{n\ell}\int {\mathbf{d}^\ell}^\dagger(\beta(r)) \Omega(r) \mathbf{J}_z {\mathbf{d}^\ell}(\beta(r))\ K(r)\ \mathrm{d} r\\ &= b_{n\ell}\int \Omega(r) (\hat{\mathbf{n}} \cdot \vec{\mathbf{J}})\ K(r)\ \mathrm{d} r\\ &= \boxed{b_{n\ell}\left(\int \vec{\mathbf{\Omega}}(r) K(r)\ \mathrm{d} r\right) \cdot \vec{\mathbf{J}}}. \end{aligned} \]
\(\implies\) For each mode, AM
matrix is
specified by usual vector addition.
\[ \small \left(\begin{bmatrix} {\color{grey}\mathbf{L}_{\pi\pi}} & \mathbf{L}_{\pi\gamma} \\ \mathbf{L}_{\pi\gamma}^T & {\color{red}\mathbf{L}_{\gamma\gamma}} \end{bmatrix} \otimes \mathbb{1}_{2\ell+1} + 2 \omega \begin{bmatrix} {\color{forestgreen}\tilde{\mathbf{R}}_\pi} & 0 \\ 0 & {\color{forestgreen}\tilde{\mathbf{R}}_\gamma}\end{bmatrix} + \omega^2 \begin{bmatrix} \mathbb{1} & \mathbf{D} \\ \mathbf{D}^T & \mathbb{1} \end{bmatrix} \otimes \mathbb{1}_{2\ell+1} \right)\mathbf{x} = 0. \]
Mode visibilities are specified by \(\mathbf{x}^\dagger(\mathbb{1} \otimes \mathbf{P})\mathbf{x}\), where \(\mathbf{P}\) is the projection matrix onto \(m=0\) in the observer’s coordinate frame.
\[ \vdots \]
New Theory
(Ong 2025)
\[\begin{array}{c} {\scriptsize\text{Power}}\\ \big\uparrow \end{array}\]
\[\longrightarrow {\scriptsize\text{Frequency}}\]
Ong 2025

From Huber et al. 2013
Ong 2025
Probing the star-planet connection, and non-canonical evolution
e.g. Ong 2025;
Ong et al. 2024;
Hon et al. (incl Ong) Nat. 2023
(also incl. Ong: Huber+ 2019, 2022; …)
Kepler-56’s core and envelope
rotate around different axes.
(Ong 2025)
Zvrk rotates too fast to
not have eaten something recently.
(Ong et al. 2024)
8 UMi b should have
been consumed, but wasn’t?
(Hon+ 2023 incl. Ong)
In the data-rich régime,
technique development and
theoretical interpretation are
rate-limiting steps to discovery.

\[ {\large\color[RGB]{0, 120, 0}\textbf{Theory} \longrightarrow} \mathrm{Single~Stars} \longrightarrow \text{Populations} {\color[RGB]{0, 120, 0}\longrightarrow} \large\text{Astrophysics} \]
Technique
Development
Interpretation
Data Volume
Statistical sample (Ong: WP 120, 128)
\(\sim 10^6\) red giants in galactic bulge & globular clusters
Long temporal baselines
(Ong: TASOC WG1, WG2, WG7)
\(\sim100 \to 10,000++\)
stars:
How do we cope with
a deluge of new data?
e.g. Hey, Huber, Ong, et al. in review;
Nielsen, Ong, et al., in review.
(incl. Ong: Cunha+ 2021; Nielsen+
2021; Campante+ 2023)

Ong: TASOC WG1, WG2; PLATO WP120, 128
?
Combining with Photometric surveys —
e.g. Ong et al. (2024); Hart+
incl. Ong (2023) 
Seismology from Extreme Precision Radial Velocities
(Li, Huber, Ong, et al. 2025; Hon+ incl. Ong 2024)
Ong: ASAS-SN; SONG WG1 & WG2; Keck Planet Finder via CPS



Coming Soon: HARPS-N; MINERVA; Rubin/LSST; G-CLEF?

Probes of mixing processes
near convective boundaries
(Lindsay, Ong, Basu 2022†, and
2024†;
Ong, Lindsay, Reyes et al. 2025;
Reyes, Stello, Ong, et al. 2025†, Nature)


†: mentoree paper
Ong: TESS Guest Investigator Cycle 7, PI
?
from Blouin et al. 2023a, b
Nonlinear evolution is
the primary obstacle to
g-mode inversions.
Will understanding this (Hoogendam, Ong, in prep.†) permit further technique development?

proxy for age \(\to\)
†: mentoree paper

(RHD simulations courtesy of Joel D. Tanner)
?
Stellar oscillations have revolutionised our
understanding of
stellar structure, dynamics, and
evolution.
The development of
theory and technique
unlocks their use as a tool (e.g.: mixed modes
\(\to\) RGB science).
My group will combine
these developments with new surveys
to make these measurements astrophysically
interesting.
The Sun as seen by SOHO:
Bedding & Kjeldsen 2005
Procyon from MOST vs. RVs:
Huber et al. 2011




Normal modes satisfy an operator eigenvalue problem
\(\hat{\mathcal{L}}\left|\xi_{n\ell
m}\right> = -\omega^2_{n\ell m}\left|\xi_{n\ell
m}\right>\).
Under a small perturbation \(\hat{\mathcal{L}} \mapsto \hat{\mathcal{L}} + \lambda \hat{\mathcal{V}}\),
\[\begin{aligned} \omega^2_{n\ell m} &\mapsto \omega^2_{n\ell m} - \lambda \left<\xi_{n\ell m}|\hat{\mathcal{V}}|\xi_{n\ell m}\right> + \mathcal{O}\left(\lambda^2\right)\\ &\approx \omega^2_{n\ell m} - \lambda \int \mathrm d^3 x \cdot \rho \left(\vec{\xi}_{n\ell m}^* \hat{\mathcal{V}} \vec{\xi}_{n\ell m}\right) + \mathcal{O}\left(\lambda^2\right). \end{aligned}\]
e.g. for (slow) rotation in particular
\[ \omega^2 \mapsto \omega^2 + 2 m\ b_{n\ell m} \omega \int \Omega(r) K_{nlm}(r) \mathrm dr \]











For slow rotation, \[\boxed{\delta\omega_{nlm} \sim m b_{nl} \int \Omega(r) K_{nl}(r) \mathrm d r}\]
Variations in \(\delta\omega_\text{rot}\) are
radial differential rotation.
\[\boxed{{\color{orange}{\delta\omega_{nlm}}} \sim {\color[RGB]{0,100,255}{m b_{nl} \sum_i {\color{black}{\Omega(r_i)}} K_{nl}(r_i) \Delta r_i}}}\] is of the form \({\color[RGB]{0,100,255}\mathbf{A}}\mathbf{x} = {\color{orange}\mathbf{b}}\).
i.e. Inferring the rotational profile
\(\Omega(r)\) is a linear
inverse problem.
Rotational Inversions
\[\scriptsize\begin{aligned} \int \Omega(r) {\color[RGB]{0,100,255}\left(\sum_i c_i K_{i}(r)\right)} \mathrm d r &\sim \sum_{i} \left(c_i \over m_i b_{i} \right) \delta\omega_{\mathrm{rot}, i} \\ & \to \boxed{\Omega({\color[RGB]{0,100,255}r_0})}\end{aligned}\]
e.g. OLA
(Backus & Gilbert 1968; Gough 1985;
Pijpers & Thompson 1992; Schunker 2016; Ong 2024;
etc.)
\(\implies\) we can make measurements of internal rotational structure!
two standard
solar models
\[\huge {{\color{orange}{\delta\omega_i \over \omega_i}} \approx {\color[RGB]{0,100,255}{\int K_{\rho|c_s^2, i}(r)}} {\delta\rho \over \rho}{\color[RGB]{0,100,255}{\mathrm d r}} + {\color[RGB]{0,100,255}{\int K_{c_s^2|\rho, i}(r)}} {\delta c_s^2 \over c_s^2}{\color[RGB]{0,100,255}{\mathrm d r}}} \]
\[{\Large{\color[RGB]{0,100,255}\mathbf{A}}\mathbf{x} = {\color{orange}\mathbf{b}}}\]
Structure perturbations
do not reorder mixed modes!
(cf. Ball+ 2018)
Mode mixing yields avoided crossings
between multiplet components of identical \(m\)
(cf. Mosser+ 2012, Ouazzani+ 2013, Deheuvels+ 2017)
pure p
pure g
mixed

Very high \(\mathrm{A(Li)}\),
but otherwise innocuous
\[\ell = 0,2?\]
But Kepler says \(\ell =
0\) have to live here!
(and theory says so too…)



Rotational splittings from seismology:
\(\implies P_\text{rot} \sim 115\ \mathrm{d}?\)

Rotational signal probably got
detrended away
by systematic corrections…
Very suggestive…
but is this real?
Probably yes!
Asteroseismology:
Photometry:
Spectroscopy:
Gaia DR3:
Gaia RV scatter rules out
large RV semiamplitudes…
\(P_\text{orb} \gg 99\ \mathrm{d}\)
cannot spin star
up to 99-day rotational period…
Remaining permissible orbits are
unstable to tidal dissipation!
\(\implies\) ENGULFMENT?
\(\forall \ell, \exists (2\ell + 1) \times
(2\ell + 1)\) matrices
\(\mathbf{J}_x^\ell\), \(\mathbf{J}_y^\ell\), \(\mathbf{J}_z^\ell\) satisfying commutation
relations
\(\left[\mathbf{J}_i, \mathbf{J}_j\right] =
-i\epsilon_{ijk}\mathbf{J}_k\), with
\(\mathbf{J}_z \hat{=} \mathrm{diag}(-\ell,
-\ell + 1 \ldots \ell - 1, \ell).\)
\[\small\mathbf{J}_x \hat{=} {1 \over \sqrt{2}}\begin{bmatrix}0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0\end{bmatrix}; \mathbf{J}_y \hat{=} {1 \over \sqrt{2}}\begin{bmatrix}0 & i & 0 \\ -i & 0 & i \\ 0 & -i & 0\end{bmatrix}; \mathbf{J}_z \hat{=} \begin{bmatrix}-1 & 0 & 0 \\ 0 &0 &0 \\ 0& 0 &1\end{bmatrix}.\]
For fixed \(m\) (to leading order):
\[\left(-\mathbf{\Omega}_0^2 + 2 \omega m \mathbf{R} + \omega^2 \mathbf{\Delta}\right)\mathbf{c} = 0\]
For fixed \(n\) (to leading order):
\[(-\omega_0^2 \mathbb{1}_{2\ell+1} + 2 \omega \mathbf{J}_z R_{n,n} + \omega^2 \mathbb{1}_{2\ell+1})\mathbf{y} = 0\]
\[\vec{\xi}_{n\ell m}(r, \theta, \varphi) = \vec{\xi}_{n\ell}(r) Y_\ell^m(\theta, \phi) \iff \underbrace{\tilde{R}_{n, n', m, m'} = R_{n,n'} J_{m,m'}}_{\text{this is a }\textbf{tensor product}!}\]
\[\implies \left(-\mathbf{\Omega_0}^2 \otimes \mathbb{1}_{2\ell+1} + 2 \omega \underbrace{\mathbf{R} \otimes \mathbf{J}_z}_{\tilde{\mathbf{R}}} + \omega^2 \mathbf{\Delta} \otimes \mathbb{1}_{2\ell+1}\right)\mathbf{x} = 0\]
Let’s associate with each mass shell at \(r\) both \(\Omega(r)\)
(as is customary), and also an axis \(\hat{\mathbf{n}}(r) =\sum_i n_i
\mathbf{e}_i\).
\[\small \begin{aligned} \mathbf{R}_{n\ell, n\ell} &= b_{n\ell}\int {\mathbf{d}^\ell}^\dagger(\beta(r)) \Omega(r) \mathbf{J}_z {\mathbf{d}^\ell}(\beta(r))\ K(r)\ \mathrm{d} r\\ &= b_{n\ell}\int \Omega(r) (\hat{\mathbf{n}} \cdot \vec{\mathbf{J}})\ K(r)\ \mathrm{d} r\\ &= \boxed{b_{n\ell}\left(\int \vec{\mathbf{\Omega}}(r) K(r)\ \mathrm{d} r\right) \cdot \vec{\mathbf{J}}}. \end{aligned} \]
\[ \boxed{\tilde{\mathbf{R}} = \int \mathbf{K}\otimes\left(\Omega_\text{rot}(r)\ \hat{\mathbf{n}}\cdot \vec{\mathbf{J}}\right)\ \mathrm{d}r} \]
\(\implies\) For each mode, AM
matrix is
specified by usual vector addition.
We only assume that the pure p- and g-mode
solutions
are separately amenable to separation of
variables;
the mixed-mode eigenfunctions need not be.
\[ \small \left(\begin{bmatrix} {\color{grey}\mathbf{L}_{\pi\pi}} & \mathbf{L}_{\pi\gamma} \\ \mathbf{L}_{\pi\gamma}^T & {\color{red}\mathbf{L}_{\gamma\gamma}} \end{bmatrix} \otimes \mathbb{1}_{2\ell+1} + 2 \omega \begin{bmatrix} {\color{forestgreen}\tilde{\mathbf{R}}_\pi} & 0 \\ 0 & {\color{forestgreen}\tilde{\mathbf{R}}_\gamma}\end{bmatrix} + \omega^2 \begin{bmatrix} \mathbb{1} & \mathbf{D} \\ \mathbf{D}^T & \mathbb{1} \end{bmatrix} \otimes \mathbb{1}_{2\ell+1} \right)\mathbf{x} = 0. \]




\(\beta = 0\)
\(\beta = {\pi\over10}\)
\(\beta = {\pi\over2}\)
\(\beta = \pi\)
