Quantization of a 3 D Nonstationary Harmonic plus an Inverse Harmonic Potential System

1Department of Material Sciences, Faculty of Science, University of M’sila, 28000 M’sila, Algeria 2Faculty of Technology, University of Ferhat Abbas Sétif-1, 19000 Sétif, Algeria 3Laboratory of Optoelectronics and Compounds (LOC), Department of Physics, Faculty of Science, University of Ferhat Abbas Sétif-1, 19000 Sétif, Algeria 4Department of Radiologic Technology, Daegu Health College, Daegu 41453, Republic of Korea


Introduction
The quantum dynamics of central potential systems that is originally developed in the context of stationary Hamiltonian systems [1][2][3] can be extended to more complicated systems that involve time-dependent parameters.Time-dependent Hamiltonians are useful for describing nonstationary quantum dynamical systems which are ubiquitous in the physics world.Among many time-dependent Hamiltonian systems, central potential systems with time-dependent parameters have attracted considerable interest in the literature because they provide solvable models for various physical systems.As a consequence, active researches have been carried out in this context, accompanying lots of reports so far [4][5][6][7][8][9][10][11][12][13][14][15][16].For instance, a dissipative harmonic oscillator perturbed by an inverse quadratic potential described with time-dependent parameters has been studied [7].Quantum properties of harmonic motion of a charged particle that have timedependent effective mass under time-dependent singular potential have been investigated [16].However, the topics of most of these researches were limited to the case of one-or two-dimensional systems due to the difficulty in developing mathematical procedures for three-or higher-dimensional systems that involve time-dependent parameters.
For a full understanding of profound physical properties of such systems, quantum analytical solutions of corresponding Schrödinger equations, as well as classical solutions, are necessary.Various mathematical tools and methods for solving Schrödinger equations of complicated time-dependent Hamiltonian systems have been suggested and developed up to now.They are the invariant method [17], the method based on path integral formulation [18], the trial function method [19], the canonical/unitary transformation approach [20], and so on.In particular, the invariant method which was devised by Lewis and Riesenfeld [21,22] is very useful for unfolding fundamental theories of quantum dynamics for the case where the Hamiltonian of the given system explicitly depends on time.Among diverse nonstationary quantum systems, the time-dependent quantum system with a harmonic plus an inverse harmonic potential in one or two dimensions is one of the interesting problems that has attracted permanent interest on account of its emerging applications in actual physical systems [9,16,23,24].For instance, exact quantum solutions for many body systems [2,3] and polyatomic molecules [23] and the quantum behavior of a charged particle in electromagnetic fields [24] can be investigated through this model.These studies for nonstationary quantum systems may be possibly extended to more complicated systems, such as the time-dependent quantum systems with spherical harmonic plus an inverse harmonic potential in three dimensions.
Stimulated by this trend for recent researches, we investigate exact quantum solutions of the time-dependent harmonic potential system perturbed by an inverse harmonic potential in three dimensions in the present work.To this end, several particular mathematical procedures will be employed in order to overcome the mathematical difficulty for unfolding quantum theory under the consideration of time-dependent parameters in three dimensions.Typically, the Schrödinger solutions or wave functions for a timedependent Hamiltonian system are represented in terms of the eigenfunctions of an invariant operator of the system.Such eigenfunctions are obtained by solving the eigenvalue equation of the invariant operator.Because the invariant operator is represented in terms of some time functions, the direct evaluation of the eigenvalue equation is somewhat difficult.This difficulty can be resolved by introducing the unitary transformation method and the Nikiforov-Uvarov (NU) method together.The eigenvalue equation can be simplified by unitary transformation with an appropriate unitary operator.The NU method enables us to solve such arbitrary-type second-order differential equations.

Hamiltonian and Invariant Operator
As is well known, the study of quantum dynamical properties of a system requires Schrödinger solutions as a preliminary step.While quantum features of time-dependent harmonic plus an inverse harmonic potential in one or two dimensions are extensively studied in the literature, the same problems in three dimensions are not extensively studied yet as far as we know.As a solvable dynamical system in three dimensions, we consider a time-dependent harmonic plus an inverse harmonic potential with a linearly increasing mass [25,26].In this case, the Hamiltonian has the form where  0 ,  0 , , , and  are real constants.Notice that the linearly increasing mass appears in several dynamical systems such as the expanding shells with mass accretion in a plasma [25] and the mass-accreting Algol-type eclipsing binaries [26].Simple classical modeling of a system with a linearly increasing mass has been carried out by Kunkel and Harrington [27].A pail-rain model of a linearly massaccreting system is described in [28,29].Here, we consider only a slowly time-varying function for convenience, which is viable under the levy of a condition  ≪ 1.
The momentum operator in three dimensions is defined as where and  is the total angular momentum of the form Notice that the position and the momentum operators satisfy the commutation relation [,   ] = ℏ.
To investigate quantum features of the system, it is necessary to solve the Schrödinger equation in the first place.Let us express the Schrödinger equation for this system in the form Because this equation is a three-dimensional form that involves time-dependent parameters, it may not be an easy task to solve it.Hence, special mathematical techniques are necessary.To solve this equation, we first use the dynamical invariant method devised for evaluating quantum solutions of time-dependent Hamiltonian systems by Lewis and Riesenfeld [22].This method can be developed through an establishment of a nontrivial Hermitian operator () which satisfies the equation In case that such invariant  exists for a system and if it does not contain time derivative operators, the solutions of the Schrödinger equation ( 5) take the form where Φ( ⃗ , ) is an eigenfunction of () that entails a timeindependent eigenvalue and () is a phase of the wave function that can be derived from Advances in Mathematical Physics 3 According to the Lewis-Riesenfeld method which is available for deriving quantum solutions, it is necessary to construct an invariant operator of the system.By solving (6), we have the following invariant operator [30]: where  is a unitary operator of the form and  0 is another invariant quantity that is given by Here,  is the solution of an auxiliary equation which is where  =  0 (1 + ) 2 and  =  0  (/2) .Now, we try to solve (12).The transformation  = (1/(1 + )) leads to the Yermakov's differential equation [31] which is Through the aid of the solution of this equation, we obtain the solution for () to be where  1 and  2 are arbitrary real constants, and () and () are given by where  1 and  2 are arbitrary constants.

Eigenvalues and Eigenfunctions of the Invariant
Let us express the eigenvalue equation for the invariant operator as where  is a time-independent eigenvalue.Because the wave functions for the time-dependent Hamiltonian system are represented in terms of the eigenfunction of the invariant operator, we now solve this equation in order to derive the eigenfunction Φ( ⃗ , ).Through the use of the unitary relation with , (17) can be mapped, in terms of  0 given in (11), into where ( ⃗ ) is the eigenfunction of  0 .Notice that the relation between ( ⃗ ) and Φ( ⃗ , ) is given by Using (11), the eigenvalue equation given in (18) can be written as where  = /,  = tan −1 (/),  = cos −1 (/), and , ,  are the Cartesian coordinates.By taking the function  to be and substituting this into (20), the differential equations for the radial and angular parts are separated as where  is the orbital momentum number.The solution of ( 23) is just the spherical harmonics: By taking the transformation  =  2 , the radial equation leads to the following differential equation: We use the NU method [32] in order to solve this equation in the next section.

The Nikiforov-Uvarov Method
The NU method was proposed to treat complicated secondorder differential equations in a convenient way by reducing it to a hypergeometric-type equation through an appropriate coordinate transformation  = ().To illustrate the general theory of this, let us consider a differential equation which is where () and σ() are polynomials that are allowed up to the second order of  and τ() is a polynomial allowed up to the first order of .If we take the factorization where () is defined through a logarithmic derivative of the form then, according to [33], the function  required for this method is given by where while the other variable () is the hypergeometric-type function whose polynomial solutions are given by the Rodrigues relation: where   are normalizing constants and the weight function () must satisfy the condition with For the case of a general form of the Schrödinger equation, which is given by the value  given in (30) can be represented as [33] where with Again on the basis of the evaluation of [33], the general form of solution for  =  1 = − 1 − 2√ 2 is given by where while, for  =  2 = − 1 + 2√ 2 , the general solution is expressed as where By comparing (34) with ( 25), we get At this stage, we should determine which kind of wave solutions is physically acceptable among the two given in ( 38) and (40), respectively.Because the wave solutions not only should be finite at  = 0 but also should converge to zero at  = ∞, we can confirm that (38) is the only physically suitable kind of wave solutions among them.Then, after a minor evaluation regarding this, the solutions of (25) are given in terms of generalized Laguerre polynomials [34] such that where  = 2√( + 1) + 1/16 + /ℏ 2 , and the eigenvalue  of the invariant is expressed in the form The complete normalized eigenfunctions of the invariant  can be constructed by taking account of normalizing (21) with the use of the following relation [35]: Then, by gathering ( 19), (24), and (43) and introducing an appropriate factor 1/ 3/2 into these functions for the purpose of ensuring the normalization condition, we obtain where the normalization constants are given by Hence, the exact solutions of the original Schrödinger equation given in (5), which is associated with the central harmonic plus an inverse harmonic potential, can now be obtained in the form Ψ (, , , ) = √ 2 +2 !  3 ℏ +1 Γ ( +  + 1)  () (   ) where the global phases are given by These wave functions are useful for investigating quantum properties of the system.For example, the uncertainty relations between a Cartesian coordinate  and its conjugate momentum are obtained through the use of these wave functions, such that Thus, the ground state uncertainty product under the timedependent situation is (ℏ/2)[1 +  2 () 2 ρ 2 ] 1/2 .If the timedependence of the Hamiltonian disappears, the wave functions reduce to the familiar ones known for the timeindependent harmonic oscillator perturbed by the inverse harmonic potential and the uncertainty relations become ΔΔ  = ℏ( + 1/2) as expected.

Conclusion
Quantum solutions of the central time-dependent harmonic plus an inverse harmonic potential system in three dimensions were investigated using several mathematical techniques.According to the theory of Lewis and Risenfeld, the wave functions of the system are represented in terms of the eigenfunctions of the invariant operator.An exact invariant is established on the basis of its fundamental definition via the use of its unitary relation with the unitary operator given in (10).The invariant operator is represented in terms of a time function (), while its overall time derivative results in zero.
The eigenvalue equation for () has been written in terms of a second-order differential equation.To derive the solutions of such equation, the unitary transformation approach and the NU method are used together.The NU method enabled us to reduce the second-order differential equation to the hypergeometric equation.From these procedures, we have obtained the eigenfunctions and they are represented in terms of the Laguerre polynomials.The eigenfunctions of the invariant operator derived in this way are necessary for representing the Schrödinger solutions, that is, the wave functions.The finally obtained wave functions are expressed in terms of these eigenfunctions and time-dependent phase factors, given in (49), that correspond to global phases.The difference of this research from the previous one, performed by Ferkous et al. [30], is that we used the NU method in order to derive exact quantum solutions of the system whereas Ferkous et al. did not.Quantum wave functions that we derived in this work were represented in terms of spherical harmonics    (, ), while those of Ferkous et al. were represented in terms of the Jacobi polynomials.Our results for the wave functions given in (48) are useful for investigating diverse fundamental quantum properties of the system.For example, they can be used to derive fluctuations of canonical variables, uncertainty relations, energy eigenvalues, Wigner functions, and so on.For a particular case where the time-dependence of parameters disappears, the wave functions (48) reduce to the familiar ones that are well known in quantum mechanics of the stationary system.