Effect of noise on the golden-mean quasiperiodicity at the chaos threshold |
The circle map
is a fundamental model describing many dissipative nonlinear systems like a self-oscillator under external periodic driving, Josephson junction in a microwave field, space charge density waves in solid-state physics, a pendulum with periodic external force. In studies concerning medical and biological problems, the circle map appears as a model of cardiac rhythms.
The circle map is not only a qualitative model, but a
representative of a class of quantitative universality
associated with a transition to chaos via quasiperiodic motion.
The most well studied case relates to the ratio of basic frequencies, or the rotation number,
equal to the golden mean constant
In real physical systems, effect of noise takes place inevitably. Let us turn to a circle map with added random force:
where x_{n} is a sequence of statistically independent random numbers with a zero mean and fixed variance s, e is a parameter of intensity of noise supposed to be small.
Below we show charts of the Lyapunov exponent for the circle map
on the parameter plane
In presence of noise, the periodic or quasiperiodic dynamics does not occur in exact sense, but the picture of characteristic domains on the Lyapunov chart remains visible, at least for small and moderate noise intensities. We can speak either on "noisy periodic regime", when the Lyapunov exponent is notably negative, on "noisy quasiperiodic", if it is close to zero, or of "noisy chaotic" as the Lyapunov exponent is positive. The Lyapunov charts allow to distinguish visually these regimes. In a sense, the effect of noise looks rather trivial: the subtle details of the picture appear to be smeared off.
The following figure shows phase portraits of the attractors on the iteration diagrams, plots of x_{n+Fk} versus x_{n}, where F_{k} is one of the Fibonacci numbers. The diagrams are presented for F_{6} = 8 at the parameter values K and r associated with the critical point GM.
The first picture corresponds to a very weak noise, the second and the third to larger intesities (see indicated values of e). Observe successive smearing of subtle and than more large scale structural details of the attractor as the noise intensity grows.
As known, without noise near the critical point GM
definite scaling regularities take place
[Shenker, Ostlund, Rand, Feigenbaum, Kadanoff], hence, one can think that
effect of noise also must be associated with some properties of
universality and scaling.
The renormalization group (RG) analysis suggested by Hamm and Graham
[Hamm and Graham, 1992] leads to a conclusion of existence
of a universal constant, which indicates how much the noise intensity has to be
decreased to observe each next level of the small-scale structure in a
vicinity of the GM critical point:
More general formulation of the scaling property is the following.
Let us suppose that with deflection of r from the value corresponding to the
golden-mean rotation number
Dr and K close to
K_{c}=1 in presence of noise
e a certain regime of behavior
takes place. Then, at
Dr decreased by factor
Let us consider some computer illustrations of the mentioned scaling. We start with a series of plots, each of which depicts a noisy attractor on the iteration diagrams in coordinates x(nF_{k}), x((n+1)F_{k}, where F_{k} is a Fibonacci number with k=6, 7, 8. From one picture to another the noise parameter is decreased by factor g. Values of K and r correspond to the critical point GM. On each plot we select a fragment near the origin shown separately with magnification. Similarity of the objects observed on these magnified pictures illustrates the expected scaling property in the system with noise.
Let us turn now to diagrams shown the rotation number
dependence on the parameter r at K=1.
The basic plot and the first inset relate to the circle map with noise
at e=0.1. The other pictures show fragments with
growing magnification along the horizontal and vertical axes, respectively,
with factors
The next diagram depicts a dependence of the Lyapunov exponent on the noise intensity at K and r corresponding to the GM critical point. As seen, presence of noise leads to decrease of the Lyapunov exponent, i.e. facilitates stabilization (in contrast to the situation of the Feigenbaum transition, for which the noise gives rise no increase of the Lyapunov exponent). In the inset a fragment of the plot is shown with rescaling along horizontal by factor g and along vertical by factor W=(5^{1/2}-1)/2. Similarity of both pictures illustrates the expected scaling property.
One more way to demonstrate scaling consists in plotting of the Lyapunov exponent versus the noise intensity in the double logarithmic scale. Observe that the dots are located in average along a straight line with the slope coefficient log_{g}W=0.5759, where W=(5^{1/2}+1)/2.
Another series of illustrations contains the Lyapunov charts
relating to the critical value of the parameter K=1.
The horizontal axis is that of parameter r,
and the vertical one of the noise intensity parameter
e. The gray tones from dark to light
code Lyapunov exponent values from minus infinity to zero, and black corresponds to
positive values.
For each next picture for illustration of the scaling property
the horizontal and vertical scales are changed by a factor, respectively,
Finally, the last series of diagrams represents the Lyapunov charts illustrating the scaling property of the picture of the Arnold tongues. For this, it is necessary to introduce certain special local coordinates (C_{1},C_{2}) on the parameter plane near the GM critical point: One coordinate axis is directed along the critical line K=1, and the second one along the curve of constant rotation number defined in absence of the noise.
The main diagram is a Lyapunov chart on the plane (r,K) at the noise intensity e=0.03. The first inset shows an interior of the selected parallelogram in the local coordinates, which are expressed wia the parameters of the circle map as follows
The other two pictures are obtained with the noise intensity parameter decreased, respectively, by g and g^{2} times. The horizontal and vertical scales correspond to magnification from one picture to another by factors d_{1} and d_{2}, with respective redefinition of the color coding rule for the Lyapunov exponent.