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]uok.ac.ir
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) = hx|Uˆ(t, t0)|x0i, (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ˆ †(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) hx|U(t)[α(t)x + p − ]|x0i = xhx|U(t)|x0i, (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) = hx|U(t)|x0i, 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 ˆ|x i = x |x i,p ˆ|x i = i~ ∂x0 |x 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 †(t)ˆp U(t), (35) so Uˆ(t)ˆp(t) =p ˆUˆ(t), leading to
hx|Uˆ(t)[mα˙ (t)ˆx + β˙(t)ˆp − γ˙ (t)]|x0i = hx|pˆUˆ(t)|x0i, (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 2 ˙ i C0(x, t) = βx + (γβ − γβ˙ )x + ψ(t), (39) 2~β m ~ where ψ(t) is an unknown function of time. The propagator now can be written as
2γ 0 i ψ(t) im βx˙ 2+αx02−2xx0+ 2 (γβ˙−γβ˙ )x− x0 K(x, t; x , 0) = e ~ e 2~β m m . (40)
By making use of the identity
† 0 0 0 hx|Uˆ(t)Uˆ (t)|x i = hx|I|x i = δ(x − x ), (41) or equivalently Z dx00 hx|Uˆ(t)|x00ihx00|Uˆ †(t)|x0i = δ(x − x0), (42) we find easily r i ψ(t) m i ϕ(t) e ~ = e ~ , (43) 2π~β
5 where ϕ(t) is a real function of time to be determined. Up to now, the quantum propagator has the following form r 0 m i ϕ(t) im αx02+βx˙ 2−2(αβ˙−αβ˙ )xx0 i (γβ˙−γβ˙ )x−γx0 K(x, t; x , 0) = e ~ e 2~β e ~β . (44) 2π~β Another well known property of the propagator is lim K(x, t; x0, 0) = δ(x − x0), (45) t→0 from Eqs.(26,27,28) we have α(t) → 1, t→0 γ(t) → 0, t→0 β(t) → t + o(t2), t→0 (46) where β(t) given in Eq.(27), has been expanded by a Taylor series around t = 0 up to the first order. Note that α(t) and γ(t) in Eq.(44) can be simply replaced by their respective limiting values but β(t) carrying the singularity, should be replaced with its asymptotic value, namely t. Therefore, i i ϕ(0) 0 2 e ~ (x − x ) lim K(x, t; x0, 0) = lim √ e2~t , t→0 t→0 2π~t = δ(x − x0). (47) Comparing this result with the well known identity
1 − A (x−x0)2 0 lim(A/πt) 2 e t = δ(x − x ), (48) t→0 √ i ϕ(0) we obtain e ~ = 1/ i. If we redefine λ(t) = ϕ(t) − ϕ(0), then for the quantum propagator we find
r 2γ 0 m i λ(t) im βx˙ 2+αx02−2xx0+ 2 (γβ˙−γβ˙ )x− x0 K(x, t; x , 0) = e ~ e 2~β m m . (49) 2πi~β To determine λ(t), we know that the propagator K(x, t; 0, 0) satisfies the Schr¨odingerequation 2 1 i ∂ K(x, t; 0, 0) = − ~ ∂2 + c(t) + e(t)x K(x, t; 0, 0), (50) ~ t 2m x 2 taking temporal and spatial differentiations, and setting x = 0, we find easily
1 γβ˙ − γβ˙ 2 λ˙ (t) = − . (51) 2m β Finally, the quantum propagator corresponding to the quadratic Lagrangian Eq.(7) is
r γβ˙−γβ˙ 2γ 0 m − i R t dt0 ( )2 im βx˙ 2+αx02−2xx0+ 2 (γβ˙−γβ˙ )x− x0 K(x, t; x , 0) = e 2m~ 0 β e 2~β m m . (52) 2πi~β
IV. A DIFFERENTIAL EQUATION FOR β(t)
Setting t0 = 0 in Eq.(17) and comparing the result with Eq.(27), we have G(t) = β(t) Θ(t). (53) From Eq.(27) it is clear that β(t) satisfies the equation d2 1 + c(t) β(t) = 0, (54) dt2 m
6 with the initial conditions
dβ(t) β(0) = 0, = 1, (55) dt t=0 which is the same result obtained by Schulman [7] using a different approach. Example 1. For a harmonic oscillator in an external classical source (e(t) 6= 0), we set c(t) = mω2 in Eq.(9) and find f(t) = cos(ωt), g(t) = sin(ωt), α(t) = cos(ωt), sin(ωt) β(t) = , ω sin[ω(t − t0)] G(t, t0) = θ(t − t0), ω Z t sin[ω(t − t0)] γ(t) = dt0 e(t0), 0 ω therefore, r 0 mω imω (x2+x02) cos(ωt)−2xx0 K(x, t; x , 0) = e 2~ sin(ωt) 2πi~ sin(ωt) − i x R t dt0 sin(ωt0)e(t0)+x0 R t dt0 sin ω(t−t0)e(t0) × e ~ sin(ωt) 0 0 0 − i R t dt0 R t ds e(t0)e(s) sin ω(t−t0) sin(ωs) × e m~ω sin(ωt) 0 0 . (56) as expected [7].
V. DENSITY MATRIX
The evolution of the density matrix ρ(t) corresponding to a system described by the Hamiltonian Eq.(9) is given by ρ(t) = U(t)ρ(0)U †(t), (57) where U(t) is the evolution operator connected to the propagator or kernel trough K(x, t; x0, 0) = hx|U(t)|x0i. There- fore, in position basis we can write ZZ ∞ 0 ∗ 0 ρ(x, x ; t) = dx1dx2 K(x, t; x1, 0) ρ(x1, x2; 0) K (x , t; x2, 0). (58) −∞ By inserting the general form Eq.(52) for the propagator into Eq.(58) we find
0 mω im [β˙(x2−x02)+ 2 (γβ˙−γβ˙ )(x−x0)] ρ(x, x ; t) = e 2~β m 2π~|β(t)| ZZ im 2 γ 2 0 γ [αx −2x1(x+ )−αx +2x2(x + )] × dx1dx2 e 2~β 1 m 2 m ρ(x1, x2; 0). (59)
To proceed, we assume the following gaussian initial state r λ − λ (x2+x2) ρ(x , x ; 0) = e 2 1 2 , (60) 1 2 π leading to s 0 λ imα(ζ−1) (x2−x02) [ iαγ(ζ−1) − iγ˙ ](x−x0) ρ(x, x ; t) = e 2~βζ e ~βζ ~ πζ(t)
− λ [(x+ γ )2+(x0+ γ )2] × e 2ζ m m , (61)
7 where for notational simplicity we have defined 1 ζ(t) = 2 2 2 . (62) 2 λ ~ β (t) α (t) + m2 The diagonal elements of the density matrix give the probability density to find the particle at position x at time t s λ − λ (x+ γ )2 P (x, t) = ρ(x, x; t) = e ζ(t) m . (63) πζ(t)
Note that Eq.(63) is a general result obtained for a system governed by the Hamiltonian Eq.(9). For a harmonic oscillator in an external classical field e(t), by setting λ = mω/~ we find ζ(t) = 1, and r 0 λ −iγ˙ (x−x0) − λ [(x+ γ )2+(x0+ γ )2] ρ(x, x ; t) = e ~ e 2 m m . (64) π If we rewrite Eq.(64) as
ρ(x, x0; t) = hx|ρ|x0i = hx|ψ(t)ihψ(t)|x0i, (65) we deduce that the ground state Eq.(60) has evolved to the pure state