Chirikov standard map

From Scholarpedia
Boris Chirikov and Dima Shepelyansky (2008), Scholarpedia, 3(3):3550. doi:10.4249/scholarpedia.3550 revision #183857 [link to/cite this article]
(Redirected from Kicked rotor)
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 Raizen group [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.Petrosky, "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).


a) 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 b) 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).


a) 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 b) 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-Kontorova 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). Advanced numerical methods for the standard map are described in (UP4),(UP5). The results of the field theory for the quantum standard map are discussed in (UP6).

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 known 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 survival Monte Carlo method proposed in [UP2a] operating here with \( N_f =1000 \) surviving trajectories; magenta and red curves show independence of results on accuracy (overlapped 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 resonance structures. On the investigated scales the Poincare decay exponent, in \( P \propto 1/t^{\alpha} \) is found to be approximately \( \alpha \approx 1.58 \) 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 [UP2c]; 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 [UP3a] (see also the section Ulam networks in Google matrix). The fractal Weyl law for the quantum standard map with absorption is analyzed in [44]. The semiclassical properties of eigenstates in quantum dissipative systems are analyzed in [UP3b],[UP3c].
  • UP4) For many important Hamiltonian maps (e.g., the standard map) it is possible to construct related mappings that (i) carry a lattice into itself; (ii) approach the original map as the lattice spacing is decreased; (iii) can be iterated exactly using integer arithmetic; and (iv) are Hamiltonian themselves [UP4]. These lattice maps are compared to maps that use floating-point arithmetic to evaluate the original map. The problems associated with roundoff error are analyzed and it is argued that lattice maps are superior to floating-point maps for the study of the long-term behaviour of Hamiltonian dynamical systems.
  • UP5) A powerful visualization method for measure-preserving dynamical systems, based on frequency analysis and Koopman operator theory is applied to the standard map, and other maps, in [UP5].
  • UP6) The field theory methods for the quantum standard map have been developed in [UP6]. 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

Various models related to the standard map are discussed here: localization for the case of linear rotational spectrum (UR1), Shnirelman peak for level spacing statistics (UR2), studies of quantum synchronization (UR3), kicked rotator as a deterministic detector (UR4), effect of two interacting particles for two coupled kicked rotators (UR5), renormalization dynamical chaos for the critical spiral mean in the frequency modulated kicked rotator (UR6), effects of nonlinearity on localization in kicked rotator (UR7), fast Arnold diffusion, chaos measure and Poincare recurrences in coupled standard maps and many-body Hamiltonian systems (UR8).

  • 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 proof 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.
  • UR2) In 1975 Shnirelman proved the theorem about asymptotic multiplicity of Laplace operator [UR2a] which implies that the eigenenergies of generic integrable 2D billiards are exponentially quasidegenrate at large level numbers thus forming pairs of quasidegenerate levels forming the Shnirelman peak in the level spacing statistics. A physical interpretation of the Shnirelman theorem about such bulk quasidegeneracy is given in [UR2b]. Conditions for the strong impact of degeneracy on quantum level statistics are formulated allowing to extend the applications of the Shnirelman theorem to a broad class of quantum systems. It is shown that in some sense the degeneracy between the states connected by time-reversal symmetry is destroyed by tunneling between the future and the past (corresponding to a double well in momentum space). The numerical tests are done with the kicked rotator model of Eq.(2) with the modified potential \( \cos \theta \rightarrow \cos \theta - 0.5 \sin 2\theta \) so that the space symmetry is broken. The generic aspect of the Shnirelman peak is confirmed by the numerical results for rough billiards [UR2c].
Figure 11: Fig. 11U. Quantum synchronization in the standard map with dissipation \( \gamma =0.25 \) and static force \( f \) at \( K=0.8\); panels show the dependence of average momentum \( P \) on \( f\) for dimensionless Planck constant \(\hbar =0.05; 0.5 \) for top; bottom panel. Top panel: synchronization remains stable in respect to quantum fluctuations; bottom panel: quantum fluctuations destroy synchronization. (from [UR3a])
  • UR3) The quantum standard map in infinite space \( x \) (kicked particle) is studied numerically [UR3a] by methods of quantum trajectories in presence of dissipation \( \gamma \) and applied static force \( f \). The model allows to analyze the effects of quantum fluctuations on synchronization and establish the regimes where the synchronization is preserved in a quantum case (quantum synchronization). Thus at small values of dimensionless Planck constant \( \hbar \) the classical devil’s staircase remains robust with respect to quantum fluctuations while at large values synchronization plateaus are destroyed (see Fig.11U). Quantum synchronization in the model has close similarities with Shapiro steps in Josephson junctions [UR3b].
Figure 12: Fig. 12U. Standard map as a detector: Husimi function in action-angle variables \( -\pi \leq p <\pi, 0 \leq x < 2\pi \) for qubit state up (left) and qubit state down (right) at \( K=4.5, \epsilon_c = 0.8\), \( \delta=0.1 \), \( \hbar =0.0123 \) at time \( t = 20 \). The initial states of kicked rotator and qubit (spin) are a Guassian packet centered at the fixed point \( p=0, x=\pi \) and spin state \( (|0>+|1>)/\sqrt{2} \); color show density with blue for zero and red for maximum. (from [UR4])
  • UR4) The properties of kicked rotator as a deterministic detector of qubit (spin or two state system) are analyzed in [UR4]. Th Hamiltonian of the whole system is a sum of kicker rotator Hamiltonian \( H_{kr} \), qubit \( H_s = \delta \sigma_x \) and coupling term \( H_{int}= \epsilon_c \sigma_z \cos \theta \sum_m \delta(t-m) \). It is shown that in the regime of quantum chaos the detector acts as a chaotic bath inducing qubit decoherence. The dependence of dephasing and relaxation rates on parameters is established. For a strong qubit-detector coupling the dephasing rate is given by the Lyapunov exponent of classical dynamics. For the strong coupling the detector performs an efficient measurement of qubit (see Fig.12U). In the case of weak coupling, due to chaos, the dynamical evolution of the detector is strongly sensitive to the state of the qubit. However, in this case it is unclear how to extract a signal from any measurement with a coarse-graining in the phase space on a size much larger than the Planck cell.
  • UR5) The model of two kicked rotators with short range and finite range interactions in the momentum space is analyzed in [UR5a], [UR5b]. It is shown that the interaction leads to a strong enhancement of localization length. The case of interaction between two particles in higher effective dimensions is considered in [42]a) for the case of frequency modulated kicked rotators with two (\( d_{eff}=3 \)) and three (\( d_{eff}=4 \)) modulation frequencies. It is shown that in such models the interactions create delocalized pairs of particle in the regime when all one-particle states are exponentially localized. While without interactions the above models have been realized with cold atoms in kicked optical lattices (see [32] for the kicked rotator and [UR5c] for the case of two frequencies, described in more detail in UE4) below, the realization of local interactions in the momentum space is rather difficult for such systems.
  • UR6) The frequency modulated kicked rotator model was introduced in [46] for the quantum case. The corresponding classical volume preserving map has the form \( \bar{y}=y - (k + \epsilon \cos z) \sin x , \) \( \bar{x} = x + \bar{y} , \) \( \bar{z} = x +2\pi r_2 \). The destruction of the spiral mean torus with the rotation numbers \( r_1= 1/\lambda^2\), \( r_2 = 1/\lambda\), with \( \lambda^3 - \lambda -1 =0\), \( \lambda =1.324718...\) is analyzed in [UR6a]. The critical torus exists along a critical curve in the plane \( (k, \epsilon) \). In a certain interval of this curve the Greene residue dynamics with the renormalization time step is not universal indicating an emergence of dynamical renormalization chaos. Such a behaviour is strikingly different from the case of the critical golden curve in the standard map where the renormalization dynamics is universal corresponding to a fix point. Further analysis of critical tori in this map is reported in [UR6b]. An approximate renormalization-group transformation for Hamiltonian systems with three degrees of freedom is constructed in [UR6c].
  • UR7) The effects of nonlinearity on localization in kicked rotator (2) are analyzed in [UR7a] by adding after each kick a nonlinear phase shift of wavefunction amplitudes \( \psi_n \) in the momentum representation \( \bar{\psi_n} = \exp(i\beta |\psi_n|^2) \psi_n \). The model was called the kicked nonlinear rotator (KNR). It is argued that there a certain critical strength of nonlinearity \( \beta_c \) below which the localization is essentially preserved, while for \( \beta > \beta_c\) a subdiffusive spreading over momentum states \( n\) takes place in time with \( n^2 \propto t^\alpha \) with the exponent \( \alpha=2/5\). The further numerical studies [UR7b] give a smaller value \( \alpha = 0.35 \pm 0.03 \) where the error bar is obtained from averaging over 10 realizations, taken from different rotation phases \( n(n+\zeta)/2\) with \( t \leq 10^8 \); one realization with \( t \leq 10^9 \) has \( \alpha = 0.35 \) (the case of random rotational phases instead of \( n(n+\xi)/2\) has the same exponent). In [UR7a] it was argued that KNR describes the nonlinear spreading over momentum harmonics in the kicked Gross-Pitaevskii equation on a ring \( (0 \leq x < 2\pi ) \) studied in [49]

\[\tag{5} i \hbar{{\partial {\psi}} \over {\partial t}} = -{{1} \over {2}}{{\partial^2 {\psi}} \over {\partial^2 {x}}} +{\beta} |\psi|^2 \psi - k \cos x \; \psi \sum_{m=-\infty}^{\infty} \delta(t-mT) \; . \]

This is confirmed in further numerical simulations [UR7c] with the numerically found exponent \( \alpha \approx 0.4 \) for \( t \leq 10^7 \). In [UR7c] the analogy between Kolmogorov energy flow from large to small spacial scales and conductivity in disordered solid state systems is proposed for model (5). It is argued that the Anderson localization can stop such an energy flow. The effects of nonlinear wave interactions on such a localization are analyzed. The results obtained for finite size systems show the existence of an effective chaos border between the Kolmogorov-Arnold-Moser (KAM) integrability at weak nonlinearity, when energy does not flow to small scales, and developed chaos regime emerging above this border with the Kolmogorov turbulent energy flow from large to small scales. Another conjecture, pushed forward in [UR7a], is that the destruction of Anderson localization in 1D disordered potential is characterized by the same subdiffusive spreading as for the KNR. This conjecture was confirmed in extensive numerical simulations with the discrete Anderson nonlinear Schrodinger equation (DANSE) (see [UR7c] and Refs. therein).

Figure 13: Fig. 13U. Dependence of chaotic layer width \( w_s \) on the adiabaticity parameter \( \lambda=1/\sqrt{K} \) is shown by broken solid lines connecting various symbols with resonance dimension \( L=N \) indicated by numbers. The rate of Arnold diffusion inside chaotic layer is \( D \sim {w_s}^2 \). Data demonstrate nonexponential (approximately algebraic) decay of \( w_s \) on \( \lambda \) (from [UR8b] with more details given there).

UR8) The model of coupled standard maps with nearest left-right neighbour couplings \( K \) is used for investigations of classical chaos properties in many-body (or many dimensional) systems [UR8a],[UR8b]. In [UR8b] a skillful numerical method is used to compute the width of chaotic layers and Arnold diffusion rate inside the layer. The method is based on the computation of unstable fix points in high-dimensional phase space, and then determination of rotational period inside the layer via a certain number of trajectories [UR8b]. This approach allows to determine very small layer width \( w_s \sim 10^{-22} \) and the related fast Arnold diffusion inside the layer \( D \sim {w_s}^2 \sim 10^{-44} \) (see Fig.13U from [UR8b]). The main result obtained in [UR8b] is approximately algebraic (nonexponential) decay of the Arnold diffusion with a decrease of perturbation parameter \(K\) or increase of the related adiabaticity parameter \( \lambda =1/\sqrt{K} \). Certain explanations are proposed in [UR8b] but the origin of this slow decay of \( D(\lambda) \) remains unclear. Further studies of this model [UR8c] show that at small coupling the measure of chaos is found to be proportional to the coupling strength with the typical maximal Lyapunov exponent being proportional to the square root of coupling. This strong chaos appears as a result of triplet resonances between nearby modes. The dynamics in such triplets remains chaotic even for \( K \rightarrow 0 \). In addition to strong chaos there is a weakly chaotic component having much smaller Lyapunov exponent, the measure of which drops approximately as a square of the coupling strength (\( K^2 \)) down to smallest couplings reached. It is argued that this weak chaos is linked to the regime of fast Arnold diffusion discussed in [UR8b]. The investigation of Poincare recurrences in this model shows that their statistics is characterized by the algebraic decay \( P \propto 1/t^\alpha \) with the Poincare exponent \( \alpha \approx 1.3 \) being independent of number of degrees of freedom [UR8d]. A conjecture is made about universal value of the Poincare exponent in systems with many degrees of freedom [UR8d]. A certain confirmation of this conjecture is given by the numerical results of Poincare recurrences in protein and DNA molecules where a similar value of the Poincare exponent is obtained from numerical simulations [UR8e],[UR8f]. Finally, the case of the many standard maps with a relatively strong kick amplitude \( K\) and weak couplings between all maps is considered in [UR8g]; the case of 4D coupled standard maps is studied in [UR8h]a) with approximate \( \alpha \approx 1.6 \) found; a similar model is studied in [UR8h]b) analyzing the influence of recurrence set choice.

Experimental realizations

Figure 14: Fig. 14U. Root mean square of occupied levels \( \Delta n\) in modulated quantum pendulum versus the normalized modulation amplitude \( \lambda\) at \( k=15\), \(\hbar_{eff}=1.58\); classical results are shown by stars connected by dashed lines, quantum results are shown by full curve and the quantum localization theory is shown by dashed curve. (from [UE1a], this Fig. is also used in [UE1b]).

Various experiments with systems related to the standard map are described here including cold atom experiments with phase modulated effective pendulum (UE1), high-order quantum resonances observed with BEC and noncondensed cold atoms (UE2), time reversal of BEC atomic waves in quantum chaos regime (UE3), cold atom experiments with frequency modulated kicked rotator and observation of Anderson transition (UE4), realizations of kicked rotator with molecules in pulsed laser field (UE5), stabilization theory of electron edge states in magnetic and microwave fields and its experimental observation (UE6), observation of the Aubry transition for cold ions in optical periodic lattice and thermoelectric properties of this system (UE7).

  • UE1) The classical dynamics of phase modulated pendulum, described by the rescaled dimensionless Hamiltonian \( H = p^2/2 - k \cos(\phi - \lambda \sin t) \) is analyzed in [1],[10]. In the regime of fast resonance crossing \( k/\lambda > 1, \lambda \gg 1\) the dynamics is approximately described by the standard map since each crossing gives a kick to momentum of particle. The chaos region is restricted to \( |p| < \lambda \) since the crossing is possible only in this region. From the map the chaos border is \( k > 0.04 \sqrt{\lambda} \) with a diffusion rate \( D \approx k^2/\lambda \). The quantum case [UE1a] is characterized by an effective dimensionless Planck constant \( \hbar_{eff} \) with the dynamical localization length \( \ell = \pi D/{\hbar_{eff}}^2 \). Thus in the classical case the fluctuations of momentum (or average populated number of quantum states \( \Delta n \)) grow linearly with \( \lambda \) while in the quantum case at they drops at large \( \lambda \) as \( \Delta n \approx \ell \propto D \propto 1/\lambda\) (see Fig.14U). It is argued that this model describes the quantum dynamics of Josephson junctions at small dissipation. In [UE1b] it is shown that the same Hamiltonian describes cold atoms in a modulated standing light wave. This physical system happens to be more accessible for experimental realization and the dynamical localization of quantum chaos in an effective modulated pendulum is observed by Raizen group in [UE1c].
  • UE2) After the realization of kicked rotator with cold atoms by Raizen group [32] the properties of this model have been studied by different experimental groups. Thus the high-order quantum resonances predicted in [28]a) are observed with BEC [UE2a] and noncondensed cold atoms [UE2b]. The cold atom experiments with a kicked rotator (particle) under an applied static field are discussed at the article Kicked cold atoms in gravity field.
Figure 15: Fig. 15U. Time reversal of Bose-Einstein condensate of Rb cold atoms in the regime of quantum chaos for \( \epsilon =1, k = K/T \approx 2\) (see text); top panel shows probability distribution in momentum (expressed in recoil units) after each kick; the reversal time is \( t_r=5\); the bottom left panel shows the distribution at the return moment \( 2t_r=10\) (experimental data are shown by the solid blue curve, numerical simulations by the dashed green curve) and the right panel shows zoom of the left panel with initial distribution shown by dashed red curve. (from [UE3d]).
  • UE3) The famous dispute between Boltzmann [UE3a] and Loschmidt [UE3b] on time reversal of moving atoms (see also [45]) remained without any experimental verification from 1876 till recently since it is rather hard to invert time for matter waves. The method for realization of time reversal of atomic waves and BEC has been proposed for cold atoms [48] and BEC [UE3c] moving in kicked optical lattice in the regime of quantum chaos for kicked rotator. This is reached by propagating the evolution described by (2) or (5) during time \( t_r \) at \( T =4\pi + \epsilon\) and then replacing the kick period \( T \rightarrow 4\pi - \epsilon\) and displacing the lattice in \( x\) by \( \pi\) (see (5)) that generates an approximate time reversal of atoms with small initial momentum at time \( 2t_r\). This theoretical proposal is realized experimentally with BEC of Rb atoms [UE3d] with the time reversal after \( t_r=5 \) kicks and return time \( 2t_r=10\), as it is shown in Fig.15U. Even if the classical dynamics of this model is deeply in the chaotic regime with \( K \approx 22 \) the quantum system returns close to the initial distribution. Of course, if would be desirable to increase \( t_r\) by a factor 10 to 20 with a larger kick amplitude \( k \sim 10 \) so that after such a time the classical trajectories simulated on a computer would not return to the origin (see Fig.7) in contrast to the quantum evolution (see Fig.6). It is possible to hope that this first experimental test for the Bolzmann-Loschmidt dispute will be extended to the above parameters to understand in a better manner the properties of time reversal for dynamics of classical and quantum chaotic atomic motion.
  • UE4) The frequency modulated kicked rotator is introduced in [46]. It is described by the Hamiltonian \(H = T n^2/2 + k (1+ \epsilon {\Pi_{m=1}}^{n_f} \cos (\omega_m t)) \cos \theta \delta_p(t-1) \), where \( n =-i \partial/\partial \theta \) and \( T =\hbar\) plays the role of effective Plank constant; \( \delta_p(t-1)\) is periodic delta-function of unit period and \(\Pi\) notes the product. Since the phases \( \theta_m = \omega_m t\) evolve linearly with time it is possible to go to extended phase space with additional actions \( n_m = -i \partial/\partial \theta_m \) so that the system will have the effective dimension \( d_{eff}=n_f+1\). The frequencies \( \omega_m \) are incommensurate between themselves and \( 2\pi\). Thus the case with \( d_{eff}=2\) is studied in [46] showing that the number of excited levels is growing exponentially with \( \epsilon, k = K/T=K/\hbar \) (see Fig.9 in [46]) corresponding to Anderson localization in two dimensions. The case of \( d_{eff}=3,4\) is studied in [42]a) , [UE4a] (the case of random phases is considered there instead of \( Tn^2/2 \)) showing that there is the Anderson transition from exponential localization of probability distribution over levels to a diffusive spreading over them above a certain delocalization border for \( k \) at given \( \epsilon >0 \). Thus the critical values are \( k_c \approx 1.8\) at \( \epsilon =0.75, d_{eff}=3\) and \( k_c \approx 1.15\) at \( \epsilon =0.9, d_{eff}=4\) [42]a). The critical exponents for the localization length and diffusion rate in the vicinity of critical point are found to be in agreement with the renormalization theory for Anderson transition [UE4a]. The results obtained in [42]a) , [46] attracted interest of cold atom experimental groups. The model with \( d_{eff}=3\) is studied experimentally by Garreau group [UR5c] with cold cesium atoms finding the transition at \( k = K/T \approx 1.88\) at \(\epsilon =0.75, T=\hbar = 2.89\) being in agreement with numerical simulations and the above value found in [42]a) (the difference in the rotational phases is not very important since for quadratic rotational phases the classical chaos border \( K \approx 0.3\) established in [UR6a] is much below the Anderson critical point). The Garreau group experiments [UR5c], [UE4b] succeeded to obtain experimentally the critical exponents in a vicinity of 3D Anderson transition that is strikingly exceptional and was never realized in the solid state experiments. The overview of the Garreau group experiments is given in [UE4b]. The experiments for 2D case demonstrated exponential growth of the localization length confirming the early results obtained in [46] (Fig.9 there). Since the numerical simulations of the frequency modulated kicked rotator are very efficient (evolution takes place only in one dimension) this model and its extensions are investigated with various physical effects including quantum Hall effect in two dimensions [UE4c], metallic phase of the quantum Hall in four dimensions with computations of critical exponents [UE4d], topological quantum phenomena with spin-half quasiperiodic quantum kicked rotators [UE4e]. However, an experimental realization of these models with cold atoms is challenging. The frequency modulated kicked rotator with up to 10 frequencies is studied numerically in [UE4f], the critical exponent for the diffusion rate in the critical point vicinity is found to be in agreement with the results of renormalization theory, however, certain deviations are found for the exponent in the localization phase at large dimensions, even if this can be related to a restricted computation times since large localization length requires times \( t_{comp} \gg l^{d_{eff}} \).
  • UE5) The kicked rotator model is realized with nitrogen molecules kicked by a periodic train of femtosecond laser pulses [UE5a], [UE5b], [UE5c]. These experiments allows to realize about 15 kicks with up to 25 rotational states. These experiments demonstrated the effects of quantum resonance and dynamical localization of quantum chaos.
  • UE6) For edge channels of electron transport in two-dimensional electron gas (2DEG) under perpendicular magnetic field and microwave field it is shown that the semiclassical electron dynamics is described by the standard map [UE6a], the dynamics of orbits touching the edge is stabilized by the map principal resonance leading to a vanishing longitudinal resistance of edge states. The signatures of this stabilization are observed with 2DEG high mobility samples [UE6b].
Figure 16: Fig. 16U. Dependence of ZT figure of merit of Wigner crystal in a periodic potential on parameters rescaled temperature \( T/K_c\) and periodic potential amplitude \( K/K_c\) with \( K_c\) being the critical amplitude at Aubry transition; ZT is proportional to color (red for maximum, blue for zero; contour curves show ZT values). (from [UE7d]).
  • UE7) It is shown in [53] that a Wigner crystal of cold ions placed in a periodic optical lattice potential has the Aubry transition when the potential amplitude \( K \) becomes larger than a certain critical value \( K_c \) (measured in units of Coulomb interaction strength on a unit length at period \( 2\pi \)). This system is locally described by the Frenkel-Kontorova model and corresponding standard map with the chaos border parameter \( K_{eff} = 0.5 K_c (2\pi/\nu)^3 \approx 1\),where \( \nu \) is the charge density per period. Thus for a given density the Aubry transition takes place at \( K_c \approx 0.034 (\nu/\nu_g)^{3} \) where for the golden mean density \( \nu_g = 1.618...\) the numerics gives a more exact, but close, value \( K_c =0.0462 \) [53]. Below \( K_c \) the ion chain can easily slide while above the transition it is pinned by the potential. This strongly affects the friction of ion chain and in this way the transition is observed with cold Yb ions by the Vuletic group [UE7a]. Even if the number of ions remains small, up to 5, due to experimental restrictions in [UE7a], other groups start to observe signatures of Aubry transition with tens of ions [UE7b]. The properties of this system are rather nontrivial and their theoretical and experimental investigations are important for understanding of physics of friction on nanoscale as discussed in [UE7c]. It is striking that in the Aubry phase the system has remarkable thermoelectric properties with large values of Seebeck coefficient and such high figure of merit values as \( ZT \approx 4 \) [UE7d] (see Fig.16U). Thus it is rather plausible that this model will allow to understand the main physical features of thermoelectricity which foundations have been done by Ioffe in far 1957 [UE7e]. Due to important technological applications of thermoelectric materials with high figure of merit various materials are actively investigated with first-principles calculations (see e.g, [UE7f]). However, according to [53] the quantum properties of this system are rather nontrivial, being similar to a dynamical version of spin glass systems with a huge quasi-degeneracy in the ground state vicinity, and hence it remains questionable if these first-principle calculations are able to describe correctly thermoelectricity in the quantum Aubry phase.

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).


G.G.Carlo, L.Ermann, A.M.F.Rivas, M.E.Spina, "Classical properties of the leading eigenstates of quantum dissipative systems", arXiv:1705.03847(quant-ph) (2017).


G.G.Carlo, L.Ermann, A.M.F.Rivas, M.E.Spina, D.Poletti "Classical counterparts of quantum attractors in generic dissipative systems", Phys. Rev. E 95: 062202 (2017).


D.J.D. Earn, S.Tremaine, "Exact numerical studies of Hamiltonian maps: Iterating without roundoff error", Physica D 56: 1 (1992).


Z.Levnajic, I.Mezic, "Ergodic theory and visualization. II. Fourier mesochronic plots visualize (quasi)periodic sets", Chaos 25: 053105 (2015).


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).


A.I.Shnirelman, "On asymptotic multiplicity of spectrum of Laplace operator", Usp. Mat. Nauk 30: 265 (1975) (in Russian); A.I.Shnirelman, addendum in V.F.Lazutkin, "KAM Theory and semiclassical approximations of eigenfunctions, Springer, Berlin (1993).


B.V.Chirikov, D.L.Shepelyansky, "Shnirelman peak in level spacing statistics", Phys. Rev. Lett. 74: 518 (1995), [26]


K.M.Frahm, D.L.Shepelyansky, "Quantum localization in rough billiards", Phys. Rev. Lett. 78: 1440 (1997).


O.V.Zhirov, D.L.Shepelyansky, "Quantum synchronization", Eur. Phys. J. D 38: 375 (2006).


S.Shapiro, "Josephson currents in superconducting tunneling: the effect of microwaves and other observations", Phys. Rev. Lett. 11: 80 (1963).


J.W.Lee, A.V.Averin, G.Benenti, D.L.Shepelyansky, "Model of a deterministic detector and dynamical decoherence", Phys. Rev. A 72: 012310 (2005).


D.L.Shepelyansky, "Coherent propagation of two interacting particles in a random potential", Phys. Rev. Lett. 73: 2607 (1994).


F.Borgonovi, D.L.Shepelyansky, "Enhancement of localization length for two interacting kicked rotators", Nonlinearity 8: 877 (1995).


J.Chabe, G.Lemarie, B.Gremaud, D.Delande, P.Szirftgiser, J.C.Garreau, "Experimental observation of the Anderson metal-insulator transition with atomic matter waves", Phys. Rev. Lett. 101: 255702 (2008).


R.Artuso, G.Casati, D.L.Shepelyansky, "Break-up of the spiral mean torus in a volume-preserving map", Chaos, Solitons & Fractals 2(2): 181 (1992).


A.M.Fox, J.D.Meiss, "Greene’s residue criterion for the breakup of invariant tori of volume-preserving maps", Physica D 243(1): 45 (2013).


C.Chandre, H.R.Jauslin, G.Benfatto, A.Celletti, "Approximate renormalization-group transformation for Hamiltonian systems with three degrees of freedom", Phys. Rev. E 60: 5412 (1999).


D.L.Shepelyansky, "Delocalization of quantum chaos by weak nonlinearity", Phys. Rev. Lett. 70: 1787 (1993).


L.Ermann, D.L.Shepelyansky, "Destruction of Anderson localization by nonlinearity in kicked rotator at different effective dimensions", J. Phys. A: Math. Theor. 47: 335101 (2014).


D.L.Shepelyansky, "Kolmogorov turbulence, Anderson localization and KAM integrability", Eur. Phys. J. B 85: 199 (2012).


K.Kaneko, T.Konishi, "Diffusion in Hamiltonian dynamical systems with many degrees of freedom", Phys. Rev. A 40(10): 6130 (1989).


B.V.Chirikov, V.V.Vecheslavov, "Arnold diffusion in large systems", JETP 85(3): 616 (1997) (Zh. Eksp. Teor. Fiz. 112: 1132 (1997), [27]


M.Mulansky, K.Ahnert, A.Pikovsky, D.L.Shepelyansky, "Strong and weak chaos in weakly nonintegrable many-body Hamiltonian systems", J. Stat. Phys. 145: 1256 (2011).


D.L.Shepelyansky, "Poincare recurrences in Hamiltonian systems with a few degrees of freedom", Phys. Rev. E 82: 055202(R) (2010).


D.E.Shaw, P.Maragakis, K.Lindorff-Larsen, S.Piana, R.O.Dror, M.P.Eastwood, J.A.Bank, J.M. Jumper, J.K.Salmon, Y.Shan, W.Wriggers, "Atomic-level characterization of the structural dynamics of proteins", Science 330: 341 (2010).


A.K.Mazur, D.L.Shepelyansky, "Algebraic statistics of Poincare recurrences in a DNA molecule", Phys. Rev. Lett. 115: 188104 (2015).


E.G.Altmann, H.Kantz, "Hypothesis of strong chaos and anomalous diffusion in coupled symplectic maps", EPL 78: 10008 (2007).


a) S.Lange, A.Backer, R.Ketzmerick, "What is the mechanism of power-law distributed Poincaré recurrences in higher-dimensional systems?", EPL 116: 30002 (2016); b) M.Sala, R.Artuso, C.Manchein, "Anomalous dynamics and the choice of Poincare recurrence set", Phys. Rev. E 94: 052222 (2016).


R.Graham, M.Schlautmann, D.L.Shepelyansky, "Dynamical localization in Josephson junctions", Phys. Rev. Lett. 67: 255 (1991).


R.Graham, M.Schlautmann, P.Zoller, "Dynamical localization of atomic-beam deflection by a mudulated standing light wave", Phys. Rev. A 45: R19 (1992).


F.L.Moore, J.C.Robinson, C.Bharucha, P.E.Williams, M.G.Raizen, "Observation of dynamical localization in atomic momentum transfer: a new testing ground for quantum chaos", Phys. Rev. Lett. 73: 2974 (1994).


C.Ryu, M.F.Anderson, 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).


J.Kanem, S.Maneshi, M.Partlow, M.Spanner, A.M.Steinberg, "Observation of high-order quantum resonances in the kicked rotor", Phys. Rev. Lett. 98: 083004 (2007).


L. Boltzmann. "Über die Beziehung eines allgemeine mechanischen Satzes zum zweiten Haupsatze der Wärmetheorie". Sitzungsberichte der Akademie der Wissenschaften, Wien, Vol. II 75, pp. 67-73 (1877).


J. Loschmidt. "Über den Zustand des Wärmegleichgewichts eines Systems von Körpern mit Rücksicht auf die Schwerkraft", Sitzungsberichte der Akademie der Wissenschaften, Wien, Vol. II 73, pp. 128-142 (1876).


J.Martin, B.Georgeot, D.L.Shepelyansky, "Time reversal of Bose-Einstein condensates", Phys. Rev. Lett. 101: 074102 (2008).


A.Ullah, M.D.Hoogerland, "Experimental observation of Loschmidt time reversal of a quantum chaotic system", Phys. Rev. E 83: 046218 (2012).


D.L.Shepelyansky, "Anderson transition in three and four effective dimensions for the frequency modulated kicked rotator", arXiv:1102.4450(cond-mat.dis-nn) (2011).


J.-C.Garreau, "Quantum simulation of disordered systems with cold atoms", Comptes Rendus Physique 18: 31 (2017).


J.P.Dahlhaus, J.M.Edge, J.Tworzydlo, C.W.J.Beenakker, "Quantum Hall effect in a one-dimensional dynamical system", Phys. Rev. B 84: 115133 (2011).


J.M.Edge, J.Tworzydlo, C.W.J.Beenakker, "Metallic phase of the quantum Hall effect in four-dimensional space", Phys. Rev. Lett. 109: 135701 (2012).


C.Tian, Y.Chen, J.Wang, "Emergence of integer quantum Hall effect from chaos", Phys. Rev. B 93: 075403 (2016).


H.S.Yamada, F.Matsui, K.S.Ikeda, "Critical phenomena of dynamical delocalization in a quantum Anderson map", Phys. Rev. E 92: 062908 (2015).


J.Flob, A.Kamalov, I.S.Averbukh, P.H.Bucksbaum, "Observation of Bloch oscillations in molecular rotation", Phys. Rev. Lett. 115: 203002 (2015).


M.Bitter, V.Milner, "Experimental observation of dynamical localization in laser-kicked molecular rotors", Phys. Rev. Lett. 117: 144104 (2016).


M.Bitter, V.Milner, "Experimental demonstration of coherent control in quantum chaotic systems", Phys. Rev. Lett. 118: 034101 (2017).


A.D.Chepelianskii, D.L.Shepelyansky, "Microwave stabilization of edge transport and zero-resistance states", Phys. Rev. B 80: 241308(R) (2009).


A.D.Levin, Z.S.Momtaz, G.M.Gusev, A.K.Bakarov, "Microwave-induced nonlocal transport in a two-dimensional electron system", Phys. Rev. B 89: 161304(R) (2014).


A.Bylonskii, D.Gangloff, I.Countis, V.Vuletic, "Observation of Aubry-type transition in finite atom chains via friction", Nat. Mater. 11: 915 (2015).


J.Kiethe, R.Nigmatullin, D.Kalincev, T.Schmirander, T.E.Mehlstaubler, "Probing nanofriction and Aubry-type signatures in a finite self-organized system", Nat. Comm. 8: 15364 (2017).


N.Manini, G.Mistura, G.Paolicelli, E.Tosatti, A.Vanossi, "Current trends in the physics of nanoscale friction" Adv. Phys. X 2(3): 569 (2017).


O.V.Zhirov, D.L.Shepelyansky, "Thermoelectricity of Wigner crystal in a periodic potential", EPL 103: 68008 (2013).


A.F.Ioffe, "Semiconductor thermoelements and thermoelectric cooling", Infosearch, London (1957).


S.Bang, D.Wee, A.Li, A.Fornari, B.Kozinsky, "Thermoelectric properties of pnictogen-substituted skutterudites with alkaline-earth fillers using first-principles calculations", Jour. Appl. Phys. 119: 205102 (2016).

See also added after 2008

Anderson localization and quantum chaos maps, Chirikov criterion, Cold atom experiments in quantum chaos, Google matrix, Kicked cold atoms in gravity field, kicked Harper model, Microwave ionization of hydrogen atoms

Personal tools

Focal areas