arxiv:nlin/ v1 [nlin.cd] 21 Sep 2005

Transcription

1 General Quantum Resonances of the Kicked Particle Itzhack Dana 1 and Dmitry L. Dorofeev 1,2 1 Minerva Center and Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel 2 Department of Physics, Voronezh State University, Voronezh , Russia arxiv:nlin/ v1 [nlin.cd] 21 Sep 2005 Abstract The quantum resonances (QRs) of the kicked particle are studied in a most general framework by also considering arbitrary kicking potentials. It is shown that QR can arise, in principle, for general rational values of the Bloch quasimomentum. This is, in fact, what is found in the case of the main QRs for a potential containing all the harmonics. In this case, which is precisely described by the linear kicked rotor, exact formulas are derived for the diffusion coefficient associated with either an incoherent mixture of plane waves or a general wave packet. The momentum probability distribution is exactly calculated for a two-harmonic potential. It clearly exhibits new resonant values of the quasimomentum and it is robust under small deviations from QR. PACS numbers: Mt, Ac, w, Gg Typeset using REVTEX 1

2 There has been recently a considerable experimental [1] and theoretical [2,3] interest in new remarkable phenomena associated with the quantum resonances of the periodically kicked particle (KP), either in the presence or in the absence of gravity. In this paper, we report about new general results concerning the quantum resonances of the KP in the absence of gravity. These results extend some of previous work [2,3] in a significant way, shedding light on new basic aspects of the problem and connecting an important case of the system with a well-known model. We start with a brief summary of previous results. The quantum KP is described by the Hamiltonian, using notation in Ref. [3], Ĥ = ˆp2 2 + kv (ˆx) t δ(t tτ), (1) where (x, p) are the position and momentum of the particle, k is a nonintegrability parameter, V (x) is a periodic potential, t takes all the integer values, t is the continuous time, and τ is the kicking period. The units are chosen so that the particle mass is 1, the Planck s constant = 1, and the period of V (x) is 2π. In essence, V (x) has been always chosen as the standard potential V (x) = cos(x). The one-period evolution operator for (1), from t = t + 0 to t = t + τ + 0, is given by Û = exp [ ikv (ˆx)]exp ( iτ ˆp 2 /2 ). (2) The translational invariance of (2) in ˆx implies the conservation of a quasimomentum β (0 β < 1): The application of Û on a Bloch function Ψ β(x) = exp(iβx)ψ β (x), where ψ β (x+2π) = ψ β (x), results in a Bloch function Ψ β (x) = exp(iβx)ψ β (x) associated with the same value of β. Here ψ β (x) is the 2π-periodic function ψ β (x) = Ûβψ β (x), where Û β = exp [ ikv (ˆx)] exp [ iτ (ˆp + β) 2 /2 ]. (3) The restriction of the operator (3) to 2π-periodic functions ψ β (x) allows one to interpret x as an angle θ and ˆp as an angular-momentum operator ˆN = id/dθ with integer eigenvalues n. One can then view (3) as the one-period evolution operator for a β-kicked rotor (β-kr). Now, an arbitrary KP wave packet Ψ(x) can be always expressed as a superposition of Bloch functions, Ψ(x) = 1 0 dβ exp(iβx)ψ β(x), where 2

3 ψ β (x) = 1 2π n Ψ(n + β) exp(inx), (4) Ψ(p) being the momentum representation of Ψ(x). One then gets the basic relation Û t Ψ(x) = 1 0 dβ exp(iβx)ût β ψ β(x) (5) for integer time t, connecting the quantum dynamics of the KP with that of β-krs. For typical irrational values of τ/(2π), β-krs are expected to feature dynamical localization in the angular momentum n [4,5], implying a similar localization of the KP wave packet Ψ(p) in p = n + β. If τ/(2π) is rational, the usual (β = 0) KR exhibits quantum resonance (QR), i.e., a ballistic (quadratic in time) growth of its kinetic energy [5]. QR in general β-krs was studied in the case of integer τ/(2π) (main QRs) and V (x) = cos(x) [2,3]. It was found [3] that QR arises in this case only for special resonant values of β, finite in number. Thus, QR is exhibited by the KP only if Ψ(p) is sharply localized on the discrete set of momenta p = n + β, with β in the finite resonant set. However, QR leaves a clear fingerprint in the evolution of either a general KP wave packet (5), involving all values of β, or an incoherent mixture of plane waves: The average kinetic energy in both cases grows diffusively (linearly) in time [2,3]. This diffusive behavior is robust under small deviations ǫ of τ/(2π) from integers, in the sense that it is still observed on time scales t kǫ 1/2 [3]. In this paper, the results above concerning the QRs of the system (1) are fully extended by also considering arbitrary periodic potentials V (x). We show that the general condition for QR in β-krs is just the rationality of both τ/(2π) and β. Thus, a resonant value of β may be a general rational number in [0, 1) and one then has a (countable) infinity of such values. This is, in fact, what we find in the case of the main QRs for a potential containing all the harmonics. In this case, a β-kr is precisely the well-known linear kicked rotor (linear KR) [7,8] and we derive exact formulas for the diffusion coefficient associated with both an incoherent mixture of plane waves and a general KP wave packet. For a mixture uniformly 3

4 distributed in β, the formula essentially coincides with one given by Berry in the context of the linear KR [8]. Multi-harmonic potentials cause the diffusion coefficient for a general wave packet to be explicitly dependent on quantum correlations reflecting the localization features of the initial wave packet in momentum space. The momentum probability distribution is exactly calculated for the potential V (x) = cos(x) + γ cos(2x). It clearly exhibits new resonant values of β due to the second harmonics and it is found numerically to be robust under small deviations of τ/(2π) from integers. The basic origin of QR is a band (continuous) quasienergy (QE) spectrum due to some translational invariance in phase space. For KR systems, this is the invariance of the evolution operator under translations ˆT q = exp( iqˆθ) by q in the angular momentum ˆN [5,6]; here q must be an integer since θ is an angle (0 θ < 2π). The translational invariance of β-krs is expressed by [Ûβ, ˆTq ] = 0, where Ûβ is given by Eq. (3) (ˆx ˆθ, ˆp ˆN). Using the fact that ˆN has integer eigenvalues, one easily finds that [Ûβ, ˆTq ] = 0 holds only if τ 2π = l q, (6) β = r l q 2 mod(1), (7) where l and r are integers. Eq. (6) is the rationality condition for τ/(2π) while Eq. (7) is a formula for the general resonant values of β. Let us write l = gl 0 and q = gq 0, where l 0 and q 0 are coprime integers and g is the greatest common factor of (l, q); the value of τ/(2π) = l 0 /q 0 will be kept fixed in what follows. It is then clear that β in Eq. (7) can take any rational value β r in [0, 1) since g can be always chosen so that r = (β r + gq 0 /2)gl 0 is integer. For given β = β r, we shall choose g as the smallest integer satisfying the latter requirement, so as to yield the minimal values of l = gl 0 and q = gq 0. In general, g > 1, so that (l, q) are not coprime. For the usual KR (β = 0), g = 1 if l 0 q 0 is even and g = 2 if l 0 q 0 is odd (compare with Ref. [6]). For general β = β r, the QE states ϕ can be chosen as simultaneous eigenstates of Ûβ and ˆT q : Û β ϕ = exp( iω)ϕ, ˆTq ϕ = exp( iqα)ϕ, where ω is 4

5 the QE and α is a quasiangle, 0 α < 2π/q. Using standard methods [5,6], it is easy to show from the latter equations that at fixed α one has precisely q QE eigenvalues ω b (α, β), b = 0,..., q 1. As α is varied continuously, these eigenvalues form q QE bands. If not all these bands are flat (have zero width), QR can arise and β = β r is indeed a resonant value. The important case of the main QRs (τ = 2πl 0, q 0 = 1) was studied in Refs. [2,3] for V (θ) = cos(θ), and it was found [3] that QR arises only for l 0 resonant values of β, given by Eq. (7) with r = 0,..., l 0 1 and g = 1 (i.e., l = l 0 and q = 1). In what follows, we study the main QRs for arbitrary V (θ). The case of τ = 2πl 0 is the only one in which the term ( ˆN + β) 2 in Eq. (3) (ˆx ˆθ, ˆp ˆN) can be replaced by the operator ˆN + 2β ˆN + β 2, linear in ˆN; this is because exp( iπl 0 n 2 ) = exp( iπl 0 n) for the integer eigenvalues n of ˆN. Then, after omitting the nonrelevant constant phase factor exp( iπl 0 β 2 ), one can express (3) as [ ] ( ) Û β = exp ikv (ˆθ) exp iτ β ˆN, (8) where τ β = πl 0 (1 + 2β). We identify Ûβ in Eq. (8) as the evolution operator for the wellknown linear KR [7,8] with Hamiltonian Ĥ = τ β ˆN+kV (ˆθ) t= δ(t t); the corresponding Schrödinger equation is exactly solvable for arbitrary potential V (θ), V (θ) = m V m exp( imθ). (9) Assuming the quantum state of the linear KR to be initially (at t = 0) an angular-momentum state, ψ β,0 (θ) = exp(in 0 θ)/ 2π, the expectation value of the kinetic energy in this state at time t is given by the exact expression [8,9]: E n0,β(t) = 1 2 ψ β,t ˆN 2 ψ β,t = n k2 m 2 V m 2sin2 (mτ β t/2) sin 2 (mτ m>0 β /2). (10) If β is a typical irrational number, so that also τ β /(2π) is such, the QE spectrum of the linear KR is pure point [7] and E n0,β(t) in Eq. (10) is a bounded, quasiperiodic function of time [8]. A rational value of β given by Eq. (7) (with l = gl 0, q = g, and g chosen as above) can be written, up to a fixed integer shift in r, as β = β r,g = r/(gl 0 ) 1/2 mod(1), 5

6 where r and g are coprime; the corresponding value of τ β /(2π) is r/g. Then, from work [7], the QE spectrum consists of g bands which are all nonflat if there exists at least one integer j 0 such that the Fourier coefficient V jg in Eq. (9) is nonzero (otherwise, all the bands are flat). Thus, if V (θ) contains all the harmonics (V m 0 for all m), all rational values of β = β r,g are resonant. For τ β /(2π) = r/g and large t, one finds from Eq. (10) a QR (ballistic) behavior, E n0,β(t) St 2 /2, with the stream constant S = k 2 g 2 j j2 V jg 2 [8]. As a simple example, let V (θ) = cos(θ)+γ cos(2θ). The only nonzero coefficients V m are V ±1 = 1/2 and V ±2 = γ/2, so that QR arises only for g = 1, 2. The resonant β values are β r,g = r/(gl 0 ) 1/2 mod(1), where r = 0, 1,..., l 0 1 for g = 1 and r = 1, 3,..., 2l 0 1 for g = 2. The values for g = 1 are just the known ones for V (θ) = cos(θ) [3] while those for g = 2 are new ones, due entirely to the second harmonics. See, however, note [10]. We now consider, as in experimental situations [1], an incoherent mixture of plane waves with initial momentum distribution f(p). The average kinetic energy of this mixture is Ē(t) = 1 0 dβ n 0 f(n 0 + β)e n 0,β (t), where E n 0,β (t) is given by Eq. (10) with n 0 replaced by p = n 0 + β. As in Ref. [3], we define f 0 (β) = n 0 f(n 0 + β) and use the relation 1 0 dβf 0 (β) sin2 [πml 0 (β + 1/2)t] sin 2 [πml 0 (β + 1/2)] t ml 0 ml 0 1 r=0 f 0 (β r,m ) [β r,m = r/(ml 0 ) 1/2 mod(1)], valid in the limit t. We then find for large t a diffusive behavior of Ē(t), Ē(t) D 0 t, where the diffusion coefficient D 0 is given by D 0 = k2 l 0 m>0 ml 0 1 m V m 2 r=0 f 0 (β r,m ). (11) The simple mixture with f(n 0 + β) = δ n0,n 0 is uniformly distributed in β, f 0(β) = 1. In this case, Ē(t) = 1 0 dβe n 0,β(t), leading to an exact equality for all times t: Ē(t) = Ē(0) + D 0t with D 0 = k 2 /(4π) 2π 0 [dv (θ)/dθ]2 dθ. The latter result was essentially obtained by Berry [8] as the kinetic energy of a linear KR whose value of τ is completely unknown. 6

7 Next, we consider a general KP wave packet (5), a coherent mixture of Bloch waves exhibiting all quasimomenta. The kinetic-energy expectation value in Ψ t (x) = Ût Ψ(x) is E t = 1 Ψt ˆp π 2 Ψt dβ dθ dψ β,t (θ) 2 2 dθ, (12) 0 where the last relation, with ψ β,t (θ) = Ût β ψ β(θ), holds for large t provided that E t is unbounded as t [2]; this will turn out to be indeed the case, see below. For τ = 2πl 0, we easily obtain from Rels. (8) and (9) that 0 ψ β,t (θ) = exp [ ik V β,t (θ) ] ψ β (θ tτ β ), (13) V β,t (θ) = t 1 s=0 V (θ sτ β ) = m V m sin(mτ β t/2) sin(mτ β /2) e im[θ (t 1)τ β/2]. (14) As shown in Ref. [2] (Appendix A) for V (θ) = cos(θ), with straightforward extension to arbitrary V (θ), the expression (13) implies that for large t the dominant contribution of dψ β,t (θ)/dθ 2 to E t in Eq. (12) is k 2 ψβ (θ tτ β )d V β,t (θ)/dθ 2. This contribution can be calculated using Eq. (14) and a relation following from Eq. (4): ψ β (θ) 2 = 1 2π where C β (m) are correlations in momentum space, C β (m) exp(imθ), (15) m C β (m) = n Ψ(m + n + β) Ψ (n + β). (16) From Eqs. (14) and (15), one can write the Fourier expansion of k 2 ψ β (θ tτ β )d V β,t (θ)/dθ 2. After inserting this expansion in Eq. (12), we get where B(m, m ; t) = E t k m,m 0 mm V m V m B(m, m ; t), (17) dβ sin(mτ βt/2) sin(m τ β t/2) sin(mτ β /2) sin(m τ β /2) C β(m m )e i(m m)(t+1)τ β /2. (18) By a tedious but straightforward calculation, whose details will be given elsewhere [11], one can show that the asymptotic behavior of the quantity (18) for t is given by 7

8 B(m, m ; t) g(m,m t )l 0 1 2l 0 mm ( m + m m m ) C βr,g (m m ), (19) where g = g(m, m ) is the greatest common factor of ( m, m ) and β r,g = r/(gl 0 ) 1/2 mod(1), r = 0,..., gl 0 1. Rels. (17) and (19) imply a diffusive behavior of E t for large t, E t Dt. The diffusion coefficient D can be expressed, after some algebra, as the sum of two terms: D = D I + D II = k2 l 0 + 2k2 l 0 m>0 m=1 m =m+1 ml 0 1 m V m 2 r=0 C βr,m (0) g(m,m )l 0 1 mre (V m Vm ) r=0 r=0 Re [ C βr,g (m m ) ]. (20) For a potential containing harmonics V m of sufficiently high order m, the term D II [on the second line of Eq. (20)] can generally lead to a sensitive dependence of D on the correlations C β (m), which reflect the profile of the initial wave packet in momentum space by Eq. (16). This term vanishes in some cases, e.g., for V (θ) = cos(θ) and/or a uniform probability distribution ψ β (θ) 2 in Eq. (15). One is then left with only the first term (D I ), which is completely analogous to the diffusion coefficient (11) for the incoherent mixture; this is because C β (0) = n Ψ(n + β) 2 from Eq. (16) and this is analogous to f 0 (β) = n f(n+β). Clearly, all the resonant values β r,g of β appear in formulas (11) and (20). In fact, the momentum probability distribution (MPD) P(p, t) = Ψ t (p) 2 must exhibit ballistic peaks at p = n+β r,g and averaging P(p, t)p 2 around each peak gives its contribution to the diffusive behavior for t. We now study P(p, t) in the case of V (θ) = cos(θ)+γ cos(2θ), already considered above. From Rel. (4), Ψ t (n + β) are just the Fourier coefficients of 2πψ β,t (θ) and, for τ = 2πl 0, ψ β,t (θ) can be calculated using Rel. (13). We assume an initial wave packet Ψ(p) satisfying Ψ(p) = 1 for 0 p < 1 and Ψ(p) = 0 otherwise. This corresponds, by Rel. (4), to a uniform ψ β (θ), ψ β (θ) = (2π) 1/2 for all β. Rel. (13) then implies that 2πψβ,t (θ) = exp [ ik V β,t (θ) ]. The Fourier coefficients of exp [ ik V β,t (θ) ] can be easily calculated from (14) since V m = 0 for m > 2. We finally obtain 8

9 P(p = n + β, t) = m= i m J n 2m [k 1 (β, t)] J m [k 2 (β, t)γ] 2, (21) where J n are ordinary Bessel functions and k j (β, t) = k sin (jτ β t/2)/sin (jτ β /2), j = 1, 2. Let us consider some behaviors of (21) at the resonant values of p = n + β r,g, where β r,g (g = 1, 2) were determined above: (a) For β r,1 = r/l 0 1/2 mod(1), r = 0, 1,..., l 0 1, one has k 1 (β r,1, t) = k 2 (β r,1, t) = kt. Then, on the time scale of t kγ 1, one can approximate (21) at p = n + β r,1 by P(p, t) Jn 2 (kt), so that the most prominent g = 1 peaks of the MPD are those with n kt. In fact, for γ = 0, the expression (21) reduces exactly to P(p, t) = J 2 n [k 1(β, t)], which is the known result for the MPD in the case of V (θ) = cos(θ) and an initial Ψ(p) localized at p = β (a plane wave) [3]; due to the uniformity of ψ β (θ), Eq. (21) gives also the MPD for such an initial Ψ(p). (b) For β r,2 = r/(2l 0 ) 1/2 mod(1), r = 1, 3,..., 2l 0 1, one has k 1 (β r,2, t) = k or 0 for t odd or even, respectively, and k 2 (β r,2, t) = kt. Then, if k 1, J n 2m [k 1 (β r,2, t)] is relatively small for n 2m > 1 and one can approximate (21) at p = n + β r,2 by P(p, t) J[n/2] 2 (kγt), where [n/2] is the integer part of n/2. The new (g = 2) peaks thus start to emerge in a significant way when t > kγ 1 and, for n 2 kγt, their magnitude should be comparable to that of the g = 1 ones. All these behaviors are illustrated in Fig. 1 for τ = 2π (l 0 = 1), t = 100, k = 0.1, γ = 0 [Fig. 1(a)] and γ = 0.2 [Fig. 1(b), with 2 kγt = 4]. We have checked numerically for many values of (t, k, γ) that the MPD is robust under sufficiently small perturbations of τ, τ = 2π + ǫ, at least within the domains of p where the principal peaks above are found. As an example, compare Fig. 1(c) with Fig. 1(b). In conclusion, the results in this paper should provide new insights into the nature of the spectra and quantum dynamics of the KP under general QR conditions and for realistic, multi-harmonic potentials. At fixed τ/(2π) = l 0 /q 0 (l 0 and q 0 are coprime integers), QR may arise for arbitrary rational values of the quasimomentum β. A resonant β is characterized by an integer pair (r, g) which determines β through Eq. (7), with l = gl 0 and q = gq 0 assuming their minimal values. The QE spectrum consists of q = gq 0 bands. By slightly varying β on 9

10 the rationals, the corresponding value of g changes very erratically, implying similar changes in the QE spectrum. The case of the main QRs (q 0 = 1) is precisely described by the linear KR. In this case, β = β r,g is resonant provided that the potential contains harmonics whose order is a multiple of g. The diffusion coefficients (11) and (20) depend explicitly on all the resonant values β r,g of β, as these values correspond to ballistic peaks in the momentum probability distribution (see Fig. 1). Formula (20) implies a new phenomenon for general multi-harmonic potentials [thus excluding the standard case of V (x) = cos(x)]: A sensitive dependence of the diffusion coefficient on the specific localization features of the initial wave packet in momentum space. Assuming the robustness of our general results under small variations of τ, which was numerically verified for a two-harmonic potential, it may be possible to observe this phenomenon in experimental realizations of the system. This work was partially supported by the Israel Science Foundation (Grant No. 118/05). DLD acknowledges partial support from the Russian Ministry of Education and Science and the US Civilian Research and Development Foundation (CRDF BRHE Program, Grants Nos. VZ and Y2-P-10-01). 10

11 REFERENCES [1] R.M. Godun et al., Phys. Rev. A 62, (2000); M.B. d Arcy et al., Phys. Rev. E 64, (2001); S. Schlunk et al., Phys. Rev. Lett. 90, (2003); ibid., (2003); Z.-Y. Ma, M.B. d Arcy, and S.A. Gardiner, Phys. Rev. Lett. 93, (2004). [2] S. Fishman, I. Guarneri, and L. Rebuzzini, Phys. Rev. Lett. 89, (2002); J. Stat. Phys. 110, 911 (2003). [3] S. Wimberger, I. Guarneri, and S. Fishman, Nonlinearity 16, 1381 (2003); Phys. Rev. Lett. 92, (2004). [4] D.R. Grempel, R.E. Prange, and S. Fishman, Phys. Rev. A 29, 1639 (1984); D.L. Shepelyansky, Physica (Amsterdam) 28D, 103 (1987). [5] F.M. Izrailev, Phys. Rep. 196, 299 (1990), and references therein. [6] S.-J. Chang and K.-J. Shi, Phys. Rev. A 34, 7 (1986). [7] D.R. Grempel, S. Fishman, and R.E. Prange, Phys. Rev. Lett. 49, 833 (1982); R.E. Prange, D.R. Grempel, and S. Fishman, Phys. Rev. B 29, 6500 (1984). [8] M.V. Berry, Physica (Amsterdam) 10D, 369 (1984). [9] In Ref. [8], the potential (9) was assumed to be even (real V m = V m ). However, the results can be straightforwardly extended to the most general V (θ). [10] If τ/(2π) is rational but noninteger (q 0 > 1), resonant values of β with g > 1 may exist also for V (θ) = cos(θ). For example, in the case of τ/(2π) = 1/2 and V (θ) = cos(θ), we found numerically that for g = 2 all the q = gq 0 = 4 QE bands for β = 1/2 [r = 1 and l = gl 0 = 2 in Eq. (7)] have nonzero width, implying QR. [11] I. Dana and D.L. Dorofeev, in preparation. 11

12 FIGURES FIG. 1. Momentum probability distribution (21) for t = 100, k = 0.1, and (a) τ = 2π, γ = 0; (b) τ = 2π, γ = 0.2; (c) τ = 2π + ǫ [ǫ = π( 5 1)/ ], γ = 0.2. Note the new resonant peaks at integer p = n for γ = 0.2 [in (b) and (c)], adding to the peaks at p = n + 1/2, the only ones for γ = 0 [in (a)]. 12

13 P(p) Figure 1 (a) (b) (c) p