Next Article in Journal
Evaluation of Transport Properties and Energy Conversion of Bacterial Cellulose Membrane Using Peusner Network Thermodynamics
Next Article in Special Issue
Superior Resilience of Non-Gaussian Entanglement against Local Gaussian Noises
Previous Article in Journal
Matter-Aggregating Low-Dimensional Nanostructures at the Edge of the Classical vs. Quantum Realm
Previous Article in Special Issue
Exact Solution of a Time-Dependent Quantum Harmonic Oscillator with Two Frequency Jumps via the Lewis–Riesenfeld Dynamical Invariant Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero

by
Viktor V. Dodonov
1,2,* and
Alexandre V. Dodonov
1,2
1
Institute of Physics, University of Brasilia, P.O. Box 04455, Brasilia 70919-970, DF, Brazil
2
International Center of Physics, University of Brasilia, Brasilia 70910-900, DF, Brazil
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(1), 2; https://doi.org/10.3390/e25010002
Submission received: 12 November 2022 / Revised: 11 December 2022 / Accepted: 14 December 2022 / Published: 20 December 2022
(This article belongs to the Special Issue Quantum Nonstationary Systems)

Abstract

:
We study the evolution of the energy of a harmonic oscillator when its frequency slowly varies with time and passes through a zero value. We consider both the classical and quantum descriptions of the system. We show that after a single frequency passage through a zero value, the famous adiabatic invariant ratio of energy to frequency (which does not hold for a zero frequency) is reestablished again, but with the proportionality coefficient dependent on the initial state. The dependence on the initial state disappears after averaging over the phases of initial states with the same energy (in particular, for the initial vacuum, the Fock and thermal quantum states). In this case, the mean proportionality coefficient is always greater than unity. The concrete value of the mean proportionality coefficient depends on the power index of the frequency dependence on a time near the zero point. In particular, the mean energy triplicates if the frequency tends to zero linearly. If the frequency attains zero more than once, the adiabatic proportionality coefficient strongly depends on the lengths of time intervals between zero points, so that the mean energy behavior becomes quasi-stochastic after many passages through a zero value. The original Born–Fock theorem does not work after the frequency passes through zero. However, its generalization is found: the initial Fock state becomes a wide superposition of many Fock states, whose weights do not depend on time in the new adiabatic regime. When the mean energy triplicates, the initial Nth Fock state becomes a superposition of, roughly speaking, 6N states, distributed nonuniformly. The initial vacuum and low-order Fock states become squeezed, as well as the initial thermal states with low values of the mean energy.

1. Introduction

One of many brilliant results of classical and quantum mechanics is the existence of adiabatic invariants in the case where parameters of a system vary slowly with time. The simplest invariant is the ratio of the energy of a harmonic oscillator E ( t ) to its time-dependent frequency ω ( t ) [1]:
E ( t ) / ω ( t ) = c o n s t if | ω ˙ | / ω 2 ( t ) 1 .
Then, suppose that the frequency returns to its initial value after some slow variations. What is the final energy of the oscillator? According to Equation (1), the answer is obvious: the final energy coincides with the initial one. However, there exists a remarkable exclusion from this result, when the frequency passes through a zero value in the process of evolution, so that the condition of validity of Equation (1) is obviously broken for any rate of the evolution. The goal of our paper is to study the dependence of the final energy on the shape of the time-dependent frequency, when this frequency slowly passes through a zero value. We perform analytic calculations for the quantum oscillator and numeric calculations for the classical oscillator. Note that various aspects of the harmonic oscillator evolution in the adiabatic regime have been studied by many authors for decades (see, e.g., papers [2,3,4,5,6,7,8]). However, the situation when the frequency passes slowly through a zero value have not been considered in the known publications. An oscillator whose frequency exponentially goes adiabatically and asymptotically to zero has been considered, e.g., in papers [9,10]. However, we are interested in the case when the frequency passes through zero and returns to its initial value.
The structure of the paper is as follows. In Section 2, we describe the results of the numeric solutions of the classical equations of motion for two kinds of time dependence of the frequency: ω 2 ( t ) | t | n and ω 2 ( t ) | tanh ( a t ) | n . Several figures show the differences between the cases of fast and slow frequency variations, when the frequency does not attain a zero value and when it passes through a zero value. In Section 3, we present general formulas describing the evolution of the mean oscillator energy in the quantum case, including adiabatic regimes without and with a zero-crossing frequency value. Exact solutions for the power profile of the frequency are derived and analyzed in Section 4. Transition rules in the case of a single frequency passage through zero are obtained in Section 5. An example of tanh-like frequency functions is considered in Section 6. Double transitions of the frequency through zero values are studied in Section 7 in the quantum and classical cases. Section 8, Section 9 and Section 10 are devoted to the energy fluctuations, the violation and generalization of the Born–Fock theorem and the appearance of squeezing, respectively. Section 11 contains a discussion of the main results.

2. Evolution of the Oscillator Energy in the Classical Case

The basic equation has the form
x ¨ + ω 2 ( t ) x = 0 .
We consider a special case when the time-dependent frequency can be written as ω 2 ( t ) = ω 0 2 f ( t / τ ) , where f ( t / τ ) is a non-negative function with the properties
f ( 1 ) = 1 , f ( 0 ) = 0 , τ > 0 .
We assume that ω 2 ( t ) = ω 0 2 = c o n s t for t τ . Since we are interested mainly in the evolution of the energy at t > τ , we consider a one-parameter family of classical trajectories with the same initial energy E 0 . Then, assuming the particle mass m = 1 , it is convenient to use the initial coordinate x 0 and initial velocity x ˙ 0 in the form
x 0 = cos ( φ ) 2 E 0 / ω 0 , x ˙ 0 = sin ( φ ) 2 E 0 , 0 φ < 2 π .
To solve Equation (2) numerically, we introduce dimensionless variables X = ω 0 x / 2 E 0 and T = t / τ , arriving at the equation for T 1
d 2 X / d T 2 + G 2 f ( T ) X = 0 , G = ω 0 τ ,
with the initial conditions
X ( 1 ) = cos ( φ ) , d X / d T | T = 1 = G sin ( φ ) .
Then, the dimensionless energy ratio R = E ( t ) / E 0 depends on the dimensionless time T and two parameters, G and φ :
R ( T ; G ; φ ) = f ( T ) X 2 ( T ) + G 2 d X / d T 2 ,
where the derivative d X / d T must be taken at instant T. The existence of the adiabatic invariant (1) implies that R ( T ; G ) = f ( T ) for T < 0 and big enough values of parameter G, independently on the values of parameter φ. In the following subsections, we study what can happen if T > 0 , for different families of functions f ( T ) .

2.1. A Power Profile of the Frequency

Our first example is the power profile f ( T ) = | T | n with n > 0 . We solved Equation (4) numerically and calculated the dimensionless ratio
R n ( T ; G ; φ ) = | T | n X 2 ( T ) + G 2 d X / d T 2 .
The adiabatic ratio must equal R n ( T ; G ) = | T | n / 2 . In Figure 1 and Figure 2, we show R n as function of φ for the fixed values T = 1 / 2 , n = 1 / 2 , n = 2 and n = 4 , comparing different behaviors when G = 1 (no adiabaticity) and G = 1000 . Figure 1 shows a strong dependence of the energy on the phase φ in the nonadiabatic regime. However, this dependence becomes negligible in Figure 2, which shows that the mean value of R n is, indeed, | T | n / 2 .
However, the situation becomes quite different for T 0 : the dependence of R n on φ does not disappear even for very big values of parameter G. This is shown in Figure 3 for R n ( 1 ; 1000 ; φ ) , with n = 1 / 2 , n = 2 and n = 4 . It is important to pay attention to the different vertical scales in the different plots. Strong oscillations are observed. The mean value of these oscillations depends on index n.

2.2. A Tanh-like Profile of the Frequency

The second example is the profile f ( T ) = | tanh ( a T ) / tanh ( a ) | n with n > 0 and a = 5 , which describes a softer transition from the constant frequency to a time-dependent one. The plots in this case turn out very similar to those of the preceding section. We show one of them in Figure 4, for the value T = 1 (when the final frequency coincides with the initial one). While the lines R ( φ ) are shifted with respect to each other in the cases of f ( T ) = | T | n and f ( T ) = | tanh ( a T ) / tanh ( a ) | n , the maximal and minimal values coincide. Moreover, the average values of the final energy as functions of index n turn out practically identical for two families of frequency profiles, as shown in Figure 5. This coincidence is explained in the following sections.

3. Evolution of the Mean Oscillator Energy in the Quantum Case

In the quantum case, one has to solve the time-dependent Schrödinger equation and use the solutions to calculate various mean values, in particular, those contributing to the mean energy. However, a simpler way is to use the Ehrenfest equations for the mean values, which are immediate consequences of the Schrödinger equation. It was shown in the seminal papers by Husimi [11], Popov and Perelomov [12], Lewis and Riesenfeld [13] and Malkin, Man’ko and Trifonov [14] that the solutions of the Schrödinger and Ehrenfest equations for the harmonic oscillator with an arbitrary time-dependent frequency depend on complex functions ε ( t ) and ε * ( t ) , satisfying Equation (2) and the initial conditions
ε ( τ ) = ω 0 1 / 2 , ε ˙ ( τ ) = i ω 0 1 / 2 .
The Wronskian identity for the solutions ε ( t ) and ε * ( t ) has the form
ε ˙ ε * ε ˙ * ε = 2 i .
Then, we can write at t τ
x ( t ) = x 0 ω 0 Re [ ε ( t ) ] + p 0 ω 0 Im [ ε ( t ) ] , p ( t ) = x 0 ω 0 Re [ ε ˙ ( t ) ] + p 0 ω 0 Im [ ε ˙ ( t ) ] .
Equation (9) holds both for the classical and quantum particles (in the Heisenberg representation in the latter case). Its immediate consequences are the following formulas for the second-order moments of the canonical operators for t τ :
x 2 t = x 2 τ ω 0 Re [ ε ( t ) ] 2 + p 2 τ ω 0 Im [ ε ( t ) ] 2 + x p + p x τ Re [ ε ( t ) ] Im [ ε ( t ) ] ,
p 2 t = x 2 τ ω 0 Re [ ε ˙ ( t ) ] 2 + p 2 τ ω 0 Im [ ε ˙ ( t ) ] 2 + x p + p x τ Re [ ε ˙ ( t ) ] Im [ ε ˙ ( t ) ] .
The time-dependent mean energy is given by the formula
E ( t ) = 1 2 p 2 t + ω 2 ( t ) x 2 t .
It is worth remembering that for systems with quadratic Hamiltonians with respect to x and p, the dynamics of the first-order mean values x and p are totally independent from the dynamics of the variances σ x = x 2 x 2 , σ p = p 2 p 2 and σ x p = x p + p x / 2 x p . This means that the equations of the same form as (10) and (11) exist for the sets ( x 2 , p 2 , x p ) and ( σ x , σ p , σ x p ) .
The adiabatic (quasi-classical) approximate complex solution to Equation (2), satisfying the initial conditions (7), has the form
ε ( t ) [ ω ( t ) ] 1 / 2 e i ϕ τ ( t ) , ε ˙ ( t ) i [ ω ( t ) ] 1 / 2 e i ϕ τ ( t ) , ϕ τ ( t ) = τ t ω ( z ) d z .
Putting the solution (13) in the Equations (10)–(12), we arrive immediately at the adiabatic invariant
E ( t ) / ω ( t ) = E ( τ ) / ω 0 ,
for arbitrary initial values at t = τ . However, the solution (13) obviously loses its sense if ω ( t ) = 0 at some time instant t * (taken as t = 0 in our paper). Nonetheless, when the frequency slowly passes through a zero value and slowly becomes not too small, the conditions of the quasi-classical approximation are reestablished again. Hence, the solution for t > 0 can be written (outside some interval near t = 0 ) in the most general quasi-classical form as follows,
ε ( t ) [ ω ( t ) ] 1 / 2 u + e i ϕ ( t ) + u e i ϕ ( t ) , ε ˙ ( t ) i [ ω ( t ) ] 1 / 2 u + e i ϕ ( t ) u e i ϕ ( t ) ,
where
ϕ ( t ) = t * t ω ( z ) d z , d ϕ ( t ) / d t = ω ( t ) .
Constant complex coefficients u ± must obey the condition
| u + | 2 | u | 2 = 1 ,
which is the consequence of Equation (8). Then, Equation (12) assumes the form
E t E τ = ω ( t ) ω 0 ( β + Δ β ) ,
where
β = | u + | 2 + | u | 2 = 1 + 2 | u | 2 ,
Δ β = ω 0 2 x 2 τ p 2 τ R e u + u + ω 0 x p + p x τ I m u + u / E τ .
Equation (18) can be interpreted as a generalized adiabatic formula for the energy after the frequency passes slowly through a zero value. It shows that the quantum mechanical mean energy is proportional to the instant frequency ω ( t ) in the adiabatic regime. However, the proportionality coefficient strongly depends on the initial conditions in the most general case. This is in agreement with the classical results shown in Figure 3 and Figure 4. For this reason, we concentrate hereafter on the important special case when
p 2 τ = ω 0 2 x 2 τ , x p + p x τ = 0 .
It includes the vacuum, thermal and Fock initial quantum states. Then, Δ β = 0 . In addition, many formulas can be simplified:
x 2 t = ω 0 x 2 τ | ε ( t ) | 2 , p 2 t = p 2 τ | ε ˙ ( t ) | 2 / ω 0 ,
E t = E τ 2 ω 0 ω 2 ( t ) | ε ( t ) | 2 + | ε ˙ ( t ) | 2 .
In principle, the choice of the initial point of integration in Equation (16), defining the phase function ϕ ( t ) in Equation (15) can be arbitrary, since it influences the phases of coefficients u ± only. However, the point t * = 0 is distinguished in our problem, because ω ( t * ) = 0 . Therefore, we assume hereafter that t * = 0 in the definition of the phase ϕ ( t ) (16).
Note that after Lewis and Riesenfeld’s paper [13], many authors working on various problems related to the harmonic oscillator with a time-dependent frequency did not use as a starting point linear Equation (2) but its nonlinear analog (known under the name “Ermakov equation”)
ρ ¨ + ω 2 ( t ) ρ = ρ 3 , ρ ( t ) | ε ( t ) | ,
which follows from (2) if one writes ε = ρ exp ( i χ ) and takes into account the condition (8). Then, one can rewrite Equation (23) as follows,
E t = E τ 2 ω 0 ω 2 ( t ) ρ 2 ( t ) + [ ρ ˙ ( t ) ] 2 + [ ρ ( t ) ] 2 E τ 2 ω 0 2 | ω ( t ) | + [ ρ ˙ ( t ) ] 2 .
Consequently, the mean energy always increases when the frequency returns to its initial value, unless the time derivative ρ ˙ ( t ) is negligibly small, i.e., if u 0 (under the condition (21)). Many references on the subjects related to the Ermakov equation can be found, e.g., in the review [15] and a recent paper [16]. However, we prefer to use linear equation (2), because the key terms that help us solve the adiabatic problem are the coefficients u ± in the asymptotic formula (15).
However, how to find these constant coefficients u ± ? The numeric results of Section 2 (especially Figure 5) indicate that the answer depends on the exponent n in the form of the frequency transition through zero: ω 2 ( t ) | t | n when t 0 (assuming that ω ( t ) = 0 at t = 0 ). It is remarkable that the explicit dependence of | u | 2 on the index n can be found analytically, as shown in the next section.

4. Exact Solutions for the Power Profile of the Frequency

It is known that Equation (2) with the time-dependent frequency ω 2 ( t ) = ω 0 2 | t / τ | n can be reduced to the Bessel equation
d 2 Z d y 2 + 1 y d Z d y + 1 ν 2 y 2 Z = 0
for t > 0 (see, e.g., papers [17,18]). The same can be done for t < 0 , as soon as the initial equation is invariant with respect to the time reflection t t . One can verify that Equation (2) goes to (26) with the aid of the following transformations:
x ( t ) = | t | Z [ y ( t ) ] , ν = 1 n + 2 , y ( t ) = g t τ γ , γ = 1 2 ν , g = 2 G ν , G = ω 0 τ .
Hence, the function ε ( t ) can be written as a superposition of the Bessel functions J ν ( y ) and J ν ( y ) , although with different coefficients in the regions of t < 0 and t > 0 :
ε ( t ) = | t | × A J ν [ y ( t ) ] + B J ν [ y ( t ) ] , t < 0 A + J ν [ y ( t ) ] + B + J ν [ y ( t ) ] , t > 0 .
Constant complex coefficients A and B can be found from the initial conditions (7). Remembering that d | t | / d t = 1 for t < 0 , one obtains the following equations:
A J ν g + B J ν g = 1 / G ,
A J ν g + B J ν g = 1 G i + 1 2 G ,
where J ± ν ( z ) means the derivative of the Bessel function J ± ν ( z ) with respect to its argument z. Using the known Wronskian [19,20]
J ν ( z ) J ν ( z ) J ν ( z ) J ν ( z ) = 2 sin ( ν π ) / ( z π ) ,
we obtain the following expressions:
A = ν π G sin ( ν π ) J ν g + i + 1 2 G J ν g ,
B = ν π G sin ( ν π ) i + 1 2 G J ν g + J ν g .
Using the known identities (see, e.g., formulas 7.2 (54) and 7.2 (55) in [20])
J ν ( z ) ± z ν J ν ( z ) = z ν J ν 1 ( z ) ,
we can simplify formulas for the coefficients A and B :
A = ν π G sin ( ν π ) J 1 ν g i J ν g , B = ν π G sin ( ν π ) i J ν g + J ν 1 g .
The time derivative of function (28) at t < 0 (when d | t | / d t = 1 ) can be written with the aid of identities (29) as follows:
d ε / d t = y 2 ν | t | B J 1 ν ( y ) A J ν 1 ( y ) , t 0 .
On the other hand,
d ε / d t = y 2 ν t A + J ν 1 ( y ) B + J 1 ν ( y ) , t 0 .
Using the leading term of the Bessel function J p ( z ) = z p / [ 2 p Γ ( p + 1 ) ] at z 0 , one can see that | t | J ν ( y ) 0 when t 0 , while the product | t | J ν ( y ) tends to a finite value in this limit. Consequently, the continuity of function ε ( t ) at t = 0 implies the condition B + = B . On the other hand, y J 1 ν ( y ) / | t | 0 at t 0 , while the product y J ν 1 ( y ) / | t | tends to a finite value in this limit. Hence, the continuity of derivative d ε / d t at t = 0 can be guaranteed under the condition A + = A . Then, one can verify that the Wronskian identity (8) is satisfied identically, both for t 0 and t 0 , in view of the identity [20]
J ν ( z ) J 1 ν ( z ) + J ν ( z ) J ν 1 ( z ) = 2 sin ( ν π ) / ( z π ) .
Using Equations (28) and (30), one can write down formula (23) for the mean energy ratio R ( t ) = E ( t ) / E ( τ ) as follows:
R ( t < 0 ) R ( t > 0 ) = 1 8 g π sin ( ν π ) 2 t τ n + 1 K ( g ) K + ( y ) + K + ( g ) K ( y ) 2 K 0 ( g ) K 0 ( y ) ,
where
K + ( z ) = J ν 2 ( z ) + J ν 1 2 ( z ) , K ( z ) = J ν 2 ( z ) + J 1 ν 2 ( z ) , K 0 ( z ) = J ν 1 ( z ) J 1 ν ( z ) J ν ( z ) J ν ( z ) .
Figure 6 shows the function R ( b ) , where 1 b t / τ 1 , for ν = 1 / 4 (i.e., n = 2 and y = g b 2 ) and three values of parameter g = 0.1 , 1.0 , 10 .
Using the known formulas
J ν ( y ) ( y / 2 ) ν / Γ ( ν + 1 ) , y 1 ,
Γ ( z ) Γ ( 1 z ) = π / sin ( π z ) ,
one can see that R ( b ) is totally symmetric in the limit g 0 (an instantaneous frequency jump through a zero value), when R ( b ) = 1 + | b | n / 2 . However, the symmetry is broken for not very small values of parameter g. The known asymptotic formula for the Bessel functions of large arguments [19,20],
J ν ( z ) 2 π z cos z ν π 2 π 4 ,
results in the following simple expressions for z 1 :
K ± ( z ) = 2 π z , K 0 ( z ) = 2 cos ( ν π ) π z .
Hence, in the adiabatic limit ( g 1 and y 1 ), we obtain
R ( t < 0 ) R ( t > 0 ) = ω ( t ) 1 cos 2 ( ν π ) ω 0 sin 2 ( ν π ) .
If t < 0 , we arrive exactly at the adiabatic formula (1) for any value of the power n. On the other hand, if t > 0 (i.e., after the frequency has passed through a zero value), we see again the linear proportionality
R ( t ) = β ω ( t ) / ω 0 , β = 1 + cos 2 ( ν π ) sin 2 ( ν π ) , ν = 1 n + 2 .
The proportionality coefficient β depends on parameter n. For example, β = 5 / 3 for n = 1 , β = 3 for n = 2 and β = 7 for n = 4 , in accordance with Figure 5. If n 1 , then β 1 , while β 2 ( n / π ) 2 for n 1 . To see the limitations on the validity of the adiabatic approximation g 1 , we plot in Figure 7 the ratio ρ R ( τ ) ,
ρ = 1 4 g π sin ( ν π ) 2 K ( g ) K + ( g ) + K 0 2 ( g ) ,
as a function of g for ν = 1 / 3 , 1 / 4 , 1 / 6 . We see that the generalized adiabatic approximation (38) has an accuracy better than 1% for g > 100 .
In view of formula (35), the only nonzero contribution to the right-hand side of Equation (34) at t 0 is given by the function K + ( y ) J ν 1 2 ( y ) ( y / 2 ) ν 1 / Γ ( ν ) 2 | t | ( n + 1 ) , whereas the contributions of K ( y ) | t | 1 and K 0 ( y ) | t | 0 are eliminated by the term | t | ( n + 1 ) . Hence,
R ( t = 0 ) = π g 2 ν 1 2 ν Γ ( ν ) sin ( π ν ) 2 , g 1 .
This means, in particular, that the adiabatic formula (38) holds under the condition ω ( t ) / ω 0 g 2 ν 1 (provided g = 2 ν ω 0 τ 1 ).

5. Transition Rules for Adiabatic Coefficients after Frequency Passes through Zero Value

Equation (38) means that parameter | u | (determining the adiabatic evolution of the mean energy after the frequency passes through a zero value according to formula (18)) has the following form:
| u | = cot ( ν π ) = cot π n + 2 .
This formula can be derived directly from Equation (28). If g 1 , Equation (30) assumes the following asymptotic form:
A ν π sin ( ν π ) exp i g + ν π 2 3 π 4 , B ν π sin ( ν π ) exp i g ν π 2 + π 4
Then, Equation (28) results in the following expressions for y 1 :
ε ( t < 0 ) [ ω ( t ) ] 1 / 2 e i ( g y ) , ε ( t > 0 ) e i g [ ω ( t ) ] 1 / 2 sin ( ν π ) e i y + i cos ( ν π ) e i y .
On the other hand, calculating the phase ϕ ( t ) according to the definition (16), we find ϕ = y for t < 0 and ϕ = y for t > 0 . Hence, omitting the common phase term e i g in Equation (43) and comparing this equation with (13) and (15), we obtain the following expressions for the coefficients u ± :
u + = [ sin ( ν π ) ] 1 , u = i cot ( ν π ) .
They satisfy exactly the identity (17) and result in formula (41). Note that coefficient u + given by Equation (44) is real. However, probably, the reality of this coefficient is due to the specific exact power shape of function ω ( t ) considered in this section. For other functions ω ( t ) with a similar behavior when ω 0 , this coefficient can be complex, although with the same absolute value. An example is given in the next section. However, the formulas for the absolute values | u ± | ,
| u + | = [ sin ( ν π ) ] 1 , | u | = cot ( ν π ) ,
seem to be universal after a single frequency passage through zero.

6. Exact Solution for the Tanh Profile of the Frequency

An interesting example of exact solutions corresponds to the time-dependent frequency (a special case of the family of Epstein–Eckart profiles [21,22])
ω 2 ( t ) = ω 0 2 tanh 2 ( κ t / 2 ) , < t < , κ > 0 .
In this case, solutions to Equation (2) can be written in terms of the Gauss hypergeometric function
F ( a , b ; c ; x ) = n = 1 ( a ) n ( b ) n x n ( c ) n n ! ,
satisfying the equation
x ( 1 x ) F + ( c ( a + b + 1 ) x ) F a b F = 0 .
The first step to come to Equation (47) is to introduce the new variable ξ = tanh ( κ t / 2 ) . Then, Equation (2) takes the form
( 1 ξ 2 ) 2 d 2 x d ξ 2 2 ξ ( 1 ξ 2 ) d x d ξ + 4 ω ˜ 0 2 ξ 2 x = 0 , ω ˜ 0 ω 0 / κ .
We wish to arrive to the function F ( a , b ; c ; y ) with y = ( 1 + ξ ) / 2 . In such a case, y = 0 when t = , so that the initial condition (7) can be easily satisfied, as soon as F ( a , b ; c ; 0 ) = 1 . On the other hand, there are many relations for the function F ( a , b ; c ; 1 ) , which arises when t . Then, the asymptotics of function ε ( t ) can be easily found. Therefore, looking for the solution in the form x ( t ) = [ y ( 1 y ) ] α f ( y ) , we obtain the equation
y 2 ( 1 y ) 2 f + y ( 1 y ) ( 1 2 y ) ( 1 + 2 α ) f + ( 2 y 1 ) 2 α 2 + ω ˜ 0 2 2 α y ( 1 y ) f = 0 .
Consequently, choosing α = i ω ˜ 0 , we arrive at the solution
ε ( t ) = ω 0 1 / 2 [ y ( 1 y ) ] i ω ˜ 0 F a + , a ; c ; y = ω 0 1 / 2 [ 2 cosh ( κ t / 2 ) ] 2 i ω ˜ 0 F a + , a ; c ; y ,
where
y = 1 2 [ 1 + tanh ( κ t / 2 ) ] = 1 + e κ t 1 ,
a ± = 1 2 + 2 i ω ˜ 0 ± r , r = 1 2 1 16 ω ˜ 0 2 , c = 1 + 2 i ω ˜ 0 .
If t , function (50) goes to ε ( t ) = ω 0 1 / 2 exp ( i ω 0 t ) . If t (and y 1 ), we can use the analytic continuation of the hypergeometric function, given, e.g., by formula 2.10(1) from [20],
F ( a , b ; c ; y ) = Γ ( c ) Γ ( c a b ) Γ ( c a ) Γ ( c b ) F ( a , b ; a + b + 1 c ; 1 y ) + Γ ( c ) Γ ( a + b c ) Γ ( a ) Γ ( b ) ( 1 y ) c a b F ( c a , c b ; c + 1 a b ; 1 y ) .
Then, function (50) assumes the asymptotic form
ε ( t ) = ω 0 1 / 2 v + e i ω 0 t + v e i ω 0 t ,
v + = Γ ( 1 + 2 i ω ˜ 0 ) Γ ( 2 i ω ˜ 0 ) Γ ( 1 / 2 + 2 i ω ˜ 0 + r ) Γ ( 1 / 2 + 2 i ω ˜ 0 r ) , v = Γ ( 1 + 2 i ω ˜ 0 ) Γ ( 2 i ω ˜ 0 ) Γ ( 1 / 2 + r ) Γ ( 1 / 2 r ) .
Using relation (36), we can simplify the expression for coefficient v :
v = i cos π 1 / 4 4 ω ˜ 0 2 sinh ( 2 π ω ˜ 0 ) .
The quantity | v | 2 increases with an increase of ω ˜ 0 . If ω ˜ 0 = 1 / 4 , then | v | 2 = [ sinh ( π / 2 ) ] 2 0.19 . In the adiabatic limit ω ˜ 0 1 , we have v i coth ( 2 π ω ˜ 0 ) , i.e., v i and | v | 1 when ω ˜ 0 . In this limit, we have r = 2 i ω ˜ 0 + O ( ω ˜ 0 1 ) . Then, we can write
v + 2 i ω ˜ 0 [ Γ ( 2 i ω ˜ 0 ) ] 2 π Γ ( 1 / 2 + 4 i ω ˜ 0 ) .
Using the asymptotic Stirling formula for the gamma function,
Γ ( z ) 2 π exp ( z 1 / 2 ) ln ( z ) z , | z | 1 .
we obtain the expression
v + 2 exp 4 i ω ˜ 0 ln ( 2 ) .
Consequently, | v + | 2 | v | 2 = 1 and E ( ) / E ( ) = 3 , in accordance with formula (38).
The asymptotic form (54) is similar to the general adiabatic solution (15). Using the definition (16) of the phase ϕ ( t ) , we obtain the formula (remember that ω ( t ) = | tanh ( κ t / 2 ) | in this section, as soon as function ω ( t ) is assumed to be non-negative in formula (15))
ϕ ( t ) = 2 ω ˜ 0 ln [ cosh ( κ t / 2 ) ] sign ( t ) .
If t ± , then ϕ ( t ) ω 0 t 2 ω ˜ 0 ln ( 2 ) sign ( t ) . This means that, according to (15), the function ω 0 1 / 2 exp i ω 0 t + 2 ω ˜ 0 ln ( 2 ) at t goes to the following superposition at t :
ω 0 1 / 2 u + exp i ω 0 t 2 ω ˜ 0 ln ( 2 ) + u exp i ω 0 t + 2 ω ˜ 0 ln ( 2 ) .
Comparing this expression with (54), we conclude that
u + = v + = 2 exp 4 i ω ˜ 0 ln ( 2 ) , u = v = i .
We see that the phases of complex coefficients u ± are sensitive to the rate of the adiabatic evolution through the term ω ˜ 0 . A strong consequence of this result is considered in the next section.

7. Double Adiabatic Passage of Frequency through Zero Value

What can happen if the frequency passes again through a zero value? Then, we have to make the transformation of function (15), using the superposition principle and two additional observations. First, the function ω 1 / 2 exp ( i ϕ ) transforms as function ε * ( t ) after the frequency passes through a zero value. Second, applying the transformation rule (15) to the second transition, we must use the phase ϕ ˜ ( t ) , where the integral over the frequency is taken from the second transition point t * * . Obviously,
ϕ ( t ) = ϕ ˜ ( t ) + Φ , Φ ( t * , t * * ) = t * t * * ω ( z ) d z .
We suppose that the transition rule through the second zero has the form
[ ω ( t ) ] 1 / 2 e i ϕ ˜ ( t ) t < t * * [ ω ( t ) ] 1 / 2 w + e i ϕ ˜ ( t ) + w e i ϕ ˜ ( t ) t > t * * , | w + | 2 | w | 2 = 1 .
Comparing the two forms of the solution ε ( t ) for t > t * * (far enough from that point), we arrive at the equality (omitting the common term ω 1 / 2 and using the notation U ± for the coefficients at t > t * * )
U + e i ϕ ˜ ( t ) + U e i ϕ ˜ ( t ) = u + e i Φ w + e i ϕ ˜ ( t ) + w e i ϕ ˜ ( t ) + u e i Φ w + * e i ϕ ˜ ( t ) + w * e i ϕ ˜ ( t ) .
Hence,
U + = w + u + e i Φ + w * u e i Φ , U = w u + e i Φ + w + * u e i Φ .
One can verify that the identity | U + | 2 | U | 2 = 1 is fulfilled exactly. The adiabatic mean energy amplification factor after the second passage through a zero frequency equals
β = 1 + 2 | U | 2 = 1 + 2 | w | 2 | u + | 2 + | w + | 2 | u | 2 + 2 R e w + w u + u * e 2 i Φ .
In the adiabatic regime, Φ 1 . Moreover, this phase is very sensitive to the form of function ω ( t ) and the distance between the zero-point instances t * and t * * . In addition, coefficients u ± and w ± can be strongly phase-sensitive, as shown in Section 6. This means that it is practically impossible to predict the energy mean value after twice crossing the zero frequency (quite differently from the single crossing). The extremal values of β are as follows,
β m i n = 1 + 2 | w + u | | w u + | 2 , β m a x = 1 + 2 | w + u | + | w u + | 2 .
In particular, if w ± = u ± , then β m i n = 1 , meaning that, in principle, the mean energy can return to the initial value after the frequency passes through a zero value two times. On the other hand, β m a x = 1 + 8 | u + u | 2 under the same conditions. If n = 2 , β m a x = 17 .

Classical Illustrations

To illustrate the effects after single and double frequency crossings through a zero value, we considered the classical motion with the frequency ω 2 ( t ) = ω 0 2 sin 2 [ π t / ( 2 τ ) ] | for t = τ and the initial conditions (5). Figure 8 shows the energy ratio R (6) at the instants T = t / τ = 1 and T = 3 , for several values of parameter G close to G = 1000 . If T = 1 (single crossing), the variations of parameter G result in shifts of the curves without changing the maximal, minimal and average values. On the other hand, the picture is totally different for T = 3 (double crossing). In this case, Φ = 4 G / π . Hence, the variation Δ Φ = 2 π , when one can expect a similar behavior, corresponds to Δ G = π 2 / 2 5 . On the other hand, a twice smaller variation Δ G 2 , yielding Δ Φ π , results in a totally different behavior, as one can see in the figure for T = 3 .

8. Energy Fluctuations

Figure 3, Figure 4 and Figure 8 show strong energy fluctuations (as functions of the initial phase) after the frequency passes through a zero value. These fluctuations can be characterized by the variance σ E = E 2 E 2 . Using the solutions (9) of the Heisenberg equations of motion, one can write σ E in terms of the fourth- and second-order moments of the canonical variables x and p and various products of functions ε ( t ) , ε ˙ ( t ) and their complex-conjugated partners. The complete formula is rather cumbersome in the most general case. For this reason, we consider here the simplest case of the initial Fock quantum state | N . This special case is probably the most interesting, because the famous adiabatic theorem in quantum mechanics was proven by Born and Fock [23] exactly for the Fock states. In this special case (as well as for arbitrary diagonal mixtures of Fock states), the nonzero statistical moments are those containing even powers of each variable, x or p. After some algebra, one can obtain the following formula (using the dimensionless variables, assuming = m = ω 0 = 1 , so that x 4 = p 4 ):
16 E 2 t = 2 x 4 τ A 2 + B 2 + x 2 p 2 + p 2 x 2 τ A 2 B 2 + ( x p + p x ) 2 τ C 2 ,
where
A ( t ) = ω 2 ( t ) | ε ( t ) | 2 + | ε ˙ ( t ) | 2 , B ( t ) = Re ω 2 ( t ) ε 2 ( t ) + ε ˙ 2 ( t ) , C ( t ) = I m ω 2 ( t ) ε 2 ( t ) + ε ˙ 2 ( t ) .
In the adiabatic regime (15), we have
A = 2 ω ( t ) u + | 2 + | u | 2 , B = 4 ω ( t ) Re u + u , C = 4 ω ( t ) I m u + u .
For the initial Fock state | N , we have
x 4 τ = 3 4 2 N 2 + 2 N + 1 , x 2 p 2 + p 2 x 2 τ = 1 2 2 N 2 + 2 N 1 ,
( x p + p x ) 2 τ = 2 N 2 + N + 1 .
Hence,
E 2 t / ω 2 ( t ) = u + | 2 + | u | 2 2 ( N + 1 / 2 ) 2 + 2 | u + u | 2 N 2 + N + 1 .
Remembering that the mean energy equals E t = ω ( t ) u + | 2 + | u | 2 ( N + 1 / 2 ) , we arrive at the unexpectedly simple formula for the energy variance:
σ E ( t ) = 2 ω 2 ( t ) | u + u | 2 N 2 + N + 1 , σ E ( t ) E t 2 = 2 | u + u | 2 N 2 + N + 1 N 2 + N + 1 / 4 .
In the absence of zero-frequency values we have u = 0 . In this case, σ E ( t ) 0 , in accordance with the Born–Fock theorem. However, this theorem is broken when the frequency passes through a zero value. For example, for the initial vacuum state ( N = 0 ) and the power index n = 2 of the single frequency transition through a zero value, we obtain σ E ( t ) / E t 2 = 16 . This ratio can be four times smaller if N 1 .
In quantum optics, fluctuations are frequently characterized by the Mandel factor [24]
Q = n ^ 2 n ^ 2 / n ^ 1 , n ^ E ^ / ( ω ) 1 / 2 .
Then, for the initial Fock state | N (having Q = 1 , which means the so called sub-Poissonian statistics), we obtain the following instantaneous values:
n ^ = N + 2 | u | 2 ( N + 1 / 2 ) , n ^ 2 n ^ 2 = σ E / ( ω ) 2 .
Consequently,
Q = | u | 2 1 + 2 | u | 2 + N 2 | u | 4 1 + 2 N 2 | u + u | 2 N + 2 | u | 2 ( N + 1 / 2 ) .
In particular, for | u | 2 = 1 , we have
Q = 3 + N + 4 N 2 3 N + 1 .
This means that the statistics become super-Poissonian after the frequency passage through a zero value. However, the super-Poissonianity is not very strong, because Q ( 4 / 3 ) N ( 4 / 9 ) n ^ for N 1 , whereas Q = n ^ for the “strongly super-Poissonian” thermal states.

9. Evolution of the Fock States

What happens with the initial Fock state | N when the frequency passes through a zero value? Obviously, it cannot survive, as soon as the mean energy and especially the energy variance increase substantially. This means that the initial Fock state becomes a superposition of many Fock states. However, what is the width of the new distribution? Is it concentrated near some distinguished states, or it is very wide and almost uniform, especially when N 1 ? To answer these questions, one can use general results concerning the quantum harmonic oscillator with a time-dependent frequency [11,12,13,14] (the details can be found, e.g., in the review [25]). Recall that the Fock states | N are eigenstates of the operator a ^ a ^ , where a ^ and a ^ are standard annihilation and creation operators. When the frequency varies with time, operators a ^ and a ^ become new operators, A ^ and A ^ , which are quantum integrals of motion. The time-dependent state | N t remains the eigenstate of operator A ^ A ^ . As soon as operators A ^ and A ^ maintain their linear form with respect to operators x ^ and p ^ , the wave function of state | N t maintains its functional form as the product of some Gaussian exponential by the Hermite polynomial. The explicit form, found in [12,14] (see also [18]), is as follows (in dimensionless units with = m = 1 ),
x | N t = N ! ε π 1 / 2 ε * 2 ε N / 2 exp i ε ˙ 2 ε x 2 H N x | ε | ,
where ε ( t ) is the solution to Equation (2) satisfying conditions (7) and (8). Transition probabilities | M | N t | 2 between the instantaneous Fock state | M (when ε ( t ) = ω 1 / 2 exp ( i ω t ) ) and exact time-dependent state | N t were calculated in different forms in [11,12,13,14]. In the generalized adiabatic regime (15), the results of [12,14] can be written in the form (symmetric with respect to M and N)
| M | N t | 2 = N < ! N > ! | u + | 1 P ( M + N ) / 2 | M N | / 2 | u + | 1 2 ,
where P j k ( z ) is the associated Legendre polynomial, N < = m i n ( M , N ) , N > = m a x ( M , N ) . Formula (72) holds provided | M N | / 2 is an integer; otherwise, the probability equals zero. Note that the probabilities do not depend on time, as soon as the adiabatic solution (15) is valid.
In some cases, it can be convenient to use the expression of the associated Legendre polynomials in terms of the Gauss hypergeometric function,
P n m ( x ) = ( 1 ) m ( n + m ) ! 2 m ( n m ) ! m ! ( 1 x 2 ) m / 2 F m n , m + n + 1 ; m + 1 ; 1 x 2 .
Then,
| M | N t | 2 = 2 ( N > ) ! | u | | M N | ( N < ) ! | M N | 2 ! 2 | 2 u + | | M N | + 1 F N < , N > + 1 ; | M N | 2 + 1 ; | u + | 1 2 | u + | 2 .
In the case of a single frequency crossing through a zero value, this formula can also be written as
| M | N t | 2 = ( N > ) ! sin ( ν π ) [ cos ( ν π ) ] | M N | ( N < ) ! | M N | 2 ! 2 2 | M N | F N < , N > + 1 ; | M N | 2 + 1 ; sin 2 [ ( 1 2 ν ) π / 4 ] 2 .
Among different special cases, we bring here two formulas:
| 2 K | 0 t | 2 = ( 2 K 1 ) ! ! | u | 2 K ( 2 K ) ! ! | u + | 2 K + 1 = ( 2 K 1 ) ! ! ( 2 K ) ! ! sin ( ν π ) [ cos ( ν π ) ] 2 K ,
| N | N t | 2 = | u + | 1 P N | u + | 1 2 = sin ( ν π ) P N sin ( ν π ) 2 ,
where P N ( z ) P N 0 ( z ) is the usual Legendre polynomial. The first equalities in Equations (76) and (77) hold in the most general adiabatic case (including multiple frequency crossings through a zero value), whereas the second equalities are valid for the single crossing. The distribution (76) (describing the evolution of the initial ground state) decreases monotonously as function of parameter K. However, the situation is totally different for other initial Fock states, especially when N 1 .
The survival probabilities p s ( N ) | N | N t | 2 rapidly diminish with the quantum number N. For example, if n = 2 ( ν = 1 / 4 ), we find the following surviving probabilities after the single frequency zero-crossing (when the mean energy triplicates):
p s ( 0 ) = 1 2 , p s ( 1 ) = 1 2 2 , p s ( 2 ) = 1 16 2 , p s ( 3 ) = 1 32 2 ,
This means that the initial Fock state | N becomes a superposition of a large number of different Fock states | M , see Figure 9. It is impressive that the probability of transition N 3 N is very small, whereas the probability p N ( M 3 N ) is about 50%. In addition, the distribution of probabilities p N ( M ) with M < 3 N looks rather irregular, whereas some regular picture is observed for M > 3 N . Unfortunately, we did not succeed to find an analytic approximation for this regular picture.

10. Squeezing Evolution

If the frequency ω does not depend on time, the evolution of the coordinate variance is given by the formula
σ x ( t ) = σ x ( 0 ) cos 2 ( ω t ) + σ p ( 0 ) ( m ω ) 2 sin 2 ( ω t ) + σ x p ( 0 ) ( m ω ) sin ( 2 ω t ) .
Minimizing this expression over time, one can write the minimal value σ m i n as (similar formulas were obtained, e.g., in [26,27,28])
σ m i n = ( m ω 2 ) 1 E E 2 ω 2 D ,
where
E = σ p / ( 2 m ) + m ω 2 σ x / 2 , D = σ x σ p σ x p 2 .
The quantity D is the simplest example of quantum universal invariants [29], which do not depend on time (although they depend on the initial state) for arbitrary quadratic Hamiltonians. On the other hand, D 2 / 4 for any (normalizable) quantum state due to the Schrödinger–Robertson uncertainty relation. The energy of quantum fluctuations E satisfies the inequality E ω / 2 . Therefore, it is convenient to use two dimensionless parameters, λ 1 and γ 1 , according to the relations E = λ ω / 2 and D = γ 2 2 / 4 . Then, normalizing the minimal value σ m i n by the variance in the vacuum state σ v a c = / ( 2 m ω ) , one can obtain the following formula for the invariant squeezing coefficient s = σ m i n / σ v a c :
s = λ λ 2 γ 2 = γ 2 λ + λ 2 γ 2 .
For the states satisfying the initial conditions (21), we have γ = λ and s = λ 1 . Moreover, parameter λ maintains it initial value in the standard adiabatic case (14). However, if the frequency passes through a zero value, in the new adiabatic regime (18), λ goes to β λ , while γ maintains its initial value. Hence, the new squeezing coefficient equals
s = λ β + β 2 1 .
Hence, the initial vacuum state becomes squeezed when the frequency passes adiabatically through a zero value. Using Equation (38), we obtain the following value of the squeezing coefficient after a single passage through zero:
s = λ tan 2 ( ν π / 2 ) .
In particular, s 0.17 λ for ν = 1 / 4 (i.e., n = 2 and β = 3 ), so the Fock states | N become squeezed for N 2 (when λ 5 ) in this special case.

11. Conclusions

The first main result of the paper is the discovery of the existence of the generalized adiabatic invariant in the form of Equation (18). In the most general case, the adiabatic proportionality coefficient in this equation depends on the initial state. This dependence disappears after averaging over parameters of families of initial states with the same energy (in particular, such averaging happens automatically for the initial vacuum, Fock and thermal states). Then, universal relations (45) exist, provided the frequency passes through zero only once. In the cases of multiple frequency passages through zero, the energy adiabatic coefficients become sensitive to the additional parameter, the phase Φ , according to Equation (63). As a consequence, the adiabatic behavior after many crossings through a zero frequency value can be quasi-chaotic. Under specific conditions, the mean energy can return to the initial value after a double frequency passage through zero. The original Born–Fock adiabatic theorem is broken after the frequency passes through a zero value. Although the functional shape of the wave function of the initial Fock state is preserved in the form of the product of a Gaussian exponential by the Hermite polynomial, the arguments of this form are not determined totally by the instantaneous frequency. However, the probability distribution over the instantaneous Fock states, determined by the adiabatic coefficients | u ± | , according to Equation (72), does not depend on time, as soon as the adiabatic regime is justified. This statement can be considered as the generalized Born–Fock theorem; it is the second main result of the paper. Note that the time-independent probability distributions can be different after each frequency passage through zero.
In view of the mean energy amplification (e.g., triplication in the most natural case of a linear frequency dependence near the zero point), any initial state becomes significantly deformed. For example, coherent states (which possess the same quadrature variances as the vacuum state) will be transformed into squeezed states. The same can be said about initial thermal states: they will become Gaussian mixed states with unequal quadrature variances (and squeezed under certain conditions), maintaining the initial value of the quantum purity.

Author Contributions

V.V.D.: conceptualization, methodology, analytical calculations and writing the paper. A.V.D.: numeric calculations, plotting figures and writing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the partial support of the Brazilian funding agency Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Landau, L.D.; Lifschitz, E.M. Mechanics; Mir: Moscow, Russia, 1966. [Google Scholar]
  2. Kulsrud, R.M. Adiabatic invariant of the harmonic oscillator. Phys. Rev. 1957, 106, 205–207. [Google Scholar] [CrossRef]
  3. Knorr, G.; Pfirsch, D. The variation of the adiabatic invariant of the harmonic oscillator. Z. Naturforschung A 1966, 21, 688–693. [Google Scholar] [CrossRef]
  4. Solimeno, S.; Di Porto, P.; Crosignani, B. Quantum harmonic oscillator with time-dependent frequency. J. Math. Phys. 1969, 10, 1922–1928. [Google Scholar] [CrossRef]
  5. 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]
  6. Keller, J.B.; Mu, Y. Changes in adiabatic invariants. Ann. Phys. 1991, 205, 219–227. [Google Scholar] [CrossRef]
  7. Robnik, M.; Romanovski, V.G. Exact analysis of adiabatic invariants in the time-dependent harmonic oscillator. J. Phys. A Math. Gen. 2006, 39, L35–L41. [Google Scholar] [CrossRef] [Green Version]
  8. Robnik, M.; Romanovski, V.G. Energy evolution in time-dependent harmonic oscillator. Open Syst. Inf. Dyn. 2006, 13, 197–222. [Google Scholar] [CrossRef] [Green Version]
  9. Zaugg, T.; Meystre, P.; Lenz, G.; Wilkens, M. Theory of adiabatic cooling in cavities. Phys. Rev. A 1994, 49, 3011–3021. [Google Scholar] [CrossRef]
  10. Möller, K.B.; Henriksen, N.E. On wave-packet dynamics in a decaying quadratic potential. Phys. Scripta 1997, 55, 542–546. [Google Scholar] [CrossRef]
  11. Husimi, K. Miscellanea in elementary quantum mechanics. II. Prog. Theor. Phys. 1953, 9, 381–402. [Google Scholar] [CrossRef]
  12. Popov, V.S.; Perelomov, A.M. Parametric excitation of a quantum oscillator. Zhurn. Eksp. Teor. Fiz. 1969, 56, 1375–1390, [Sov. Phys. JETP 1969, 29, 738–745]. [Google Scholar]
  13. 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]
  14. Malkin, I.A.; Man’ko, V.I.; Trifonov, D.A. Coherent states and transition probabilities in a time-dependent electromagnetic field. Phys. Rev. D 1970, 2, 1371–1385. [Google Scholar] [CrossRef]
  15. Cordero-Soto, R.; Suazo, E.; Suslov, S.K. Quantum integrals of motion for variable quadratic Hamiltonians. Ann. Phys. 2010, 325, 1884–1912. [Google Scholar] [CrossRef] [Green Version]
  16. Ramos-Prieto, I.; Urzúa-Pineda, A.R.; Soto-Eguibar, F.; Moya-Cessa, H.M. KvN mechanics approach to the time-dependent frequency harmonic oscillator. Sci. Rep. 2018, 8, 8401. [Google Scholar] [CrossRef] [Green Version]
  17. 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]
  18. Kim, S.P. A class of exactly solved time-dependent quantum harmonic oscillators. J. Phys. A Math. Gen. 1994, 27, 3927–3936. [Google Scholar] [CrossRef]
  19. Gradshtein, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic: Amsterdam, The Netherlands, 2007. [Google Scholar]
  20. Erdélyi, A. (Ed.) Bateman Manuscript Project: Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  21. Eckart, C. The penetration of a potential barrier by electrons. Phys. Rev. 1930, 35, 1303–1309. [Google Scholar] [CrossRef]
  22. Epstein, P.S. Reflection of waves in an inhomogeneous absorbing medium. Proc. Nat. Acad. Sci. USA 1930, 16, 627–637. [Google Scholar] [CrossRef] [Green Version]
  23. Born, M.; Fock, V. Beweis des Adiabatensatzes. Z. Phys. 1928, 51, 165–180. [Google Scholar] [CrossRef]
  24. Mandel, L. Sub-Poissonian photon statistics in resonance fluorescence. Opt. Lett. 1979, 4, 205–207. [Google Scholar] [CrossRef] [PubMed]
  25. Dodonov, V.V.; Man’ko, V.I. Invariants and correlated states of nonstationary quantum systems. In Invariants and the Evolution of Nonstationary Quantum Systems; Proceedings of Lebedev Physics, Institute; Markov, M.A., Ed.; Nova Science: Commack, NY, USA, 1989; Volume 183, pp. 103–261. [Google Scholar]
  26. Lukš, A.; Peřinová, V.; Peřina, J. Principal squeezing of vacuum fluctuations. Opt. Commun. 1988, 67, 149–151. [Google Scholar] [CrossRef]
  27. Lukš, A.; Peřinová, V.; Hradil, Z. Principal squeezing. Acta Phys. Polon. A 1988, 74, 713–721. [Google Scholar]
  28. Dodonov, V.V.; Man’ko, V.I.; Polynkin, P.G. Geometrical squeezed states of a charged particle in a time-dependent magnetic field. Phys. Lett. A 1994, 188, 232–238. [Google Scholar] [CrossRef]
  29. Dodonov, V.V. Universal integrals of motion and universal invariants of quantum systems. J. Phys. A Math. Gen. 2000, 33, 7721–7738. [Google Scholar] [CrossRef]
Figure 1. The dimensionless energy of a classical particle at the dimensionless instant T = 1 / 2 in the nonadiabatic regime ( G = 1 ).
Figure 1. The dimensionless energy of a classical particle at the dimensionless instant T = 1 / 2 in the nonadiabatic regime ( G = 1 ).
Entropy 25 00002 g001
Figure 2. The dimensionless energy of a classical particle at the dimensionless instant T = 1 / 2 in the adiabatic regime ( G = 1000 ).
Figure 2. The dimensionless energy of a classical particle at the dimensionless instant T = 1 / 2 in the adiabatic regime ( G = 1000 ).
Entropy 25 00002 g002
Figure 3. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ), for the profile f ( T ) = | T | n .
Figure 3. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ), for the profile f ( T ) = | T | n .
Entropy 25 00002 g003
Figure 4. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ) for the profile f ( T ) = | tanh ( a T ) / tanh ( a ) | n a = 5 .
Figure 4. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ) for the profile f ( T ) = | tanh ( a T ) / tanh ( a ) | n a = 5 .
Entropy 25 00002 g004
Figure 5. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ), averaged over the initial phase φ , as function of index n. Left: for the profile f ( T ) = | T | n . Right: for the profile f ( T ) = | tanh ( a T ) / tanh ( a ) | n , a = 5 .
Figure 5. The dimensionless energy of a classical particle at the dimensionless instant T = 1 in the case of slow evolution ( G = 1000 ), averaged over the initial phase φ , as function of index n. Left: for the profile f ( T ) = | T | n . Right: for the profile f ( T ) = | tanh ( a T ) / tanh ( a ) | n , a = 5 .
Entropy 25 00002 g005
Figure 6. The function R ( b ) for g = 0.1 , 1.0 , 10 and ν = 1 / 4 .
Figure 6. The function R ( b ) for g = 0.1 , 1.0 , 10 and ν = 1 / 4 .
Entropy 25 00002 g006
Figure 7. The function ρ ( g ) (39) for g > 1 and ν = 1 / 3 , 1 / 4 , 1 / 6 .
Figure 7. The function ρ ( g ) (39) for g > 1 and ν = 1 / 3 , 1 / 4 , 1 / 6 .
Entropy 25 00002 g007
Figure 8. The dimensionless energy of a classical particle at the dimensionless instants T = 1 (top) and T = 3 (bottom) for several different values of the adiabatic parameter G = ω 0 τ , shown near the respective lines. The frequency profile is ω 2 ( t ) = ω 0 2 sin 2 [ π t / ( 2 τ ) ] | .
Figure 8. The dimensionless energy of a classical particle at the dimensionless instants T = 1 (top) and T = 3 (bottom) for several different values of the adiabatic parameter G = ω 0 τ , shown near the respective lines. The frequency profile is ω 2 ( t ) = ω 0 2 sin 2 [ π t / ( 2 τ ) ] | .
Entropy 25 00002 g008
Figure 9. The probability p ( M ) (given by Equation (74)) of finding the initial Fock state | N in the Fock state | M after the frequency slowly passes through a zero value, in the case of | u | = 1 and | u + | = 2 .
Figure 9. The probability p ( M ) (given by Equation (74)) of finding the initial Fock state | N in the Fock state | M after the frequency slowly passes through a zero value, in the case of | u | = 1 and | u + | = 2 .
Entropy 25 00002 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dodonov, V.V.; Dodonov, A.V. Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero. Entropy 2023, 25, 2. https://doi.org/10.3390/e25010002

AMA Style

Dodonov VV, Dodonov AV. Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero. Entropy. 2023; 25(1):2. https://doi.org/10.3390/e25010002

Chicago/Turabian Style

Dodonov, Viktor V., and Alexandre V. Dodonov. 2023. "Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero" Entropy 25, no. 1: 2. https://doi.org/10.3390/e25010002

APA Style

Dodonov, V. V., & Dodonov, A. V. (2023). Adiabatic Amplification of the Harmonic Oscillator Energy When the Frequency Passes through Zero. Entropy, 25(1), 2. https://doi.org/10.3390/e25010002

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