ABSTRACT

We propose an approach based on Floquet theorem combined with the resonating averages method (RAM), to solve the time-dependent Schrödinger equation with a time-periodic Hamiltonian. This approach provides an alternative way to determine directly the evolution operator, and then we deduct the wave functions and the corresponding quasi-energies, of quantum systems. An application is operated for the driven cubic or/and quatric anharmonic as well as for the Morse potential. Comparisons of our results with those of other authors are discussed, and numerical evaluations are performed, to determine the dissociation energy of (HCl) and (CO) molecules.

Keywords:

Floquet Theorem, Resonating Averages, Anharmonic Potential

1. Introduction

Several physical systems have been modeled by a simple harmonic oscillator as a first approximation devoted to the construction of elegant formalism, to understand their dynamics and are dealing with the description of the nature of various physical phenomena under certain conditions. Moreover, it is possible to compare solutions derived from these studies to experimental results obtained by studying laser-matter interaction problems in many branches of physics, ranging from black body radiation to vibrations of crystal lattices [1] [2] [3] [4]. It is now generally realized that interesting connections of numerous real classical or quantum systems are a consequence of the anharmonic non-linear character of evolution in these systems, such as vibrations of the real molecules, graded alloys and irregular semiconductor structures.

In this regard, several approaches have been used, among them the well-known perturbation theory [1] [2] [5] [6], the variational method based on the squeezed states [7], Lie algebraic method using canonical transformations [8] and multiple scale method [9] which have been developed and introduced in some special cases of the time-dependent harmonic and anharmonic oscillators in the presence of an external excitation force, to obtain the evolution operator of these systems which in turn enables to find the explicit expressions of quasi-energies and wave functions. In the regime of strong field-matter interaction, a non-perturbative approach, based on the quantization version of the Floquet formalism [10] [11] [12], has been developed and applied to describe physical systems with time-periodic Hamiltonians. This was used to give explanation of various phenomena, including multiple high order harmonic generation in intense laser field [12], multi-photon ionisation [13], and to describe the motion of a charged particle in an oscillating electric field [14], etc. Our motivation consists especially in solving the problem of a quantum system driven by a periodic time-varying force. The key feature of the application of Floquet theorem is that it permits the reduction of the time-dependent Schrödinger equation of the system to an equivalent conservative eigenvalue problem. It is an efficiency tool in the case of a system submitted to strong time-periodic field, and consequently for which the usual stationary states, solutions to the time-independent Schrödinger equation when no external field acts, do not exist [4]. These derived quantum steady states, lead to set correspondences with stationary states of conservative systems [12] [13].

In previous works, we have applied the above mentioned approach, to the case of harmonic oscillator with time-periodic frequency and to the simple forced harmonic oscillator [15]. Thus, some properties of the Floquet states for this system were elucidated. Moreover, we used resonating averages method (RAM) [16] which provides a useful tool for constructing the evolution operators in a whole resonance zone, which enabled us to obtain readily steady states and the associated quasi-energies of the time-periodic Hamiltonian of these systems. In this work, we apply our method, to first and second ameliorated order approximations, to some driven anharmonic potentials. A comparison of our analytical results with those of the literature was made and numerical evaluations, allowed calculation of the maximum vibrational quantum number to estimate the dissociation energy for the (HCl) and (CO) molecules.

The paper is organized as follows. In Section 2, we review the basic formulation of our approach. Section 3 consists of its application to the driven cubic and quatric anharmonic oscillators, and to the Morse potentail expansion. In Section 4, some comparisons of our results with those of other published works are presented and discussed and an example of numerical evaluations was performed for the (HCl) and (CO) molecules. Concluding remarks are given in Section 5.

2. Basic of the Proposed Approach [15]

A quantum system that is submitted to a perturbation may be described by the following hamiltonian;

$H\left(t\right)=H_0+\lambda H^{\prime }\left(t\right)$ (1)

where $H_0$ is the Hamiltonian of the unperturbed system, and $H^{\prime }\left(t\right)$ is the interaction Hamiltonian, which amplitude $\lambda$ is taken as being very small.

In the case of a periodical perturbation, according to the Floquet theorem there exists a couple of operators (R, T(t)), so that the time-evolution operator can be written in the following form [10] [11] [12]

$U\left(t\right)=T\left(t\right){\text{e}}^{\frac{-iRt}{\hslash }}$ (2)

and $U\left(0\right)=1$

where T(t) is a periodic unitary operator of the same period as $H^{\prime }\left(t\right)$, and R is a constant hermitian operator.

In the interaction picture, $U\left(t\right)$ satisfies the following differential equation

$i\hslash \frac{\text{d}{U}_I\left(t\right)}{\text{d}t}=\lambda H_I\left(t\right)U_I\left(t\right)$ (3)

where

$H_I\left(t\right)={\text{e}}^{i{H}_{0}t/\hslash }{H}^{\prime }\left(t\right){\text{e}}^{-i{H}_{0}t/\hslash }$ (4)

A unitary transformation T(t) may be applied to Equation (3) to obtain the so-called reduced equation of the system such as,

$i\hslash \frac{\text{d}|{\varphi }_n\left(t\right)\rangle }{\text{d}t}=R|{\varphi }_n\left(t\right)\rangle$ (5a)

$|{\varphi }_n\left(t\right)\rangle ={\text{e}}^{-i{E}_{n}t/\hslash }|n\rangle$ (5b)

where $|{\varphi }_n\left(t\right)\rangle$ are the eigenstates of the Floquet operator R, corresponding to the eigenvalues $E_n$, and where $|n\rangle$ are the states of the unperturbed system. Consequently, the Floquet states (or steady states), solutions of the time-dependent Schrödinger equation, in the Floquet representation are defined such as,

$|{\psi }_n\left(t\right)\rangle =T\left(t\right)|{\varphi }_n\left(t\right)\rangle$ (6)

The obtained steady states form a complete set of time-dependent solutions in the extended Hilbert space, and do not depend on the choice of the couple (R, T(t)).

Search of Floquet operators is based on the resonating averages method (RAM) [15] [16], which consists in the separation of the perturbed Hamiltonian $H_I\left(t\right)$, written in the interaction picture Equation (4), into an averaging part, ${\bar{H}}_I\left(t\right)$, and an oscillating part, ${\tilde{H}}_I\left(t\right)$, such as:

$H_I\left(t\right)={\bar{H}}_I\left(t\right)+\frac{\text{d}{\tilde{H}}_I\left(t\right)}{\text{d}t}$ (7)

The application of the RAM to Equation (3) gives rise to the following solutions to first and second order in $\lambda$

${}^{\left(1a\right)}U{}_I\left(t\right)=\left[1-\frac{i\lambda }{\hslash }{\tilde{H}}_I\left(t\right)\right]{}^{\left(1\right)}V{}_I\left(t\right)$ (8)

${}^{\left(2a\right)}U{}_I\left(t\right)=\left[1-\frac{i\lambda }{\hslash }{\tilde{H}}_I\left(t\right)+{\lambda }^2{A}_2\left(t\right)\right]G\left(t\right)$ (9)

The determination of the solutions ${}^{\left(1a\right)}U{}_I\left(t\right)$ and ${}^{\left(2a\right)}U{}_I\left(t\right)$ Equations ((8), (9)) and a comparison with Floquet representation of Equation (2), enabled us to obtain the first and second order couples $\left({}^{\left(1a\right)}R,{}^{\left(1a\right)}T\left(t\right)\right)$,$\left({}^{\left(2a\right)}R,{}^{\left(2a\right)}T\left(t\right)\right)$, and

thence the quasi-energies ${}^{\left(1a\right)}E{}_n$,${}^{\left(2a\right)}E{}_n$ and Floquet-states $|{}^{\left(1a\right)}\psi {}_n\left(t\right)\rangle$,$|{}^{\left(2a\right)}\psi {}_n\left(t\right)\rangle$, respectively.

3. Applications

3.1. Driven Cubic Anharmonic Oscillator

The Hamiltonian of the considered quantum system is given by the following expression

$H_1\left(t\right)=\frac{p^2}{2m}+\frac{1}{2}m{\omega }_0^2{q}^2+{\mu }_1\hslash {\omega }_0{\hat{q}}^3+\mu \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}$ (10)

where $\hat{q}=\sqrt{\frac{m{\omega }_0}{\hslash }}q$

${\mu }_1\hslash {\omega }_0{\hat{q}}^3$ and $\mu \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}$ are the conservative perturbation and the time-dependent perturbation due to an external force, respectively.

m and ${\omega }_{\text{0}}$ are the mass and the frequency of the simple harmonic oscillator, respectively, $\nu$ is the driven oscillation frequency, and ${\mu }_{\text{1}}$,$\mu$, are very weak amplitudes of the perturbations.

Adjusting Equation (10) with the RAM formulation needs the following variables changes,

${\mu }_1=\lambda {\gamma }_1$;$\mu =\lambda \gamma$ and $\lambda \ll 1$

Thus, we can write $H_1\left(t\right)$ in the form,

$H_1\left(t\right)=\frac{p^2}{2m}+\frac{1}{2}m{\omega }_0^2{q}^2+\lambda \left({\gamma }_1\hslash {\omega }_0{\hat{q}}^3+\gamma \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}\right)$ (11)

Introduction of the creation and annihilation operators $a^{+}$ and a of the unperturbed Hamiltonian [15] allows to write $H^{\prime }\left(t\right)$ given in Equation (1) as,

$H^{\prime }\left(t\right)=\frac{\hslash {\omega }_0}{\sqrt{2}}\left[\frac{{\gamma }_1}{2}\left(a^{+3}+a^3+3a^{+}a{a}^{+}+3\left(a^{+}a+1\right)a\right)+\gamma \left(a^{+}+a\right)\sin \left(\nu t\right)\right]$ (12)

The RAM applied to the interaction picture form of $H_1\left(t\right)$ (Equations ((4), (7))) gives

${\bar{H}}_I\left(t\right)=0$ (13)

$\begin{array}{c}{\tilde{H}}_I\left(t\right)=\frac{i{\gamma }_1\hslash }{6\sqrt{2}}\left[-{\text{e}}^{3i{\omega }_{0}t}{a}^{+3}+{\text{e}}^{-3i{\omega }_{0}t}{a}^3-9{\text{e}}^{i{\omega }_{0}t}{a}^{+}a{a}^{+}+9{\text{e}}^{-i{\omega }_{0}t}\left(a^{+}a+1\right)a\right]\\ +\frac{\gamma \hslash {\omega }_0}{\sqrt{2}}\left({\alpha }_0^{\star }\left(t\right){\text{e}}^{i{\omega }_{0}t}{a}^{+}+{\alpha }_0\left(t\right){\text{e}}^{-i{\omega }_{0}t}a\right)\end{array}$ (14)

where

${\alpha }_0\left(t\right)=-\frac{\nu \cos \left(\nu t\right)+i{\omega }_0\sin \left(\nu t\right)}{{\nu }^2-{\omega }_0^2}$ (15)

${\alpha }_0^{\star }\left(t\right)$ being the complex conjugate of ${\alpha }_0(t)$

From Equation (8), we obtain the evolution operator to the first ameliorated order such as

$\begin{array}{c}{}^{\left(1a\right)}U{}_1\left(t\right)=[1+\frac{{\mu }_1}{6\sqrt{2}}\left(-a^{+3}+a^3-9a^{+}a{a}^{+}+9\left(a^{+}a+1\right)a\right)\\ -i\frac{\mu {\omega }_0}{\sqrt{2}}\left({\alpha }_0^{\star }\left(t\right)a^{+}+{\alpha }_0\left(t\right)a\right)]{\text{e}}^{\frac{-i{H}_{0}t}{\hslash }}\end{array}$ (16)

By comparison with the formulation of Equation (2), we deduce the first ameliorate order Floquet operators (R, T(t)), such as:

${}^{\left(1a\right)}R{}_1=\hslash {\omega }_0\left(a^{+}a+\frac{1}{2}\right)$ (17a)

$\begin{array}{c}{}^{\left(1a\right)}T{}_1\left(t\right)=1+\frac{{\mu }_1}{6\sqrt{2}}\left(-a^{+3}+a^3-9a^{+}a{a}^{+}+9\left(a^{+}a+1\right)a\right)\\ -i\frac{\mu {\omega }_0}{\sqrt{2}}\left({\alpha }_0^{\star }\left(t\right)a^{+}+{\alpha }_0\left(t\right)a\right)\end{array}$ (17b)

Thus, the quasi-energies, Floquet states and wave functions developed to the first order respectively are:

${}^{\left(1a\right)}E{}_n=\hslash {\omega }_0\left(n+\frac{1}{2}\right)$ (18)

$|{}^{\left(1a\right)}\psi {}_n\left(t\right)\rangle ={\text{e}}^{\frac{-i{}^{\left(1a\right)}E{}_{n}t}{\hslash }}\left[c_{-3}|n-3\rangle +c_{-1}|n-1\rangle +|n\rangle +c_{+1}|n+1\rangle +c_{+3}|n+3\rangle \right]$ (19)

${}^{\left(1a\right)}\psi {}_n\left(q,t\right)={\text{e}}^{\frac{-i{}^{\left(1a\right)}E{}_{n}t}{\hslash }}\left[c_{-3}{\phi }_{n-3}\left(q\right)+c_{-1}{\phi }_{n-1}\left(q\right)+{\phi }_n\left(q\right)+c_{+1}{\phi }_{n+1}\left(q\right)+c_{+3}{\phi }_{n+3}\left(q\right)\right]$ (20)

where the coefficients $c_{\pm 1}$ and $c_{\pm 3}$ are given by

$c_{-3}=\frac{{\mu }_1}{6\sqrt{2}}\sqrt{n\left(n-1\right)\left(n-2\right)}$ (21a)

$c_{-1}=\left[\frac{3{\mu }_1}{2\sqrt{2}}n-i\frac{\mu {\omega }_0}{\sqrt{2}}{\alpha }_0\left(t\right)\right]\sqrt{n}$ (21b)

$c_{+1}=-\left[\frac{3{\mu }_1}{2\sqrt{2}}\left(n+1\right)+i\frac{\mu {\omega }_0}{\sqrt{2}}{\alpha }_0^{\star }\left(t\right)\right]\sqrt{n+1}$ (21c)

$c_{+3}=-\frac{{\mu }_1}{6\sqrt{2}}\sqrt{\left(n+1\right)\left(n+2\right)\left(n+3\right)}$ (21d)

and

${\phi }_n\left(q\right)={\left(\frac{\alpha }{\pi }\right)}^{\frac{1}{4}}\frac{{\text{e}}^{\frac{-\alpha q^2}{2}}}{\sqrt{2^{n}n!}}{H}_n\left(\sqrt{\alpha }q\right)$ (22)

${\phi }_n\left(q\right)$ is the wave function of the simple oscillator, the parameter $\alpha =\frac{m{\omega }_0}{\hslash }$ and $H_n\left(\sqrt{\alpha }q\right)$ are the Hermite polynomials.

Using Equation (9) and with the help of Equations ((13), (14)) one can write the second ameliorated order evolution operator and deduce the Floquet operators ${}^{\left(2a\right)}R{}_1$ in the following form

${}^{\left(2a\right)}R{}_1=H_0-{\mu }_1^2\frac{15\hslash {\omega }_0}{4}\left({\left(a^{+}a\right)}^2+a^{+}a+\frac{11}{30}\right)+\frac{{\mu }^2\hslash {\omega }_0^3}{4\left({\nu }^2-w_0^2\right)}$ (23)

Thus, the quasi-energies of the system up to second ameliorated order is given by

${}^{\left(2a\right)}E{}_n=\hslash {\omega }_0\left(n+\frac{1}{2}\right)-{\mu }_1^2\frac{15\hslash {\omega }_0}{4}\left(n^2+n+\frac{11}{30}\right)+\frac{{\mu }^2\hslash {\omega }_0^3}{4\left({\nu }^2-w_0^2\right)}$ (24)

We note that the correction effects on the quasi-energies to second order approximation, and the Floquet shift levels depend on the amplitudes $\left({\mu }_1,\mu \right)$, of the perturbations and the quantum number n.

3.2. Driven Quatric Anharmonic Oscillator

We consider the system which Hamiltonian is given by,

$H_2\left(t\right)=\frac{p^2}{2m}+\frac{1}{2}m{\omega }_0^2{q}^2+{\mu }_2\hslash {\omega }_0{\hat{q}}^4+\mu \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}$ (25)

where

${\mu }_2\hslash {\omega }_0{\hat{q}}^4$ and $\mu \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}$ are the quatric anharmonic perturbation with amplitude ${\mu }_{\text{2}}$ and the external time-dependent perturbation with amplitude $\mu$ respectively.

We operate the change variable ${\mu }_2=\lambda {\gamma }_2$ and $\mu =\lambda \gamma$ on Equation (25), then we have,

$H_2\left(t\right)=\frac{p^2}{2m}+\frac{1}{2}m{\omega }_0^2{q}^2+\lambda \left({\gamma }_2\hslash {\omega }_0{\hat{q}}^4+\gamma \hslash {\omega }_0\sin \left(\nu t\right)\hat{q}\right)$ (26)

Using the usual creation and annihilation operators yields to write $H^{\prime }\left(t\right)$ in Equation (1) as

$\begin{array}{c}{H}^{\prime }\left(t\right)={\gamma }_2\frac{\hslash {\omega }_0}{4}[a^{+4}+a^4+4a^{+}a{a}^{+2}-2a^{+2}+4\left(a^{+}a+1\right)a^2+2a^2\\ +3{\left(a^{+}a\right)}^2+3{\left(a^{+}a+1\right)}^2]+\gamma \frac{\hslash {\omega }_0}{\sqrt{2}}\left(a^{+}+a\right)\sin \left(\nu t\right)\end{array}$ (27)

The RAM applied to the interaction picture form of $H_2\left(t\right)$ (Equations ((4), (7))), gives

${\bar{H}}_I\left(t\right)=\frac{3}{4}{\gamma }_2\hslash {\omega }_0\left(2{\left(a^{+}a\right)}^2+2a^{+}a+1\right)$ (28)

$\begin{array}{c}{\tilde{H}}_I\left(t\right)={\gamma }_2\frac{i\hslash }{16}[-{\text{e}}^{4i{\omega }_{0}t}{a}^{+4}+{\text{e}}^{-4i{\omega }_{0}t}{a}^4-4{\text{e}}^{2i{\omega }_{0}t}\left(2a^{+}a-1\right)a^{+2}\\ +4{\text{e}}^{-2i{\omega }_{0}t}\left(2a^{+}a+3\right)a^2]+\gamma \frac{\hslash {\omega }_0}{\sqrt{2}}\left({\alpha }_0^{\star }\left(t\right){\text{e}}^{i{\omega }_{0}t}{a}^{+}+{\alpha }_0\left(t\right){\text{e}}^{-i{\omega }_0}a\right)\end{array}$ (29)

where ${\alpha }_0\left(t\right)$ is given by Equation (15).

Following the previous procedure given in the subsection 3.1, we obtain the time evolution operator, then the quasi-energies and Floquet states to the first ameliorated order, for this system, respectively such as,

$\begin{array}{c}{}^{\left(1a\right)}U{}_2\left(t\right)=[1-\frac{i\lambda }{\hslash }[{\gamma }_2\frac{i\hslash }{16}\left[-a^{+4}+a^4-4\left(2a^{+}a-1\right)a^{+2}+4\left(2a^{+}a+3\right)a^2\right]\\ +\gamma \frac{\hslash {\omega }_0}{\sqrt{2}}\left({\alpha }_0^{\star }\left(t\right)a^{+}+{\alpha }_0\left(t\right)a\right)]]\\ \times \exp \left(\frac{-it}{\hslash }\left[\frac{3}{4}{\mu }_2\hslash {\omega }_0\left(2{\left(a^{+}a\right)}^2+2a^{+}a+1\right)+H_0\right]\right)\end{array}$ (30)

${}^{\left(1a\right)}E{}_n=\hslash {\omega }_0\left(n+\frac{1}{2}\right)+\frac{3}{4}{\mu }_2\hslash {\omega }_0\left(2n^2+2n+1\right)$ (31)

$\begin{array}{c}|{}^{\left(1a\right)}\psi {}_n\left(t\right)\rangle ={\text{e}}^{\frac{-i{}^{\left(1a\right)}E{}_{n}t}{\hslash }}[k_{-4}|n-4\rangle +k_{-2}|n-2\rangle +k_{-1}|n-1\rangle \\ +|n\rangle +k_{+1}|n+1\rangle +k_{+2}|n+2\rangle +k_{+4}|n+4\rangle ]\end{array}$ (32)

where

$k_{-4}=\frac{{\mu }_2}{16}\sqrt{n\left(n-1\right)\left(n-2\right)\left(n-3\right)}$ (33a)

$k_{-2}=\frac{{\mu }_2}{4}\left(2n-1\right)\sqrt{n\left(n-1\right)}$ (33b)

$k_{-1}=-i\frac{\mu {\omega }_0}{\sqrt{2}}{\alpha }_0\left(t\right)\sqrt{n}$ (33c)

$k_{+1}=-i\frac{\mu {\omega }_0}{\sqrt{2}}{\alpha }_0^{\star }\left(t\right)\sqrt{n+1}$ (33d)

$k_{+2}=-\frac{{\mu }_2}{4}\left(2n+3\right)\sqrt{\left(n+1\right)\left(n+2\right)}$ (33e)

$k_{+4}=-\frac{{\mu }_2}{16}\sqrt{\left(n+1\right)\left(n+2\right)\left(n+3\right)\left(n+4\right)}$ (33f)

Application of the RAM Equation (9) to second ameliorated order, gives the expressions of the Floquet operator and quasi-energies respectively, such as:

$\begin{array}{c}{}^{\left(2a\right)}R{}_2=\hslash {\omega }_0\left(a^{+}a+\frac{1}{2}\right)+{\mu }_2\frac{3}{4}\hslash {\omega }_0\left(2{\left(a^{+}a\right)}^2+2a^{+}a+1\right)\\ -{\mu }_2^2\frac{\hslash {\omega }_0}{8}\left(34{\left(a^{+}a\right)}^3+51{\left(a^{+}a\right)}^2+59a^{+}a+21\right)+\frac{{\mu }^2\hslash {\omega }_0^3}{4\left({\nu }^2-w_0^2\right)}\end{array}$ (34)

$\begin{array}{c}{}^{\left(2a\right)}E{}_n=\hslash {\omega }_0\left(n+\frac{1}{2}\right)+{\mu }_2\frac{3}{4}\hslash {\omega }_0\left(2n^2+2n+1\right)\\ -{\mu }_2^2\frac{\hslash {\omega }_0}{8}\left(34n^3+51n^2+59n+21\right)+\frac{{\mu }^2\hslash {\omega }_0^3}{4\left({\nu }^2-w_0^2\right)}\end{array}$ (35)

We note that the correction on quasi-energies of this system exist to first and second orders, and the Floquet shift levels depends on the parameters $\left({\mu }_2,\mu \right)$ and the quantum number n, which means that these energies levels are not equidistants.

We also note that in the absence of the cubic and quatric anharmonic perturbations ( ${\mu }_1={\mu }_2=0$ ) we find the Floquet states and quasi-energies of the simple forced harmonic oscillator [15].

3.3. Morse Potential Expansion

The Morse potential $V\left(q\right)$, is the agreed model for diatomic molecules and is given by [6] [7],

$V\left(q\right)=D_e{\left(1-{\text{e}}^{-\rho q}\right)}^2$ (36)

where $\rho$ is a parameter that controls the width of attraction (or measure the curvature at the bottom of the well), $D_e$ is a parameter that controls the depth of the attracting well (or the depth of potential well), and q is the variation of interatomic distance with respect to the equilibrium distance.

Let us consider the Taylor development of exponential term in $V\left(q\right)$ to forth order, and collecting terms in this development gives,

$V\left(q\right)=D_e{\rho }^2{q}^2-D_e{\rho }^3{q}^3+\frac{7}{12}{D}_e{\rho }^4{q}^4$ (37)

Then we make the following changes to the notations:

${\omega }_0=\rho \sqrt{\frac{2D_e}{m}}$ (38a)

${\mu }_1=-\frac{D_e{\rho }^3}{m{\omega }_0^2}\sqrt{\frac{\hslash }{m{w}_0}}$ (38b)

${\mu }_2=\frac{7}{12}\frac{D_e{\rho }^4\hslash }{m^2{\omega }_0^3}$ (38c)

where ${\omega }_0$ is the vibrational constant with the reduced mass m of the diatomic molecule.

With the help of Equations ((36), (37)), we obtain the similar situation given, by the cubic together with quatric anharmonic oscillator, in the absence of the time dependent perturbation ( $\mu =0$ ), for which the Hamiltonian is given by,

$H_3\left(t\right)=\frac{p^2}{2m}+D_e{\rho }^2{q}^2-D_e{\rho }^3{q}^3+\frac{7}{12}{D}_e{\rho }^4{q}^4$ (39)

Using Equations (24) and (35), one can obtain the full quasi-energy to second order as:

$\begin{array}{c}{}^{\left(2a\right)}E{}_n=\hslash {\omega }_0\left(n+\frac{1}{2}\right)+\frac{7}{16}\frac{D_e{\rho }^4{\hslash }^2}{m^2{w}_0^2}\left(2n^2+2n+1\right)\\ -\frac{7}{16}\frac{D_e^2{\rho }^6{\hslash }^2}{m^3{w}_0^4}\left(n^2+n+\frac{11}{30}\right)\\ -\frac{49}{1152}\frac{D_e^2{\rho }^8{\hslash }^3}{m^4{\omega }_0^5}\left(34n^3+51n^2+59n+21\right)\end{array}$ (40)

The difference between two adjacent Floquet levels for the cubic and quatric anharmonic oscillators are given by,

$\begin{array}{c}\Delta E={}^{\left(2a\right)}E{}_{n+1}-{}^{\left(2a\right)}E{}_n\\ =\hslash {\omega }_0\left[1+3{\mu }_2\left(n+1\right)-{\mu }_1^2\frac{15}{2}\left(n+1\right)-{\mu }_2^2\frac{51}{4}\left({\left(n+1\right)}^2+\frac{7}{17}\right)\right]\end{array}$ (41)

We note that when n increases $\Delta E$ decreases, until becoming equal to zero when the energy level reached the dissociation energy of diatomic molecular system, then the quantum number takes the maximum value $n_{\max }$.

4. Comparisons and Numerical Evaluations

In the previous paragraphs we have developed calculations to first and second orders, and presented a number of results of the quantum anharmonic oscillator (Floquet states, quasi-energies). Table 1 compares our results with the works, of some other published works that used perturbation theory method.

The coefficients $A\left(n\right)$,$B\left(n\right)$,$C\left(n\right)$,$D\left(n\right)$,$E\left(n\right)$,$F\left(n\right)$,$G\left(n\right)$,$H\left(n\right)$,$A^{\left(2\right)}\left(n\right)$,$B^{\left(1\right)}\left(n\right)$ and $B^{\left(2\right)}\left(n\right)$, are given by Equation (15) in page (77) of the reference [5].

We observe in Table 1 that our first order states of the anharmonic (cubic and quatric) oscillator are similar to those obtained by the application of the stationary perturbation theory given in the work of Wang et al. [5]. We also note that our quasi-energy expression to second order for the anharmonic oscillator agreed with those obtained by Wang et al. [5] for this system.

We express ${}^{\left(2a\right)}E{}_n$ of the Equation (40) as a function of $\left(n+\frac{1}{2}\right)$ and since the experimental data show that the terms in ${\left(n+\frac{1}{2}\right)}^3$ are negligible [6], we obtain

${}^{\left(2a\right)}E{}_n\approx \hslash \rho \sqrt{\frac{2D_e}{m}}\left(1-\frac{3283}{18432}\frac{{\hslash }^2{\rho }^2}{m{D}_e}\right)\left(n+\frac{1}{2}\right)-\frac{{\hslash }^2{\rho }^2}{2m}{\left(n+\frac{1}{2}\right)}^2$ (42)

Notations: $m\to m^{\star }$,$\mu =0$,${\mu }_1=\frac{{\alpha }_w}{\hslash {\omega }_0}$,${\mu }_2=\frac{{\beta }_w}{\hslash {\omega }_0}$,$\hat{q}=\frac{z}{L}$.

Let us consider the values of the parameters corresponding of (HCl) and (CO) diatomic molecules given in Table 2,

The derivative of Equation (42) with respect to n and setting it equal to zero gives the vibrational level associated with the dissociation limit $n_{\max }$, which allowed us to find the values of the theoretical dissociation energy corresponding to the equilibrium point ( $D_e^{th}$ ) and the dissociation energy with respect to the zero-point level ( $D_0=D_e^{th}-E_0$ ) which are given in Table 3.

5. Conclusions

The Floquet theory is one of the most useful tools which provides an alternative way for solving the Schrödinger equation of quantum systems with time-periodic Hamiltonian. In this paper, this approach was applied to driven quantum anharmonic oscillators. We have given the Floquet operators, solutions of the Schrödinger evolution equation, with the help of the RAM applied to first and to second ameliorated orders approximation and then we have calculated the Floquet states and the corresponding quasi-energies as well as the wave functions. Indeed, the approach used in our study determined, in a natural way, the explicit expressions for the time-dependent states of the anharmonic potential systems. It can be noticed that when we switch off the time-perturbation, we obtain the conservative energies of the cubic and the quatric anharmonic potential, and that the energy levels spacing decrease with increasing values of n and allows us to estimate the dissociation energy of the molecule.

The comparisons of our expressions with published works by other authors, which have used different methods [1] [2] [5] revealed a good concordance, and the numerical evaluation carried out for the values of the parameters of (HCl) and (CO) diatomic molecules, illustrated clearly our results. This approach can be a useful tool to solve the Schrödinger equation of other types of driven time-dependent quantum systems. Therefore, it can be applied to investigate transitions between excited states and evaluate the dissociation energy of diatomic and polyatomic molecules. The goal of future work will be the application of the established approach to the driven Mathieu oscillator.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Idrissi, M.J., Fedoul, A., Sayouri, S. and Amila, I. (2020) Anharmonic Potentials Analysis through the Floquet Representation. Journal of Applied Mathematics and Physics, 8, 184-195. https://doi.org/10.4236/jamp.2020.81014

References

Appendix

$H\left(t\right)$ : Time-dependent Hamiltonian.

$H_0$ : Unperturbed Hamiltonian.

$H^{\prime }\left(t\right)$ : Interaction Hamiltonian.

$\lambda$ : Amplitude of the perturbation.

$\hslash$ : Reduced Planck's constant.

$|{\varphi }_n\left(t\right)\rangle$ : Eigenstates of the operator R.

$E_n$ : Eigenvalues of R (Quasi-energies).

$|n\rangle$ : Sates of the unperturbed system.

$|{\psi }_n\left(t\right)\rangle$ : Floquet states.

${}^{\left(2a\right)}U\left(t\right)$ : Second order ameliorated evolution operator.

p: Impulsion operator.

q: Position operator.

m: Mass of the system

${\omega }_{\text{0}}$ : Unperturbed oscillator frequency.

${\mu }_{\text{1}}$ : Amplitude of the cubic anharmonic oscillator.

$\mu$ : Amplitude of the time-dependent perturbation.

$\nu$ : Frequency of the time-dependent perturbation.

a: Annihilation operator.

$a^{+}$ : Creation operator.

${\mu }_{\text{2}}$ : Amplitude of the quatric anharmonic oscillator.

$D_e$ : Depth of potential well.

$\rho$ : Parameter that controls the width of attraction.

$D_e^{th}$ : Theoretical dissociation energy corresponding to the equilibrium point.

$E_0$ : Zero-point energy (n = 0).

$D_0$ : Dissociation energy with respect to the zero-point level.