Issue |
A&A
Volume 496, Number 1, March II 2009
|
|
---|---|---|
Page(s) | 217 - 222 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/200811168 | |
Published online | 30 January 2009 |
False periods in complex chaotic systems
V. Votruba1,2 - P. Koubský1 - D. Korcáková1 - F. Hroch2
1 - Astronomický ústav, Akademie ved Ceské republiky v.v.i., 251 65 Ondrejov, Czech Republic
2 - Ústav teoretické fyziky a astrofyziky, Prirodovecká fakulta, Masarykova univerzita, Kotlárská 2, 611 37 Brno, Czech Republic
Received 16 October 2008 / Accepted 19 December 2008
Abstract
Aims. Astrophysical objects frequently exhibit some irregularities or complex behaviour in their light curves. We focus primarily on hot stars, where both radial and non-radial pulsations are observed. One of the primary research goals is to determine physical parameters of stellar pulsations by analyzing their light curves or spectra, focusing on periodic or quasiperiodic behaviour.
Methods. We analyse the feasibility of classical methods for period searches in a nonlinear chaotic system, such as the Rössler system, where a period does not exist at all. As an astrophysical application of the chaotic system, we utilize a simple model of stellar pulsation with two different sets of parameters corresponding to periodic and chaotic behaviour. For both models we create a synthetic signal, and then apply widely used methods for period finding, such as the phase dispersion method and periodograms. For comparison, a quasi-periodic signal is employed as well.
Results. The period analysis indicates periods even for the chaotic signal. Such periods are apparently spurious. This implies that it is very problematic to distinguish chaotic and quasiperiodic process by such an analysis only.
Key words: stars: oscillations - stars: variables: Cepheids - chaos
1 Introduction
Variability of astronomical sources like stars is an important aspect of their observation. We can detect variability in their light and/or spectra on different timescales. Dynamical pulsations, radial or non-radial, may be one of the main sources of intrinsic variability in many types of stars. As was shown by Eddington (1919), stars may exhibit pulsation due to the unstable layers in the stellar envelope. The first theoretical modelling of pulsations predicted regular oscillations and such modelling was succesfully applied to objects like classical Cepheids. But later, observations indicated other types of pulsating stars, e.g. Population II Cepheids, which in some cases exhibit irregular patterns in observational data.
There was no reliable explanation of such variability until the pioneering work of Baker et al. (1966), which showed that a simple one-zone model of stellar oscillators with a nonlinear adiabatic term can in some specific conditions produce irregular pulsation. Similar results were also obtained by Moore & Spiegel (1966), Auvergne & Baglin (1985) and Auvergne (1988). The disadvantage of these models is their dynamically unstable equilibrium state, which seems to be nonrealistic for stars. Takeuti (1987) and later Saitou & Takeuti (1989) showed that similar, yet dynamically stable and more precise one-zone model also produced chaotic solutions, with a strange attractor similar to the Rössler one (see Rössler 1976). Detailed, more realistic calculations of the state equation for helium ionization in a one zone model, which was performed by Buchler & Regev (1982), confirmed that pulsation instability may lead to chaotic behaviour.
Goupil et al. (1986) studied light curves of the pulsating white dwarf PG 1351+489 and
found strong evidence of period doubling. Period doubling is a
scenario where a system that primarily oscillates with a fundamental frequency
,
after change of the control parameter oscillates with frequency
(see Gurel & Rössler 1979; Helleman 1980). Cascades of period doubling is a
typical transition routes to chaotic behaviour. Such an effect is typical of chaotic systems.
Serre et al. (1996) analysed long time-series of photometric observations of the variable star R Scuti, one
possible candidate pulsating star which exhibits a strange irregular light
curve. By using the method of global flow reconstruction (see more about global
time series predictions in Brown 1993; Abarbanel et al. 1994), they were able to
reproduce the observation with the chaotic dynamics very well. More details and examples of chaotic pulsations can be found in Buchler (1987).
These clearly show that there is both theoretical and experimental evidence of quasiperiodic or chaotic behaviour in variable stars. The problem is how to distinguish between these two possibilities, if only limited time series are available, often very unevenly sampled. In other words, what will be the result, if classical astronomical tools for the period searches, such as periodograms, the phase dispersion minimalization method or more sophisticated techniques are applied to data originating from quasiperiodic or chaotic sources. In the next paragraphs we tackle the problem using synthetic data.
![]() |
Figure 1: Left panel: chaotic attractor constructed from a numerical solution of the ODE system (1) with the parameters from Table 1. - chaotic case. Right panel: two-limit cycle attractor constructed from numerical solution of (1) with the parameters from Table 1. - periodic case. |
Open with DEXTER |
2 Generation of synthetic time series
As an example of irregular stellar variability due to pulsation we use the model of Tanaka & Takeuti (1988). In this model, the set of ordinary differential equations (ODE) for a stellar oscillator with nonlinear, nonadiabatic terms, reads
where the first two equations represent equations of motion and the last one is a state equation with a nonadiabatic term. We can adjust the behaviour of the model using different values of

Table 1: Parameters of the stellar pulsation model.
Table 2: Parameters for the quasiperiodic signal.
![]() |
Figure 2:
Upper left panel: x(t) solution of the ODE system (1) for the periodic case with a two-limit cycle. Upper right panel: x(t) solution of Eq. (2) for the quasiperiodic case. Lower left panel: x(t) solution of the ODE system (1) for the chaotic case. Lower right panel: Gaussian noise (3) with parameters
|
Open with DEXTER |
We use the initial conditions (
x0 = 0.37,
y = -0.249,
z = -0.174) for the chaotic solution.
The initial conditions are very close to the orbit of the attractor in phase space. We obtain a time-dependence of the state variables
x(t),y(t),z(t) in the form of a time series with an equidistant time step. We choose one state variable, namely x(t), as a representation of the measured variable. To simulate a synthetic signal from the quasiperiodic process, a different approach is employed. The signal is modeled by the following equation
where a1,a2,a3 are amplitudes of incommensurate modes with frequencies f1,f2,f3. Their values are summarized in Tables 1 and 2. These parameters are carefully chosen in order to obtain similar dependencies of x(t) and

where



Data for synthetic time series satisfy the condition t>300, i.e. reasonably far away from the starting point (t=0). It is necessary to make sure that the solution is relaxed away from the influence of the initial conditions.
It can be readily seen that chaotic and quasiperiodic oscillations (see Fig. 2) have very similar patterns, different to the simple periodic signal (two-limit cycle).
![]() |
Figure 3:
Correlation integrals (4) for different lengthscales |
Open with DEXTER |
![]() |
Figure 4: Upper panels: power spectra constructed for quasiperiodic oscillation ( on the left) described by (2) and the chaotic case of stellar oscillation ( on the right) described by (1). Bottom panels: periodogram constructed using the phase dispersion minimalization method for the quasiperiodic case ( on the left) and chaotic case ( on the right). For the frequency we used arbitrary units. |
Open with DEXTER |
![]() |
Figure 5:
Left upper panel: phase diagram for the quasiperiodic case and the dominant frequency
|
Open with DEXTER |
3 Analysis of the signal
3.1 Classical period analysis
The most common techniques used for period searches are the power spectral analysis (PWS), the modified Lomb periodogram (LPD), the phase dispersion minimalization (PDM), and the string length method (SLM). In the first case, discrete Fourier transformation is applied to a dataset resulting in a power spectrum (for more details see e.g. Scargle 1982). For unevenly sampled series it is better to use a modification of PWS, method LPD (see for more details Lomb 1976). The third method binds data in a grid of phases and minimizes phase dispersion statistics (see for more details Stellingwerf 1978; Lafler& Kinman 1965). The SLM method is similar to the previous one (see Dworetsky 1983), and the LPD method behaves like the PWS. So we limit on the two methods, PWS and PDM only. Obviously, PDM and PWS techniques applied on a periodic signal detect corresponding frequencies without any doubt, but, as we will show, a different situation arises for quasiperiodic and chaotic regimes.
- 1. Quasiperiodic signal. The PWS analysis of the quasiperiodic signal
produced by known incommensurate frequencies (see Table 2)
detects the frequencies used. Moreover, it is possible to recognize the leading
frequency, which is in our case
. The noise in the periodogram in Fig. 4 is caused by the simulated Gaussian noise of the signal. The PDM method yields a similar estimate of the input frequencies, and the results also show some other frequencies that are subharmonics or harmonics of the frequency set. The phase diagram in Fig. 5 (on the left) is folded with the leading frequency f2. Clearly, the structure is not a simple sine wave. It is a quasiperiodic, and hence it is composed of a set of harmonics with slightly different phase. After accounting for Gaussian noise, this difference is smeared out.
- 2. Chaotic signal. The results become much more interesting for the chaotic signal. Again, both PWS and PDM methods detect periods, as can be clearly seen in Fig. 4. After some checks we can recognize the two longest independent periods.The power spectrum of the chaotic signal shows the most significant frequencies together with the frequency noise. Nevertheless, low-level noise is present even if Gaussian noise is not applied. It is the consequence of the chaotic behaviour.
The above results clearly show that it is very difficult to distinguish between chaotic and quasiperiodic cases using period analysis only. Thus, we must look for other methods, which can help to resolve both processes.
3.2 Nonlinear analysis
For detailed analysis of the signal that can be generated by nonlinear
chaotic processes, we can use methods of nonlinear time series analysis. These
methods were developed for classification and identification of chaotic
phenomenon in time series of the experimental data (see review Abarbanel et al. 1993). As a indication of chaotic behaviour we can use topological invariants, namely the correlation dimension and Lyapunov exponents. Algorithms for numerical estimates of such invariants from time series are described by Abarbanel et al. (1993) and by Kantz & Schreiber (2004). In our case the indicator of chaotic behaviour is the correlation integral defined (see more detail Kantz & Schreiber 2004) as
where m is the embedding dimension,




![]() |
(5) |
This value agrees with the typical value of the Rössler type attractor. In the case of the quasiperiodic signal we cannot identify the plateau region (right panel Fig. 3). This indicates that there is no self-similar geometry for this type of the signal.
4 Summary
Analysis of synthetic signals simulating quasiperiodic and chaotic data led to the following findings:
- 1)
- In the realistic case (finite number of observations with Gaussian noise), it is impossible to distinguish between quasiperiodic and chaotic behaviour by classical methods of period analysis.
- 2)
- When the data contain large amounts of noise, the chaotic or quasiperiodic character of the signal can be overlooked and a simple periodic solution with the strongest frequency can be determined from the period analysis.
- 3)
- In order to definitely rule out presence of chaos it is neccessary to use a nonlinear time series analysis. Thus, we must reconstruct the topologically similar phase portrait of the system using the technique described by Abarbanel et al. (1993). Then we must determine the invariants, such as the correlation dimension or Lyapunov coefficients from the reconstructed time series, and from their values we can decide if deterministic chaos is present in the data or not. However, both these procedures need better quality of data (noise, length of time series).
- 4)
- We would like to draw attention to the fact that irregular chaotic behaviour may be more usual in astronomical data than thought before, and that some systems that were identified as quasiperiodic may actually be chaotic ones.
Acknowledgements
We thank S. Shore and T. Dytrych for helpful comments on our analysis. This research has made use of NASA's Astrophysics Data System. This work was supported by the grant CZ 205/09/P476. The Astronomical Institute Ondrejov is supported by project AV0 Z10030501.
References
- Abarbanel, H. D. I., Brown, R., Sidorowich, J. J., & Tsimring, L. S. 1993, Rev. Mod. Phys., 65, 1331 [NASA ADS] [CrossRef] (In the text)
- Abarbanel, H. D. I., Carroll, T. A., Pecora, L. M., Sidorowich, J. J., & Tsimring, L. S. 1994, Phys. Rev. E, 49, 1840 [NASA ADS] [CrossRef] (In the text)
- Auvergne, M., & Baglin, A. 1985, A&A, 142, 388 [NASA ADS] (In the text)
- Auvergne, M. 1988, A&A, 204, 341 [NASA ADS] (In the text)
- Baker, N. H., Moore, D. W., & Spiegel, E. A. 1966, ApJ, 71, 844 (In the text)
- Brown, R. 1993, Phys. Rev. E, 47, 3962 [NASA ADS] [CrossRef] (In the text)
- Buchler, J. R., & Regev, O. 1982, ApJ, 263, 312 [NASA ADS] [CrossRef] (In the text)
- Buchler, J. R. 1987, Chaotic Phenomena in Astrophysics, ed J. R. Buchler, & H. Eichhorn, Annals New York Academy of Sciences, 37 (In the text)
- Dworetsky, M. M. 1983, MNRAS, 203, 917 [NASA ADS] (In the text)
- Eddington, A. S. 1919, MNRAS, 79, 2 (In the text)
- Goupil, M. J., Auvergne, A., & Baglin, A. 1986, A&A, 196, L13 [NASA ADS] (In the text)
- Gurel, O., & Rössler, O. E. 1979, Bifurcation Theory and Applications to Scientific Disciplines, Annals New York Academy of Science, 316
- Hamming, R. W. 1962, Numerical methods fo Scientist and Engineers, (McGraw-Hill) (In the text)
- Helleman, R. ed. 1980, Nonlinear dynamics, Annals New York Academy of Science, 357
- Kantz, H., & Schreiber, T. 2004, Nonlinear Time Series Analysis (Cambridge University Press) (In the text)
- Lafler, J., & Kinman, T. D. 1965, ApJS, 11, 216 [NASA ADS] [CrossRef]
- Lomb, N. R. 1976, Ap&SS, 39, 447 [NASA ADS] [CrossRef] (In the text)
- Moore, D. W., & Spiegel, E. A. 1966, ApJ, 143, 871 [NASA ADS] [CrossRef] (In the text)
- Omodei, N., Bellazzini, J., & Montangero, S. 2004, A&A, 414, 1177 [NASA ADS] [CrossRef] [EDP Sciences]
- Regev, O. 2006, Chaos and Complexity in Astrophysics (Cambridge University Press) (In the text)
- Rössler, O. E. 1976, , 57, 397 (In the text)
- Saitou, M., & Takeuti, M. 1989, PASJ, 41, 297 [NASA ADS] (In the text)
- Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] (In the text)
- Serre, T., Kolláth, Z., & Buchler, J. R. 1996, A&A, 311, 833 [NASA ADS] (In the text)
- Schreiber, T. 1993, Phys. Rev. E, 47, 2401 [NASA ADS] [CrossRef] (In the text)
- Stellingwerf, R. F. 1978, ApJ, 224, 953 [NASA ADS] [CrossRef]
- Takeuti, M. 1987, Ap&SS, 136, 129 [NASA ADS] [CrossRef] (In the text)
- Tanaka, Y., & Takeuti, M. 1988, Ap&SS, 148, 229 [NASA ADS] [CrossRef] (In the text)
All Tables
Table 1: Parameters of the stellar pulsation model.
Table 2: Parameters for the quasiperiodic signal.
All Figures
![]() |
Figure 1: Left panel: chaotic attractor constructed from a numerical solution of the ODE system (1) with the parameters from Table 1. - chaotic case. Right panel: two-limit cycle attractor constructed from numerical solution of (1) with the parameters from Table 1. - periodic case. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Upper left panel: x(t) solution of the ODE system (1) for the periodic case with a two-limit cycle. Upper right panel: x(t) solution of Eq. (2) for the quasiperiodic case. Lower left panel: x(t) solution of the ODE system (1) for the chaotic case. Lower right panel: Gaussian noise (3) with parameters
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Correlation integrals (4) for different lengthscales |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Upper panels: power spectra constructed for quasiperiodic oscillation ( on the left) described by (2) and the chaotic case of stellar oscillation ( on the right) described by (1). Bottom panels: periodogram constructed using the phase dispersion minimalization method for the quasiperiodic case ( on the left) and chaotic case ( on the right). For the frequency we used arbitrary units. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Left upper panel: phase diagram for the quasiperiodic case and the dominant frequency
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.