Quantum chaotic attractor in a dissipative system W. Vincent Liu and William C. Schieve Department of Physics and Center for Statistical Mechanics & Complex System, The University of Texas, Austin, Texas 78712 fluence functional approach of Feynman and Vernon [9] A dissipative quantum system is treated here by coupling and its application to Brownian motion [10], the quantum it with a heat bath of harmonic oscillators. Through quan- Langevin equations [11], the master equations [5–7,11], tum Langevin equations and Ehrenfest’s theorem, we estab- lish explicitly the quantum Duffing equations with a double- etc. [12]. Among these, the quantum Langevin formal- well potential chosen. A quantum noise term appears the only ism is more interesting for us since the dynamics of quan- driving force in dynamics. Numerical studies show that the tum operators for the system can be explicitly given in a chaotic attractor exists in this system while chaos is certainly fashion quite similar to the classical one. Most of above forbidden in the classical counterpart. theories considered a generic model of the system-bath interaction. Recently, Pattanayak and Schieve [3] pro- PACS numbers: 05.45.+b,03.65.Sq,05.40.+j,42.50.Lc posed a derivation of semiclassical dynamics directly from Heisenberg equations of motion via Ehrenfest’s theorem where ostensible (apparent) quantum chaos was found Quantum chaos of Hamiltonian systems has been stud- for the conservative (Hamiltonian) system with a double- ied extensively [1–4]. By contrast, very little work has well potential. In this paper, we will study a dissipative been done in looking at quantum chaos of dissipative sys- quantum model similar to that of Refs. [9–11] and write tems. Many quantum mechanical systems (e.g., SQUID explicitly the quantum Duffing equations for a system with Josephson junction, an atom in a cavity of electro- of double-well potential by means of Langevin formal- magnetic fields, NMR quantum measurements), however, ism and Ehrenfest’s theorem. Numerical results are now are neither isolated nor Hamiltonian. These interact possible to obtain and show that chaotic attractors ro- with their environment and thus are open in general and bustly exist in this model whereas chaos is impossible noisy and dissipative. Dissipation is relatively difficult for the corresponding classical Duffing system due to the to treat in a quantum system since it seems inevitable to absence of an external driving force. deal with stochastic processes via, for instance, the com- To do this, let us consider a particle of unit mass mov- monly used master equation that is extremely difficult to ing in a one-dimensional time-independent bounded po- solve numerically. Recently, Spiller and Ralph [5] stud- tential V (Qˆ) with its Hamiltonian, H = 1 Pˆ2 + V (Qˆ). ied a damped and driven non-linear oscillator by using sys 2 To introduce dissipation, we may take the system H the quantum state diffusion model based on the assumed sys linearly interacting with an external “heat bath” of many master equation of Lindblad form [6]. They simulated degrees of freedom, which here is assumed an assembly of the behavior of one member of the ensemble with a sin- harmonic oscillators [9–11]. The familiar example of this gle environment operator and found that the quantum model is a system of an atom interacting with a bath of noise “kicks” the motion between the chaotic and the equilibrium photons. Then the complete Hamiltonian is periodic behavior and so smears out any fractal struc- [9–11], ture. In general, the quantum state diffusion method arXiv:chao-dyn/9703012v1 21 Mar 1997 often demands an approximation of the effects of the Pˆ2 pˆ2 1 C Qˆ complicated environment by simple operators and always H = + V (Qˆ)+ n + mω2 (ˆx − n )2 , (1) 2 2m 2 n n mω2 has a problem regarding the dimension of the environ- n ( n ) X ment [7]. Brun [8] has tried to derive the quantum ver- sion of the forced and damped Duffing oscillator through in which we have assumed that all harmonic oscillators the decoherence approach which involves the path inte- possess the same mass m but may have different fre- gral. As he noticed, performing calculations in the low- quency ωn and that Cn is the coupling constant be- temperature limit is extremely difficult. For the high tween the system and the nth oscillator. The equal times commutation relations implicit in (1) are [Q,ˆ Pˆ]= temperature limit, Brun also suggested that the quantum ˆ ˆ ˆ ˆ noise “smears out” the quantum maps based on Wigner i¯h, [Q, xˆn] = [Q, pˆn] = [P, xˆn] = [P, pˆn]=0, [ˆxl, xˆn] = distributions and that numerical computation cannot be [ˆpl, pˆn]=0and[ˆxl, pˆn]= i¯hδln. Now it is straightforward efficiently done because of the necessity of enumerating to write down the Heisenberg equations of motion for all the possible histories and elements of the decoherence the complete system. Then we may apply the standard functional. procedure of Ref. [11] to derive the quantum Langevin The dissipative quantum system has been studied by equations. Assuming a continuous frequency distribution means of several different theories, for instance, the in- g(ω) of harmonic oscillators and using the first Markov approximation 1 2 ˙ g(ω)C (ω) 2γ hQˆi = hPˆi, (6a) 2 ≡ (2) mω π ˙ hPˆi = −hV ′(Qˆ)i− γhPˆi + η(t). (6b) where γ is assumed constant, we then have ˙ We now expand the equations around the centroid by Qˆ = P,ˆ (3a) ˆ (n) ˜n (n) using the identity F (Q)= n F Q /n!, where F = ˙ ′ n n Pˆ = −V (Qˆ) − γPˆ + ξˆ(t), (3b) ∂ F (Qˆ)/∂Qˆ | ˆ ˆ , and for a double-well potential of Q=hQi P V (Qˆ)= − 1 aQˆ2 + 1 bQˆ4 obtain the closed system of the where ξˆ(t) is the quantum noise operator due to the 2 4 stochastic differential equations heat bath and γ defined by (2) clearly represents the constant damping coefficient. Suppose that the system Q˙ = P, (7a) and the bath are initially independent (at t → −∞) so P˙ = aQ − bQ3 − 3bQ¯hµ − γP + η(t), (7b) that the complete density operator can be written into ρ = ρsys ⊗ ρb, and that the bath is initially thermal, µ˙ = α, (7c) ρ ∼ exp(−H /k T ). Then, ξˆ(t) has the properties (see, 1+ α2 b b B α˙ = +2µ(a − 3bQ2) − 6b¯hµ2 − γα. (7d) for example, [11]) 2µ hξˆ(t)i =0, (4a)Here we have written Q, P for hQˆi, hPˆi. We remind the ∞ ′ 2γ¯h ¯hω ′ reader that all physical observables must be obtained by h[ξˆ(t ), ξˆ(t)]+i = dωω coth( )cos ω(t − t), (4b) π 0 2kBT averaging over stochastic noise η(t). When damping con- Z stant γ is set to zero, the system is Hamiltonian with 2 2 in which [···]+ denotes an anticommutator. The average P π a 2 2 1 b 4 4 Hextended = + − (Q + ρ )+ 2 + (Q +3ρ + h· · ·i is over all bath variables. The operator nature of 2 2 2 8ρ 4 6¯hQ2ρ2), for which µ → ρ2 and α → 2ρπ in (7). The ξˆ(t) can be reduced by using the strategy of the adjoint above equations of motion can be reduced to those of commutative representation [13]. We can define a new Ref. [3], in which Hamiltonian semiquantal chaos was re- ′ 1 ˆ ′ operator η(t) by η(t)ρ(t ) ≡ 2 [ξ(t),ρ(t )]+ for all t and ′ ′ ported. t , which yields [η(t), η(t )] = 0. This means that η(t) is Also, in the classical limit ¯h → 0 the first two equations a c-number function of time. Let us replace the opera- ˆ of (7) decouple from the fluctuation variables (µ, α) and tor ξ(t) by the c-number η(t) in (3). Its 1- and 2-point the noise term vanishes as T → 0 [19]. The well-known correlation functions are given by (4). In the following classical Duffing equations without external driving are we will restrict ourselves in low temperature limit (not recovered (for T → 0) discussed in [8]), T → 0, in which quantum effects are most important. In this limit the noise η(t) is given by Q˙ = P, (8a) ˙ 3 hη(t)i =0, (5a) P = aQ − bQ − γP. (8b) γ¯h ∞ hη(t)η(t′)i = dωω cos ω(t′ − t) This shows that (7) are indeed the quantum analog of π Z0 Duffing equations with the presence of a quantum noise γ¯h 1 = − . (5b) term, which serves as a Langevin driving force. Now a π (t − t′)2 question perhaps rises upon how effective our equations (7) could be to give the description of exact quantum dy- To obtain this, an exponential cutoff (e−ǫω) was used, namics. The procedure of using Gaussian wavepacket to letting ǫ → 0 after integration [14]. describe the motion of a particle may be subject to large While the dependence of bath operators has been elim- errors and will even break down after a certain time scale inated in it, the quantum Langevin equations (3) are [20]. This has yet to be decided. It has been argued, still operator equations and are as difficult as that in for example, by Heller [17] who has extensively studied the deterministic quantum mechanics. To make fur- quantum chaos without noise that the generalized Gaus- ther progress, we make an additional assumption that sian approximation works extremely well. Also, recently the wave packet of the system can be described by Ashkenazy et al [4] have computed the time development the squeezed coherent state [15]. Then we have the of the wave function in the presence of a potential bar- relations [3] hQ˜2mi = (2m)!(¯hµ)m/m!2m, hQ˜2m+1i = rier in a bounded well for a long time.