# On the Canonical Transformation of Time-Dependent Harmonic Oscillator

Authors: Akpan N. Ikot; Louis E. Akpabio; Ita O. Akpan; Michael I. Umo; Oladunjoye A. Awoga; Alice B. Nyong
Journal: Physics Research International
 (2010)

---

## Abstract

We performed a two-variable canonical transformation on the time momentum operator, and without loss of generality we carried out a three-variable transformation on the coordinate and momentum space operators to trivialize the Hamiltonian operator of the system. Fortunately, this operation separates the time-coordinate and space coordinate naturally, and the wave function of the time-dependent Harmonic Oscillator is evaluated via the generator.


---

## Body

## 1. Introduction

The method of Canonical transformations (CT) has proved to be a fruitful approach in treating quantum systems [1]. It is often used in describing systems with Hamiltonian that are quadratic either in co-ordinate and momentum or equivalently in Boson creation and annihilation operators. One major advantage of the method of (CT) consists in reducing the Hamiltonian H of the system S to a Hamiltonian H′ of some simple system S′ with known solution, that is, H=CH′C-1. These two systems are canonically equivalent since their Hamiltonian is related by means of a CT. However, analytical   approaches   to isolate conserved quantities for a given physical system are a major objective in the realm of Hamiltonian theory [2]. For autonomous, where the Hamiltonian is time independent, one conserved quantity is immediately found: the Hamiltonian which represents the total energy of the system is a constant of motion.Rather unfortunately the Hamiltonian of most real physical systems is explicitly time-dependent and does not provide direct conserved quantity.  In classical mechanics, the equation of motion governed the dynamical behavior of the system and this second-order differential equation is solved directly by the trivialization of the conjugate variables through CT [3]. Park  [3] asserted that because of the non-commutability of operators, CT has not been fully realized in quantum mechanics.  However, CT is usually defined as a change of the noncommutating variables that preserved the commutation relations in quantum mechanics [4].In literature [4, 5] three types of transformations have been identified such as interchange, similarity, and point transformations. Different models of oscillators have been evaluated via the CT [6, 7]. The applications of quantum canonical transformation to a harmonic oscillator in which the angular frequency and equilibrium position are time dependent have been reviewed [3]. In this paper we extend the work of [3] by performing a three-variable transformation on both the space co-ordinate and momentum co-ordinate operators.


## 2. Time-Dependent Harmonic Oscillator

The Time-Dependent Schrodinger Wave equations have been solved explicitly by the method of Quantum Canonical Transformation (QCT). The QCT is also used in solving time-dependent frequency [3]. However, the solutions of time-dependent harmonic oscillator have been obtained through various methods including invariant operator [8], Path Integral [9, 10], and the space-time transformation [11, 12]. We write the Schrödinger equation to be solved as
(1)Ĥ1Ψ1(q,t)=0,
where(2)H1=Pt+P2+Ω2(t)(q-x0(t)),
with Pt=-𝒾∂/∂t, and P=-𝒾∂/∂q,  as the momentum operators in time and in the co-ordinate. It is worthy to note that the angular frequency Ω(t) and the equilibrium position x0(t) are all functions of time.We make time-variable similarity transformation on the momentum time-coordinate and three-variable transformations on space momentum operators as(3)Pt→Pt-𝒾2(∂f∂t)q2,P→P-𝒾fq-𝒾2(∂f∂q).
And the Hamiltonian Operator becomes
(4)H2=Pt-𝒾2(∂f∂t)q2+(P-𝒾fq-𝒾2(∂f∂q))2+Ω2(t)(q-x0(t))2=Pt+P2+(Ω2(t)-𝒾2(∂f∂t)-f2)q2-(𝒾P∂f∂q+fq∂f∂q+14(∂f∂q)2)-2Ω2(t)x0(t)q+Ω2(t)x0(t)-f.
The condition for the differential in (4) to be cancelled is
(5)4∂2f∂q2+4fq(∂f∂q)+(∂f∂q)2=0.
And we solve (5) by series method as(6)f(q)=∑kCkqk+β,
where Ck is the expansion coefficient and K and β are arbitrary constant.Substituting (6) into (5) yields
(7)f(q)=q2+2𝒾2q+8.
Imposing boundary condition on (7), we obtain f(0)=8 and (7) becomes(8)f(q)=q2+2𝒾2q+f(0).
And we recaste (8) in the form
(9)f(q)=q2+2𝒾2q,
where f(0) has been absorbed into f(q). Now multiplying the first term of the LHS of (9) by (-1/2)mω2 and replacing the magnitude 2𝒾2 by f(q) and replacing   f→L-(1/2)mq̇2, (9) becomes(10)L=12mq̇2-12mω2q2+F(q)q,
and (10) is the Lagrangian of the system whose classical equation is
(11)mq̈+mω2q-∂V∂q=0.
The result of (10) shows that the cancellation of the differential in (4) generates the Lagrangian describing the harmonic oscillator of the system. Another condition of (4) is the cancellation of quadratic term in co-ordinate and these yields
(12)𝒾∂f∂t+2f2=2Ω2(t),
however (12) is a first-order nonlinear equation and can be linearized by setting
(13)f=𝒾2∂∂tlnϕ(t),
and (12) becomes(14)∂2∂t2ϕ(t)+4Ω2(t)ϕ(t)=0.
Solving (14) gives
(15)ϕ(t)=ϕ(0)sin(ωt+δ),
where δ is the phase factor. Putting (15) into (13) leads to(16)f(t)=𝒾ω2tan(ωt+δ).
With (16) and (9), we can determine the generatore-f(q)e-f(t)(q2/2).The second Canonical transformation is carried out this time to remove the quadratic term in momentum and the associated differential in co-ordinate:(17)Pt→Pt-𝒾2(∂g∂t)P2,q→q-𝒾gP+𝒾2(∂g∂q).
Substituting (17) into (4) yields
(18)H3=Pt+(2+4gf-𝒾∂g∂t)P22-(𝒾f∂2g∂q2+𝒾Ω2(t)∂g∂q)-2𝒾fgp-2𝒾Ω2x0(t)gP-2Ω2x0(t)P+Ω2x02-f.
We eliminate the quadratic term in momentum by setting(19)2+4gf-𝒾∂g∂t=0,
and solving (19), we obtain
(20)g(t)=-2𝒾e-4𝒾∫f(t)dt∫e4𝒾∫f(t)dt+Ke-4𝒾∫f(t)dtdt.
Similarly, the differential in co-ordinate is eliminated as(21)∂2g∂q2+Ω2(t)f∂g∂q=0,
with f being known, and we determined (21) as(22)g(t)=D+Fe-α2(q)ω2q,
where D and F are arbitrary constants and α(q)=1/f(q).The next transformation is carried out to eliminate the linear term in positionq and the associated differential in co-ordinate:(23)Pt→Pt-𝒾(∂h∂t)q,q→q-𝒾h+𝒾2(∂h∂q),
and the Hamiltonian becomes(24)H4=Pt-(𝒾∂h∂t+2fh+2Ω2x0(t))q-(fq∂h∂q+Ω2x0g∂h∂q)-2𝒾fqp-2𝒾Ω2x0gp-2Ω2x0gh+Ω2x02-f.
The condition for the linear term in q and the differential to be cancelled leads to the solution of h(t) as(25)h(t)=2𝒾e2𝒾∫f(t)dt∫Ω2x0e-2𝒾∫f(t)dtdt,
and the h(q)as(26)h(q)=Ce(∫f(g)[ln(q)]-1/gΩ2)dq,
where C is an arbitrary constant.Another two- andthree-variable transformation is carried out on the Hamiltonian to eliminate the linear term in momentum and the associated differential in co-ordinate:(27)Pt→Pt-𝒾2(∂r∂t)P,q→q-𝒾r+𝒾2(∂r∂q),
and the Hamiltonian becomes(28)H5=Pt-(𝒾2∂r∂t+2𝒾Ω2x0g-2fr)P-2𝒾fqp+𝒾f∂2r∂q2-2Ω2x0gh+Ω2x02-f.
The linear term in P in (28) is eliminated if
(29)𝒾2(∂r∂t)+2𝒾Ω2x0g-2fr=0,
and solving (29) yields(30)r(t)=2e2𝒾∫f(t)dt∫Ω2x0ge-2𝒾∫f(t)dtdt+Ce-2𝒾∫f(t)dt,
and the condition for the differential in co-ordinate to cancelled yields(31)∂2r(q)∂q2=0.
However (31) is a second-order differential Equation and we rewrite (31) in the form(32)𝒾fq2ddq{q2∂r∂q}=0
or(33)∂2r∂q2+1q∂r∂q=0.
Multiplying (33) by q2 gives(34)q2∂2r∂q2+q∂2r∂q=0,
and comparing (34) with (35) we have(35)q2y′′+xy′-λ2y=0.
This shows that (34) is a Bessel equation with λ=0. The general solution of (33) becomes(36)r(q)=A1Jo(mq)+A2Yo(mq),
where m is an integer, A1, A2 are arbitrary constant, and Jo(q) and Yo(q), respectively, represent Bessel functions of first and second kind.Again a point transformation is done to eliminate the co-ordinate in the Hamiltonian of (28) [3]:(37)P→e-qP,q→eq.
Substituting (37) into (28) will be(38)H6=Pt-2𝒾fp-2Ω2x0gh+Ω2x02-f.
In order to trivialize the Hamiltonian of (38), we carried out a two-variable similarity transformation:
(39)Pt→Pt+2𝒾fp+2Ω2x0gh-Ω2x02+f,
and (38) becomes(40)H7=Pt.
The Schrodinger equation associated with the trivialized Hamiltonian of (40) is(41)PtΨ7(q)=H7Ψ7(q),
where Ψ7(q) being the time-dependent Wave function originates from the differential in the canonical transformed Hamiltonian Operators.The complete wave functionΨ1(q,t) is obtained as(42)Ψ1(q,t)=Ψ(q)Ψ(t)=ef(q)eg(p)eh(p)eγ(p)ef(t)q2/2eg(t)p2/2eh(t)Peγ(t)P×Plnge-2P∫f(t)dt+𝒾∫(2Ω2x0gh-Ω2x02+f)dt.
The first four terms of (42) are the times Independent Wave Function arising naturally from the cancellation of the differential in co-ordinate in the Hamiltonian Operator while the other terms are the time-dependent wave function gotten as the consequences of the cancellation in the co-ordinates and momentum.Solving (42) completely yields the result as [3] without the ambiguity of Ψ1(q) which they encountered because Ψ1(q) is obtained directly from the cancellation of the differentials in the Hamiltonian.We obtained after a simple algebra the unnormalized wave function from (11) as(43)Ψ1(q,t)=Ng3n/2e-(1/2)gh2exp(𝒾∫[2ω2x0gh-ω2x02+(2n+1)f]dt)×e(1/2)fq2+hqHn(g2(q-𝒾γg-h)),
where N is the normalization constant and Hn(g(q)/2(q-𝒾γ(q)/Ω-h(q))) is the Hermite polynomial of order n.


## 3. Conclusion

The introduction of a three-variable CT on the co-ordinate and the momentum space removed the ambiguity encountered in [3] in evaluating the time-independent Schrödinger wave function. In (41), we obtain the generalized wave function of the time dependent harmonic oscillator with the first four terms arising from the introduction of the three canonical variables. These terms were simply given as Ψ1 in [3] without its explicit values as we have obtained above. Rather fortunately the three-variable transformations separate naturally the time-dependent and time-independent function explicitly. This separation aids us to solve the Schrödinger equation completely by trivializing the time-dependent Hamiltonian Operator. In addition our result can be viewed as a generalized and explicit case of that obtained in [3].The extension of this work to different time dependent systems with different potential is possible. Thus canonical transformation is a good tool in solving time-dependent systems.



---

*Source: 103538-2009-12-02.xml*