Asteroseismic Analysis of WFST J053009.62+594557.0

• The paper demonstrates robust asteroseismic modeling of the DAV white dwarf using high-precision photometry, spectroscopy, and advanced WDEC simulations.
• It identifies three pulsation modes with precise frequency and amplitude measurements, resulting in tightly constrained global and interior stellar parameters.
• The study highlights the significance of thin hydrogen and helium envelopes in DAV evolution, supported by consistent Gaia astrometry and SED analyses.

WFST J053009.62+594557.0 is a newly identified faint DA white dwarf variable (ZZ Ceti or DAV) discovered by the Wide Field Survey Telescope (WFST). Its asteroseismological properties were characterized through high-precision g-band photometric monitoring, spectroscopy, and detailed theoretical modeling using the White Dwarf Evolution Code (WDEC). This object is notable for its low Gaia G-band magnitude (G ≈ 19.13) and its placement among the thinnest-envelope DAVs studied, with exceptionally well-constrained global and interior parameters. Multi-faceted observational constraints, including Gaia astrometry, spectral fitting, and broadband SED analysis, robustly converge on the asteroseismological solution.

1. Observational Campaign and Photometric Analysis

WFST J053009.62+594557.0 was monitored with the 2.5 m WFST situated on Saishiteng Mountain (Lenghu), utilizing a 9k × 9k CCD and a 6.5 deg² field of view in the SDSS-like g band. The observational cadence was set at 20 s per exposure, with an effective temporal resolution of ~50 s after accounting for overhead. Two observation campaigns (>2 h each) occurred on 2023-09-18 and 2023-11-16.

Photometric reduction was processed via an LSST-derived pipeline, with PSF photometry anchored to Gaia DR3 astrometry and zero-points from Pan-STARRS DR2. Detected periodic signals exhibited high signal-to-noise ratios (SNR ≥ 5.1); three independent pulsational modes exceeded SNR = 4 (the standard detection threshold).

Mode	Frequency (μHz)	δ Frequency (μHz)	Period (s)	δ Period (s)	Amplitude (ppt)	SNR
f₁	2496.005	0.030	400.640	0.005	24.198	6.33
f₂	3247.170	0.008	308.435	0.001	20.737	6.99
f₃	3467.206	0.013	288.417	0.001	15.156	5.12

Each mode’s frequency and amplitude were determined via Period04 analysis with Monte Carlo estimation of uncertainties.

2. Spectroscopic Classification and External Constraints

Follow-up spectroscopy was conducted with the P200/DBSP instrument on the 200-inch Hale Telescope. The dichroic D55 configuration was used: blue arm with a 600 l/mm grism (blaze 3780 Å), red arm with a 316 l/mm grism (blaze 7150 Å), and a 1.5″ slit for a total exposure of 1800 s. Data reduction followed standard PypeIt and IRAF long-slit pipelines.

Balmer-line fits yielded:

• $T_{\rm eff} = 11\,609 \pm 605$ K
• $\log g = 8.05 \pm 0.36$
• Mass $M_* = 0.63 \pm 0.22\,M_\odot$

These atmospheric parameters are typical for ZZ Ceti variables and are consistent with WD evolutionary expectations for the corresponding mass/temperature regime.

Additional constraints derive from Gaia DR3 parallax ($\pi = 3.4973 \pm 0.2282$ mas, corresponding to $d = 285.94 \pm 18.66$ pc), Gaia XP spectroscopy ($T_{\rm eff} = 11\,657 \pm 505$ K, $M = 0.619 \pm 0.070\,M_\odot$), Gaia color–magnitude diagram ($T_{\rm eff} = 12\,110 \pm 770$ K, $M = 0.63 \pm 0.08\,M_\odot$), and SED fitting ($T_{\rm eff} \simeq 11\,000$ K, $M = 0.56\,M_\odot$).

3. Asteroseismological Modeling with WDEC

To probe the inner structure, a grid of $\sim 1.45 \times 10^7$ WD models was generated with the WDEC. Model input physics included:

• Equation of state and opacities from MESA tables
• Convection via Mixing-Length Theory ($\alpha=0.6$)
• Helium diffusion coefficients set to $0.12$ at the envelope and He-layer base

The parameter space encompassed:

• $T_{\rm eff}$: $10\,600$–$12\,600$ K
• $M_*$: $0.50$–$0.85\,M_\odot$
• $-\log(M_{\rm env}/M_*)$: $1.50$–$3.00$
• $-\log(M_{\rm He}/M_*)$: $2.00$–$5.00$
• $-\log(M_{\rm H}/M_*)$: $4.00$–$10.00$
• Helium mass fraction in the mixed region: $0.10$–$0.90$
• Central oxygen mass fraction and transition widths (C/O profile parameters $h_1$–$h_3$, $w_1$–$w_3$)

Model optimization minimized the root-mean-square (RMS) period difference for the three observed periods using

$$\chi^2 = \frac{1}{N}\sum_{i=1}^N (P_{\rm obs,\,i} - P_{\rm calc,\,i})^2$$

with $N=3$.

4. Mode Identification and Model Fit

The best-fitting model (referred to as Model 2) precisely reproduced the observed periods ($<0.02\,$s deviation for each mode). Identified modes and matches:

$P_{\rm obs}$ (s)	$P_{\rm calc}$ (s)	$|{\Delta P}|$ (s)	$\ell$	$k$
288.417	288.436	0.019	1	3
308.435	308.432	0.003	2	7
400.640	400.642	0.002	1	5

Key model parameters ($\pm1\,\sigma$ from $1/\chi^2$ width):

• $T_{\rm eff} = 11\,850 \pm 10$ K
• $M_* = 0.600 \pm 0.005\,M_\odot$
• $\log g = 8.02 \pm 0.01$
• $-\log(M_{\rm env}/M_*) = 1.50 \pm 0.01$
• $$-\log(M_{\rm He}/M_*) = 4.00 \pm 0.01\,\,\Rightarrow\, M_{\rm He} / M_* \approx 10^{-4}$$
• $$-\log(M_{\rm H}/M_*) = 7.00 \pm 0.01\,\,\Rightarrow\, M_{\rm H} / M_* \approx 10^{-7}$$
• Central oxygen mass fraction: $0.61 \pm 0.01$
• Minimum fit RMS: $\chi^2_{\rm min} = 3.58 \times 10^{-4}\,\mathrm{s}$

These parameters are in close agreement with those derived from independent Gaia and SED analyses.

5. Synthesis: Independent Consistency and Physical Interpretation

The convergence between asteroseismic, photometric, spectroscopic, and astrometric solutions is robust. The asteroseismic distance, deduced via model luminosity and bolometric corrections,

$d_{ast} = 281.84 \pm 1.88\,\mathrm{pc}$

is within $1.45$% of the Gaia parallax distance, establishing external consistency.

Structurally, WFST J053009.62+594557.0 is characterized by:

• An unusually thin helium layer ($M_{\rm He}/M_* \approx 10^{-4}$) and an even thinner hydrogen envelope ($M_{\rm H}/M_* \approx 10^{-7}$), placing it among the thinnest-envelope DAVs studied.
• A core oxygen mass fraction ($\approx 0.61$) typical of post-main-sequence progenitors near $2.9\,M_\odot$. A pure-carbon buffer overlays the mixed C/O interior for $-\log[1-M_r/M_*] \approx 1.6$–$3.5$.
• The evolutionary pathway is well-matched by a MESA track from a $2.9\,M_\odot$ ZAMS star, culminating in a $0.60\,M_\odot$ WD at $T_{\rm eff}\simeq11\,850$ K, reinforcing the consistency of the asteroseismic model.

6. Theoretical Framework and Governing Equations

The asteroseismological modeling is underpinned by the theory of non-radial g-mode oscillations in white dwarfs. The adiabatic pulsation equation can be expressed (schematically) as:

$$\frac{d^2\xi_r}{dr^2} + \cdots + N^2\,\xi_r = \omega^2\,\xi_r$$

where $N$ is the Brunt–Väisälä (buoyancy) frequency:

$$N^2 = \frac{g}{H_P}\left(\nabla_{\rm ad} - \nabla + B\right)$$

with $B$ capturing the contribution from composition gradients.

For high-order g-modes, the asymptotic period spacing is:

$$\Delta P_\ell \approx \frac{\Pi_0}{\sqrt{\ell(\ell+1)}}$$

with

$$\Pi_0 = 2\pi^2\left[\int_{r_1}^{r_2}\frac{N}{r}dr\right]^{-1}$$

The parameter optimization minimized:

$$\chi^2 = \frac{1}{N}\sum_{i=1}^N (P_{\rm obs,\,i} - P_{\rm calc,\,i})^2$$

7. Significance within White Dwarf Pulsation Studies

WFST J053009.62+594557.0 represents a rare example of a very faint DAV ($G \approx 19.13$ mag), for which a multi-modal observational campaign and exhaustive theoretical search yield a convergent and tightly constrained asteroseismic solution. The extremely thin H/He layers and well-determined core O mass fraction provide key constraints on DAV envelope evolution and core composition, informing scenarios of post-AGB mass loss and chemical stratification. The exceptional agreement between the asteroseismic and Gaia parallax distances exemplifies the current precision attainable in white dwarf structural inference, supporting both the modeling methodology and the input physics underlying the WDEC grid (Yonghui et al., 2 Jan 2026).