Received Date: June 06, 2016; Accepted Date: September 23, 2016; Published Date: November 04, 2016
Citation: Kenfack SC, Nguimeya GP, Talla PK, Fotue AJ et al. (2016) Decoherence of Driven Coupled Harmonic Oscillator. J Nanosci Curr Res 1: 104.
Copyright: © 2016 Kenfack SC, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Nanosciences: Current Research
We consider a pair of linearly coupled harmonic oscillators to explore the decoherence phenomenon induces by interaction of a quantum system with a classical environment using Feynman path’s integral method. We determine the DCHO propagator afterwards thermodynamics parameters associated to the system. We show numerically that one can reduce decoherence of a system by coupling this one to a driven harmonic oscillator or by carrying this system to the resonance.
Decoherence; Propagator; Coupled oscillator
The study of the phenomenon of decoherence has attracted increasing attention in recent years . In fact, it has been recognized that decoherence is of fundamental importance in the comprehension of the nature of the boundary between the quantum and classical worlds. The nature of this boundary has been meticulously examined from both the theoretical and experimental points of view [2,3]. Decoherence is closely related to one of the greatest mysteries of physics which up to now remains unanswered, that is, the quantum measurement problem. This problem was clearly formulated by Erwin Schrodinger in 1935 in his experiment of thought known as “Schrodinger’s cat paradox” [4-6] pointing out the interaction problem between quantum systems where quantum effects dominate and classical systems which operate in the classical limit. Two linearly coupled oscillators provide an ideal test model for exploring this aspect of the quantum-classical limit interaction . In this work, we use the Feynman path integral formalism which enables an evaluation of the sum over all possible paths in a consistent way. This formulation which appears to be appropriate for a correct description of the reality is exact, equivalent to the Hilbert space formalism and leads to the same results for the Schrödinger equation. It is an alternative formulation of the quantum mechanics, which introduces only functions with numerical values and does completely without operators. This paper is organized as follows: In section 2, we present the model and determine the coupled harmonic oscillator (CHO) propagator from which we deduce the driven coupled harmonic oscillator(DCHO) propagator. In section 3, we derive the thermodynamic parameters associated to the system which is entropy. In section 4, we present the numerical results and then, end with the conclusion in section 5.
The DCHO propagator
The DCHO is a system constituted of two harmonic oscillators coupled by a spring and whose extremities are subjected to external driving forces represented by and . Here, we assume that the two oscillators have equal mass m and coupling constant k .
The CHO propagator
The determination of this propagator is of great interest because it will allow us later by comparison to establish the explicit form of the DCHO propagator. The simplest version of a pair of coupled quantum mechanical simple harmonic oscillators with linear coupling consists of two identical oscillators, with equal masses, spring constants, and frequencies, plus a connecting spring with its own spring constant. Considering the system described above, the Hamiltonian reads:
where x1, x2, p1, p2 are respectively the displacements and momenta of each oscillator, and ? their oscillation frequency.
From the formula , the Lagrangian of the system reads
The differentials equations that govern the evolution of the system are given by
When setting f1(t)=f2(t)=0, the system of Equation(3) becomes
The Lagrangian is transform to
Then, the action is
by making the following change of variables,
the action splits into two terms
is the classical action and yi the deviation between xi(t) and theirs classical path xci which is given by yi = xi - xci,(i = 1,2).
The propagator is thus transformed as follows:
is a smooth function in the limit
the normalization factor F(t) of the system becomes
This factor can be written as the product of the normalization factor of each oscillator as
Thus, it yields that 
Therefore, the propagator takes the form
We now focus on the classical action. Performing the integration of Equation(11) and taking into account Equation(4), one finds 
To obtain the exact form of this action, we have to solve the system of Equation(4). By setting
This system becomes
of which the solutions are
When returning to initials variables, one gets
thus, we have
From the system of Equation(30) we obtain
and by performing the computations x1+x2 and x1-x2, we get respectively
By using Eqs. (27) - (33) one finds from Equation(23) that 
substituting this result into Equation(22), we obtain the CHO propagator which is 
Deducing of the DCHO propagator: Setting f1(t) = f2(t) = 0, the DCHO is reducing to CHO. Since the Lagrangian is quadratic, then we can assume the propagator for the DCHO to have the form 
Where a(t),b(t),c(t),d(t),g(t),h(t) are the functions of time which have to be determined.
Taking and the Hamiltonian of Equation(1) takes the form
The Hamiltonian of Equation(37) satisfies the Schrödinger equation 
From Equation(37), when using Equation (35), we get
From the two first relations of above system, we get, taking this relation into account
we have from Equation(41) the following system
When solving each equation of the system of Equation(43), one obtains
From Equation (43) and using Equation (41), we get
One can see that Equation (44) does not include the excitation forces f1(t) and f2(t). Thus setting f1(t) = f2(t) = 0 the functions a(t) and b(t) can be equalized with the coefficients of x12 and x22 of equation (34). The comparison between these two equations shows that ?1 and ?2 must be equal to zero. And so
By replacing a(t) = c(t) into the forth equation of system in Equation (39) and making the following change
Solving these two differentials equations gives
Where ?1 and ?2 are the constants to be determine. When taking into account Equation (48), then, Equation(46) gives 
When combining Eqs. (45) and (49), we find 
Where ϑ is a constant to be determined.
To determine the constants Λ1, Λ2 and ϑ , we make a comparison between Equations (34) and (35) for f1(t) = f2(t) = 0. From this comparison, it comes that
When inserting Equations (45) and (49)-(51) into Equation(35), gathering the similar terms and replacing by we work out the DCHO propagator which reads 
Thermodynamics Parameters for the DCHO
To determine thermodynamic parameters associated to the system, the knowledge of the time-dependent wave function is necessary as well as the wave function at initial time.
Wave function at initial time: The considered system here can be modeled as follows (Figure 1)
The potential energy of this system is:
Where both k and kc are the elastic constant and kc characterize the coupling of particles at point x`=0 with each other.
The kinetic energy of this system is written as:
When introduce the coordinate of the center of mass Y and the relative coordinate X.
Equations (53) and (54) are transformed to
Thus, the Hamiltonian reads
From the Hamilton canonical equations, we get
Then the Schrodinger equation can be written as:
By separation of variables while taking E= E1 + E2 this equation yields to
The solutions of system in Equation(62) are looking for under the forms
In the fundamental state n1=0 and, n2=0, what leads to:
When returning to the initials variables, we obtain
with and given by Equations. (52) and (69) respectively.
As a reminder, and x1, x2 are the initial and final coordinates of the oscillators respectively. In a much more compact form easy to manipulate, the propagator in Equation (54) can be written as follows:
When changing I1, I2, I3 and I4 into
this propagator can finally be expressed as follows:
When putting Equation(77) into Equation(71) and doing some arrangement, we then get
When changing x1 and x2 into
then, Equation(78) becomes
Equation (81) is easily integrate using Gaussian integral, then, one gets the expression of the time-dependent wave function which is
Where we recall that
Probability density: In term of the time-dependent wave function, the probability density is define by
Carrying out this product using Gaussian integrals and taking into account Equations (83), we obtain after few arrangements the probability density of the system, which reads
Entropy: By definiton, the Gibbs-Shannon entropy reads
and proceeding to the changes
Integrating Equation(92) using Gaussian integrals yields to the entropy of the system which is
We notice in figure 2 that the entropy decreases with time and tends towards a limit.
Therefore, decoherence decrease (Figure 2).
This result is in accordance with that obtained by M.A. de Ponte et al. for two coupled dissipative oscillators. Physically, this means that the whole system is organizes itself and evolves towards a coherent state in which information will be preserved.
From figure 3, we note that the entropy of the system oscillates by keeping a constant amplitude which means that the exchange between the system and its surrounding are quasi-perfect so that the system keeps a certain coherence. The system is slowly affected by external forces. Physically, this coherence can be interpreted as being due to the fact that information remains preserved in time. This result is in accordance with that obtained by Tabue et al.  for the pure states of a damped harmonic oscillator.
Figure 4 shows that entropy increases in time, and thus decoherence increases. This means that, the whole system collects environment which accentuates decoherence and so the system disorganizes itself. One could expect such a result which is in accordance with that obtained by Tabue et al.  for a damped harmonic oscillator.
Using the Feynman path integral method we have determined the DCHO propagator after which, considering the particular case where the external driving forces are sinusoidal, we have determined the entropy associated to the system. We reveal numerically that, when two harmonic oscillators subjected to an external driving force are coupled, one compensates the effects of the other so that the whole system organizes itself and evolves towards an equilibrium position in which it will conserve certain coherence. That is explained by the decrease of the entropy which thereafter tends towards a limit. We also reveal that at resonance the system is in a coherent state which is explained by the conservation, with time of the oscillation amplitudes of the entropy of the system. Thus, we have shown that to reduce decoherence in a system one can just couple such a system with a driven harmonic oscillator or the system in question to resonance.