Statistical properties of a modified standard map in quantum and classical regimes

We present a model—a modified standard map. This model has interesting properties that allow quantum–classical correspondences to be studied. For some range of parameters in the classical phase space of this model, there exist large accelerator modes. We can create a family of maps that have large accelerator modes.

δ-function: n δ(t − nT ), then it is possible to integrate it over one period in time and, in the end, obtain a discrete map [1][2][3]. One of the most famous models of this type was introduced by Chirikov [4], and it is known as a standard or Chirikov map (SM). In SM, perturbation is realized by periodic kicks in time and by a sine function of the position. This model is crucial from a theoretical point of view, as a quite simple model that exhibits Hamiltonian chaos. It has found many applications in accelerator physics, plasma physics, and solid-state physics. In this article, we investigate properties of the modified standard map (MSM) introduced in [5] and its quantum counterpart. In this model, a two-dimensional rotator is also perturbed in time by force modeled by a periodic Dirac δ-function, but we changed the functional dependency of this perturbation on position. The Hamiltonian function of our model has the following form: where P is the momentum, I is moment of inertia, t is time, T is the period of the force, and V (Q) is a periodic function-external potential. We can write the Hamilton equations for the dynamics of our model as where f (Q) and V (Q) are related by the standard formula: In the case of SM, perturbation depends on position Q according to the formula V (Q) = cos Q.
For the models considered in this article, the socalled accelerating modes (AMs) emerge [4,6,7] in phase space for a wide range of parameters. There are areas of the phase space where momentum increases (or decreases) approximately linearly in time. In the case of SM, there are two kinds of AM. However, areas of AM in SM are quite small. Using our modifications, we can create AMs that are several times larger and parametrically control their areas. We can also propose modification when only one kind of AMs exist with quite large areas. Simultaneously, trajectories of the area outside the accelerator modes can "stick" to their boundaries and follow a linear evolution for a very long time. This phenomenon leads to statistically interesting results in the form of power-like tails in the momentum distribution function [8].
In this paper, we study the classical and quantum dynamics of a system, for a particular range of parameters for which AMs exist. Because AMs have an important practical interpretation [9][10][11][12], we hope that our results may be relevant in this area.

Modified standard map
Equations (2) can be integrated over time, and this leads to classical discrete dynamical system, which is a modification of SM: where Q n is a position at the nth step, P n represents momentum in the nth step, and f (Q) is a piecewise linear continuous periodic function: Let us recall that the iteration steps of the map correspond to equal time increments. It is also assumed that the mean value of the momentum increment P n+1 − P n is equal to zero. We additionally require that there exists a stable AM with a center at Q = Θ 1 /2, and this implies the following equation: These assumptions allow one to express parameters a i , b i , i = 1, 2, 3 as a linear function of one parameter, say a 1 . Moreover, one can easily check that the solution (Q n , P n ) = (Θ 1 /2, 2nπ) is an elliptic point for a 1 ∈ (− 4, 0), and therefore, only for this range of values, there exists a stable AM. With this modification, it is much easier to investigate the quantum counterpart of the classical map, especially effects like tunneling from the AM or the stickiness of its boundary (see Fig. 1). Using a proper choice of parameters, the area of AM in the phase space is relatively large (see Fig. 2). It should be noted that MSM as well as SMs is obtained from Hamilton equations in accordance with the assumption that the time dependence of the acting force is a periodic Dirac δfunction; thus, we are considering the so-called kicked model. This also allows the exact form of the quantum map to be obtained. However, in the case of continuous and smooth time dependence in the quantum description, one can apply the well-known Floquet Hamiltonian.

Quantum map
In this section, we construct a quantum counterpart of Equations (4). With this aim, we start from Schrödinger's equation with the Hamilton operator in the following form: where the second term represents periodic kicks with period T . Due to the form of the Hamiltonian, the evolution operator from time t = 0 to t = T may be written as a product of two simple operators: Points belonging to AM escape to the large value of momentum whereĤ 0 is the operator of free rotation and we assume that we start from free rotation. It should be noted that each operator in the above equation is a simple multiplication in the corresponding representation:Û Q is a multiplication in position representation, andÛ P is a multiplication in momentum representation. In our case, the phase space is a cylinder with periodicity in variable Q; therefore, we have discrete momentum representation with P m = mh, m = . . . , −2, −1, 0, 1, 2, . . .. In numerical computations, we take into account only a finite number of Pvalues say P m = −Nh, . . . , (N − 1)h and conse-quently a finite number of discrete Q-values: The transitions between both representations are given by the discrete Fourier transform and its inversion: In the classical map (4), the area , the classical mapping is invariant with respect to the translation of 2kπ in momentum. The number of quantum states that belong to such an elementary cell is strictly determined by the value of constanth. However, we can treath as a parameter that may take any positive value. Then, the number of quantum states belonging to the elementary cell can be calculated as 2 [π/h] + 1, where symbol [·] means an integer part of the fraction π/h or as 2π/h, when multiples ofh are equal to π .

Classical and quantum correspondence
In the studies of relationships between classical and quantum dynamics, it is useful to represent a quantum state in phase space. Such possibility gives us a Wigner function and Husimi-Kano function. In our case, the phase space has a cylindrical geometry and position is a cyclic variable. For that geometry of phase space, the Wigner function can be written in the following way [13]: (11) and the Husimi function on cylinder is defined by the following formula [14,15]: We analyze the quantum evolution of a state given by: where D ∈ (real number), Q 0 and P 0 denote the mean value of position and momentum, respectively. When we choose the numbers Q 0 and P 0 , we can localize the maximum of the Husimi function in various regions of the classical phase space both in accelerator modes and in the chaotic part of phase space. We can also transform initial Gaussian states to obtain a rotated Husimi or Wigner function. Such a situation is presented in Fig. 3 where both functions are drawn as contours inside the AM which has the shape of a rotated ellipse.
This initial Gaussian state can be obtained by rotating the Wigner function (13) by the angle ϕ and finding its wave function ψ ϕ (Q). In the case of the plane as a phase space, one obtains the following form of this function: where We can see that the transformation from ψ(Q) to ψ ϕ (Q) is equivalent to replacing D with a complex parameter. We can say that the function ψ ϕ (Q) is a product of Gaussian and linear chirp. The formula above can also be used in the cylinder case if the Gaussian function is sufficiently narrow, i.e., if its width is much smaller than 2π .
In Fig. 4, we present the Husimi function for the initial state (rotated Gaussian) localized in AM in the cell described by the Cartesian product: [0, 2π) × [−π, π) and five of its iterations for n = 1, 2, 10, 100, 1000. We can observe that the Husimi functions of these iterations are localized in AMs in cells of phase space described by [0, 2π ] × [(2n − 1)π, (2n + 1)π ], where n is the iteration number. Moreover, the form of the Husimi function changes slightly. However, there exists tunneling from AM to the remaining phase space. In order to investigate the properties of this tunneling, we calculate the probability W (n) = i |ψ(P i )| 2 , where P i ∈ [(2n − 1)π, (2n + 1)π ] at each iteration of the initial state. When proceeding in this way, however, we are calculating rather the escape rate from the cell with AM than from AM itself.
The results of numerical calculations for the chosen value of D andh are presented in Figs. 5 and 6. In the first one, we present Wigner functions of the initial state . We can see that the state is still well localized in AM as W (4000) ≈ 0.996. These two pictures give an alternative presentation of states to the Husimi representation.
The left side of Fig. 6 presents a graph of function ln(W (n)) and suggests that it can be described by the formula: with β = 9.46×10 −7 . Such a good fit was obtained for the same initial state as in the previous figure. On the right side of Fig. 6, where the initial state has larger values of parameter D, we obtain a good fit to the straight line for results with n > 5000. A value of β = 9.44 × 10 −7 was estimated for iterations with n > 30000. Such a calculated value of the β parameter is practically the same as in the previous figure (D = 0.358). Similar simulations can be performed when the Husimi function is well localized outside of the AM, in the chaotic part of phase space. In contrast to the situation presented in Fig. 4, the initial state gradually fills the chaotic part of phase space leaving intact regions where AMs are localized, which is clearly visible for the iteration n = 3 in Fig. 7.
We have also studied the asymptotics of numerically obtained densities, where a maximum of Husimi functions lay outside the accelerator mode. In Fig. 8, we present an evolution of a wave function amplitude |Ψ (P, t)| in momentum representation. The amplitude has a maximum, which moves in the direction of the negative P-values with constant velocity. At the same time, it possesses a long tail for positive P-values with the front which moves with constant velocity as well (the front jumps by 2π at every iteration). In order to  and P 0 = 0. For n = 1, 2, 3, we have the next quantum iterative steps. We use following parameters: I = 1.0000018, K = 1, a 1 = −2, 1, θ 1 = π , θ 2 = 3 2 π ,h = π 64 , T = 1 and the momentum range [−8π, 8π ]  Fig. 8 Amplitude |Ψ (n, P)| of a wave function in momentum representation for selected iteration steps n verify the nature of this tail, we consider a state after a large number of iterations as is presented in Fig. 9 on the left side. The value of the wave function amplitude is drawn in two scales, which can be seen on the right side of this figure. The upper figure presents a linear fit to the results in log-log scale that has a high coefficient of determination. It is a strong numerical argument that the asymptotics of the fragment of |Ψ (P)| under consideration has an inverse power character as in the classical case (see Fig. 2 in [5]).
Taking into account the results of our numerical calculations, we can state that basic properties of quantum distribution retain the properties of the classical model (despite quite a large Planck constant). The AMs still have an influence on quantum evolution as in the classical model. Both in the classical picture and the quantum one, densities have a pronounced maximum for some values of momentum while for larger values of momentum they have inverse power asymptotics.

Further modifications of the standard map
Replacing the function f (Q) with a four-part piecewise linear continuous periodic function for which P n+1 − P n = 0. For certain parameter values, we can observe the emergence of two kinds of AM in phase space. In one type of AM, the P value grows approximately linearly during successive iterations. In the other type of AM, the movement goes in the opposite direction. Both types of AM are presented in Fig. 10. We can also modify function (5) requiring that it be differentiable, i.e., that it be a function of class C n . It can be realized by joining its adjacent linear parts by a polynomial of order n + 1 on some small intervals of length D. All coefficients of the linear and polynomial parts can be uniquely determined by periodicity, continuity of functions, and their derivatives. The proper choice of model parameters allows us to obtain AM in the phase space (see Fig. 11). In the limit D → 0, this function becomes f (Q) defined by formula (5).

Conclusions
We have presented a model that is a modification of the SM. The kicking function is a smooth piecewise linear function parametrized by one parameter (besides position of three intervals). For properly chosen parameters, there exist large AMs. This feature simplifies the numerical analysis of the quantum mapping. Our analysis shows that AM and its stickiness property are retained in the quantum picture. The quantum evolution of states localized in AM follows classical evolution, i.e., the quantum map translates the localized state in AM to the next AM, practically without changing its Husimi function, with the exception of very small tunneling to the rest of phase space. Quantum probability density in momentum representation preserves many essential features of classical momentum density. Numerical results show that states which are localized in the chaotic part of phase space delocalize under the influence of quantum mapping in a manner similar to the classical case. For a sufficiently long time (a number of a map iterations), density has an inverse power asymptotic for large momentum with a characteristic front moving with constant velocity. This model is relatively easy to analyze when instead of kick, i.e., δ-time dependence, we use a continuous periodic function. In such a case, one has to apply Floquet theory.

Compliance with ethical standards
Conflict of interest The authors declare that there is no conflict of interests regarding the publication of this paper.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.