Brought to you by:

Articles

PRECISION ASTEROSEISMOLOGY OF THE PULSATING WHITE DWARF GD 1212 USING A TWO-WHEEL-CONTROLLED KEPLER SPACECRAFT

, , , , , , , , , , , , , , , , and

Published 2014 June 17 © 2014. The American Astronomical Society. All rights reserved.
, , Citation J. J. Hermes et al 2014 ApJ 789 85 DOI 10.1088/0004-637X/789/1/85

0004-637X/789/1/85

ABSTRACT

We present a preliminary analysis of the cool pulsating white dwarf (WD) GD 1212, enabled by more than 11.5 days of space-based photometry obtained during an engineering test of the two-reaction-wheel-controlled Kepler spacecraft. We detect at least 19 independent pulsation modes, ranging from 828.2–1220.8 s, and at least 17 nonlinear combination frequencies of those independent pulsations. Our longest uninterrupted light curve, 9.0 days in length, evidences coherent difference frequencies at periods inaccessible from the ground, up to 14.5 hr, the longest-period signals ever detected in a pulsating WD. These results mark some of the first science to come from a two-wheel-controlled Kepler spacecraft, proving the capability for unprecedented discoveries afforded by extending Kepler observations to the ecliptic.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The endpoints of stellar evolution, white dwarf (WD) stars provide important boundary conditions on the fate of all stars with masses ⩽8 M, as is the case for more than 97% of all stars in our Galaxy, including the Sun. When a WD cools to the appropriate effective temperature to foster a hydrogen partial-ionization zone (∼12, 000 K, although there is a slight dependence on WD mass), global oscillations driven as non-radial g-modes become unstable and reach observable amplitudes (see Winget & Kepler 2008; Fontaine & Brassard 2008).

These hydrogen-atmosphere pulsating WDs (so-called DAVs or ZZ Ceti stars) have spent hundreds of Myr passively cooling before reaching this evolutionary state. Global oscillations provide a unique window below the thin photosphere and deep into the interior of these stars, enabled by matching the observed periods to theoretical periods generated by adiabatic pulsation calculations.

Given the number of free parameters for full asteroseismic fits, the most reliable results require securing a large number of pulsation periods and uniquely identifying the oscillation modes. However, with typical g-mode periods ranging from 100 to 1400 s, ground-based photometry suffers from frequent gaps in coverage, frustrating efforts to disentangle multiperiodic signals and alias patterns.

Multi-site campaigns coordinated across the globe via the Whole Earth Telescope (WET; Nather et al. 1990) have proved the richness of well-resolved WD pulsation spectra. For example, less than a week of nearly continuous observations of the helium-atmosphere (DBV) GD 358 revealed more than 180 significant periodicities in the power spectrum, providing exquisite constraint on the helium-envelope mass, (2.0 ± 1.0) × 10−6 M, the overall mass, 0.61 ± 0.03 M, and the magnetic field strength, 1300  ±  300 G (Winget et al. 1994). Similarly, multiple WET campaigns on the pre-WD PG 1159−035 have revealed 198 periodicities, accurately constraining the mass, rotation rate and magnetic field of this DOV (Wing et al. 1991; Costa et al. 2008).

DAVs have also been extensively studied by more than 10 WET campaigns, with varying results. In part, this is a result of how pulsation modes excited in DAVs are characteristically influenced by the WD effective temperature: hotter DAVs tend to have fewer modes, lower amplitudes and shorter-period pulsations, while cooler DAVs driven by substantially deeper convection zonest end to have more modes at higher amplitude and longer periods (e.g., Mukadam et al. 2006). WET campaigns have borne this out. More than 5 days of nearly continuous monitoring of the hot DAV G226−29 revealed just one significant triplet pulsation (Kepler et al. 1995), whereas the cooler DAV G29−38 has more than a dozen modes of higher amplitude (Kleinman et al. 1994).

In fact, G29−38 illustrates the challenges faced to performing asteroseismology of cooler DAVs: although the WD exhibits at least 19 independent oscillation frequencies, there is significant amplitude and phase modulation of these modes, which change dramatically from year to year (Wing et al. 1990; Kleinman et al. 1998). Another excellent example of this complex behavior is the cool DAV HL Tau 76, which shows 34 independent periodicities along with many oscillation frequencies at linear combinations of the mode frequencies (Dolez et al. 2006). The complex mode amplitude and frequency variations are likely the result of longer-period pulsations having much shorter linear growth times, increasing the prevalence of amplitude and phase changes in cooler DAVs with longer periods (e.g., Goldreich & Wu 1999).

The Kepler mission has already uniquely contributed to long-term distinctions between the handful of hot and cool DAVs eventually found in the original pointing. The longest-studied by Kepler, the cool DAV KIC 4552982 (10, 860 K, log g =8.16) discovered from ground-based photometry (Hermes et al. 2011), shows considerable frequency modulation in the long-period modes present between 829–1292 s (K. J. Bell et al. 2014, in preparation). A much hotter DAV was also observed for six months, KIC 11911480 (12, 160 K, log g =7.94), which shows at least six independent pulsation modes from 172.9–324.5 s that are incredibly stable and evidence consistent splitting from a 3.5 ± 0.5 day rotation rate (Greiss et al. 2014).

After the failure of a second reaction wheel in 2013 May, the Kepler spacecraft has demonstrated a mission concept using two-wheel control, observing fields in the direction of the ecliptic. This mission concept aims to obtain uninterrupted observations for approximately 75 days. As part of an initial test to monitor the two-wheel-controlled pointing behavior on long timescales, short-cadence photometry was collected every minute on the cool DAV GD 1212 during a preliminary engineering run in 2014 January and February.

GD 1212 (V = 13.3 mag) was discovered to pulsate by Gianninas et al. (2006), with roughly 0.5% relative amplitude photometric variability dominant at 1160.7 s. The most recent model atmosphere fits to spectroscopy of GD 1212 find this WD has a Teff =10, 970  ±  170 K and log g =8.03  ±  0.05 (Gianninas et al. 2011), which corresponds to a mass of 0.62  ±  0.03 M after incorporating the latest three-dimensional-model corrections of Tremblay et al. (2013). This puts the WD at a distance of roughly 18 pc, although GD 1212 has the lowest proper motion of any WD within 25 pc of the Sun, 33.6 ± 1.0 mas yr−1 (Subasavage et al. 2009).

In this paper we provide a preliminary analysis of the unique two-wheel-controlled Kepler observations of GD 1212. In Section 2 we outline the observations and reductions. We analyze the independent pulsation modes and nonlinear combination frequencies in Sections 3 and 4, respectively. We conclude in Section 5 with a preliminary asteroseismic interpretation and a discussion of these results in the context of a two-wheel-controlled Kepler mission observing into the ecliptic.

2. OBSERVATIONS AND REDUCTIONS

We observed GD 1212 (GJ 4355, WD 2336−079, WD J233850.74−074119.9) for a total of 264.5 hr using the Kepler spacecraft in two-wheel mode. With only two working reaction wheels, the spacecraft pointing cannot be stabilized in three axes. By pointing the bore-sight close to the plane of the spacecraft's orbit (within about a degree of the ecliptic, Howell et al. 2014), the unconstrained spacecraft roll is placed in equilibrium with respect to the solar pressure, the dominant external force exerted on the spacecraft.

The unconstrained spacecraft roll is in unstable equilibrium, and the pointing needs to be corrected every 3–6 hr by firing the thrusters. No data is lost during a thruster firing (known as a reset); the change in attitude causes a minor change in the pointing, which is accounted for by moving the aperture centroid. The spacecraft rolls by no more than 25 arcsec between resets, which corresponds to motion of a star at the extreme edge of the field of view by no more than 1 pixel. The spacecraft can point at a single field for up to 80 days before the angle of the Sun with respect to the solar panels exceeds allowed limits. Putnam & Wiemer (2014) provides further details of the capabilities and limitations of two-wheel operation.

We collected data on GD 1212 from 2014 January 17 to 2014 February 13 in short-cadence mode, where each exposure is 58.8 s. After the first 2.6 days, the observations were interrupted for 15.1 days by a safe-mode event and subsequent engineering fault analysis. The final 9.0 day light curve can be found in Figure 1. We have removed all points falling more than 4σ from the light-curve mean (112 points), resulting in 2.6 and 9.0 day observations with a duty cycle of 98.2% and 98.9%, respectively.

Figure 1.

Figure 1. Portion of the minute-cadence photometry collected on the pulsating white dwarf GD 1212 in 2014 February. This 9.0 day light curve (shown without any smoothing) establishes the capabilities of a two-wheel-controlled Kepler spacecraft on this Kp = 13.3 mag white dwarf. An additional 2.6 days of data, not shown, were collected on GD 1212 in 2014 January and included in our analysis.

Standard image High-resolution image

In contrast to the primary mission, where only small masks were placed around each star, the data on GD 1212 were collected using a 50 × 50 pixel "super aperture." We expect the aperture size used in two-wheel mode will decrease as confidence in the spacecraft pointing ability increases. The observed pixels were processed through the CAL module of the Kepler pipeline (Quintana et al. 2010) to produce target pixel files (TPFs); light-curve files were not produced for this engineering data. TPF data for all stars observed during this engineering run are available at the MAST Kepler archive.13 GD 1212 is object no. 60017836 in this engineering data.

We estimated a local background for each image using a section of pixels absent of stars and subtracted this. We then created a 50 × 50 pixel image summed in time, and fit a two-dimensional Gaussian to this summed image, which defined a standard deviation in the semi-major axis, a standard deviation in the semi-minor axis, and a rotation angle. We held these constant to fix an elliptical aperture but allowed it to move to account for motion of the star. We extracted all pixels within five times the standard deviation, weighted by how much of the area of each pixel covered by the ellipse.

Our final 9.0 day data set is nearly continuous and has a formal frequency resolution of 1.29 μHz. The median noise level in the FT near 500 and 1500 μHz for this 9.0 day run is roughly 0.0036% (36 ppm). For the entire 26.7 day (42.8% duty cycle) data set on the Kp = 13.3 mag GD 1212, the median noise level is roughly 0.0028% (28 ppm).

3. INDEPENDENT PULSATION MODES

As expected given the relatively cool spectroscopically determined temperature, the pulsation periods excited in GD 1212 are relatively long, ranging from 828.2–1220.8 s. A cool effective temperature is also borne out from model atmosphere fits to the photometry of this WD, which find Teff =10, 940 ± 320 K and log g =8.25 ± 0.03 (Giammichele et al. 2012).

For our discussion, we will adopt the more precise parameters derived from spectroscopy: Teff =10, 970 ± 170 K and log g =8.03 ± 0.05. The spectroscopy of GD 1212 is notable in that much of it was collected in order to identify the then-dominant mode at 1160.7 s, by measuring the pulsation amplitude as a function of wavelength with time-resolved spectroscopy from the Bok 2.3 m telescope at Kitt Peak National Observatory. However, the identification was inconclusive given the variable, low amplitude of the pulsations (Desgranges 2008).

A Fourier transform (FT) of our entire data set is shown in Figure 2, which orients us to the dominant frequencies of variability. An FT of our entire data set finds significant variability ranging from 19.2 μHz (14.5 hr) down to 2703.9 μHz (369.8 s), which is shown in full in the top panel of Figure 2. The only periodicity detected with marginal significance at longer frequencies occurs at 4531.8 μHz (0.12 ppt), which is an instrumental artifact sampling the long-cadence exposures of 29.4 minutes that often appears in short-cadence Kepler data (Gilliland et al. 2010). We do not include this signal in our analysis.

Figure 2.

Figure 2. Fourier transforms (FTs) of our two-wheel-controlled Kepler observations of GD 1212. The top panel shows the FT of the entire data set out to 2800 μHz; there are no significant periodicities at longer frequencies. The middle and bottom panels show the independent pulsation modes in more detail. The blue FT represents only the final 9.0 days of data, while the black FT was calculated from the entire 26.7 day data set. Additionally, we have split the data into four 2.6 day subsets, and we average the amplitudes of the FTs of each of these subsets, shown in gray for the bottom two panels. We mark the pulsation modes adopted, detailed in Table 1, with red points and associated uncertainties at the top of each panel. Additionally, the magenta lines in the top panel mark nonlinear combination frequencies also detected in our data, shown in more detail in Figure 4. Amplitudes are expressed in ppt, where 1 ppt = 0.1% relative amplitude.

Standard image High-resolution image

The highest-amplitude variability in GD 1212 occurs in the region between 810–1210 μHz (826–1234 s). This region evidences at least 19 independent pulsation modes, which we show in more detail in the bottom panels of Figure 2. However, the amplitudes of variability in this region are not stable over our 26.7 days of observations, so we additionally display the FT of just our final 9.0 day light curve. We also find a number of nonlinear combination frequencies of these highest-amplitude pulsations, which we discuss further in Section 4.

The amplitude and frequency variability of the pulsations in the region between 810–1210 μHz are best shown by a running FT of our 9.0 day light curve using a 4.0 day sliding window, shown in Figure 3. For example, the highest-amplitude peak in the FT near 840.0 μHz is essentially monotonically increasing in amplitude over these last 9.0 days, and is broadly unstable in frequency. Conversely, the running FT shows that the second-highest peak in the FT near 910.5 μHz is decreasing in amplitude, and appears further to bifurcate.

Figure 3.

Figure 3. Running FT of the range of independent pulsation modes detected in our 9.0 day, high-duty-cycle two-wheel-controlled Kepler light curve of GD 1212. This figure uses a four-day sliding window, and darker grayscale corresponds to higher amplitudes. This demonstrates the lack of amplitude and frequency stability for the pulsations in GD 1212, even over less than a week.

Standard image High-resolution image

A common cause of amplitude modulation in pulsating stars is the presence of additional nearby frequencies that are unresolved over the course of observations. The pulsational variability in DAVs is the result of non-radial stellar oscillations, so rotation acts to break the spherical symmetry of a pulsation mode and generates multiplets about the m = 0 component, spaced by an amount proportional to the rotation rate (e.g., Dolez & Vauclair 1981). WDs with identified pulsations all vary in ℓ = 1, 2 g-modes, and thus DAVs which have measured rotational splittings range from 2.5–19.0 μHz (see Kawaler 2004; Fontaine & Brassard 2008 and references therein).

The DAV G226−29 is an excellent example of a closely spaced multiplet affecting pulsation amplitudes in a WD (Kepler et al. 1983). In this WD, the 109 s pulsation is actually a superposition of three signals at 109.08648 s, 109.27929 s, and 109.47242 s, which are separated by roughly 16.1 μHz and require more than 17.2 hr of observations to resolve the beat period of the closely spaced signals. Shorter observations show the signal sinusoidally varying at the beat period of 17.2 hr.

However, rotational splittings are unlikely to explain the incoherent changes in the running FT of GD 1212, nor the disappearance of f5 near 847.2 μHz, which was strong during the first 2.6 days and was virtually unseen in our final 9.0 days of data. Instead, we are likely observing genuine amplitude and frequency variability internal to the star, as has been witnessed during long campaigns on other cool DAVs (e.g., G29−38; Kleinman et al. 1998) and DBVs (e.g., GD 358; Provencal et al. 2009).

The frequency, amplitude, and possibly phase variability observed in GD 1212 cause difficulty in defining the true pulsation periods. For hotter DAVs, which show exceptionally stable oscillations, pre-whitening the light curve by the highest peaks in an FT removes most power in that region and allows for a relatively simple extraction of the periods present in the star (e.g., Greiss et al. 2014). However, none of the peaks in either the 26.7 day FT or the 9.0 day FT of GD 1212 can be smoothly pre-whitened in the standard way, because fitting a pure sine wave does not completely remove the pulse shape.

One can always decompose the broad peaks in an FT with a linear combination of sine waves, but not all of these sine waves represent physical modes. For example, using just the final 9.0 days of data, the highest-amplitude signal near 840.0 μHz requires five sine waves to reproduce the broad peak: 840.21 μHz (2.48 ppt), 839.18 μHz (1.54 ppt), 842.02 μHz (1.10 ppt), 838.00 μHz (0.82 ppt), and 836.83 μHz (0.38 ppt). This series of pre-whitened peaks in all likelihood corresponds to only one independent pulsation mode.

The pre-whitening method does provide a quantitative test for significance. We can iteratively fit and pre-whiten the highest peaks in the FT until there are none above some threshold; we estimate this locally for each frequency by calculating the median value, σ, of the FT in a ±200 μHz sliding window, and mark as significant the signals that exceed 5σ. Since the 15.1 day gap in our whole data set strongly degrades the window function, we have adopted this approach for our final 9.0 days of data and our initial 2.6 day light curve. This provides an input list of 48 significant signals.

We do not, however, use the pre-whitening method to define adopted periods and the associated period uncertainties of the independent pulsation modes, which form the input for asteroseismic modeling. Performing a linear least-squares fit of the periods determined from the iterative pre-whitening technique would greatly underestimate the period uncertainties.

Instead, we have fit a Lorentzian profile to a ±5 μHz range of the regions of power defined significant using our pre-whitening method, and use the central peak, half-width-half-maximum, and intensity of this Lorentzian to define the adopted periods, period uncertainties, and amplitudes of the broad peaks in the FT. These values provide the red points and uncertainties shown in Figure 2 that conservatively represent the periods of excited variability, and are detailed in Table 1. Lorentzian profiles are frequently used to describe stochastic oscillations driven by turbulent convection in cool main-sequence and red giant stars, but our use of these profiles is primarily out of convenience and does not suggest that stochastic driving is undergoing in GD 1212 (a discussion of driving in DAVs can be found in Brickhill 1991).

Table 1. Frequency Solution for GD 1212

ID Period Frequency Amplitude
(s) (μHz) (ppt)
Independent Pulsation Modes
f14 1220.75 ± 7.15 819.17 ± 4.80 0.23
f1 1190.53 ± 2.28 839.96 ± 1.61 2.11
f4 1180.40 ± 4.02 847.17 ± 2.89 0.85
f18 1166.67 ± 4.81 857.14 ± 3.54 0.21
f7 1147.12 ± 3.19 871.75 ± 2.42 0.52
f17 1125.37 ± 3.01 888.60 ± 2.38 0.22
f2 1098.36 ± 1.65 910.45 ± 1.37 1.04
f12 1086.12 ± 3.27 920.71 ± 2.78 0.24
f5 1063.08 ± 4.13 940.66 ± 3.66 0.63
f20 1048.19 ± 4.01 954.02 ± 3.65 0.13
f9 1025.31 ± 2.26 975.31 ± 2.15 0.49
f11 1008.07 ± 1.20 991.99 ± 1.18 0.31
f3 987.00 ± 3.73 1013.17 ± 3.82 0.95
f16 956.87 ± 4.91 1045.08 ± 5.36 0.22
f15 871.06 ± 2.13 1148.02 ± 2.81 0.22
f8 857.51 ± 0.86 1166.17 ± 1.17 0.49
f13 849.13 ± 0.76 1177.68 ± 1.05 0.23
f10 842.96 ± 1.02 1186.30 ± 1.43 0.46
f6 828.19 ± 1.79 1207.45 ± 2.61 0.53
f21 371.05 ± 0.36 2695.04 ± 2.63 0.07
f19 369.83 ± 0.17 2703.91 ± 1.24 0.15
Nonlinear Combination Frequencies
f10f8 50,500 ± 10,600 19.82 ± 4.16 0.24
f7f1 31,540 ± 3580 31.70 ± 3.60 0.27
f2f7 25,990 ± 2870 38.47 ± 4.24 0.35
f18f1 20,920 ± 2280 47.80 ± 5.21 0.30
f2f1 14,210 ± 539 70.37 ± 2.67 0.21
f5f1 10,463 ± 428 95.57 ± 3.91 0.21
f3f2 9967 ± 274 100.33 ± 2.75 0.25
f3f1 5910 ± 162 169.19 ± 4.66 0.22
2f1 595.45 ± 1.29 1679.41 ± 3.65 0.08
f1 + f2 571.18 ± 1.30 1750.75 ± 4.00 0.08
f5 + f7 552.36 ± 0.68 1810.41 ± 2.23 0.08
f1 + f3 540.23 ± 0.64 1851.07 ± 2.18 0.09
f2 + f3 519.94 ± 0.88 1923.29 ± 3.27 0.09
f1 + f8 498.36 ± 0.88 2006.60 ± 3.54 0.11
f4 + f8 497.13 ± 0.99 2011.53 ± 4.01 0.08
f5 + f8 474.56 ± 0.45 2107.22 ± 1.99 0.12
f3 + f8 458.89 ± 0.49 2179.19 ± 2.35 0.12

Download table as:  ASCIITypeset image

This method strongly underestimates the adopted pulsation amplitudes, but theoretical pulsation calculations do not account for mode amplitudes, and our results reasonably represent the two quantities requisite for a seismic analysis: the pulsation periods and associated uncertainties. We italicize the signals at 371.1 s and 369.8 s in Table 1 because they may not in fact be independent pulsation modes (see discussion at the end of Section 4).

When establishing the coherence of variability over long timescales, it is often instructive to calculate FTs from multiple different subgroups of the data of identical resolution and then average the amplitudes of the FTs of each of these subsets. This effectively throws away the phase information of the different data sets, and identifies the most coherent regions of the FT. Figure 2 shows such an average FT for four equal-length (2.6 day) subsets of our data on GD 1212. We find that smaller 2.6 day subsets are not sufficient to resolve some closely spaced variability that can be resolved within longer subsets, such as the FT of our 9.0 day light curve.

It is possible that this frequency broadening for the averaged FTs is the result of unresolved rotational multiplets embedded within the natural frequency variability of the underlying pulsation modes. The average half-width-half-maximum of the profiles fit to the highest seven peaks of the averaged FT is 3.14 μHz, compared to 1.12 μHz for the 9.0 day FT. This difference could be explained by rotational splittings arising from a WD rotation rate between 1.9–5.2 days if these are ℓ = 1 modes and between 6.1–17.2 days if these are ℓ = 2 modes.

An inferred rotation rate of 1.9–17.2 days provides one of the first constraints on the rotation rate of a cool DAV, albeit a very coarse estimate. Hotter, stable DAVs show detected rotation rates from rotational splittings between 0.4–2.3 days (Kawaler 2004; Fontaine & Brassard 2008).

4. NONLINEAR COMBINATION FREQUENCIES

In addition to the 19 significant independent pulsation modes we identify in Section 3, we see a number of nonlinear combination frequencies in our data that arise at summed and difference frequencies of the independent modes. For example, the peak at 2107.3 μHz is, within the uncertainties, a linear sum of the two independent pulsations f4 = 940.3 μHz and f7 = 1166.2 μHz.

Combination frequencies are a common feature in the pulsation spectra of variable WDs, which are created by a nonlinear distortion of an underlying linear oscillation signal. Work by Brickhill (1992) attributed these nonlinear distortions to the changing thickness of the convection zone of a DAV, a result of local surface temperature variations from the underlying global stellar oscillations. The convection zone depth is extremely sensitive to temperature (∝T−90; Montgomery 2005), distorting the observed emergent flux. For the hottest DAVs, the nonlinear response of the flux to a temperature perturbation (i.e., the "T 4" nonlinearity) may also play a role in creating these signals (Brassard et al. 1995).

Observations of the rate of period change of combination frequencies show that they match identically the rates of their parent modes, establishing that these signals are not independently excited pulsations but rather nonlinear distortions (Hermes et al. 2013). Thus, these signals provide no additional asteroseismic constraints on the interior structure of the star. However, the amplitude ratios of these combination frequencies to their parent modes can yield insight into the depth of the convection zone, which is responsible for driving the underlying independent pulsations (e.g., Montgomery et al. 2010).

We detail the observed nonlinear combination frequencies in GD 1212 and identify their likely parent modes in Table 1. Their amplitudes and frequencies were determined in an identical manner to the independent pulsation modes, as discussed in Section 3. However, we have used a less stringent test for significance, and include as significant the signals that exceed 4σ. We mark these adopted nonlinear combination frequencies as magenta points in our FT of all the data, shown in the top panel of Figure 2. We show a more detailed view of the regions of these combination frequencies in Figure 4.

Figure 4.

Figure 4. Fourier transforms (FTs) of the nonlinear combination frequencies detected in our observations of GD 1212. As in Figure 2, the blue FT represents only the final 9.0 days of data, while the black FT was calculated from the entire 11.6 days of data spread over 26.7 days. Additionally, the average FT of four 2.6 day subsets is shown in gray. We mark the adopted combination frequencies, detailed in Table 1, with magenta points and associated uncertainties at the top of each panel.

Standard image High-resolution image

Good evidence that these are in fact combination frequencies comes qualitatively from the region in the FT near 2006 μHz (bottom panel of Figure 4) as compared to the region near 847 μHz (middle panel of Figure 2). There appear two combination frequencies, 2006.6 μHz and 2011.5 μHz, which are combination of f1 + f7 and f5 + f7, respectively. The power for the f5 + f7 combination frequency is nearly nonexistent in our final 9.0 days of data, exactly as the power for f5 near 847.2 μHz is diminished in our final 9.0 days of data.

Unfortunately, the underlying amplitude and frequency variability of the independent parent modes extends to these combination frequencies, broadening their power in the FTs. Therefore, these nonlinearities do not provide much assistance in refining the periods of the parent modes. Still, the eight significant low-frequency difference frequencies, ranging from 19.2–170.0 μHz (1.6–14.5 hr), are the longest-period signals ever detected in a pulsating WD, accessible only because of the exceptionally long and uninterrupted space-based observations provided by the Kepler spacecraft.

We preliminarily identify the variability at f21 = 371.1 s and f19 = 369.8 s as independent pulsation modes, but they arise at substantially shorter periods, far isolated from the other independent pulsations, which are all longer than 828.2 s. The signals f21 and f19 may instead be more complicated nonlinear combination frequencies. For example, f21 may be the sum 2f1 + f2. However, no obvious frequency combination can produce the higher-amplitude f19.

It is also possible that the signal near 370 s may be a new two-wheel-controlled Kepler instrumental artifact. However, we do not observe a significant peak near 2700 μHz in any other short-cadence targets observed during the two-wheel-controlled Kepler engineering test run.

5. DISCUSSION AND CONCLUSIONS

We report a detailed census of pulsation modes excited in the cool DAV GD 1212 using more than 11.5 days of data collected during a two-wheel concept engineering test of the repurposed Kepler spacecraft. Our results suggest the presence of at least 19 independent pulsation modes excited to observable amplitudes, with up to 2.5% peak-to-peak variability.

Observed frequency and amplitude modulation complicates a precise determination of the pulsation periods in GD 1212, likely the result of these relatively long-period oscillations (828.2–1220.8 s) having short linear growth times. We have defined a conservative estimate of the excited pulsations by fitting Lorentzian functions to the broad bands of power in an FT of our entire data set, yielding an average uncertainty of 3.0 s (<0.5%) for each adopted period listed in Table 1.

A complete asteroseismic solution of the data in hand would be premature until we have more secure identification of the spherical degree (ℓ values) of the independent pulsation modes, the subject of a future paper. We can make a first attempt at these identifications given that such long-period modes likely have high radial order (k > 17). This is near the asymptotic limit, in which patterns in consecutive period spacing can yield insight into the spherical degree of the modes.

There appear to be two regions of semi-evenly spaced independent pulsations, one at longer periods (involving f1-f7-f2-f5-f9-f3) that have a mean period spacing of 41.5 ± 2.5 s. Additionally, four shorter-period modes (f15-f8-f10-f6) have a spacing of 14.4  ±  1.5 s. The asymptotic period spacing for ℓ = 1 modes of a cool (11, 147 K), 0.837 MWD with a 10−5.408 M hydrogen layer mass is roughly 40.3 s, according to the models of Romero et al. (2012). However, this model does not accurately predict the other period spacings, as they calculate ℓ = 2 modes to have a 23.3 s mean period spacing.

The asymptotic period spacing for higher spherical degrees decreases as $\sqrt{\ell (\ell +1)}$, so it is possible that the shorter-period sequence is composed of a series of ℓ = 3 modes. However, ℓ = 3 modes have never been previously inferred in WDs, although ℓ = 4 modes have been proposed in G 185-32 (Thompson et al. 2004). Limb-darkening and geometric effects greatly reduce the observed amplitudes of ℓ = 1, 2 modes compared to ℓ = 3 modes at shorter (near-UV) wavelengths (e.g., Kepler et al. 2000), so multicolor observations of GD 1212 could conclusively test this hypothesis.

It is perhaps more likely that the two seemingly evenly spaced regions are not consecutive radial order modes of the same spherical degree at all, but rather a coincidence of overlapping ℓ = 1 and ℓ = 2 modes. A coarse grid search using the models described in Castanheira & Kepler (2008), allowing the highest-amplitude periodicities to be ℓ = 1, 2, can modestly reproduce the observed periods if f8, f3, f5, f7, f1 are ℓ = 1 modes and if f6, f10, f15, f2 are ℓ = 2 modes, if the WD is 10, 700 K, 0.81 M, and has a 10−5.5 M hydrogen layer.

Another interpretation of the period-spacing patterns suggests that all modes with periods in excess of 956 s are ℓ = 1 modes, with overtones spaced by about 37 s. Each overtone would appear as an m = ±1 doublet split in total by roughly 16.5 μHz. Modes with shorter periods roughly follow a series spaced by about 21 s (consistent with l = 2 given the l = 1 spacing for the longer-period modes), with each overtone expressed as a doublet with roughly 19.6 μHz splittings. It is not uncommon to observe only the ±1 components of an l = 1 triplet and not the m = 0 central component, due to inclination or excitation effects (e.g., Wing et al. 1991).

More quantitative answers will require a full asteroseismic analysis, matching all 19 independent pulsation modes to the periods expected from theoretical models, which should explicitly reveal the overall mass, temperature, and hydrogen-layer mass of this cool pulsating WD.

As DAVs cool and approach the red edge of the ZZ Ceti instability strip, the observed pulsation periods generally increase, often exhibiting large-amplitude, nonlinear pulse shapes (Kanaan et al. 2002). The periods observed in the cool (10, 970 K) GD 1212 fit this trend: the weighted mean period of the observed pulsations is roughly 1061 s, one of the longest observed for any DAV (Mukadam et al. 2006). However, the pulsation amplitudes and nonlinearities observed in GD 1212 are notably weak, suggesting this WD is a somewhat rare low-amplitude DAV near the red edge.

Our final 9.0 day light curve displayed in Figure 1 is, to date, the highest-quality continuous look at a pulsating WD, unmolested by differential atmospheric effects, clouds, or sunrises. Much of the long-term variability about the mean of the light curve is real, caused by nonlinear combination frequencies of the excited pulsation modes at periods extending up to 14.5 hr. Such long-period signals have never been detected from ground-based campaigns.

This relatively short engineering test proves the impact that a two-wheel-controlled Kepler spacecraft, observing in the ecliptic, can have in making new discoveries in the exploration of stellar interiors and asteroseismology. There are at least 36 known pulsating WDs brighter than g < 18.0 mag within 12° of the first nine proposed K2 fields (see Howell et al. 2014). We expect that such a mission will afford us unprecedented insight into the heretofore elusive cause of why pulsations in DAVs shut down at the red edge, as well as uncovering additional low-amplitude pulsations simply inaccessible from ground-based observations that will further constrain asteroseismic modeling of these degenerate stars.

We acknowledge Working Group 11 of the Kepler Asteroseismic Science Consortium for the eager discussion and analysis that led to this rapid result. We thank the referee, S.O. Kepler, for helpful comments that improved aspects of the discussion. J.J.H. acknowledges funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement No. 320964 (WDTracer). M.H.M. and D.E.W. gratefully acknowledge the support of the NSF under grants AST-0909107 and AST-1312678 and the Norman Hackerman Advanced Research Program under grant 003658-0252-2009. M.H.M. acknowledges the support of NASA under grant NNX12AC96G, and D.E.W. acknowledges the support of NASA under grant NNX13AC23G. S.B. is supported by the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO). D.H. acknowledges support by NASA under grant NNX14AB92G issued through the Kepler Participating Scientist Program.

Facilities: Kepler - The Kepler Mission

Footnotes

Please wait… references are loading.
10.1088/0004-637X/789/1/85