In mathematics, particularly in dynamical systems, Arnold tongues (named after Vladimir Arnold)[1] [2] are a pictorial phenomenon that occur when visualizing how the rotation number of a dynamical system, or other related invariant property thereof, changes according to two or more of its parameters. The regions of constant rotation number have been observed, for some dynamical systems, to form geometric shapes that resemble tongues, in which case they are called Arnold tongues.[3]
Arnold tongues are observed in a large variety of natural phenomena that involve oscillating quantities, such as concentration of enzymes and substrates in biological processes[4] and cardiac electric waves. Sometimes the frequency of oscillation depends on, or is constrained (i.e., phase-locked or mode-locked, in some contexts) based on some quantity, and it is often of interest to study this relation. For instance, the outset of a tumor triggers in the area a series of substance (mainly proteins) oscillations that interact with each other; simulations show that these interactions cause Arnold tongues to appear, that is, the frequency of some oscillations constrain the others, and this can be used to control tumor growth.
Other examples where Arnold tongues can be found include the inharmonicity of musical instruments, orbital resonance and tidal locking of orbiting moons, mode-locking in fiber optics and phase-locked loops and other electronic oscillators, as well as in cardiac rhythms, heart arrhythmias and cell cycle.[5]
One of the simplest physical models that exhibits mode-locking consists of two rotating disks connected by a weak spring. One disk is allowed to spin freely, and the other is driven by a motor. Mode locking occurs when the freely-spinning disk turns at a frequency that is a rational multiple of that of the driven rotator.
The simplest mathematical model that exhibits mode-locking is the circle map, which attempts to capture the motion of the spinning disks at discrete time intervals.
Arnold tongues appear most frequently when studying the interaction between oscillators, particularly in the case where one oscillator drives another. That is, one oscillator depends on the other but not the other way around, so they do not mutually influence each other as happens in Kuramoto models, for example. This is a particular case of driven oscillators, with a driving force that has a periodic behaviour. As a practical example, heart cells (the external oscillator) produce periodic electric signals to stimulate heart contractions (the driven oscillator); here, it could be useful to determine the relation between the frequency of the oscillators, possibly to design better artificial pacemakers. The family of circle maps serves as a useful mathematical model for this biological phenomenon, as well as many others.[6]
The family of circle maps are functions (or endomorphisms) of the circle to itself. It is mathematically simpler to consider a point in the circle as being a point
x
2\pi
2\pi
[0,2\pi]
\thetai+1=g(\thetai)+\Omega
where
\Omega
g
g(\theta)=\theta
\theta
\Omega
\Omega
The particular circle map originally studied by Arnold,[8] and which continues to prove useful even nowadays, is:
\thetai+1=\thetai+\Omega+
K | |
2\pi |
\sin(2\pi\thetai)
where
K
\thetai
1
K
\Omega
\Omega=1/3
K
Another way to view the circle map is as follows. Consider a function
y(t)
a
z(t)=c+b\sin(2\pit)
\{tn\}
This model tells us that at time
tn-1
y(tn-1)=c+b\sin(2\pitn-1)
y
tn
y
\begin{align} 0&=y(tn-1)-a ⋅ (tn-tn-1)\\[0.5em] 0&=\left[c+b\sin(2\pitn-1)\right]-atn+atn-1\\[0.5em] tn&=
1 | |
a |
\left[c+b\sin(2\pitn-1)\right]+tn-1\\[0.5em] tn&=tn-1+
c | |
a |
+
b | |
a |
\sin(2\pitn-1) \end{align}
\Omega=c/a
K=2\pib/a
tn=tn-1+\Omega+
K | |
2\pi |
\sin(2\pitn-1).
argues that this simple model is applicable to some biological systems, such as regulation of substance concentration in cells or blood, with
y(t)
In this model, a phase-locking of
N:M
y(t)
N
M
z(t)
N/M
Consider the general family of circle endomorphisms:
\thetai+1=g(\thetai)+\Omega
g(\theta)=\theta+(K/2\pi)\sin(2\pi\theta)
f(\theta)
\thetai+1=f(\thetai)=\thetai+\Omega+
K | |
2\pi |
\sin(2\pi\thetai).
P1.
f
K<1
K
\thetai
f
f'(\theta)=1+K\cos(2\pi\theta)
K<1
P2. When expanding the recurrence relation, one obtains a formula for
\thetan
\thetan=\theta0+n\Omega+
K | |
2\pi |
n-1 | |
\sum | |
i=0 |
\sin(2\pi\thetai).
P3. Suppose that
\thetan=\theta0\bmod1
n
\thetai
M=(\thetan-\theta0) ⋅ 1
n:M
P4. For any
p\inN
f(\theta+p)=f(\theta)+p
f(\theta+p)=f(\theta)\bmod{1}
\thetai
1
P5 (translational symmetry).[9] Suppose that for a given
\Omega
n:M
\Omega'=\Omega+p
p
n:(M+np)
\theta0,...,\thetan
\Omega
\Omega'=\Omega+p,p\inN
P6. For
K=0
\Omega
\Omega=p/q\inQ
q:p
For small to intermediate values of K (that is, in the range of K = 0 to about K = 1), and certain values of Ω, the map exhibits a phenomenon called mode locking or phase locking. In a phase-locked region, the values θn advance essentially as a rational multiple of n, although they may do so chaotically on the small scale.
The limiting behavior in the mode-locked regions is given by the rotation number.
\omega=\limn\toinfty
\thetan | |
n |
.
which is also sometimes referred to as the map winding number.
The phase-locked regions, or Arnold tongues, are illustrated in yellow in the figure to the right. Each such V-shaped region touches down to a rational value Ω = in the limit of K → 0. The values of (K,Ω) in one of these regions will all result in a motion such that the rotation number ω = . For example, all values of (K,Ω) in the large V-shaped region in the bottom-center of the figure correspond to a rotation number of ω = . One reason the term "locking" is used is that the individual values θn can be perturbed by rather large random disturbances (up to the width of the tongue, for a given value of K), without disturbing the limiting rotation number. That is, the sequence stays "locked on" to the signal, despite the addition of significant noise to the series θn. This ability to "lock on" in the presence of noise is central to the utility of the phase-locked loop electronic circuit.
There is a mode-locked region for every rational number . It is sometimes said that the circle map maps the rationals, a set of measure zero at K = 0, to a set of non-zero measure for K ≠ 0. The largest tongues, ordered by size, occur at the Farey fractions. Fixing K and taking a cross-section through this image, so that ω is plotted as a function of Ω, gives the "Devil's staircase", a shape that is generically similar to the Cantor function.One can show that for K<1, the circle map is a diffeomorphism, there exist only one stable solution. However, as K>1 this holds no longer, and one can find regions of two overlapping locking regions. For the circle map it can be shown that in this region, no more than two stable mode locking regions can overlap, but if there is any limit to the number of overlapping Arnold tongues for general synchronised systems is not known.
The circle map also exhibits subharmonic routes to chaos, that is, period doubling of the form 3, 6, 12, 24,....
The Chirikov standard map is related to the circle map, having similar recurrence relations, which may be written as
\begin{align} \thetan+1&=\thetan+pn+{
K | |
2\pi |
Arnold tongues have been applied to the study of