Pre-Main Sequence Solar-Like Oscillators

• Pre-Main Sequence solar-like oscillators are young, contracting stars with convective envelopes that show stochastic, low-amplitude p-mode oscillations similar to the Sun.
• Asteroseismic scaling relations and Bayesian detection methods quantify oscillation parameters, with νₘₐₓ values >500 μHz and Δν scaling tied to stellar density.
• High-precision observations from missions like PLATO enable robust constraints on internal structure and accretion histories, refining early stellar evolution models.

Pre-main sequence (PMS) solar-like oscillators are young, contracting stars—typically with convective envelopes—showing stochastically excited, low-amplitude p-mode oscillations analogous to those observed in the Sun and solar-type main-sequence stars. These oscillations encode diagnostics of the stars’ internal structure, mean density, and accretion history, providing access to stellar physics during the earliest, most active evolutionary phases. Direct detection has long been precluded by strong surface magnetic activity, accretion, and photometric noise, but recent advances in photometric analysis and asteroseismic modeling have identified robust candidates among T Tauri stars and have established the theoretical basis for high-precision asteroseismology in the PMS phase (Müllner et al., 2020, Jørgensen et al., 4 Jan 2026).

1. Asteroseismic Scaling Relations in PMS Stars

Solar-like oscillations in PMS stars peak at the frequency of maximum power, $\nu_\mathrm{max}$, determined by the star’s surface gravity and effective temperature. The canonical scaling relation is

$$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$

with $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$. In theoretical models of accreting PMS stars (0.7–1.6 $M_\odot$), $\nu_\mathrm{max}$ is always greater than $500\,\mu\mathrm{Hz}$, reflecting their compact, high-gravity structures—even at ages well before the zero-age main sequence (ZAMS). As contraction proceeds, $\nu_\mathrm{max}$ rises, reaching values up to $4000\,\mu\mathrm{Hz}$ near the ZAMS for solar-mass stars (Jørgensen et al., 4 Jan 2026).

The large frequency separation, $\Delta\nu$, scales as the square root of mean stellar density: $$\Delta\nu \propto \sqrt{\frac{M}{R^3}}, \qquad \frac{\Delta\nu}{\Delta\nu_\odot} = \left(\frac{M}{M_\odot}\right)^{1/2} \left(\frac{R}{R_\odot}\right)^{-3/2}$$ with $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$0. Empirically, $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$1 in standard asteroseismic practice, and in PMS grids, a power law $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$2 provides an accurate fit (Müllner et al., 2020, Jørgensen et al., 4 Jan 2026).

2. Internal Structure and Accretion Imprints

PMS evolution diverges from the classical Hayashi–Henyey description when accretion is accounted for. Accreting models exhibit distinct signatures in the Brunt–Väisälä frequency profile: $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$3 where $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$4 is local gravity, $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$5 density, $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$6 pressure, $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$7 the adiabatic temperature gradient, $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$8 the actual gradient, and $$\frac{\nu_\mathrm{max}}{\nu_{\mathrm{max},\odot}} = \frac{M/M_\odot}{(R/R_\odot)^2 \sqrt{T_\mathrm{eff}/T_{\mathrm{eff},\odot}}}$$9 the chemical composition gradient. Heat and deuterium delivered by episodic accretion events induce a pronounced inner “bump” in $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$0 at small radii, contrasting with the smooth rise seen in classical PMS structure. This affects the acoustic Lamb frequency ($$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$1), producing a kink in mode propagation cavities. Across 35 simulated accretion histories, the amplitude and location of these interior features vary, introducing measurable changes in mode frequencies up to $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$220\,$$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$3Hz at the ZAMS (Jørgensen et al., 4 Jan 2026).

3. Observational Methodologies: Bayesian Detection and APOLLO

Detection of PMS solar-like oscillations is achieved by exploiting high-quality light curves, advanced filtering, and robust Bayesian model comparison pipelines. The APOLLO software package exemplifies this workflow:

• Light curve pre-processing: KASOC filtering or Everest detrending, $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$4 outlier clipping, gap interpolation, and Lomb–Scargle periodogram generation (Parseval normalization).
• Background modeling: The power spectral density (PSD) is fit with two models: a noise-only background (white noise + multi-component granulation) or a background plus a Gaussian-shaped p-mode excess. The latter includes a sum of three Harvey-like granulation terms and a Gaussian mode envelope centered at $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$5.
• Statistical detection: Bayesian nested sampling (Diamonds) computes the odds ratio (Bayes factor, $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$6) between models: ln $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$7 constitutes “strong evidence” for oscillations per standard statistical criteria.
• Parameter inference: Posterior distributions of all model terms are derived, and if oscillation power is detected, $$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$8 is estimated from autocorrelation analysis in the p-mode envelope region. Priors are informed by solar-calibrated scaling relations and granulation parameter trends (Müllner et al., 2020).

A crucial pre-selection employs the FliPer metric ($$\nu_{\mathrm{max},\odot} \approx 3090\,\mu\mathrm{Hz}$$9) to identify candidate stars with gravitational parameters consistent with expectations from scaling laws.

4. Empirical Confirmation: Key Candidate Detections

Application of these methods to Kepler K2 Campaign 2 Upper Scorpius members (T Tauri stars, $M_\odot$010 Myr) yielded several candidates. EPIC 205375290 displayed a clear p-mode envelope at $M_\odot$1Hz, with ln $M_\odot$2, indicating strong Bayesian evidence for genuine oscillations. Fundamental parameters include $M_\odot$3 K, $M_\odot$4, $M_\odot$5 km s$M_\odot$6, and solar metallicity. Spectroscopic features confirm its youth (Li I absorption, H$M_\odot$7 and Ca II emission lines, weak accretion). Evolutionary models place it below the Kepler LC Nyquist line and confirm its status as a pre-main sequence object. $M_\odot$8 could not be robustly measured due to spectral truncation (Müllner et al., 2020).

5. Acoustic Diagnostics and Frequency Separations

Asteroseismic diagnostics leverage large and small separations: $M_\odot$9 with averages computed for modes within $\nu_\mathrm{max}$0 of $\nu_\mathrm{max}$1. In PMS stars, $\nu_\mathrm{max}$2 increases as the star contracts toward the ZAMS; $\nu_\mathrm{max}$3 concurrently rises but remains lower by $\nu_\mathrm{max}$4–$\nu_\mathrm{max}$5 of $\nu_\mathrm{max}$6, as confirmed in grids spanning 0.7–1.6 $\nu_\mathrm{max}$7. Frequency differences across accretion histories can reach $\nu_\mathrm{max}$8Hz at high radial orders at the ZAMS, while mean and interquartile ranges are at the $\nu_\mathrm{max}$9Hz level (Jørgensen et al., 4 Jan 2026).

Empirical and theoretical tracks in the C–D diagram (plot of $500\,\mu\mathrm{Hz}$0 vs $500\,\mu\mathrm{Hz}$1) are well separated by mass and sensitive, at few-percent level, to individual accretion histories.

6. Detection Prospects and Implications for Stellar Evolution

Observational strategies are shaped by the amplitude, frequency regime, and noise characteristics of PMS oscillators. Predicted mode amplitudes are at the ppm level but may be damped by surface activity and accretion-driven variability. For typical $500\,\mu\mathrm{Hz}$2Hz, oscillations lie beyond the Nyquist frequency of the TESS 30-min cadence but are fully resolvable in PLATO’s fast cadence (6.25–25 s), with anticipated two-year baselines allowing frequency precision better than $500\,\mu\mathrm{Hz}$3Hz.

The high $500\,\mu\mathrm{Hz}$4 and $500\,\mu\mathrm{Hz}$5 of young PMS stars distinguish them from main-sequence stars of equivalent mass. Accretion imprints—buoyancy frequency bumps, small separations, and order-dependent shifts—persist to the ZAMS and are observable as systematic frequency offsets. Upcoming missions (PLATO, K2, TESS) are thus positioned to deliver statistically significant samples of PMS solar-like oscillators.

Observational detection of stochastic, convectively driven p-modes in PMS stars enables direct asteroseismic probes of interior structure, core properties, and convective efficiency during the Hayashi track. This confirms theoretical predictions of mode excitation in magnetically active, accreting young stars (Müllner et al., 2020, Jørgensen et al., 4 Jan 2026). A plausible implication is that measurement of frequency separations and absolute modes in PMS stars could retrospectively constrain individual accretion histories, thereby anchoring evolutionary models in previously inaccessible early stellar phases.