The Blazhko behaviour of RR~Geminorum II

The Blazhko behaviour of RR~Geminorum Ii

long-term photometric results
Á. Sódor Konkoly Observatory of the Hungarian Academy of Sciences. P.O. Box 67, H-1525 Budapest, Hungary
11email: sodor,szeidl,jurcsik@konkoly.hu
   B. Szeidl Konkoly Observatory of the Hungarian Academy of Sciences. P.O. Box 67, H-1525 Budapest, Hungary
11email: sodor,szeidl,jurcsik@konkoly.hu
   Konkoly Observatory of the Hungarian Academy of Sciences. P.O. Box 67, H-1525 Budapest, Hungary
11email: sodor,szeidl,jurcsik@konkoly.hu
   J. Jurcsik Konkoly Observatory of the Hungarian Academy of Sciences. P.O. Box 67, H-1525 Budapest, Hungary
11email: sodor,szeidl,jurcsik@konkoly.hu
Received ; accepted
Key Words.:
Stars: individual: RR~Gem – Stars: variables: RR Lyr – Stars: oscillations – Stars: horizontal-branch – Techniques: photometric –
thanks: Tables are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/
Abstract

Context:RR~Gem is one of the few Blazhko RR Lyrae that has photometric observations available extended enough to study the long-term courses of its pulsation and modulation properties in detail.

Aims:We investigate the pulsation and modulation properties and the relations between them in RR~Gem using photometric observations from the past 70 years in order to gain further insight into the nature of the Blazhko modulation.

Methods:We studied the photographic, photoelectric, and CCD light curves obtained at the Konkoly Observatory and other authors’ published maxima observations. Detailed analysis of the light curves, maximum brightness, and data are carried out.

Results:RR~Gem showed modulation most of the time it was observed. The modulation amplitude showed strong variations from the undetectable level (less than 0.04 mag in maximum brightness) to about 0.20 mag. The amplitudes of the amplitude and phase modulations showed parallel changes, thus the total ‘power’ of the modulation have changed during the past 70 years. Parallel changes in the pulsation and modulation periods occur with a ratio. We also detected  mag changes in the mean maximum brightness and mean pulsation amplitude.

Conclusions:

1 Introduction

There is no satisfactory and commonly accepted explanation yet for the phenomenon of the light curve modulation of RR Lyrae stars, the century’s old Blazhko effect (Blazhko 1907). Less than a dozen Blazhko RR Lyrae have photometric observations available that are extended enough to study their long-term modulation properties in detail. RR~Gem (, , J2000) is one of these stars that has been observed at the Konkoly Observatory from time to time since 1935. The observations were always obtained with the most up-to-date photometric technique available at the time: between 1935 and 1953 photographic, from 1954 to 1983 photoelectric, and in 2004 and 2005 CCD observations. (The photographic, photoelectric, and the 2005 CCD observations have not been published yet.)

RR~Gem was claimed to be Blazhko-modulated with a modulation period of 48 days based on the photographic observations taken in the 1930s at the Konkoly Observatory (Balázs 1960). Although the modulation period was dubious, the fact of modulation itself seemed to be unambiguous. In contrast, in the 1950s the modulation seemed to cease (Detre 1970). We continued to observe RR~Gem in 2004 and 2005 and found it to be modulated but with very low amplitude and with the shortest period (7.216 d) known at that time (Jurcsik et al. 2005a; Sódor 2006).

In the first part of this series of papers, the results of the 2004 CCD observations of RR~Gem were reported (Jurcsik et al. 2005a, Paper1 hereafter). This second part reports our results of the reexamination of all the available photometric data of RR~Gem in order to reveal the long-term behaviour of the modulation.

2 Data

We utilise photometric observations made by different instruments of the Konkoly Observatory. Most of these data are previously unpublished archive observations. We also use maxima timing observations of the past 114 years from the literature.

begin end telescope  detector filters time res. data observed comp.
[JD] date [JD] date [min]  points nights/maxima star 
2 427 833 1935 Jan 2 434 457 1953 Mar 16 cm pg unf. 668 37/19
2 435 052 1954 Nov 2 436 227 1958 Jan 60 cm pe unf. 604 22/19 A
2 436 229 1958 Jan 2 443 931 1979 Feb 60 cm pe 1 754 37/28 A
2 441 679 1972 Dec 2 445 673 1983 Dec 50 cm pe 440 4/4 A
2 453 019 2004 Jan 2 453 440 2005 Mar 60 cm CCD 13 355 63/31 B
See details about the telescopes in Sections .
Time resolution means the typical cycle times in the case of multicolour observations.
Comparison stars used are: A - GSC 02452-01847; B - BD~$+31°1547$22footnotemark: 2.
Table 1: Log of the Konkoly observations of RR~Gem.

2.1 Photometric data

The light curves we use were obtained at the Konkoly Observatory between 1935 and 2005. Only a limited number of photoelectric data have been published by other authors (Fitch et al. 1966; Epstein 1969; Stepien 1972; Liu & Janes 1989; ESA 1997). These data are, however, few in number and too sparse for studying the modulation behaviour.

The log of the Konkoly observations is found in Table 1. The photographic, photoelectric, and 2005 CCD data, listed in Tables , are only available electronically at the CDS (http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/). In each table, column 1 lists the HJD of the observations. In Table 4, column 2 gives the photographic magnitudes, while in Tables  differential magnitudes of RR~Gem are given with respect to GSC 02452-01847 and BD~$+31°1547$22footnotemark: 2 for the photoelectric and CCD data, respectively.

2.1.1 Photographic data

The photographic observations were made with the 16 cm astrograph of the Konkoly Observatory at Budapest, Svábhegy. Between 1935 and 1953, about 1000 exposures were taken on about 80 photographic plates. Unfortunately, many of the plates were lost or damaged during the past 70 years. We succeeded in recovering 598 measurable exposures on 56 plates.

The photographic plates were digitized on a Umax PowerLook 3000 flatbed transparency scanner with a spatial resolution corresponding to about 1.8 ”/pixel. Digital aperture photometry was applied to the images using standard IRAF111IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. packages. In this way, we determined the photographic blackenings of RR~Gem and of the photographic comparison stars.

A series of comparison stars were chosen from the 20 arcmin vicinity of RR~Gem to cover its brightness variation range. Altogether, 15 appropriate comparison stars were found. We measured the instrumental magnitudes of the comparison stars during the course of the 2005 CCD observations. The differential magnitudes were determined with respect to the magnitude of BD~$+31°1547$222This star was used as comparison for the CCD observations. In the first paper of this series, the wrong BD number was erroneously given (BD ) for this star. () we had derived earlier (see Paper1). The magnitudes of RR~Gem were determined by a order polynomial blackening curve fit for the comparison stars of each exposure. The most deviating points were omitted from the fits.

To check the consistency between our instrumental CCD photometry and the photographic measurements, we also automatically computed the magnitudes of the comparison stars for each of the exposures by the aforementioned blackening curve-fitting process using all the other comparison stars except for the most deviating ones. The magnitudes determined in this way showed good agreement with the accepted instrumental CCD magnitudes for all the 15 comparison stars. The differences have a mean of 0.00 mag with a standard deviation of 0.05 mag. Neither colour nor brightness dependency of the residuals was found according to our CCD measurements.

A record of the evaluation of only a part of the photographic observations has remained available. It contains 331 photometric measurements of RR~Gem, 70 of them from damaged or lost plates. These 70 points were brought onto our magnitude scale with a transformation that was derived from a linear regression of the 261 common data points. The photographic light curve was complemented with these 70 measurements.

2.1.2 Photoelectric data

The photoelectric observations were made with different unrefrigerated photometers between 1954 and 1983. The photometers in the Newton focus of the 60 cm Heyde telescope at Budapest, Svábhegy, first employed an RCA 1P21 and later an EMI 9052 B photomultiplier tube. In the Cassegrain focus of the 50 cm telescope at the Piszkéstető Mountain Station, a photometer was used with an EMI 9058 QB multiplier. All the instruments were equipped with conventional Schott filters ( - UG2 2mm, - BG12 1mm + GG13 2mm, - GG11 2mm), except for the first observations, which were made without filter. A standard photoelectric reduction procedure was applied. No correction was applied to the differential extinction, as the comparison star is as near as 8’ to the variable. The colour extinction was also neglected because the colour of the variable and the comparison match at maximum brightness. The filtered measurements were transformed to the standard Johnson system according to the regularly determined telescope constants from standard star measurements. Measurements in the band were only made on 3 nights in the past few years of the photoelectric observations. These data are also published electronically but not used in the analysis. For the photoelectric observations, the star GSC 02452-01847 was used as comparison. We derived the standard magnitude of of this star from the CCD measurements.

Figure 1: The diagram of RR~Gem over the past 115 years according to the average pulsation period (upper line of points). The lower line of points represents the data of the Konkoly light curves, calculated using the same period. These points are shifted downwards by 0.2 d for better visibility (see the axis on the right side). The intervals of Konkoly observations are indicated.

2.1.3 CCD data

The 2004 CCD observations and results have been published in Paper1. In order to refine the pulsation and modulation periods and to measure magnitudes of stars in a larger area to calibrate the photographic data, further CCD observations were obtained in 2005 (Sódor 2006) with the same instrumentation (the 60 cm Heyde telescope with a 1152 x 770 Wright CCD detector in its Newton focus, which gives a field of view of 24’ x 17’). These data were reduced and transformed to the standard system in the same way as for the 2004 observations (see Paper1 for details).

2.2 The data

The changes in pulsation period and the phase modulation properties are studied using the values of the maximum brightness times and the times of the brightness on the rising branch, referred to as and , respectively. The and timings of the maxima and rising branches were determined from the Konkoly observations. The values were calculated according to the corresponding ephemerides, which are given in the subsequent sections.

2.2.1

Most of the published times of maxima were collected in the GEOS333http://dbrr.ast.obs-mip.fr/ database (Graff 1922; Waterfield 1927; Luizet 1930; Ivanov 1930a, b; Chudovichev 1930; Dubiago 1930; Blazhko 1935; Lange 1935, 1969; Kleissen 1939; Batyrev 1950, 1962; Alania 1956; Guriev 1958; Mandel 1960; Fitch et al. 1966; Epstein 1969; Tsesevich 1969; Kanischeva & Lange 1971; Lange et al. 1976; Braune et al. 1977; Liu & Janes 1989; ESA 1997; Vandenbroere 1997, 1998, 1999, 2001, 2003, 2005; Gensler 1998; Agerer & Huebscher 2000, 2002, 2003; Huebscher 2000, 2003, 2005; Agerer et al. 2001; Huebscher et al. 2005; Le Borgne 2004; Le Borgne et al. 2005, 2006; Wils et al. 2006). A few additional maxima times of Waterfield (1927) and Graff (1922) were also considered. We have left out some outlying visual points where no cause of the deviation was found. Altogether, 289 reliable data points were used. Further 101 points were determined from the Konkoly light curves by order polynomial fits to the data points around the maxima. The times of the maxima were averaged for the , , , and bands in the case of multicolour observations, since the times of maximum brightness do not differ significantly in these bands. Times of maxima in the band were not taken into account because of the larger uncertainty of these data. The fits also yielded the magnitudes of maximum brightnesses, which were utilised in studying the modulation behaviour.

2.2.2

The phase modulation can be studied through the data of a particular point of the rising branch or of the maximum brightness. Because the time of maximum brightness can be determined with less accuracy than the time of a certain magnitude phase on the rising branch, we decided to investigate the modulation of the 11.3 mag brightness phase on the rising branch in the band. This brightness approximately bisects the range from the middle of the rising branch to the mean maximum light.

The variation measures, in fact, a combination of the amplitude and phase modulations. It shows oscillations due to the changing slope of the rising branch when the light curve change can be described as only a scaling in magnitude (exact amplitude modulation), as well as due to exact phase modulation. Just to estimate the contribution of the exact amplitude modulation in the amplitude of the , test data were generated using the mean CCD light curve of RR~Gem and assuming exact amplitude modulation. The modulation properties of the test data are compared with the observations of RR~Gem in Sect. 5.1.

The phase modulation can be studied in the data only if the effect of the long-term period changes are eliminated. Therefore we used different ephemerides for the different parts of the observations, as the period changes require. The ephemerides applied to construct the data for the different parts of the observations are given in the corresponding sections.

The maximum brightness times, magnitudes, and data are given in Tables  for the photographic, photoelectric, and CCD observations, respectively. These tables are only available electronically at the CDS (http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/). In all 3 Tables, columns  give the data, their errors, and the times of the maximum light, respectively. In Table 13, columns 4 and 5 contain the magnitudes and their errors, respectively. In Table 14, columns 4 and 6 list the photoelectric and magnitudes, and columns 5 and 7 give their errors, respectively. In Table 15, columns 4 and 5 list the CCD magnitudes and their errors.

3 The diagram

The variation in the pulsation period during the past 115 years can be followed in the diagrams shown in Fig. 1. The diagrams were constructed using the ephemerides

The period corresponds to the average pulsation period of RR~Gem over the past 115 years. The data calculated with this single ephemeris shows similar changes to the data.

The diagram indicates steady and abrupt pulsation period changes. The period of the pulsation varied between the extrema of 0.397253 d and 0.3973177 d. The transition between these values occurred abruptly just in the middle of the photographic observations, around JD 2 428 800. After the period jump, the period was decreasing during the next 20000 days (see Fig. 1 in Sódor 2006) and seem to remain stable in the past ten years. Details of the period change rates and observed period values are given in the next sections.

Note that the shape of the diagram, in particular the aforementioned sudden period change, strongly resembles that of the Blazhko RR Lyrae stars, XZ~Cyg (Baldwin & Samolyk 2003) and RW~Dra (Balázs & Detre 1962).

4 Light curve analysis

The Konkoly photometric data were analysed in order to detect the modulation if existed, to follow any changes in the pulsation and modulation characteristics, and to find possible connections between them.

The observations were divided into several parts whenever the different observing methods, the data distribution in time, and the changes in the pulsation period allowed. When the phase coverage and data quality made it possible, Fourier analysis was performed to search for modulation components in the spectra. The modulation was also examined through the maximum brightness and the data points calculated using the actual mean pulsation values. In this way, the maximum brightness variation caused by amplitude modulation and the phase modulation of the rising branch were studied.

Data analysis was performed using the different applications of the MUFRAN package (Kolláth 1990) and the linear and nonlinear curve fitting abilities of gnuplot444http://www.gnuplot.info/ and Mathematica555Mathematica is a registered trademark of Wolfram Research Inc. (Wolfram 1996).

4.1 Results from photographic observations

Because of the sudden and remarkable period change that took place at the middle of the photographic observations, this data set was divided into two parts separated at JD 2 428 800 (denoted as pg1 and pg2, respectively). In this way, we can follow whether the physical processes that gave rise to the sudden and significant change in the pulsation period have any effect on the modulation properties.

The pg2 data set has a much less favourable distribution in time than that of the pg1 data. There are 9 light maxima and 9 rising branches from two consecutive seasons in pg1, whereas we have 10 observed maxima and 15 rising branches from the 15-year long pg2 interval.

The pulsation period was determined by Fourier fit to the light curves taking 8 harmonics into account. The pg2 data set does not cover each phase of the pulsation, since the observations were focused on the rising branch and maximum light. The gap in the middle of the descending branch was bypassed with artificial points based on the CCD mean light curve, shifted and scaled appropriately. These points were used only to stabilize the Fourier fit so they had no effect on any of the conclusions. The data show that a slight period decrease occurred during the pg2 interval, therefore a linear period change was taken into account. The pulsation periods derived for the two parts of the photographic data are

where

The light curves of the two intervals folded with the appropriate pulsation periods are plotted in Fig. 2. The fitted mean light curves are also shown in the plots. The plotted pg2 light curve was transformed to constant period by a HJD transformation (see details in Jurcsik et al. 2001).

Figure 2: Konkoly photographic light curves for the pg1 and pg2 intervals. Solid lines show the mean light curves that take 8 harmonics into account. Grey lines show the mean curve of the other data set to help compare the pulsation amplitudes. The times of the plotted pg2 light curve were transformed in order to eliminate the period change.

From the times (listed in Table 13), values were calculated according to the following ephemerides:

for the two intervals, respectively.

The light curves shown in Fig. 2 indicate that the mean pulsation amplitude had increased for the second interval perceptibly. No differences in the plate material and in the data evaluation of the pg1 and pg2 observations account for these discrepancies.

4.1.1 Modulation properties of the Pg1 data (JD )

The rising branch of the pg1 light curve shows greater scatter around the mean curve than that of the pg2 data. Possible explanations of this phenomenon are incessant period change or phase modulation during the pg1 interval. Because the values change sign more than twice during the pg1 interval, these changes cannot be explained with a unidirectional period change. Recently Chadid & Chapellier (2006) have reported an irregular behaviour of the Blazhko modulation of RR~Lyr, and they explain this phenomenon by the shock wave passage across the atmosphere. According to Chadid (2000), there are important irregularities in the atmosphere of RR~Lyr that occur during the rising branch of successive pulsation cycles. This effect can also be responsible for the widening of the rising branch in the pg1 data set of RR~Gem. The irregular behaviour on the rising branch of RR~Lyr pointed out by Chadid & Chapellier (2006) is, however, connected with the modulation. Therefore, as steady period change can be excluded, the larger scatter on the rising branch of RR~Gem indicates the presence of modulation.

The modulation of the pg1 interval is first investigated in the rising branch phase data. These data, calculated from the times (see electronic Table 13) according to the ephemeris in Sect. 4.1, were Fourier-analysed. The plausible modulation frequency range of were examined. The Fourier amplitude spectrum of the pg1 data is plotted in the top panel of Fig. 3. The most significant peak of the spectrum is at 0.139 c/d frequency. Prewhitening the data with this frequency yields the residual spectrum shown in the bottom panel of Fig. 3. This spectrum only shows low residual noise. The best sine curve fit to the data gives , , and . denotes peak-to-peak amplitude throughout this paper. The data phased with this period and the fitted harmonic curve are shown in the top left panel of Fig. 4. The order harmonic curve fit the data with very small scatter, giving convincing evidence of phase modulation.

Due to the large scatter and significant nightly systematic errors, the maximum brightness magnitudes from the pg1 interval show the modulation period with uncertainty. These data cannot be used to refine the modulation period valid for this interval, but by accepting the period of the modulation from the data, we can estimate the extent of the amplitude modulation. The middle left panel of Fig. 4 shows the maximum brightness data folded with the 7.20 d modulation period. A sine curve fit with this period gives and  mag values for the modulation amplitude and the average maximum brightness, respectively. The fact that the maximum brightness values can be phased with the period derived from the data strengthens the finding of modulation during this interval.

Figure 3: Fourier amplitude spectrum of the data of the pg1 observations. The insert shows the spectral window function. The dotted line shows the significance level. An arrow points to the peak at . Bottom panel shows the residual spectrum after prewhitening with this frequency.
Figure 4: The modulation in (top), in maximum brightness (middle) and the relation between them (bottom) for the studied light curves. The and maximum brightness points are folded with the determined modulation period for the pg1, pe1, and CCD data, and are plotted versus HJD for the pg2 and pe2 data. The fitted harmonic curves are also shown. In the second bottom panels, each point represents a particular pulsation cycle for which both and maximum brightness could be determined. These plots are independent of the modulation period and reflect the amplitude relation and phase connection of the amplitude and phase modulation components. The plotted ellipses are the combinations of the fitted harmonic curves. In the bottom right panel we show the relation valid for test data of exact amplitude modulation.

4.1.2 Modulation properties of the Pg2 data (JD )

The unfavourable data distribution makes the Fourier amplitude spectra of the pg2 maximum brightness and data uninterpretable because of the severe aliasing effects and the possible change in the modulation period.

Though the period of the modulation cannot be derived from the pg2 data, we are able to estimate the amplitudes of the phase and amplitude modulations. In the bottom panels of Fig. 4 the measured versus maximum brightness data of the individual pulsation cycles are plotted. These plots are independent of the actual value of the modulation period. Any correlation or connection between these data is the sign that the measured scatters of the and maximum brightness values are due to modulation rather than observational inaccuracy or any irregular behaviour. The bottom panel of the pg2 data in Fig. 4 shows that this is indeed the case, since the variations in and in maximum brightness are not independent. The possible extent of the amplitudes of the modulation can also be read from this plot: and . If we assume that the pg1 and pg2 data have similar noise statistics and that the ratio of the modulation amplitudes and the scatters of the data are identical for the two data sets, then we can estimate the modulation amplitudes of the pg2 data simply from the observed scatter of the maximum brightness and values (see Table 2). The modulation amplitudes estimated in this way are in very good agreement with those read from Fig. 4, confirming that the derived amplitudes of the modulation are reliable.

The average maximum brightness was  mag according to the pg2 data, which was  mag brighter than observed during the pg1 period. Although the scatter around minimum light made the estimation of the mean pulsation amplitude somewhat uncertain, it had increased with about  mag from the pg1 to the pg2 interval.

To check whether the observed change of the mean pulsation amplitude and/or the maximum brightness is real or is only an instrumental effect, we compared the mean photographic magnitudes of each of the comparisons stars separately for the pg1 and pg2 data. The magnitudes show no systematic variations between the two intervals, and the magnitude differences are less than 0.04 mag for each star. Consequently, we regard the observed 0.14 mag increase in the average maximum brightness as real, and it is most probably a sign of real change in the mean pulsation amplitude as well.

4.2 Results from photoelectric observations

The photoelectric observations made without any filter are quite inhomogeneous as they were obtained during the test phase of the new photometer with several different instrument settings. Therefore, this data set is only suitable for determining the data but not for investigating the modulation.

The data of the unfiltered observations, the times, and magnitudes of all the observed and light maxima, as well as the timing data of the photoelectric light curves, are listed in Table 14.

The pulsation period changed with a constant rate between JD 2 435 000 and JD 2 445 000. The change rate has increased with about an order of magnitude in the years after JD 2 445 000. The time distribution of the observations are uneven as nearly the half of the data points were obtained during the first 2 seasons (pe1 data, before JD 2 436 629). Later, the observations were very sparse: nights per annum (pe2 data). The pe2 data contain the observations made between JD 2 536 956 and JD 2 543 931, as the period change rate increased later.

The period of the pulsation during the pe1 and pe2 observations was

where

The period and its change rate were determined by a order polynomial fit to the Konkoly photoelectric data.

The values were calculated from the pe1 and pe2 times (listed in Table 14) according to the following ephemeris:

The pe1 data folded with the average pe1 pulsation period and the fitted mean light curve taking 13 harmonics into account are plotted in Fig. 5.

Figure 5: The pe1 photoelectric observations phased with the pulsation period and the fitted mean light curve.

4.2.1 Modulation properties of the Pe1 data (JD )

The light curve in Fig. 5 shows no sign of phase modulation on the rising branch; however, the scatter in magnitude around the light maximum is higher than the uncertainty of the observations that indicates variation in amplitude, either in a regular or an irregular way. The effect of period change is negligible during the pe1 interval. The Fourier amplitude spectra of the pe1 and maximum brightness data are shown in Fig. 6. Both spectra show two peaks around the expected 0.139 c/d frequency value. The frequency that fits both the and data sets the best is (i.e. ). The distance between this peak and its neighbour is 0.003 c/d. The cycle/year alias frequency results in a fit with significantly higher rms. The pe1 and maximum brightness data phased with are plotted in Fig. 7. The amplitudes of the fitted harmonic curves are and for the and data, respectively. The average maximum brightness in is . According to the fitted mean light curve, the mean pulsation amplitude in band is .

Figure 6: Fourier amplitude spectra of the pe1 and maximum brightness data. Dotted lines show the significance level. Arrows point to the peaks at . Panels marked with ‘pw’ show the residual spectra after prewhitening with this frequency. Bottom panel shows the spectral window function.
Figure 7: Photoelectric and maximum brightnesses of the pe1 data folded with the 7.28 d modulation period.

4.2.2 Modulation properties of the Pe2 data (JD )

The sparse data distribution in time and the pulsation period change during this interval makes these data hard to interpret. Therefore, the pe2 data are treated similarly to those of pg2. Upper limits of the amplitudes of the possible modulation are estimated from the plots shown in Fig. 4 and from the comparison of the parameters of the pe1 and pe2 data listed in Table 2. The estimations were made assuming identical ratios for the two photographic data sets and also for the two photoelectric data sets.

The scatter of the maximum brightness magnitudes of the pe2 data is more than 3 times greater than the average error of the individual points. It suggests that modulation at this time might also exist, but with even slightly smaller amplitude than in the first part of the photoelectric observations. The data show no sign of phase modulation. The average maximum brightness of the pe2 light curve was and the mean pulsation amplitude was the same as for the pe1 data within the error range.

data set
[mag] [mag] [d] [d]
pg1 0.065 0.19 0.0034 0.0094
pg2 0.052 0.15 0.0011 0.0030
pe1 0.022 0.05 0.0005 0.0005
pe2 0.017 0.04 0.0004 0.0004
Table 2: Measured scatter of the maximum brightness and data and the measured and estimated (italics) amplitudes of their variations.

4.3 Results from CCD observations

A detailed study of the modulation properties of RR Gem during the 2004 observations was given in Paper1. Here we only summarise the global parameters of the modulation in order to compare them with those observed during the photographic and photoelectric observations. The and maximum brightness data determined for CCD light curves of the 2004 and 2005 seasons are listed in Table 15. The values were calculated according to the ephemeris:

where the pulsation period is

The and maximum brightness data of light curves are shown in the right panels of Fig. 4. These plots were constructed using the modulation period valid for the two seasons: (Sódor 2006). The amplitude of the modulation in and maximum brightness are and , respectively (see right panels of Fig. 4). The average maximum brightness of the CCD data was , and the mean pulsation amplitude was .

5 Discussion of long-term changes

The long-term changes in the pulsation and modulation of RR~Gem can be followed in Figs. 4 and 9. The properties of the modulation and pulsation at the studied intervals are summarised in Table 3.

Each of the measured parameters (pulsation and modulation periods and amplitudes, phase difference between maximum brightness and data, and mean maximum brightness) show much larger changes than their uncertainties allow.

5.1 Changes in the modulation properties

In Fig. 4 the modulation in and in maximum brightness, as well as the relation between them, are plotted for all the studied intervals. The top and middle panels show the and maximum brightness values folded with the derived modulation periods of the pg1, pe1, and CCD data and the fitted harmonic curves. Because the modulation period of the pg2 and pe2 data could not be determined due to imperfect data sampling, these data are plotted versus the HJD of the measurements. In the bottom panels each point represents a particular pulsation cycle for which both and maximum brightness could be determined reliably. These plots are independent of the period of the modulation, and reflect the amplitude relation and phase connection of the amplitude and phase modulation components. The plotted ellipses are the combinations of the fitted sine curves of the maximum brightness and data for the pe1, pg1, and CCD data.

It is not easy to make an exact distinction between amplitude and phase modulation of RR Lyrae stars, because the shape of the light curve (i.e. the relative amplitudes and phase differences of the Fourier components) always change more or less during the modulation cycle. Even if a node on the rising branch exists, a slight oscillation in the times of maxima could happen (this is only the case of RR~Gem at the CCD observations; see Fig. 5 in Jurcsik et al. 2006). Especially problematic is how to correctly measure phase modulation. As described in Sect. 2.2.1, the variation measures a combination of the phase and amplitude modulations. For comparison, the bottom panel of the CCD data in Fig. 4 line shows the  maximum brightness connection predicted from test data of the exact amplitude modulation (see Sect. 2.2.1). The actual slope of the observed vs. maximum-brightness plots of the larger amplitude modulations are flatter than for the test data. During the photoelectric observations, when the amplitude modulation was the smallest, if there was any, no modulation in was detected.

The bottom panels in Fig. 4 show that significant changes in the modulation behaviour of RR~Gem occurred during the past 70 years. The most significant change in the nature of the modulation happened at around JD 2 428 800 in conjunction with the sudden change in the pulsation period. Prior to this event, a pronounced phase modulation was accompanied by the modulation in amplitude. After the sudden period increase the character of the modulation had changed. The phase of the middle of the rising branch show a small variation in excess, if any, to what is expected from exact amplitude modulation, while the amplitude of the amplitude modulation changes in the range of an order of magnitude. The bottom panels of Fig. 4 show that, after the sudden pulsation period change, only the strength of the modulation shows significant changes, but its character mostly remains the same.

The parameters plotted in Fig. 9 show that the changes in the amplitudes of the and maximum-brightness variations are the most tightly connected values. The correlated variations of these amplitudes point to the total ‘power’ of the modulation changing with time. Only a few Blazhko stars have extended enough photometric data to follow the changes in the modulation properties. The best-studied case is RR~Lyr itself, where the amplitude of the modulation shows cyclic changes on a 4-year timescale (Szeidl 1976). In other Blazhko stars XZ~Cyg, AR~Her (Almár 1961), AH~Cam (Smith et al. 1999), RW~Dra (Balázs & Detre 1962), there are definite observations showing that the character of their modulation can change considerably, although no detailed study comprising their entire photometric data has been published.

5.2 Changes in the pulsation properties

The most important and surprising finding of this study of long-term changes in RR~Gem is the definite change in the mean maximum brightness and pulsation amplitude in the middle of the photographic observations. We have detected a difference of about 0.12 mag between the mean maximum brightnesses of the pg1 and pg2 data, while the uncertainties allow only less than 0.05 mag, according to the comparison stars photometry (see details in Sect. 2.1.1).

It cannot be excluded that the changes in these parameters between the photographic and photoelectric intervals are also real, but because of the defects of the calibration of the magnitudes, they most probably arise from the differences between the photographic and photoelectric photometries.

The explanation of the Blazhko phenomenon suggested by Dziembowski & Mizerski (2004) involves energy transfer from the fundamental radial to a nonradial mode. It predicts that the larger the amplitude of the modulation, the smaller the pulsation mean amplitude of Blazhko stars. When the modulation was the strongest (pg1 data), the pulsation amplitude of RR~Gem was about 0.1 mag smaller than at any time later, which agrees with the model predictions of Dziembowski & Mizerski (2004).

There are very few long-term Blazhko star observations available that are homogeneous, extended, and accurate enough to investigate similar variations in the modulation and pulsation amplitudes. Extended long-term observations in the same band have been summarised and/or published for RS~Boo (Nagy 1998, and references therein), XZ~Cyg (LaCluyzé et al. 2004, and references therein), RW~Dra (Szeidl et al. 2001, and references therein), XZ~Dra (Jurcsik et al. 2002, and references therein), AR~Her (Almár 1961, and references therein), RR~Lyr (Kolenberg et al. 2006; Smith et al. 2003; Szeidl et al. 1997, and references therein), and RV~UMa (Kovács 1995, and references therein). Nevertheless, only two of them (RR~Lyr and AR~Her) are appropriate for investigation of modulation amplitude and mean maximum-brightness changes.

Almár (1961) showed that in AR~Her the brightest maxima remained at the same brightness level (see Fig. 12 in Almár 1961), although the amplitude of the modulation had changed considerably. This behaviour indicates that the mean brightness was also fainter in AR~Her when the modulation amplitude was large, similar to what we have found in RR~Gem. It supports the correlation of the modulation amplitude and the mean maximum brightness of Blazhko variables.

The behaviour of the modulation of RR~Lyr is, however, different. The mean maximum brightness seems to remain unaltered during the 4-year cycle, while the modulation amplitude changes considerably (see Fig. 6 in Szeidl 1976; Preston 1965, Fig. 2). Nonetheless, the cyclic behaviour of the variation in the modulation amplitude of this star also differs from the irregular changes of RR~Gem and AR~Her, therefore it can have a different physical origin.

The possible differences between the pg and pe data does not allow us to make any conclusions about further changes in the mean maximum brightness and pulsation amplitude values between these observations. During the pe and CCD observations, though the amplitude of the modulation has changed these properties of the pulsation remained the same, indicating that the above-mentioned correlation is not strictly valid.

5.3 Connections between the pulsation and modulation properties

The period changes of the pulsation and modulation indicate a positive correlation between the periods as shown in Fig. 8, although they can be compared only on the bases of three epochs of data. The data define for RR~Gem, which completes the list of Blazhko stars with measured period change of both their periods (LaCluyzé et al. 2004, Table 8) with a new item.

It was shown in Jurcsik et al. (2005b) that the largest possible amplitude of the modulation depends on the pulsation frequency of Blazhko stars. However, at a given pulsation frequency, the modulation has very different amplitudes, from very low values to the extent of the possible maximum value. That those Blazhko variables with long-term photometric observations show significant changes in their modulation amplitudes proves that the measured amplitude of the modulation is just a temporal parameter. It may explain why the pulsation frequency gives only a limit for the possibly highest value of the modulation amplitude and does not determine its actual value.

Most probably the ‘instability’ of the modulation is a general property of Blazhko stars, and any plausible explanation of the phenomenon should give a reason for the irregularity of the observed light curve modulation (see also Chadid & Chapellier 2006).

What the modulation properties (amplitude, phase relation between the amplitude, phase modulation, etc.) tell us about the physics of the phenomenon, if they show considerable changes on larger timescales, is unclear. During the time of the photoelectric observations, RR~Gem showed hardly any modulation. This fact warns us that the modulation can only be temporarily detected in some cases. In contrast to previous expectations, the amplitude of the modulation can also be small (Jurcsik et al. 2006), below the detectability limit of most of the extended surveys. Therefore, the most plausible conclusion is that the modulation is a common, intrinsic property of RR Lyrae stars.

Figure 8: Modulation period versus pulsation period.
data set (error) (error) (error) (error)  (error) (error)  phase(error) 
[d] [d] [mag] [mag] [d] [mag] [rad]
pg1 0.397253(5) 7.20(2) 1.55(2) 10.83(2) 0.0094(4) 0.19(7) -1.8(3)
pg2 0.397317(2)  - 1.77(3) 10.71(5) 0.003 0.15 0.6(7)
pe unf. 0.3973156(35) - - - - - -
pe1 0.3973148(2) 7.28(2) 1.61(2) 10.77(1) 0.0005(4) 0.05(1) -
pe2 0.3973160(2)  - 1.61(2) 10.76(2) 0.0004 0.04 -
CCD 0.3972893(3) 7.216(3) 1.62(1) 10.76(1) 0.0018(2) 0.114(4) 0.0(1)
The errors of the mean maximum magnitudes were determined in two different ways. When the modulation period was found, it corresponds to the zero point uncertainty of the fitted sine curve. For the pg2 and pe2 data, the given errors are the scatters of the maximum points, which overestimate the true errors, if this scatter originates partially from modulation.
 phase denotes the difference between the phases of the maximum magnitude and variation.
As the pulsation period changed during the pg2 and pe2 observations, the given values correspond to the middle of these intervals. Details of the period change are given in Sects. 4.1 and 4.2.
Table 3: Modulation and pulsation properties derived from the observations.
Figure 9: Variation of the pulsation and modulation properties of RR~Gem. The plotted quantities from top to bottom are , pulsation period, pulsation amplitude, average maximum brightness, modulation period, modulation amplitude in , and modulation amplitude in maximum brightness.

6 Conclusions

The rise in atmospheric shock waves during the pulsation of RR Lyrae stars is well-established both observationally and theoretically. Preston (1965) proposed that the depth of the formation of the main shock wave during rising light shifts periodically in the course of the Blazhko cycle. Chadid & Chapellier (2006) report irregular changes in the atmosphere of RR~Lyr during the Blazhko cycle. They find that the residual scatter of the radial velocity curve is largest when the nonlinear effects are the most intense, which is during the shock wave passage across the atmosphere.

The uncertainty of our photographic data and the unfavourable data distribution of both the photographic and photoelectric data sets do not allow us to distinguish between observational errors and irregularities. However, the detected systematic variations give clear evidence of the light curve modulation. Analysis of the 70-year long photometric observations of RR~Gem has revealed some new and important properties of the Blazhko modulation.

RR~Gem is the first Blazhko RR Lyrae star where the amplitude of the pulsation showed definite changes during the observation times. An unambiguous amplitude increase happened together with a sudden pulsation period increase around JD 2 428 800. Before this event, a clear and strong modulation was observable with an amplitude of nearly 0.2 mag in maximum brightness and 14 min in phase on the rising branch. After the sudden pulsation period change, the strength of the modulation decreased, along with an increase in the mean pulsation amplitude. The photoelectric observations of show clear but weak modulation only in amplitude, with a slightly longer Blazhko period. The CCD observations show modulation with an amplitude of about 0.1 mag in band and with very weak phase modulation.

Cyclic changes in the modulation and pulsation properties are known for two Blazhko variables, RR~Lyr (Szeidl 1976) and XZ~Dra (Jurcsik et al. 2002). Supposing that the changes observed in RR Gem are also a part of a cyclic behaviour on a century-long timescale, both the pulsation and modulation periods are expected to decrease, while the amplitude of the modulation is expected to increase in the future. If a sudden jump in the pulsation and modulation properties of RR Gem, similar to that observed in 1937, would happen in the future that would give us a unique opportunity to study the connections of the pulsation and modulation properties of Blazhko variables in detail. The extremely short modulation period (7.2 d) makes it easier to check the modulation amplitude of RR Gem annually.

Acknowledgements.
We would like to thank the referee, M. Chadid, for her useful suggestions on this paper. This research has made use of the SIMBAD database, operated at the CDS Strasbourg, France, and the GEOS RR Lyrae database. The financial support of OTKA grants T-043504, T-046207, and T-048961 is acknowledged.

References

  • Agerer & Huebscher (2000) Agerer, F., & Huebscher, J. 2000, IBVS, 4912
  • Agerer & Huebscher (2002) Agerer, F., & Huebscher, J. 2002, IBVS, 5296
  • Agerer & Huebscher (2003) Agerer, F., & Huebscher, J. 2003, IBVS, 5485
  • Agerer et al. (2001) Agerer, F., Dahm, M., & Huebscher, J. 2001, IBVS, 5017
  • Alania (1956) Alania, I. F. 1956, AC, 173, 20
  • Almár (1961) Almár, I. 1961, Communications of the Konkoly Observatory, No. 51
  • Balázs (1960) Balázs-Detre, J. 1960, Kleine Veröffentlichungen der Remeis-Sternwarte Bamberg, Nr. 27, p. 26
  • Balázs & Detre (1962) Balázs-Detre, J., & Detre, L. 1962, Kleine Veröffentlichungen der Remeis-Sternwarte Bamberg, Nr. 34, p. 90
  • Baldwin & Samolyk (2003) Baldwin, M. E., & Samolyk, G. 2003, Observed Maxima Timings of RR Lyrae Stars, Number 1 (Cambridge: AAVSO)
  • Batyrev (1950) Batyrev, A. A. 1950, VS, 7, 243
  • Batyrev (1962) Batyrev, A. A. 1962, VS, 14, 230
  • Blazhko (1907) Blazhko, S. 1907, Astron. Nachr., 173, 325
  • Blazhko (1935) Blazhko, S. 1935, Tadzhik Eph., 1935
  • Braune et al. (1977) Braune, W., Huebscher, J., & Mundry, E. 1977, AN, 298, 121
  • Chadid (2000) Chadid, M. 2000, A&A, 359, 991
  • Chadid & Chapellier (2006) Chadid, M., Chapellier, E. 2006, A&A, 456, 305
  • Chudovichev (1930) Chudovichev, N. I. 1930, VS, 2, 58
  • Detre (1970) Detre, L. 1970, Transactions of the IAU XIV A. Reports on Astronomy, p.259
  • Dubiago (1930) Dubiago, A. 1930, Astron. Nachr., 239, 15
  • Dziembowski & Mizerski (2004) Dziembowski, W. A., Mizerski, T. 2004, Acta Astron., 54, 363
  • Epstein (1969) Epstein, I. 1969, AJ, 74, 1131
  • ESA (1997) ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP-1200
  • Fitch et al. (1966) Fitch, W. S., Wisniewski, W. Z., & Johnson, H. L. 1966, Comm. Lunar Plan. Lab., 5, 71
  • Gensler (1998) Gensler, J. 1998, BAV Mitt., 113
  • Guriev (1958) Guriev, N. I. 1958, Tadzhik Trudy, 7
  • Graff (1922) Graff, K. 1922, Astron. Nachr., 217, 305
  • Huebscher (2000) Huebscher, J. 2000, BAV Mitt., 131
  • Huebscher (2003) Huebscher, J. 2003, BAV Mitt., 157
  • Huebscher (2005) Huebscher, J. 2005, IBVS, 5643
  • Huebscher et al. (2005) Huebscher, J., Paschke, A., & Walter, F. 2005, IBVS, 5657
  • Ivanov (1930a) Ivanov, N. 1930a, Astron. Nachr., 237, 291
  • Ivanov (1930b) Ivanov, N. 1930b, Proceedings of II, III, IV Astronomical Congresses (1920-1928), Leningrad, p151
  • Jurcsik et al. (2001) Jurcsik, J., Clement, C., Geyer, E. H., & Domsa, I. 2001, AJ, 121, 951
  • Jurcsik et al. (2002) Jurcsik, J., Benkő, J. M., & Szeidl, B. 2002, A&A, 396, 539
  • Jurcsik et al. (2005a) Jurcsik, J., Sódor, Á., Váradi, M. et al. 2005a, A&A, 430, 1049 (Paper1)
  • Jurcsik et al. (2005b) Jurcsik, J., Sódor, Á, Váradi, M. 2005b, IBVS, 5666
  • Jurcsik et al. (2006) Jurcsik, J., Szeidl, B., Sódor, Á. et al. 2006, AJ, 132, 61
  • Kanischeva & Lange (1971) Kanishcheva, R. K., & Lange, G. A. 1971, VS supp. 1, 107
  • Kleissen (1939) Kleissen, E. 1939, AN, 267, 137
  • Kolláth (1990) Kolláth, Z. 1990, Occ. Techn. Notes Konkoly Obs., No. 1, http://www.konkoly.hu/staff/kollath/mufran.html
  • Kovács (1995) Kovács, G. 1995, A&A, 295, 693
  • Kolenberg et al. (2006) Kolenberg, K., Smith, H. A., Gazeas, K. D. et al. 2006, A&A, 459, 577
  • LaCluyzé et al. (2004) LaCluyzé, A., Smith, H. A., Gill, E.-M., Hedden, A., Kinemuchi, K. et al. 2004, AJ, 127, 1653
  • Lange (1935) Lange, G. A. 1935, Tadzhik Eph. 1935, 23
  • Lange (1969) Lange, G. A. 1969, AC, 503, 7
  • Lange et al. (1976) Lange, G. A., Motrich, V. D., Firmanyuk, D. N., & Tsesevich, V. P. 1976, AC, 900, 5
  • Le Borgne (2004) Le Borgne, J. F. 2004, GEOS Circ., RR22
  • Le Borgne et al. (2005) Le Borgne, J. F., Klotz, A., & Boër, M. 2005, IBVS 5650
  • Le Borgne et al. (2006) Le Borgne, J. F., Klotz, A., & Boër, M. 2006, IBVS 5717
  • Liu & Janes (1989) Liu, T., & Janes, K. A. 1989, ApJS, 69, 593
  • Luizet (1930) Luizet, M. 1930, Lyon Bull., 12, 94
  • Mandel (1960) Mandel, O. E. 1960, AC, 209, 28
  • Nagy (1998) Nagy, A. 1998, A&A, 339, 440
  • Preston (1965) Preston, G. W., Smak, J., Paczynski, B. 1965, ApJS, 12, 99
  • Smith et al. (1999) Smith, H. A., Barnett, M., Silbermann, N. A., & Gay, P. 1999, AJ, 118, 572
  • Smith et al. (2003) Smith, H. A., Church, J. A., Fournier, J. et al. 2003, PASP, 115, 43
  • Sódor (2006) Sódor, Á. 2006, PADEU, 17, 115
  • Stepien (1972) Stepien, K. 1972, Acta Astron., 22, 175
  • Szeidl (1976) Szeidl, B. 1976, in IAU Colloq. 29, Multiple Periodic Variable Stars, ed. W. S. Fitch (Dordrecht: Reidel), 133
  • Szeidl et al. (1997) Szeidl, B., Guinan, E. F., Oláh, K., & Szabados, L. 1997, Communications of the Konkoly Observatory, No. 99
  • Szeidl et al. (2001) Szeidl, B., Oláh, K., Barlai, K., & Szabados, L. 2001, Communications of the Konkoly Observatory, No. 102
  • Tsesevich (1969) Tsesevich, V. P. 1969, ”RR Lyrae stars”, Israel program for scientific translation. 1966 for the original version in russian, Izn. Naukova Dumka, Kiev
  • Vandenbroere (1997) Vandenbroere, J. 1997, GEOS Circular, 14
  • Vandenbroere (1998) Vandenbroere, J. 1998, GEOS Circular, 15
  • Vandenbroere (1999) Vandenbroere, J. 1999, GEOS Circular, 16
  • Vandenbroere (2001) Vandenbroere, J. 2001, GEOS Circular, 17
  • Vandenbroere (2003) Vandenbroere, J. 2003, GEOS Circular, 19
  • Vandenbroere (2005) Vandenbroere, J. 2005, GEOS Circular, 23
  • Waterfield (1927) Waterfield, R. L. 1927, Bulletin of the Harvard College Observatory, 848, 25
  • Wils et al. (2006) Wils, P., Lloyd, Ch., & Bernhard, K. 2006, MNRAS, 368, 1757
  • Wolfram (1996) Wolfram, S. 1996, The Mathematica Book (Wolfram Media – Cambridge Univ. Press, Champaign – Cambridge)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
60179
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description