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
to its time-dependent frequency
[
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:
and
. 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
We consider a special case when the time-dependent frequency can be written as
, where
is a non-negative function with the properties
We assume that
for
. Since we are interested mainly in the evolution of the energy at
, we consider a one-parameter family of classical trajectories with the same initial energy
. Then, assuming the particle mass
, it is convenient to use the initial coordinate
and initial velocity
in the form
To solve Equation (
2) numerically, we introduce dimensionless variables
and
, arriving at the equation for
with the initial conditions
Then, the dimensionless energy ratio
depends on the dimensionless time
T and two parameters,
G and
:
where the derivative
must be taken at instant
T. The existence of the adiabatic invariant (
1) implies that
for
and big enough values of parameter
G,
independently on the values of parameter φ. In the following subsections, we study what can happen if
, for different families of functions
.
2.1. A Power Profile of the Frequency
Our first example is the power profile
with
. We solved Equation (
4) numerically and calculated the dimensionless ratio
The adiabatic ratio must equal
. In
Figure 1 and
Figure 2, we show
as function of
for the fixed values
,
,
and
, comparing different behaviors when
(no adiabaticity) and
.
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
is, indeed,
.
However, the situation becomes quite different for
: the dependence of
on
does not disappear even for very big values of parameter
G. This is shown in
Figure 3 for
, with
,
and
. 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
with
and
, 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
(when the final frequency coincides with the initial one). While the lines
are shifted with respect to each other in the cases of
and
, 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
and
, satisfying Equation (
2) and the initial conditions
The Wronskian identity for the solutions
and
has the form
Then, we can write at
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
:
The time-dependent mean energy is given by the formula
It is worth remembering that for systems with
quadratic Hamiltonians with respect to
x and
p, the dynamics of the first-order mean values
and
are
totally independent from the dynamics of the variances
,
and
. This means that the equations of the same form as (
10) and (
11) exist for the sets
and
.
The adiabatic (quasi-classical) approximate complex solution to Equation (
2), satisfying the initial conditions (
7), has the form
Putting the solution (
13) in the Equations (
10)–(
12), we arrive immediately at the adiabatic invariant
for arbitrary initial values at
. However, the solution (
13) obviously loses its sense if
at some time instant
(taken as
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
can be written (outside some interval near
) in the most general quasi-classical form as follows,
where
Constant complex coefficients
must obey the condition
which is the consequence of Equation (
8). Then, Equation (
12) assumes the form
where
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
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
It includes the vacuum, thermal and Fock initial quantum states. Then,
. In addition, many formulas can be simplified:
In principle, the choice of the initial point of integration in Equation (
16), defining the phase function
in Equation (
15) can be arbitrary, since it influences the phases of coefficients
only. However, the point
is distinguished in our problem, because
. Therefore, we assume hereafter that
in the definition of the phase
(
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”)
which follows from (
2) if one writes
and takes into account the condition (
8). Then, one can rewrite Equation (
23) as follows,
Consequently, the mean energy always
increases when the frequency returns to its initial value, unless the time derivative
is negligibly small, i.e., if
(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
in the asymptotic formula (
15).
However, how to find these constant coefficients
? 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:
when
(assuming that
at
). It is remarkable that the explicit dependence of
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
can be reduced to the Bessel equation
for
(see, e.g., papers [
17,
18]). The same can be done for
, as soon as the initial equation is invariant with respect to the time reflection
. One can verify that Equation (
2) goes to (
26) with the aid of the following transformations:
Hence, the function
can be written as a superposition of the Bessel functions
and
, although with different coefficients in the regions of
and
:
Constant complex coefficients
and
can be found from the initial conditions (
7). Remembering that
for
, one obtains the following equations:
where
means the derivative of the Bessel function
with respect to its argument
z. Using the known Wronskian [
19,
20]
we obtain the following expressions:
Using the known identities (see, e.g., formulas 7.2 (54) and 7.2 (55) in [
20])
we can simplify formulas for the coefficients
and
:
The time derivative of function (
28) at
(when
) can be written with the aid of identities (
29) as follows:
On the other hand,
Using the leading term of the Bessel function
at
, one can see that
when
, while the product
tends to a finite value in this limit. Consequently, the continuity of function
at
implies the condition
. On the other hand,
at
, while the product
tends to a finite value in this limit. Hence, the continuity of derivative
at
can be guaranteed under the condition
. Then, one can verify that the Wronskian identity (
8) is satisfied identically, both for
and
, in view of the identity [
20]
Using Equations (
28) and (
30), one can write down formula (
23) for the mean energy ratio
as follows:
where
Figure 6 shows the function
, where
, for
(i.e.,
and
) and three values of parameter
.
Using the known formulas
one can see that
is totally symmetric in the limit
(an instantaneous frequency jump through a zero value), when
. 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],
results in the following simple expressions for
:
Hence, in the adiabatic limit (
and
), we obtain
If
, we arrive exactly at the adiabatic formula (
1) for any value of the power
n. On the other hand, if
(i.e., after the frequency has passed through a zero value), we see again the linear proportionality
The proportionality coefficient
depends on parameter
n. For example,
for
,
for
and
for
, in accordance with
Figure 5. If
, then
, while
for
. To see the limitations on the validity of the adiabatic approximation
, we plot in
Figure 7 the ratio
,
as a function of
g for
. We see that the generalized adiabatic approximation (
38) has an accuracy better than 1% for
.
In view of formula (
35), the only nonzero contribution to the right-hand side of Equation (
34) at
is given by the function
, whereas the contributions of
and
are eliminated by the term
. Hence,
This means, in particular, that the adiabatic formula (
38) holds under the condition
(provided
).
5. Transition Rules for Adiabatic Coefficients after Frequency Passes through Zero Value
Equation (
38) means that parameter
(determining the adiabatic evolution of the mean energy after the frequency passes through a zero value according to formula (
18)) has the following form:
This formula can be derived directly from Equation (
28). If
, Equation (
30) assumes the following asymptotic form:
Then, Equation (
28) results in the following expressions for
:
On the other hand, calculating the phase
according to the definition (
16), we find
for
and
for
. Hence, omitting the common phase term
in Equation (
43) and comparing this equation with (
13) and (
15), we obtain the following expressions for the coefficients
:
They satisfy exactly the identity (
17) and result in formula (
41). Note that coefficient
given by Equation (
44) is real. However, probably, the reality of this coefficient is due to the specific exact power shape of function
considered in this section. For other functions
with a similar behavior when
, 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
,
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])
In this case, solutions to Equation (
2) can be written in terms of the Gauss hypergeometric function
satisfying the equation
The first step to come to Equation (
47) is to introduce the new variable
. Then, Equation (
2) takes the form
We wish to arrive to the function
with
. In such a case,
when
, so that the initial condition (
7) can be easily satisfied, as soon as
. On the other hand, there are many relations for the function
, which arises when
. Then, the asymptotics of function
can be easily found. Therefore, looking for the solution in the form
, we obtain the equation
Consequently, choosing
, we arrive at the solution
where
If
, function (
50) goes to
. If
(and
), we can use the analytic continuation of the hypergeometric function, given, e.g., by formula 2.10(1) from [
20],
Then, function (
50) assumes the asymptotic form
Using relation (
36), we can simplify the expression for coefficient
:
The quantity
increases with an increase of
. If
, then
. In the adiabatic limit
, we have
, i.e.,
and
when
. In this limit, we have
. Then, we can write
Using the asymptotic Stirling formula for the gamma function,
we obtain the expression
Consequently,
and
, in accordance with formula (
38).
The asymptotic form (
54) is similar to the general adiabatic solution (
15). Using the definition (
16) of the phase
, we obtain the formula (remember that
in this section, as soon as function
is assumed to be non-negative in formula (
15))
If
, then
. This means that, according to (
15), the function
at
goes to the following superposition at
:
Comparing this expression with (
54), we conclude that
We see that the phases of complex coefficients
are sensitive to the rate of the adiabatic evolution through the term
. 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
transforms as function
after the frequency passes through a zero value. Second, applying the transformation rule (
15) to the second transition, we must use the phase
, where the integral over the frequency is taken from the second transition point
. Obviously,
We suppose that the transition rule through the second zero has the form
Comparing the two forms of the solution
for
(far enough from that point), we arrive at the equality (omitting the common term
and using the notation
for the coefficients at
)
Hence,
One can verify that the identity
is fulfilled exactly. The adiabatic mean energy amplification factor after the second passage through a zero frequency equals
In the adiabatic regime,
. Moreover, this phase is very sensitive to the form of function
and the distance between the zero-point instances
and
. In addition, coefficients
and
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,
In particular, if
, then
, 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,
under the same conditions. If
,
.
Classical Illustrations
To illustrate the effects after single and double frequency crossings through a zero value, we considered the classical motion with the frequency
for
and the initial conditions (
5).
Figure 8 shows the energy ratio
R (
6) at the instants
and
, for several values of parameter
G close to
. If
(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
(double crossing). In this case,
. Hence, the variation
, when one can expect a similar behavior, corresponds to
. On the other hand, a twice smaller variation
, yielding
, results in a totally different behavior, as one can see in the figure for
.
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
. Using the solutions (
9) of the Heisenberg equations of motion, one can write
in terms of the fourth- and second-order moments of the canonical variables
x and
p and various products of functions
,
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
. 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
, so that
):
where
In the adiabatic regime (
15), we have
For the initial Fock state
, we have
Hence,
Remembering that the mean energy equals
, we arrive at the unexpectedly simple formula for the energy variance:
In the absence of zero-frequency values we have
. In this case,
, 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 (
) and the power index
of the single frequency transition through a zero value, we obtain
. This ratio can be four times smaller if
.
In quantum optics, fluctuations are frequently characterized by the
Mandel factor [
24]
Then, for the initial Fock state
(having
, which means the so called sub-Poissonian statistics), we obtain the following instantaneous values:
Consequently,
In particular, for
, we have
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
for
, whereas
for the “strongly super-Poissonian” thermal states.
9. Evolution of the Fock States
What happens with the initial Fock state
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
? 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
are eigenstates of the operator
, where
and
are standard annihilation and creation operators. When the frequency varies with time, operators
and
become new operators,
and
, which are quantum integrals of motion. The time-dependent state
remains the eigenstate of operator
. As soon as operators
and
maintain their
linear form with respect to operators
and
, the wave function of state
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
),
where
is the solution to Equation (
2) satisfying conditions (
7) and (
8). Transition probabilities
between the instantaneous Fock state
(when
) and exact time-dependent state
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)
where
is the associated Legendre polynomial,
,
. Formula (
72) holds provided
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,
Then,
In the case of a single frequency crossing through a zero value, this formula can also be written as
Among different special cases, we bring here two formulas:
where
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
.
The survival probabilities
rapidly diminish with the quantum number
N. For example, if
(
), we find the following surviving probabilities after the single frequency zero-crossing (when the mean energy triplicates):
This means that the initial Fock state
becomes a superposition of a large number of different Fock states
, see
Figure 9. It is impressive that the probability of transition
is very small, whereas the probability
is about 50%. In addition, the distribution of probabilities
with
looks rather irregular, whereas some regular picture is observed for
. 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
Minimizing this expression over time, one can write the minimal value
as (similar formulas were obtained, e.g., in [
26,
27,
28])
where
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,
for any (normalizable) quantum state due to the Schrödinger–Robertson uncertainty relation. The energy of quantum fluctuations
E satisfies the inequality
. Therefore, it is convenient to use two dimensionless parameters,
and
, according to the relations
and
. Then, normalizing the minimal value
by the variance in the vacuum state
, one can obtain the following formula for the
invariant squeezing coefficient:
For the states satisfying the initial conditions (
21), we have
and
. 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
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:
In particular,
for
(i.e.,
and
), so the Fock states
become squeezed for
(when
) 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
, 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.