Phase-Space Description of the Coherent State Dynamics in a Small One-Dimensional System
Total Page:16
File Type:pdf, Size:1020Kb
Open Phys. 2016; 14:354–359 Research Article Open Access Urszula Kaczor, Bogusław Klimas, Dominik Szydłowski, Maciej Wołoszyn, and Bartłomiej J. Spisak* Phase-space description of the coherent state dynamics in a small one-dimensional system DOI 10.1515/phys-2016-0036 space variables. In this case, the product of any functions Received Jun 21, 2016; accepted Jul 28, 2016 on the phase space is noncommutative. The classical alge- bra of observables (commutative) is recovered in the limit Abstract: The Wigner-Moyal approach is applied to investi- of the reduced Planck constant approaching zero ( ! 0). gate the dynamics of the Gaussian wave packet moving in a ~ Hence the presented approach is called the phase space double-well potential in the ‘Mexican hat’ form. Quantum quantum mechanics, or sometimes is referred to as the de- trajectories in the phase space are computed for different formation quantization [6–8]. Currently, the deformation kinetic energies of the initial wave packet in the Wigner theory has found applications in many fields of modern form. The results are compared with the classical trajecto- physics, such as quantum gravity, string and M-theory, nu- ries. Some additional information on the dynamics of the clear physics, quantum optics, condensed matter physics, wave packet in the phase space is extracted from the analy- and quantum field theory. sis of the cross-correlation of the Wigner distribution func- In the present contribution, we examine the dynam- tion with itself at different points in time. ics of the initially coherent state in the double-well poten- Keywords: Wigner distribution, wave packet, quantum tial using the phase space formulation of the quantum me- trajectory chanics. This autonomous path to quantization of the clas- sical theory is based on the Wigner distribution function PACS: 03.65.Sq,03.65.Ca (WDF), ϱW , which is defined by the Weyl transform of the density operator [4, 9, 10], Z (︂ )︂ 1 1 1 −(i/~)pX 1 Introduction ϱW (x, p; t) = dX ρ x − X, x + X; t e . 2π~ 2 2 (1) Description of the quantum dynamical processes in the The WDF is normalized and real, but not always posi- phase space can be realized on the Wigner-Moyal alge- tive function of position and momentum over the phase bra which is the non-commutative algebra of observables space [11]. The negative values of the WDF in some regions based on the symplectic symmetry of the phase-space [1]. of the space are a consequence of the Heisenberg’s uncer- In this approach, the state of the quantum dynamical sys- tainty principle and they can be regarded as an indicator tem cannot be longer represented by the Dirac measure of the non-classicality of the state [12]. Despite the fact that localized at a point in the phase space, but rather by an the WDF can be negative, the expectation value of any dy- appropriate quasi-distribution function [2–4], because the namical variable can be computed with respect to the WDF Heisenberg’s uncertainty principle prohibits precise deter- in the same way as the average value in the classical sta- mination of momentum and position at the same time [5]. tistical mechanics, namely [4, 11] The uncertainty also tampers with the algebra of observ- Z ables which are represented by functions of the phase- hF(t)i = dxdp FW (x, p)ϱW (x, p; t), (2) where FW (x, p) is the Weyl symbol of quantum- *Corresponding Author: Bartłomiej J. Spisak: AGH University mechanical operator F in the position representation [13], of Science and Technology, Faculty of Physics and Applied Com- puter Science, Al. A. Mickiewicza 30, 30-059 Krakow, Poland; Email: [email protected] Z (︂ )︂ 1 1 1 −(i/~)pX Urszula Kaczor, Bogusław Klimas, Dominik Szydłowski, Maciej FW (x, p) = dX F x − X, x + X e . (3) 2π~ 2 2 Wołoszyn: AGH University of Science and Technology, Faculty of Physics and Applied Computer Science, Al. A. Mickiewicza 30, 30- 059 Krakow, Poland © 2016 U. Kaczor et al., published by De Gruyter Open. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivs 3.0 License. Phase-space description of the coherent state dynamics Ë 355 2 Theory 100 ω2=2, λ=0.1 ω2=1.5, λ=0.02 The unitary time evolution of the WDF is generated by 80 ω2=3, λ=0.03 ω2 λ the Moyal equation which exhibits the non-local nature =4, =0.08 of the quantum dynamics. The equation of motion has a 60 Liouville-like form, i.e., 40 ∂ ϱW (x, p; t) = LϱW (x, p; t), (4) ∂t potential (a.u.) 20 where the right hand side of the equation (4) is expressed by the Moyal bracket. Its explicit form is given by the fol- 0 lowing formula [1, 8]: −20 −8 −6 −4 −2 0 2 4 6 8 LϱW (x, p; t) = fHW (x), ϱW (x, p; t)g* {︂ [︂ (︂−! − −! −)︂]︂ position (a.u.) 1 ~ ∂ ∂ ∂ ∂ = HW (x, p) exp − i~ 2 ∂p ∂x ∂x ∂p }︂ Figure 1: The potential used in the calculations. All values are given in atomic units. × ϱW (x, p; t) = HW (x, p) * ϱW (x, p; t), (5) form (presented in Fig. 1). The system may serve as a pro- where HW (x, p) is the Weyl symbol of the small system totype of an electronic nanosystem exhibiting nonequilib- Hamiltonian, and the arrows indicate in which direction rium dynamics in the presence of the external gradients or the derivatives act. We note that the lowest order of the fields [22]. The Weyl symbol of the small system Hamilto- series expansion of the time evolution generator for the nian takes the form 2 Moyal equation with respect to the deformation parameter p 1 2 2 4 HW = − mω x + λx , (7) ~/2 corresponds to the Poisson bracket, and therefore the 2m 2 Moyal equation is reduced to the exact form of the classi- where m is the mass of the electron, while ω and λ are the cal Liouville equation [1]. Alternatively the classical limit parameters of the potential. For the present calculations of the dynamics generated by the Moyal equation can be their values are taken to be equal ω = 2 a.u. (atomic units deduced on the basis of the Egorov’s theorem [14]. Follow- of e = ~ = m = 1) and λ = 0.08 a.u. ing this theorem, the quantum dynamics is reduced to the classical as the semiclassical parameter ϵ which effectiv- elly replaces ~ determined by the physical scales involved 3 Results in the system goes to zero corresponding to ϵ2 [15]. Some general similarities between the propagation of the classical and coherent states in the phase space based We have determined the propagation of the WDF in on the careful studies of the semiclassical limit [16] al- the small system by means of the method presented in low incorporating the split-operator technique to deter- Ref. [23]. For the present calculations the followingp val- δ mine the time evolution of the WDF [17–19]. In the present ues of the wave packet are assumed: x = 1/ 2 a.u., x p calculations, the initial condition for the Moyal equation 0 = −3.5 a.u., and three different initial momenta: 0 = is taken in the coherent state form which is represented −10, −3.5 and 5 a.u., which correspond to different kinetic by the Wigner form of the Gaussian wave packet centered energies of the coherent state. The phase-space calcula- tions are performed on the computational grid with Nx = around some point (x0, p0) in the phase space [20], 512 mesh points for the x-coordinate, and Np = 512 mesh {︂ 2 2 2 }︂ 1 2δx (p − p0) (x − x0) x ϱW (x, p; 0) = exp − − , (6) points for the -component of momentum. The quantum π 2 δ2 ~ ~ 2 x trajectories are determined by simultaneous calculations where δx is the initial half-width of the wave packet. The of the expectation values of the position and momentum form represents the most classical quantum state permit- variables according to the formula (2) with the WDF being ted by the uncertainty principle [21]. a solution of the Moyal equation (4). On the other hand, We apply this formalism to the description of the co- the classical trajectories are computed on the basis of the herent state dynamics in the small system which is mod- Hamilton’s equations [24]. The results of performed com- eled by the double well potential in the “Mexican hat” puter simulations are presented in Figs. 2, 3, and 4, which 356 Ë U. Kaczor et al. display not only the snapshots of the WDF at different times, but also the quantum and classical trajectories. A comparison of the both types of the trajectories seems to be extremely informative, because it displays the main differences between the dynamics of the classical and quantum states in the phase space, generated by the Poisson and Moyal brackets, respectively. In the case of the classical dynamics of the particle which moves under the influence of a conservative force (in the present case, the force is F = mω2x − 4λx3), the trajectory does not have in- tersections and is always closed, although its geometrical center as well as the shape strongly depends on the initial momentum. Due to the knowledge of the classical trajecto- ries, the Liouville formulation of the classical dynamics in the phase space [25] can be simply reconstructed, namely the Dirac measure, ρ(x, p, t), (︀ )︀ (︀ )︀ ρ(x, p, t) = δ x − x(t) δ p − p(t) , (8) which imitates the distribution function for the single par- ticle, moves along the classical trajectory.