Kepler Discovery of GW Vir Pulsations of the Central Star of Planetary Nebula Kn 61

Technical Summary: Kepler Discovery of GW Vir Pulsations of the Central Star of Planetary Nebula Kn 61

Problem and Context
PG 1159 stars represent a critical, yet poorly understood, transitional evolutionary phase between the post-asymptotic giant branch (post-AGB) and the white dwarf cooling stage. A subset of these stars, known as GW Vir stars, exhibit non-radial g-mode pulsations driven by the $\kappa-\gamma$ mechanism in partial ionization zones of carbon and oxygen. A significant open question in the field is the "nitrogen dichotomy," a hypothesis suggesting that all nitrogen-rich (N-rich) PG 1159 stars are pulsators, while nitrogen-poor (N-poor) stars are not. Furthermore, the central star of the planetary nebula Kn 61 (KIC 3231337) had previously shown sporadic, long-duration brightening events in Kepler Long Cadence data, but its pulsational nature remained unconfirmed due to the limitations of ground-based observations.

Methodology
The authors analyzed Kepler Short Cadence (SC) data from Quarter 12.3 (one month of observations) to search for high-frequency pulsations in Kn 61. The data were processed using the KASOC pipeline to remove transits and trends. Frequency analysis was performed using the Period04 software, employing Fourier amplitude spectra up to the Nyquist frequency. The authors applied a signal-to-noise (S/N) threshold of $\ge 5$ to identify significant peaks, alongside a Rayleigh resolution criterion.

To characterize the pulsations, the team searched for rotational splitting and mean period spacing ($\Delta\Pi$). They utilized the Fourier transform method, the Inverse Variance method, and the Kolmogorov-Smirnov (K-S) statistical test to identify period spacings typical of GW Vir stars. Once a sequence of modes was identified, they applied asymptotic theory to determine the spherical degree ($\ell$) and radial orders ($k$).

The asteroseismic mass was derived using the relation between the constant period spacing ($\Pi_0$) and stellar mass, luminosity, and helium fraction, as established by Kawaler & Bradley (1994). This was compared against an evolutionary mass derived by interpolating the star's position on the Hertzsprung-Russell diagram against post-AGB evolutionary tracks (Miller Bertolami & Althaus 2006).

Additionally, the authors analyzed Kepler Long Cadence (LC) data from Quarters 10–13 to characterize the sporadic brightening events previously noted by De Marco et al. (2015). They flattened the light curve using a Tukey biweight algorithm and identified events where the smoothed flux exceeded three standard deviations. They calculated the equivalent duration (ED) and estimated the energy of these events by modeling the relationship between flux increases in the Kepler bandpass and bolometric luminosity changes, assuming the brightening was caused by an increase in effective temperature.

Key Results

1. Pulsation Discovery: The analysis detected four significant pulsation peaks (S/N $\ge 5$) in the frequency range of 87–111 d$^{-1}$ (periods of 812–991 s), consistent with g-modes in GW Vir stars. Two additional peaks with S/N = 4.1 were also noted.
2. Mode Identification and Period Spacing: The authors identified a sequence of three modes ($f_1, f_2, f_3$) with a mean period spacing of $\Delta\Pi = 21.526(6)$ s. This spacing is consistent with dipole ($\ell=1$) modes. The authors assigned radial orders $k \approx 36, 42,$ and $46$ to these modes.
3. Rotational Splitting: A potential doublet ($f_5, f_6$) with a separation of $0.2256$d$^{-1}$ was identified. If interpreted as components of an $\ell=1$ triplet, this implies a rotation period of approximately 2.22 days or 4.43 days. However, the authors note this identification is tentative due to the sparse spectrum.
4. Mass Determination: Based on the $\ell=1$ period spacing, the asteroseismic mass was calculated to be $M_{\text{astro}} = 0.551(6) M_{\odot}$. This is consistent with the evolutionary mass of $M_{\text{evol}} = 0.575(19) M_{\odot}$ derived from evolutionary tracks, validating the asteroseismic analysis within $1.2\sigma$.
5. Brightening Events: The study characterized 49 sporadic brightening events in the LC data. These events last 1–3 days, increase flux by 2–5%, and recur with a median delay of 5 days. The median equivalent duration is 26 minutes. Assuming the brightening results from an increase in effective temperature, the estimated energy of these events is $\sim 4 \times 10^{40}$ erg.

Significance and Claims
The paper claims that the discovery of multiperiodic pulsations in Kn 61 categorizes it as a GW Vir pulsator. As Kn 61 is spectroscopically confirmed as an N-rich PG 1159 star, this finding supports the nitrogen dichotomy hypothesis that all N-rich PG 1159 stars are pulsators.

The detection of a clear mean period spacing allows for a robust asteroseismic mass determination, making Kn 61 the 12th pulsating PG 1159 star with such a mass measurement. The consistency between the asteroseismic and evolutionary masses reinforces the validity of the underlying stellar evolution models for this phase.

Regarding the brightening events, the authors note their similarity to outbursts in pulsating DA white dwarfs (DAVs) but highlight a massive energy difference ($\sim 10^{40}$ erg vs. $\sim 10^{34}$ erg for DAVs). They explicitly state that the mechanism for these events is unknown and that their energy estimation relies on the unverified assumption that the brightness increase is solely due to a temperature rise. Consequently, they do not claim a definitive physical link to DAV outbursts but note the intriguing fact that the only PG 1159 star observed with such behavior is a pulsator.

The authors conclude that despite the sparse pulsation spectrum, the derived period spacing provides sufficient constraints to determine fundamental stellar parameters, adding Kn 61 as a valuable target for future detailed asteroseismic modeling.