High-Time-Resolution Photometry of AR Scorpii: Confirmation of the White Dwarf’s Spin-Down
The unique binary AR Scorpii consists of an asynchronously rotating, magnetized white dwarf (WD) that interacts with its red-dwarf companion to produce a large-amplitude, highly coherent pulsation every 1.97 minutes. Over the course of two years, we obtained thirty-nine hours of time-resolved, optical photometry of AR Sco at a typical cadence of 5 seconds to study this pulsation. We find that it undergoes significant changes across the binary orbital period and that its amplitude, phase, and waveform all vary as a function of orbital phase. We show that these variations can be explained by constructive and destructive interference between two periodic, double-peaked signals: the spin-orbit beat pulse, and a weaker WD spin pulse. Modelling of the light curve indicates that in the optical, the amplitude of the primary spin pulse is 50% of the primary beat amplitude, while the secondary maxima of the beat and spin pulses have similar amplitudes. Finally, we use our timings of the beat pulses to confirm the presence of the disputed spin-down of the WD. We measure a beat-frequency derivative of Hz s and show that this is attributable to the spin-down of the WD. This value is approximately twice as large as the estimate from Marsh et al. (2016) but is nevertheless consistent with the constraints established in Potter & Buckley (2018). Our precise measurement of the spin-down rate confirms that the decaying rotational energy of the magnetized white dwarf is sufficient to power the excess electromagnetic radiation emitted by the binary.
0000-0002-7181-2554]Robert A. Stiller \move@AU\move@AF\@affiliationDepartment of Physics, University of Notre Dame, Notre Dame, IN 46556, USA
0000-0001-7746-5795]Colin Littlefield \move@AU\move@AF\@affiliationDepartment of Physics, University of Notre Dame, Notre Dame, IN 46556, USA
0000-0003-4069-2817]Peter Garnavich \move@AU\move@AF\@affiliationDepartment of Physics, University of Notre Dame, Notre Dame, IN 46556, USA
0000-0003-4773-4602]Charlotte Wood \move@AU\move@AF\@affiliationDepartment of Physics, University of Notre Dame, Notre Dame, IN 46556, USA
AAVSO / Vereniging Voor Sterrenkunde (VVS), Brugge, Belgium \move@AU\move@AF\@affiliationBundesdeutsche Arbeitsgemeinschaft für Veränderliche Sterne e.V. (BAV), Berlin, Germany \move@AU\move@AF\@affiliationCBA Mol, Belgium
0000-0002-9810-0506]Gordon Myers \move@AU\move@AF\@affiliationCenter for Backyard Astrophysics, San Mateo \move@AU\move@AF\@affiliationAAVSO, 5 Inverness Way, Hillsborough, CA, USA
The cataclysmic variable AR Scorpii (AR Sco) shows large-amplitude, highly periodic pulsations across the electromagnetic spectrum every 1.97 minutes, superimposed upon a strong waveform at the system’s 3.56-h orbital period (Marsh et al., 2016). The system’s low X-ray luminosity rules out the presence of significant accretion by the white dwarf (WD) primary from its M-dwarf companion, so neither of these signals is powered by accretion (Marsh et al., 2016; Takata et al., 2018). Instead, AR Sco has been called a white-dwarf pulsar because its pulsations consist of synchrotron radiation and are apparently powered by the spin-down of its highly magnetized ( 500 MG) WD, similar to neutron-star pulsars (Marsh et al., 2016; Buckley et al., 2017).
The spin-down of the WD is a foundational conclusion from Marsh et al. (2016), who calculated that the spin-down rate they detected is large enough to power AR Sco’s pulsations. However, Potter & Buckley (2018) contested the significance of this spin-down after finding that the Marsh et al. (2016) spin-down ephemeris did not accurately predict the frequencies of the spin and orbital periods in their optical photometry. Although Potter & Buckley (2018) concluded that a linear spin ephemeris accurately described their data, they were also careful to note that their result constrained—but did not rule out—the WD spin-down. An unambiguous detection of the slowing spin rate, they wrote, would require additional observations.
AR Sco’s light curve contains a number of remarkable features at different timescales. The orbital waveform is brightest at phase 0.4 and has a peak-to-peak amplitude of 1.5-2 mag in the optical, depending on the bandpass. Katz (2017) proposes two alternative models to explain the orbital modulation and why it does not peak at superior conjunction. To provide observational constraints for these models, Littlefield et al. (2017) analyzed archival, ground-based photometry as well as 79 days of continuous photometry by the Kepler K2 mission. They reported that while the system’s overall brightness has remained relatively stable since 2005, the orbital waveform peaked at a different phase and had a slightly lower amplitude between 2005-2007. Moreover, the K2 photometry showed aperiodic brightness fluctuations at the level of a few percent on a timescale of days (Littlefield et al., 2017).
The 1.97-min pulsations are arguably AR Sco’s defining observational characteristic and are remarkable for their speed, amplitude (a factor of in the optical), and phase coherence across a wide range of wavelengths, including radio, near-infrared, optical, ultraviolet (Marsh et al., 2016; Stanway et al., 2018) and even soft X-rays (Takata et al., 2018). Their period corresponds with the beat period (i.e., the orbital sideband) between the binary orbital period and 1.95-min WD spin period, and they are thought to originate on the inner hemisphere of the M5-class companion star (Marsh et al., 2016; Takata et al., 2018). Geng et al. (2016) proposes that the WD’s magnetic axis is inclined with respect to its rotational axis and that the pulses are caused by the interaction of the WD’s magnetosphere with the secondary’s wind.
The 1.95-min spin period of the white dwarf in AR Sco is extremely short when compared with the system’s orbit. White dwarfs are not born spinning so rapidly, and it is thought that a phase of high accretion powered AR Sco’s rapid spin-up, followed by the current epoch of little or no mass transfer (Marsh et al., 2016; Buckley et al., 2017).
Here, we present high-cadence photometry with the twin objectives of (1) searching for a spin-down and (2) disentangling the spin and beat periods.
2.1 SLKT observations
We obtained 39 hours of high-time-resolution photometry of AR Sco using the 80-cm Sarah L. Krizmanich Telescope (SLKT) and an unfiltered Santa Barbara Instrument Group STL-1001 CCD camera at the University of Notre Dame in 2016, 2017, and 2018. As indicated in Table 2.2, which lists details of each time series, the exposure time was 2 s, and factoring in the overhead between images, the typical cadence was 5-6 s, and each time series usually spanned 1-3 h. Fig. 1 plots a representative, 2.5-hour-long light curve and zooms in on one segment during which the pulsations were especially prominent.
The system’s optical pulsations are extremely fast and well-defined, so CCD data must be well-timed in order to be useful. Consequently, we synchronized the clock of the CCD control computer to an atomic clock prior to the start of each time series. Additionally, we measured the shutter lag of the detector (the time offset between the shutter actuation and the timestamp recorded in the FITS header), found it to be stable from night-to-night, and applied an appropriate correction to the image timestamps. Finally, we applied a BJD correction to all observations using routines in Astropy (Astropy Collaboration et al., 2013).
We used aperture photometry to extract AR Sco’s light curve, but there were no optimal comparison stars within the field of view. According to APASS photometry, each of the field stars is quite red, probably because there is a dark nebula associated with the Ophiuchi complex along the line of sight to AR Sco. Faced with this paucity of choices, we selected UCAC4 336-082341 ( = 16h 22m 06.481s, 53’ 27.84”, ) as our comparison star. Since the spectral properties of such a reddened star are likely a poor match for those of AR Sco, any attempt to infer standard magnitudes from the unfiltered photometry would probably suffer from serious systematic errors. Consequently, we do not place our photometry on a standard magnitude scale.
2.2 AAVSO observations
As indicated in Table 2.2, coauthors FJH and GM observed AR Sco and submitted their observations to the AAVSO International Database333https://www.aavso.org under AAVSO observer codes HMB and MGW, respectively. Their observations had cadences that ranged from 14-35 s.
Table 0. \Hy@raisedlink\hyper@@anchor\@currentHrefLog of SLKT observations
|UT Start Date||Length (hr)||Cadence (s)|
Table 0. \Hy@raisedlink\hyper@@anchor\@currentHrefLog of AAVSO observations. The “Obs.” column contains the AAVSO code of the observer.
|Observer||UT Start Date||Length (hr)||Cadence (s)|
Lomb-Scargle power spectra of the SLKT data, presented in Fig. 2.1, show that the measured spin and beat frequencies agree with those reported in Potter & Buckley (2018). However, the predicted spin and beat frequencies from the Marsh et al. (2016) spin-down ephemeris are in poor agreement with the peaks observed in our data, confirming the result in Potter & Buckley (2018). We explore this issue in depth in Sec. 4.
After eliminating a small number of SLKT observations with a signal-to-noise ratio of less than 5, we display our photometry in Fig. 2.1. The data have been phased to the orbital period based on the ephemeris in Marsh et al. (2016). The zero point of the orbital phase is defined as the moment that the red secondary star is at inferior conjunction. The overall light curve shape is similar to the slower cadence data analyzed by Littlefield et al. (2017). The light curve is asymmetric, with a quick rise to maximum brightness at an orbital phase of 0.4 and a second peak at an orbital phase of 0.75.
We assume that the observed light curve can be represented as the sum of a slowly varying orbital modulation and the high-frequency, pulsed emission. To remove the variations associated with the orbit, we divided the orbit into phase bins and identified the faintest 5% of the points in each bin. Because contamination by the pulsed emission inflates the amplitude of the orbital waveform in a simple Lomb-Scargle power spectrum, the faintest observations in these bins more accurately describe the underlying orbital modulation, and we selected the exact threshold from trial-and-error. We then represented these points with a Fourier series so that we could predict the strength of the orbital modulation as a function of orbital phase. As seen in Fig. 2.1, the orbital variation is moderately asymmetric, with the peak brightness just before phase 0.5 and the minimum at phase 0.0.
Subtracting the orbital-modulation function from the phased orbital light curve yields the pulsed light curve as a function of orbital phase (bottom panel in Fig. 2.1). The pulse amplitude displays a strong orbital-phase dependence, peaking at orbital phase 0.25, reaching minimum amplitude near phase 0.6, and rebounding around phase 0.75.
To more closely investigate the changes in the pulse amplitude and morphology, we split the data into ten equally-sized, non-overlapping orbital bins and phased each bin to the beat period. The resulting plots (Fig. 2.2) show that the beat pulse has two unequal maxima per cycle, as has been noted by Marsh et al. (2016). Our beat phaseplots show that the shape, amplitude, and phase of both maxima vary as a function of orbital phase. At orbital phase 0.0, the major pulse is broad and symmetric, but between orbital phases 0.1-0.3, its shape changes, with the rise to beat maximum becoming longer than the decline. This is also the orbital phase when the major pulse has the largest amplitude. Around orbital phase 0.35, the major peak is again symmetric but soon skews in the other direction. Eventually, the major peak nearly disappears around phase 0.55. For the second half of the orbit, the major pulse remains symmetrical, but is weaker than in the first half of the orbit. The minor pulse reaches its largest amplitude between orbital phases 0.3 and 0.4, significantly later than phase 0.25 for the major pulse.
The beat pulse undergoes a significant phase shift, depending on the orbital phase. This is most easily seen in a heatmap of the pulse brightness versus both the orbital and beat phase (see Fig. 3). This figure clearly shows that between orbital phases 0.1 and 0.5, the major peak shifts by 10% in beat phase. The shift in the beat pulse over the second half of the orbit is smaller than the first half, but the pulse is generally seen to be broader and fainter. Takata et al. (2018) reported similar findings from their analysis of 39 ks of data obtained on 2016 September 19 with the XMM Newton satellite’s Optical/UV Monitor Telescope, and our results imply that this behavior is stable on timescales of years.
The amplitude and phase shifts of the beat pulse are consistent with the addition of second periodic signal with a slightly different frequency (e.g., as observed in FO Aqr; Osborne & Mukai, 1989). For example, if the spin pulse were to be isolated and plotted in the heat map in Fig. 3, it would run diagonally since the figure phases the data to the beat period. To explore this possibility, we modeled the full light curve as the superposition of three periodic signals: the spin, beat, and orbital periods (similar to the light curve model for FO Aqr in its low state (Littlefield et al., 2016)). We found the best-fit trigonometric function at each of those frequencies by a simple least-squares fit, with each term consisting of three harmonics. We also attempted to fit additional frequencies detected in the power spectrum, but adding these terms did not significantly improve the quality of the fit.
The results of this fit (Fig. 3) reveal that the beat and spin models are both double-peaked, with the beat pulse having two unequal maxima. The maxima of the spin pulse, by contrast, are roughly equal. Our model predicts that the amplitude of the spin pulse is about 50% of the amplitude of the beat pulse, while the minor pulses are comparable in brightness. However, one limitation of our model is that it underestimates the amplitude of the highest-amplitude beat pulses. This might be a consequence of the fact that the model does not account for the changing visibility of the secondary’s inner hemisphere across the orbital period.
As a test of this model, we computed the OC for the blended spin-beat pulse, and as shown in Fig. 3, it offers a reasonably accurate prediction of the actual OC. The simulated OC includes a correction for light-travel delays caused by the secondary’s orbital motion. The light-travel delay for pulsations originating from material orbiting in the plane of the binary is given by
where is the distance of the emission from the binary center of mass, is the orbital inclination, is the speed of light, and is orbital phase of inferior conjunction for the emitting material. Given the requirement that the emission originate on the donor (Marsh et al., 2016), the smallest-possible light-travel delay would occur for emission arising at the first Lagrangian point (L). If we adopt M = 0.8 M and M = 0.3 M, as did Marsh et al. (2016), and assume an orbital inclination of , then the semi-amplitude of the light-travel delay would be 0.6 s.
Although the effect of light-travel delays is relatively minor at our 5-s cadence, a higher time resolution might potentially be able to discern light-travel delays from emission from different regions of the secondary (i.e., across a range of values of and ).
Although there might be a light-travel delay associated with the spin pulse, the region in which the optical spin pulse is generated is unknown, making it impossible to calculate its light-travel delay. Buckley et al. (2017) provide evidence that it might originate in a “striped wind” outside the light cylinder of the magnetosphere, which is almost an order of magnitude larger than the binary orbital separation.
4 Spin-down ephemeris
4.1 Measuring the spin-down
Marsh et al. (2016) detected a significant spin-frequency derivative of Hz s, implying a spin-down of sufficient magnitude to power the optical pulsations. However, Potter & Buckley (2018) found that the Marsh et al. (2016) spin-down ephemeris did not accurately predict the frequencies observed in their power spectra and concluded their photometry was consistent with a constant spin period of 0.008538220(3) Hz, where the number in parentheses is the uncertainty on the final digit. While Potter & Buckley (2018) showed that the spin-down ephemeris from Marsh et al. (2016) was inaccurate, they noted that their results could still be reconciled with a nonzero and that a longer baseline of observations was necessary to investigate this possibility.
Our results in Sec. 3 show that the beat pulse is more readily measured than the lower-amplitude spin pulse, so we use OC measurements of the beat pulse to search for a change in the spin period. Assuming that any change in the orbital period is small over the baseline of observations, the derivative in the beat frequency will be a direct measure of the derivative of the spin frequency.
We measured 1,077 beat-pulse timings444These timings are available as an online table, and in Table 4.1, we provide a sample of them to illustrate the format of the data. from our dataset by fitting a Gaussian to each well-observed beat pulse. The timings, which span two years and three observing seasons, have a sufficiently long baseline to search for . We applied a correction to all pulse timings to compensate for the orbital-phase dependence of their arrival times, using the empirical fit from Fig. 3. We then tested the Potter & Buckley (2018) linear ephemeris by computing an OC diagram, using the beat period implied by their spin period. To improve the signal-to-noise ratio, we averaged the OC residuals for each night and used the standard error of their mean as the 1 uncertainty for each night. The residuals in this plot (Fig. 4.1) show a rising trend, suggesting that the beat period inferred from Potter & Buckley (2018) is not a good match to our data. Further, there appears a significant curvature in the measurements consistent with the presence of a period derivative.
Table 0. \Hy@raisedlink\hyper@@anchor\@currentHrefBeat-pulse timings. The full table is available online as a machine-readable table.
Relative to Eq. 2. Raw pulse timings, uncorrected for orbital-phase dependence of arrival times. Pulse timings corrected for orbital phase.
We employed two independent fitting procedures—a bootstrap fit and an affine-invariant Markov Chain Monte Carlo algorithm (Foreman-Mackey et al., 2013)—to generate linear and quadratic ephemerides for the pulse maxima. We inspected the residuals from the linear and quadratic ephemerides (Fig. 4.1) to determine which best described the data. While the measurements for each season systematically evade the linear fit, this trend vanished with the inclusion of a quadratic term, and the statistic dropped from 13.1 to 2.6. Based on the quadratic fit, we calculate a beat ephemeris of
The quadratic coefficient is equivalent to where is the average beat period, yielding a unitless period derivative of To convert this to a frequency derivative, we start with the definition and differentiate with respect to , obtaining The time derivative of this relation is
Thus, our is equivalent to Hz s.
4.2 Testing the spin-down with AAVSO photometry
Fig. 4.1 includes OC values from beat pulses in AAVSO photometry, but in order to provide an independent test of the spin-down ephemeris, these values were not included in the calculation of Eq. 2. Because the cadences of these time series were too slow for OC analysis of individual beat maxima, we took each AAVSO light curve and phased it to the beat period using our linear and quadratic ephemerides. The resulting beat-phase plots, an example of which is shown in Fig. 4.1, can be used to measure the average phase of the beat-pulse maxima provided that enough beat cycles were observed. We filtered the AAVSO photometry to exclude any time series with a cadence slower than 35 s or a duration shorter than 1 hr, and we applied a correction for the orbital-phase dependence of the pulse arrival times. We further excluded any time series afflicted by aliasing between the beat period and the observing cadence.
With the benefit of the extended baseline of observations, it is immediately obvious in Fig. 4.1 that the linear ephemeris leads to strong curvature in the residuals. By contrast, the quadratic residuals do not show any systematic trend.
No shutter-lag correction was applied to the AAVSO data, and the raw AAVSO quadratic residuals showed an offset of -0.7 s with respect to contemporaneous SLKT residuals. An uncorrected shutter lag will cause timestamps to be earlier than the data that they describe, imparting a small, negative OC to the pulse timings. We therefore attribute the constant, negative offset of the AAVSO residuals to uncorrected shutter lag, and in Fig. 4.1, we add an offset of 0.7 s to each AAVSO residual to compensate for this effect. Regardless of the cause of this offset, the lack of a systematic trend in the AAVSO quadratic residuals supports our measurement of the spin-down.
4.3 The nature of
Because the beat frequency is defined as , its time derivative is , meaning that the observed could be caused either by a decrease of the spin frequency or an increase in the orbital frequency. If the orbital period were changing so rapidly, there would be detectable consequences in long-term photometry. To explore this possibility, we used phase-dispersion minimization (PDM; Stellingwerf, 1978) to determine whether an increasing orbital frequency could reduce the scatter in orbital phase plots using 13 years of survey photometry from the Catalina Real-Time Transient Survey (CRTS; Drake et al., 2009) and All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al., 2014; Kochanek et al., 2017). For a range of between Hz s and Hz s, we phased the survey photometry using a quadratic orbital ephemeris in which we assumed an initial orbital period of d (Marsh et al., 2016) at an epoch of JD = 2455000.5 (Potter & Buckley, 2018).555This is the epoch of the CRTS photometry used by Marsh et al. (2016) to measure the orbital period. It is different than the epoch in the orbital ephemeris in Marsh et al. (2016), which reported a spectroscopically determined epoch and assumed—correctly—that the orbital period did not change appreciably in the intervening 5 years since the CRTS epoch. Phase plots using the CRTS epoch will have a uniform horizontal offset with respect to the Marsh et al. (2016) epoch, but this does not impact our PDM analysis. For each quadratic phase plot, we computed the PDM statistic and normalized it to the the PDM statistic for the phase plot constructed with the linear ephemeris from Marsh et al. (2016). The results of this analysis are shown in Fig. 4.3. For Hz s, the phase plots contain an obvious increased scatter with respect to the linear ephemeris; below that value, the effect of becomes negligible, so we constrain Hz s. To verify this constraint, we created a simulated orbital light curve of AR Sco, injected an artificial , and matched the sampling of the simulated data to the actual observations. As Fig. 4.3 shows, our algorithm successfully recovered values of in excess of Hz s, but below that threshold, it could not discern the simulated , bolstering our constraint.
The non-detection of is consistent with the measurements of the orbital phase of maximum light from Littlefield et al. (2017). An increasing orbital frequency (i.e., ) would induce concave-down curvature in their Fig. 4, but this is not seen. We conclude, therefore, that any contributes negligibly to , such that to within our measurement uncertainty. Thus, the subscript notation becomes unnecessary for .
The disagreement between the Marsh et al. (2016) spin-down ephemeris and the measured spin frequency in Potter & Buckley (2018) provides a stringent test of our spin-down ephemeris. In Fig. 4.4, we plot the extrapolated beat period from Eq. 2 as a function of beat-cycle count, overlaying the beat periods from both Marsh et al. (2016) and Potter & Buckley (2018) at the appropriate epochs.666Potter & Buckley (2018) reported only , so we used from Marsh et al. (2016) to calculate the corresponding Even though our beat ephemeris was calculated without regard to either of these two measurements, the projected beat period is in excellent agreement with both measurements. Furthermore, the difference between the Potter & Buckley (2018) and Marsh et al. (2016) spin frequencies, when divided by the difference in epochs yields Hz s, consistent with both our measurement ( Hz s) and the constraints from Potter & Buckley (2018).
This suggests that the inability of the Marsh et al. (2016) ephemeris to correctly predict the spin and beat frequencies in Potter & Buckley (2018) is a consequence of an underestimated and can be remedied by using our measurement of the spin-down. Despite the error in their estimate of , it is remarkable, given the sparse sampling and comparatively low time resolution of the survey photometry available to them, that they were able to accurately measure , , and while also estimating to within a factor of 2.
4.5 Spin-down luminosity
Our precise measurement of the frequency decay rate allows us to improve the estimate of the spin power available for conversion into the observed electromagnetic (EM) energy. The spin-down luminosity, is given by where is the WD’s moment of inertia (Marsh et al., 2016). As did Marsh et al. (2016), we assume a 0.8 M WD with a radius of 0.01 R. The mass-radius relation for non-relativistic WD stars means that the moment of inertia changes as , and is rather insensitive to variations in the assumed mass. In the non-relativistic regime, WDs are predicted to have a density structure like that of a polytrope with an index of 1.5, and we calculate a moment of inertia of kgm. However, the strong magnetic fields and the rapid spin rate of the WD may have an effect on the precise value of the moment of inertia (Franzon & Schramm, 2017).
Thus, we find the spin-down power is W. The Gaia DR2 parallax (Gaia Collaboration et al., 2016, 2018) to AR Sco ( mas) is equivalent to a distance of pc, which is close enough to the Marsh et al. (2016) value that it is unnecessary to correct their measurement of the average EM power of the pulsations. We find the efficiency of converting the spin energy to detected EM emission to be 4%.
The complex morphological variations seen in the optical light curve of AR Sco are explained as the superposition of the spin and beat pulses, as well as a slower orbital modulation. As noted in earlier studies, a beat cycle consists of a major and minor pulse. Here, we show the minor pulse has half the amplitude of the brighter signal after removing contamination from the spin pulse. We find the spin cycle also consists of two distinct pulses, but their amplitudes are comparable. The major spin pulse is half the amplitude of the major beat pulse and therefore, about the same amplitude as the minor beat pulse.
Over an orbit, the major beat and spin pulses add constructively between phases 0.2 and 0.3, resulting in the highest-amplitude optical variations observed in the system. Half an orbit later, the major beat and minor spin pulses add together along with the minor beat and major spin pulses. This combination results in lower amplitude peaks when compared with orbital phases around 0.25. The smallest amplitudes are seen when the beat and spin pulses are out of phase and destructively interfere around orbital phases between 0.5 and 0.6. This model provides a good overall fit to the rapid variations seen in AR Sco. We also show that this model explains the OC variations seen in the timings of the beat pulses.
Our simple model assumes constant beat and spin amplitudes over an orbit, and this does not fully match the largest amplitudes fluctuations observed around orbital phases 0.2 to 0.4. A more complete model would take into account the changing viewing angle of the secondary star thought to be the origin of the beat emission.
The major beat and spin pulses coincide at orbital phase 0.3, resulting in the highest amplitude peaks in the AR Sco light curve. The beat pulse likely comes from near the surface of the red secondary star (Marsh et al., 2016), while the spin pulse is likely originating in the magnetic field of the WD. If we assume that beat maxima occur when one of the WD magnetic poles is pointing toward the secondary, then the geometry of the system suggests that we see the peak spin pulse when the WD magnetic pole is nearly perpendicular to our line-of-sight.
Perhaps most importantly, our results establish that the WD is indeed spinning down (as originally proposed by Marsh et al., 2016), but the frequency derivative that we measure is almost twice as large as their estimate. Our updated spin-down ephemeris successfully passes two tests: it accurately predicts the evolution of the beat period between the Marsh et al. (2016) and Potter & Buckley (2018) epochs, and it also accounts for the pulse-arrival times in AAVSO photometry. Furthermore, it comports with the constraints on the spin-down rate established in Potter & Buckley (2018). Our measurement of the spin-down of the WD confirms the conclusion by Marsh et al. (2016) that the pulsed EM emission from AR Sco can be powered by the rotational energy of the WD.
The Sarah L. Krizmanich Telescope is a generous donation by the Krizmanich family to the University of Notre Dame. It is named in honor of their daughter.
We thank David Buckley for showing us several figures from an early draft of the Potter & Buckley (2018) manuscript. We avoided using any of this advance information in our paper.
We are grateful to our anonymous referee for urging us to measure the spin-down rate of the white dwarf using the SLKT data.
This research uses observations of AR Sco from the AAVSO International Database contributed by GM and FJH, both of whom also participate in the Center for Backyard Astrophysics collaboration (Patterson, 2012).
- Astropy Collaboration et al. (2013) Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33
- Buckley et al. (2017) Buckley, D. A. H., Meintjes, P. J., Potter, S. B., Marsh, T. R., & Gänsicke, B. T. 2017, Nature Astronomy, 1, 29
- Drake et al. (2009) Drake, A. J., Djorgovski, S. G., Mahabal, A., et al. 2009, ApJ, 696, 870
- Foreman-Mackey et al. (2013) Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306
- Franzon & Schramm (2017) Franzon, B. & Schramm, S. 2017, MNRAS, 467, 4484
- Gaia Collaboration et al. (2016) Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1
- Gaia Collaboration et al. (2018) Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, arXiv:1804.09365
- Geng et al. (2016) Geng, J.-J., Zhang, B., & Huang, Y.-F. 2016, ApJL, 831, L10
- Katz (2017) Katz, J. I. 2017, ApJ, 835, 15
- Kochanek et al. (2017) Kochanek, C. S., Shappee, B. J., Stanek, K. Z., et al. 2017, PASP, 129, 104502
- Littlefield et al. (2016) Littlefield, C., Garnavich, P., Kennedy, M. R., et al. 2016, ApJ, 833, 93
- Littlefield et al. (2017) Littlefield, C., Garnavich, P., Kennedy, M., et al. 2017, ApJL, 845, L7
- Osborne & Mukai (1989) Osborne, J., & Mukai, K. 1989, MNRAS, 238, 1233
- Patterson (2012) Patterson, J. 2012, JAAVSO, 40, 240
- Potter & Buckley (2018) Potter, S. B. & Buckley, D. A. H. 2018, MNRAS, 478, L78
- Marsh et al. (2016) Marsh, T. R., Gänsicke, B. T., Hümmerich, S., et al. 2016, Nature, 537, 374
- Shappee et al. (2014) Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48.
- Stanway et al. (2018) Stanway, E. R., Marsh, T. R., Chote, P., et al. 2018, arXiv:1801.07258
- Stellingwerf (1978) Stellingwerf, R. F. 1978, ApJ, 224, 953
- Takata et al. (2018) Takata, J., Hu, C.-P., Lin, L. C. C., et al. 2018, ApJ, 853, 106