Chirikov standard map

From Scholarpedia
Boris Chirikov and Dima Shepelyansky (2008), Scholarpedia, 3(3):3550. doi:10.4249/scholarpedia.3550 revision #183629 [link to/cite this article]
Jump to: navigation, search
Post-publication activity

Curator: Dima Shepelyansky

The Chirikov standard map [1], [2] is an area-preserving map for two canonical dynamical variables, i.e., momentum and coordinate \( (p,x)\). It is described by the equations:

\[\tag{1} \begin{array}{lcr} \bar{p} = p+K\sin x \\ \bar{x} = x+\bar{p} \end{array} \]

where the bars indicate the new values of variables after one map iteration and \(K\) is a dimensionless parameter that influences the degree of chaos. Due to the periodicity of \( \sin x \) the dynamics can be considered on a cylinder (by taking \( x \!\!\! \mod{2\pi} \)) or on a torus (by taking both \( x,p \!\!\! \mod{2\pi} \)). The map is generated by the time dependent Hamiltonian \( H(p,x,t)= p^2/2 +K \cos(x) \, \delta_1(t) \), where \( \delta_1(t) \) is a periodic \( \delta- \)function with period 1 in time. The dynamics is given by a sequence of free propagations interleaved with periodic kicks.

Examples of the Poincare sections of the standard map on a torus are shown in the following Figs. 1,2,3.

Figure 1: K=0.5
Figure 2: K=0.971635
Figure 3: K=5

Below the critical parameter \( K < K_c \) (Fig.1) the invariant Kolmogorov-Arnold-Moser (KAM) curves restrict the variation of momentum \( p \) to be bounded. The golden KAM curve with the rotation number \[ r=r_g=(\sqrt{5}-1)/2 =0.618033... \] is destroyed at \(K=K_g=0.971635... \) [3], [4] (Fig.2). This Fig. shows a generic phase space structure typical for various area-preserving maps with smooth generating functions: stability islands are embedded in a chaotic sea, similar structure appears on smaller and smaller scales. In a vicinity of a critical invariant curve, with a golden tail in a continued fraction expansion of \( r\), the phase space structure is universal for all smooth maps [4]. Above the critical value \( K > K_c \) (see Fig.3 showing a chaotic component and visible islands of stability) the variation of \(p\) becomes unbounded and is characterized by a diffusive growth \( p^2 \sim D_0 t\) with number of map iterations \( t \). Here \( D_0 \) is a diffusion rate with \( D_0 \approx (K-K_c)^3/3 \) for \( K_c < K < 4 \) and \( D_0 \approx D_{ql}=K^2/2 \) for \( 4 < K \) [2], [5]. There are strong arguments in favor of the equality \( K_c = K_g \) but rigorously it is only proven that there are no KAM curves for \( K > 63/64 = 0.984375 \) [6]. With the numerical results [3], [4] this implies inequality for the global chaos border, \( K_g \leq K_c < 63/64 \).

A simple analytical criterion proposed in 1959 and now known as the Chirikov resonance-overlap criterion [7] gives the chaos border \( K_c = \pi^2/4\) [1] and after some improvements leads to \( K_c \approx 1.2\) [2],[8]. This accuracy is not so impressive compared to modern numerical methods but still up to now this criterion remains the only simple analytical tool for determining the chaos border in various Hamiltonian dynamical systems.

The Kolmogorov-Sinai entropy of the map is well described by relation \( h \approx \ln(K/2)\) valid for \( K > 4 \)[1], [2].


Universality and Applications

The map (1) describes a situation when nonlinear resonances are equidistant in phase space that corresponds to a local description of dynamical chaos. Due to this property various dynamical systems and maps can be locally reduced to the standard map and due to this reason the term standard map was coined in [2]. Thus, the standard map describes a universal, generic behavior of area-preserving maps with divided phase space when integrable islands of stability are surrounded by a chaotic component. A short list of systems reducible to the standard map is given below:

  • chaotic layer around separatrix of a nonlinear resonance induced by a monochromatic force (the whisker map) [2]
  • charged particle confinement in mirror magnetic traps [1], [2], [7], [9]
  • fast crossing of nonlinear resonance [1], [10]
  • particle dynamics in accelerators [11]
  • comet dynamics in solar system [12] with a rather similar map for the comet Halley [13]
  • microwave ionization of Rydberg atoms (linked to the Kepler map) [14] and autoionization of molecular Rydberg states [15]
  • electron magnetotransport in a resonant tunneling diode [16]

Open Problems

  • In spite of fundamental advances in ergodic theory [17], a rigorous proof of the existence of a set of positive measure of orbits with positive entropy is still missing, even for specific values of \( K \) (see e.g. [18]).
  • What are the fractal properties of critical chaos parameter \( K_c(r) \) as a function of arithmetic properties of the rotation number \( r \) of KAM curve? do local maxima correspond only to a golden tail of continuous fraction expansion [3], [4] or they may have tails with Markov numbers as it is conjectured in [19]? (see also [20])
  • Due to trajectory sticking around stability islands the statistics of Poincare recurrences in Hamiltonian systems with divided phase space (see e.g. Fig.2 with a critical golden KAM curve) is characterized by an algebraic decay \( P(\tau) \propto 1/\tau^\alpha \) with \( \alpha \approx 1.5 \) while a theory based on the universality in a vicinity of critical golden curve gives \( \alpha \approx 3 \ ;\) this difference persists up to 1013 map iterations; as a result correlation functions decay rather slowly \( C(\tau) \sim \tau P(\tau) \propto 1/\tau^{\alpha-1} \) that can lead to a divergence of diffusion rate \( D \sim \tau C(\tau) \) (see [21] and Refs. therein)

Quantum Map

Figure 4: Dependence of rescaled rotator energy \( E/(k^2/4) \) on time \( t \) for \( K=kT=5, \hbar=0.25 (k=20, T=0.25) \); the full curve shows numerical data and the straight line gives the diffusive energy growth in the classical case (from [23]).

The quantization of the standard map is obtained by considering variables in (1) as the Heisenberg operators with the commutation relation \( [p,x] = -i \hbar \), where \( \hbar \) is an effective dimensionless Planck constant. In a same way it is possible to use the Schrödinger equation with the Hamiltonian \( H(\hat{p},\hat{x},t) \) given above and \( \hat{p}=-i\hbar \partial/\partial x \). Integration on one period gives the quantum map for the wave function \( \psi \):

\[\tag{2} \bar{\psi} = \hat{U} \psi = e^{-i{\hat p}^2/2\hbar} e^{-i K/\hbar \cos {\hat x}} \psi \]

where bar marks the new value of \( \psi \) after one map iteration. Due to space periodicity of the Hamiltonian the momentum can be presented in the form \( p=\hbar (n + \beta) \), where \( n \) is an integer and \( \beta \) is a quasimomentum preserved by the evolution operator \( \hat U \) . The case with \( \beta =0 \) corresponds to a periodic boundary conditions with \( \psi(x+2\pi) =\psi(x) \) and is known as the kicked rotator introduced in [22].

Other notations with \( \hbar \rightarrow T\), \( K/\hbar \rightarrow k\) are also used to mark the dependence on the period \( T\) between kicks, then \( K = k T\). The diffusion rate over quantum levels \( n \) is \( D=D_0/\hbar^2= n^2/t \approx K^2/2\hbar^2 =k^2/2\), thus the rotator energy \( E = <n^2>/2\) grows linearly with time. Quantum interference effects lead to a suppression of this semiclassical diffusion [22] on the diffusive time scale \( t_D \) so that the quantum probability spreads effectively only on a finite number of states \( \Delta n \sim \sqrt{D t_D} \) (Fig.4). According to the analytical estimates obtained in [23]:

\[\tag{3} t_D \sim \Delta n \sim D \sim k^2 \sim D_0/\hbar^2 . \]

This diffusive time scale is much larger than the Ehrenfest time scale [23], [24] \( t_E \sim \ln(1/\hbar)/2h \) after which a minimal coherent wave packet spreads over the whole phase space due to the exponential instability of classical dynamics. For \( t < t_E \) a quantum wave packet follows the chaotic dynamics of a classical trajectory as it is guaranteed by the Ehrenfest theorem [23]. For the case of Fig.4 the Kolmogorov-Sinai entropy \( h \approx 1 \) and the Ehrenfest time \( t_E \sim 1 \) is extremely short comparing to the diffusive time \( t_D \sim D \sim 200 \). The quantum suppression of chaotic diffusion is similar to the Anderson localization in disordered systems if to consider the level number as an effective site number in a disordered lattice, such an analogy has been established in [25]. However, in contrast to a disordered potential for the case of Anderson localization, in the quantum map (2) diffusion and chaos have a pure deterministic origin appearing as a result of dynamical chaos in the classical limit.

Figure 5: Dependence of the localization length \( l \) on the quantum parameter of chaos \( K \rightarrow K_q=2k \sin T/2 \). The circles and the curve are, respectively, the numerical data and the theory for the classical diffusion \( D(K) \) (see [8]). The quantum data for \( l \) are shown by \( + \) (for \( 0<T<\pi \)) and by \( \times \) (for \( \pi<T<2\pi \)); here \( k=30; D_{ql}=k^2/2 \) (from [27]).

Due to that this phenomenon is called the dynamical localization. The eigenstates of the unitary evolution operator \( \hat U \) are exponentially localized over momentum states \( \psi_m(n) \sim \exp(-|n-m|/l)/\sqrt{l} \) with the localization length \( l \sim \Delta n \sim t_D \) given by the relation [26], [27] \[\tag{4} l=D(K)/2 =D_0(K)/2\hbar^2, \]

where \( D \) is the semiclassical diffusion expressed via a square number of levels per period of perturbation. For \( \hbar = T > 1 \) the chaos parameter \( K \) in the dependence \( D(K) \) should be replaced by its quantum value \( K \rightarrow K_q = 2 k \sin T/2 \) [27]. The quantum localization length \( l \) repeats the characteristic oscillations of the classical diffusion as it is shown in Fig.5. The relation (4) assumes that \( T/4\pi \) is a typical irrational number while for rational values of this ratio the phenomenon of quantum resonance takes place and the energy grows quadratically with time for rational values of quasimomentum [28]. The derivations of the relation (4) based on the field theory methods applied to dynamical systems with chaotic diffusion can be find in [29], [30] (see also Refs. therein).

If the quantum map (2) is taken on a torus with \( N \) levels then the level spacing statistics is described by the Poisson law for \( N \gg l \) and by the Wigner-Dyson law of the random matrix theory for \( N \ll l \) [24],[31]. In the later case the quantum eigenstates are ergodic on a torus in agreement with the Shnirelman theorem and the level spacing statistics agrees with the Bohigas-Giannoni-Schmit conjecture (see books on quantum chaos in Recommended Reading).

The quantum map (2) was built up experimentally with cold atoms in a kicked optical lattice by the group of M.Raizen [32]. Such a case corresponds to a particle in an infinite periodic lattice with averaging over many various \( \beta \). The quantum resonances at \( \beta \approx 0 \) were also experimentally observed with the Bose-Einstein condensate (BEC) in [33]. Quantum accelerator modes for kicked atoms falling in the gravitational field were found and analyzed in [34].

Extensions and Related Quantum Systems

Due to universal properties of the standard map its quantum version also finds applications for various systems and various physical effects:

  • dynamical localization for ionization of excited hydrogen atoms in a microwave field

was theoretically predicted in [35] and was experimentally observed by the group of P.Koch [36] (see more details in [14],[37],[38])

  • quantum particle in a triangular well and monochromatic field with a quantum delocalization transition [39]
  • the kicked Harper model where in contrast to the relation (4) the quantum delocalization can take place due to quasi-periodicity of unperturbed spectrum (see [40], [41] and Refs. therein)
  • 3D Anderson transition in kicked rotator with modulated kick strength and quantum transport in mesoscopic conductors (see [42] and Refs. therein)
  • dissipative quantum chaos [43]
  • fractal Weyl law for the quantum standard map with absorption (see [44] and Refs. therein)
Figure 6: Dependence of rescaled energy \( E/(k^2/4) \) on time in the classical map (1) at \( K=5 \); time reversal is performed at \( t=150 \); numerical simulations are done on BESM-6 with relative accuracy \( \epsilon \approx 10^{-12} \) (from [46]).

Time Reversibility and Boltzmann - Loschmidt Dispute

Figure 7: Same as in Fig.6 but for the quantum map (2) with \( K=5, \hbar=0.25 \), the straight line shows the classical diffusion; time reversal is performed at the moment \( t=150 \) marked by the vertical line, numerical simulations are done on the same computer BESM-6, in addition random quantum phases \( 0<\Delta \phi <0.1 \) are added for quantum amplitudes in momentum representation at the moment of time reversal (from [46]).

The statistical theory of gases developed by Boltzmann leads to macroscopic irreversibility and entropy growth even if dynamical equations of motion are time reversible. This contradiction was pointed out by Loschmidt and is now known as the Loschmidt paradox. The reply of Boltzmann relied on the technical difficulty of velocity reversal for material particles: a story tells that he simply said "then go and do it" [45]. The modern resolution of this famous dispute, which took place around 1876 in Wien, came with the development of the theory of dynamical chaos (see e.g. [8], [17]). Indeed, for chaotic dynamics the Kolmogorov-Sinai entropy is positive and small perturbations grow exponentially with time, making the motion practically irreversible. This fact is convenient to illustrate on the example of the standard map which dynamics is time reversible, e.g. by inverting all velocities at the middle of free propagation between two kicks (see Fig.6). This explanation is valid for classical dynamics, while the case of quantum dynamics requires special consideration. Indeed, in the quantum case the exponential growth takes place only during the rather short Ehrenfest time, and the quantum evolution remains stable and reversible in presence of small perturbations [46] (see Fig.7). Quantum reversibility in presence of various perturbations has been actively studied in recent years and is now described through the Loschmidt echo (see [47] and Refs. therein). A method of approximate time reversal of matter waves for ultracold atoms in the regime of quantum chaos, like those in [32], [33], is proposed in [48]. In this method a large fraction of the atoms returns back even if the time reversal is not perfect. This fraction of the atoms exhibits Loschmidt cooling which can decrease their temperature by several orders of magnitude. At the same time a kicked BEC of attractive atoms (soliton) described by the Gross-Pitaevskii equation demonstrates a truly chaotic dynamics for which the exponential instability breaks the time reversibility [49]. However, since a number of atoms in BEC is finite and since BEC is a really quantum object one should expect that the Ehrenfest time is still very short and hence the time reversibility should be preserved in presence of small errors if the second quantization is taken into account.

Links to Other Physical Topics

Frenkel-Kontorova Model

The Frenkel-Kontorova model describes a one-dimensional chain of atoms/particles with harmonic couplings placed in a periodic potential [50]. This model was introduced with the aim to study crystal dislocations but it also successfully applies for the description of commensurate-incommensurate phase transitions, epitaxial monolayers on the crystal surface, ionic conductors, glassy materials, charge-density waves and dry friction [51]. The Hamiltonian of the model is \( H= \sum_i \left({P_i^2 \over 2} + {(x_i -x_{i-1})^2 \over 2}- K \cos x_i \right)\), where \( P_i, x_i \) are momentum and position of atom \( i \). At the equilibrium the momenta \( P_i =0\) and \( \partial H/\partial x_i =0\) so that the positions of atoms are described by the map (1) with \( p_{i+1} = x_{i+1}- x_i , \; p_{i+1}= p_i +K\sin x_i\). The density of atoms corresponds to the rotation number \( r \) of an invariant KAM curve. For the golden density with \( r =r_g\) the chain slides in the periodic potential for \( K < K_g \) (KAM curve regime) while for \( K > K_g \) the transition by the breaking of analyticity, or Aubry transition, takes place, the chain becomes pinned and atoms form an invariant Cantor set called cantorus (see [52] and Aubry-Mather theory). In this regime the phonon spectrum has a gap so that the phonon excitations are suppressed at low temperature. The mathematical Aubry-Mather theory guarantees that the ground state of the chain exists and is unique. However there exist exponentially many static equilibrium configurations which are exponentially close to the energy of the ground state. The energies of these configurations form a fractal quasi-degenerate band structure and become mixed at any physically realistic temperature. Thus, such configurations can be viewed as a dynamical spin glass. For a case of Coulomb interactions between particles (e.g. ions or electrons) one obtains a problem of Wigner crystal in a periodic potential which again is locally described by the Frenkel-Kontorova model since the map (1) gives the local description of the dynamics. For the quantum Frenkel-Kontorova model the dynamics of atoms (ions) in the chain is quantum. In this case the quantum vacuum fluctuations and instanton tunneling lead to a quantum melting of pinned phase: above a certain effective Planck constant a quantum phase transition takes place from pinned instanton glass to sliding phonon gas (see [53] and Refs. therein).

Quantum Computing

One iteration of maps (1) and (2) can be simulated on a quantum computer in a polynomial number of quantum gates for an exponentially large vector representing a Liouville density distribution or a quantum state. The quantum algorithm of such a quantum computation is described in [54], effects of quantum errors are analyzed in [55] (see also Refs. therein).

Historical Notes

The standard map (1) in a form of recursive relation for atoms in a periodic potential appears already in the works of Kontorova and Frenkel [50]. As a dynamical map it first appeared as a description of electron dynamics in a new relativistic accelerator proposed by V.I.Veksler (Dokl. Akad. Nauk SSSR 43: 346 (1944)). The regime of a stable regular acceleration was studied later also by A.A.Kolomensky (Zh. Tekh. Fiz. 30: 1347 (1960)) and S.P.Kapitsa, V.N.Melekhin ("Microtron", Nauka, Moscow (1969) in Russian). Among the early researchers of model (1) was also British physicist J.B.Taylor (unpublished reports). The description of chaos in map (1) and its main properties, including chaos border, diffusion rate and positive entropy, was given in [1]. The term "standard map" appeared in [2], "Chirikov-Taylor map" [8] and "Chirikov standard map" [16] are also used, the quantum standard map or kicked rotator was first considered in [22]. Appearance of other terms: Kolmogorov-Arnold-Moser theory [1], Arnold diffusion [1], Kolmogorov-Sinai entropy [2], Ehrenfest time [24].

Recommended Reading

B.V.Chirikov, "Research concerning the theory of nonlinear resonance and stochasticity", Preprint N 267, Institute of Nuclear Physics, Novosibirsk (1969), (Engl. Trans., CERN Trans. 71-40 (1971))

B.V.Chirikov, "A universal instability of many-dimensional oscillator systems", Phys. Rep. 52: 263 (1979).

B.V.Chirikov, "Time-dependent quantum systems" in "Chaos and quantum mechanics", Les Houches Lecture Series, Vol. 52, pp.443-545, Eds. M.-J.Giannoni, A.Voros, J.Zinn-Justin, Elsevier Sci. Publ., Amsterdam (1991)

A.J.Lichtenberg, M.A.Lieberman, "Regular and chaotic dynamics", Springer, Berlin (1992).

F.Haake, "Quantum signatures of chaos", Springer, Berlin (2001).

L.E.Reichl, "The Transition to chaos in conservative classical systems and quantum manifestations", Springer, Berlin (2004).

Internal references

  • Martin Gutzwiller (2007) Quantum chaos. Scholarpedia, 2(12):3146.

External Links

Selected publications of Boris Chirikov [1]

Sputnik of Chaos [2]

Google query for "standard map" [3]



B.V.Chirikov, "Research concerning the theory of nonlinear resonance and stochasticity", Preprint N 267, Institute of Nuclear Physics, Novosibirsk (1969) [4], (Engl. Trans., CERN Trans. 71-40 (1971)) [5].


B.V.Chirikov, "A universal instability of many-dimensional oscillator systems", Phys. Rep. 52: 263 (1979) [6].


J.M.Greene, "Method for determining a stochastic transition", J. Math. Phys. 20(6): 1183 (1979).


R.S.MacKay, "A renormalization approach to invariant circles in area-preserving maps", Physica D 7(1-3): 283 (1983).


R.S.MacKay, J.D.Meiss, I.C.Percival, "Transport in Hamiltonian systems", Physica D 13(1-2): 55 (1984).


R.S.MacKay, I.C.Percival, "Converse KAM - theory and practice", Comm. Math. Phys. 94(4): 469 (1985).


B.V.Chirikov, "Resonance processes in magnetic traps", At. Energ. 6: 630 (1959) (in Russian [7]) (Engl. Transl., J. Nucl. Energy Part C: Plasma Phys. 1: 253 (1960) [8]).


A.J.Lichtenberg, M.A.Lieberman, "Regular and chaotic dynamics", Springer, Berlin (1992).


B.V.Chirikov, "Particle confinement and adiabatic invariance", Proc. R. Soc. Lond. A 413: 145 (1987) [9].


B.V.Chirikov, D.L.Shepelyanskii, "Diffusion during multiple passage through a nonlinear resonance", Sov. Phys. Tech. Phys. 27(2): 156 (1982) [10] (in Russian [11])


F.M.Izraelev, "Nearly linear mappings and their applications", Physica D 1(3): 243 (1980).


T.Y.Petrowsky, "Chaos and cometary clouds in the solar system", Phys. Lett. A 117(7): 328 (1986).


B.V.Chirikov, V.V.Vecheslavov, "Chaotic dynamics of comet Halley", Astron. Astrophys. 221: 146 (1989) [12].


G.Casati, I.Guarneri, D.L.Shepelyansky, "Hydrogen atom in monochromatic field: chaos and dynamical photonic localization", IEEE J. of Quant. Elect. 24: 1420 (1988).


F.Benvenuto, G.Casati, D.L.Shepelyansky, "Chaotic autoionization of molecular Rydberg states", Phys. Rev. Lett. 72: 1818 (1994).


D.L.Shepelyansky, A.D.Stone, "Chaotic Landau level mixing in classical and quantum wells", Phys. Rev. Lett. 74: 2098 (1995).


I.P.Cornfeld, S.V.Fomin, Ya.G.Sinai, "Ergodic theory", Springer, Berlin (1982).


A.Giorgilli, V.F.Lazutkin, "Some remarks on the problem of ergodicity of the standard map", Phys. Lett. A 272: 359 (2000).


B.V.Chirikov, D.L.Shepelyansky, "Chaos border and statistical anomalies", Eds. D.V.Shirkov, D.I.Kazakov and A.A.Vladimirov, World Sci. Publ., Singapore, "Renormalization Group" p.221 (1988) [13].


J.M.Greene, R.S.MacKay, J.Stark, "Boundary circles for area-preserving maps", Physica D 21(2-3): 267 (1986).


B.V.Chirikov, D.L.Shepelyansky, "Asymptotic statistics of Poincare recurrences in Hamiltonian systems with Divided Phase Space", Phys. Rev. Lett. 82: 528 (1999) [14]; 89: 239402 (2002) [15] .


F.M.Izrailev, G.Casati, J.Ford, B.V.Chirikov, "Stochastic behavior of a quantum pendulum under a periodic perturbation", Preprint 78-46, Institute of Nuclear Physics, Novosibirsk (1978) (extended version in Russian) [16]; G.Casati, B.V.Chirikov, F.M.Izrailev, J.Ford, Lecture Notes in Physics, Springer, Berlin, 93: 334 (1979) [17].


B.V.Chirikov, F.M.Izrailev, D.L.Shepelyansky, "Dynamical stochasticity in classical and quantum mechanics", Sov. Scient. Rev. C 2: 209 (1981) (Section C - Mathematical Physics Reviews, Ed. S.P.Novikov vol.2, Harwood Acad. Publ., Chur, Switzerland (1981)) [18].


B.V.Chirikov, F.M.Izrailev, D.L.Shepelyansky, "Quantum chaos: localization vs. ergodicity", Physica D 33: 77 (1988) [19].


S.Fishman, D.R.Grempel, R.E.Prange, "Chaos, quantum recurrences, and Anderson localization", Phys. Rev. Lett. 49: 509 (1982).


B.V.Chirikov, D.L.Shepelyanskii, "Localization of dynamical chaos in quantum systems", Izv. Vyssh. Ucheb. Zaved. Radiofizika 29(9): 1041 (1986) (in Russian [20]); (English Trans. Plenum Publ. [21] ).


D.L.Shepelyansky, "Localization of diffusive excitation in multi-level systems", Physica D 28: 103 (1987).


F.M.Izrailev, D.L.Shepelyanskii, "Quantum resonance for a rotator in a nonlinear periodic field", Theor. Math. Phys. 43(3): 553 (1980); see also I. Dana and D.L. Dorofeev, "General quantum resonances of kicked particle", Phys. Rev. E 73: 026206 (2006).


K.M.Frahm, "Localization in a rough billiard: a sigma model formulation", Phys. Rev. B 55: 8626(R) (1997).


C.Tian, A.Kamenev, A.Larkin, "Weak dynamical localization in periodically kicked cold atomic gases", Phys. Rev. Lett. 93: 124101 (2004).


F.M.Izrailev, "Simple models of quantum chaos: spectrum and eigenfunctions", Phys. Rep. 196: 299 (1990).


F.L.Moore, J.C.Robinson, C.F.Bharucha, B.Sundaram, M.G.Raizen, "Atom optics realization of the quantum \( \delta \)-kicked rotor", Phys. Rev. Lett. 75: 4598 (1995).


C.Ryu, M.F.Anderen, A.Vaziri, M.B.d'Arcy, J.M.Grossman, K.Helmerson, W.D.Phillips, "High-order quantum resonances observed in a periodically kicked Bose-Einstein condensate", Phys. Rev. Lett. 96: 160403 (2006).


A.Buchleitner, M.B.d'Arcy, S.Fishman, S.A.Gardiner, I.Guarneri, Z.-Y.Ma, L.Rebuzzini, G.S.Summy, "Quantum accelerator modes from the Farey tree", Phys. Rev. Lett. 96: 164101 (2006).


G.Casati, B.V.Chirikov, D.L.Shepelyansky, ""Quantum limitations for chaotic excitation of the hydrogen atom in a monochromatic field", Phys. Rev. Lett. 53: 2525 (1984) [22].


E.J.Galvez, B.E.Sauer, L.Moorman, P.M.Koch, D.Richards, "Microwave ionization of H atoms: breakdown of classical dynamics for high frequencies", Phys. Rev. Lett. 61: 2011 (1988).


G.Casati, B.V.Chirikov, D.L.Shepelyansky, I.Guarneri, ""Relevance of classical chaos in quantum mechanics: the hydrogen atom in a monochromatic field", Phys. Rep. 154: 77 (1987) [23].


P.M.Koch, K.A.H. van Leeuwen, "The importance of resonances in microwave “ionization” of excited hydrogen atoms", Phys. Rep. 255: 289 (1995).


F.Benvenuto, G.Casati, I.Guarneri, D.L.Shepelyansky, "A quantum transition from localized to extended states in a classically chaotic system", Z.Phys.B - Cond. Matt. 84: 159 (1991).


R.Lima, D.L.Shepelyansky, "Fast delocalization in a model of quantum kicked rotator", Phys. Rev. Lett. 67: 1377 (1991).


T.Prosen, I.I.Saija, N.Shah, "Dimer decimation and intricately nested localized-ballistic phases of a kicked Harper model", Phys. Rev. Lett. 87: 066601 (2001).


F.Borgonovi, D.L.Shepelyansky, "Two interacting particles in an effective 2-3-d random potential", J. de Physique I France 6: 287 (1996) and F. Borgonovi, I.Guarneri and L. Rebuzzini, "Chaotic diffusion and statistics of universal scattering fluctuations", Phys. Rev. Lett. 72: 1463 (1994).


G.G.Carlo, G.Benenti, D.L.Shepelyansky, "Dissipative quantum chaos: transition from wave packet collapse to explosion", Phys. Rev. Lett. 95: 164101 (2005).


D.L.Shepelyansky, "Fractal Weyl law for quantum fractal eigenstates", Phys. Rev. E 77: 015202(R) (2008).


J.E.Mayer, M.Goppert-Mayer, "Statistical mechanics", John Wiley & Sons, N.Y. (1977).


D.L.Shepelyansky, "Some statistical properties of simple classically stochastic quantum systems", Physica D 8: 208 (1983).


T.Gorin, T.Prosen, T.H.Seligman, M.Znidaric, "Dynamics of Loschmidt echos and fidelity decay", Phys. Rep. 435: 33 (2006).


J.Martin, B.Georgeot, D.L.Shepelyansky, "Loschmidt cooling by time reversal of atomic matter waves", arxiv:0710.4860[cond-mat]), Phys. Rev. Lett. 100: 044106 (2008).


F.Benvenuto, G.Casati, A.S.Pikovsky, D.L.Shepelyansky, "Manifestations of classical and quantum chaos in nonlinear wave propagation", Phys. Rev. A 44: 3423(R) (1991).


T.A.Kontorova, Ya.I.Frenkel, "On the theory of plastic deformation and doubling", Zh. Eksp. Teor. Fiz 8: 89 (1938); 8: 1340 (1938); 8: 1359 (1938) (in Russian).


O.M. Braun, Yu.S. Kivshar, "The Frenkel-Kontorova Model: Concepts, Methods, and Applications", Springer, Berlin (2004).


S.Aubry, "The twist map, the extended Frenkel-Kontorova model and the devil's staircase", Physica D 7(1-3): 240 (1983).


I.Garcia-Mata, O.V.Zhirov, D.L.Shepelyansky, "Frenkel-Kontorova model with cold trapped ions", Eur. Phys. J. D 41: 325 (2007).


B.Georgeot, D.L.Shepelyansky, "Exponential gain in quantum computing of quantum chaos and localization", Phys. Rev. Lett. 86: 2890 (2001).


K.M.Frahm, R.Fleckinger, D.L.Shepelyansky, "Quantum chaos and random matrix theory for fidelity decay in quantum computations with static imperfections", Eur. Phys. J. D 29: 139 (2004).

See also

Hamiltonian systems, Mapping, Chaos, Kolmogorov-Arnold-Moser Theory, Kolmogorov-Sinai entropy, Aubry-Mather theory, Quantum chaos

Updates added after 2008

Here some additional features of the standard map, which appeared after 2008 or were omitted in the 2008 article edition, are presented.

Mathematical aspects

Mathematical results are presented here for the classical (UM1-UM3) and quantum (UM4) standard map.

  • UM1) A rigorous proof is given in the standard map for the existence of chaotic trajectories with unbounded momenta for large enough coupling constant \( K > K_0 \), where \( K_0 \) depends on a coding representation of a trajectory. The obtained chaotic trajectories correspond to stationary configurations of the Frenkel-Kontorowa model with a finite (non-zero) photon gap. The concept of anti-integrability emerges from the theorems presented in [UM1].
  • UM2) The large basic sets, which fill in the torus as the parameter runs to infinity, are constructed. It is proven that, for a residual set of large parameters, these basic sets accumulated by elliptic periodic islands. It is shown there exists a \( K_0 > 0 \) and a dense set of parameters \( K_0 \leq K < \infty \) for which the standard map exhibits homoclinic tangencies [UM2].
  • UM3) It is proven that stochastic sea of the standard map has full Hausdorff dimension

for sufficiently large topologically generic parameters [UM3].

  • UM4) For the quantum standard map with a generic quadratic rotational spectrum the localization is proven for small kick amplitudes [UM4].

Physical aspects and numerical results

Results for the Ulam method for the standard map are presented in [UP1], Poincare recurrences in [UP2] and the fractal Weyl law for Perron-Frobenius operators in [UP3]. The results of the field theory for the quantum standard map are discussed in [UP4].

Figure 8: Fig. 8U. Amplitude of the eigenstate of the Ulam approximate of Perron-Frobenius operator of the standard map at \( K=0.971635406\); the number of cells and matrix size are \( N_d = 127282\) and eigenvalue is \( \lambda_2 = 0.99878108 \); amplitude is proportional to color with maximum for red and zero for blue; upper part of phase plane is shown for the range \( 0 \leq x/2\pi \leq 1; 0 \leq p/2\pi \leq 0.5 \) .(from [UP1b])
Figure 9: Fig. 9U. The complex plane of spectrum of eigenvalues \( \lambda \) of UFPO \( S\) of the standard map at \( K=0.971635406\); the number of cells and matrix size are \( N_d = 16609\); right panel shows zoom of the global spectrum of left panel. (from [UP1b])
  • UP1) In 1960 Ulam proposed a method, known now as the Ulam method, for construction a finite size matrix approximate for the Perron-Frobenius operator of a dynamical system in a continuous phase space [UP1a]. The method allows to construct numerically a matrix of Markov transitions between cells in a discretized phase space with fully chaotic dynamics. The method is knwon to be converging to the continuous limit of Perron-Frobenius operator when the phase space is fully chaotic. However, for systems with a divided phase space an effective noise induced by a finite cell size breaks the convergence leading to a destruction of the invariant KAM curve. This problem was resolved in [UP1b] by a generalized Ulam method in which the Markov transitions between cells are generated by one chaotic trajectory starting inside a chaotic component. The extensive numerical studies based on the Arnoldi method show that the Ulam approximate of the Perron-Frobenius operator \( S \) (UPFO) on a chaotic component converges to a continuous limit. Typically, in this regime the spectrum of relaxation modes is characterized by a power law decay for small relaxation rates. The numerical results show that the exponent of this decay is approximately equal to the exponent of Poincare recurrences in such systems. The eigenmodes, or eigenstates, show links with trajectories sticking around stability islands. An example of such an eigenstate is shown in Fig.8U. The spectrum of UFPO \( S \) is shown in Fig.9U.
Figure 10: Fig. 10U. The left panel shows the statistics of Poincare recurrences \( P(t) \) averaged over 10 random realizations and obtained by the surviaval Manonte Carlo method proposed in [UP2a] operating here with \( N_f =1000 \) surviving trajectories; magenta and red curves show independence of results on accuracy (overlaped curves), green curve shows data at \( N_f =100 \) with larger fluctuations on large times \( t \); blue curve shows data from [UP2b] at shorter times; the dashed line shows the slope \( P \propto 1/t^{1.5}\). The right panel shows the results of red curve of left panel with the results obtained with the generalized Ulam method with the number of cells \( M^2/2\) in the upper half plane of Fig.8U. (from [UP2a])
  • UP2) The numerical studies of the Poincare recurrences in the standard map with the critical golden curve have been performed in [UP2a] with a new survival Monte Carlo method which allows to study recurrences on times changing by ten orders of magnitude (see Fig.10U). The comparison is done with the results of generalized Ulam method and localization properties of eigenstates of the Ulam matrix are analyzed. The recurrences at long times are determined by trajectory sticking in a vicinity of the critical golden curve and secondary resonance structures. On the investigated scales the Poincare decay exponent, in <math P \propto 1/t^{\alpha} </math> is found to be approximately </math> \alpha \approx 1.58 </math> in a satisfactory agreement with early and more recent studies of various symplectic 2D maps [19], [20], [UP2c], [UP2d] thus indicating the universality of the exponent. The detailed theoretical explanation of the algebraic decay of Poincare recurrences and the exponent value is still lucking. It is interesting to add a few notes: A) in the standard map with dissipation (e.g. the right hand side of upper line Eq.(1) is multiplied by a coefficient less than unity) in the regime of a strange attractor the decay of Poincare recurrences is exponential [UP2c]; B) in the symplectic case additional noise leads to a diffusive type decay \( P \propto 1/\sqrt{t} \) on a long time scale due to diffusion inside stability islands, when they are present; C) in the quantum standard map the quantum Poincare recurrences are characterized by a decay \( P \propto 1/t \) due to quantum tunneling inside the stability islands [UP2e].
  • UP3) For the standard map with absorption or dissipation in a chaotic regime it is show that the Ulam approximate of Perron-Frobenius operator is characterized by the fractal Weyl law with the exponent given by a half of the fractal dimension of related chaotic repeller or strange attractor [UP3] (see also the section Ulam networks in Google matrix). The fractal Weyl law for the quantum standard map with absorption is analyzed in [44].
  • UP4) The field theory methods for the quantum standard map have been developed in [UP4]. It is shown that the effective theory describing the long wave length physics of the system is precisely the supersymmetric nonlinear sigma-model for quasi one-dimensional metallic wires. It is shown that the localization length is given by Eq.(4). This proves that the analogy between chaotic systems with dynamical localization and disordered metals can indeed be exact, as claimed by the authors. However, this approach misses certain properties of quantum evolution, thus it gives the finite localization length for the quantum kicked Harper model in the chaotic regime while the numerical results show the existence of delocalized quantum phase.

Related models and systems

  • UR1) The properties of quantum kicked rotator with rotational phases depending linearly on level number in Eq.(2) and generalized to any number of dimensions are considered in [UR1a],[UR1b], [UR1c]. The mathematical prove of localization of all eigenstates is given for small [UR1b] and arbitrary kick amplitudes [UR1c]. This result is rather clear from the view point of classical dynamics where linear dependence of Hamiltonian on actions (linear spectrum) leads to complete integrability of motion.

References added after 2008


S.Aubry, G.Abramovici, "Chaotic trajectories in the standard map. The concept of anti-integrability", Physica D 43: 199 (1990).


P.Duarte, "Plenty of elliptic islands for the standard family of area preserving maps", Annales de l'I.H.P. Sec. C 11(4): 359 (1994).


A.Gorodentski, "On stochastic sea of the standard map", Comm. Math. Phys. 309(1): 155 (2012).


J.Bourgain, "Estimates on Green's functions, localization and the quantum kicked rotor model", Ann. Math. 156: 249 (2002).


S.M.Ulam, "A Collection of mathematical problems, Interscience tracs in pure and applied mathematics", Interscience, New York 8: 73 (1960).


K.M.Frahm, D.L.Shepelyansky, "Ulam method for the Chirikov standard map", Eur. Phys. J. B 76: 57 (2010).


K.M.Frahm, D.L.Shepelyansky, "Poincare recurrences and Ulam method for the Chirikov standard map", Eur. Phys. J. B 86: 322 (2013).


M.Weissm L.Hufnagel, R.Ketzmerick, "Universal power-law decay in Hamiltonian systems?", Phys. Rev. Lett. 100: 184101 (2008).


B.V.Chirikov, D.L.Shepelyansky, "Correlation properties of dynamical chaos in Hamiltonian systems", Physica D 13: 395 (1984), [24], detailed BINP preprint in RUssian [25]


G.Cristadoro, R.Ketzmerick, "Universality of algebraic decays in Hamiltonian systems", Phys. Rev. Lett. 100: 184191 (2008).


G.Casati, G.Maspero, D.L.Shepelyansky, "Quantum poincare recurrences", Phys. Rev. Lett. 82: 524 (1999).


L.Ermann, D.L.Shepelyansky, "Ulam method and fractal Weyl law for Perron-Frobenius operators", Eur. J. B 75: 299 (2010).


A.Altland, M.Zinbauer, "Field theory of the quantum kicked rotor", Phys. Rev. Lett. 77: 4536 (1996).


S.Fishman, D.R.Grenpel, R.E.Prange, "Localization in an incommensurate potential: an exactly solvable model", Phys. Rev. Lett. 49: 833 (1982).


J.Bellissard, R.Lima, E.Scoppola, "Localization in v-dimensional incommensurate structures", Comm. Math. Phys. 88: 465 (1983).


A.L.Figotin, L.A.Pastur, "An exactly solvable model of a multidimensional incommensurate structure localization for linear spectrum", Comm. Math. Phys. 95: 401 (1984).

See also added after 2008

Anderson localization and quantum chaos maps, Cold atom experiments in quantum chaos, Google matrix, kicked Harper model, Microwave ionization of hydrogen atoms

Personal tools

Focal areas