Issue 
A&A
Volume 599, March 2017



Article Number  A116  
Number of page(s)  8  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201629797  
Published online  09 March 2017 
Nonradial and radial period changes of the δ Scuti star 4 CVn
II. Systematic behavior over 40 years
^{1} Institut für Astrophysik der Universität Wien, Türkenschanzstr. 17, 1180 Wien, Austria
email: breger@astro.as.utexas.edu
^{2} Department of Astronomy, University of Texas, Austin, TX 78712, USA
^{3} N. Copernicus Astronomical Center, Bartycka 18, 00716 Warsaw, Poland
Received: 27 September 2016
Accepted: 12 December 2016
Aims. Radial and nonradial pulsators on and near the main sequence show period and amplitude changes that are too large to be the product of stellar evolution. The multiperiodic δ Sct stars are well suited to study this, as the period changes of different modes excited in the same star can be compared. This requires a very large amount of photometric data covering years and decades as well as mode identifications.
Methods. We have examined over 800 nights of highprecision photometry of the multiperiodic pulsator 4 CVn obtained from 1966 through 2012. Because most of the data were obtained in adjacent observing seasons, it is possible to derive very accurate period values for a number of the excited pulsation modes and to study their systematic changes from 1974 to 2012.
Results. Most pulsation modes show systematic significant period and amplitude changes on a timescale of decades. For the wellstudied modes, around 1986 a general reversal of the directions of both the positive and negative period changes occurred. Furthermore, the period changes between the different modes are strongly correlated, although they differ in size and sign. For the modes with known values of the spherical degree and azimuthal order, we find a correlation between the direction of the period changes and the identified azimuthal order, m. The associated amplitude changes generally have similar timescales of years or decades, but show little systematic or correlated behavior from mode to mode.
Conclusions. A natural explanation for the opposite behavior of the prograde and retrograde modes is that their period changes are driven by a changing rotation profile. The changes in the rotation profile could in turn be driven by processes, perhaps the pulsations themselves, that redistribute angular momentum within the star. In general, different modes have different rotation kernels, so this will produce period shifts of varying magnitude for different modes.
Key words: stars: variables: δ Scuti / stars: oscillations / stars: individual: 4 CVn / techniques: photometric
© ESO, 2017
1. Introduction
Period and amplitude changes are common among the pulsating δ Sct stars. These changes are not caused by stellar evolution, for example, Breger & Pamyatnykh (1998) showed that for Pop. I radial pulsators, the measured (1 /P)dP/ dt values are around 10^{7} yr^{1} with equal distribution between period increases and decreases. The evolutionary models, on the other hand, predict that the vast majority should show increasing periods with values a factor of ten smaller than observed. For nonradial δ Sct pulsators of Pop. I, the discrepancies are even larger. The behavior suggests that for these relatively unevolved stars the rate of evolution cannot be deduced from the period changes. Since that study, additional attempts to match the observed changes to evolutionary models have also failed. A recent example is that of Murphy et al. (2012), who attempt to explain the observed amplitude changes in the star KIC 3429637 (ρ Puppis subclass) by evolutionary changes in the pulsation driving.
It is, therefore, essential to study the period and amplitude changes in detail, in order to separate the evolutionary from the other effects. These “other effects” are at present poorly understood and presumably temporary changes in the structure of the star. Such investigations require extensive data covering decades. We refer to promising results for the cepheid Polaris, where the longterm changes of 4.5 s yr^{1}, covering 150 yr, are evolutionary, while smaller changes are superposed, for example, around 1965 (Turner et al. 2005).
In δ Sct stars, such studies are made more difficult by the presence of hundreds of simultaneously excited pulsation modes of generally small amplitudes. Even when the studies concentrate on the largeamplitude modes only, these stars often show close frequencies from modes of different spherical degree, ℓ, caused by mode trapping in the envelope (Breger et al. 2009). However, the presence of radial as well as nonradial modes can also be advantageous to distinguish between different astrophysical explanations. Pulsation modes of different degree probe different parts of the star and their period changes show changes in the stellar structure. Changes in stellar radius and/or rotation lead to changes in the observed rotational splitting and thereby a dependence of the period changes on the azimuthal order, m, of a pulsation mode. Furthermore, magnetic cycles also lead to period changes.
Photometric observations used.
One of the most promising and beststudied stars is the evolved δ Sct variable 4 CVn (HR 4715 = HD 107904 = AI CVn, F3IIIIV). Its variability was discovered by Jones & Haslam (1966). This led to a large number of photometric studies, of which some were extensive multisite campaigns. Most of the early photometry obtained prior to 2000 have been summarized in Table 1 of Breger (2000). Since then, additional photometry from 1991 was made available by R. J. Dukes (Breger et al. 2008).
The largest photometric study of 4 CVn was obtained in two different passbands using a dedicated automatic photoelectric telescope. It covered 702 nights from 2005 through 2012 (Breger 2016, hereafter Paper I). 64 pulsation frequencies could be detected at a high significance level. During these eight years, the periods of the wellstudied dominant nonradial modes were found to be variable with rates of change ranging from (1 /P)dP/ dt = −16 × 10^{6} to + 13 × 10^{6} yr^{1}. These values are one or two orders of magnitude larger than expected from stellar evolution. Furthermore, both increasing and decreasing periods were found. However, for each of the more dominant modes studied, the rates of period changes were constant from month to month as well as from year to year.
The question arises as to whether these rates of period change are constant over longer time scales than the previously studied eight years. Due to the large sizes of the observed period changes, we would expect sign reversals at some point. Do such reversals occur simultaneously in all pulsation modes, irrespective of the sign of the period changes? Such a situation would suggest a change in the internal structure of the star as the origin of the period changes.
The available photometric data for 4 CVn cover four decades. In this paper we examine the period changes of the main pulsation modes of 4 CVn over these decades and correlates the result with the type of pulsation mode.
2. Photometric frequency determinations
In this investigation we have examined systematic frequency changes on the order of 10^{5} yr^{1} or smaller. This means that the data covering many decades need to subdivided into a number of separate time groups, each of which is suitable for precise multifrequency solutions. Each group needs to contain a large number of highprecision measurements to yield accurate multifrequency solutions and span a sufficiently long timebase for high frequency resolution. These requirements suggest that each group should contain two successive observing seasons, if possible.
We are fortunate that more than 800 nights of observation are available and that they naturally divide into several groups. However, a small amount of data could not be used for the present analysis. Earlier photometry covering 1966–1970 (see Table 1 in Breger 2000) was not extensive enough to determine the frequencies of 4 CVn with sufficient accuracy. Furthermore, some additional data from 1997 (Stankov et al. 2000) were not included because of relatively lower internal consistency. We note though that inclusion of these data would not change the conclusions of the present paper.
Table 1 lists the photometric data used in the present analysis. The pulsation frequency determinations were performed with a package of computer programs (PERIOD04, Lenz & Breger 2005). PERIOD04 utilizes both Fourier and multipleleastsquares algorithms, which do not rely on sequential prewhitening. It is especially suited for extensive data containing large numbers of simultaneously excited frequencies.
In order to examine the period changes in stars, the standard (O−C) phaseshift technique is usually applied. As part of this technique, an optimum frequency is determined. This frequency is then used to derive the differences between the actually observed and calculated phases as a function of time. However, for 4 CVn the (O−C) plots revealed the common difficulty with the phase shift method: the unknown number of elapsed cycles between the different groups of years with sufficient data. While the photometric coverage of 4 CVn is excellent for a large number of individual groups of years, the data also show large gaps without coverage between these groups of years.
Consequently, for the present analysis we have not used the (O−C) phaseshift method, but actually determine the optimum values of the main frequencies for each individual data set. This standalone approach also eliminates the potential wavelengthdependent phaseshift problem of the extensive 1974 data, which were obtained through the Strömgren b bandpass, as opposed to the Johnson V or Strömgren y bandpasses used for the other data sets.
The quantity and spacing of these data prior to the year 2000 immediately suggest the following five groups: the 1974 data, combined 1976–1978 data, multisite 1983 and 1984 campaigns, the 1991 and 1992 data, and the 1996 multisite campaign combined with the extensive 1997 singlesite data. For the 1976–1978 solution, one frequency value was unreliable and had to be omitted: during 1976, the amplitude of the 7.376 cycle d^{1} mode decreased to become essentially undetectable before it started to increase its amplitude again. For the 2005–2012 seasons we have adopted the following approach to obtain the desired frequency resolution: for each solution we combined two adjacent observing seasons, covering approximately 18 months. The eight observing seasons from 2005 to 2012 therefore lead to seven solutions.
To determine the period changes, we have calculated new multifrequency fits for these data sets. Additional details on the multiple data sets and their properties can be found in the references listed in Table 1.
Schmid et al. (2014) also discovered that 4 CVn is an eccentric binary system with an orbital period of 124.4 ± 0.03 d. This raises the question whether the measured pulsation periods and period changes are seriously affected by the lighttime effects originating in the binary orbit. Consequently, we repeated our analyses without the orbital corrections. No significant changes in the derived frequencies and amplitudes presented in this paper were found, although the residuals of the analysis without the orbital lighttime corrections were slightly higher. Also, the orbital period can be deduced from an analysis of the small residuals of the uncorrected photometry. This remarkable effect is examined in more detail in a later section. We emphasize that ignoring the binary nature would not change the results and conclusions of the present paper. The numbers presented in this paper are based on the data after applying orbital corrections.
3. Pulsation mode identifications for 4 CVn
In order to understand observed intrinsic period changes, it is important to identify the individual observed pulsation modes in terms of the pulsational degree, ℓ, and the azimuthal order, m. For δ Sct stars, spectroscopic and photometric methods complement each other for mode identification. In a pioneering study, Schmid et al. (2014) analyzed the lineprofile variations of 4 CVn from over 2000 highdispersion spectra obtained from 2008 to 2011. This study was an extension of an earlier spectroscopic project by Castanheira et al. (2008). For seven pulsation modes, Schmid et al. determined the m values.
For A and F stars, the brightness variations measured through two different bandpasses show small phase shifts and different amplitudes, depending on the type of pulsation mode. These phase shifts and amplitude ratios are mainly a function of the spherical degree, ℓ. Consequently, the spectroscopic and photometric mode identification techniques can be combined to determine the spherical degree and azimuthal order for the pulsation modes.
The 702 nights of highaccuracy photometry, obtained from 2005 to 2012 (see Paper I), were obtained with two Strömgren v,y filters. The bandpasses of these filters are ideal to determine the phase shifts and amplitude ratios of stars of spectral type A and F. In order to minimize any effects due to amplitude and frequency changes, we have treated the data from each annual observing season separately.
We determined the values of the phase shifts and amplitude ratios for each observing season from 2005 to 2012. The annual values are in excellent agreement with each other. The results and their uncertainties are listed in Table 2.
Phase differences, amplitude ratios and mode identifications.
Fig. 1 Amplitude ratios and phase differences in the Strömgren v and y passbands for different pulsation modes. The error bars are determined by comparing results from eight different years. In A/F stars the observed phase differences are strong indicators of the spherical degree. The results from a theoretical model (see text) covering the frequency range from 4.5 to 9.5 cycle d^{1} are also shown. 
Fig. 2 Frequency variations of the six dominant modes from 1974 to 2012. The errors bars shown are usually the size of the symbols. The diagram shows the systematic difference between prograde and retrograde modes. Furthermore, around ~1986 the sign of the frequency changes reversed systematically. 
Frequency values changing over four decades.
We can now compare the observed values with those predicted by theoretical models for different pulsation modes. Our calculation of amplitude ratios and phase differences for the Strömgren v and y passbands relies on the formalism presented in DaszyńskaDaskiewicz et al. (2003). An evolved stellar model with a mass of 2.0 solar masses was used as input for our computation. The model parameters log L/L_{⊙} = 1.486, T_{eff} = 6845 K, log g = 3.49 match the effective temperature and luminosity determined by Schmid et al. (2014) very well (i.e., log L/L_{⊙} = 1.47 ± 0.05, T_{eff} = 6875 ± 120 K). Moreover, solar metallicity was assumed.
The comparison between the observed and the corresponding computed phase shifts and amplitude ratios is shown in Fig. 1. Detailed inspection of the figure leads to the following conclusions: while the theoretically predicted phase differences agree well with the observations, the amplitude ratios show a larger deviation. This is a well known problem in δ Sct stars and is caused by the imperfect theoretical description of the outer layers in these rather cool stars. Therefore, only the phase difference is commonly used as a discriminator of the spherical mode degree for δ Sct stars. As we can see from Fig. 1, the modes can be separated into three groups, each corresponding to ℓ = 0, 1 or 2.
We note that the derived ℓ values agree remarkably well with preliminary results given in Lenz et al. (2010), which used only the data from the much more limited 1996 multisite campaign. We have combined the new results with the spectroscopically determined m values. The final mode identifications are also listed in Table 2.
An interesting detail is revealed when comparing the detected azimuthal orders of the modes to their amplitude ratios. The three modes with positive azimuthal order have lower v/y amplitude ratios than their negative counterparts. This suggests an effect of rotation on the amplitude ratios. However, this intriguing result must be regarded as somewhat uncertain due to smallnumber statistics.
4. Frequency variations: 1974–2012
Let us examine the results for the six dominant modes with known mode identifications. Figure 2 and Table 3 show the derived frequencies with their statistical uncertainties for each of the data groups discussed in Sect. 2. We note that, except for the 1974 data, the statistical uncertainties are very small and often smaller than the size of the symbols in Fig. 2.
Fig. 3 Confirmation of the binary orbit from a comparison of the 2005–2012 photometric data without and with the orbital lighttime corrections. Top left panel: Fourier amplitudes of the observed data near the f_{1} frequency. Bottom left panel: Fourier amplitudes of the residuals after prewhitening the best 64frequency solutions. Note the changed amplitude scale. The uncorrected data (dotted curve) shows the side peaks from the 124 (~1/ 0.008) d orbit, while the data corrected for the spectroscopically detected binary orbit eliminates these peaks. The difference of the f_{1} side peaks from the lower left panel is demonstrated in the top right panel. A similar analysis for f_{2} (bottom right panel) confirms the detection. 
The pulsation frequencies of 4 CVn are systematically changing on timescales of decades. The photometric amplitudes of the dominant pulsation modes are also variable and range from 2 to 24 millimag.
While the timescale of the amplitude variations are also long, the amplitude variations showed no obvious patterns, no correlations from mode to mode, and no correlations with the systematic frequency changes. Consequently, they will not be considered further in the present discussion.
The study of the frequency variations of the dominant modes show the following remarkable patterns:

i.
The frequency of the (axisymmetric) radial mode shows verylittle, if any, variation from 1974 to 2012.

ii.
The frequencies of the nonradial modes vary in a continuous manner on a timescale of decades. From 2005 to 2012, the rates of frequency variations were essentially constant.

iii.
For the nonradial modes studied, we found a simultaneous sign reversal of the rates of frequency change around 1986.

iv.
Prograde modes show increasing frequencies before ~1986 and decreasing frequencies after ~1986.

v
Retrograde modes show decreasing frequencies before ~1986 and increasing frequencies after ~1986.

vi.
The (axisymmetric) radial mode shows only small or zero frequency changes before and after ~1986.
5. Can the binary orbit be detected photometrically?
The longterm period changes reported in this paper were calculated with and without orbital lighttime corrections. We found that the derived (longterm) frequency changes presented in this paper are not significantly affected by the binary nature of 4 CVn. This is not surprising since the orbital period is about 1/3 yr, implying that each year covers almost three full orbits. This means that combining data from two adjacent years covers the orbital phases in a similar way. Also, the annual observing seasons cover at least one full orbit for most years. For the 702 nights of photometry presented in Paper I (2005–2012 data), correction of the lighttime effects in the eccentric binary orbit leads to a small reduction of the residuals of the 64frequency solution, providing support for the spectroscopic binary detection.
In a series of recent papers (e.g., Shibahashi & Kurtz 2012; Shibahashi et al. 2015; Murphy et al. 2016), it has been pointed out that the binary nature can also be detected from the observed frequency modulations and phase shifts of the pulsation modes caused by the orbital motion. This has been successfully applied to pulsators measured by the Kepler spacecraft. If the observational time span of each data group is long enough, we expect that the Fourier spectrum of the pulsation modes are triplets with small side lobes separated from the central component by 1/P_{orb}. For 4 CVn we would expect side lobes separated from the intrinsic pulsation frequencies by 0.008 cycles d^{1} and an amplitude near half a millimag for the dominant f_{1} mode.
This does not mean that such side lobes must necessarily be detectable in the groundbased photometry of 4 CVn: the detection of the orbital side lobes is severely affected by the known longterm frequency changes and amplitude variability of the pulsation modes, gaps in the observational coverage, as well as the large number of undetected pulsation modes beyond the adopted 64 frequencies.
In fact, the Fourier frequency spectrum of the observed 2005–2012 data does not clearly reveal the peaks, even after comparison with the spectral window. This is illustrated for the dominant mode, f_{1}, in the top left panel of Fig. 3, where the small orbital side lobes are hidden.
In order to search for such side lobes, we have used the data from the 702 nights of observations without applying the orbital corrections. For each year, an optimum 64frequency solution was calculated and prewhitened. This approach eliminated most the effects of amplitude and longterm period variability, although the annual (as opposed to biannual) solutions may lead to some overinterpretation of this variability. We then calculated a Fourier frequency spectrum from the residuals.
The bottom left panel of Fig. 3 shows the Fourier spectrum of the residuals near the frequency of f_{1}. Two side lobes separated from the prewhitened f_{1} frequency by 0.008 cycles d^{1} are clearly seen. We have repeated the analysis after correcting for the known spectroscopically determined orbit. The top right panel depicts the difference between the two Fourier spectra. The orbital side lobes are clearly seen together with their expected annual side lobes (±0.003 cycles d^{1}). The same result is obtained for f_{2}.
We conclude that the observed photometry confirms the spectroscopic binary detection and the orbital period near 125 (=1/0.008) d.
6. Astrophysical interpretations
In this section we explore a number of possible astrophysical explanations for the observed systematic frequency changes.
6.1. Binary system effects
Could both components of 4 CVn be pulsating with the systematically different changes of the pulsation frequencies corresponding to the two stellar components? This has been seen before: in the δ Sct star θ^{2} Tau (Breger 2005), the observed pulsation modes divide into two groups with opposite signs of period changes. For θ^{2} Tau the binary orbit is known. The predicted orbital lighttime effects of the brighter component fit the observed period changes of the modes with frequencies smaller than 20 cycles d^{1}, while the fainter component fit the changes of the frequencies higher than 20 cycles d^{1}. In principle, for 4 CVn such a scenario might even explain the general reversal of all the signs of period change around 1986, which would be caused by the geometry of the orbital motion of both components. However, in 4 CVn the known orbital period of 124 d does not fit the timescale of the observed frequency changes covering decades. Furthermore, in the binarysystem model, all pulsation modes should show similar lighttime shifts. The opposite is observed. Consequently, we have to reject the binarysystem explanation for the observed period changes.
6.2. Nonlinear mode interaction
Another theoretical explanation for large period changes is nonlinear mode interaction, for example, as proposed for XX Pyx by Handler et al. (2000). Such mode interactions between three modes were described by Moskalik (1985). As noted by Blake et al. (2003), the large number of modes excited in most δ Sct stars makes theoretical predictions from models of such mode interactions difficult. In addition, it does not seem to naturally provide an explanation for the opposite behavior of prograde and retrograde modes.
6.3. Rotation
In the limit of slow rotation, the frequencies of stellar pulsation modes are split according to the following formula: (1)where n is the radial overtone number, l is the order, and m is the azimuthal order of the mode; f_{rot} = Ω/2π is the solidbody rotation rate, and C is the Coriolis coefficient, which is usually much less than 1 for p modes. If there is differential rotation, then the term for the splitting is replaced by an integral over the star of the local rotation rate, Ω(r,θ), times the rotation kernel of a given mode. As Eq. (1) illustrates, changes in Ω would produce opposite effects for prograde (m > 0) and retrograde (m < 0) modes. Thus, this interpretation could offer a natural explanation for these changes.
6.3.1. Change in rotation due to a radius change
As the simplest possible model, let us assume that the frequency changes are driven by a change in the radius of the star. Assuming solid body rotation and conservation of angular momentum, this implies that R^{2}f_{rot} is constant, so a nonradial mode will have its frequency, f_{NR}, perturbed by (2)The frequency of a radial mode, f_{R}, is ∝(G ⟨ ρ ⟩ )^{1/2}. Consequently, we find that (3)Thus, we expect that (4)Since the measured rotational velocity for 4 CVn is greater than 100 km s^{1}, the rotation period is of the order of one day or smaller. Using f_{rot} ~ 1 d^{1} and f_{R} ~ 7 d^{1} we find that Δf_{NR}/ Δf_{R} ~ 0.1. However, the measured values of Δf_{NR} and Δf_{R} are ~6 × 10^{5} d^{1} and ~0.5 × 10^{5} d^{1}, respectively, so empirically we have Δf_{NR}/ Δf_{R} ~ 10. Therefore, a simple radius change fails to adequately differentiate between the behaviors of radial and nonradial modes.
6.3.2. Differential rotation
We now leave aside the origin of a possible change in rotation rate and simply explore what change in rotation would be required to explain the frequency changes seen in the nonradial modes. In Fig. 4 we show the frequency changes of the different pulsation modes relative to their 1986 values (colored points). We note that we have plotted the negative of the frequency changes for the retrograde modes; after accounting for the sign of m, we see that the prograde and retrograde modes contribute morphologically similar frequency changes.
The solid curve in Fig. 4 shows the average of these frequency changes; it is the average change in rotation frequency implied by these shifts, assuming solid body rotation. This interpretation makes sense qualitatively, explaining the overall behavior of the prograde, retrograde, and radial modes. Unfortunately, there is a problem: since C_{nl} ≪ 1 in δ Sct models, Eq. (1) implies that Δf_{nlm} ≈ m Δf_{rot}, that is, for a given change in rotation rate the magnitude of the frequency changes of all the nonradial modes should be nearly the same. As can be seen in Fig. 4, the frequency changes of the modes differ from each other by factors of up to two or three.
Fig. 4 Frequency changes of the different pulsation modes relative to their 1986 values (symbols). As indicated in the legend, we have plotted the negative of the frequency changes for the retrograde modes; the solid curve shows the average of these frequency changes. After accounting for the sign of m, we see that the prograde and retrograde modes contribute morphologically similar frequency changes. 
To salvage this explanation we must assume three things: (1) the star is differentially rotating; (2) different modes sample different effective rotation rates; and (3) the frequencies of the modes are more affected by rotation than is usually seen in pmode pulsations. Given that f_{rot}/f_{pulse} ≳ 0.2 is not small, the modes in this star are not in the slow rotation limit. Since Eq. (1) was derived assuming slow rotation, it may not be applicable to this star. In particular, the eigenmodes we observe could be more sensitive to changes in Ω(r,θ) than is expected in a slowly rotating star.
This interpretation naturally explains the opposite behavior of the prograde and retrograde modes, so we find it the most compelling possibility. Even so, it may require a revision of our understanding of the effects of rotation on the pulsations of δ Sct stars in the moderate rotation regime.
As a final point, we mention that the rotation profile we hypothesize is not changing due to an external torque on the star. Rather, we imagine that there are processes, perhaps the pulsations themselves, that redistribute angular momentum within the star. Since the modes in general have different rotation kernels, this will produce different period changes for different modes.
7. Conclusions
We have examined over 800 nights of highprecision photometry of the evolved multiperiodic pulsator 4 CVn obtained from 1966 through 2012. Most of the data were obtained in adjacent observing seasons. This made it possible to derive very accurate period values for a number of the excited pulsation modes and to study their systematic changes from 1974 to 2012.
Most pulsation modes of 4 CVn show systematic significant period and amplitude changes on a timescale of decades. Around 1986, a general reversal of the directions of both the positive and negative period changes was observed in the wellstudied modes. Furthermore, the period changes between the different modes are strongly correlated, although they differ in size and sign. An important discovery is that for the modes with known (ℓ,m) values, there exists a correlation between the direction of the period changes and the identified azimuthal order, m. The associated amplitude changes generally have similar timescales of years or decades, but show little systematic or correlated behavior from mode to mode.
To explain these observations, we have examined several potential physical explanations. Binary system effects, nonlinear mode interaction and changes in rotation due to a radius change were found to be unsatisfactory. However, a natural explanation for the opposite behavior of the prograde and retrograde modes is that their period changes are driven by a changing rotation profile. The changes in the rotation profile could in turn be driven by processes, perhaps the pulsations themselves, that redistribute angular momentum within the star. Our data show that the magnitudes of the period changes of the modes are different, which are not expected for pmode pulsations in the slow rotation limit. Given that f_{rot}/f_{pulse} ≳ 0.2 is not small, we see the variation in magnitude of period change as evidence that the pulsations are strongly modified by rotation.
Acknowledgments
M.B. is grateful to the Austrian Fonds zur Förderung der wissenschaftlichen Forschung for supporting this investigation through project P 21830N16. M.H.M. is grateful for the support of the NSF under grant AST1312983. A.A.P. acknowledges partial financial support from the Polish NCN grant 2011/01/B/ST9/05448.
References
 Blake, C., Fox, D. W., Park, H. S., & Williams, G. G. 2003, A&A, 399, 365 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Breger, M. 2000, MNRAS, 313, 129 [CrossRef] [Google Scholar]
 Breger, M. 2005, ASP Conf. Ser., 335, 85 [NASA ADS] [Google Scholar]
 Breger, M. 2016, A&A, 592, A97 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Breger, M., & Hiesberger, F. 1999, A&AS, 135, 547 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Breger, M., & Pamyatnykh, A. A. 1998, A&A, 332, 958 [NASA ADS] [Google Scholar]
 Breger, M., McNamara, B. J., Kerschbaum, F., et al. 1990, A&A, 231, 56 [NASA ADS] [Google Scholar]
 Breger, M., Handler, G., Garrido, R., et al. 1997, A&A, 324, 566 [NASA ADS] [Google Scholar]
 Breger, M., Handler, G., Garrido, R., et al. 1999, A&A, 349, 225 [NASA ADS] [Google Scholar]
 Breger, M., Davis, K. A., & Dukes, R. J. 2008, Comm. Asteroseismol., 153, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Breger, M., Lenz, P., & Pamyatnykh, A. A. 2009, MNRAS, 396, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Castanheira, B., Breger, M., Beck, P. G., et al. 2008, Comm. Asteroseismol. 157, 124 [Google Scholar]
 DaszyńskaDaskiewicz, J., Dziembowski, W. A., & Pamyatnykh, A. A. 2003, A&A, 407, 999 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fitch, W. S. 1980, Lect. Notes Phys., 125, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Handler, G., Arentoft, T., Shobbrook, R. R., et al. 2000, MNRAS, 318, 511 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, D. H. P., & Haslam, C. M. 1966, The Observatory, 86, 34 [NASA ADS] [Google Scholar]
 Lenz, P., & Breger, M. 2005, CoAst, 146, 53 [Google Scholar]
 Lenz, P., Pamyatnykh, A. A., & Breger, M. 2010, in Variable Stars, the Galactic halo and Galaxy Formation, Proceedings of an international conference held in Zvenigorod, Russia, 12–16 October 2009, published by Sternberg Astronomical Institute of Moscow University, Russia, 83 [Google Scholar]
 Loumos, G. 1980, Ph.D. dissertation, unpublished (Univ. of Texas Austin) [Google Scholar]
 Moskalik, P. 1985, Acta Astron., 35, 229 [NASA ADS] [Google Scholar]
 Murphy, S. J., Grigahcène, A., Niemczura, E., Kurtz, D. W., & Uytterhoeven, K. 2012, MNRAS, 427, 1418 [NASA ADS] [CrossRef] [Google Scholar]
 Murphy, S. J., Shibahashi, H., & Bedding, T. R. 2016, MNRAS, 461, 4215 [Google Scholar]
 Schmid, V., Themessl, N., Breger, M., et al. 2014, A&A, 570, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shibahashi, H., & Kurtz, D. W. 2012, MNRAS, 422, 738 [NASA ADS] [CrossRef] [Google Scholar]
 Shibahashi, H., Kurtz, D. W., & Murphy, S. J. 2015, MNRAS, 450, 3999 [NASA ADS] [CrossRef] [Google Scholar]
 Stankov, A., Handler, G., Mkrtichian, D. E., et al. 2000, J. Astron. Data, 6, 5 [Google Scholar]
 Turner, D. G., Savoy, J., Derrah, J., AbdelSabour AbdelLatif, M., & Berdnikov, L. 2005, PASP, 117, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Warman, J., Pena, J. H., & ArrelanoFerro, A. 1979, AJ, 84, 109 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Amplitude ratios and phase differences in the Strömgren v and y passbands for different pulsation modes. The error bars are determined by comparing results from eight different years. In A/F stars the observed phase differences are strong indicators of the spherical degree. The results from a theoretical model (see text) covering the frequency range from 4.5 to 9.5 cycle d^{1} are also shown. 

In the text 
Fig. 2 Frequency variations of the six dominant modes from 1974 to 2012. The errors bars shown are usually the size of the symbols. The diagram shows the systematic difference between prograde and retrograde modes. Furthermore, around ~1986 the sign of the frequency changes reversed systematically. 

In the text 
Fig. 3 Confirmation of the binary orbit from a comparison of the 2005–2012 photometric data without and with the orbital lighttime corrections. Top left panel: Fourier amplitudes of the observed data near the f_{1} frequency. Bottom left panel: Fourier amplitudes of the residuals after prewhitening the best 64frequency solutions. Note the changed amplitude scale. The uncorrected data (dotted curve) shows the side peaks from the 124 (~1/ 0.008) d orbit, while the data corrected for the spectroscopically detected binary orbit eliminates these peaks. The difference of the f_{1} side peaks from the lower left panel is demonstrated in the top right panel. A similar analysis for f_{2} (bottom right panel) confirms the detection. 

In the text 
Fig. 4 Frequency changes of the different pulsation modes relative to their 1986 values (symbols). As indicated in the legend, we have plotted the negative of the frequency changes for the retrograde modes; the solid curve shows the average of these frequency changes. After accounting for the sign of m, we see that the prograde and retrograde modes contribute morphologically similar frequency changes. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.