Next Article in Journal
Entropy Interpretation of Hadamard Type Fractional Operators: Fractional Cumulative Entropy
Next Article in Special Issue
Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero
Previous Article in Journal
Universal Non-Extensive Statistical Physics Temporal Pattern of Major Subduction Zone Aftershock Sequences
Previous Article in Special Issue
Extreme Electron Acceleration with Fixed Radiation Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method

1
Faculdade de Física, Universidade Federal do Pará, Belém 66075-110, PA, Brazil
2
Centro de Física, Universidade do Minho, 4710-057 Braga, Portugal
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(12), 1851; https://doi.org/10.3390/e24121851
Submission received: 30 November 2022 / Revised: 14 December 2022 / Accepted: 16 December 2022 / Published: 19 December 2022
(This article belongs to the Special Issue Quantum Nonstationary Systems)

Abstract

:
Harmonic oscillators with multiple abrupt jumps in their frequencies have been investigated by several authors during the last decades. We investigate the dynamics of a quantum harmonic oscillator with initial frequency ω 0 , which undergoes a sudden jump to a frequency ω 1 and, after a certain time interval, suddenly returns to its initial frequency. Using the Lewis–Riesenfeld method of dynamical invariants, we present expressions for the mean energy value, the mean number of excitations, and the transition probabilities, considering the initial state different from the fundamental. We show that the mean energy of the oscillator, after the jumps, is equal or greater than the one before the jumps, even when ω 1 < ω 0 . We also show that, for particular values of the time interval between the jumps, the oscillator returns to the same initial state.

1. Introduction

The quantum harmonic oscillator potential with time-dependent parameters is relevant in modeling several problems in physics, and it has been investigated [1,2,3,4,5,6,7,8,9]. For example, the interaction between a spinless charged quantum particle and a time-dependent external classical electromagnetic field can be studied through a harmonic potential whose frequency depends explicitly on time [4,10,11,12,13], and this is used to model the quantum motion of this particle in a trap [14,15,16,17,18,19]. In the context of quantum electrodynamics, this potential is useful, for instance, to describe the free electromagnetic field in nonstationary media [8,9,20]. In the context of shortcuts to adiabaticity, time-dependent quantum oscillators have also been considered [21,22,23,24,25,26]. Other applications are found in relativistic quantum mechanics, quantum field theory, dynamical Casimir effect, and gravitation [27,28,29,30,31,32,33].
A particular case of a quantum harmonic oscillator with time-dependent parameters that shows sudden frequency jumps is investigated, for instance, in refs. [21,23,24,34,35,36,37,38,39,40,41]. Under such jumps (or any time dependence in the parameters), a classical oscillator in its ground state remains in the same state, whereas a quantum oscillator can become excited [35]. Moreover, the wave functions of quantum harmonic oscillators with time-dependent parameters describe squeezed states [5,9,42], For example, a sudden change in the oscillation frequency of 85 Rb atoms in the vibrational fundamental state of a one-dimensional optical lattice generates squeezed states [43]. The description of squeezed states is relevant, for instance, in the implementation of schemes for noise minimization in quantum sensors, which increases their sensitivity (see, for instance, ref. [44] and references therein). Subtle points involving the squeezed states for the model of two frequency jumps were investigated, for instance, by Tibaduiza et al. in ref. [40], where the solution for this case was obtained via algebraic method.
In the present paper, we investigate the dynamics of a quantum harmonic oscillator with initial frequency ω 0 that undergoes a sudden jump to a frequency ω 1 and, after a certain time interval, suddenly returns to its initial frequency. Instead of using the algebraic method used in ref. [40], here, we use the Lewis–Riesenfeld (LR) method of dynamical invariants. The LR method [2,3,4] enables the calculation of the exact wave function of a system subjected, for instance, to a harmonic oscillator potential with time-dependent parameters, such as mass and frequency [5,42,45]. Using this method, we show that the results for the squeeze parameters, the quantum fluctuations of the position and momentum operators, and the probability amplitude of a transition from the fundamental state to an arbitrary energy eigenstate coincide with those found in ref. [40]. In addition, using the same LR method, we also obtain expressions for the mean energy value and for the mean number of excitations (which were not calculated in ref. [40]), as well as for the transition probabilities considering the initial state different from the fundamental (which generalizes the formula found in ref. [40]).
The paper is organized as follows. In Section 2.1, we review some results of the application of the LR method to the quantum harmonic oscillator with time-dependent frequency. In Section 2.2, we define the squeezing parameters and, from these and the oscillator wave function obtained via the LR method, we determine the quantum fluctuations of the position, momentum and Hamiltonian operators, the mean number of excitations, and the transition probabilities between different states. In Section 3, we apply the results of previous sections to the model of ref. [40] and analyze their physical implications. In Section 4, we present our final remarks.

2. Analytical Method

2.1. The Wave Function of the Harmonic Oscillator via Lewis–Riesenfeld Method

Let us consider the one-dimensional Schrödinger equation for a system whose Hamiltonian H ^ ( t ) explicitly depends on time [46,47,48],
i Ψ ( x , t ) t = H ^ ( t ) Ψ ( x , t ) .
According to the LR method [2,3,4,5,9,42], given an invariant Hermitian operator I ^ ( t ) which satisfies
I ^ ( t ) t + 1 i I ^ ( t ) , H ^ ( t ) = 0 ,
a particular solution Ψ n ( x , t ) of Equation (1) is
Ψ n ( x , t ) = exp i α n ( t ) Φ n ( x , t ) ,
in which Φ n ( x , t ) are the eigenfunctions of I ^ ( t ) , found from
I ^ ( t ) Φ n ( x , t ) = λ n Φ n ( x , t ) ,
with λ n being time-independent eigenvalues of I ^ ( t ) and α n ( t ) phase functions, obtained from the equation
d α n ( t ) d t = + d x Φ n * ( x , t ) i t 1 H ^ ( t ) Φ n ( x , t ) .
The general solution Ψ ( x , t ) of Equation (1) is
Ψ ( x , t ) = n = 0 C n Ψ n ( x , t ) ,
where the time-independent coefficients C n depend only on the initial conditions.
Specifically, for a time-dependent one-dimensional harmonic oscillator with mass m 0 , whose time dependence is contained purely in its oscillation frequency ω ( t ) , the Hamiltonian is given by
H ^ ( t ) = p ^ 2 2 m 0 + 1 2 m 0 ω ( t ) 2 x ^ 2 ,
where x ^ and p ^ are position and momentum operators, respectively, with x ^ , p ^ = i . An operator I ^ ( t ) associated with Equation (7) is [4,5,9,42]
I ^ ( t ) = 1 2 x ^ ρ ( t ) 2 + ρ ( t ) p ^ m 0 ρ ˙ ( t ) x ^ 2 ,
wherein ρ ( t ) is a real parameter which is the solution of the Ermakov–Pinney equation [49,50,51,52]
ρ ¨ ( t ) + ω ( t ) 2 ρ ( t ) = 1 m 0 2 ρ ( t ) 3 .
The eigenfunctions of I ^ ( t ) , given by Equation (8), are
Φ n ( x , t ) = 1 2 n n ! Φ 0 ( x , t ) H n x 1 2 ρ ( t ) ,
where H n are the Hermite polynomials of order n [53], λ n = n + 1 / 2 , and
Φ 0 ( x , t ) = 1 π ρ ( t ) 2 1 4 exp i m 0 2 ρ ˙ ( t ) ρ ( t ) + i m 0 ρ ( t ) 2 x 2 .
From Equation (5), the functions α n ( t ) are given by
α n ( t ) = 1 m 0 n + 1 2 0 t d t ρ t 2 .
Thus, from Equations (3), (10) and (12), the wave function Ψ n ( x , t ) associated with the Hamiltonian (7) is
Ψ n ( x , t ) = 1 2 n n ! exp i m 0 n + 1 2 0 t d t ρ ( t ) 2 Φ 0 ( x , t ) H n x 1 2 ρ ( t ) .
For the case in which the frequency is always constant ( ω ( t ) = ω 0 ), the solution of Equation (9) is ρ ( t ) = ρ 0 , where [4,6,42]
ρ 0 = 1 m 0 ω 0 .
Therefore, Equation (13) falls back to the wave function of a harmonic oscillator with time-independent mass and frequency, Ψ n ( 0 ) ( x , t ) , given by [46,47,48]
Ψ n ( 0 ) ( x , t ) = 1 2 n n ! m 0 ω 0 π 1 4 exp i n + 1 2 ω 0 t m 0 ω 0 x 2 2 H n m 0 ω 0 1 2 x .

2.2. Squeeze Parameters, Quantum Fluctuations, Mean Number of Excitations, and Transition Probability

As discussed in refs. [5,9,42,54], the quantum states of the time-dependent oscillator, characterized by the wave function Ψ n ( x , t ) (Equation (13)), are squeezed. Thus, we can define the squeeze parameter r ( t ) and the squeeze phase ϕ ( t ) , which specify the squeezed state, in terms of the parameter ρ ( t ) [55]:
r ( t ) = cosh 1 m 0 2 ρ ˙ ( t ) 2 + ρ ( t ) 2 + 2 m 0 ω 0 + m 0 2 ω 0 2 ρ ( t ) 2 4 m 0 ω 0 1 2 ,
ϕ ( t ) = cos 1 1 + m 0 ω 0 ρ ( t ) 2 2 cosh 2 [ r ( t ) ] 2 sinh [ r ( t ) ] cosh [ r ( t ) ] ,
with r ( t ) 0 and 0 ϕ ( t ) 2 π . From Equation (13), one can also obtain the expected value of a given observable O ^ ( t ) in the state Ψ n ( x , t ) as
O ^ ( t ) ( n , t ) = + d x Ψ n * ( x , t ) O ^ ( t ) Ψ n ( x , t ) ,
which, from Equations (16) and (17), can be written in terms of r ( t ) and ϕ ( t ) . For the operators x ^ and p ^ , one has [56]:
x ^ ( n , t ) = p ^ ( n , t ) = 0 ,
x ^ 2 ( n , t ) = n + 1 2 m 0 ω 0 { cosh 2 r ( t ) + sinh 2 r ( t ) + 2 sinh r ( t ) × cosh r ( t ) cos ϕ ( t ) } ,
p ^ 2 ( n , t ) = n + 1 2 m 0 ω 0 { cosh 2 r ( t ) + sinh 2 r ( t ) 2 sinh r ( t ) × cosh r ( t ) cos ϕ ( t ) } ,
where it follows, from Equations (7), (20), and (21), that
H ^ ( t ) ( n , t ) = p ^ 2 ( n , t ) 2 m 0 + 1 2 m 0 ω ( t ) 2 x ^ 2 ( n , t ) .
From Equations (19)–(21), one finds the variances of the operators x ^ :
[ Δ x ^ ] 2 ( n , t ) = x ^ 2 ( n , t ) [ x ^ ( n , t ) ] 2 ,
and p ^ :
[ Δ p ^ ] 2 ( n , t ) = p ^ 2 ( n , t ) [ p ^ ( n , t ) ] 2 ,
which implies the uncertainty relationship
[ Δ x ^ ] 2 ( n , t ) [ Δ p ^ ] 2 ( n , t ) n + 1 2 2 2 { cosh 4 [ r ( t ) ] + sinh 4 [ r ( t ) ] 2 sinh 2 [ r ( t ) ] × cosh 2 [ r ( t ) ] cos [ 2 ϕ ( t ) ] } .
Due to the time dependence of the frequency, one can also determine the mean number of excitations N ^ ( n , t ) that a system, subjected to this potential, can undergo. This is given by [57,58,59]
N ^ ( n , t ) = n + 2 n + 1 sinh 2 r ( t ) .
For the fundamental state n = 0 , one finds N ^ ( 0 , t ) = sinh 2 r ( t ) , a result that agrees with refs. [33,56] for vacuum squeezed states. This means that a system, even in the fundamental state, could be excited due to the temporal variations in its frequency. The system subjected to the time-dependent harmonic potential can also make transitions between different states, since time-dependent potentials induce quantum systems to make transitions [46,47,48]. Let us consider that the system is initially at a stationary state Ψ m ( 0 ) ( x , t = 0 ) with frequency ω 0 and, due to a modification in its frequency from ω 0 to ω ( t ) , it evolves to a new state Ψ m ( x , t ) (Equation (13)). In this way, the probability to find the system in the state Ψ n ( 0 ) ( x , t ) (Equation (15)), is given by [46,47,48]
P ( t ) m n = + d x Ψ n * ( 0 ) ( x , t ) Ψ m ( x , t ) 2 .
Using Equations (13) and (15), one can find that P ( t ) m n = 0 for odd values of | n m | , and [57,60]
P ( t ) m n = 2 m + n min ( m , n ) ! 2 sinh [ r ( t ) ] | n m | m ! n ! cosh [ r ( t ) ] × k = | n m | 2 n + m 2 n + m 2 k n + m + 2 k 2 4 n + m 2 k ! k | n m | 2 ! cosh k [ r ( t ) ] 2 ,
for even values of | n m | , where min ( m , n ) is the smallest value between m and n. Note that the fact that Equation (28) is nonzero only for even values of | n m | is related to the parity of the harmonic potential [61]. From Equations (26) (making n = 0 ) and (28), we can relate the probability P ( t ) m n to the mean number of excitations in the fundamental state N ^ ( 0 , t ) :
P ( t ) m n = 2 m + n min ( m , n ) ! 2 N ^ ( 0 , t ) | n m | 2 m ! n ! N ^ ( 0 , t ) + 1 1 2 × k = | n m | 2 n + m 2 n + m 2 k n + m + 2 k 2 4 n + m 2 k ! k | n m | 2 ! N ^ ( 0 , t ) + 1 k 2 2 .
It follows that if the mean number of excitations in the fundamental state is nonzero, then the oscillator will have nonzero probabilities of making transitions between different energy levels. When we consider m = n = 0 in Equation (28) (or in Equation (29)), one has the probability of persistence in the fundamental state, and, as a consequence, one can also obtain the probability of excitation, given by 1 P ( t ) 0 0 [40].

3. Oscillator with Two Frequency Jumps

Now, we apply the formulas shown in Section 2 to investigate the model discussed in ref. [40], namely an oscillator with
ω ( t ) = ω 0 , t < 0 , ω 1 , 0 < t < τ , ω 0 , t > τ ,
in which ω 0 and ω 1 are constant frequencies, and τ is the length of the time interval between the frequency jumps.

3.1. Solution and General Behavior of the ρ ( t ) Parameter

Due to the form of Equation (30), the ρ ( t ) parameter can be written as
ρ ( t ) = ρ 0 , t < 0 , ρ 1 ( t ) , 0 < t < τ , ρ 2 ( t ) , t > τ ,
where ρ 0 is given in Equation (14), and ρ 1 ( t ) and ρ 2 ( t ) are calculated next.

3.1.1. Interval 0 < t < τ

For the interval 0 < t < τ , the equation to be solved is
ρ ¨ 1 ( t ) + ω 1 2 ρ 1 ( t ) = 1 m 0 2 ρ 1 ( t ) 3 ,
with the conditions [6]
ρ 1 t = 0 = 1 m 0 ω 0 , ρ ˙ 1 t = 0 = 0 .
The general solution for ρ 1 ( t ) is of the form [50,51]
ρ 1 ( t ) = A 1 sin 2 ω 1 t + B 1 cos 2 ω 1 t + 2 C 1 sin ω 1 t cos ω 1 t 1 2 ,
and the relationship between the constants A 1 , B 1 , and C 1 is
A 1 B 1 C 1 2 = 1 m 0 2 ω 1 2 .
Then, applying conditions (33) to Equation (34) and using relation (35), we obtain
ρ 1 ( t ) = ω 0 sin 2 ω 1 t m 0 ω 1 2 + cos 2 ω 1 t m 0 ω 0 1 2 .

3.1.2. Interval t > τ

In the interval t > τ , the Ermakov–Pinney equation has the form
ρ ¨ 2 ( t ) + ω 0 2 ρ 2 ( t ) = 1 m 0 2 ρ 2 ( t ) 3 .
The general solution of Equation (37) is
ρ 2 ( t ) = A 2 sin 2 ω 0 t + B 2 cos 2 ω 0 t + 2 C 2 sin ω 0 t cos ω 0 t 1 2 ,
with the constants A 2 , B 2 , and C 2 determined from the relationship
A 2 B 2 C 2 2 = 1 m 0 2 ω 0 2 ,
and the conditions for continuity
ρ 1 ( t = τ ) = ρ 2 ( t = τ ) , ρ ˙ 1 ( t = τ ) = ρ ˙ 2 ( t = τ ) .
Using Equations (36), (38), (39), and (40) results in
A 2 = 1 m 0 ω 0 3 ω 1 2 { ω 0 2 ω 1 2 + ( ω 0 4 ω 1 4 ) sin 2 ( ω 0 τ ) ω 0 2 ω 1 2 + ω 1 4 sin 2 ( ω 1 τ ) + 2 ω 0 ω 1 ( ω 0 ω 1 ) ( ω 0 + ω 1 ) sin ( ω 0 τ ) cos ( ω 0 τ ) sin ( ω 1 τ ) cos ( ω 1 τ ) } ,
B 2 = 1 m 0 ω 0 3 ω 1 2 { ω 0 2 ω 1 2 + ( ω 1 4 ω 0 4 ) sin 2 ω 0 τ ω 0 2 ω 1 2 + ω 0 4 sin 2 ω 1 τ 2 ω 0 ω 1 ω 0 ω 1 ω 0 + ω 1 sin ω 0 τ cos ω 0 τ sin ω 1 τ cos ω 1 τ } ,
C 2 = 1 m 0 ω 0 3 ω 1 2 { [ ( ω 0 2 + ω 1 2 ) sin ω 0 τ cos ω 0 τ sin ω 1 τ 2 ω 0 ω 1 ( sin 2 ω 0 τ 1 / 2 ) cos ω 1 τ ] ω 0 ω 1 ω 0 + ω 1 sin ω 1 τ } .

3.1.3. General Behavior

The general solution for the ρ ( t ) parameter is given by Equation (31), with ρ 0 , ρ 1 ( t ) and ρ 2 ( t ) given by Equations (14), (36), (38), (41), (42), and (43). From these equations, it can be seen that the ρ ( t ) parameter is a periodic function of time. Moreover, even when the frequency returns to its initial value ω 0 , this parameter will still, in general, be a periodic function of time. However, if we define τ u = u π / ω 1 ( u > 0 ), and make τ = τ l , where l N , the ρ ( t ) parameter returns to ρ 0 (Equation (14)), which means that although the oscillator feels the effect of the change in its frequency when it jumps from ω 0 to ω 1 , if the frequency returns to ω 0 at τ = τ l , for t > τ l , the oscillator behaves as if nothing happened. In other words, if τ = τ l , the abrupt change in the frequency is imperceptible to the oscillator when t > τ l . On the other hand, when τ = τ l + 1 / 2 , the ρ ( t ) parameter reaches its maximum value. The behavior of ρ ( t ) is shown in Figure 1.

3.2. Squeeze Parameters

Because the frequency varies abruptly, squeezing occurs in the system [35,36,40]. Thus, now, we calculate the parameters r ( t ) and ϕ ( t ) associated with the model in Equation (30). We show that our results for these parameters agree with those found in ref. [40] via an exact algebraic method.

3.2.1. Parameter r ( t )

The parameter r ( t ) for any time interval is given by
r ( t ) = 0 , t < 0 , r 1 ( t ) , 0 < t < τ , r 2 ( t ) , t > τ .
Note that r ( t < 0 ) = 0 because the frequency of the oscillator is time-independent in this interval. Using Equation (36) in Equation (16), we obtain, for the interval 0 < t < τ , the squeezing parameter r 1 ( t ) , where
r 1 ( t ) = cosh 1 1 + ω 1 2 ω 0 2 2 ω 0 ω 1 2 sin 2 ( ω 1 t ) ,
which is a result that agrees with the one found in ref. [40]. For the interval t > τ , using Equation (38) in Equation (16), we find that r 2 ( t ) , with
r 2 ( t ) = r 1 ( τ ) ,
this also agrees with ref. [40]. Note that for τ = τ l , one has r 2 ( t ) = 0 . In Figure 2 (also found in ref. [40]), one can see the behavior of r ( t ) for some values of τ .

3.2.2. Parameter ϕ ( t )

The squeeze phase ϕ ( t ) for any time interval has the form
ϕ ( t ) = undefined , t < 0 , ϕ 1 ( t ) , 0 < t < τ , ϕ 2 ( t ) , t > τ ,
where ϕ ( t ) for t < 0 is undefined because there is no squeeze in this interval. Using Equations (36) and (45) in Equation (17), we obtain that the squeeze phase for the interval 0 < t < τ , ϕ 1 ( t ) is given by
ϕ 1 ( t ) = cos 1 ω 0 4 ω 1 4 sin 2 ( ω 1 t ) 4 ω 0 2 ω 1 2 + ω 1 2 ω 0 2 2 sin 2 ( ω 1 t ) ω 1 2 ω 0 2 2 sin 2 ( ω 1 t ) ,
which also agrees with ref. [40]. To calculate the squeeze phase in the interval t > τ , the reasoning is analogous, simply substituting Equations (38) and (46) into Equation (17).
From Figure 3, it can be seen that the squeezing phase will continue to vary in time, even for t > τ .
Due of this time dependence, the fluctuations of the x ^ and p ^ operators will continue to depend on time in this interval, as we see later in Section 3.3 (see Equations (23) and (24)). Another point to be observed in Figure 3 concerns the behavior of the squeezing phase in the interval t > τ , when τ = τ l (in the specific case of Figure 3, τ l = π ). Since the squeeze parameter r 2 ( t ) (Equation (46)) is zero in this case, the system is no longer squeezed. Consequently, the squeeze phase is undefined for τ l . Therefore, its effect on the system will be negligible because there will be no more squeezing.

3.3. Quantum Fluctuations

The variance of the x ^ operator for any time interval is given by
[ Δ x ^ ] 2 ( n , t ) = [ Δ x ^ ] 0 2 ( n ) , t < 0 , [ Δ x ^ ] 1 2 ( n , t ) , 0 < t < τ , [ Δ x ^ ] 2 2 ( n , t ) , t > τ ,
where [46]
[ Δ x ^ ] 0 2 ( n ) = n + 1 2 m 0 ω 0 .
Substituting Equations (45) and (48) into (23), we find, for the x ^ operator in the interval 0 < t < τ ,
[ Δ x ^ ] 1 2 ( n , t ) = ω 0 2 ω 1 2 sin 2 ( ω 1 t ) + cos 2 ( ω 1 t ) [ Δ x ^ ] 0 2 ( n ) ,
which is in agreement with refs. [35,40]. For the interval t > τ , the procedure is analogous, using Equations (16), (17), and (38) in Equation (23). In Figure 4, we show the behavior of [ Δ x ^ ] 2 ( n , t ) for n = 0 .
Similarly, for the variance of the p ^ operator, we have
[ Δ p ^ ] 2 ( n , t ) = [ Δ p ^ ] 0 2 ( n ) , t < 0 , [ Δ p ^ ] 1 2 ( n , t ) , 0 < t < τ , [ Δ p ^ ] 2 2 ( n , t ) , t > τ ,
being [46]
[ Δ p ^ ] 0 2 ( n ) = n + 1 2 m 0 ω 0 .
Through Equations (24), (45), and (48), we find, for the interval 0 < t < τ , the expression [35,40]
[ Δ p ^ ] 1 2 ( n , t ) = ω 1 2 ω 0 2 sin 2 ( ω 1 t ) + cos 2 ( ω 1 t ) [ Δ p ^ ] 0 2 ( n ) .
To calculate the variance of the operator p ^ in the interval t > τ , we must use Equations (16), (17), and (38) in Equation (24). The general behavior of [ Δ p ^ ] 2 ( n , t ) , when n = 0 , is schematized in Figure 5.
Thus, it is direct to see that the uncertainty relation between these operators in the interval 0 < t < τ has the form
[ Δ x ^ ] 1 2 ( n , t ) [ Δ p ^ ] 1 2 ( n , t ) [ Δ x ^ ] 0 2 ( n ) [ Δ p ^ ] 0 2 ( n ) { 1 + [ ω 1 2 ω 0 2 ω 0 ω 1 sin ( ω 1 t ) × cos ( ω 1 t ) ] 2 } ,
and the uncertainty relation for the interval t > τ is obtained in a similar way. Clearly, when ω 1 = ω 0 , the uncertainty relation (Equation (55)) falls back to the uncertainty relation of a time-independent oscillator [46]. Furthermore, the uncertainty relation for an oscillator with time-independent frequency is also reobtained when τ = τ l , as shown in Figure 6.
The results found by us (via the LR method) given in Equations (45), (46), (48), (51), and (54) are in agreement with those found in ref. [40]. Hereafter, we use the LR method to obtain new results concerning the model given in Equation (30).

3.4. Mean Energy

The expected value of the Hamiltonian operator is identified as the mean energy of the system, that is, E n , t = H ^ ( t ) ( n , t ) . The mean energy for the model in Equation (30) can be written as
E ( n , t ) = E 0 ( n ) , t < 0 , E 1 ( n , t ) , 0 < t < τ , E 2 ( n , t ) , t > τ ,
wherein [46]
E 0 ( n ) = n + 1 2 ω 0 .
For the interval 0 < t < τ , through Equations (22), (51), (54), and (57), we find that the mean energy is time-independent, given by
E 1 ( n , t ) = 1 2 1 + ω 1 2 ω 0 2 E 0 ( n ) .
When ω 1 = ω 0 , Equation (58) reduces to E 1 ( n , t ) = E 0 ( n ) , as expected. Note that for ω 1 / ω 0 < 1 , E 1 ( n , t ) < E 0 ( n ) , whereas for ω 1 / ω 0 > 1 , we have E 1 ( n , t ) > E 0 ( n ) . When ω 1 = 0 , which means that the system is free in interval 0 < t < τ , we have E 1 ( n , t ) = E 0 ( n ) / 2 . This can also be obtained by making ω 0 ω 1 , which leads to E 1 ( n , t ) E 0 ( n ) / 2 .
For the interval t > τ , from Equations (16), (17), (22), (38), and (57), we have that the mean energy E 2 ( n , t ) is given by
E 2 ( n , t ) = 1 + 1 2 ω 1 2 ω 0 2 ω 0 ω 1 2 sin 2 ( ω 1 τ ) E 0 ( n ) .
Note that Equation (59) is independent of t, and E 2 ( n , t ) E 0 ( n ) , even for ω 1 < ω 0 . The behavior of the ratio E 2 ( n , t ) / E 0 ( n ) is illustrated in Figure 7.
For ω 1 = ω 0 (dashed line in Figure 7), Equation (59) recovers E 2 ( n , t ) = E 0 ( n ) , as expected. Furthermore, for τ = τ l (dotted lines in Figure 7), Equation (59) also gives E 2 ( n , t ) = E 0 ( n ) . In particular, the energy of this system is maximized when τ = τ l + 1 / 2 (dot–dashed lines in Figure 7). For ω 1 / ω 0 0 , we obtain E 2 ( n , t ) = 1 + ω 0 2 τ 2 / 2 E 0 ( n ) . Moreover, from Equations (46) and (59), we obtain
E 2 ( n , t ) = 2 cosh 2 r 2 ( t ) 1 E 0 ( n ) .
Thus, while there is squeeze, E 2 ( n , t ) > E 0 ( n ) . Therefore, the squeezing caused by the frequency jumps results in an increase in the mean energy of the oscillator with respect to its initial energy E 0 ( n ) .
It is interesting to investigate the behavior of E 2 ( n , t ) / E 1 ( n , t ) . Unlike the ratio E 2 ( n , t ) / E 0 ( n ) , which is such that E 2 ( n , t ) / E 0 ( n ) 1 , the ratio E 2 ( n , t ) / E 1 ( n , t ) can be lesser than, equal to, or greater than one, as shown in Figure 8. More specifically, for ω 1 / ω 0 < 1 , we have E 2 ( n , t ) / E 1 ( n , t ) > 1 , and when ω 1 / ω 0 0 , we find E 2 ( n , t ) = 2 + ω 0 2 τ 2 E 1 ( n , t ) . For ω 1 / ω 0 > 1 , the ratio E 2 ( n , t ) / E 1 ( n , t ) oscillates between zero and one (see Figure 8). We highlight that, besides the trivial case ω 1 = ω 0 , there are other values of the ratio ω 1 / ω 0 that result in E 2 ( n , t ) = E 1 ( n , t ) .

3.5. Mean Number of Excitations

The mean number of excitations, N ^ ( n , t ) , for the model in Equation (30), is given by
N ^ ( n , t ) = N ^ 0 ( n ) , t < 0 , N ^ 1 ( n , t ) , 0 < t < τ , N ^ 2 ( n , t ) , t > τ ,
where [46]
N ^ 0 ( n ) = n .
Given this, for the interval 0 < t < τ , by means of Equations (26) and (45), we have
N ^ 1 ( n , t ) = n + n + 1 2 1 2 ω 1 2 ω 0 2 ω 0 ω 1 2 sin 2 ( ω 1 t ) .
We remark the time dependence in N ^ 1 ( n , t ) , whereas the mean energy in the interval 0 < t < τ (Equation (58)) is time-independent.
For the interval t > τ , through Equations (26) and (46), we obtain N ^ 2 ( n , t ) , given by
N ^ 2 ( n , t ) = N ^ 1 ( n , τ ) .
We also remark that for τ = τ l , the behavior of the system returns to that of the time-independent oscillator found before the frequency jumps, i.e, N ^ 2 ( n , t ) = N ^ 0 ( n ) . We highlight that there is excitation even for n = 0 , which means that, under jumps in its frequency, a quantum oscillator initially in its ground state can become excited (a classical oscillator in its ground state would remain in the same state). In addition, excitation can also occur when ω 1 / ω 0 < 1 , as shown in Figure 9.
Using Equations (59) and (64), we can also write E 2 ( n , t ) as
E 2 ( n , t ) = N ^ 2 ( n , t ) + 1 2 ω 0 ,
which has the same structure as the expression for the energy eigenvalues of a time-independent oscillator (see Equation (57)). The time evolution of N ^ ( 0 , t ) , given by Equation (61) for n = 0 , is shown in Figure 10.

3.6. Transition Probability

The general transition probability, P ( t ) m n , for the model in Equation (30), is given by
P ( t ) m n = δ m , n , t < 0 , P 1 ( t ) m n , 0 < t < τ , P 2 ( t ) m n , t > τ ,
with δ m , n being the Kronecker delta. Using Equations (28) and (45) (or Equations (29) and (63) with n = 0 ), we find, for the interval 0 < t < τ , P 1 ( t ) m n , whose expression is
P 1 ( t ) m n = 2 m + n min ( m , n ) ! 2 ω 1 2 ω 0 2 2 ω 0 ω 1 sin ( ω 1 t ) | n m | m ! n ! 1 + ω 1 2 ω 0 2 2 ω 0 ω 1 2 sin 2 ( ω 1 t ) 1 2 × k = | n m | 2 n + m 2 n + m 2 k n + m + 2 k 2 4 n + m 2 k ! k | n m | 2 ! 1 + ω 1 2 ω 0 2 2 ω 0 ω 1 2 sin 2 ( ω 1 t ) k 2 2 ,
for even values of | n m | , and P ( t ) m n = 0 for odd values of | n m | . Since the r 1 ( t ) parameter (Equation (45)) in this interval is explicitly time-dependent, consequently, the transition probability m n will also depend.
For the interval t > τ , using Equation (46) in Equation (28) (or Equations (29) and (64) with n = 0 ), we see that the result is
P 2 ( t ) m n = P 1 ( τ ) m n .
Thus, even when the frequency returns to ω 0 after an instant τ , one can find a nonzero m n transition probability, depending on the value of τ . It is noticeable from Equation (67) that P ( t ) m n = P ( t ) n m , with such symmetry being a consequence of the parity of the potential in Equation (7) [61]. The behavior of P ( t ) 1 n is illustrated in Figure 11 (for τ = 3 τ 1 / 2 ) and Figure 12 (for τ = τ 1 ). In Figure 11, as n increases, P ( t ) 1 n decreases. We highlight that for t > τ = τ 1 in Figure 12, P ( t ) 1 1 = 1 or, in other words, the oscillator remains in its same initial state.
In addition, for m = 0 , Equation (68) gives
P 2 ( t ) 0 n = n ! ω 1 2 ω 0 2 2 ω 0 ω 1 n sin n ( ω 1 τ ) 2 n n 2 ! 2 1 + ω 1 2 ω 0 2 2 ω 0 ω 1 2 sin 2 ( ω 1 τ ) n + 1 2 ,
where n = 0 , 2 , 4 , , recovering one of the results found in ref. [40]. Thus, Equation (68) generalizes the result for the transition probability found in ref. [40]. By also making n = 0 in Equation (69), we find the probability of the oscillator of persisting in the fundamental state, and the probability of the oscillator being excited, after the frequency returns to ω 0 , which is given by 1 P 2 ( t ) 0 0 .

4. Final Remarks

Using the Lewis–Riesenfeld method, we investigated the dynamics of a quantum harmonic oscillator that undergoes two abrupt jumps in its frequency (Equation (30)). We reobtained the analytical formulas of ref. [40] for the squeeze parameters (Equations (45), (46), and (48)), the quantum fluctuations of the position (Equation (51)) and momentum (Equation (54)) operators, and the probability amplitude of a transition from the fundamental state to an arbitrary energy eigenstate (Equation (69)). We also obtained expressions for the mean energy value (Equations (58) and (59)), the mean number of excitations (Equations (63) and (64)) (which were not calculated in ref. [40]), and for the transition probabilities considering the initial state different from the fundamental (Equations (67) and (68)) (which generalizes the formula found in ref. [40]).
We found that, as expected, the mean energy of the system is independent of time in each one of the intervals: t < 0 , 0 < t < τ , and t > τ . Moreover, we showed that the mean energy of the oscillator after the jumps is equal or greater than that before these jumps, even when ω 1 < ω 0 . We also obtained, for t > τ τ l , a non-null value for the mean number of excitations when the oscillator starts in the fundamental state (Equations (58) and (59) with n = 0 ), which means that, under the jumps in its frequency, a quantum oscillator, initially in the ground state, can become excited. We showed that transitions between arbitrary m and n states only occur if | n m | is an even number. We highlighted that, for t > τ τ l and a fixed value of m, as n increases, P ( t ) m n decreases. Finally, we showed that, for t > τ = τ l , P ( t ) m n = δ m , n , so that the oscillator returns to the same initial state (this generalizes, for any initial state m, the result found in ref. [40] for m = 0 ).

Author Contributions

Conceptualization, S.S.C., L.Q. and D.T.A.; investigation, S.S.C., L.Q. and D.T.A.; writing—original draft preparation, S.S.C.; writing—review and editing, S.S.C., L.Q. and D.T.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)—Brazil, by the program PIBIC/CNPq, project 144456/2020-6, and Fundação Amazônia de Amparo a Estudos e Pesquisas (Fapespa) - Brazil, by the program PIBIC/Fapespa, and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)—Brazil, Finance Code 001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Alexandre Costa, Andreson Rego, Edson Nogueira, and Van Sérgio Alves for valuable discussions, as well as Adolfo del Campo, Bogdan M. Mihalcea, Daniel Tibaduiza, and Victor Dodonov for their valuable suggestions to this paper. S.S.C. was partially supported by CNPq—Brazil (program PIBIC/CNPq—project 144456/2020-6), Fapespa—Brazil (program PIBIC/Fapespa), and CAPES—Brazil (Finance Code 001). L.Q. was also supported by CAPES, Finance Code 001.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Husimi, K. Miscellanea in Elementary Quantum Mechanics, II. Prog. Theor. Phys. 1953, 9, 381–402. [Google Scholar] [CrossRef]
  2. Lewis, H.R. Classical and Quantum Systems with Time-Dependent Harmonic-Oscillator-Type Hamiltonians. Phys. Rev. Lett. 1967, 18, 510–512. [Google Scholar] [CrossRef]
  3. Lewis, H.R. Class of Exact Invariants for Classical and Quantum Time-Dependent Harmonic Oscillators. J. Math. Phys. 1968, 9, 1976–1986. [Google Scholar] [CrossRef]
  4. Lewis, H.R.; 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]
  5. Pedrosa, I.A. Exact wave functions of a harmonic oscillator with time-dependent mass and frequency. Phys. Rev. A 1997, 55, 3219–3221. [Google Scholar] [CrossRef]
  6. Ciftja, O. A simple derivation of the exact wavefunction of a harmonic oscillator with time-dependent mass and frequency. J. Phys. A. Math. Gen. 1999, 32, 6385–6389. [Google Scholar] [CrossRef]
  7. Guasti, M.F.; Moya-Cessa, H. Solution of the Schrödinger equation for time-dependent 1D harmonic oscillators using the orthogonal functions invariant. J. Phys. A. Math. Gen. 2003, 36, 2069–2076. [Google Scholar] [CrossRef]
  8. Pedrosa, I.A.; Rosas, A. Electromagnetic Field Quantization in Time-Dependent Linear Media. Phys. Rev. Lett. 2009, 103, 010402. [Google Scholar] [CrossRef]
  9. Pedrosa, I.A. Quantum electromagnetic waves in nonstationary linear media. Phys. Rev. A 2011, 83, 032108. [Google Scholar] [CrossRef]
  10. Dodonov, V.; Man’ko, V.; Polynkin, P. Geometrical squeezed states of a charged particle in a time-dependent magnetic field. Phys. Lett. A 1994, 188, 232–238. [Google Scholar] [CrossRef]
  11. Xu, X.-W.; Ren, T.-Q.; Liu, S.-D. Analytic solution for one-dimensional quantum oscillator with a variable frequency. Acta Phys. Sin. 1999, 8, 641–645. [Google Scholar] [CrossRef]
  12. Aguiar, V.; Guedes, I. Entropy and information of a spinless charged particle in time-varying magnetic fields. J. Math. Phys. 2016, 57, 092103. [Google Scholar] [CrossRef]
  13. Dodonov, V.V.; Horovits, M.B. Squeezing of Relative and Center-of-Orbit Coordinates of a Charged Particle by Step-Wise Variations of a Uniform Magnetic Field with an Arbitrary Linear Vector Potential. J. Russ. Laser Res. 2018, 39, 389–400. [Google Scholar] [CrossRef]
  14. Brown, L.S. Quantum motion in a Paul trap. Phys. Rev. Lett. 1991, 66, 527–529. [Google Scholar] [CrossRef]
  15. Agarwal, G.S.; Kumar, S.A. Exact quantum-statistical dynamics of an oscillator with time-dependent frequency and generation of nonclassical states. Phys. Rev. Lett. 1991, 67, 3665–3668. [Google Scholar] [CrossRef]
  16. Mihalcea, B.M. A quantum parametric oscillator in a radiofrequency trap. Phys. Scr. 2009, T135, 014006. [Google Scholar] [CrossRef] [Green Version]
  17. Aguiar, V.; Nascimento, J.; Guedes, I. Exact wave functions and uncertainties for a spinless charged particle in a time-dependent Penning trap. Int. J. Mass Spectrom. 2016, 409, 21–28. [Google Scholar] [CrossRef]
  18. Menicucci, N.C.; Milburn, G.J. Single trapped ion as a time-dependent harmonic oscillator. Phys. Rev. A 2007, 76, 052105. [Google Scholar] [CrossRef] [Green Version]
  19. Pedrosa, I.A. On the Quantization of the London Superconductor. Braz. J. Phys. 2021, 51, 401–405. [Google Scholar] [CrossRef]
  20. Choi, J.R. Interpreting quantum states of electromagnetic field in time-dependent linear media. Phys. Rev. A 2010, 82, 055803. [Google Scholar] [CrossRef]
  21. Salamon, P.; Hoffmann, K.H.; Rezek, Y.; Kosloff, R. Maximum work in minimum time from a conservative quantum system. Phys. Chem. Chem. Phys. 2009, 11, 1027–1032. [Google Scholar] [CrossRef] [PubMed]
  22. Schaff, J.F.; Song, X.L.; Vignolo, P.; Labeyrie, G. Fast optimal transition between two equilibrium states. Phys. Rev. A 2010, 82, 033430. [Google Scholar] [CrossRef] [Green Version]
  23. Chen, X.; Ruschhaupt, A.; Schmidt, S.; del Campo, A.; Guéry-Odelin, D.; Muga, J.G. Fast Optimal Frictionless Atom Cooling in Harmonic Traps: Shortcut to Adiabaticity. Phys. Rev. Lett. 2010, 104, 063002. [Google Scholar] [CrossRef]
  24. Stefanatos, D.; Ruths, J.; Li, J.S. Frictionless atom cooling in harmonic traps: A time-optimal approach. Phys. Rev. A 2010, 82, 063422. [Google Scholar] [CrossRef] [Green Version]
  25. Dupays, L.; Spierings, D.C.; Steinberg, A.M.; del Campo, A. Delta-kick cooling, time-optimal control of scale-invariant dynamics, and shortcuts to adiabaticity assisted by kicks. Phys. Rev. Res. 2021, 3, 033261. [Google Scholar] [CrossRef]
  26. Martínez-Tibaduiza, D.; Pires, L.; Farina, C. Time-dependent quantum harmonic oscillator: A continuous route from adiabatic to sudden changes. J. Phys. B At. Mol. Opt. Phys. 2021, 54, 205401. [Google Scholar] [CrossRef]
  27. Landim, R.R.; Guedes, I. Wave functions for a Dirac particle in a time-dependent potential. Phys. Rev. A 2000, 61, 054101. [Google Scholar] [CrossRef]
  28. Gao, X.C.; Fu, J.; Li, X.H.; Gao, J. Invariant formulation and exact solutions for the relativistic charged Klein-Gordon field in a time-dependent spatially homogeneous electric field. Phys. Rev. A 1998, 57, 753–761. [Google Scholar] [CrossRef]
  29. Dodonov, V.V.; Klimov, A.B. Generation and detection of photons in a cavity with a resonantly oscillating boundary. Phys. Rev. A 1996, 53, 2664–2682. [Google Scholar] [CrossRef]
  30. Dodonov, V.; Klimov, A.; Man’ko, V. Generation of squeezed states in a resonator with a moving wall. Phys. Lett. A 1990, 149, 225–228. [Google Scholar] [CrossRef]
  31. Pedrosa, I.A.; Guedes, I. Exact quantum states of an inverted pendulum under time-dependent gravitation. Int. J. Mod. Phys. A 2004, 19, 4165–4172. [Google Scholar] [CrossRef]
  32. Carvalho, A.M.d.M.; Furtado, C.; Pedrosa, I.A. Scalar fields and exact invariants in a Friedmann-Robertson-Walker spacetime. Phys. Rev. D 2004, 70, 123523. [Google Scholar] [CrossRef]
  33. Greenwood, E. Time-dependent particle production and particle number in cosmological de Sitter space. Int. J. Mod. Phys. D 2015, 24, 1550031. [Google Scholar] [CrossRef] [Green Version]
  34. Janszky, J.; Yushin, Y. Squeezing via frequency jump. Opt. Commun. 1986, 59, 151–154. [Google Scholar] [CrossRef]
  35. Janszky, J.; Adam, P. Strong squeezing by repeated frequency jumps. Phys. Rev. A 1992, 46, 6091–6092. [Google Scholar] [CrossRef] [PubMed]
  36. Kiss, T.; Janszky, J.; Adam, P. Time evolution of harmonic oscillators with time-dependent parameters: A step-function approximation. Phys. Rev. A 1994, 49, 4935–4942. [Google Scholar] [CrossRef] [PubMed]
  37. Moya-Cessa, H.; Fernández Guasti, M. Coherent states for the time dependent harmonic oscillator: The step function. Phys. Lett. A 2003, 311, 1–5. [Google Scholar] [CrossRef] [Green Version]
  38. Stefanatos, D. Minimum-Time Transitions between Thermal and Fixed Average Energy States of the Quantum Parametric Oscillator. SIAM J. Control Optim. 2017, 55, 1429–1451. [Google Scholar] [CrossRef] [Green Version]
  39. Stefanatos, D. Minimum-Time Transitions Between Thermal Equilibrium States of the Quantum Parametric Oscillator. IEEE Trans. Automat. Contr. 2017, 62, 4290–4297. [Google Scholar] [CrossRef] [Green Version]
  40. Tibaduiza, D.M.; Pires, L.; Szilard, D.; Zarro, C.A.D.; Farina, C.; Rego, A.L.C. A Time-Dependent Harmonic Oscillator with Two Frequency Jumps: An Exact Algebraic Solution. Braz. J. Phys. 2020, 50, 634–646. [Google Scholar] [CrossRef]
  41. Tibaduiza, D.M.; Pires, L.; Rego, A.L.C.; Szilard, D.; Zarro, C.; Farina, C. Efficient algebraic solution for a time-dependent quantum harmonic oscillator. Phys. Scr. 2020, 95, 105102. [Google Scholar] [CrossRef]
  42. Pedrosa, I.A.; Serra, G.P.; Guedes, I. Wave functions of a time-dependent harmonic oscillator with and without a singular perturbation. Phys. Rev. A 1997, 56, 4300–4303. [Google Scholar] [CrossRef]
  43. Xin, M.; Leong, W.S.; Chen, Z.; Wang, Y.; Lan, S.Y. Rapid Quantum Squeezing by Jumping the Harmonic Oscillator Frequency. Phys. Rev. Lett. 2021, 127, 183602. [Google Scholar] [CrossRef] [PubMed]
  44. Wolf, F.; Shi, C.; Heip, J.C.; Gessner, M.; Pezzè, L.; Smerzi, A.; Schulte, M.; Hammerer, K.; Schmidt, P.O. Motional Fock states for quantum-enhanced amplitude and phase measurements with trapped ions. Nat. Commun. 2019, 10, 2929. [Google Scholar] [CrossRef] [Green Version]
  45. Choi, J.R. The dependency on the squeezing parameter for the uncertainty relation in the squeezed states of the time-dependent oscillator. Int. J. Mod. Phys. B 2004, 18, 2307–2324. [Google Scholar] [CrossRef]
  46. Sakurai, J.J.; Napolitano, J. Modern Quantum Mechanics, 3rd ed.; Cambridge University Press: Cambridge, UK, 2020; pp. 83–88. [Google Scholar]
  47. Griffiths, D.J. Introduction to Quantum Mechanics, 3rd ed.; Cambridge University Press: Cambridge, UK, 2018; pp. 39–55. [Google Scholar]
  48. Cohen-Tannoudji, C.; Diu, B.; Laloe, F. Quantum Mechanics, Volume 1: Basic Concepts, Tools, and Applications, 2nd ed.; Wiley-VCH: Weinheim, Germany, 2019; pp. 502–521. [Google Scholar]
  49. Prykarpatskyy, Y. Steen–Ermakov–Pinney Equation and Integrable Nonlinear Deformation of the One-Dimensional Dirac Equation. J. Math. Sci. 2018, 231, 820–826. [Google Scholar] [CrossRef] [Green Version]
  50. Pinney, E. The nonlinear differential equation y″ + p(x)y + cy−3 = 0. Proc. Am. Math. Soc. 1950, 1, 681. [Google Scholar] [CrossRef]
  51. de Lima, A.L.; Rosas, A.; Pedrosa, I. Quantum dynamics of a particle trapped by oscillating fields. J. Mod. Opt. 2009, 56, 75–80. [Google Scholar] [CrossRef]
  52. Cariñena, J.F.; de Lucas, J. Applications of Lie systems in dissipative Milne-Pinney equations. Int. J. Geom. Methods Mod. Phys. 2009, 06, 683–699. [Google Scholar] [CrossRef] [Green Version]
  53. Weber, H.J.; Arfken, G.B. Essential Mathematical Methods for Physicists, 6th ed.; Academic Press: San Diego, CA, USA, 2003; pp. 638–642. [Google Scholar]
  54. Pedrosa, I.A. Comment on “Coherent states for the time-dependent harmonic oscillator”. Phys. Rev. D 1987, 36, 1279–1280. [Google Scholar] [CrossRef]
  55. Daneshmand, R.; Tavassoly, M.K. Dynamics of Nonclassicality of Time- and Conductivity-Dependent Squeezed States and Excited Even/Odd Coherent States. Commun. Theor. Phys. 2017, 67, 365–376. [Google Scholar] [CrossRef]
  56. Guerry, C.C.; Knight, P.L. Introductory Quantum Optics, 1st ed.; Cambridge University Press: Cambridge, UK, 2005; pp. 150–165. [Google Scholar]
  57. Kim, M.S.; de Oliveira, F.A.M.; Knight, P.L. Properties of squeezed number states and squeezed thermal states. Phys. Rev. A 1989, 40, 2494–2503. [Google Scholar] [CrossRef] [PubMed]
  58. Marian, P. Higher-order squeezing and photon statistics for squeezed thermal states. Phys. Rev. A 1992, 45, 2044–2051. [Google Scholar] [CrossRef]
  59. Moeckel, M.; Kehrein, S. Real-time evolution for weak interaction quenches in quantum systems. Ann. Phys. (N. Y). 2009, 324, 2146–2178. [Google Scholar] [CrossRef] [Green Version]
  60. Kim, M.; de Oliveira, F.; Knight, P. Photon number distributions for squeezed number states and squeezed thermal states. Opt. Commun. 1989, 72, 99–103. [Google Scholar] [CrossRef]
  61. Popov, V.; Perelomov, A. Parametric Excitation of a Quantum Oscillator. Sov. J. Exp. Theor. Phys. 1969, 30, 1375–1390. [Google Scholar]
Figure 1. General behavior of ρ ( t ) , as a function of ω 0 t , with ω 1 = 3 ω 0 and different values of τ (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units). For τ = [ τ l + 1 / 2 ] l = 2 = 5 π / 6 , we note that the amplitude of oscillation for the ρ ( t ) parameter is maximum. For τ = 59 π / 62 , the amplitude of oscillation is intermediate. Finally, when τ = [ τ l ] l = 3 = π , there is no oscillation at all, and it follows that ρ ( t ) becomes time independent.
Figure 1. General behavior of ρ ( t ) , as a function of ω 0 t , with ω 1 = 3 ω 0 and different values of τ (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units). For τ = [ τ l + 1 / 2 ] l = 2 = 5 π / 6 , we note that the amplitude of oscillation for the ρ ( t ) parameter is maximum. For τ = 59 π / 62 , the amplitude of oscillation is intermediate. Finally, when τ = [ τ l ] l = 3 = π , there is no oscillation at all, and it follows that ρ ( t ) becomes time independent.
Entropy 24 01851 g001
Figure 2. Behavior of the squeeze parameter r ( t ) as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Figure 2. Behavior of the squeeze parameter r ( t ) as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Entropy 24 01851 g002
Figure 3. Behavior of the squeeze phase ϕ ( t ) as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Figure 3. Behavior of the squeeze phase ϕ ( t ) as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Entropy 24 01851 g003
Figure 4. Behavior of the [ Δ x ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Figure 4. Behavior of the [ Δ x ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Entropy 24 01851 g004
Figure 5. Behavior of the [ Δ p ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Figure 5. Behavior of the [ Δ p ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Entropy 24 01851 g005
Figure 6. Behavior of the [ Δ x ^ ] 2 ( n , t ) [ Δ p ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Figure 6. Behavior of the [ Δ x ^ ] 2 ( n , t ) [ Δ p ^ ] 2 ( n , t ) for n = 0 as a function of ω 0 t , where ω 1 = 3 ω 0 (we consider, for simplicity, m 0 = ω 0 = 1 and = 1 in arbitrary units).
Entropy 24 01851 g006
Figure 7. Ratio E 2 ( n , t ) / E 0 ( n ) as a function of ω 0 τ and ω 1 / ω 0 . The dashed line corresponds to ω 1 = ω 0 . The dotted lines correspond to τ = τ l . The dot–dashed lines correspond to τ = τ l + 1 / 2 .
Figure 7. Ratio E 2 ( n , t ) / E 0 ( n ) as a function of ω 0 τ and ω 1 / ω 0 . The dashed line corresponds to ω 1 = ω 0 . The dotted lines correspond to τ = τ l . The dot–dashed lines correspond to τ = τ l + 1 / 2 .
Entropy 24 01851 g007
Figure 8. Ratio E 2 ( n , t ) / E 1 ( n , t ) as a function of ω 1 / ω 0 . For ω 1 / ω 0 < 1 , E 2 ( n , t ) / E 1 ( n , t ) > 1 , whereas for ω 1 / ω 0 > 1 , E 2 ( n , t ) / E 1 ( n , t ) 1 .
Figure 8. Ratio E 2 ( n , t ) / E 1 ( n , t ) as a function of ω 1 / ω 0 . For ω 1 / ω 0 < 1 , E 2 ( n , t ) / E 1 ( n , t ) > 1 , whereas for ω 1 / ω 0 > 1 , E 2 ( n , t ) / E 1 ( n , t ) 1 .
Entropy 24 01851 g008
Figure 9. Some examples of mean number of excitations N ^ 2 ( 0 , t ) that an oscillator could undergo as a function of ω 1 / ω 0 for different values of ω 0 τ .
Figure 9. Some examples of mean number of excitations N ^ 2 ( 0 , t ) that an oscillator could undergo as a function of ω 1 / ω 0 for different values of ω 0 τ .
Entropy 24 01851 g009
Figure 10. Mean number of excitations N ^ ( 0 , t ) as a function of ω 0 t for different values of τ , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Figure 10. Mean number of excitations N ^ ( 0 , t ) as a function of ω 0 t for different values of τ , where ω 1 = 3 ω 0 (we consider ω 0 = 1 in arbitrary units).
Entropy 24 01851 g010
Figure 11. Behavior of P ( t ) 1 n , as a function of ω 0 t , with ω 1 = 2 ω 0 and τ = 3 τ 1 / 2 = 3 π / 4 (we consider ω 0 = 1 in arbitrary units).
Figure 11. Behavior of P ( t ) 1 n , as a function of ω 0 t , with ω 1 = 2 ω 0 and τ = 3 τ 1 / 2 = 3 π / 4 (we consider ω 0 = 1 in arbitrary units).
Entropy 24 01851 g011
Figure 12. Behavior of P ( t ) 1 n , as a function of ω 0 t , with ω 1 = 2 ω 0 and τ = τ 1 = π / 2 (we consider ω 0 = 1 in arbitrary units).
Figure 12. Behavior of P ( t ) 1 n , as a function of ω 0 t , with ω 1 = 2 ω 0 and τ = τ 1 = π / 2 (we consider ω 0 = 1 in arbitrary units).
Entropy 24 01851 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Coelho, S.S.; Queiroz, L.; Alves, D.T. Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method. Entropy 2022, 24, 1851. https://doi.org/10.3390/e24121851

AMA Style

Coelho SS, Queiroz L, Alves DT. Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method. Entropy. 2022; 24(12):1851. https://doi.org/10.3390/e24121851

Chicago/Turabian Style

Coelho, Stanley S., Lucas Queiroz, and Danilo T. Alves. 2022. "Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method" Entropy 24, no. 12: 1851. https://doi.org/10.3390/e24121851

APA Style

Coelho, S. S., Queiroz, L., & Alves, D. T. (2022). Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method. Entropy, 24(12), 1851. https://doi.org/10.3390/e24121851

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