%Numero de paginas en dvi = 6 \documentclass{rmf-d} \usepackage{nopageno,rmfbib,multicol,times,epsf,amsmath,amssymb,cite} \usepackage[latin1]{inputenc} \usepackage[]{caption2} \usepackage{graphics} % \def\rmfcornisa{RESEARCH \hfill\rmf\ {\bf 64} (2018) 30--35 \hfill JANUARY-FEBRUARY 2018} \newcommand{\ssc}{\scriptscriptstyle} % %\numero{1} % \newcommand{\nt}{\notag \\ & \times} \newcommand{\na}{\notag \\ & } \newcommand{\no}{\notag \\} % \def\rmfcintilla{{\it Rev.\ Mex.\ Fis.\/} {\bf 64} (2018) 30--35} \clearpage \rmfcaptionstyle \pagestyle{myheadings} \setcounter{page}{30} \markboth{Surarit Pepore } {Integrals of the motion and green functions for time-dependent mass harmonic oscillators} \begin{document} \title{Integrals of the motion and Green functions for time-dependent\\ mass harmonic oscillators \vspace{-6pt}} \author{Surarit Pepore } \address{Department of Physics, Faculty of Science and Technology, \\ Rajamangala University of Technology Thanyaburi, \\ Rangsit-Nakornayok Road, Pathumthani 12110, Thailand. \\ e-mail: surapepore@gmail.com } \maketitle \recibido{31 July 2017}{25 October 2017 \vspace{-12pt}} \begin{abstract} The application of the integrals of the motion of a quantum system in deriving Green function or propagator is established. The Green function is shown to be the eigenfunction of the integrals of the motion which described initial points of the system trajectory in the phase space. The explicit expressions for the Green functions of the damped harmonic oscillator, the harmonic oscillator with strongly pulsating mass, and the harmonic oscillator with mass growing with time are obtained in co-ordinate representations. The connection between the integrals of the motion method and other method such as Feynman path integral and Schwinger method are also discussed. \end{abstract} \keys{ Integrals of the motion;\ Green function;\ Time-dependent mass harmonic oscillators \vspace{-4pt}} \pacs{03.65.-w \vspace{-4pt}} \begin{multicols}{2} \section{Introduction} %1 In non-relativistic quantum mechanics, the propagator is represented as the transition probability amplitude for a particle to motion from initial space-time configuration to final space-time configuration. The Feynman path integral [1] and the Schwinger action principle [2] are the well-known methods in calculating the propagator. The aim of this paper is to present the connection between the integrals of the motion of a quantum system and its Green function or propagator. As reveal by V.V. Dodonov \textit{et al}. [3] that the Green function is the eigenfunction of the integrals of the motion describing initial points of the system trajectory in the phase space. D.B. Lemeshevskiy and V.I. Man'ko [4] constructed the Green functions for the driven harmonic oscillator with the aid of integrals of the motion. In the present paper we want to calculate the Green functions or propagators for the damped harmonic oscillator [5-7], the harmonic oscillator with strongly pulsating mass, [8] and the harmonic oscillator with mass growing with time [9] by the method developed by V.V. Dodonov \textit{et al}. [3] This paper is organized as follows. In Sec. 2, the Green function for the damped harmonic oscillator is derived. In Section 3, the calculation of the Green function for the harmonic oscillator with strongly pulsating mass is presented. The Green function for the harmonic oscillator with mass growing with time is evaluated in Sec. 4. Finally, the conclusion is given in Sec. 5. \section{The Green function for a damped harmonic oscillator} %2 The Hamiltonian operator for a damped harmonic oscillator is described by [5-7] \begin {align} \hat{H} (t) = e ^{-rt}\frac{\hat{p}^2}{2m} + \frac{1}{2} m \omega^2 e^{rt} \hat{q}^2, \end {align} where $r$ is the damping constant coefficient. The aim of this section is to drive the Green function $G(x,x',t)$ of the Schrodinger equation by the method of integrals of motion [3-4]. The classical correspondence of the Hamiltonian operator in Eq. (1) is \begin {align} H (q,p,t) = e ^{-rt}\frac{{p}^2}{2m} + \frac{1}{2} m \omega^2 e^{rt} {q}^2, \end {align} The Hamilton equation of motion for position and momentum are [10] \begin {align} \dot{q} = \frac{p}{m} e^{-rt}, \qquad \dot{p}= - -m\omega^2 e^{rt} q . \end {align} The classical paths in the phase space under the initial conditions $q(0)=q_0$ and $p(0)=p_0$ are given by \begin {align} q(t) &= q_0 \left(e^{-rt}\cos \Omega t + \frac{re^{-rt/2}}{2 \Omega}\right) \na + p_0 \left(\frac{re^{-rt/2}}{m \Omega}\sin\Omega t\right), \\ p(t) &= p_0 \left(e^{rt/2}\cos \Omega t - \frac{re^{rt/2}}{2 \Omega}\sin \Omega t \right) \na - q_0 \left(\frac{m \omega^2 }{ \Omega}re^{-rt/2} \sin\Omega t\right), \end {align} where $\Omega^2 = \omega^2 - r^2 / 4$. Now we consider the system of Eqs.~(4) and (5) as an algebraic system for unknown initial position\hspace{1pt} $q_0$\hspace{1pt} and\hspace{1pt} momentum\hspace{1pt} $p_0$,\hspace{1pt} respectively.\hspace{1pt} The\hspace{1pt} variables \pagebreak \noindent $q$, $p$, and $t$ are taken as the parameters. The solution of this system are given as \begin {align} q_0 (q,p,t)& = q \left( e^{rt/2}\cos \Omega t - \frac{r}{2\Omega} e^{rt/2} \sin \Omega t \right) \na -p \left( \frac{e^{- {rt}/{2}}}{m\Omega} \sin \Omega t \right), \\ p_0 (q,p,t)& = q \left( \frac{m \omega^2 }{\Omega } e^{rt/2}\sin \Omega t \right) \na + p \left(e^{-{rt}/{2}}\cos \Omega t+ \frac{r e^{-{rt}/{2}} }{2 \Omega} \sin \Omega t \right). \end {align} We define operators acting in the Hilbert space as follows \begin {align} \hat{q}_0 (\hat{q},\hat{p},t)& = \hat{q} \left( e^{rt/2}\cos \Omega t - \frac{r}{2\Omega} e^{rt/2} \sin \Omega t \right) \na - \hat{p} \left( \frac{e^{-{rt}/{2}}}{m\Omega} \sin \Omega t \right), \\ \hat{p}_0 (\hat{q},\hat{p},t)& = \hat{q} \left( \frac{m \omega^2 }{\Omega } e^{rt/2}\sin \Omega t \right) \na + \hat{p} \left(e^{-{rt}/{2}}\cos \Omega t+ \frac{r e^{-{rt}/{2}} }{2 \Omega} \sin \Omega t \right). \end {align} Calculating the total derivative of the operator $\hat{q}_0 (\hat{q},\hat{p},t)$ with respect to time $t$, we obtain \begin {align} \frac{d \hat{q}_0}{dt} = \frac{\partial \hat{q}_0}{\partial t } + \frac{i}{ \hbar } [\hat{H}, \hat{q}_0]. \end {align} Similarly, the total time-derivative of the operator $\hat{p}_0 (\hat{q},\hat{p},t)$ is \begin {align} \frac{d \hat{p}_0}{dt} = \frac{\partial \hat{p}_0}{\partial t } + \frac{i}{ \hbar } [\hat{H}, \hat{p}_0]. \end {align} Thus, operators in Eqs. (8) and (9) are integrals of the motion and correspond to the initial position and momentum. Then these operators must satisfy equations for the Green function $G(x,x',t)$, [3-4] \begin {align} \hat{q}_0 (x) G(x,x',t) = \hat{q} (x')G(x,x',t), \\ \hat{p}_0 (x) G(x,x',t) = -\hat{p} (x')G(x,x',t), \end {align} where the operators on the left-hand sides of the equations act on variable $x$, and on the right- hand sides, on $x'$. Now we write Eqs. (12) and (13) explicitly, \begin {align} &\bigg( x \left( e^{rt/2} \cos \Omega t - \frac{r}{2 \Omega} e^{rt/2} \sin \Omega t \right) \na + \frac{i\hbar}{m \Omega} e^{- rt/2} \sin \Omega t \frac{\partial}{\partial x} \bigg)G(x,x',t) = x' G(x,x',t), \\ &\bigg( x \left(\frac{m\omega^2}{\Omega} e^{rt/2} \sin \Omega t\right) - i \hbar \bigg( e^{- rt/2}\cos \Omega t \na {+}\frac{re^{- rt/2}}{2 \Omega}\sin \Omega t \bigg) %\nt \frac{\partial}{\partial x}\bigg) G(x,x',t) {=} i\hbar \frac{\partial G(x,x',t)}{\partial x'}. \end {align} By modifying Eqs. (14) and (15), the system of equations for deriving the Green function $G(x,x',t)$ are \begin {align} & \frac{\partial G(x,x',t) }{\partial x } = - \frac{im \Omega }{ \hbar} \bigg( \frac{e^{rt/2}}{\sin \Omega t } x' \na - \big( e^{rt} \cot \Omega t - \frac{re^{rt}}{ 2 \Omega} \big) x \bigg) G(x,x',t), \\ & \frac{\partial G(x,x',t) }{\partial x' } = - \frac{i }{ \hbar} \bigg( \frac{m \Omega e^{rt/2}}{\sin \Omega t } x \na - \Big( m\Omega \cot \Omega t + \frac{mr}{2} \Big) x' \bigg) G(x,x',t). \end {align} Now one can integrate Eq. (16) with respect to the variable $x$ to obtain \begin {align} G(x,x',t) &= C(x',t)\exp \bigg( \frac{i}{ \hbar } \bigg \{ \bigg ( \frac{m \Omega}{2} e^{rt} \cot \Omega t \na - \frac{mr}{4} e^{rt} \bigg) x^2-\frac{m \Omega}{\sin \Omega t} e^{rt/2} xx'\bigg\} \bigg), \end {align} where $C(x',t)$ is the function of $x'$ and $t$. Substituting Eq. (18) into Eq. (17), we obtain the differential equation for $C(x',t)$ as \begin {align} \frac{\partial C(x',t)}{\partial x'} = \frac{i}{ \hbar} \big( m\Omega \cot \Omega t + \frac{mr}{2} \big) x' C(x', t). \end {align} Solving Eq. (19), the function $C(x',t)$ can be expressed as \begin {align} C(x',t) = C(t) \exp \bigg( \frac{i}{ \hbar} \bigg(\frac{m\Omega}{2}\cot \Omega t + \frac{mr}{4} \bigg) x'^2 \bigg), \end {align} where $C(t)$ is the pure function of time. So, the Green function in Eq. (18) can be written as \begin {align} G(x,& x',t)= C(t) \exp \bigg( \frac{i}{\hbar} \bigg\{ \bigg( \frac{m\Omega}{2} e^{rt}\cot \Omega t- \frac{mr}{4} e^{rt} \bigg)x^2 \na + \bigg( \frac{m\Omega}{2} \cot \Omega t + \frac{mr}{4} \bigg) x'^2 - \frac{m\Omega e^{rt/2}}{\sin \Omega t} xx'\bigg\} \bigg). \end {align} To find $C(t)$, we must substitute the Green function of Eq. (21) into the Schrodinger equation \begin {align} & i\hbar \frac{\partial G(x,x', t) }{\partial t} = - \frac{\hbar^2}{2m} e^{-rt} \frac{\partial^2 G(x,x', t)}{\partial x^2 } \na +\frac{1}{2}m \omega^ 2 e^{rt}x^2 G(x,x', t). \end {align} After some algebra, we obtain an equation that does not contain the variables $x$ and $x'$, \begin {align} \frac{dC(t)}{dt} = C(t) \left( \frac {r} {2} - \frac {\Omega \cot \Omega t} {2} \right) \end {align} Eq. (23) can be simply integrated with respect to time $t$, and one obtains \begin {align} C(t) = \frac {c} {\sqrt{\sin \Omega t}} e^{rt/4} , \end {align} where $C$ is a constant. Substituting Eq. (24) into Eq. (21) and applying the initial condition \begin {align} G(x,x',t=0)=\delta (x-x'), \end {align} we get \begin {align} C= \sqrt{\frac{m \Omega }{2 \pi i \hbar}} \, . \end {align} So, the Green function or propagator for a damped harmonic oscillator can be written as \begin {align} G(x,& x',t ) = \sqrt{\frac{m \Omega e^{rt/2}}{2\pi i \hbar \sin \Omega t}} \exp \bigg( \frac {i} {\hbar} \bigg\{ \bigg( \frac{m \Omega}{2} e^{rt}\cot \Omega t \na - \frac{mr e^{rt}}{4} \bigg)x^2 + \bigg( \frac{m \Omega}{2} \cot \Omega t %\na +\frac{mr }{4} \bigg)x'^2 \na - \frac{m \Omega e^{rt/2}}{ \sin \Omega t} xx' \bigg \} \bigg), \end {align} which is the same form as the result of S. Pepore \textit{et al}. [5] calculating from Feynman path integral. \section{The Green function for a harmonic oscillator with strongly pulsating mass} %3 The Hamiltonian operator for a harmonic oscillator with strongly pulsating mass can be expressed as [9] \begin {align} \hat{H} (t) = \frac{p^2}{2m \cos^2 vt} + \frac{1}{2} m \cos^2 vt \omega^2 \hat{q}^2 , \end {align} where $ v $ is the frequency of mass. The classical analog of the Hamiltonian operator in Eq. (28) is \begin {align} H(q,p,t) = \frac{p^2}{2m \cos^2 vt} +\frac{1}{2} m\cos^2 vt \omega^2 q^2 . \end {align} The classical equations of motion determining the oscillator position and momentum are \begin {align} \ddot{q} - 2v \tan v t\dot{q} + \omega^2 q =0 \end {align} The classical trajectories in the phase space under the initial conditions $q(0)=q_0$ and $p(0)= p_0$ can be written as \begin {align} q(t) & =q_0 \sec vt \cos \Omega t+ p_0/m \Omega \sec v t \sin \Omega t, \\ p(t) & =q_0 (mv \cos vt \tan vt \cos \Omega t -m \Omega \cos vt \sin \Omega t) \na +p_0 (\cos vt \cos \Omega t+ \frac{ v}{\Omega } \cos vt \tan vt sin \Omega t), \end {align} where $ \Omega^2 = \omega^2 + v^2 $. By eliminating $p_0$ in Eq. (31) and $q_0$ in Eq. (32), the solutions are \begin {align} q_0 (q,p,t)& = q \bigg(\cos vt \cos \Omega t+ \frac{v}{\Omega} \sin vt \sin \Omega t\bigg) \na -p\bigg( \frac{ \sec vt \sin \Omega t}{m \Omega} \bigg), \\ p_0 (q,p,t) & =q (m \Omega \cos vt \sin \Omega t \na -m v \sin vt \cos \Omega t)+p(\sec vt \cos \Omega t). \end {align} The Hilbert space operators of $q_0$ and $p_0$ are \begin {align} \hat{q}_0 (\hat{q}, \hat{p},t)& = \hat{q} \bigg(\cos vt \cos \Omega t+ \frac{v}{\Omega} \sin vt \sin \Omega t\bigg) \na - \hat{p} \left(\frac{\sec vt \sin \Omega t }{ m \Omega} \right), \\ \hat{p} (\hat{q} ,\hat{p},t)& =\hat{q}(m\Omega \cos vt \sin \Omega t-mv \sin vt \cos \Omega t) \na + \hat{p}(\sec vt \cos \Omega t). \end {align} We can determine that $\hat{q}_0$ and $\hat{p}_0$ are integrals of the motion by finding total time derivatives of \begin {align} \frac{d\hat{q}_0 }{dt }&=\frac{\partial \hat{q}_0 }{\partial t }+ \frac{i }{ \hbar} [\hat{H}, \hat{q}_ 0]=0, \\ \frac{d\hat{p}_0 }{dt }&=\frac{\partial \hat{p}_0 }{\partial t }+ \frac{i }{ \hbar} [\hat{H}, \hat{p}_ 0]=0. \end {align} Then these operators must satisfy the equations for the Green function $G(x,x',t)$ [3-4] \begin {align} \hat{q}_0 (x)G(x,x',t) & =\hat{q}(x' )G(x,x',t), \\ \hat{p}_0 (x)G(x,x',t) & =-\hat{p}(x' )G(x,x',t). \end {align} Now we write Eqs. (39) and (40) explicitly, \begin {align} & \bigg(x \bigg( \cos vt \cos \Omega t + \frac{v}{\Omega}\sin vt \sin \Omega t \bigg) \na + \frac{i \hbar \sec vt \sin \Omega t }{m \Omega } \frac{ \partial }{ \partial x}\bigg) G(x,x',t) = x' G (x,x',t) , \\& \bigg( x \big(m\Omega \cos vt \sin \Omega t - mv \sin vt \cos \Omega t \big) \na - i\hbar \sec vt \cos \Omega t \frac{\partial}{\partial x} \bigg) G(x,x',t) \na = i \hbar \frac{ \partial}{\partial x'} G (x,x',t) \end {align} The system of equations for defining the Green function $G(x,x',t)$ are \vspace{-0pt} \begin {align} \frac{\partial G(x,x',t) }{\partial x } &= - \frac{i}{\hbar} [ x' (m \Omega \cos vt \csc \Omega t ) \na -x(m\Omega \cos^2 vt \cot \Omega t \na + mv \sin vt \cos vt )] G(x,x',t), \\ \frac{\partial G(x,x',t) }{\partial x' } &= - \frac{i}{\hbar} [ x (m \Omega \cos vt \csc \Omega t ) \na -x'(m\Omega \cot \Omega t )] G(x,x',t). \end {align} Now we can integrate Eq. (43) with respect to the variable $x$ to get \begin {align} G(x,x',t) &= C(x', t ) \exp \bigg [ \frac{i}{ \hbar} \bigg(\bigg( \frac{1}{2} m \Omega \cos^2 vt \cot \Omega t \na +\frac{1}{2}mv \sin vt \cos vt \bigg)x^2 \na - m\Omega \cos vt \csc \Omega t xx' \bigg) \bigg]. \end {align} Substituting Eq. (45) into Eq. (44), we obtain the differential equation for $C(x',t)$ as \begin {align} \frac{\partial C(x',t) }{\partial x' } = \frac{i}{\hbar} (m\Omega \cot \Omega t) x' C(x',t) \end {align} Solving Eq.(46), we obtain \begin {align} C(x',t)=C(t) \exp \left(\frac{i}{2\hbar} m\Omega \cot \Omega t x'^2 \right). \end {align} After substituting Eq. (47) into Eq. (45), we arrive at \begin {align} G(x,x',t) &= C( t ) \exp \bigg [ \frac{i}{ \hbar} \bigg(\bigg( \frac{1}{2} m \Omega \cos^2 vt \cot \Omega t \na +\frac{1}{2}mv \sin vt \cos vt \bigg)x^2 \na \frac{1}{2} m\Omega \cot \Omega t x'^2 - m\Omega \cos vt \csc \Omega t xx' \bigg) \bigg]. \end {align} To find $C(t)$, we must substitute the Green function of Eq. (48) into the Schr\"odinger equation \begin {align} i\hbar \frac{ \partial }{ \partial t} G(x,x',t) &= - \frac{\hbar^2 }{2 m \cos^2 vt} \frac{\partial^2 G(x,x',t) }{\partial x^2 } \na + \frac{1}{2} m \cos^2 vt \omega^2 x^2 G(x,x',t). \end {align} After some algebra, we get \begin {align} \frac{ dC(t)}{dt} = - \frac{1}{2} (\Omega \cot \Omega t + v \tan vt )C(t). \end {align} Integrating Eq. (50) with respect to time, we obtain \begin {align} C(t) = C \left( \frac{\cos vt }{\sin \Omega t } \right)^{1/2} . \end {align} Substituting Eq. (51) into Eq. (48) and applying the initial condition in Eq. (25), the constant $C$ is \begin {align} C= \sqrt{\frac{m \Omega}{2 \pi i \hbar}} \end {align} Thus, the Green function for a harmonic oscillator with strongly pulsating mass can be expressed as \begin {align} G(x,x',t) &= \sqrt{\frac{ m \Omega \cos vt}{ 2 \pi i \hbar \sin \Omega t}} \exp \bigg[\frac{i}{ \hbar } \bigg( \bigg( \frac{1}{2} m \Omega \cos^2 vt \cot \Omega t \na + \frac{1}{2} mv \sin vt \cos vt \bigg)x^2 + \frac{1}{2} m\Omega \cot \Omega tx'^2 \na - m\Omega \cos vt \csc \Omega t xx' \bigg) \bigg], \end {align} which is the same result as M. Sabir and S. Rajagopalan [9] by Feynman path integral method. \section{The Green function for a harmonic oscillator with mass growing with time} %4 The Hamiltonian operator for a harmonic oscillator with mass growing with time can be written as \begin {align} \hat{H} (t) = \frac{ \hat{p}^2 }{2m (1+\alpha t)^2} + \frac{1}{2} m (1 + \alpha t)^2 \omega^2 \hat{q}^2, \end {align} where $ \alpha $ is a constant. \begin {align} H(q,p,t) = \frac{ {p}^2 }{2m (1+\alpha t)^2} + \frac{1}{2} m (1 + \alpha t)^2 \omega^2 {q}^2 \end {align} The equation of motion for this oscillator is \begin {align} \ddot{q} + \frac{2\alpha}{(1+ \alpha t)}\dot{q} + \omega^2 q = 0 . \end {align} The classical paths in the space under the initial conditions $q(0)=q_0$ and $p(0)=p_0$ can be expressed as \begin {align} q(t) & = q_0 \left( \frac {\alpha \sin \omega t + \omega \cos \omega t} {\omega (1+ \alpha t)} \right) \na + p_0 \left( \frac {\sin \omega t} {m \omega (1+ \alpha t)} \right), \\ p(t) & {=} q_0 \bigg(m \alpha^2 t \cos \omega t{-} m\omega (1- \alpha t) \sin \omega t {-} \frac{m \alpha^2}{ \omega } \sin \omega t \bigg) \na {+} p_0 \bigg((1{+}\alpha t )\cos \omega t {-} \frac{ \alpha}{\omega} \sin \omega t \bigg). \end {align} We can express $q_0$ and $p_0$ in terms of $q$, $p$, and $t$ by \begin {align} q_0 (q,p,t) &= q \bigg( (1+ \alpha t ) \cos \omega t - \frac {\alpha} {\omega} \sin \omega t \bigg ) \na - p \left( \frac {\sin \omega t} {m \omega (1+\alpha t)} \right), \\ p_0 (q,p,t) &{=} p \left( \frac {\alpha \sin \omega t {+} \omega \cos \omega t} {\omega(1{+}\alpha t)} \right) \na {-} q \bigg( m \alpha^2 t \cos \omega t {-} m \omega (1+\alpha t) \sin \omega t \na - \frac{m \alpha^2 \sin \omega t}{ \omega } \bigg) . \end {align} We define operators acting in the Hilbert space as follows \begin {align} \hat{q}_0 (\hat{q},\hat{p}, t)& = \hat{q} \bigg( (1-\alpha t) \cos \omega t - \frac{\alpha}{\omega} \sin \omega t \bigg) \na - \hat{p} \bigg( \frac{\sin \omega t }{m \omega (1+ \alpha t)} \bigg), \\ \hat{q}_0 (\hat{q},\hat{p}, t) &= \hat{q} \left( \frac {\alpha \sin \omega t + \omega \cos \omega t} {\omega (1+\alpha t)} \right) \na \hat{q} \bigg( m\alpha^2 t \cos \omega t - m \omega(1+\alpha t) \sin \omega t \na - \frac{m \alpha^2 \sin \omega t }{ \omega} \bigg) . \end {align} Calculating the total derivatives of the operators $\hat{q}_0$ and $\hat{q}_0$ with respect to time, we obtain \begin {align} \frac{d\hat{q}_0 }{dt } & = \frac{\partial \hat{q}_0 }{\partial t } + \frac{i}{ \hbar} [\hat{H}, \hat{q}_0 ] =0, \\ \frac{d\hat{p}_0 }{dt } & = \frac{\partial \hat{p}_0 }{\partial t } + \frac{i}{ \hbar} [\hat{H}, \hat{p}_0 ] =0. \end {align} Hence, operators in Eqs. (61) and (62) are integrals of motion and correspond to the initial position and momentum. Then these operators must satisfy the equations for the Green function $G(x,x',t)$ [3-4] \begin {align} \hat{q}_0 (x)G(x,x',t) & = \hat{q}(x' )G(x,x',t), \\ \hat{p}_0 (x)G(x,x',t) & = -\hat{p}(x' )G(x,x',t). \end {align} Writing Eqs. (65) and (66) explicitly, it can be shown that \begin {align} \bigg( & x \bigg[ (1+\alpha t) \cos \omega t - \frac{\alpha}{\omega} \sin \omega t \bigg] + \frac{i\hbar \sin \omega t}{m \omega (1+ \alpha t) } \frac{\partial }{\partial x }\bigg) \nt G(x,x',t) = x' G(x,x',t), \\ \bigg( & -i \hbar \bigg(\frac{\alpha \sin \omega t + \omega \cos \omega t }{\omega (1+\alpha t)} \bigg) \frac{ \partial}{ \partial x} - x \bigg(m\alpha^2 t \cos \omega t \na - m\omega (1+ \alpha t) \sin \omega t- \frac{m \alpha^2 \sin \omega t }{\omega} \bigg) G (x,x',t) \na = i\hbar \frac{\partial G (x,x',t)}{\partial x' } \end {align} The system of equations for calculating the Green function $G(x,x',t)$ are \begin {align} \frac{\partial G(x,x',t) }{\partial x } &=\bigg(-\frac{im\omega(1+\alpha t)x' }{\sin \omega t } +\frac{ix }{\hbar} [m\omega (1+ \alpha t)^2 \nt \cot \omega t -m\alpha (1+\alpha t)]\bigg)G(x,x',t), \\ \frac{\partial G(x,x',t) }{\partial x' } &=\bigg(\frac{i m }{\hbar } (\alpha + \omega \cot \omega t) x' - \frac{ im \omega(1+\alpha t)}{ \hbar \sin \omega t} x \bigg) \nt G(x,x',t). \end {align} Now we can integrate Eq. (70) with respect to the variable $x$ to obtain \begin {align} G(x,x',t) & = C(x', t)\exp \bigg[ \frac{i}{2 \hbar} \bigg( (m \omega (1+\alpha t)^2 \cot \omega t \na {-} m\alpha (1+\alpha t) ) x^2 {-} \frac{2 m \omega (1+\alpha t) }{\sin \omega t} xx' \bigg) \bigg]. \end {align} Substituting Eq. (71) into Eq. (70), we obtain the differential equation for $C(x',t)$ as \begin {align} \frac{\partial C(x', t) }{ \partial x'}= \frac{im }{\hbar}(\alpha + \omega \cot \omega t )x' C(x', t). \end {align} Solving Eq. (72), we obtain \begin {align} C(x', t)= C(t) \exp \bigg( \frac{i}{2\hbar} (m \alpha + m\omega \cot \omega t)x'^2 \bigg) \end {align} After substituting Eq. (73) into Eq. (71), we obtain \begin {align} G(x,x',t) & = C(x', t)\exp \bigg( \frac{i}{2 \hbar} \bigg[ (m \omega (1+\alpha t)^2 \cot \omega t \na {-} m\alpha (1+\alpha t) ) x^2 {+} (m\alpha + m\omega \cot \omega t)x'^2 \na - \frac{2 m \omega (1+\alpha t) }{\sin \omega t} xx' \bigg] \bigg). \end {align} To get $C(t)$, we must substitute the Green function of Eq. (74) into the Schr\"odinger equation \begin {align} i\hbar \frac{\partial G(x,x',t) }{ \partial t} &= - \frac{\hbar^2 }{2m(1+\alpha t)^2 } \frac{\partial^2 G(x,x',t) }{\partial x^2 } \na + \frac{1}{2} m (1+\alpha t)^2 \omega^2 G(x,x',t). \end {align} After some algebra, we obtain \begin {align} \frac{dC(t) }{dt} = \bigg( \frac{\alpha }{2(1+\alpha t) } - \frac{1}{2}\omega \cot \omega t \bigg)C(t). \end {align} Integrating Eq. (76) with respect to time, we get \begin {align} C(t) = C \bigg(\frac{1+ \alpha t}{\sin \omega t} \bigg)^{1/2}. \end {align} Substituting Eq. (77) into Eq. (74) and employing the initial condition in Eq. (25), the constant $C$ becomes \begin {align} C= \sqrt{\frac{m \omega }{2 \pi i \hbar}} . \end {align} So, the Green function for a harmonic oscillator with mass growing with time can be written as \begin {align} G(x,x',t) & = \bigg(\frac{m \omega (1+\alpha t) }{2\pi i \hbar \sin \omega t } \bigg)^{1/2} \exp \bigg( \frac{i}{2 \hbar} \bigg[ (m \omega (1+\alpha t)^2 \nt \cot \omega t {-} m\alpha (1+\alpha t) ) x^2 {+} (m\alpha + m\omega \cot \omega t)x'^2 \na - \frac{2 m \omega (1+\alpha t) }{\sin \omega t} xx' \bigg] \bigg), \end {align} which is agreement with the result of S. Pepore and B. Sukbot~[11] calculating by Schwinger method. \section{Conclusion} %5 The method for deriving the Green functions with the helping of integrals of the motion presented in this paper can be successfully applied in solving time-dependent mass harmonic oscillator problems. This method has the important steps in finding the constant of motions $q_0$ and $p_0$ and implying that the Green functions $G(x,x',t)$ is the eigenfunctions of the operators $\hat{q}_0 (x)$ and $\hat{p}_0 (x)$. In fact, this method has many common features with the Schwinger method, [11-14] but the Schwinger method requires the operator $\hat{q}_0 (x)$ and $\hat{p}_0 (x)$ in calculating the matrix element of Hamiltonian operator in the Green function \begin {align} G(x,x',t) & = C(x,x') \exp \bigg\{ - \frac{i}{\hbar} \nt \int\limits_{0}^{t} \frac{ \langle x(t) \mid \hat{H} ( \hat{x}(t),\hat{x} (0) \mid x' (0) \rangle \mid \rangle }{\langle x(t) \mid x'(0) \rangle} dt \bigg\}\, . \end {align} In Feynman path integral, the pre-exponential function $C(t)$ comes from sum over all historical paths that depend on the calculation of functional integration while in the integrals of motion method this term appears in solving Schrodinger equation of Green function. In my opinion the method in this article seems to be more simple from the viewpoint of calculation. \end{multicols} \medline \begin{multicols}{2} \begin{thebibliography}{99} %\bibitem{c1} %1. \bibitem{c1} R.P. Feynman and A.R Hibbs, \textit{Quantum Mechanics and Path Integral}, (McGraw-Hill, New York, 1965). %2. \bibitem{c1} J. Schwinger, \textit{Phys. Rev.} \textbf{82} (1951) 664. %3. \bibitem{c1} V.V. Dodonov, I.A. Malkin, and V. I. Man'ko, \textit{Int. J. Theor. Phys.} \textbf{14} (1975) 37. %4. \bibitem{c1} D.B. Lemeshevskiy and V.I. Man'ko, \textit{Journal of Russian Laser Research}, \textbf{33} (2012) 166. %5. \bibitem{c1} S. Pepore, P. Winotai, T. Osotchan, and U. Robkob, \textit{Science Asia.} \textbf{32} (2006) 173. %6. \bibitem{c1} P. Caldirola, \textit{Nuovo Cim.} \textbf{18} (1941) 393. %7. \bibitem{c1} E. Kanai, \textit{Prog. Theor. Phys.} \textbf{3} (1948) 440. %8. \bibitem{c1} R.K. Colegrave and M.S. Abdalla, \textit{J. Phys. A.} \textbf{19} (1982) 1549. %9. \bibitem{c1} M. Sabir and S. Rajagopalan, \textit{J. Phys. A.} \textbf{37} (1991) 253. %10. \bibitem{c1} H. Goldstein, C. Poole, and J. Safko, \textit{Classical Mechanics}, (Addison Wesley, San Francisco, 2000). %11. \bibitem{c1} S. Pepore and B. Sukbot, \textit{Chinese. J. Phys.} \textbf{53} (2015) 060004. %12. \bibitem{c1} S. Pepore and B. Sukbot, \textit{Chinese. J. Phys.} \textbf{47} (2009) 753. %13. \bibitem{c1} S. Pepore and B. Sukbot, \textit{Chinese. J. Phys.} \textbf{53} (2015) 100002. %14. \bibitem{c1} S. Pepore and B. Sukbot, \textit{Chinese. J. Phys.} \textbf{53} (2015) 120004. \end{thebibliography} \end{multicols} \end{document} % %Notas: % %\tabletopline\vspace{2pt}\lilahf{\sc Table I.\ {\rm Table caption}} %\begin{center} %\small{\renewcommand{\arraystretch}{1.3} %\renewcommand{\tabcolsep}{1.35pc} %\begin{tabular}{cccc} %\hline %\hline %\end{tabular}} %\end{center} %