Next Article in Journal / Special Issue
Solutions of the Multivariate Inverse Frobenius–Perron Problem
Previous Article in Journal
Modified Hilbert Curve for Rectangles and Cuboids and Its Application in Entropy Coding for Image and Video Compression
Previous Article in Special Issue
Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators

1
Laborarory of Chemical Materials, Department of Physics, Faculty of Sciences, University of Mohamed Boudiaf M’sila, M’sila 28000, Algeria
2
Laboratory of Optoelectronics and Compounds (LOC), Departement of Physics, Faculty of Science, University of Ferhat Abbas Setif 1, Setif 19000, Algeria
3
Department of Nanoengineering, Kyonggi University, Yeongtong-gu, Suwon 16227, Korea
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(7), 837; https://doi.org/10.3390/e23070837
Submission received: 30 May 2021 / Revised: 21 June 2021 / Accepted: 26 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Dynamical Systems, Differential Equations and Applications)

Abstract

:
We study the dynamical invariant for dissipative three coupled oscillators mainly from the quantum mechanical point of view. It is known that there are many advantages of the invariant quantity in elucidating mechanical properties of the system. We use such a property of the invariant operator in quantizing the system in this work. To this end, we first transform the invariant operator to a simple one by using a unitary operator in order that we can easily manage it. The invariant operator is further simplified through its diagonalization via three-dimensional rotations parameterized by three Euler angles. The coupling terms in the quantum invariant are eventually eliminated thanks to such a diagonalization. As a consequence, transformed quantum invariant is represented in terms of three independent simple harmonic oscillators which have unit masses. Starting from the wave functions in the transformed system, we have derived the full wave functions in the original system with the help of the unitary operators.

1. Introduction

Description and interpretation of coupled systems are of particular interest in physics because the interaction caused by coupling is responsible for novel quantum effects such as entanglement [1,2] and quadrature squeezing [3]. Coupled oscillatory quantum motions are found in almost all areas of physical sciences, ranging from nanotechnology to biology [4,5,6,7,8]. Coupled oscillatory systems can be used as a model to describe the interactions between atoms in a one-dimensional crystal with spring-like forces under white noise excitations [9,10]. In this regard, oscillatory motions of three masses coupled by four springs were studied analytically by Ndikilar et al. [9]. The dynamics of mixedness and entanglement for three coupled oscillators with arbitrary time-dependent frequencies and coupling parameters has also been investigated [11]. Entanglement and its control in coupled quantum oscillators are crucial for the process of information in quantum cryptography, quantum teleportation, and quantum computing [12].
In particular, coupled oscillator model for nano-optomechanical systems is important because many novel state-of-the-art techniques can be realized by utilizing it. For instance, couplings of photonic systems with mechanical resonators provide fundamental platform for quantum technologies, such as slow/fast-light generation [13], cooling of nanomechanical resonators [14,15], frequency conversion [16,17], and phononic-structure preparation [18]. You can see an example of nano-optomechanical three coupled oscillators from Figure 1.
In this contribution, we will study quantum dynamical invariant for dissipative three coupled oscillators based on exact quantum description of the system. Such invariant can be used in analyzing various mechanical properties of the coupled oscillatory systems [19,20]. For example, quantum mechanics of optomechanical systems can be unfolded by means of the theory of dynamical invariants [20].
There are several methods to calculate dynamical invariants based on algebraic approaches [21,22,23,24,25,26,27]. In the classical case, the Lutzky’s approach [23,24], which relies on the Noether’s theorem is applicable. Some years ago, Bertin et al. developed a new way to calculate dynamical invariants [25], which utilizes the combinations of the classical equations of motion. We also note that there are Lewis-Riesenfeld method for time-dependent Hamiltonian systems [26,27]. In particular, the last method can be flexibly used in both classical and quantum systems. Many research papers were devoted to the construction and applications of dynamical invariants in nonconservative systems [28,29,30]. The method of invariants was also used in the study of three coupled oscillators [31,32,33].
Our paper is structured as follows: In Section 2, we will represent the method for treating three coupled oscillators from preliminary level of mechanics. The results of this work and the related discussion will be placed in Section 3. The classical and quantum invariant quantities for the three coupled oscillators will be derived based on the fundamental Hamiltonian dynamics. The quantum invariant operator will be simplified by a unitary transformation together with a diagonalization through a rotational unitary transformation parameterized by Euler angles. Then, we evaluate the eigenvalues and eigenfunctions of the original invariant operator by taking advantage of the simple expression of the diagonalized invariant operator. We will also derive solutions of the Schrödinger equation of the oscillatory systems using their close relationship with the eigenstates of the invariant operator. Finally, we will conclude our research in Section 4 with some remarks.

2. Methods

We start our study by briefly representing how to describe the dynamics of a simple dissipative mechanical oscillatory system. For the linearly damped motion of a 1D oscillatory system subjected to a force field V ( X ) , the Newton’s equation can be written as
m X ¨ + γ X ˙ = V ( X ) .
The Lagrangian that produces this motion is of the form
L = exp ( δ t ) 1 2 m X ˙ 2 V ( X ) , δ = γ m ,
where δ is the damping coefficient. From the well-known relation between the Lagrangian and the Hamiltonian, we can easily have the corresponding Hamiltonian such that
H = P 2 2 m exp δ t + V ( X ) exp ( δ t ) .
This is the famous Caldirola–Kanai model [34,35], which celebrates 80 years now. This Hamiltonian is dependent on time due to the damping of the system, even if the mass is independent of time.
The above simple mechanical description can be readily extended to coupled oscillators. We consider a system of dissipative three coupled oscillators which have different effective masses ( m 1 , m 2 , m 3 ), where they are parameterized by the three coordinates X 1 , X 2 , X 3 [36]. The Hamiltonian of this system can be represented as
H ( t ) = 1 2 i = 1 3 P i 2 e δ t m i + e δ t k i X i 2 + 1 2 e δ t k 12 X 1 X 2 2 + k 13 X 1 X 3 2 + k 23 X 2 X 3 2 ,
where the parameters k i ( i = 1 , 2 , 3 ) , k 12 , k 13 , and k 23 are constants; the convention for the lower subscript i (including j and k) given here will be applied throughout the paper. The Hamiltonian in Equation (4) is actually a generalization of 1D single mechanical system expressed by Equation (3). We assumed that damping coefficients for the three oscillators are the same as each other for simplicity; this means that
γ 1 / m 1 = γ 2 / m 2 = γ 3 / m 3 δ .
Let us denote the Poisson bracket for two arbitrary observables μ i and ν j as { μ i , ν j } . Then, for canonical local coordinates, we have
{ X i , P j } = δ i j , X i , X j = P i , P j = 0 .
From the use of the Hamilton’s equations, X ˙ i = H / P i and P ˙ i = H / X i , we confirm that the classical equations of motion are given by
X ¨ 1 X ¨ 2 X ¨ 3 = δ X ˙ 1 δ X ˙ 2 δ X ˙ 3 + k 1 + k 12 + k 13 m 1 k 12 / m 1 k 13 / m 1 k 12 / m 2 k 2 + k 12 + k 23 m 2 k 23 / m 2 k 13 / m 3 k 23 / m 3 k 3 + k 13 + k 23 m 3 X 1 X 2 X 3 .
These equations reveal that coordinates of the system are intricate due to the coupling between mechanical oscillators. Hence, the investigation of mechanical properties of the system is not an easy task. Nevertheless, it may be possible to overcome such a knotty situation by finding a classical invariant quantity of the system. An invariant quantity is an important tool for analyzing mechanical properties of dynamical systems [37,38,39,40]. Let us assume that the formula of the invariant quantity for our system is of the form
I ( t ) = 1 2 i = 1 3 α i t P i 2 + γ i ( t ) X i 2 + β i t X i P i + 1 2 η 12 ( t ) X 1 X 2 + η 13 t X 1 X 3 + η 23 ( t ) X 2 X 3 ,
where α i ( t ) , β i ( t ) , γ i ( t ) , and η i j ( t ) are real differentiable functions of time, of which formulae will be evaluated later on.
In fact, the search for such an invariant (a constant of motion) is somewhat difficult for a complicated system. To obtain the invariant quantity by determining the time functions in Equation (8), we use the Liouville’s theorem which is that the phase-space distribution function is constant along the trajectories of the dynamical systems. Via the analysis of the invariant quantity based on this theorem, we can confirm the time evolution of a given dynamical system. Usually, there exist constants of motion in addition to the energy for an integrable dynamical system. These constants for a time-dependent Hamiltonian system do not commute with the Hamiltonian under the Poisson bracket.

3. Results and Discussion

3.1. Classical Analysis

If we insert Equation (8) with Equation (4) into the Liouville equation
d I d t = I t + i = 1 3 I X i H P i I P i H X i = 0 ,
we obtain the coupled differential equations for α i ( t ) , β i ( t ) , γ i ( t ) , and η i j ( t ) as
α ˙ i ( t ) = 2 β i ( t ) m i e δ t ,
β ˙ 1 ( t ) = k 1 + k 12 + k 13 e δ t α 1 ( t ) γ 1 ( t ) m 1 e δ t ,
β ˙ 2 ( t ) = k 2 + k 12 + k 23 e δ t α 2 ( t ) γ 2 ( t ) m 2 e δ t ,
β ˙ 3 ( t ) = k 3 + k 13 + k 23 e δ t α 3 ( t ) γ 3 ( t ) m 3 e δ t ,
γ ˙ 1 ( t ) = 2 k 1 + k 12 + k 13 e δ t β 1 ( t ) ,
γ ˙ 2 ( t ) = 2 k 2 + k 12 + k 23 e δ t β 2 ( t ) ,
γ ˙ 3 ( t ) = 2 k 3 + k 13 + k 23 e δ t β 3 ( t ) ,
η ˙ 12 ( t ) = 2 k 12 ( t ) e δ t β 1 t + β 2 ( t ) ,
η ˙ 13 ( t ) = 2 k 13 ( t ) e δ t β 1 t + β 3 ( t ) ,
η ˙ 23 ( t ) = 2 k 23 ( t ) e δ t β 2 t + β 3 ( t ) .
By solving the above coupled equations, we have
α i ( t ) = 1 m i e δ t , β i ( t ) = δ 2 ,
γ 1 ( t ) = k 1 + k 12 + k 13 e δ t , γ 2 ( t ) = k 2 + k 12 + k 23 e δ t , γ 3 ( t ) = k 3 + k 13 + k 23 e δ t ,
η 12 ( t ) = 2 k 12 e δ t , η 13 ( t ) = 2 k 13 e δ t , η 23 ( t ) = 2 k 23 e δ t .
From the substitution of these outcomes into Equation (8), we have the classical invariant quantity as
I ( t ) = H ( t ) + δ 2 i = 1 3 X i P i .
This can be used in analyzing dynamical properties of the oscillatory systems.
Besides the invariant approach that we are interested here, there are other treatments with Hamiltonian for dissipative harmonic oscillators. Bateman carried out a variational approach for the dissipative oscillators with the Langrangian that gives suitable equations of motion of the system [41]. Lemos made the Hamiltonian of the dissipative system simpler through a canonical transformation based on the Hamilton–Jacobi method [42]. Later on, McDonald also proposed another linear (canonical) transformation of position and momentum variables for the damped oscillator, which holds the Liouville’s theorem [43]. These researches may serve as a contrast to the current treatment and would also highlight potential novelties along this line.
In the next section, we will extend this invariant to quantum mechanics, i.e., we will consider its counterpart quantum description in order to unfold the associated quantum theory based on the dynamical invariant.

3.2. Quantum Analysis

Due to the analogy of quantum mechanics with the classical mechanics, we can still use the similar notions of an observable and the associated physical state related to our previous description of the system even in a quantum domain. By replacing the canonical variables in the classical invariant with corresponding quantum operators, one gets the quantum invariant operator, such that
I ^ ( t ) = H ^ ( t ) + δ 4 i = 1 3 ( X ^ i P ^ i + P ^ i X ^ i ) ,
where the momentum coordinate is given by P ^ i = i X i . Through the re-interpretation of the classical variables as counterpart quantum operators, Poisson brackets are replaced by commutators [ X ^ i , P ^ j ] = i δ i j . Then, the analog of Liouville equation in quantum mechanics is the von Neumann one, which is of the form
d I ^ ( t ) d t = I ^ ( t ) t + 1 i [ I ^ ( t ) , H ^ ( t ) ] = 0 .
This can be used to elucidate the time evolution of a quantum state of the system. This consequence stems from the fact that the canonical quantization of the system is possible on the basis of the theorem related to the invariant operator [26,27]. However, for the practical use of the invariant operator for such a purpose, it may be advantageous to transform it to a simple form.
To simplify the invariant, Equation (24), we use the unitary transformation approach. As a first step, we introduce the following transformation of it:
I ^ = U ^ 1 ( t ) I ^ ( t ) U ^ t ,
where the unitary operator U ^ ( t ) is given by [36]
U ^ ( t ) = i = 1 3 exp i 2 ( P ^ i X ^ i + X ^ i P ^ i ) ln m i + δ 2 t × exp i δ 4 i = 1 3 X ^ i 2 .
Then, after an algebraic evaluation, it is possible to get the transformed invariant as
I ^ = 1 2 P ^ 1 2 + P ^ 2 2 + P ^ 3 2 + 1 2 k 1 + k 12 + k 13 m 1 δ 4 2 X ^ 1 2 + 1 2 k 2 + k 12 + k 23 m 2 δ 4 2 X ^ 2 2 + 1 2 k 3 + k 13 + k 23 m 3 δ 4 2 X ^ 3 2 + 1 2 2 k 12 m 1 m 2 X ^ 1 X ^ 2 + 2 k 13 m 1 m 3 X ^ 1 X ^ 3 + 2 k 23 m 2 m 3 X ^ 2 X ^ 3 .
We can confirm that, through this transformation, the invariant operator has been simplified in a way that its P ^ i 2 terms are represented in terms of the unit mass. However, the coupling terms which involve X i X j still remain. Due to this, it is still difficult to investigate the basic quantum features of the system in a straightforward way relying on the invariant.
In the next transformation, we will eliminate the coupling terms by diagonalizing the invariant. To this end, we write the invariant in a matrix form by introducing vectors X = ( X 1 , X 2 , X 3 ) T and P = ( P 1 , P 2 , P 3 ) T such that
I ^ = 1 2 i , j = 1 3 P ^ i δ i j P ^ j + 1 2 i , j = 1 3 X ^ i k i j X ^ j ,
where k i j are ith row and jth column elements of the matrix
k = ϖ 1 2 1 2 K 12 1 2 K 13 1 2 K 12 ϖ 2 2 1 2 K 23 1 2 K 13 1 2 K 23 ϖ 3 2 ,
while the involved parameters are of the form
ϖ 1 = ( k 1 + k 12 + k 13 ) / m 1 δ / 4 2 ,
ϖ 2 = ( k 2 + k 12 + k 23 ) / m 2 δ / 4 2 ,
ϖ 3 = ( k 3 + k 13 + k 23 ) / m 3 δ / 4 2 ,
K 12 = 2 k 12 m 1 m 2 , K 13 = 2 k 13 m 1 m 3 , K 23 = 2 k 23 m 2 m 3 .
In the next section, we will diagonalize the resultant invariant operator, Equation (29), by eliminating the coupling terms.

3.3. Rotation Matrix and Diagonalization of Invariant Operator

The diagonalization of the invariant operator I ^ can be done by making use of an algebraic approach based on another unitary transformation that corresponds to a 3D rotation. The unitary operator Λ ^ which will be used to perform a 3D rotation is parameterized by three Euler angles ( ϕ , θ , φ ) .
We consider a rotation matrix of the form
R = R X 1 ( ϕ ) R X 2 ( θ ) R X 3 ( φ ) ,
where
R X 1 ( ϕ ) = 1 0 0 0 c ϕ s ϕ 0 s ϕ c ϕ , R X 2 ( θ ) = c θ 0 s θ 0 1 0 s θ 0 c θ , R X 3 ( φ ) = c φ s φ 0 s φ c φ 0 0 0 1 .
In Equation (36), we have abbreviated notations as c ζ , s ζ cos ζ , sin ζ for convenience, while ζ ϕ , θ , φ . Note that the rotation matrix R is a real ( 3 × 3 ) orthogonal matrix where its determinant is unity; this can be expressed by an analogous unitary operator (rotation operator) Λ ^ ( t ) in quantum mechanics, such that
R Λ ^ = exp ( i ϕ J ^ 3 ) exp ( i θ J ^ 2 ) exp ( i φ J ^ 3 ) ,
where the three operators, J ^ 1 , J ^ 2 , and J ^ 3 , are angular momentum generators which are calculated from the definition of the angular momentum J = X × P . Therefore, from the infinitesimal version of Equation (37), we derive the three basic matrices
J ^ 1 = 0 0 0 0 0 i 0 i 0 , J ^ 2 = 0 0 i 0 0 0 i 0 0 , J ^ 3 = 0 i 0 i 0 0 0 0 0 .
These can also be expressed by a single formula
( J ^ k ) i j = i ϵ i j k ,
where ϵ i j k is the Levi-Civita symbol, which is an antisymmetric tensor. Note that the three matrices J ^ i satisfy the commutation relation
[ J ^ i , J ^ j ] = i ϵ i j k J ^ k .
The switching relationship given in the above equation shows that the operators, J ^ 1 , J ^ 2 , and J ^ 3 , constitute the generators of a Lie algebra with the structure of constants i ϵ i j k . The associated Lie group is in fact the group of rotation SO(3), which displays the relation between angular momentum operators and the rotation.
We will now see how to diagonalize the invariant I ^ by using the matrix representation of the rotation operator. Let us start by writing the matrix k in terms of the new diagonal matrix
k = R diag Ω 1 2 , Ω 2 2 , Ω 3 2 R 1 .
Then, it is possible to verify the relation R 1 k R = D , where
D = diag Ω 1 2 , Ω 2 2 , Ω 3 2 .
The new frequencies Ω i 2 , which represent the eigenvalues of the matrix k , are expressed as [44]
Ω 1 2 = 1 3 ϖ 1 2 + ϖ 2 2 + ϖ 3 2 + Ω 2 3 2 cos α ,
Ω 2 2 = 1 3 ϖ 1 2 + ϖ 2 2 + ϖ 3 2 + Ω 2 3 2 cos α + 2 π 3 ,
Ω 3 2 = 1 3 ϖ 1 2 + ϖ 2 2 + ϖ 3 2 + Ω 2 3 2 cos α 2 π 3 ,
where
Ω 2 = ϖ 1 2 ϖ 2 2 2 + ϖ 1 2 ϖ 3 2 2 + ϖ 2 2 ϖ 3 2 2 + 3 K 12 2 + K 13 2 + K 23 2 1 / 2 ,
α = arccos A 2 B 3 / 2 ,
while
A = 3 ϖ 1 2 + ϖ 2 2 ϖ 1 2 + ϖ 3 2 ϖ 2 2 + ϖ 3 2 27 1 4 ϖ 1 2 K 23 2 + 1 4 ϖ 2 2 K 13 2 + 1 4 ϖ 3 2 K 12 2 + 2 ϖ 1 6 + ϖ 2 6 + ϖ 3 6 + 18 ϖ 1 2 ϖ 2 2 ϖ 3 2 + 3 8 K 12 K 13 K 23 + 9 ϖ 1 2 + ϖ 2 2 + ϖ 3 2 1 4 K 12 2 + 1 4 K 13 2 + 1 4 K 23 2 ,
B = 1 2 ϖ 1 2 ϖ 2 2 2 + ϖ 1 2 ϖ 3 2 2 + ϖ 2 2 ϖ 3 2 2 + 3 2 1 4 K 12 2 + 1 4 K 13 2 + 1 4 K 23 2 .
By taking into account the commutation relation [ P 2 , J ^ i ] = 0 , we can confirm that the expressions of the old conjugate momenta P ^ i 2 appeared in the invariant operator I ^ remain unchanged. However, the new expressions of the coordinates are given by
x 1 = X 1 c θ c ϕ X 2 s θ s φ + c θ c φ s ϕ + X 3 c θ s ϕ s φ s θ s φ ,
x 2 = X 1 s ϕ + X 2 c ϕ c φ X 3 c ϕ s φ ,
x 3 = X 1 c ϕ s θ + X 2 c θ s φ s θ c φ s ϕ + X 3 c θ c φ + s θ s φ s ϕ .
p i also take similar expressions. From an algebraic evaluation considering these relations, we can show that I ^ takes the form
I ^ = 1 2 i = 1 3 p ^ i 2 + Ω i 2 x ^ i 2 .
Hence, the invariant I ^ corresponds to a sum of three simple harmonic oscillators with unit masses and constant frequencies. Although this simplified invariant has been obtained under the assumption given in Equation (5), it may also be possible to find a simple form of the invariant operator (as in Equation (53)) for the case where the damping coefficients of the three coupled oscillators are not the same each other. Thanks to this simple formula of the transformed invariant, the theory of invariant that we have developed is substantially useful in analyzing the mechanical properties of the system. In subsequent sections, we will show that quantum mechanical analysis of the system is possible by taking advantage of the simplified invariant operator, Equation (53).

3.4. Eigenfunctions of the Invariant Operator

As previously mentioned, the invariant quantity is helpful in understanding the quantum dynamics of the system as well as the classical dynamics. To analyze quantum mechanical characteristics of the system, let us see the eigenvalues and eigenfunctions of the invariant operator. In order to obtain them, we introduce creation and annihilation operators in the diagonalized system such that
b ^ i = Ω i / 2 x ^ i + i 2 Ω i p ^ i ,
b ^ i = Ω i / 2 x ^ i i 2 Ω i p ^ i .
We easily confirm that these operators obey the usual properties of ladder operators, including the boson canonical commutation rule [ b ^ i , b ^ i ] = 1 . We can also express Equation (53) in terms of b ^ i and b ^ i as
I ^ = i = 1 3 Ω i b ^ i b ^ i + 1 2 .
Now, let us turn our attention to the invariant operator I ^ ( t ) in the original system. The formula of I ^ ( t ) can be obtained from the inverse unitary transformation, which is
I ^ ( t ) = U ^ ( t ) I ^ U ^ 1 ( t ) .
We confirm that the canonical variables are changed through this transformation in a way that
X ^ i U ^ ( t ) X ^ i U ^ 1 ( t ) = m i e δ t 1 / 2 X ^ ,
P ^ i U ^ ( t ) P ^ i U ^ 1 ( t ) = m i e δ t 1 / 2 P ^ i + δ 2 m i e δ t X ^ i .
Regarding this, the invariant operator I ^ ( t ) , Equation (24), can be readily evaluated to be
I ^ ( t ) = i = 1 3 Ω i a ^ i ( t ) a ^ i ( t ) + 1 2 ,
where a ^ i ( t ) and a ^ i ( t ) are time-dependent canonical annihilation and creation operators that are defined as
a ^ i ( t ) = U ^ ( t ) b ^ i U ^ 1 ( t ) ,
a ^ i ( t ) = U ^ ( t ) b ^ i U ^ 1 ( t ) .
It is possible to obtain the full expressions of a ^ i ( t ) and a ^ i ( t ) using the relations in Equations (58) and (59): we have provided them in Appendix A for convenience.
According to conventional quantum mechanics, the number operators are given by a ^ i a ^ i . The eigenvalue equations for them can be written in the form
a ^ i a ^ i n i , t = n i n i , t ,
where the eigenvalues, n i , are three integer numbers. Because Equation (60) is represented in terms of a ^ i a ^ i , the eigenstates of the invariant operator are the same as those of a ^ i a ^ i , which are n i , t . Therefore, its normalized eigenstates can be written as a product of the three eigenfunctions of which formulae are derived from the independent zero-point states:
n 1 , n 2 , n 3 , t = n 1 , t n 2 , t n 3 , t = 1 n 1 ! n 2 ! n 3 ! ( a ^ 1 ) n 1 ( a ^ 2 ) n 2 ( a ^ 3 ) n 3 0 , 0 , 0 , t ,
whereas the corresponding eigenvalues are given by
λ n 1 , n 2 , n 3 = i = 1 3 Ω i n i + 1 2 .
On account of the normalization condition, the states in Equation (64) obey
n 1 , n 2 , n 3 , t n 1 , n 2 , n 3 , t = δ n 1 , n 1 δ n 2 , n 2 δ n 3 , n 3 .
In the configuration space, the normalized eigenstates are easily obtained by solving the eigenvalue equation of I ^ and they are given by
X 1 , X 2 , X 3 n 1 , n 2 , n 3 , t = i = 1 3 Ω i m i e δ t / 2 π 1 / 2 n i ! 2 n i 1 / 2 H n i ( Y i ) × exp i = 1 3 Ω i 2 + i δ 4 j = 1 3 R i j m j e δ t / 2 X j 2 ,
where H n are nth order Hermite polynomials, R i j are elements of R that correspond to ith row and jth column, and the functions Y i are represented as
Y 1 = Ω 1 e δ t 1 / 2 [ m 1 X 1 c θ c ϕ m 2 X 2 s θ s φ + c θ c φ s ϕ m 3 X 3 s θ c φ c θ s ϕ s φ ] ,
Y 2 = Ω 2 e δ t 1 / 2 m 1 X 1 s ϕ + m 2 X 2 c ϕ c φ m 3 X 3 c ϕ s φ ,
Y 3 = Ω 3 e δ t 1 / 2 [ m 1 X 1 c ϕ s θ + m 2 X 2 c θ s φ s θ c φ s ϕ + m 3 X 3 c θ c φ + s θ s φ s ϕ ] .
In the next section, we will derive Schrödinger solutions of the system by taking advantage of the eigenstates of I ^ , which we have obtained here.

3.5. The Schrödinger Equation and Its Solutions

The invariant operator and its eigenstates Equation (67) are useful in analyzing the dynamical properties of the system. We will derive the solutions of the Schrödinger equation of the system using their close relationship with the eigenstates of the invariant operator. If we write the time-dependent Schrödinger equation as
i t ψ n 1 , n 2 , n 3 ( t ) = H ^ ψ n 1 , n 2 , n 3 ( t ) ,
its solutions (wave functions) are written in the form
ψ n 1 , n 2 , n 3 ( t ) = e i ζ n 1 , n 2 , n 3 ( t ) n 1 , n 2 , n 3 , t ,
where ζ n 1 , n 2 , n 3 ( t ) are time-varying phases. From the substitution of Equation (72) into Equation (71), we see that the phases ζ n 1 , n 2 , n 3 ( t ) satisfy the equation
t ζ n 1 , n 2 , n 3 ( t ) = 1 n 1 , n 2 , n 3 , t t H ^ n 1 , n 2 , n 3 , t .
A minor evaluation for this equation leads to
ζ n 1 , n 2 , n 3 ( t ) = i = 1 3 Ω i n i + 1 2 t .
Finally, the wave functions in the configuration space are expressed as
X 1 , X 2 , X 3 ψ n 1 , n 2 , n 3 ( t ) = X 1 , X 2 , X 3 n 1 , n 2 , n 3 , t × exp i i = 1 3 Ω i n i + 1 2 t ,
where the eigenstates X 1 , X 2 , X 3 n 1 , n 2 , n 3 , t are given in Equation (67). Although the wave functions, Equation (75), in the original system are somewhat complicated, they are complete. These wave functions are necessary in evaluating quantum mechanical expectation values of various observables such as position and momentum quadratures, quantum energy, etc. The propagator, Wigner function, and entanglement of the system can also be investigated on the basis of them.

4. Conclusions

The dynamical invariant and its diagonalization for dissipative three coupled oscillators were investigated. We also treated the application of the dynamical invariant on quantization of the system. We have constructed the classical invariant quantity at first. Then, it was extended to a quantum one, i.e., we have obtained a rigorous form of the quantum invariant operator using the Liouville-von Neumann equation.
By a unitary transformation, the quantum invariant was transformed to that of a simple coupled oscillatory system which has unit masses. The resultant invariant operator was diagonalized eventually by further transformation using a rotation matrix parameterized by the Euler angles. During the diagonalization procedure, the coupling terms X ^ i X ^ j was eliminated and, as a consequence, the invariant reduced to the form of three uncoupled oscillators. The simplicity of the diagonalized quantum invariant is very advantageous in utilizing it in the analysis of the dynamical properties of the system.
The eigenfunctions of the quantum invariant operator were derived in the Fock state by solving its eigenvalue equation. By taking advantage of such eigenfunctions, the wave functions satisfying the Schrödinger equation have been obtained as shown in Equation (75). Such wave functions can be used to characterizing the quantum properties of various coupled oscillatory systems such as nano-optomecanical systems from a fundamental level.
Diagonalizing the invariant operator may have the similar complexity as the diagonalization of the Hamiltonian of the system. Nevertheless, the reason why we are interested in the invariant operator and its diagonalization is that the quantum wave functions associated with a time-dependent Hamiltonian such as Equation (4) are described in terms of the eigenstates of the invariant operator. In our case, the wave functions given in Equation (75) are represented in terms of the eigenstates of I ^ given in Equation (67).
We note that the invariant formalism for dynamical systems of which Hamiltonians depend on time admits to obtaining exact classical and quantum solutions without resorting to variational techniques or other approximation manipulations. As a matter of fact, we did not use any approximation in the derivation of (quantum) solutions from the diagonalization of the invariant in this research. This the merit of the approach of complicated dynamical systems based on an invariant, which distinguishes it from other methods in this field.

Author Contributions

S.M. (Salim Medjber), S.M. (Salah Menouar) and J.R.C. co-wrote the paper. S.M. (Salah Menouar) carried out final check of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea(NRF) grant funded by the Korea government(MSIT) (No.: NRF-2021R1F1A1062849).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Annihilation and Creation Operators in the Original System

From straightforward evaluations of Equations (61) and (62) using the formulae of b ^ i ( t ) and b ^ i ( t ) given in Equations (54) and (55), we have the expressions of the annihilation and creation operators as
a ^ 1 ( t ) = Ω 1 e δ t 2 1 / 2 [ m 1 X 1 c θ c ϕ m 2 X 2 s θ s φ + c θ c φ s ϕ + m 3 X 3 c θ s ϕ s φ s θ s φ ] + i e δ t 2 Ω 1 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 c θ c ϕ 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 s θ s φ + c θ c φ s ϕ + 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c θ s ϕ s φ s θ s φ ] ,
a ^ 1 ( t ) = Ω 1 e δ t 2 1 / 2 [ m 1 X 1 c θ c ϕ m 2 X 2 s θ s φ + c θ c φ s ϕ + m 3 X 3 c θ s ϕ s φ s θ s φ ] i e δ t 2 Ω 1 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 c θ c ϕ 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 s θ s φ + c θ c φ s ϕ + 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c θ s ϕ s φ s θ s φ ] ,
a ^ 2 ( t ) = Ω 2 e δ t 2 1 / 2 m 1 X 1 s ϕ + m 2 X 2 c ϕ c φ m 3 X 3 c ϕ s φ + i e δ t 2 Ω 2 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 s ϕ + 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 c ϕ c φ 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c ϕ s φ ] ,
a ^ 2 ( t ) = Ω 2 e δ t 2 1 / 2 m 1 X 1 s ϕ + m 2 X 2 c ϕ c φ m 3 X 3 c ϕ s φ i e δ t 2 Ω 2 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 s ϕ + 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 c ϕ c φ 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c ϕ s φ ] ,
a ^ 3 ( t ) = Ω 3 e δ t 2 1 / 2 [ m 1 X 1 c ϕ s θ + m 2 X 2 c θ s φ s θ c φ s ϕ + m 3 X 3 c θ c φ + s θ s φ s ϕ ] + i e δ t 2 Ω 3 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 c ϕ s θ + 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 c θ s φ s θ c φ s ϕ + 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c θ c φ + s θ s φ s ϕ ] ,
a ^ 3 ( t ) = Ω 3 e δ t 2 1 / 2 [ m 1 X 1 c ϕ s θ + m 2 X 2 c θ s φ s θ c φ s ϕ + m 3 X 3 c θ c φ + s θ s φ s ϕ ] i e δ t 2 Ω 3 1 / 2 [ 1 m 1 P ^ 1 + δ 2 m 1 e δ t X ^ 1 c ϕ s θ + 1 m 2 P ^ 2 + δ 2 m 2 e δ t X ^ 2 c θ s φ s θ c φ s ϕ + 1 m 3 P ^ 3 + δ 2 m 3 e δ t X ^ 3 c θ c φ + s θ s φ s ϕ ] .

References

  1. Zhou, J.-Y.; Zhou, Y.-H.; Yin, X.-L.; Huang, J.-F.; Liao, J.-Q. Quantum entanglement maintained by virtual excitations in an ultrastrongly-coupled-oscillator system. Sci. Rep. 2020, 10, 12557. [Google Scholar] [CrossRef] [PubMed]
  2. Kao, J.-Y.; Chou, C.-H. Quantum entanglement in coupled harmonic oscillator systems: From micro to macro. New J. Phys. 2016, 18, 073001. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, S.; Zhu, A.-D.; Jin, Z.; Yeon, K.-H.; Um, C.-I. Quantum squeezed effect of a mesoscopic capacitance coupled circuit with mutual inductance. J. Korean Phys. Soc. 2003, 43, 773–777. [Google Scholar] [CrossRef]
  4. Yeo, I.; de Assis, P.-L.; Gloppe, A.; Dupont-Ferrier, E.; Verlot, P.; Malik, N.S.; Dupuy, E.; Claudon, J.; Gérard, J.-M.; Auffèves, A.; et al. Strain-mediated coupling in a quantum dot–mechanical oscillator hybrid system. Nat. Nanotechnol. 2014, 9, 106–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Treutlein, P. Optomechanics: A strained couple. Nat. Nanotechnol. 2014, 9, 99–100. [Google Scholar] [CrossRef] [PubMed]
  6. Wilson, C.J.; Callaway, J.C. Coupled oscillator model of the dopaminergic neuron of the substantia nigra. J. Neurophysiol. 2000, 83, 3084–3100. [Google Scholar] [CrossRef] [PubMed]
  7. Frank, T.D.; Daffertshofer, A.; Peper, C.E.; Beek, P.J.; Haken, H. Towards a comprehensive theory of brain activity: Coupled oscillator systems under external forces. Physica D 2000, 144, 62–86. [Google Scholar] [CrossRef]
  8. Merdaci, A.; Jellal, A. Entanglement in three coupled harmonic oscillators. Phys. Lett. A 2019, 384, 126134. [Google Scholar] [CrossRef] [Green Version]
  9. Ndikilar, C.E.; Saleh, S.I.; Hafeez, H.Y.; Taura, L.S. Analytical calculation of power flow in three coupled oscillators and its application to one dimensional crystal. Eur. J. Appl. Phys. 2020, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]
  10. Deymier, P.; Runge, K. One-dimensional mass-spring chains supporting elastic waves with non-conventional topology. Crystals 2016, 6, 44. [Google Scholar] [CrossRef] [Green Version]
  11. Park, D. Dynamics of entanglement in three coupled harmonic oscillator system with arbitrary time-dependent frequency and coupling constants. Quantum Inf. Process. 2019, 18, 282. [Google Scholar] [CrossRef] [Green Version]
  12. Gonzalez-Henao, J.C.; Pugliese, E.; Euzzor, S.; Meucci, R.; Roversi, J.A.; Arecchi, F.T. Control of entanglement dynamics in a system of three coupled quantum oscillators. Sci. Rep. 2017, 7, 9957. [Google Scholar] [CrossRef] [Green Version]
  13. Agarwal, G.S.; Huang, S.M. Electromagnetically induced transparency in mechanical effects of light. Phys. Rev. A 2010, 81, 041803. [Google Scholar] [CrossRef] [Green Version]
  14. Naik, A.; Buu, O.; LaHaye, M.D.; Armour, A.D.; Clerk, A.A.; Blencowe, M.; Schwab, K.C. Cooling a nanomechanical resonator with quantum back-action. Nature 2006, 443, 193–196. [Google Scholar] [CrossRef] [Green Version]
  15. Chan, J.; Alegre, T.P.M.; Safavi-Naeini, A.H.; Hill, J.T.; Krause, A.; Groblacher, S.; Asp, M. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 2011, 478, 89–92. [Google Scholar] [CrossRef] [Green Version]
  16. Lin, Q.; Rosenberg, J.; Chang, D.; Camacho, R.; Eichenfield, M.; Vahala, K.J.; Painter, O. Coherent mixing of mechanical excitations in nano-optomechanical structures. Nat. Photon. 2010, 4, 236–242. [Google Scholar] [CrossRef] [Green Version]
  17. Shkarin, A.B.; Flowers-Jacobs, N.E.; Hoch, S.W.; Kashkanova, A.D.; Deutsch, C.; Reichel, J.; Harris, J.G.E. Optically mediated hybridization between two mechanical modes. Phys. Rev. Lett. 2014, 112, 013602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Hossein-Zadeh, M.; Vahala, K.J. Observation of optical spring effect in a microtoroidal optomechanical resonator. Opt. Lett. 2007, 32, 1611–1613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Kaushal, R.S.; Gupta, S. Construction of exact dynamical invariants in coupled oscillator problems. J. Phys. A Math. Gen. 2001, 34, 9879–9892. [Google Scholar] [CrossRef]
  20. Abdalla, M.S.; Leach, P.G.L. Wigner functions for time-dependent coupled linear oscillators via linear and quadratic invariant processes. J. Phys. A Math. Gen. 2005, 38, 881–893. [Google Scholar] [CrossRef]
  21. Korsch, H.J. Dynamical invariants and time-dependent harmonic systems. Phys. Lett. A 1979, 74, 294–296. [Google Scholar] [CrossRef]
  22. Kaushal, R.S.; Korsch, H.J. Dynamical Noether invariants for time-dependent nonlinear systems. J. Math. Phys. 1981, 22, 1904–1908. [Google Scholar] [CrossRef]
  23. Lutzky, M. Symmetry groups and conserved quantities for the harmonic oscillator. J. Phys. A Math. Gen. 1978, 11, 249–258. [Google Scholar] [CrossRef]
  24. Lutzky, M. Noether’s theorem and the time-dependent harmonic oscillator. Phys. Lett. A 1978, 68, 3–4. [Google Scholar] [CrossRef]
  25. Bertin, M.C.; Pimentel, B.M.; Ramirez, J.A. Construction of time-dependent dynamical invariants: A new approach. J. Math. Phys. 2012, 53, 042104. [Google Scholar] [CrossRef] [Green Version]
  26. Lewis, H.R., Jr.; Riesenfeld, W.B. An exact quantum theory of the time-dependent harmonic oscillator and of a charged particle in a time-dependent electromagnetic field. J. Math. Phys. 1969, 10, 1458–1473. [Google Scholar] [CrossRef]
  27. Lewis, H.R., Jr.; Riesenfeld, W.B. Class of exact invariants for classical and quantum time-dependent harmonic oscillators. J. Math. Phys. 1968, 9, 1976–1986. [Google Scholar] [CrossRef]
  28. Malkin, I.A.; Man’ko, V.I.; Trifonov, D.A. Invariants and evolution of coherent states for charged particle in time-dependent magnetic field. Phys. Lett. A 1969, 30, 414–416. [Google Scholar] [CrossRef]
  29. Malkin, I.A.; Man’ko, V.I.; Trifonov, D.A. Linear adiabatic invariants and coherent states. J. Math. Phys. 1973, 14, 576–582. [Google Scholar] [CrossRef]
  30. Dodonov, V.V.; Malkin, I.A.; Man’ko, V.I. Integrals of the motion, Green functions and coherent states of dynamical systems. Int. J. Theor. Phys. 1975, 14, 37–54. [Google Scholar] [CrossRef]
  31. De Castro, A.S.M.; Peruzzo, J.F.; Dodonov, V.V. Quantum state exchange between indirectly coupled modes. Phys. Rev. A 2005, 71, 032319. [Google Scholar] [CrossRef]
  32. Fan, H.-Y. Squeezing in the triatomic linear molecule model revealed by virtue of IWOP technique. J. Phys. A 1993, 26, 151–158. [Google Scholar]
  33. El-Orany, F.A.A.; Perina, J.; Abdalla, M.S. Statistical properties of three quantized interacting oscillators. Phys. Scr. 2001, 63, 128–140. [Google Scholar] [CrossRef]
  34. Caldirola, P. Porze non conservative nella meccanica quantistica. Nuovo Cimento 1941, 18, 393–400. [Google Scholar] [CrossRef]
  35. Kanai, E. On the quantization of dissipative systems. Prog. Theor. Phys. 1948, 3, 440–442. [Google Scholar] [CrossRef]
  36. Choi, J.R.; Menouar, S. Quantum dynamics of nano-optomechanical three coupled oscillators. 2021; unpublished. [Google Scholar]
  37. Qiao, Y.-F.; Li, R.-J.; Zhao, S.-H. Symmetry and invariant in generalized mechanical systems in the high- dimensional extended phase space. Acta Phys. Sin. 2001, 50, 814–815. [Google Scholar]
  38. Mei, F.X.; Wu, H.B.; Zhang, Y.F. Symmetries and conserved quantities of constrained mechanical systems. Int. J. Dyn. Control 2014, 2, 285–303. [Google Scholar] [CrossRef] [Green Version]
  39. Mei, F.X. Lie symmetries and conserved quantities of constrained mechanical systems. Acta Mech. 2000, 141, 135–148. [Google Scholar] [CrossRef]
  40. Peng, W. Conformal invariance and conserved quantities of mechanical system with unilateral constraints. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 463–471. [Google Scholar]
  41. Bateman, H. On dissipative systems and related variational principles. Phys. Rev. 1931, 38, 815–819. [Google Scholar] [CrossRef]
  42. Lemos, N.A. Canonical approach to the damped harmonic oscillator. Am. J. Phys. 1979, 47, 857–858. [Google Scholar] [CrossRef]
  43. McDonald, K. A Damped Oscillator as a Hamiltonian System. 2015. Preprint. Available online: http://kirkmcd.princeton.edu/examples/damped.pdf (accessed on 16 May 2021).
  44. Kronenburg, M.J. A method for fast diagonalization of a 2 × 2 or 3 × 3 real symmetric matrix. arXiv 2013, arXiv:1306.6291. [Google Scholar]
Figure 1. Schematic of a kind of nano-optomechanical three coupled oscillators on which our work can be applied potentially. ω 1 and ω 2 are frequencies associated with optical modes in cavities, whereas ω 3 is a frequency of a mechanical oscillation for a flexible nano membrane. This figure shows the interaction of cavity fields with the nano-mechanical oscillator via the force of radiation pressure.
Figure 1. Schematic of a kind of nano-optomechanical three coupled oscillators on which our work can be applied potentially. ω 1 and ω 2 are frequencies associated with optical modes in cavities, whereas ω 3 is a frequency of a mechanical oscillation for a flexible nano membrane. This figure shows the interaction of cavity fields with the nano-mechanical oscillator via the force of radiation pressure.
Entropy 23 00837 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Medjber, S.; Menouar, S.; Choi, J.R. Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators. Entropy 2021, 23, 837. https://doi.org/10.3390/e23070837

AMA Style

Medjber S, Menouar S, Choi JR. Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators. Entropy. 2021; 23(7):837. https://doi.org/10.3390/e23070837

Chicago/Turabian Style

Medjber, Salim, Salah Menouar, and Jeong Ryeol Choi. 2021. "Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators" Entropy 23, no. 7: 837. https://doi.org/10.3390/e23070837

APA Style

Medjber, S., Menouar, S., & Choi, J. R. (2021). Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators. Entropy, 23(7), 837. https://doi.org/10.3390/e23070837

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop