Quantization of the 1-D forced harmonic oscillator in the space ($x, v$)

The quantization of the forced harmonic oscillator is studied with the quantum variable ($x,\hat v$), with the commutation relation $[x,\hat v]=i\hbar/m$, and using a Shr\"odinger's like equation on these variable, and associating a linear operator to a constant of motion $K(x,v,t)$ of the classical system, The comparison with the quantization in the space ($x,p$) is done with the usual Schr\"odinger's equation for the Hamiltonian $H(x,p,t)$, and with the commutation relation $[x,\hat p]=i\hbar$. It is found that for the non resonant case, both forms of quantization brings about the same result. However, for the resonant case, both forms of quantization are different, and the probability for the system to be in the exited state for the ($x,\hat v$) quantization has less oscillations than the ($x,\hat p$) quantization, the average energy of the system is higher in ($x,\hat p$) quantization than on the $(x,\hat v$) quantization, and the Boltzmann-Shannon entropy on the ($x,\hat p$) quantization is higher than on the ($x,\hat v$) quantization.


Introduction
The usual quantum mechanics formulation is done in the space (x,p) [1], wherep = −i ∂/∂x is the linear operator associated to the classical generalized linear momentum of the motion of a particle of mass "m,"where the commutation relation [x,p] = i [2] is satisfied. A linear operator is associated to the classical Hamiltonian, H(x,p, t), to form the so called Schrödinger's equation [3] i ∂Ψ ∂t = H(x,p, t)Ψ, where Ψ = Ψ(x, t) is the wave function. This formulation has have enormous success to explain and to predict many microscopic behavior of the nature [4]. However, despite this enormous success, Hamiltonian-Lagrangian mathematical formulation has some details, even for 1-D problem where one knows that the Lagrangian (therefore the Hamiltonian) always exists [5]. First, from the expression to obtain the generalized linear momentum given the Lagragian, L(x,ẋ, t) for the system, it is not always possible to obtain explicitlyẋ =ẋ(x, p, t) to be able to get the explicit expression for the Hamiltonian from the Legrandre's transformation [6], H(x, p, t) =ẋ(x, p, t)p − L(x,ẋ(x, p, t), t).
Second, when one is dealing with classical dissipative systems [7], either it is not possible to find its Hamiltonian, or two different Hamiltonians are possible to find for the system [8][9][10][11][12]. Last one, for those problems of variable mass systems, which are not invariant under Galileo's transformations and Sommerfeld modification is not consistent, to find the Hamiltonian for this system [13] requires to start from the "Inverse Problem of the Mechanics." Therefore, one has the necessity to find some extension of the known quantization arised from the Hamilton-Lagrangian approach. In this way, there is already a proposition [14,15] of using a function K(x, v, t) that could be a constant of motion of the classical system, and to associate a linear operator to the velocity of the formv such that [x,v] = i /m, and to associate a linear operator which can be used to form the Shrödinger's like equation In this paper, approaches (1) and (8) are used to study the 1-D forced harmonic oscillator to determine whether or not there is a difference on the quantization, and hopefully to see if the approach (8) could have with these result and experimental verification.
2. Analytical Approach for K(x, v, t) The forced harmonic oscillator is classically characterized by Newton's equation where "m" is the mass of the particle, ω 0 is the natural frequency of oscillation (when α = 0), and α is the amplitude of the forced force. The well known solution of this problem is where one has the non resonant case (ω = ω 0 ) and the resonant case (ω = ω 0 ). The velocity is known by making the differentiation of (10) with respect the time, and the constants C 1 and C 2 are determined by the initial condition (x(0), v(0)). For the non resonant case, these constants are and For the resonant case (ω = ω 0 ), one has and Now, by choosing a constant of motion of the form where 'nr"means non resonant and "r"means resonant, it follows that which represents the usual energy of the harmonic oscillator, independently of the non resonant case or resonant case. This constant of motion can be written as where K 0 and W (nr,r) α are defined as and where one has made the definitions and (19b) To solve equation (8), one observes that the eigenvalues problem for the operator K 0 , has exactly the same solution of the that one given by the Hamiltonian problem, ( (21b) Using Dirac's notation [16], where Φ n (x) = x|n , with |n characterizing the nth-state, and then one has the eigenvalue problem written as Therefore, one can propose the solution of the Shrödinger's equation (8) with the operator constant of motion K, of the form Taking into consideration (22), the orthogonality of the states ( m|n = δ mn ), one obtains the following equation for the coefficients where W (nr,r) mn (t) represents the matrix element The equation (25) can be simplified using the new variable The equations for these new coefficients are where C k (0) = D k (0) and the probability to find the system in the state |k is |C k (t)| 2 = |D k (t)| 2 . Matrix elements are much easier to calculate by using the non Hermitian ascent 'a † " and descent "a" operators, with the knows properties [17] [a, a] = [a † , a † ] = 0, [a, a † ] = 1, and For non resonant case (nr), after calculating the matrix elements , using the orthogonality of the states, and making some rearrangements, one gets the equations for the real and imaginary parts of the coefficients, where a 1 , b 1 , c and d have been defined as to eliminate the quadratic time dependence appearing in the expression (18) and (19b). Note that D k (0) = D k (0) and |D k (t)| 2 = | D k (t)| 2 . Doing the same as it was done above, the real and imaginary parts of these new coefficients, D k = X k + i Y k , obey the equationṡ where the functions f , h, and g have been defined as and The dynamical systems (31) and (34) are solved by Runge-Kutta method a 4th-order.

Analytical Approach for H(x, p, t)
The Hamiltonian of the forced harmonic oscillator is [18] H(x, p, t) = H 0 (x, p) + αx cos(ωt + ϕ), where H 0 is given by The solution of the eigenvalue problem is well known [19], and its solution is the same as (21). Therefore, to solve the Shrödinger's equation (1), one proposes a solution of the form which, after substituting in the Shrödinger's equation, using the eigenvalues and the orthogonality between any two states, and making some rearranging , the following dynamical systems is brought about for the real and imaginary parts of the coefficients, D k (t) = x k (t) + iy k (t), where the constant λ has been defined as These equations are also solved by using Runge-Kutta method at 4th-order.

Boltzmann-Shannon entropy and energy
Besides the probability to find the system in the state |n at the time "t", |D k (t)| 2 , for the analysis of the dynamics of the system in the spaces (x,v) and (x,p), one can also consider the Boltzmann-Shannon entropy, and its average over a evolution time "T",S as parameter which characterize the quantum dynamics of the system. This parameter gives us an indication of how many states enter in the dynamics evolution of the system. Therefore, it gives an indication of the information lost in the dynamics due to the increasing of the entropy in the quantum system.
In addition, one can also consider the expectation value of the energy and its average value over the evolution time of the system, This parameter gives information about how the energy is distributed among the states and how many of them are involved in the quantum dynamics.
In this way, solving the dynamical systems (31), (34, and (40), the evolution of the probabilities |D k (t)| 2 's are gotten. Thus, the Boltzmann-Shannon entropy (42), the expectation value of the energy (44) and their average values (43) and (45) can be calculates and can be compared for the quantization in the spaces (x, v) and (x, p).

Results
One considers a proton with mass m = 1,6726219×10 −27 Kg oscillating with a frequency ω 0 = 2π ×10 9 Hz on a one dimensional line, and interacting with a periodic force force of amplitude α = 10 −13 N ewtons with frequency ω and phase ϕ = 0. The initial conditions of the system are that is, the system is on the ground state, and one selects ten exited possible state of the system.
For the non resonant case (ω = ω 0 ), the resulting dynamics from expressions (1) and (8) are exactly the same. There is not excitation of the system at all since the system remains in the ground state in both cases.
For the resonant case (ω = ω 0 ), Figure 1 shows the probabilities of having the system on the ground state (k = 0) and on the first excited state (k = 1) for the quantization on the space (x, v), solid lines, and the quantization on the space (x, p), dotted lines. As one can see, for the Hamiltonian quantization approach (H) there are much more oscillations of the probabilities than the quantization of the constant of motion approach (K), that is, there are more transitions per unit time in the H-approach case than in the K-approach case. One must note that the probability to have the system in the first excited state for the K-approach case is totally different from the H-approach case.
Figura 1: Ground state and first excited state evolution Figure 2 shows the average value of the energy as a function of the strength of the forced force (α), and as one can see, this value is always higher for H-approach case than for the K-approach case. Although the difference is quite small and maybe out of experimental resolution.  Figure 3 shows the average value of the Boltzmann-Shannon entropy as a function of the strength of the forced force (α). Having total number of 11 states, the possible maximum entropy is 2.398. As the previous case. This parameter is always higher for the H-approach case than for the K-approach case. However, this difference is not so small and maybe could be used as a good parameter for experimental proposes. This difference means that the H-approach case brings about more complex behavior in the quatum dynamics than the K-approach case, and that the H-approach case losses more information than the K-approach case.

Conclusion
The quantization of the 1-D forced harmonic oscillator was carried out with the operators (x,v) using the assigned linear operator to a constant of motion K(x, v, t) of the classical case. The restriction imposed on this constant was that it must reduced to the known energy expression when the forced force is zero. This quantization was compared with the usual quantization with the operators (x,p) and the associated Hamiltonian H(x, p, t) of the classical case. It was shown that the probabilities to find the system in the state |n , |D n (t)| 2 , has less oscillations in the K-quantization than in the H-quantization. In addition, the average values of the energy and the average value of the Boltzmann-Shannon entropy are lower in the K-quantization than in the H-quantization. Since the difference in the average value of the energy is quite small, this parameter does not look good to measure it experimentally. However, the difference in the entropy is significant and it represents a good parameter to look experimentally.