## Revisiting CoRoT RRLyrae stars: detection of period doubling and temporal variation of additional frequencies ? ?? ??? ' '

R. Szab´ o 1 , J. M. Benk˝ o 1 , M. Papar´ o 1 , E. Chapellier 2 , E. Poretti 3 , A. Baglin 4 , W. W. Weiss 5 , K. Kolenberg 6 7 ; , E. Guggenberger 8 9 ; , and J.-F. Le Borgne 10 11 ;

- 1 Konkoly Observatory, MTA CSFK, Konkoly Thege Mikl´ os ´ ut 15-17. H-1121 Budapest, Hungary e-mail: rszabo@konkoly.hu
- 2 Laboratoire Lagrange, Universit´ e Nice Sophia-Antipolis, UMR 7293, Observatoire de la Cˆte d'Azur 06300, Nice, France o
- 3 INAF - Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807 Merate (LC), Italy
- 4 LESIA, Universit´ e Pierre et Marie Curie, Universit´ e Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
- 5 Institute of Astronomy, University of Vienna, T¨rkenschanzstrasse 17, 1180 Vienna, Austria u
- 6 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA
- 7 Instituut voor Sterrenkunde, K.U. Leuven, Celestijnenlaan 200D, B-3001 Heverlee, Belgium
- 8 Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, 37077 G¨ ottingen, Germany
- 9 Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
- 10 Universit´ e de Toulouse, UPS-OMP, IRAP, Toulouse, France
- 11 CNRS, IRAP, 14, avenue Edouard Belin, F-31400 Toulouse, France

received; accepted

## ABSTRACT

Context. High-precision, space-based photometric missions like CoRoT and Kepler revealed new and surprising phenomena in classical variable stars. Such discoveries were the period doubling in RR Lyrae stars and the frequent occurrence of additional periodicities some of which can be explained by radial overtone modes, but others are discordant with the radial eigenfrequency spectrum.

Aims. We search for signs of period doubling in CoRoT RR Lyrae stars. The occurrence of this dynamical e GLYPH&lt;11&gt; ect in modulated RRLyrae stars might help us to gain more information about the mysterious Blazhko e GLYPH&lt;11&gt; ect. The temporal variability of the additional frequencies in representatives of all subtypes of RR Lyrae stars is also investigated.

Methods. We pre-process CoRoT light curves by applying trend and jump correction and outlier removal. Standard Fourier technique is used to analyze the frequency content of our targets and follow the time dependent phenomena.

Results. The most comprehensive collection of CoRoT RR Lyrae stars, including new discoveries is presented and analyzed. We found alternating maxima and in some cases half-integer frequencies in four CoRoT Blazhko RR Lyrae stars, as clear signs of the presence of period doubling. This reinforces that period doubling is an important ingredient to understand the Blazhko e GLYPH&lt;11&gt; ect - a premise we derived previously from the Kepler RRLyrae sample. As expected, period doubling is detectable only for short time intervals in most modulated RRab stars. Our results show that the temporal variability of the additional frequencies in all RR Lyrae sub-types is ubiquitous. The ephemeral nature and the highly variable amplitude of these variations suggest a complex underlying dynamics of and an intricate interplay between radial and possibly nonradial modes in RR Lyrae stars. The omnipresence of additional modes in all types of RR Lyrae - except in non-modulated RRab stars - implies that asteroseismology of these objects should be feasible in the near future.

Key words. Stars: variables: RR Lyrae - stars: oscillations - stars: interiors - techniques: photometric

## 1. Introduction

MOST (Walker et al. 2003), CoRoT (Baglin et al. 2006) and Kepler (Borucki et al. 2010).

The advent of space photometry has opened up new vistas in investigating stellar pulsations and oscillations. Pulsating variable stars located in the classical instability strip are no exception. RR Lyrae stars in particular have benefited from the continuous and ultra-precise space photometric data delivered by

? The CoRoT space mission was developed and operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany and Spain.

?? Follow-up observations were obtained at Piszk´ es-tet˝, o the Mountain Station of Konkoly Observatory.

??? Table 3 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http: // cdsweb.ustrasbg.fr cgi-bin qcat?J A / / / + A /

One of the surprising findings was the period doubling phenomenon (hereafter PD) in the ultra-precise Kepler RR Lyrae light curves (Kolenberg et al. 2010; Szab´ o et al. 2010), which manifests as the alternating height of maximum brightness and alternating light curve shape from cycle to cycle in the photometric light curve and half-integer frequencies (HIFs) between the dominant pulsation mode and its harmonics in the frequency domain. So far only the Kepler stars were investigated to describe the phenomenon, such a detailed analysis is still missing for CoRoT RRLyrae stars. This prompted us to reexamine the published CoRoT RR Lyrae light curves and extend the study by adding new ones, as well.

PD is found only in case of modulated (Blazhko) RR Lyrae stars, though a thorough search was done on non-modulated

Table 1. Basic parameters of the CoRoT RR Lyrae stars. The columns contain the CoRoT ID number, the coordinates, V magnitude, the CoRoT run, RR Lyrae sub-type, the amplitude of the dominant mode(s), the Blazhko period in case of modulated stars, the pulsation period and references.

| CoRoT ID    | R.A. (2000)   | Dec. (2000)   | V [mag]   | Run   | Type    | f 0 ; 1 ampl. [mag]   | P Bl [days]   |   puls. period [days] | Ref.      |
|-------------|---------------|---------------|-----------|-------|---------|-----------------------|---------------|-----------------------|-----------|
| 0101370 131 | 19 28 14.40   | + 0 06 02.27  | 15.60     | LRc01 | RRab    | 0 : 257               | -             |              0.619332 | (1)       |
| 0101315 488 | 19 27 47.40   | + 0 58 35.62  | 16.15     | LRc01 | RRab    | 0 : 029 b             | -             |              0.485303 | this work |
| 0103800 818 | 18 31 23.80   | + 9 10 10.45  | 13.70     | LRc04 | RRab    | 0 : 345 c             | -             |              0.465935 | this work |
| 0104315 804 | 18 34 28.92   | + 8 57 00.43  | 15.94     | LRc04 | RRab    | 0 : 130               | -             |              0.721822 | this work |
| 0100689 962 | 19 24 00.10   | + 1 41 48.70  | 14.96     | LRc01 | RRab Bl | 0.201                 | 26.88         |              0.355997 | (2)       |
| 0101128 793 | 19 26 37.32   | + 1 13 34.90  | 15.93     | LRc01 | RRab Bl | 0.235                 | 17.86         |              0.47193  | (3)       |
| 0100881 648 | 19 25 05.43   | + 1 39 23.83  | 16.16     | LRc01 | RRab Bl | 0 : 033 b             | 59.77         |              0.607186 | this work |
| 0101503 544 | 19 29 10.12   | + 0 43 47.14  | 14.79     | LRc01 | RRab Bl | 0 : 009 b ; c         | 25.60         |              0.605087 | this work |
| 0105288 363 | 18 39 30.86   | + 7 26 53.95  | 15.32     | LRc02 | RRab Bl | 0 : 179               | 35.6          |              0.567441 | (4,5,6)   |
| 0103922 434 | 18 32 08.55   | + 8 32 40.78  | 15.84     | LRc04 | RRab Bl | 0 : 280               | 54.5          |              0.541383 | (7)       |
| 0105036 241 | 18 38 09.60   | + 7 43 56.68  | 15.58     | LRc02 | RRc     | 0 : 196               | -             |              0.372921 | this work |
| 0105735 652 | 18 42 10.13   | + 6 33 05.15  | 15.01     | LRc02 | RRc     | 0 : 204               | -             |              0.27916  | this work |
| 0101368 812 | 19 28 13.61   | + 0 40 42.46  | 15.86     | LRc01 | RRd f 0 | 0 : 053               | -             |              0.488041 | (8)       |
|             |               |               |           |       | RRd f 1 | 0 : 143               | -             |              0.363602 |           |

Notes. The superscript 'b' denotes blended pulsators. Superscript 'c' denotes a CoRoT target with colors. V magnitude values were taken from the ExoDat catalog. No good candidates classified as RR Lyrae stars were found in the third runs. References: (1) : Papar´ et al. (2009), o (2) : Chadid et al. (2010), (3) : Poretti et al. (2010), (4) : Guggenberger et al. (2011), (5) : Chadid et al. (2011), (6) : Guggenberger et al. (2012), (7) : Poretti et al., in prep., (8) : Chadid (2012)

Kepler RRab stars, too (Szab´ o et al. 2010; Nemec et al. 2011). Therefore, it is plausible to assume that either PD as a nonlinear dynamical phenomenon plays an important role in causing the Blazhko-modulation (see Smolec et al. 2012 for a hydrodynamic example of PD causing modulation in BL Herculis stars) or at least it shows up as a frequent companion e GLYPH&lt;11&gt; ect of the modulation. In either case, it is important to find well-documented cases, establish the frequency of the PD occurrence, investigate its temporal behavior and other characteristics. In addition, the physical explanation of PD itself was unambiguously traced back to a 9:2 resonance between the fundamental pulsational mode and a high-order radial (strange) overtone (Koll´ ath et al. 2011) which opens a new way to study the dynamics of these high-amplitude variable stars belonging to the horizontal branch.

find out how frequent this phenomenon is, and further strengthen its role in and connection to the mysterious Blazhko-e GLYPH&lt;11&gt; ect. Our second aim is to investigate the recently found new periodicities in all types of RR Lyrae stars. This kind of research is only possible with the photometric precision of CoRoT for these stars. We are especially interested in the temporal stability of these additional frequencies.

In recent years a new picture has started to emerge regarding the frequency spectrum of RR Lyrae stars. Thanks to dedicated telescopes (Jurcsik et al. 2009a) and space-based photometry, additional frequencies were found in many stars that were observed with at least millimagnitude precision: AQ Leo (RRd) (Gruberbauer et al. 2007), MW Lyr (Jurcsik et al. 2008), CoRoT (Poretti et al. 2010; Chadid et al. 2010) and Kepler stars (Benk˝ o et al. 2010).

These additional periodicities do not fit the series of radial pulsational modes (nor the occasional half-integer series due to PD) and present a low-amplitude variability in each case. Presently, the best explanation for their presence is the excitation of nonradial modes. If this turns out to be true, we may have another handle on the interior of horizontal branch stars to understand better the stability, pulsation and evolution of these standard candles. In some cases we unexpectedly found frequencies with low-amplitudes at or near the radial overtones (e.g. Benk˝ o et al. 2010) that might correspond to radial modes excited by resonances or nonradial modes with frequencies being in 1:1 resonance with the corresponding radial mode (Dziembowski &amp; Mizerski 2004; Van Hoolst et al. 1998). The era of exploiting the power of nonlinear seismology using radial modes in RR Lyrae stars is imminent (Moln´ ar et al. 2012).

In this work we embarked on investigating all the known CoRoT Blazhko RRLyrae stars up to LRc04 in order to detect period-doubling. The detection or non-detection would help to

The structure of this paper is the following: in Sec. 2 we introduce the original space-borne and ground-based follow-up data we use in this analysis. In Sec. 3 we discuss the methods to find period doubling and to detect the variations of low amplitude periodicities. In Sec. 4 we present our results, then in Sec. 5 we discuss them and draw our conclusions followed by a short summary in Sec. 6. Hitherto unpublished frequency tables and other complementary results can be found in the Appendices.

Throughout this work we use the following notations: by f 0 , f 1 , f 2 , ... we denote the frequency of the radial pulsation modes (fundamental, first overtone, second overtone, etc.), while f 0 , f 00 stand for additional, independent frequencies. We refer to the CoRoT stars with the last three digits of their ID number for short notation, which provides a unique identification, e.g. , 962 for CoRoT 0100689 962 .

## 2. Observations

## 2.1. CoRoT observations

CoRoT has conducted long-duration, continuous, very-highprecision relative photometry on a few very bright stars, and a large number of faint ones (Baglin et al. 2006). As an attempt to review the CoRoT RR Lyrae stars we investigate 13 stars listed in Table 1. Table 1 gives the CoRoT ID, the coordinates, the V magnitude, the CoRoT run in which the object was observed and discovered, its pulsation and modulation parameters (if applicable), as well as references if the star was already analyzed.

The brightness of the targets that were observed in the exofield of CoRoT are measured through a prism which creates images that contain red, green and blue fluxes. The color fluxes are not available for all stars, in some cases only the sum of

Table 2. Observations of the CoRoT LRc01 RR Lyrae stars at the Konkoly Observatory that were used in this paper. The CoRoT ID, the night of the observations, total number of the scientific frames (Nf), the ID of the primary comparison stars from the USNO-A2.0 are given.

|   CoRoT ID | DATE                                                                                            | N f                | Comp.         |
|------------|-------------------------------------------------------------------------------------------------|--------------------|---------------|
|  100881648 | 2008-07-28 / 29 2008-07-29 / 30                                                                 | 55 30              | 0900-14969350 |
|  101503544 | 2008-06-20 / 21 2008-06-21 / 22                                                                 | 49 68              | 0900-15291694 |
|  100689962 | 2008-06-09 / 10 2008-06-19 / 20 2008-07-08 / 09                                                 | 57 58 38           | 0900-14903871 |
|  101370131 | 2008-07-11 / 12 2008-06-22 / 23 2008-06-25 / 26 2008-07-10 / 11 2008-07-11 / 12 2008-07-12 / 13 | 57 110 33 85 36 38 | 0900-15209129 |
|  101128793 | 2008-07-27 / 28 2009-06-29 / 30                                                                 | 63 47              | 0900-15077950 |

them (white light) is retained. Though these fluxes are not related to any photometric systems, they are still useful for example to compare the amplitudes of variability in some of our cases. The interested reader can find an example for the use of CoRoT color fluxes in Papar´ et al. (2011). o

The CoRoT RR Lyrae group analyzed those stars that were classified as such by the CoRoT Variability Classifier (CVC) (Debosscher et al. 2009) up to the fourth pointing (LRc04). We note that one of these stars, CoRoT 1027817 750 , observed during LRa01, originally classified as an RR Lyrae star with a 67% probability is found to be non-RR Lyrae by Papar´ et al. (2011). o Wetook those targets that were classified as RRab pulsators with non-zero probability and found eight bona fide RRab stars out of 14. In addition to that sample we found CoRoT RR Lyrae star 648 early in the mission. Furthermore we added one more object to the classified RR Lyrae stars which was found by A GLYPH&lt;11&gt; er et al. (2012) while investigating stellar rotation periods in CoRoT light curves in runs LRc01 and LRa01. Upon inspecting the published light curves manually, we confirmed that 488 is a new RRab variable. Light curve characteristics, frequencies, epoch and more details on this target can be found in Appendix A.

RRc classification is more problematic, as eclipsing binaries ( e.g. , WUMa stars) can exhibit similar light curve shapes, thus the contamination of a sample originating from automated classification is much higher. Therefore we chose the &gt; 80% probability level assigned by CVC and found only two genuine RRc stars out of 17 candidates. We note that decreasing the probability limit to 50% or 10% did not help to find more RRc stars, but increased the number of eclipsing binaries and other types enormously.

Interestingly, we have RR Lyrae stars only in the direction of the Galactic Center. This obviously might be a selection e GLYPH&lt;11&gt; ect due to the fact that that we observe mostly discpopulation stars when the satellite points toward the anticenter direction. The length of observations varies from 158 days (LRc01) and 150 days (LRc02) to 88 days (LRc04), in all cases the duty cycle being over 90%. Time is given in Heliocentric Julian Date throughout the paper, however when showing light curves we use CoRoT Julian Date (CJD). The relation between the Heliocentric Julian Date and the CJD is

$$\text{from the } \ H J D = C J D + 2 4 5 1 5 4 5. 0$$

We used the N2 level calibrated CoRoT light curves. We applied trend and jump filtering and outlier removal as described in detail in Chadid et al. (2010). As we usually deal with relatively faint targets, only two RR Lyrae in our sample have CoRoT colors ( 544 and 818 ), but we use only the integrated (white) light to increase the signal-to-noise ratio. The nominal sampling is 512 sec for all the stars, except the above mentioned two targets, where the initial part of the light curves were observed with the 512 sec mode while the rest of the light curve was observed with the much denser short cadence mode (32 sec).

## 2.2. Ground-based follow-up observations

Ground-based multi-color observations on some of the CoRoT RRLyrae stars discussed in this paper were collected in 2008 and 2009 with the 1-m Ritchey-Chr´ etien-Coud´ e (RCC) telescope mounted at Piszk´ es-tet˝ Mountain Station of the Konkoly o Observatory. These observations were especially helpful in determining the true pulsational nature in the cases where our targets are heavily blended in the CoRoT apertures.

AVersarray 1300B camera with an UV-enhancement coating constructed by Princeton Instruments was used. This device contains a back illuminated EEV CCD36-40 1340 GLYPH&lt;2&gt; 1300 chip that corresponds to a 6 6 : 0 GLYPH&lt;2&gt; 6 8 field of view (FOV) with 0.303 : 0 00 / pixel scale. Standard Johnson BV and Kron-Cousins R C filters were used in the observations. The typical exposure times were 300, 180, 100 s for bands B V , and R C, respectively. Each night dome flats, bias and dark frames were taken as main calibration images.

We used the iraf/ccdred 1 package for the standard reduction procedures: bias, dark and flat field correction. The brightness of stars was determined by using the aperture photometry task daophot/phot of iraf . We carried out di GLYPH&lt;11&gt; erential photometry with carefully selected comparison stars. These were found to be constant and fit the target stars both in their brightness and colours. The telescope constants were obtained using the standard stars of the open cluster M67 (Chevalier &amp; Ilovaisky 1991). The typical errors of the individual observations were between 0.01-0.02 mag in BVR C. The log of the observations is presented in Table 2, while the ground-based BVR C light curves are available in electronic form via the website of the journal. A sample of this file is shown in Table 3 for guidance.

## 3. Methodology

We applied standard packages such as M u F A r n (Koll´ ath 1990) and P eriod 04 (Lenz &amp; Breger 2005) to perform Fourier-analysis.

## 3.1. Detection of period doubling

There is no obvious, strongly alternating pattern in the CoRoT RRLyrae maxima that could have triggered a thorough analysis

1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation.

amplitude [mmag]

amplitude [mmag]

freq. [c/d]

3.2

4.0

3.0

2.0

1.0

0.0

3.0

2.0

1.0

0.0

3040

3.2

3060

3.2

3080

3.2

2.8

3100

2.8

3120

2.8

3140

2.8

3.2

3160

CoRoT JD

Fig. 1. Test results for the detection of the temporal variation of low amplitude frequencies. Left panel : The amplitude of the constant, 3 mmag amplitude frequency injected at f test = 3 0 : d GLYPH&lt;0&gt; 1 into the light curve of 241 , a CoRoT RRc star (see Sec. 4.3). Right panel : A more complex case: we added a constant (3 mmag) amplitude frequency to the light curve of the strongly modulated Blazhko RR Lyrae, 962 . We chose f test = 3 4 d : GLYPH&lt;0&gt; 1 . Horizontal lines show the amplitude of the injected signal, points with error bars come from our method. The upper row shows the shape and the vicinity of the recovered frequency in the applied eight bins.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-1d34fc79-baea-4abd-8bad-a841e992470c.jpg)

Table 3. Ground-based multi-color observations of CoRoT RRLyrae stars taken with the 1m RCC telescope of the Piszk´ estet˝ o Mountain Station of the Konkoly Observatory. This table is published in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content.

|   CoRoT ID |         HJD |   magnitude | filter   |
|------------|-------------|-------------|----------|
|  100881648 | 2.45468e+06 |      16.903 | B        |
|  100881648 | 2.45468e+06 |      16.919 | B        |
|  100881648 | 2.45468e+06 |      16.913 | B        |
|  100881648 | 2.45468e+06 |      16.933 | B        |
|  100881648 | 2.45468e+06 |      16.913 | B        |

earlier, though in some cases a detailed inspection do show period doubling persisting through a few pulsational cycles as we show later in this work.

CoRoT light curves are shorter than the Kepler ones and are of lower precision at the same apparent brightness. The typical error of the individual data points for a 16 th magnitude RR Lyrae is 0.006 - 0.012 mag for the CoRoT observations, and 0.0008 - 0.003 mag for an RR Lyrae of the same brightness observed by Kepler . In addition, since period doubling is usually only temporarily noticeable in the light curves (but see Le Borgne et al. 2014), we do not expect to see PD signs in both the light curve and the frequency spectrum in each cases. The strength of PD is not correlated with the Blazhko phase, as we saw earlier on the example of RR Lyrae itself and other modulated RRab stars. Therefore here we need to relax our strict criteria applied in Szab´ o et al. (2010), where we required the presence of PD signs in both the time and the frequency domain. Consequently, in this work we investigate both the Fourier-spectrum between the dominant mode and its harmonics aiming at finding HIFs at exactly between two consecutive harmonics and the light curves themselves looking for suspicious alternating pattern. We report tentative detection in those cases where PD is present in either or both domains.

For the detection of PD in the Fourier spectrum we set a conservative limit, i.e. , 3 GLYPH&lt;27&gt; detection for the HIFs. The rms for half-integer frequencies ((2 k + 1) = 2 GLYPH&lt;1&gt; f 0 ) was computed in the intervals between consecutive harmonics of the dominant pulsational mode [ k GLYPH&lt;1&gt; f 0 ; ( k + 1) GLYPH&lt;1&gt; f 0 ] after pre-whitening with significant additional frequencies and the Blazhko side peaks. The [ f 0 ;2 f 0 ] frequency interval is the most interesting one, because all the known stars undergoing the period doubling bifurcation found in the Kepler sample show the highest half-integer frequency amplitude in this particular interval (Szab´ et al. 2010). o

In order to investigate the alternating nature of the pulsational cycles, the photometric maxima were fitted with 7th or 9th order polynomials (Chadid et al. 2010) to get rid of the detrimental e GLYPH&lt;11&gt; ects of missing photometric points or (in the case of blended Blazhko stars) large scatter. For more details of this process we refer to Chadid et al. (2010).

## 3.2. Frequency analysis and time-dependent frequencies

To get the time-dependence of the properties of a frequency peak we cut the photometric data into equal-length chunks, and the respective amplitudes and frequencies and their uncertainties were obtained from these shorter data sections with P eriod 04. We tried to cut the light curve into two-four-eight etc. pieces, and stopped where the uncertainties (primarily due to the shortening of the data sections) precludes the derivation of any meaningful information. In most of the cases, for uniformity, eight bins were used. This is a compromise between time resolution of the variabilities and the frequency resolution permitted by the lengths of the data chunks, which could not be deteriorated arbitrarily, because often close-by frequency peaks occur. We performed simulations to make sure that any arbitrary time shift ( i.e. , the exact starting epoch of the subsections) does not have significant e GLYPH&lt;11&gt; ect on the amplitudes of the studied frequencies.

In all cases the main pulsation mode(s), their harmonics, the Blazhko frequency, its harmonics and the modulation side peaks (where applicable) were subtracted before we embarked on computing the time variable amplitude of the additional peaks. Error bars were derived by a Monte Carlo method available in P eriod 04.

2.8

2.8

2.8

3.2

3.2

3.2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-807f49a6-da65-4bf3-b459-3fcff95d560c.jpg)

0.0

phase

0.0

phase phase

Fig. 2. Jump- and trend filtered phased light curves of CoRoT non-Blazhko RRab stars. Left panel : The newly found RR Lyrae, 488 . Note the small amplitude and the relatively large scatter, indicating a blended CoRoT target. Middle panel: The new, nonmodulated RR Lyrae, 818 . Having over 177,000 data points, we plot only every 10 th data points for visualization purposes. Right panel: Light curve of the new, non-modulated RR Lyrae, 804 .

In order to check whether the successive pre-whitening, the strong Blazhko modulation and other time-dependent features have an appreciable e GLYPH&lt;11&gt; ect on our method used to detect temporal variation of low amplitude frequency peaks, we designed a series of tests. First we added a 3.0 mmag, constant amplitude sinusoidal signal with a frequency of 3 0 d : GLYPH&lt;0&gt; 1 to the light curve of 241 , an RRc star (see Sec. 4.3). The vicinity of 3 0 d : GLYPH&lt;0&gt; 1 in the Fourier spectrum is relatively clean, i.e. , free of contaminating frequencies. This light curve was subjected to the same procedure as all the CoRoT targets, namely the light curve was cut into eight bins, and the main pulsation frequency and its harmonics were subtracted in each bin. As we see in the left panel of Fig. 1, the injected constant amplitude signal was preserved, the error introduced by the process is a few percent, and does not exceed 10% in any case. This means that in a sparse frequency spectrum we can confidently detect the time-dependence of lowamplitude signals without introducing amplitude variability with our methods.

The right panel of Fig. 1 shows a more complex case. This time we chose 962 , a strongly modulated Blazhko star. The frequency spectrum of this star contains more frequencies (Chadid et al. 2010) making the analysis more complicated. Now we squeezed an f test = 3 4 d : GLYPH&lt;0&gt; 1 signal with 3-mmag constant amplitude in between other frequencies. As suspected, in this case we got a higher scatter around the constant amplitude. Ten percent deviation from the nominal value is not unusual, and even larger excursions occur in some cases. The computed error bars more or less capture this scatter, the deviation rarely exceeds the assigned sigma. Evidently, in case of a dense frequency spectrum with evidence of strong temporal variations (such as a Blazhkomodulation and period-doubling) one has to exercise extra care when interpreting amplitude variability of low-amplitude signals. It is reassuring, however, that the structure of the frequency peaks does not deform due to the procedure, as long as there are no close immediate frequency peaks in the vicinity of our test signal.

## 4. Results

In the next subsections we discuss our targets one by one focusing on those stars and features that have not been discussed earlier. In case of previously published results or those in preparation we refer to the corresponding paper(s). However, for the sake of completeness we decided to list all CoRoT objects reliably classified as RR Lyrae stars up to LRc04. First we describe non-modulated RRab stars followed by Blazhko stars, then RRc

Fig. 3. Left: 2.6' x 2.6' DSS image showing the neighborhood of 488 . Right: A larger surroundings of the star 488 showing the applied aperture as well. Several faint stars are shown within the aperture. North is up and East is to the left in both images.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-a6f64aeb-f5c8-4944-a9fe-fe25e8ee640d.jpg)

stars are listed, finally the only double-mode (RRd) star is discussed.

## 4.1. Non-Blazhko RRab stars

CoRoT0101370 131 : This mono-periodic RRab star was analyzed in detail in Papar´ et al. (2009). In the course of this work o we checked again that no additional frequencies, including signs of HIFs can be found in the frequency spectrum. We give an upper limit of 0.30 mmag for the latter (Table 4). We also checked whether the maxima show any periodicity or pattern, as a sign of a so far hidden, low-amplitude Blazhko-e GLYPH&lt;11&gt; ect, but failed to find any significant periodicities.

CoRoT0101315 488 : A er et al. GLYPH&lt;11&gt; (2012) found 169 pulsating stars while investigating stellar rotation periods and ages in CoRoT light curves of spotted stars in runs LRc01 and LRa01. Upon inspecting the published light curves manually we encountered a new RR Lyrae like object: 488 , besides finding other, already known RR Lyrae variables.

As this is a new finding, we discuss the properties of this object in detail. The coordinates of the new variable are RA = 19 h 27 m 47 s : 40, Dec = 0 GLYPH&lt;14&gt; 58' 35 652. For the brightness of the : ' main target the ExoDat database gives the following values: B = 17.163, V = 16.152, R = 15.740, I = 15.072. Based on the contamination index (0.222), blending can be considered to be severe.

mag

∆

-0.8

-0.6

-0.4

-0.2

0.0

0.2

0.4

-0.4

-0.2

0.2

0.4

mag

∆

-0.20

-0.15

-0.10

-0.05

0.00

0.05

0.10

0.15

0.20

-0.4

-0.2

0.2

0.4

The light curve shape is of a typical RRab, but the total amplitude is rather small (0 : m 1). Hence, we suspect that this variable star is indeed blended. The relatively high scatter also corroborates this argument (see Fig. 2, left panel). The pulsation period is constant during the CoRoT run.

Fig. 3 shows the DSS and the EXODAT images of our target. EXODAT lists six close-by stars of magnitude 17.5 - 20.5, within the CoRoT aperture. It is conceivable that one of these might be the RR Lyrae, but based on the available data we cannot tag the exact source of variation. The pulsation parameters, frequency table, epoch and discussion of the frequency spectrum can be found in Appendix A. We note here only that no PD was found in the spectrum with an upper limit of 0.10 mmag for any half-integer frequencies.

CoRoT0103800 818 : This star was observed by CoRoT in LRc04, and it is described here for the first time. It started to be observed with nominal cadence (512 s) but soon after it was switched to short cadence (32 s) observations in color mode. To enhance the signal-to-noise here we use the integrated (white) light curve only. With a period of 0.4659348 days it is a nonmodulated RRab star. We did not detect period doubling with an the upper limit of 0.10 mmag set for the half-integer frequencies. The phase-folded light curve is shown in the middle panel of Fig. 2.

Due to the early switch to high cadence observations, fiftysix harmonics has been found in the frequency spectrum. So far this is the largest number we are aware of. The amplitude distribution is very similar to 131 (Papar´ o et al. 2009), i.e. , there is a local minimum in the amplitudes of the harmonics around the 15 th harmonics (32 d GLYPH&lt;0&gt; 1 ). The harmonics in the higher frequency range come with a di GLYPH&lt;11&gt; erent slope after this trough in the distribution. While we can not confirm this behavior with the other two non-modulated CoRoT RRab stars ( 488 and 804 ) because of the lower cadence observations, it would be interesting to see the amplitude distribution in the case of other RRab stars observed from space with high cadence.

Apart from that, only the CoRoT orbital frequencies and some residuals around the main frequency peaks are visible in the frequency spectrum. We could identify neither regular frequency side peaks that would signal the presence of the Blazhko e GLYPH&lt;11&gt; ect, nor half-integer or any other frequencies. The frequency table is given in Appendix B.

CoRoT0104315 804 : This star was also observed in LRc04, and it has not been published, yet. It was monitored with nominal cadence (512s) and only in white light. It is a non-Blazhko star with a period of 0.7218221 days. We could not detect period doubling with an upper limit of 0.30 mmag. The frequency spectrum is also devoid of statistically significant additional periodicities. The frequency table is given in Appendix C, while the phase-folded light curve is plotted in the right panel of Fig. 2.

## 4.2. RRab stars showing the Blazhko effect

CoRoT0100689 962 is identical with V1127Aql. It has a strongly modulated light curve, with a high pulsational amplitude and a pronounced phase modulation as is easily seen in the left panel of Fig. 4. The modulation is not sinusoidal (Benk˝ o et al. 2011), and is asymmetric in the sense that the maximum of the maxima and the minimum of the minima shows considerable shift. A detailed study has been published on this object by Chadid et al. (2010).

Fig. 7. Same as the left panel of Fig. 6, but for CoRoT Blazhko RRab 793 . The dominant pulsation frequency ( f 0 ), its harmonics, the Blazhko-frequency ( f B ) and 2 f B (Poretti et al. 2010) were pre-whitened. The HIFs are significant up to 17 = 2 f 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-3cb5e66e-85ec-440e-b053-d01f1ebb6330.jpg)

We found alternating cycles in the light curve of 962 (see the right panel of Fig. 5.) Half-integer frequencies between the fundamental mode and its harmonics up to 11 = 2 f 0 are also found, demonstrated in the left panel of Fig. 5. A forest of peaks can be found around the HIFs. This may be a consequence of the excitation of additional (nonradial) modes, or due to the interplay of the time-variable PD and the modulation itself (see a simulation suggesting this scenario in Szab´ et al. (2010)), or both. The amo plitude distribution of the HIFs is very similar to other cases (see Szab´ o et al. 2010): the highest peak is at 3 = 2 f 0 , then there is a monotone decrease with the order, however at 9 = 2 f 0 there is a local maximum (a standstill in the case of 962 ), which is the consequence of the 9:2 resonance with the 9th overtone (Koll´ ath et al. 2011). We could identify HIFs up to 11 = 2 f 0 above the 3 GLYPH&lt;27&gt; significance level. The alternating cycles in the light curve and the presence of the half-integer frequencies in the frequency domain constitute an unambiguous detection of the period doubling phenomenon in 962 . It is evident that the PD is not always present with the same strength throughout the 142-day CoRoT observations. Instead, we see it in the light curve for short time intervals, e.g. , CJD[2750:2754] as shown in Fig. 5 or CJD[2794:2798] (not shown).

In the frequency spectrum of 962 there is another, independent frequency peak, f 0 = 4 03265d : GLYPH&lt;0&gt; 1 . This is quite close to the 3 2 / f 0 frequency value, and there is a dense region of frequencies around them (Fig. 7 of Chadid et al. 2010), which makes the detection of the temporal variation of these frequencies extremely challenging. We made a series of tests to take into account the e GLYPH&lt;11&gt; ects of the decreased frequency resolution due the splitting of the original data set, and we found variation in the amplitude in both the 3 2 / f 0 HIF and f 0 (upper left panel of Fig. 6), although this later conclusion is quite weak and may be a consequence of contamination from other close-by frequencies.

CoRoT0101128 793 is a CoRoT RRLyrae star with a low Blazhko modulation amplitude (the amplitude of the highest side peak due to the modulation is 4.6 mmag) and a Blazhko period of 18.66 d (see Fig. 4). This star shows one of the lowest modulation amplitudes detected so far (Jurcsik et al. 2006; Skarka 2014). Poretti et al. (2010) studied this object in detail. In the course of this work we find clear sign of period doubling in the frequency spectrum in the form of half-integer frequencies

Fig. 4. Phase-folded CoRoT Blazhko light curves. Left : folded by the Blazhko-period, right : folded by the pulsational period. Representative averaged light curves are shown in two extreme modulation phases in each case in the right panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-660ae705-5588-4e4d-8130-e4274b3e6762.jpg)

Fig. 5. Left panel: Low-frequency part of the Fourier-spectrum of 962 . The dominant pulsation frequency ( f 0 ), its harmonics and the modulation triplets were pre-whitened. Black lines show the location of f 0 and its harmonics, while the blue dash-dotted lines show the location of the half-integer frequencies. To visualize the distribution of the amplitudes of the HIFs we plotted black points on top of the frequencies originating from the PD e GLYPH&lt;11&gt; ect. The HIFs are significant up to 11 = 2 f 0 . Right panel: Small part of the CoRoT light curve of 962 clearly showing alternating maxima, i.e. , period doubling. Fitted maxima are plotted to guide the eye.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-1f09b442-0684-40ed-bd96-100117b4285a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-15eb2f6e-5bdf-4ef0-949c-56935683375a.jpg)

(see Fig. 7). Poretti et al. (2010) interpreted f 2 = 3 159 d : GLYPH&lt;0&gt; 1 and f 0 + f 2 = 5 279 d : GLYPH&lt;0&gt; 1 as combination frequencies, but in our context the period doubling interpretation seems to be more plausible, especially that it is found in all possible places from 1 = 2 f 0 to 17 = 2 f 0 2 . In our interpretation these frequencies are located closer than 1% to 3 = 2 f 0 = 3 178 d : GLYPH&lt;0&gt; 1 and 5 = 2 f 0 = 5 297 d : GLYPH&lt;0&gt; 1 , thus we identify them as HIFs.

et al. 2010). Its combinations with the fundamental mode frequency and its harmonics ( f 1 + k f 0 ; k = 0 1 2 ; ; ::: ) are clearly visible in Fig. 7 with high amplitudes in between the main harmonics and the HIFs. As in the previous cases, temporal variability is easily seen for both the 3 = 2 f 0 half-integer frequency and the additional frequency f 1 = 3 63088 d : GLYPH&lt;0&gt; 1 (upper right panel of Fig. 6).

We note that there is another additional frequency, f 1 = 3 63088 d : GLYPH&lt;0&gt; 1 , which might be the second radial overtone (Poretti

2 We make an exception here and retain the original designation of f 1 and f 2 from Poretti et al. (2010)

CoRoT0100881 648 is a Blazhko star exhibiting a typical asymmetric RR Lyrae light variation with a pulsation period of 0.607186 d and a Blazhko cycle of 59.77 d (Szab´ o et al. 2009). The modulation is sinusoidal and symmetric (Fig. 4).

Fig. 6. Temporal variation of the amplitude of additional frequencies in CoRoT Blazhko RR Lyrae stars. The CoRoT photometric light curves were split into eight chunks, the dominant pulsational frequency, its harmonics and the modulation side peaks were prewhitened. The Fourier spectrum in the vicinity of the additional frequencies (denoted by small verticals lines) are shown in the upper parts of the sub-figures. Upper left: Temporal variation of the half-integer frequency at 3 = 2 f 0 = 4 1916 d : GLYPH&lt;0&gt; 1 and the independent frequency f 0 = 4 03265 d : GLYPH&lt;0&gt; 1 in the Blazhko RR Lyrae 962 . Upper right: The same for the CoRoT Blazhko RR Lyrae star 793 . The frequencies are 3 = 2 f 0 = 3 159 d : GLYPH&lt;0&gt; 1 and f 1 = 3 63088 d : GLYPH&lt;0&gt; 1 . Lower left: The variation of f 0 = 2 389287 d : GLYPH&lt;0&gt; 1 additional frequency in the blended CoRoT Blazhko RR Lyrae star 544 . Lower right: The same for the Blazhko star 363 , where f 0 = 2 98400 d : GLYPH&lt;0&gt; 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-a7465c2a-2ca2-4b9a-b329-79ce8872f1f3.jpg)

The typical RRab light variation is diluted by a string of three close stellar companions to W-NW and a fainter star to the North (Fig. 8). The flux from these stars are contained by the CoRoT aperture together with that of the star 648 , and their presence is revealed by our ground-based observations that easily resolved the near-by contaminating sources. These follow-up observations allow us to make a crude estimate on the intrinsic light variation amplitude of this Blazhko star. One V image was taken in July 2008 (HJD 2454676.44742) with the Konkoly Observatory 1m RCC telescope which corresponds to a pulsational minimum and is slightly past the Blazhko-maximum ( = 0 06, where : denotes the Blazhko phase). Based on the fluxes of the contaminating stars (assumed to be constant) and the variable in its minimumandusing the apparent total amplitude we calculate the true light variation to be 0 : m 36 in the Blazhko maximum and 0 m : 30 in the Blazhko minimum. This is still a bit lower than the amplitude of the light variation of a normal, non-blended Blazhko RRab star, but to reconstruct the undiluted light variation more precisely, much longer follow-up observations would be necessary. Other fainter stars have negligible e GLYPH&lt;11&gt; ect on our estimate.

We detected 10 harmonics of the 1 64694d : GLYPH&lt;0&gt; 1 frequency and a triplet around f 0 and 2 f 0 . The right peaks are much higher and are detected around all the harmonics. No additional frequency is detected above the 3 GLYPH&lt;27&gt; level. Also, we could not detect any signs of the period doubling phenomenon. The frequency table is available at Appendix D.

CoRoT0101503 544 is a Blazhko star with an asymmetric light variation. Its pulsation period is 0.605087 d and a Blazhkovariation of 25.60 d is clearly seen (Fig. 4). The color light curves are preserved for this object. The star exhibits a typical skewed RRab light curve, with its amplitude diminished to 0 : m 021 in red, 0 : m 045 in green and 0 : m 052 in blue passbands. Based on the CoRoT data this object is most probably an intrinsically highamplitude RR Lyrae star. This conclusion is supported by the relatively high modulation (e.g. 0 : m 028 in the green light compared to the overall light variation amplitude 0 : m 045). We note that the color amplitude ratios are typical for an RR Lyrae variable. The L1 (level-1) contamination value given by the ExoDat catalog is 0 13871 : GLYPH&lt;6&gt; 0 00663, which might be a too low estimate. :

Fig. 8. Left: EXODAT image showing the aperture applied to 648 (upper panel) and 544 (lower panel). Right: 1' x 1' V image of the targets taken by the 1m RCC telescope at Piszk´ es-tet˝. o North is to the right and East is to the bottom. The variable star is located between the thick lines. We transformed our frames to match the EXODAT images except a slight (30 degrees) rotation. 648 was caught around pulsational minimum, while 544 around pulsational maximum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-3b86113d-2687-46c3-a691-0f49d4390c8e.jpg)

CCD frames taken by one of us (J. Benk˝ o) with the Konkoly Observatory 1m RCC telescope and a Johnson V filter allowed us to reconstruct the true light variation of this heavily blended object. The contaminating star is brighter by 2 : m 0 in V compared to the brightest phase of 544 . One image was taken in June 2008 (HJD 2454639.49092) which corresponds to a pulsational maximum and is slightly past the Blazhko-maximum ( = 0 10). : Using a similar process that was applied in the case of 648 (and assuming that the Johnson V filter is reasonably close to the CoRoT green passband) we calculate the true light variation to be 0 : m 39 in the Blazhko maximum and 0 m 23 in the Blazhko : minimum in the V passband. The frequency content of 544 can be found in Appendix E.

After a few days of normal sampling (8-min) the observations were switched to the oversampling mode (32-sec). Before embarking on the frequency analysis we re-sampled the oversampled data to normal sampling to get a more tractable number of points (23 915 vs. 351 086). We checked that no information was lost in the relevant low-frequency range which is the main focus of this study.

High-precision space-based observations frequently show additional frequencies in RR Lyrae stars. This Blazhko star is no exception, it also has additional periodicities in the frequency spectrum. One of them is f 0 = 2 389287 d : GLYPH&lt;0&gt; 1 . The lower left panel of Fig. 6 shows the temporal variation of the f 0 . Again, one can easily recognize the time-variability of the appearance of this frequency peak during the CoRoT observations.

We have not found any signatures pointing to period doubling in the frequency spectrum neither in the white light, nor in the color observations, which is not surprising, given the blended nature of our target. However, a close inspection of the light curve reveals that there are sections where the alternating pulsating cycles are clearly visible (Fig. 9). We checked that the frequency spectrum does not contain any other additional frequency peaks with high enough amplitude to cause the observed variation in the maxima from cycle to cycle. Clearly, the 0.15 mmag f 0 can not cause fluctuation exceeding several mmag as seen in Fig. 9. Interestingly, CoRoT's detection capability allowed us to reveal a temporary nature of the period doubling phenomenon in this heavily blended CoRoT Blazhko RR Lyrae star.

CoRoT0105288 363 is an RRLyrae star showing a strong Blazhko modulation with a period of 35.06 d. The lower right panel of Fig. 4 shows that the maxima of the pulsational cycles do not move considerably in sharp contrast with 962 . The amplitude of the Blazhko modulation is not constant. A detailed study was presented by (Guggenberger et al. 2011). After reanalyzing the CoRoT light curve, we have not found any sign of period doubling either in the light curve or in the frequency spectrum down to 0.35 mmag in accordance with Guggenberger et al. (2011). There is an additional frequency in the star ( f 0 = 2 98400 d : GLYPH&lt;0&gt; 1 ) with a frequency ratio of 0.591, which can be identified as the second radial overtone. Despite its low amplitude and the cluttered frequency spectrum we tried to analyze its temporal variation (Fig. 6, lower right panel). Although there are hints that the amplitude of this frequency peak varies in time as well, we can not draw firm conclusions about its temporal variability due to the di GLYPH&lt;14&gt; culties mentioned. Also, there is a non-negligible chance that this periodicity comes from a di GLYPH&lt;11&gt; erent, nearby source (Guggenberger et al. 2011). According to Guggenberger et al. (2012), the star also shows a frequency around its first radial overtone ( f 1 = 2 3793 d : GLYPH&lt;0&gt; 1 , frequency ratio 0.741) and another one which is most probably a nonradial mode 3 ( f nr = 2 4422 d : GLYPH&lt;0&gt; 1 , f 0 = f nr = 0 722). Due to their low am-: plitude we did not try to analyze their temporal behavior in this study.

CoRoT0103922 434 : This is a new Blazhko RRab star (also known as V922 Oph) that will be discussed in detail in Poretti et al. (in prep). It's period is 0.5413828 days and shows a Blazhkomodulation of roughly 54.5 days, thus CoRoT LRc04 observations cover two modulation cycles.

Here we only note that despite the fact that we do not see half-integer frequencies in the frequency spectrum of this star, there are hints that the period doubling might be present temporarily. Namely, we find 6 - 10 consecutive pulsational cycles in the light curve that are alternating (Fig. 10). One is found around CJD 3490 and lasts only for 10 pulsational cycles. The other one is at the end of the CoRoT observations (starting from CJD 3548). Such regularity can not be seen in any of the other RRab stars that do not show PD. The amplitudes of the additional frequencies are not su GLYPH&lt;14&gt; cient to cause as large as many hundredths of a magnitude alternation from maximum to maximumeven in the most favorable (constructive interference) case. In addition, finding alternating maxima by chance more than once has a very low probability. 4 We also note that in RR Lyrae

3 Denoted by f N in the original paper.

4 To be more specific: in a very simplistic approach, the probability of n jump (up or down) is 1 = (2 n GLYPH&lt;0&gt; 1), since the random variable X of the height of the pulsational maximum follows a binomial distribution assuming an equal probability ( p = 1 2) of an up or down jump. Thus, = the chance configuration of n = 8 switches in the left panel of Fig. 10 is 1 127 = GLYPH&lt;25&gt; 0 0079. :

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-a23689b8-cbb1-431e-ba7a-ff1db0b1c202.jpg)

magnitude

∆

CJD

Fig. 9. Sections of light curve where period doubling is present in the blended Blazhko CoRoT RRab star 544 . Even and odd maxima are marked with di GLYPH&lt;11&gt; erent symbols. Original data points are denoted by small black points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-0d4907af-c64b-409a-a895-9a3e5bc63871.jpg)

and other Kepler PD-stars the temporal nature of the period doubling is also clearly demonstrated (Szab´ et al. 2010). Thus, we o conclude that 434 is a Blazhko RRab star showing the period doubling temporarily. One can notice that there are switches between the di GLYPH&lt;11&gt; erent branches connecting even and odd maxima (Fig. 10). The switch between the two branches of the lower and higher maximum strings is a definite sign of a more complex dynamical behavior besides the high-order resonance causing the period doubling. It points to the presence of additional resonances (see e.g. Moln´ ar et al. (2012)) or low-dimensional chaos (Plachy et al. 2013).

Although there are indications that the second radial overtone is excited in this star with f 2 = 3 165859 d : GLYPH&lt;0&gt; 1 , we examined only the strongest additional frequency f 0 = 2 612196 d : GLYPH&lt;0&gt; 1 and the result is presented in the upper right panel of Fig.11. There is a clear amplitude variation (lower panel). The figure was made after pre-whitening with the fundamental mode f 0 , its harmonics and the modulation side peaks. Large number of combination frequencies ( k f 0 + nfm ) are present as well, suggesting that these signals come from the star itself. The frequency ratio is f 0 = f 0 = 0 707091, thus : f 0 can not be a radial mode.

## 4.3. RRc stars

CoRoT0105036 241 : We have found two RR Lyrae stars predominantly pulsating in the first radial overtone (RRc) by CoRoT, one of them is 241 . The phased light curve can be seen in the left panel of Fig. 12. Besides the first overtone mode ( f 1 = 2 68153 d : GLYPH&lt;0&gt; 1 ), this RRc star has an additional frequency f 0 = 4 37783 d : GLYPH&lt;0&gt; 1 with a frequency ratio f 1 = f 0 = 0.613. This frequency ratio occurs surprisingly often in overtone RR Lyrae and Cepheid stars, see e.g. , Moskalik (2013). 241 represents another member of this class of classical pulsators. If - as it is conceivable f 1 belongs to the first radial overtone then f 0 can not be a radial mode, it is most probably a nonradial one.

The Fourier spectrum and sequential pre-whitening with the highest amplitude frequencies in 241 are displayed in Fig. 13. After successively removing the dominant pulsation frequency, its harmonics, and the additional f 0 frequency, two sets of side peaks appear. We denote the corresponding frequency di GLYPH&lt;11&gt; erences with f m and f b in Fig. 13 and in Appendix F. The first one ( f m) may be a long-period (longer than the length of our data set) Blazhko modulation, although more data would be necessary to establish its exact nature, while the second one ( f b ) most prob- ably appears because a clear amplitude variability is seen in the additional frequency, f 0 (see below).

After removing all the previously mentioned frequencies there are still a lot of remaining peaks in the spectrum, the highest among them is f 00 = 2 345174 d : GLYPH&lt;0&gt; 1 with an amplitude of 0.57 mmag. However, instead of (over)interpreting the data we prefer to stop at this level. The average value of the 'grass' of the residual frequency spectrum is 0.09 millimagnitude. A detailed discussion of the frequency spectrum can be found in Appendix F.

We computed the temporal variation of the amplitude of f 0 by using eight bins. The lower left panel of Fig. 11 shows that the amplitude of f 0 very clearly varies on long time-scales. In the same figure we also plotted the close neighborhood of f 0 in the frequency spectrum corresponding to the bins for which the amplitude was calculated (upper panel). It is interesting to see the structural variation of the peak. Sometimes the frequency peak is split into two separate peaks. The variation over the 145 day-interval is striking.

CoRoT0105735 652 is similar to 241 in many aspects. Its phased light curve is shown in the right panel of Fig. 12. However, this star has a B5 II spectrum in the CoRoT observation log Release 13 (COROT Team 2014). This and the presence of many low amplitude frequencies (that may be nonradial modes) might imply a GLYPH&lt;12&gt; Cep scenario. On the other hand, the light curve shape resembles that of an RRc star with the characteristic bump just before maximum light, though some GLYPH&lt;12&gt; Cep variables show a similar feature (Sterken et al. 1987). The CoRoT classification algorithm (Debosscher et al. 2009) assigns a probability of more than 80% for an RRc variable, while gives less than 10% chance for the GLYPH&lt;12&gt; Cep variation. In addition, the frequency ratio f 1 = f 0 = 0.615 is seen frequently in recent space-based photometric observations (e.g. Moskalik 2013), which makes it likely that the star is a first-overtone RR Lyrae star, hence we stick to this scenario.

The temporal variation of the amplitude of the f 0 frequency is plotted in lower right panel of Fig. 11. The f 0 frequency shows similar structural variations as its sibling, the other CoRoT RRc star, 241 that we discussed above. Although this star does not show as large amplitude variation as 241 , the change in the structure of this secondary frequency is obvious. Note that while the timescale of the variation is similar, the variations in the two stars are not identical, and other frequencies do not show the same effect, so we conclude that they can not come from instrumental e GLYPH&lt;11&gt; ects. We also checked the well-behaved spectral window to

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-ce5be52d-7fd4-49dd-a1ba-97f48f884ee5.jpg)

CJD

Fig. 10. Same as Fig. 9, but for the Blazhko CoRoT RRab star 434 . The switch back and forth between the two branches in the right panel is a clear signal points to more complex dynamical behavior.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-89f9576e-d29a-469b-b9d2-9a3559f4f156.jpg)

see whether the first side lobe is able to cause any trouble in combination with f 0 , but this possibility can be safely excluded. The frequency table of the star is shown in Appendix G.

## 4.4. CoRoT RRd star

presence of period doubling has not been recognized. These objects are all Blazhko RRab stars showing either characteristic alternating maxima (or pulsating cycles) in the time domain, or half-integer frequencies in between the dominant pulsation mode and its harmonics (Szab´ o et al. 2010), or both (see Table 4 for a summary).

CoRoT0101368 812 : While the original Kepler field does not contain any classical double-mode RR Lyrae stars, this RRd star has been found in CoRoT's LRc01 run and analyzed by Chadid (2012). The period of the fundamental mode is 0.4880408 d, while that of the first overtone is 0.3636016 d. The period ratio is 0.7450 and the amplitude ratio A 1 = A 0 is 2.7055, the first overtone being the dominant mode. We note in passing that RR Lyrae, the eponym of its class is an example of a modulated RR Lyrae pulsating in the fundamental mode and exhibiting the first overtone mode in addition (Moln´ ar et al. 2012), with an exceedingly tiny amplitude (non-classical double-mode RR Lyrae).

It is of interest to investigate whether classical double-mode RRLyrae show the period doubling, because of the lack of RRds in the original Kepler field. Based on the CoRoT data we see no half-integer frequencies in the Fourier-spectrum of 812 down to 0.07 mmag. We mention here that MOST observations of AQLeo, another classical RRd star do not show PD frequencies, either (Gruberbauer et al. 2007).

We examined the additional frequency f 0 = 4 4786159 d : GLYPH&lt;0&gt; 1 (Chadid 2012) ( f 2 in the original paper) and the result is presented in the upper right panel of Fig.11. The figure was made after pre-whitening with the fundamental mode f 0 (and its harmonics) the first overtone f 1 (along with all harmonics) and all the visible linear combinations of these. The f 0 periodicity shows a clear amplitude variation. We tested that the same variation is present in the original (un-prewhitened) spectrum as well, so the e GLYPH&lt;11&gt; ect is not caused by the pre-whitening process. We also checked that the observed change is not caused by any variation of the f 0 or f 1 themselves, or an interaction with the spectral window side-lobes.

## 4.5. Brief summary of the results

After discussing individual objects it is worth summarizing the most important results:

We discovered signs of period doubling in four CoRoT RRLyrae stars by re-analyzing their light curves. These are 962 , 793 544 , and 434 . Although these stars were known before, the

In Tab. 4 we list additional frequencies reported here for the first time, as well as a few previously found ones. These frequencies can be interpreted as

- -the second radial overtone in Blazhko RRab stars (or at least frequencies around its expected location with a frequency ratio GLYPH&lt;25&gt; 0 58) in case of : 962 793 , and 434 . Here we investigated only the one belonging to 793 in detail,
- -the well-known frequency with a frequency ratio of 0 61 : most probably corresponding to a nonradial mode in RRc and RRd stars: 241 , 652 , 812 (we use the notation f X throughout this work) , and
- -other nonradial modes in 962 544 , and 434 .

We emphasize that the list in Tab. 4 is by no means exhaustive or complete, since we do not investigate and discuss all the significant additional frequencies here, only the most prominent ones. Interestingly, based on this work and Nemec et al. (2011), nonBlazhko RRab stars do not show additional periodicities down to the exquisitely low amplitude limits provided of CoRoT and Kepler .

Another interesting conclusion of this work is that in many cases the amplitudes of additional frequencies vary over time on shorter time scales than the typical length of the CoRoT runs, irrespective of their nature (HIF, second radial overtone, or hypothesized nonradial modes). We found clear temporal variability of the amplitudes of a frequency peak close to the second radial overtone of 793 , the probable nonradial modes with a frequency ratio of 0.61 in 241 (RRc) and 812 (RRd), and other possibly nonradial modes in 544 and 434 . Variability is also obvious in the structure of the frequency peak f 0 = 5 82484 d : GLYPH&lt;0&gt; 1 in 652 (another RRc, f 0 = f = 0 615). Only marginal evidence : is found for temporal variability of the amplitudes of a possibly nonradial mode and the second radial overtone in the Blazhkomodulated RRab stars 962 and 363 , respectively.

Fig. 11. Temporal variability of additional frequencies in RRab, RRd and RRc stars. Upper left: Lower panel: Amplitude variation of the additional f 0 = 2 61219586 d : GLYPH&lt;0&gt; 1 frequency of CoRoT Blazhko RR Lyrae star 434 . The fundamental mode frequency, its harmonics and the side peaks related to the modulation (and the modulation frequency itself) were removed. Upper panels: the vicinity of f 0 in the frequency spectrum, lower panels: amplitude variation. Upper right: The same for the temporal variability of f 0 = 4 4786159 d : GLYPH&lt;0&gt; 1 CoRoT RRd star 812 . The first overtone frequency and its harmonics were removed. Lower left: Temporal variation of the f 0 = 4 37783 d : GLYPH&lt;0&gt; 1 frequency of CoRoT RRc star 241 . The first overtone frequency and its harmonics were removed. Lower right: The same as for the RRc star 241 , but for 652 . Note that the amplitudes are smaller than in 241 , but the variation of the structure of the f 0 = 5 82484 d : GLYPH&lt;0&gt; 1 peak is strikingly similar.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-19b0e181-09dc-4709-8375-5a0808b7b609.jpg)

Fig. 12. Phased light curves of two CoRoT RRc stars. Left: 241 . Right: 652 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-7e43ac32-056f-46de-8604-dc5b47a720c7.jpg)

## 5. Discussion and conclusions

In this work we added four non-modulated RRab stars to the list of RR Lyrae not showing PD with CoRoT's precision.

## 5.1. Period doubling

We confirm earlier results (Szab´ o et al. 2010) that period doubling occurs only in Blazhko-modulated RR Lyrae stars, and none of the non-modulated RRab stars shows this phenomenon.

Although the amplitude of the HIFs is empirically expected to be highest in the [ f 0 ;2 f 0 ] frequency interval, we nevertheless checked the whole frequency spectrum of each objects when searching for the PD. Indeed, the HIFs are significant up to 11 2 = f 0 in case of 962 and 17 2 = f 0 in 793 . In two cases ( 544 and

Fig. 13. Frequency spectrum of the star 241 with sequential pre-whitening. Some characteristic frequencies are labeled with their identification. Sat refers to frequencies connected to the orbit of the CoRoT satellite. a.) The original spectrum dominated by the first overtone pulsation mode ( f 1 ). b.) After pre-whitening with f 1 and its harmonics. c.) The result of removing f 0 and all combination terms involving f 1 and f 0 . d.) Pre-whitening with all f 1 GLYPH&lt;6&gt; k GLYPH&lt;3&gt; f m and f 0 GLYPH&lt;6&gt; k GLYPH&lt;3&gt; f m frequencies. e.) After removing all significant f 1 GLYPH&lt;6&gt; k GLYPH&lt;3&gt; f b and f 0 GLYPH&lt;6&gt; k GLYPH&lt;3&gt; f b frequencies. f.) Removing another 13 frequencies including those originating from the satellite and third order linear combinations of f f 1 ; f 0 ; f m ; f b g .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0653v1-figure-67f9eb4b-8753-4f35-bd03-53deb0df81c6.jpg)

434 ) only the alternating maxima betray the presence of PD. In 434 the reason may be twofold: the PD is rather weak and it is present only for a very brief time period. 544 is heavily blended, and we had not expected to see anything in the frequency spectrum beyond a few harmonics of the fundamental mode and some modulation side peaks. On the contrary, with CoRoT data we were able to uncover PD in this blended Blazhko RR Lyrae star (its total amplitude is only about 0 : m 02), nicely demonstrating the capabilities of high duty-cycle, extreme-precision space photometry. We checked the frequency spectrum of 544 and 434 to see whether any single additional frequency or a combination of those can cause the conspicuous alternation of consecutive maxima. It turned out that the amplitudes of these frequencies are much lower than what would be required to cause the observed e GLYPH&lt;11&gt; ect, thus we conclude that the hypothesis of the presence of PD is plausible in these cases.

As seen in the Kepler RRLyrae sample, the strength of period doubling is changing with time in the CoRoT PD stars as well, and in most cases it crops up only temporarily for a few pulsational cycles. For the estimation of the occurrence rate of the period doubling e GLYPH&lt;11&gt; ect in Blazhko and possible non-Blazhko stars it is very important to monitor RR Lyrae stars with highprecision photometry (preferably from space) for a long time (several Blazhko-modulation periods, meaning several months or years) uninterrupted because of the time variability of the phenomenon. CoRoT and Kepler data sets are ideal for this purpose, because both missions delivered many months- and years-long high-precision, uninterrupted observations, respectively. The four stars with PD out of six modulated RRab targets is entirely consistent with the results of (Benk˝ et al. 2014) who o found PD in nine out of 15 Blazhko-modulated Kepler stars. This can be augmented to 10 16, if RR Lyrae, the prototype is / accounted for (Szab´ o et al. 2010). This means that roughly two

Table 4. Summary of the main results: period doubling search and additional frequencies that were investigated in detail along with their most probable identification. ' &lt; ' denotes upper amplitude limit. For the only RRd star in the sample half-integer frequencies corresponding to both radial pulsation modes were searched for.

| CoRoT ID                                                    | Ampl. of 1 : 5 f 0 [mmag]   | alternating maxima   | Type    | PD yes / no   | f 0 or f 1 [d GLYPH<0> 1 ]   | additional freq. f 0 [d GLYPH<0> 1 ]   | frequency ratio f = f 0   | identification   |
|-------------------------------------------------------------|-----------------------------|----------------------|---------|---------------|------------------------------|----------------------------------------|---------------------------|------------------|
| 0101370 131 0101315 488 0103800 818 0104315 804 0100689 962 | < 0 : 30                    | n                    | RRab    | n             | 1 : 61464                    |                                        |                           |                  |
|                                                             | < 0 : 10                    | n                    | RRab    | n             | 2 : 06057                    |                                        |                           |                  |
|                                                             | < 0 : 07                    | n                    | RRab    | n             | 2 : 14622                    |                                        |                           |                  |
|                                                             | < 0 : 15                    | n                    | RRab    | n             | 1 : 38538                    |                                        |                           |                  |
|                                                             | 3 : 21                      | y                    | RRab Bl | y             | 2 : 80902                    | 4.03265                                | 0.697                     | nr               |
| 0101128 793                                                 | 2 : 44                      | y                    | RRab Bl | y             | 2 : 11895                    | 3.63088                                | 0.584                     | O2               |
| 0100881 648                                                 | < 0 : 10                    | n                    | RRab Bl | n b           | 1 : 64694                    |                                        |                           |                  |
| 0101503 544                                                 | < 0 : 05                    | y                    | RRab Bl | y b           | 1 : 65266                    | 2.38929*                               | 0.692                     | nr               |
| 0105288 363                                                 | < 0 : 35                    | n                    | RRab Bl | n             | 1 : 76230                    | 2.98400                                | 0.591                     | O2               |
| 0103922 434                                                 | < 0 : 30                    | y                    | RRab Bl | y             | 1 : 84712                    | 2.61220*                               | 0.707                     | nr               |
| 0105036 241                                                 | < 0 : 20                    | n                    | RRc     | n             | 2 : 68153                    | 4.37783*                               | 0.613                     | f X              |
| 0105735 652                                                 | < 0 : 13                    | n                    | RRc     | n             | 3 : 58218                    | 5.82484*                               | 0.615                     | f X              |
| 0101368 812                                                 | < 0 : 20                    | n                    | RRd f 0 | n             | 2 : 04901                    | -                                      | -                         | -                |
|                                                             | < 0 : 15                    | n                    | RRd f 1 | n             | 2 : 75026                    | 4.47862                                | 0.614                     | f X              |

Notes. The superscript 'b' denotes blended pulsators, a configuration that prevents detection of low amplitude features. Identifications: O2: second radial overtone, nr: nonradial mode, fx: ubiquitous frequency in RRc stars with a frequency ratio GLYPH&lt;25&gt; 0.61 (Moskalik 2014). '*' denotes new discoveries, other frequencies were reported in previous publications, see Tab. 1 for the relevant references.

Blazhko-modulated RRab stars out of three show the perioddoubling phenomenon. One is tempted to speculate that because of the finite length of space data (60-150 days in case of CoRoT, 4 years for Kepler ) and the temporary presence of PD it is possible that the occurrence rate can be higher. To take this argument to the extreme, it is possible that all Blazhko-modulated stars would show period doubling if long enough times series data were taken.

The connection between the Blazhko modulation and the period doubling is obvious. The high occurrence rate of PD in Blazhko-modulated stars stresses the possible underlying physical connection of the two dynamical phenomena. Indeed, the 9:2 resonance between the fundamental mode and the 9 th overtone which explains the period doubling (Szab´ o et al. 2010; Koll´ ath et al. 2011), might be the culprit for causing the Blazhko modulation itself, as well (Buchler &amp; Koll´ ath 2011). Given that the resonance paradigm is partly based on hydrodynamic models and partly on the successful and simple amplitude equation formalism, in our opinion currently this is the most plausible explanation for the century-old Blazhko enigma. It is backed up by the results of this paper.

found (Gruberbauer et al. 2007), but soon with the advent of regular space photometric observations - Blazhko-modulated RRab (Chadid et al. 2010; Guggenberger et al. 2012) and RRc stars (Moskalik 2013, 2014) followed. It is especially interesting that most of the RRc stars that were observed from space show periodicities with a frequency ratio of f 1 = f 0 GLYPH&lt;25&gt; 0 61 with : the first radial overtone mode (Moskalik 2014). These additional frequencies seem to be present in all stars in our sample showing normal amplitude first overtone pulsation, i.e. , RRc and RRd stars (see Tab. 4, where we identified them with f X ). The most plausible explanation for their origin is nonradial modes (Dziembowski 2012). Interestingly, modulated RRab stars also show various, low-amplitude additional periodicities (Benk˝ o et al. 2010; Chadid et al. 2010; Moln´r et al. 2012; Guggenberger a et al. 2012).

Interestingly, alternating cycles are found in other pulsating star types, such as RV Tauri of which the alternating deep and shallow minima are characteristic and also in BL Her stars (Smolec et al. 2012). Surprisingly, white dwarfs may show alternating cycles, as well (Papar´ et al. 2013). Although in some o cases this may be a result of an interaction of several independent modes in the frequency spectrum, at least in one case a period doubling bifurcation event seems to be well-documented (Goupil et al. 1988), showing close resemblance of RR Lyrae star pulsational dynamics (Koll´ ath et al. 2011).

## 5.2. Time-dependent additional frequencies

With the thorough analysis of CoRoT RR Lyrae data we confirm earlier emerging trends that additional periodicities are ubiquitous in RR Lyrae stars. RRd stars were the first type of objects where frequencies not fitting in the radial eigenspectrum were

In this work we attempted to unveil the temporal behavior of the ubiquitously seen additional, low-amplitude frequencies that are usually attributed to nonradial oscillations by the virtue of the extended coverage and high duty cycle provided by CoRoT. We found that in almost all cases, where the brightness of the star, the data coverage and the crowdedness of the frequency spectrum allowed a detailed analysis, these frequencies showed amplitude variation over time, followed by a variation of the shape of the frequency peak. The structure of the peaks already suggests temporal variability of the amplitude and or frequencies, / since amplitude and or frequency modulation manifests itself as / side peaks around the corresponding frequency peaks (Benk˝ o et al. 2011). This is exactly what we see around the HIFs in many cases (Szab´ o et al. 2010), see also Fig. 5 in this work. Maybe this same variation occur in other additional frequencies, as well. It is possible that the amplitude and the structure of these additional frequency peaks vary because there are close, unresolved frequencies around them. We consider this possibility unlikely based on our experiences with the Kepler RRLyrae data that have better frequency resolution.

Finding a physical explanation for ubiquitous time variability of the additional frequencies, which at the same time does not a GLYPH&lt;11&gt; ect the dominant pulsational mode is challenging. We are tempted to think that the additional frequencies might be non-

radial modes. (i) Maybe these modes are not self-excited, and are only present when they experience resonances. That would explain nonstationarity, if resonance conditions governed by the stellar structure are not always met. This might be the case in Blazhko-modulated stars, where the period and stellar structure also varies (quasi)periodically (Jurcsik et al. 2009b; S´dor o et al. 2009; Szab´ o et al. 2010; Koll´ ath et al. 2011). The typical timescale of the variation we found is several tens of days, although the variation may be seriously under-sampled because of the attainable frequency resolution. This might be compared to typical theoretical mode growth rates or timescale of interactions, as we see in stars showing period doubling (Koll´ ath et al. 2011). Anyhow, the variation is most probably a consequence of the interaction with the large-scale, long-term modulation, which may be a plausible explanation for the observed phenomena. However, this explanation is challenged by the fact that the amplitude and hence non-linearity in RRc stars is much lower than in RRab stars, still the same variation occur in these stars, as well. (ii) If not only one, but several closely-spaced nonradial modes are excited (Van Hoolst et al. 1998), they can interact causing a complex behavior in the frequency space. In this case not only the time scale of the excitation, but also that of the nonlinear interaction becomes relevant. While the found variations are quite diverse, the magnitude and the timescale of the newly found variations of the additional frequencies suggest a common origin. (iii) Rotational splitting is less likely, since we would not expect temporal variation unless other mechanism is at work. Invoking some other modes might come to the rescue if further nonlinear interactions are assumed, hence creating a complex, unresolvable (even with continuous CoRoT observations) pattern in the frequency domain. Detailed investigation of these mechanisms is beyond the scope of this paper, but we may conclude that the discovery of these variations leads us to an unexplored territory of fine details of RR Lyrae pulsation.

## 5.3. On the interpretation of the frequency spectra of Blazhko RRab stars

In Chadid et al. (2010) and Poretti et al. (2010) the frequency spectrum of the Blazhko-modulated 962 and 793 were interpreted as the result of the presence of independent, additional frequencies (presumably nonradial modes) and their combinations besides the well-known pulsational frequency, harmonics, modulation multiplets, the modulation frequency itself and its harmonics. Here, we o GLYPH&lt;11&gt; er a simpler explanation, since the presence of half-integer frequencies, as a result of the period doubling removes one independent frequency together with all the related combination terms.

In many Blazhko RR Lyrae stars, peaks appear around the expected frequencies of the second (and first) radial overtone (Benk˝ o et al. 2010). This is the case for 962 and 793 , as well. Here we identify the frequencies with a frequency ratio of 0.58 with the radial second overtone, while keeping in mind that the excitation of nonradial modes in the vicinity of radial overtones is also predicted and plausible (Dziembowski 1977; Van Hoolst et al. 1998).

Other significant frequencies were found by previous studies and this work (see in Tab. 4. denoted by 'nr') which do not fit the above described picture. Their frequency ratio ( f 0 = f 0 ) are about 0.7. As we mentioned it in Benk˝ o &amp; Szab´ o (2014) these frequencies can also be interpreted as linear combination: 2( f 2 GLYPH&lt;0&gt; f 0 ). Indeed, the components and the simplest combination ( f 0 , f 2 and f 2 GLYPH&lt;0&gt; f 0 ) of these frequencies are detectable for all Blazhko RRLyrae stars except for the heavily blended 544 and 648 .

The linear combination assumption simplifies the mathematical description, but is it a plausible physical explanation? Do these frequencies belong to radial modes? Not necessarily, since the amplitudes seem to contradict this scenario. For example, the amplitudes of f 2 GLYPH&lt;0&gt; f 0 and f 2 in 962 are about three times smaller than the amplitude of 2( f 2 GLYPH&lt;0&gt; f 0 ). Such a behaviour is highly unlikely for simple linear combinations. It is possible however, that f 2 and f 0 excite a nonradial mode ( f nr ) through a three-mode resonance f nr GLYPH&lt;25&gt; 2( f 2 GLYPH&lt;0&gt; f 0 ), in which case the amplitude of the excited mode can be higher. Similar e GLYPH&lt;11&gt; ect have been detected for GLYPH&lt;14&gt; Scuti stars, B,A,F stars and a peculiar roAp star (Balona et al. 2013; Breger 2014; Breger &amp; Montgomery 2014). All of these papers suggest the presence of coupled nonradial modes showing quasi-periodic amplitude and frequency variations similarly to our presented results. Further developments of multidimension hydrocodes, such as Geroux &amp; Deupree (2011) and Mundprecht et al. (2013) will eventually make possible to test this scenario.

## 6. Summary

The main results of this work can be summarized as follows:

- -The most comprehensive collection of CoRoT RRLyrae variables is presented to date, including new discoveries. The sample consists of all RR Lyrae sub-types: Blazhko and non-modulated RRab stars, two RRc and one RRd star. We publish epochs, periods, frequency tables and phased light curves for those object that have not been analyzed earlier.
- -Thanks to the CoRoT high cadence observations, we could detect 56 harmonics of the pulsation frequency corresponding the fundamental mode in the case of non-modulated 818 RRab star. To our knowledge this is the highest number of observed Fourier-terms describing an RR Lyrae light curve, hence the most precise one that will serve as a benchmark for model computations.
- -Period doubling is detected in CoRoT Blazhko RRab stars for the first time. We discovered brief sections of alternating maxima typical of PD e GLYPH&lt;11&gt; ect in 4 CoRoT RR Lyrae stars. It means that two out of three modulated RR Lyrae shows this dynamical phenomenon at least temporarily, in accordance with the Kepler RRLyrae statistics. Given the usually short time intervals where PD is detectable, the percentage can be even higher. The strong correlation of the PD occurrence with the Blazhko phenomenon and the fact that no PD was detected in non-modulated CoRoT and Kepler RRLyrae suggests a causal relation, such as the 'resonance paradigm' proposed by Buchler &amp; Koll´ ath (2011). In some cases the presence of PD o GLYPH&lt;11&gt; ers a simpler explanation of the frequency spectrum of Blazhko-modulated RR Lyrae.
- -Our work corroborates those recently found trends, that additional frequencies (most probably higher radial overtones and nonradial modes) are ubiquitous in all sub-types of RRLyrae stars (RRc, RRd, Blazhko RRab), except the non-modulated RRab pulsators. This is an extremely strict rule, since none of the non-Blazhko RRab stars observed by Kepler and CoRoT shows any additional frequency peaks beyond the dominant pulsational mode and its harmonics, while all the other types - except some blended objects with high-precision space-based photometric observations do show this feature. If the additional frequencies proved to be nonradial modes, then we anticipate that asteroseismology of RR Lyrae stars should be feasible in the near future.

- -We analyzed the temporal variability of additional frequencies for the first time in all sub-types of RR Lyrae stars based on the CoRoT sample. The amplitude or the shape of these frequencies vary in time in most cases where we could draw firm conclusions. This variability can be connected to the Blazhko-cycle in modulated stars, much like the half-integer frequencies signaling the period doubling are strongly variable. A physical explanation in this case may be the changes of the stellar structure during the Blazhko cycle, and the consequent close or near-miss encounters with di GLYPH&lt;11&gt; erent resonances between radial and or nonradial modes. That would / explain the temporal excitation of nonradial modes (Koll´ ath et al. 2011). In non-modulated stars, such as RRd and RRc stars, however, a di GLYPH&lt;11&gt; erent mechanism should be at work, and it is not clear at this point whether a common mechanism can explain the temporal variability of additional frequencies in all RR Lyrae sub-types. Through investigations similar to this work we may get closer to the understanding of the excitation mechanism and origin of these periodicities.

Future high-precision photometric missions will multiply the number of interesting RR Lyrae stars to be investigated in detail. NASA's continuing Kepler Mission (dubbed K2), (Howell et al. 2014), TESS (Ricker et al. 2014), and PLATO (Rauer et al. 2013) will provide hundreds-to-thousands continuous of RRLyrae light curves spanning from a few weeks (TESS), to couple of months (K2), to several years (PLATO) coverage. Prospective data sets from upcoming missions will shed new light on the occurrence of the Blazhko e GLYPH&lt;11&gt; ect, period doubling, additional radial and nonradial modes and other dynamical phenomena as a function of a broad range of stellar parameters. In the light of these prospects, we are entering a golden era of classical variable stars, and our observational data presented in this work will provide ample examples for detailed theoretical analysis.

Acknowledgements. This research has made use of the ExoDat database, operated at LAM-OAMP, Marseille, France, on behalf of the CoRoT Exoplanet pro-/ gramme. RSz, MP, and JMB acknowledge the support of the ESA PECS project No. 4000103541 11 NL KML. This project has been supported by the Hungarian / / / OTKA grant K83790 and the European Community's Seventh Framework Programme (FP7 2007-2013) under grant agreements no. 312844 (SPACEINN), / no. 269194 (IRSES ASK) and ERC grant agreement no. 338251 (StellarAges). / RSz wishes to thank the support from the J´ anos Bolyai Research Scholarship of the Hungarian Academy of Sciences. WW was supported by the Austrian Science Fonds (FWF P22691-N16). The authors thank Akos Gy˝ or ´ GLYPH&lt;11&gt; y, P´ eter P´ apics and L´szl´ Moln´ ar for their help with the ground-based observations. a o

## References

A er, L., Micela, G., Favata, F., &amp; Flaccomio, E. 2012, MNRAS, 424, 11 GLYPH&lt;11&gt;

Baglin,

A.,

Auvergne,

M.,

Boisnard,

L., et

al.

2006, in

COSPAR Meeting,

Vol. 36, 36th COSPAR Scientific Assembly, 3749

- Balona, L. A., Catanzaro, G., Crause, L., et al. 2013, MNRAS, 432, 2808
- Benk˝ o, J. M., Kolenberg, K., Szab´, R., et al. 2010, MNRAS, 409, 1585 o
- Benk˝ o, J. M., Plachy, E., Szab´, R., Moln´ ar, L., &amp; Koll´th, Z. 2014, ApJS, 213, o a 31
- Benk˝ o, J. M. &amp; Szab´, R. 2014, in IAU Symposium, Vol. 301, IAU Symposium, o
- ed. J. A. Guzik, W. J. Chaplin, G. Handler, &amp; A. Pigulski, 383-384

Benk˝ o, J. M., Szab´, R., &amp; Papar´, M. 2011, MNRAS, 417, 974

o

o

- Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977
- Breger, M. 2014, in IAU Symposium, Vol. 301, IAU Symposium, ed. J. A. Guzik,
- W. J. Chaplin, G. Handler, &amp; A. Pigulski, 93-100
- Breger, M. &amp; Montgomery, M. H. 2014, ApJ, 783, 89
- Buchler, J. R. &amp; Koll´ ath, Z. 2011, ApJ, 731, 24
- Chadid, M. 2012, A&amp;A, 540, A68
- Chadid, M., Benk˝ o, J. M., Szab´, R., et al. 2010, A&amp;A, 510, A39 o
- Chadid, M., Perini, C., Bono, G., et al. 2011, A&amp;A, 527, A146
- Chevalier, C. &amp; Ilovaisky, S. A. 1991, A&amp;AS, 90, 225
- COROT Team. 2014, VizieR Online Data Catalog, 1, 2028
- Debosscher, J., Sarro, L. M., L´pez, M., et al. 2009, A&amp;A, 506, 519 o Dziembowski, W. 1977, Acta Astron., 27, 95 Dziembowski, W. A. 2012, Acta Astron., 62, 323 Dziembowski, W. A. &amp; Mizerski, T. 2004, Acta Astron., 54, 363 Geroux, C. M. &amp; Deupree, R. G. 2011, ApJ, 731, 18 Goupil, M. J., Auvergne, M., &amp; Baglin, A. 1988, A&amp;A, 196, L13 Gruberbauer, M., Kolenberg, K., Rowe, J. F., et al. 2007, MNRAS, 379, 1498 Guggenberger, E., Kolenberg, K., Chapellier, E., et al. 2011, MNRAS, 415, 1577 Guggenberger, E., Kolenberg, K., Nemec, J. M., et al. 2012, MNRAS, 424, 649 Howell, S. B., Sobeck, C., Haas, M., et al. 2014, ArXiv e-prints Jurcsik, J., S´dor, A., Hurta, Z., et al. 2008, MNRAS, 391, 164 o ´ Jurcsik, J., S´dor, A., Szeidl, B., et al. 2009a, MNRAS, 400, 1006 o ´ Jurcsik, J., S´dor, A., Szeidl, B., et al. 2009b, MNRAS, 393, 1553 o ´ Jurcsik, J., Szeidl, B., S´dor, A., et al. 2006, AJ, 132, 61 o ´ Kolenberg, K., Szab´ o, R., Kurtz, D. W., et al. 2010, ApJ, 713, L198 Koll´ ath, Z. 1990, Konkoly Observatory Occasional Technical Notes, 1, 1 Koll´ ath, Z., Moln´r, L., &amp; Szab´, R. 2011, MNRAS, 414, 1111 a o Le Borgne, J. F., Poretti, E., Klotz, A., et al. 2014, MNRAS, 441, 1435 Lenz, P. &amp; Breger, M. 2005, Communications in Asteroseismology, 146, 53 Moln´ ar, L., Koll´th, Z., Szab´, R., et al. 2012, ApJ, 757, L13 a o Moskalik, P. 2013, in Astrophysics and Space Science Proceedings, Vol. 31, Stellar Pulsations: Impact of New Instrumentation and New Insights, ed. J. C. Su´ arez, R. Garrido, L. A. Balona, &amp; J. Christensen-Dalsgaard, 103 Moskalik, P. 2014, in IAU Symposium, Vol. 301, IAU Symposium, ed. J. A. Guzik, W. J. Chaplin, G. Handler, &amp; A. Pigulski, 249-256 Mundprecht, E., Muthsam, H. J., &amp; Kupka, F. 2013, MNRAS, 435, 3191 Nemec, J. M., Smolec, R., Benk˝ o, J. M., et al. 2011, MNRAS, 417, 1022 Papar´ o, M., Bogn´ ar, Z., Plachy, E., Moln´r, L., &amp; Bradley, P. A. 2013, MNRAS, a 432, 598 Papar´ o, M., Chadid, M., Chapellier, E., et al. 2011, A&amp;A, 531, A135 Papar´ o, M., Szab´, R., Benk˝, J. M., et al. 2009, in American Institute of Physics o o Conference Series, Vol. 1170, American Institute of Physics Conference Series, ed. J. A. Guzik &amp; P. A. Bradley, 240-244 Plachy, E., Koll´ ath, Z., &amp; Moln´r, L. 2013, MNRAS, 433, 3590 a Poretti, E., Papar´, M., Deleuil, M., et al. 2010, A&amp;A, 520, A108 o Rauer, H., Catala, C., Aerts, C., et al. 2013, ArXiv e-prints Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, ArXiv e-prints Skarka, M. 2014, A&amp;A, 562, A90 Smolec, R., Soszy´ nski, I., Moskalik, P., et al. 2012, MNRAS, 419, 2407 S´ odor, A., Jurcsik, J., &amp; Szeidl, B. 2009, MNRAS, 394, 261 ´ Sterken, C., Young, A., &amp; Furenlid, I. 1987, A&amp;A, 177, 150 Szab´ o, R., Koll´th, Z., Moln´r, L., et al. 2010, MNRAS, 409, 1244 a a Szab´ o, R., Papar´, M., Benk˝, J. M., et al. 2009, in American Institute of Physics o o Conference Series, Vol. 1170, American Institute of Physics Conference Series, ed. J. A. Guzik &amp; P. A. Bradley, 291-293 Van Hoolst, T., Dziembowski, W. A., &amp; Kawaler, S. D. 1998, MNRAS, 297, 536 Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023

Table A.1. Frequency table of the newly discovered CoRoT RRab star, 488 .

| ID     | freq.           | ampl.   | phase [rad]   |
|--------|-----------------|---------|---------------|
|        | [d GLYPH<0> 1 ] | [mag]   |               |
| f 0    | 2.060567        | 0.02902 | 5.282         |
| 2 f 0  | 4.121135        | 0.01676 | 0.301         |
| 3 f 0  | 6.181702        | 0.01121 | 1.921         |
| 4 f 0  | 8.242270        | 0.00783 | 3.485         |
| 5 f 0  | 10.302837       | 0.00521 | 5.147         |
| 6 f 0  | 12.363404       | 0.00340 | 0.345         |
| 7 f 0  | 14.423972       | 0.00258 | 1.832         |
| 8 f 0  | 16.484539       | 0.00195 | 3.453         |
| 9 f 0  | 18.545106       | 0.00132 | 5.020         |
| 10 f 0 | 20.605674       | 0.00082 | 0.236         |

## Appendix A: Frequency content of the new CoRoT RR Lyrae star: 0101315 488

Table A.1 contains the frequencies that we found during the frequency analysis. The star pulsates with a period of 0.485299 days. The following epoch for maxima was found:

$$2 4 5 4 2 3 6. 8 7 0 9 \, H J D + 0 \cdot 4 8 5 3 0 3 3 ( 1 5 ) \cdot E & & ( A. 1 ) & & \text{Des}$$

Digits in parentheses denote the uncertainties. Nine harmonics can be found in the frequency spectrum. No additional frequencies were found.

## Appendix B: Frequency table of the new CoRoT RRab star: 0103800 818

Table B.1 enumerates the frequencies of the non-modulated RR Lyrae, 818 . The period of this star is 0.4659348 days. No modulation was found in this RRab star. We found the following epoch for maxima:

$$\dot { \stackrel {. } { 2 4 5 5 0 2 9. 3 0 4 9 } H J D + 0 } { 0 } { 1 } { 4 6 5 9 3 4 8 } ( 6 ) \cdot E & & ( B. 1 ) & \text{ et al. ( 2 0 } }$$

After pre-whitening, 'forests' of peaks remain around the harmonics. These do not show any obvious modulation pattern, so we decided not to list them in Table B.1.

## Appendix C: Frequency table of the new CoRoT RR Lyrae star: 0104315 804

In Table C.1 we give the frequency content of 804 . The star pulsates with a period of 0.7218221 days. No signs of modulation or additional frequencies were found in this object. The following epoch for maxima was obtained:

$$2 4 5 5 0 1 9. 8 5 2 4 \, H \text{JD} + 0 ^ { d } 7 2 1 8 2 2 1 ( 3 6 ) \cdot E & & ( C. 1 ) & \text{CJD} \left [ 3 0 &$$

## Appendix D: Frequencies of 0100881648

In Table. D.1 we present the result of the Fourier analysis of the blended CoRoT Blazhko RRab star, 648 . The following epoch for maxima was found:

$$2 4 5 4 2 3 6. 9 4 7 2 \, H J D + 0 ^ { d }. 6 0 7 1 8 6 3 ( 4 8 ) \cdot E & & ( D. 1 ) & \underset { \text{ denote } t | } { \pmod { \cdot } }$$

This is a heavily blended Blazhko RRab star. The spectrum shows the frequency corresponding to the fundamental mode pulsation ( f 0 ) and nine harmonics. In addition, the right components of the modulation triplets are seen prominently around the harmonics after pre-whitening. In a few cases the left side peaks are also detected. The remaining spectrum consists of the known orbital frequencies of CoRoT, their linear combinations with the sideral day, many peaks due to low-frequency variations in the frequency interval 0 1 : GLYPH&lt;0&gt; 0 7 d : GLYPH&lt;0&gt; 1 , some residuals around the main frequency and around the low-order harmonics. The remaining grass is at the level of 73 GLYPH&lt;22&gt; mag. The amplitude and frequency variation of the star due to the Blazhko-modulation was already presented in Szab´ et al. (2009). o

## Appendix E: Frequencies of the blended Blazhko RRab star 0101370544

Table. E.1 lists the result of the Fourier analysis of the second blended CoRoT Blazhko RRab star, 544 . This object was observed in color mode of CoRoT, but we chose to present the co-added (white) light frequencies, because this data set is superior compared to the individual color observations. The following epoch for maxima was found:

$$2 4 5 4 2 3 7. 1 0 1 0 \, H J D + 0 ^ { d }. 6 0 5 0 8 7 0 ( 4 3 ) \cdot E & & ( E. 1 ) & & 5. 8 2 4 8 4 &$$

Despite the heavy blending, the spectrum shows the frequency corresponding to the fundamental mode pulsation ( f 0 ), and ten harmonics. In addition, the triplet components of the modulation found around f 0 are revealed around most of the harmonics after pre-whitening. The right side-lobes ( k GLYPH&lt;3&gt; f 0 + f m ) are present with higher amplitudes than their left counterparts. The remaining spectrum consists of the known orbital frequencies of CoRoT, and their linear combinations with the sideral day, many peaks due to low-frequency variations in the frequency interval 0 2 : GLYPH&lt;0&gt; 0 8 d : GLYPH&lt;0&gt; 1 , some residuals around the main frequency and low-order harmonics, and a few remaining peaks close to or below the significance level between f 0 and 2 f 0 . The remaining grass is at the level of 28 GLYPH&lt;22&gt; mag. The amplitude and frequency variation due to the Blazhko-modulation was presented in Szab´ o et al. (2009).

## Appendix F: Frequencies of the CoRoT RRc star 0105036241

Besides the dominant first overtone frequency ( f 1 = 2 68153 d : GLYPH&lt;0&gt; 1 ) and its harmonics, we see f 0 in the spectrum with a characteristic 0.613 frequency ratio with the first overtone radial pulsation (Table F.1). In addition, a bunch of frequencies were found in the [0.5;1.5] d GLYPH&lt;0&gt; 1 frequency range. Upon inspection of the data we found that their origin can be traced back to two remaining discontinuities. Therefore these portions of the data set CJD [3046.0-3048.0] and [3156.5-3157.5] were removed. In the following we analyze the remaining data set.

We found high left side peaks around the main frequency and its harmonics. If we suppose that f m = 0 00585 d : GLYPH&lt;0&gt; 1 is a modulation frequency, the period of the modulation would be longer than the data set. According to that this star may show a long-period Blazhko-modulation, but more data would be necessary to confirm this finding. Another set of modulation-like frequency di GLYPH&lt;11&gt; erence appears in the data set (Fig. 13, Table F.1). We denote the corresponding frequency f b . Neither f m , nor f b can be found in the frequency spectrum, they appear only through combination frequencies. Even combination frequencies involving both f m and f b can be identified. We note here that as we demonstrated in Sec. 4.3 f 0 has temporal amplitude variation, and this is the most probable culprit to cause the appearance of the f b modulation frequencies. The strongest argument favoring this scenario is that f b appears only close to and in combination with f 0 , and is not seen around the main pulsation frequency, f 1 . In addition, frequencies associated with the orbital period of the satellite and its daily aliases are seen at f = 13 967924 14 974027 12 969585 d : ; : ; : GLYPH&lt;0&gt; 1 as usual in CoRoT data. We omitted these peaks from Table F.1. The following epoch for maxima was found:

$$\hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat{ \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ s a r c a y } \hat { \ t h e F. 1 }$$

or taking into account a gradual period change:

$$2 4 5 4 5 7 2. 6 3 0 0 ( 7 ) H J D + 0 ^ { d } 3 7 2 9 2 1 4 ( 2 ) \cdot E - 1 ^ { d } 4 ( 7 ) \cdot 1 0 ^ { - 7 } \cdot E ^ { 2 } ( F. 2 )$$

## Appendix G: Frequencies of the CoRoT RRc star 0105735652

We detect the main frequency f 1 = 3 58218 d : GLYPH&lt;0&gt; 1 and its harmonics, but also many other frequencies with lower amplitude. Among them we found a highly significant frequency at f 0 = 5 82484 d : GLYPH&lt;0&gt; 1 with several peaks around it, then additional peaks

around f 1 . The frequencies are available in Table G.1. The following epoch for maxima and period were found:

$$2 4 5 4 5 7 2. 7 3 2 3 \, H J D + 0 ^ { d } 2 7 9 1 5 9 6 ( 3 8 ) \cdot E \quad \quad ( G. 1 ) \quad \text{ics.}$$

After pre-whitening with the frequencies enumerated in Table G.1, a dense forest of frequencies remains around f 1 . We also see similar residual power around the second and third harmonics, frequencies around 5 9 9 4 16 5 20 1 d : ; : ; : ; : GLYPH&lt;0&gt; 1 , and frequencies related to the orbital frequency of CoRoT.

The large number of side frequencies seen around frequencies f 0 and 9 470874 d : GLYPH&lt;0&gt; 1 may be the result of their nonstationary nature (modulation). We gave an example in Szab´ o et al. (2010) where the frequency forest found around the the half-integer frequencies was modeled and explained by the varying amplitude of these frequencies. We see a very similar situation here.

Table B.1. Frequency table of the non-modulated CoRoT RRab star, 818 . We refrain from listing the CoRoT orbital frequencies and a large number of residual frequencies around the harmonics.

| ID            | freq.                 | ampl.           | phase       |
|---------------|-----------------------|-----------------|-------------|
|               | [d GLYPH<0> 1 ]       | [mag]           | [rad]       |
| f 0           | 2.146223              | 0.34509         | 3.195       |
| 2 f 0         | 4.292446              | 0.17020         | 2.449       |
| 3 f 0         | 6.438669              | 0.12905         | 1.966       |
| 4 f 0         | 8.584893              | 0.08415         | 1.563       |
| 5 f 0         | 10.731116             | 0.06027         | 1.171       |
| 6 f 0         | 12.877339             | 0.03956         | 0.803       |
| 7 f 0         | 15.023562             | 0.02415         | 0.333       |
| 8 f 0         | 17.169785             | 0.01659         | 6.033       |
| 9 f 0         | 19.316008             | 0.01289         | 5.518       |
| 10 f 0        | 21.462232             | 0.00980         | 5.128       |
| 11 f 0        | 23.608468             | 0.00639         | 4.784       |
| 12 f 0        | 25.754666             | 0.00353         | 4.380       |
| 13 f 0        | 27.900991             | 0.00184         | 3.722       |
| 14 f 0        | 30.047044             | 0.00102         | 3.146       |
| 15 f 0        | 32.193767             | 0.00030         | 2.846       |
| 16 f 0        | 34.339850             | 0.00051         | 4.970       |
| 17 f 0        | 36.485549             | 0.00110         | 4.983       |
| 18 f 0        | 38.631887             | 0.00139         | 4.690       |
| 19 f 0        | 40.777791             | 0.00147         | 4.351       |
| 20 f 0        | 42.924600             | 0.00152         | 3.732       |
| 21 f 0        | 45.070642             | 0.00158         | 3.280       |
| 22 f 0        | 47.216924             | 0.00163         | 2.834       |
| 23 f 0        | 49.363091             | 0.00161         | 2.429       |
| 24 f 0        | 51.509352             | 0.00156         | 1.989       |
| 25 f 0        | 53.655410             | 0.00146         | 1.558       |
| 26 f 0        | 55.801814             | 0.00139         | 1.085       |
| 27 f 0        | 57.948115             | 0.00132         | 0.649       |
| 28 f 0 29 f   | 60.094089             | 0.00120         | 0.284       |
| 0             | 62.240271             | 0.00110         | 6.150       |
| 30 f 0        | 64.386660             | 0.00103         | 5.672       |
| 31 f 0        | 66.532893             | 0.00093         | 5.188       |
| 32 f 0 33 f 0 | 68.679352 70.825334   | 0.00086 0.00075 | 4.709 4.368 |
| 34 f 0        | 72.971610             |                 |             |
|               |                       | 0.00072         | 3.945       |
| 35 f 0        | 75.117668             | 0.00063         | 3.607       |
| 36 f 0 37 f   | 77.264339 79.410227   | 0.00056 0.00052 | 2.980 2.709 |
| 0             |                       | 0.00048         | 2.243       |
| 38 f 0        | 81.556746             |                 |             |
| 39 f 0        | 83.703001             | 0.00041         | 1.847       |
| 40 f 0 41 f 0 | 85.849075 87.995504   | 0.00039 0.00034 | 1.373 0.962 |
| f             |                       |                 | 0.570       |
| 42 0 43 f     | 90.141755             | 0.00031         | 0.102       |
| 0 44 f 0      | 92.287719 94.433877   | 0.00030 0.00027 | 5.992       |
| 45 f 0        | 96.580393             | 0.00025         | 5.741       |
| 46 f          | 98.727051             | 0.00022         | 4.954       |
| 0             |                       |                 | 4.753       |
| 47 f 0 48 f 0 | 100.872457 103.020019 | 0.00020         | 3.993       |
| 49 f 0        | 105.164947            | 0.00016 0.00015 | 4.051       |
|               |                       |                 | 3.692       |
| 50 f 0        | 107.311097            | 0.00015         | 3.666       |
| 51 f 0        | 109.456036            | 0.00013         |             |
| 52 f 0        | 111.604225            | 0.00011         | 2.343 2.633 |
| 53 f 0        | 113.749783            | 0.00008         | 1.962       |
| 54 f 0 55 f   | 115.896231            | 0.00011 0.00010 | 1.920       |
| 0 56 f 0      | 118.041414 120.189737 | 0.00010         | 0.893       |
| 57 f 0        | 122.334629            | 0.00007         | 1.414       |

Table C.1. Frequency table of the non-modulated CoRoT RRab star, 804 . We list the dominant fundamental mode pulsation frequency and its harmonics only.

| ID     | freq.           | ampl. [mag]   | phase [rad]   |
|--------|-----------------|---------------|---------------|
|        | [d GLYPH<0> 1 ] |               |               |
| f 0    | 1.385383        | 0.12961       | 2.944         |
| 2 f 0  | 2.770767        | 0.05177       | 2.608         |
| 3 f 0  | 4.156150        | 0.02263       | 2.573         |
| 4 f 0  | 5.541534        | 0.00878       | 2.853         |
| 5 f 0  | 6.926917        | 0.00592       | 3.223         |
| 6 f 0  | 8.312300        | 0.00524       | 3.217         |
| 7 f 0  | 9.697684        | 0.00423       | 3.031         |
| 8 f 0  | 11.083067       | 0.00330       | 2.748         |
| 9 f 0  | 12.468451       | 0.00242       | 2.489         |
| 10 f 0 | 13.853834       | 0.00166       | 2.207         |
| 11 f 0 | 15.239217       | 0.00139       | 1.930         |
| 12 f 0 | 16.624601       | 0.00105       | 1.759         |
| 13 f 0 | 18.009984       | 0.00070       | 1.513         |
| 14 f 0 | 19.395368       | 0.00055       | 1.244         |
| 15 f 0 | 20.780751       | 0.00042       | 0.810         |
| 16 f 0 | 22.166134       | 0.00030       | 0.915         |
| 17 f 0 | 23.551518       | 0.00033       | 0.438         |
| 18 f 0 | 24.936901       | 0.00019       | 0.808         |
| 19 f 0 | 26.322285       | 0.00021       | 6.203         |

Table D.1. Frequency table of the blended Blazhko RRab star, 648 .

| ID                 | freq.           | ampl.   | phase   |
|--------------------|-----------------|---------|---------|
|                    | [d GLYPH<0> 1 ] | [mag]   | [rad]   |
| f 0                | 1.646941        | 0.03280 | 0.911   |
| 2 f 0              | 3.293434        | 0.01208 | 6.160   |
| 3 f 0              | 4.940375        | 0.00504 | 3.919   |
| 4 f 0              | 6.586869        | 0.00163 | 3.594   |
| 5 f 0              | 8.233809        | 0.00094 | 2.063   |
| 6 f 0              | 9.880303        | 0.00084 | 1.410   |
| 7 f 0              | 11.527244       | 0.00067 | 5.489   |
| 8 f 0              | 13.173737       | 0.00048 | 4.474   |
| 9 f 0              | 14.821125       | 0.00032 | 0.821   |
| 10 f 0             | 16.467619       | 0.00020 | 6.190   |
| f 0 GLYPH<0> f m   | 1.633088        | 0.00065 | 0.490   |
| f 0 + f m          | 1.663030        | 0.00437 | 0.759   |
| 2 f 0 GLYPH<0> f m | 3.275110        | 0.00033 | 3.656   |
| 2 f 0 + f m        | 3.309971        | 0.00215 | 3.767   |
| 3 f 0 + f m        | 4.956911        | 0.00168 | 1.290   |
| 4 f 0 + f m        | 6.603405        | 0.00062 | 0.828   |
| 5 f 0 + f m        | 8.250346        | 0.00021 | 5.675   |
| 6 f 0 + f m        | 9.898180        | 0.00029 | 0.681   |
| 7 f 0 + f m        | 11.543333       | 0.00028 | 4.654   |
| 8 f 0 + f m        | 13.190274       | 0.00027 | 2.231   |
| 9 f 0 + f m        | 14.836320       | 0.00019 | 2.538   |
| 10 f 0 + f m       | 16.482814       | 0.00014 | 1.552   |

Table E.1. Frequency table of the white band light curve of the blended Blazhko RRab star, 544 .

| ID                  | freq.           | ampl.   | phase   |
|---------------------|-----------------|---------|---------|
|                     | [d GLYPH<0> 1 ] | [mag]   | [rad]   |
| f 0                 | 1.652655        | 0.00930 | 5.137   |
| 2 f 0               | 3.305310        | 0.00323 | 0.437   |
| 3 f 0               | 4.957965        | 0.00162 | 2.316   |
| 4 f 0               | 6.610620        | 0.00051 | 4.400   |
| 5 f 0               | 8.263275        | 0.00016 | 0.517   |
| 6 f 0               | 9.915930        | 0.00017 | 3.025   |
| 7 f 0               | 11.568585       | 0.00013 | 5.008   |
| 8 f 0               | 13.221240       | 0.00013 | 0.456   |
| 9 f 0               | 14.873895       | 0.00010 | 2.244   |
| 10 f 0              | 16.526550       | 0.00007 | 4.090   |
| 11 f 0              | 18.179205       | 0.00007 | 5.635   |
| f 0 + f m           | 1.691938        | 0.00163 | 0.917   |
| 2 f 0 + f m         | 3.344600        | 0.00083 | 1.996   |
| 3 f 0 + f m         | 4.997249        | 0.00067 | 3.691   |
| 4 f 0 + f m         | 6.649602        | 0.00033 | 5.857   |
| 5 f 0 + f m         | 8.302697        | 0.00014 | 1.593   |
| 6 f 0 + f m         | 9.955316        | 0.00009 | 4.047   |
| 7 f 0 + f m         | 11.608022       | 0.00013 | 6.247   |
| 8 f 0 + f m         | 13.260818       | 0.00011 | 1.749   |
| 9 f 0 + f m         | 14.913797       | 0.00010 | 3.478   |
| 10 f 0 + f m        | 16.566175       | 0.00008 | 5.297   |
| 11 f 0 + f m        | 18.219099       | 0.00006 | 0.945   |
| f 0 GLYPH<0> f m    | 1.613977        | 0.00028 | 2.873   |
| 2 f 0 GLYPH<0> f m  | 3.266438        | 0.00015 | 4.474   |
| 3 f 0 GLYPH<0> f m  | 4.919052        | 0.00015 | 5.969   |
| 4 f 0 GLYPH<0> f m  | 6.571116        | 0.00013 | 1.999   |
| 5 f 0 GLYPH<0> f m  | 8.224252        | 0.00007 | 3.496   |
| 6 f 0 GLYPH<0> f m  | 9.877031        | 0.00007 | 5.558   |
| 7 f 0 GLYPH<0> f m  | 11.529681       | 0.00004 | 0.594   |
| 8 f 0 GLYPH<0> f m  | 13.182138       | 0.00004 | 3.375   |
| 10 f 0 GLYPH<0> f m | 16.487827       | 0.00004 | 0.997   |
| f 0                 | 2.389287        | 0.00015 | 4.079   |

Table F.1. Frequency table of the CoRoT RRc star, 241 .

| ID                            | freq.           | ampl.   | phase       |
|-------------------------------|-----------------|---------|-------------|
|                               | [d GLYPH<0> 1 ] | [mag]   | [rad]       |
| f 1                           | 2.68153         | 0.19610 | 0.337       |
| 2 f 1                         | 5.36271         | 0.01186 | 4.704       |
| 3 f 1                         | 8.04458         | 0.01409 | 2.349       |
| 4 f 1                         | 10.72611        | 0.01120 | 5.662       |
| 5 f 1                         | 13.40729        | 0.00728 | 2.976       |
| 6 f 1                         | 16.08882        | 0.00454 | 6.118       |
| 7 f 1                         | 18.77035        | 0.00288 | 2.850       |
| 8 f 1                         | 21.45187        | 0.00177 | 5.771       |
| 9 f 1                         | 24.13375        | 0.00105 | 1.973       |
| 10 f 1                        | 26.81493        | 0.00065 | 5.299       |
| 11 f 1                        | 29.49611        | 0.00044 | 2.292       |
| 12 f 1                        | 32.17764        | 0.00027 | 5.439       |
| f 0                           | 4.37783         | 0.00332 | 2.953       |
| f 0 GLYPH<0> f 1              | 1.69630         | 0.00052 | 3.322       |
| f 0 + f 1                     | 7.05901         | 0.00304 | 5.886       |
| 2 f 1 + f 0                   | 9.74054         | 0.00147 | 1.980       |
| 2 f 1 GLYPH<0> f 0            | 0.98488         | 0.00039 | 5.187       |
| 3 f 1 + f 0                   | 12.42241        | 0.00035 | 4.613       |
| f 1 GLYPH<0> f m              | 2.67601         | 0.00415 | 0.369       |
| f 1 + f m                     | 2.69084         | 0.00151 | 5.743       |
| f 1 GLYPH<0> 2 f m            | 2.67015         | 0.00147 | 4.535       |
| f 1 + 2 f m                   | 2.69946         | 0.00055 | 5.690       |
| f 1 GLYPH<0> 3 f m            | 2.66256         | 0.00026 | 3.728       |
| 2 f 1 + f m                   | 5.35823         | 0.00060 | 2.256       |
| 2 f 1 GLYPH<0> f m            | 5.36823         | 0.00036 | 4.349 1.887 |
| 3 f 1 + f m                   | 8.03906         | 0.00054 |             |
| 3 f 1 GLYPH<0> f m            | 8.05286         | 0.00025 | 2.841       |
| 4 f 1 + f m                   | 10.72025        | 0.00063 | 5.756       |
| 4 f 1 GLYPH<0> f m            | 10.73232        | 0.00025 | 1.842       |
| 5 f 1 GLYPH<0> f m            | 13.40143        | 0.00043 | 2.244       |
| 6 f 1 GLYPH<0> f m            | 16.08296        | 0.00029 | 5.635       |
| 7 f 1 GLYPH<0> f m            | 18.76414        | 0.00026 | 2.625       |
| 8 f 1 GLYPH<0> f m            | 20.44610        | 0.00032 | 5.533       |
| f 0 GLYPH<0> f m              | 4.37266         | 0.00035 | 1.050       |
| f 0 + f m                     | 4.38472         | 0.00186 | 2.572       |
| f 1 GLYPH<0> f b              | 2.64290         | 0.00044 | 1.131       |
| f 1 GLYPH<0> 2 f b            | 2.60014         | 0.00038 | 3.784       |
| 2 f 1 GLYPH<0> f b            | 9.69881         | 0.00039 | 4.512       |
| f 0 GLYPH<0> f b              | 4.33852         | 0.00172 | 1.475       |
| f 0 + f b                     | 4.41852         | 0.00026 | 2.532       |
| f 0 GLYPH<0> 2 f b            | 4.30265         | 0.00074 | 5.346       |
| f 0 GLYPH<0> 3 f b            | 4.25989         | 0.00024 | 5.469       |
| f 0 + f 1 + f m               | 7.06625         | 0.00159 | 5.015       |
| f 0 + f 1 GLYPH<0> f m        | 7.05418         | 0.00062 | 3.355       |
| 2 f 1 + f 0 + f m             | 9.74743         | 0.00084 | 1.622       |
| f 0 + f 1 GLYPH<0> 2 f b      | 6.98383         | 0.00036 | 2.248       |
| f 0 + f 1 + 2 f b             | 7.14379         | 0.00043 | 3.320       |
| f 0 GLYPH<0> f b GLYPH<0> f m | 4.33369         | 0.00069 | 2.981       |
| f 0 + f b GLYPH<0> f m        | 4.34472         | 0.00043 | 4.744       |
| f 0 + f b + f m               | 4.42611         | 0.00053 | 4.089       |
| f 0 + 2 f b + f m             | 4.46197         | 0.00098 | 5.216       |

Table G.1. Frequency table of the CoRoT RRc star, 652 .

| ID                 | freq.           | ampl.   | phase   |
|--------------------|-----------------|---------|---------|
|                    | [d GLYPH<0> 1 ] | [mag]   | [rad]   |
| f 1                | 3.58218         | 0.20403 | 0.548   |
| 2 f 1              | 7.16435         | 0.03209 | 4.309   |
| 3 f 1              | 10.74653        | 0.01666 | 1.496   |
| 4 f 1              | 14.32871        | 0.01356 | 5.530   |
| 5 f 1              | 17.91088        | 0.00992 | 2.685   |
| 6 f 1              | 21.49306        | 0.00671 | 5.825   |
| 7 f 1              | 25.07541        | 0.00420 | 2.693   |
| 8 f 1              | 28.65754        | 0.00240 | 5.670   |
| 9 f 1              | 32.23960        | 0.00145 | 2.291   |
| 10 f 1             | 35.82177        | 0.00074 | 5.136   |
| 11 f 1             | 39.40394        | 0.00054 | 1.507   |
| 12 f 1             | 42.98612        | 0.00039 | 4.327   |
| 13 f 1             | 46.56830        | 0.00033 | 0.575   |
| 14 f 1             | 50.15047        | 0.00023 | 0.712   |
| f 0                | 5.82484         | 0.00216 | 5.731   |
| f 0 + f m          | 5.88925         | 0.00144 | 5.090   |
| f 0 GLYPH<0> f m   | 5.76587         | 0.00084 | 3.699   |
| f 0 + 2 f m        | 6.00000         | 0.00073 | 2.904   |
| f 0 GLYPH<0> 2 f m | 5.71282         | 0.00035 | 1.072   |
|                    | 5.93520         | 0.00043 | 6.179   |
|                    | 5.87227         | 0.00054 | 2.259   |
|                    | 5.77659         | 0.00075 | 0.001   |
|                    | 9.47087         | 0.00089 | 2.479   |
|                    | 9.40747         | 0.00064 | 0.961   |
|                    | 9.58200         | 0.00048 | 5.982   |
|                    | 9.34758         | 0.00037 | 0.292   |
|                    | 9.29643         | 0.00023 | 3.506   |
| f 1 + f m          | 9.62932         | 0.00017 | 1.606   |
|                    | 3.64669         | 0.00087 | 2.133   |
|                    | 2.66379         | 0.00064 | 3.466   |
R. Szabó J. M. Benkő M. Paparó E. Chapellier E. Poretti A. Baglin W. W. Weiss K. Kolenberg E. Guggenberger J. -F. Le Borgne
Revisiting CoRoT RR Lyrae stars: detection of period doubling and temporal variation of additional frequencies
We search for signs of period doubling in CoRoT RR Lyrae stars. The
occurrence of this dynamical effect in modulated RR Lyrae stars might help us
to gain more information about the mysterious Blazhko effect. The temporal
variability of the additional frequencies in representatives of all subtypes of
RR Lyrae stars is also investigated. We pre-process CoRoT light curves by
applying trend and jump correction and outlier removal. Standard Fourier
technique is used to analyze the frequency content of our targets and follow
the time dependent phenomena. The most comprehensive collection of CoRoT RR
Lyrae stars, including new discoveries is presented and analyzed. We found
alternating maxima and in some cases half-integer frequencies in four CoRoT
Blazhko RR Lyrae stars, as clear signs of the presence of period doubling. This
reinforces that period doubling is an important ingredient to understand the
Blazhko effect - a premise we derived previously from the Kepler RR Lyrae
sample. As expected, period doubling is detectable only for short time
intervals in most modulated RRab stars. Our results show that the temporal
variability of the additional frequencies in all RR Lyrae sub-types is
ubiquitous. The ephemeral nature and the highly variable amplitude of these
variations suggest a complex underlying dynamics of and an intricate interplay
between radial and possibly nonradial modes in RR Lyrae stars. The omnipresence
of additional modes in all types of RR Lyrae - except in non-modulated RRab
stars - implies that asteroseismology of these objects should be feasible in
the near future (Abridged).