Research

Solar-like oscillations lie at the heart of my research. These are normal modes driven by convective turbulence, associated with a characteristic excitation frequency \(\nu_\mathrm{max}\). We measure these normal-mode frequencies using time-series observations of stellar surfaces.

Time series observations of stellar surfaces (§1 — data shown are from \(\ell=0\) MDI solar dopplergrams) yield Fourier power spectra (§2) showing many peaks, each of which corresponds to a normal mode of oscillation. The normal modes of stars like the Sun show a characteristic overtone spacing \(\Delta\nu\), demonstrating that they are pressure waves, and an excitation frequency \(\nu_\text{max}\), a feature of convectively-driven pulsations. We may phase-fold the power spectrum by the overtone spacing \(\Delta\nu\) (§§3-5) to obtain so-called “échelle” diagrams — each ridge in such a diagram (§§6-7) corresponds to modes of a different spherical harmonic degree \(\ell\).

Today, asteroseismology not only provides our most precise available estimates for stellar masses, radii, and ages, but is also the only observational technique capable of directly inspecting the interior structure and dynamics of stars. These measurements have become foundational for our physical understanding of stars and planets — informing us regarding their internal physics, compositions, formation, and evolution — and for calibrating secondary methods for determining these properties (e.g. traditional spectroscopy, isochrone fitting, or gyrochronology). As a result of this usefulness, the 2020 Astronomy Decadal Survey has singled out continued investment in asteroseismology in particular as being essential for catering to the future needs of our field.

Theoretical Asteroseismology

In stars like the sun, these normal modes are pressure waves (p-modes). The shapes of these normal modes can be described by spherical harmonics, indexed by quantum numbers \(l\) and \(m\), in addition to their radial order \(n\). Observationally, the p-mode frequencies of stars like the Sun are such that modes of each \(\ell\) are separated by a roughly uniform overtone spacing \(\Delta\nu\): \[\nu_{n\ell m} = \Delta\nu\left(n + {\ell \over 2} + \epsilon_{\ell m}(\nu)\right),\] where \(\Delta\nu\) is called the “large frequency separation”. This is what we would expect for sound waves (“p-modes”, because pressure provides the restoring force) in a spherically-shaped resonant mode cavity, where a more or less constant overtone spacing emerges, given by the sound-crossing time. \(\epsilon_{\ell m}(\nu)\) is a function of frequency that specifies the eigenvalues for a given pair of \(\ell, m\), and containes information about how the star’s internal structure and dynamics differ from being a homogenous and static spherical mode cavity.

On the other hand, other stars possess normal-mode frequencies that satisfy a converse eigenvalue relation, \[{1 \over \nu_{n\ell m}} = \Delta\Pi_\ell\left(n + {\ell \over 2} + \epsilon_{\ell m}(\nu)\right),\] where it is the mode periods, rather than frequencies, which have a characteristic “undertone” period spacing relative to the fundamental mode. Waves of this kind are buoyancy waves (“g-modes”, because gravity provides the restoring force), not unlike those which you might see on the surface of a body of water, or reflected in the undulatory motions of clouds in the atmosphere.


Multicavity oscillations are well-approximated as linear combinations of pure p- and g-modes (above), in the same way that molecular orbitals are well-approximated as linear combinations of atomic orbitals (below).

Much of my research pertains to understanding, hopefully analytically, how these two kinds of normal modes relate to the internal structures of stars, react to each other, and respond to perturbations not captured in our computational modelling (e.g. as caused by global magnetic activity, radiative turbulence, or deficiencies in our descriptions of their internal structure). Recent highlights include:

Selected Publications:

  • Ong, J. M. J., & Basu, S. 2020, ApJ, 898, 127 — Hybridisation of standing waves
  • Ong, J. M. J., et al. 2021, ApJ, 920, 8 — Avoided crossings and perturbation theory
  • Ong, J. M. J., et al. 2022, ApJ, 940, 18 — Avoided crossings in rotating stars
  • Ong, J. M. J., & Gehan, C. 2023, ApJ, 946, 92 — Translating between coupling matrices and JWKB analysis

Asteroseismology as a Tool

Astronomers are often more interested in the inverse problem: that is, in inferring properties of a star given a set of observed oscillation frequencies. For example, we might infer masses, radii, and ages by comparing these mode frequencies against computational models of stellar structure. Alternatively, we might infer rotation rates and magnetic field strengths, if we are able to examine how modes of different \(m\) pulsate at different frequencies. Unfortunately, accurate measurement of the normal-mode frequencies requires relatively long observational campaigns. Often, the only quantities that are directly accessible are average values of \(\Delta\nu\), \(\epsilon\), and \(\nu_\text{max}\).

Échelle power diagram of Zvrk, a rapidly-rotating red giant, showing rotationally-split multiplets

Some of my recent work involves using these diagnostics of internal rotation to make interesting inferences about interactions between stars and their interactions with companions (e.g. planets). Highlights of this include a chemically anomalous, rapidly-rotating, planet-engulfment candidate, as well as a red giant whose core and envelope rotate separately around axes pointing in different directions.

Échelle diagram for ν Indi

For the past few years, I’ve also been working on results from the NASA Transiting Exoplanet Survey Satellite (TESS) mission through the TESS Asteroseismic Science Consortium (TASC) (in working groups 1 and 2), helping to constrain the global properties and interior structures of these stars in this manner. For this purpose, I am both a methodologist and a practitioner: I develop and investigate methods used in deriving these constraints, as well as actually perform the calculations for actual stars. I’m glad to be able to make contributions to such a large collaboration. Interesting stars I’ve studied include:

Selected Publications:

  • Ong, J. M. J., et al. 2024, ApJ, 966, 42 — A rapidly-rotating planet-engulfment candidate
  • Ong, J. M. J., ApJ. 980, 40 — Kepler-56, whose core and envelope rotate around different axes

Extreme Precision Radial Velocity (EPRV) Instrumentation

The tale of exoplanet detection and discovery is long and storied, but the oldest and still most reliable way of detecting planets, and measuring their masses, is through the radial velocity method, which requires spectroscopic characterisation. As planets orbit their stars, they exert an opposite gravitational pull on the star that causes it to also orbit (much more slowly) their common centre of mass, moving towards and then away from us over time. From Earth, we are able to measure this time-varying reflex velocity by the Doppler shift it induces on stellar spectra.

EXPRES cross-correlation functions for HD 75732 (55 Cancri) computed with respect to different Gauss-Hermite window functions

These measurements are made on high-resolution spectrographs. One of them, the EXtreme PREcision Spectrograph, was built at Yale. I wrote parts of its calibration and analysis software, including its template cross-correlation-function (CCF) radial velocity code, and various infrastructural components of its automated optimal-extraction data reduction pipeline.

While the ultimate goal of EXPRES is, of course, to discover long-period, low-mass planets with the radial-velocity method that would otherwise not be accessible to short-baseline space photometry, the high-resolution spectra returned from this process are very useful for other purposes. One example of this is Doppler imaging of active stellar surface, which can be further refined by way of joint constraints with contemporaneous photometry.

More recently, I’ve become a regular observer on the Keck Planet Finder instrument on Maunakea, Hawaiʻi, and have conducted some of the earliest science using it.

Comparison of power spectra from σ Draconis (HD 185144) as obtained from space photometry by the TESS satellite (faint blue, with a smoothed version shown in black), vs. from ground-based radial velocities measured using KPF. The frequency of maximum power \(\nu_\text{max}\) exhibited by the velocity power spectrum differs from that exhibited by photometry, shown with the gray line (with TESS modes identified with vertical dotted lines). I took these radial velocities personally!

Selected Publications:

  • Blackman, Ong, et al. 2019, AJ, 158, 40 — wavelength-dependent barycentric corrections
  • Petersburg, Ong, et al. 2020, AJ, 159, 187 — The EXPRES data reduction pipeline
  • Hon et al. incl. Ong 2024 — σ Dra, the first K dwarf with asteroseismology from KPF
  • Li, Huber, Ong et al. 2025, submitted to ApJ — HD 219134, a radius-inflated K dwarf