A novel derivation of quantum propagator useful for time-dependent trapping and control Fardin Kheirandish1, ∗ 1Department of Physics, University of Kurdistan, P.O.Box 66177-15175, Sanandaj, Iran Abstract A novel derivation of quantum propagator of a system described by a general quadratic Lagrangian is presented in the framework of Heisenberg equations of motion. The general corresponding density matrix is obtained for a derived quantum harmonic oscillator and a particle confined in a one dimensional Paul trap. Total mean energy, work and absorbed heat, Wigner function and excitation probabilities are found explicitly. The method presented here is based on the Heisenberg representation of position and momentum operators and can be generalized to a system consisting of a set of linearly interacting harmonic oscillators straightforwardly. PACS numbers: 37.10.Ty, 03.65.Ca, 03.65.-w, 05.30.-d Keywords: Trapped particle; Propagator; Green's function; Density matrix arXiv:1612.07674v3 [quant-ph] 14 Jun 2018 ∗Electronic address: [email protected] 1 I. INTRODUCTION The propagators are extensively used in many branches of Physics like quantum electrodynamics, quantum statistical mechanics, condensed matter Physics, polymer Physics, astrophysics and economics [1]. They are the main ingredients of theoretical basis of science and engineering. In non-relativistic quantum mechanics, propagators can be calculated in three different approaches [2]: (i) Using the explicit form of the energy eigenstates satisfying predefined boundary conditions. (ii) Solving the Heisenberg equations of motion for the problem in question. (iii) Using path integral techniques. The third approach from a pedagogical standpoint is interesting since it points out to the analogy with classical mechanics trough action principle. Due to the formal equivalence of the time-evolution operator in quantum mechanics and the density operator in quantum statistical mechanics, it turns out that the propagator formalism is also useful in calculating partition functions or thermodynamical quantities in quantum statistical mechanics. In non-relativistic quantum mechanics the propagator gives the probability amplitude for a particle to travel from one spatial point at one time to another spatial point at a later time. This means that, if a system has Hamiltonian H, then the appropriate propagator is a function [3] 1 G(x; t; x0; t0) = Θ(t − t0)K(x; t; x0; t0); (1) i~ satisfying @ (i − H )G(x; t; x0; t0) = δ(x − x0)Θ(t − t0); (2) ~@t x where Θ(t) is the Heaviside step function and K(x; t; x0; t0) is the kernel of the differential operator in question, often referred to as the propagator instead of G in this context. The propagator can also be written as K(x; t; x0; t0) = hxjU^(t; t0)jx0i; (3) where U^(t; t0) is the unitary time-evolution operator for the system taking states at time t to states at time t0. The general quadratic Lagrangian Eq.(5) (see Sec.II), can be transformed to the Lagrangian (7) with the corre- sponding Hamiltonian Eq.(9). As a special case, by setting c(t) = m!2(a − 2q cos(Ωt)) and e(t) = 0, we find the Hamiltonian p2 1 H = + m!2(a − 2q cos(Ωt)) x2; (4) 2m 2 describing the quantum dynamics of a particle in a Paul trap [4{6]. Here we have obtained the quantum propagator for a system described by the general quadratic Hamiltonian Eq.(9) which can have applications in quantum control and information processing. A novel derivation is presented in the framework of Heisenberg equations of motion. From the connection between propagator and density matrix, we find the corresponding density matrix for a derived quantum harmonic oscillator and a particle confined in a one dimensional Paul trap. In the following, we will find the total mean energy, work and absorbed heat, Wigner function and excitation probabilities. In the presence of external source (e(t) 6= 0), we propose to express the rate of work done by the external source or oscillator and also the absorbed heat by the oscillator, trough expressions that are different from the usual definitions of those quantities in literature. The method can be generalized to a system consisting of a set of linearly interacting harmonic oscillators, straightforwardly. II. LAGRANGIAN The general form of a classical quadratic Lagrangian is given by 1 1 L = mx_ 2 + a (t)xx_ + a (t)x2 + a (t)_x + a (t)x; (5) 2 1 2 2 3 4 regarding the fact that a total time derivative can be removed from Lagrangian, we rewrite the second and fourth terms as d 1 1 a (t)xx_ = a (t)x2 − a_ (t)x2; 1 dt 2 1 2 1 d a (t)_x = a (t)x − a_ (t)x: (6) 3 dt 3 3 2 By inserting Eqs.(6) into Eq.(5) and defining new functions c(t) =a _ 1(t) − a2(t) and e(t) =a _ 3(t) − a4(t), we can rewrite the Lagrangian given in Eq.(5) as [7] 1 1 L = mx_ 2 − c(t)x2 − e(t)x: (7) 2 2 The conjugate momentum corresponding to the position x is defined by @L p = = mx;_ (8) @x_ and the corresponding Hamiltonian H = px_ − L, is p2 1 H = + c(t)x2 + e(t)x: (9) 2m 2 The system is quantized by promoting classical variables (x; p) to quantum partners (^x; p^) and imposing canonical quantization rule [^x; p^] = i~. From Heisenberg equations for position and momentum, we find c(t) e(t) x^¨ + x^ = − : (10) m m Equation (10) can be solved formally as Z t 1 0 0 0 x^(t) =x ^H (t) − dt G(t; t ) e(t ); (11) m 0 where G(t; t0) is the Green's function of Eq.(10) satisfying d2 1 + c(t) G(t; t0) = δ(t − t0); (12) dt2 m andx ^H (t) is the homogeneous solution c(t) x^¨ (t) + x^ (t) = 0: (13) H m H Let f(t) and g(t) be classical independent solutions of Eq.(13), then c(t) f¨(t) + f(t) = 0; m c(t) g¨(t) + g(t) = 0; m x^H (t) = Af^ (t) + Bg^ (t); (14) where A^ and B^ are operator constants to be determined later. We can find the retarded Green's function in terms of the independent functions f(t) and g(t). The retarded Green's function is defined by 0; t < t0; G(t; t0) = (15) af(t) + bg(t); t > t0; From the continuity and also the discontinuity of the time derivative of Green's function at t = t0, we have af(t0) + bg(t0) = 0; dG dG − = 1; (16) dt t=t0+ dt t=t0− by inserting Eq.(15) into Eq.(16), we find the retarded Green's function as g(t0)f(t) − f(t0)g(t) G(t; t0) = θ(t − t0); (17) W (t0) 3 where Θ(t − t0) is the Heaviside step function and W (t) is the Wronskian defined by g(t) f(t) W (t) = det : (18) g_(t) f_(t) From Eqs.(8,11,17), we will find x^(0) =x ^H (0) ≡ x;^ (19) p^(0) = mx_(0) = mx_ H (0) ≡ p:^ (20) The operator constants A;^ B^ can be expressed in terms of the initial conditions or Schr¨odingeroperatorsx ^,p ^ as g(0) g_(0) A^ = p^ − x;^ mW (0) W (0) f_(0) f(0) B^ = x^ − p:^ (21) W (0) mW (0) Now by inserting Eqs.(21) into Eq.(14) and making use of Eq.(11), we find the position operator as g(t)f_(0) − f(t)_g(0) f(t)g(0) − g(t)f(0) p^ 1 Z t x^(t) = x^ + − dt0 G(t; t0)e(t0): (22) W (0) W (0) m m 0 III. PROPAGATOR In Heisenberg picture, the time evolution of the position operator is given by x^(t) = U^ y(t)^x(0)U^(t); (23) or equivalently U^(t)^x(t) =x ^(0)U^(t): (24) Inserting Eq.(22) into Eq.(24) and taking the matrix elements in position space, we find β(t) γ(t) hxjU(t)[α(t)x + p − ]jx0i = xhxjU(t)jx0i; (25) m m where g(t)f_(0) − f(t)_g(0) α(t) = ; (26) W (0) f(t)g(0) − g(t)f(0) β(t) = ; (27) W (0) Z t γ(t) = dt0 G(t; t0) e(t0): (28) 0 From Eq.(25) and the definition K(x; t; x0; 0) = hxjU(t)jx0i, we have β(t) @ γ(t) x0α(t) + i − K(x; t; x0; 0) = xK(x; t; x0; 0); (29) ~ m @x0 m 0 0 0 0 @ 0 where we usedx ^jx i = x jx i,p ^jx i = i~ @x0 jx i. Now Eq.(29) can be rewritten as @ i m γ(t) ln K(x; t; x0; 0) = − [x − x0α(t) + ]: (30) @x0 ~ β(t) m 4 Since the right side of Eq.(30) is linear in x0, we can write 1 ln K(x; t; x0; 0) = C + C x0 + C x02; (31) 0 1 2 2 where the unknown coefficients C0, C1 and C2 may depend on x and t. Inserting Eq.(31) into Eq.(30), leads to 0 im 0 γ (C1 + C2x ) = − x − α(t) x + ; (32) ~β(t) m and by comparing the coefficients of x0 on both sides we obtain im γ imα C1 = − (x + );C2 = : (33) ~β m ~β Therefore, im 0 imα 02 iγ 0 0 C0(x;t) − xx x − x K(x; t; x ; 0) = e e ~β e 2~β e ~β : (34) Now for the momentum operator we have p^(t) = mx^_ = mα_ (t)^x + β_(t)^p − γ_ (t) = U y(t)^p U(t); (35) so U^(t)^p(t) =p ^U^(t), leading to hxjU^(t)[mα_ (t)^x + β_(t)^p − γ_ (t)]jx0i = hxjp^U^(t)jx0i; (36) or equivalently @ @ (mαx_ 0 + i β_ − γ_ )K(x; t; x0; 0) = −i K(x; t : x0; 0): (37) ~ @x0 ~@x By inserting Eq.(34) into Eq.(37), we find @C (x; t) i m mαβ_ mβ_ γβ_ − γβ_ 0 = [(mα_ + − )x0 + x + ]: (38) @x ~ β β β β Since the Wronskian does not depend on time (W (t) = W (0)), we can easily proof the identity αβ_ − αβ_ = 1, and accordingly, Eq.(38) can be solved as im 2 i _ 2 _ C0(x; t) = βx + (γβ − γβ_ )x + (t); (39) 2~β m ~ where (t) is an unknown function of time.