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

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.


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 timedependent frequency ω(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(at)| 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. Sections 8,9 and 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.

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 ω 2 (t) = ω 2 0 f (t/τ), where f (t/τ) is a non-negative function with the properties We assume that ω 2 (t) = ω 2 0 = const 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ẋ 0 in the form x 0 = cos(ϕ) 2E 0 /ω 0 ,ẋ 0 = sin(ϕ) 2E 0 , 0 ≤ ϕ < 2π.
To solve Equation (2) numerically, we introduce dimensionless variables X = ω 0 x/ √ 2E 0 and T = t/τ, arriving at the equation for T ≥ −1 with the initial conditions Then, the dimensionless energy ratio R = E(t)/E 0 depends on the dimensionless time T and two parameters, G and ϕ: where the derivative dX/dT 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).

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 The adiabatic ratio must equal R n (T; G) = |T| n/2 . In Figures 1 and 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.

A Tanh-like Profile of the Frequency
The second example is the profile f (T) = | tanh(aT)/ 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(aT)/ 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.

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 The Wronskian identity for the solutions ε(t) and ε * (t) has the forṁ Then, we can write at 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 ≥ −τ: 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 x and p are totally independent from the dynamics of the variances σ x = x 2 − x 2 , σ p = p 2 − p 2 and σ xp = xp + px /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 , σ xp ).
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 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, where Constant complex coefficients u ± 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 ω(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 Figures 3 and 4. For this reason, we concentrate hereafter on the important special case when It includes the vacuum, thermal and Fock initial quantum states. Then, ∆β = 0. In addition, many formulas can be simplified: 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")ρ which follows from (2) if one writes ε = ρ exp(iχ) 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ρ(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.

Exact Solutions for the Power Profile of the Frequency
It is known that Equation (2) with the time-dependent frequency ω 2 (t) = ω 2 0 |t/τ| n can be reduced to the Bessel equation 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: 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: Constant complex coefficients A − and B − can be found from the initial conditions (7). Remembering that d|t|/dt = −1 for t < 0, one obtains the following equations: where J ±ν (z) means the derivative of the Bessel function J ±ν (z) 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 A − and B − : The time derivative of function (28) at t < 0 (when d|t|/dt = −1) can be written with the aid of identities (29) as follows: On the other hand, 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, yJ 1−ν (y)/ |t| → 0 at t → 0, while the product yJ ν−1 (y)/ |t| tends to a finite value in this limit. Hence, the continuity of derivative dε/dt 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 Using Equations (28) and (30), one can write down formula (23) for the mean energy ratio R(t) = E (t)/E (−τ) as follows: Figure 6 shows the function R(b), where −1 ≤ b ≡ t/τ ≤ 1, for ν = 1/4 (i.e., n = 2 and y = gb 2 ) and three values of parameter g = 0.1, 1.0, 10. Using the known formulas 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], results in the following simple expressions for z 1: πz .
Hence, in the adiabatic limit (g 1 and y 1), we obtain 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 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(τ), 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.

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: This formula can be derived directly from Equation (28). If g 1, Equation (30) assumes the following asymptotic form: Then, Equation (28) results in the following expressions for y 1: 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 ig in Equation (43) and comparing this equation with (13) and (15), we obtain the following expressions for the coefficients u ± : 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 ± |, seem to be universal after a single frequency passage through zero.

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 ξ = tanh(κt/2). Then, Equation (2) takes the form 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.

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, We suppose that the transition rule through the second zero has the form 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 * * ) Hence, 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 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, In particular, if w ± = u ± , then β min = 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, β max = 1 + 8|u + u − | 2 under the same conditions. If n = 2, β max = 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) = ω 2 0 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, Φ = 4G/π. 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.

Energy Fluctuations
Figures 3, 4 and 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, assuminḡ h = m = ω 0 = 1, so that x 4 = p 4 ): where A(t) = ω 2 (t)|ε(t)| 2 + |ε(t)| 2 , B(t) = Re ω 2 (t)ε 2 (t) +ε 2 (t) , C(t) = Im ω 2 (t)ε 2 (t) +ε 2 (t) .
In the adiabatic regime (15), we have For the initial Fock state |N , we have Hence, 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: 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 2 t = 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 ≡Ê/(hω) − 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: Consequently, In particular, for |u − | 2 = 1, 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 Q ≈ (4/3)N ≈ (4/9) n for N 1, whereas Q = n for the "strongly super-Poissonian" thermal states.

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â †â , 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 |N t remains the eigenstate of operatorÂ †Â . As soon as operatorsÂ andÂ † maintain their linear form with respect to operatorsx andp, 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 withh = m = 1), 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) where P k j (z) is the associated Legendre polynomial, N < = min(M, N), N > = max(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, 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 P N (z) ≡ P 0 N (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): 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 → 3N is very small, whereas the probability p N (M ≥ 3N) is about 50%. In addition, the distribution of probabilities p N (M) with M < 3N looks rather irregular, whereas some regular picture is observed for M > 3N. Unfortunately, we did not succeed to find an analytic approximation for this regular picture.

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 σ min as (similar formulas were obtained, e.g., in [26][27][28]) 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 ≥h 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 ≥hω/2. Therefore, it is convenient to use two dimensionless parameters, λ ≥ 1 and γ ≥ 1, according to the relations E = λhω/2 and D = γ 2h2 /4. Then, normalizing the minimal value σ min by the variance in the vacuum state σ vac =h/(2mω), one can obtain the following formula for the invariant squeezing coefficient s = σ min /σ vac : 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 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).

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.