Time-dependent quantum harmonic oscillator: a continuous route from adiabatic to sudden changes

D. Mart´ınez-Tibaduiza∗ Instituto de F´ısica, Universidade Federal Fluminense, Avenida Litorˆanea, 24210-346 Niteroi, RJ, Brazil

L. Pires Institut de Science et d’Ing´enierieSupramol´eculaires, CNRS, Universit´ede Strasbourg, UMR 7006, F-67000 Strasbourg, France

C. Farina Instituto de F´ısica, Universidade Federal do Rio de Janeiro, 21941-972 Rio de Janeiro, RJ, Brazil

In this work, we provide an answer to the question: how sudden or adiabatic is a change in the frequency of a quantum harmonic oscillator (HO)? To do this, we investigate the behavior of a HO, initially in its fundamental state, by making a frequency transition that we can control how fast it occurs. The resulting state of the system is shown to be a vacuum squeezed state in two bases related by Bogoliubov transformations. We characterize the time evolution of the squeezing parameter in both bases and discuss its relation with adiabaticity by changing the transition rate from sudden to adiabatic. Finally, we obtain an analytical approximate expression that relates squeezing to the transition rate as well as the initial and final frequencies. Our results shed some light on subtleties and common inaccuracies in the literature related to the interpretation of the adiabatic theorem for this system.

I. INTRODUCTION models to the dynamical Casimir effect [32–39] and spin states [40–42], relevant in optical clocks [43]. The main The harmonic oscillator (HO) is undoubtedly one of property of these states is to reduce the value of one of the most important systems in physics since it can be the quadrature variances (the variance of the orthogonal used to model a great variety of physical situations both quadrature is increased accordingly) in relation to coher- in classical and quantum contexts. In the former case, ent states, which already saturate the Heisenberg rela- almost all movements of physical systems with small am- tion [44–46]. This property is known as squeezing. An- plitude around a stable equilibrium configuration are de- other context where the study of the TDHO is relevant scribed by harmonic motions which explain, among many is shortcuts to adiabaticity. There, it is used in the de- other things, why the classical theory of dispersion in di- termination of optimal non-adiabatic protocols connect- electrics works so well [1, 2]. In the quantum context, ing two equilibrium states with different frequencies, in we can mention many examples, from quantum optics, such a way that the populations of the final and initial where for many purposes the quantized modes of the states are identical, a result naturally obtained through electromagnetic field behave like quantum HOs [3], to a perfectly adiabatic transformation [47–50]. It is worth quantum chemistry, where the main characteristics of the emphasizing that the study of the TDHO is also relevant London-van der Waals forces can be understood by con- in the context of quantum thermodynamics, where it can sidering two atoms as fluctuating dipoles modeled by os- be considered as a working medium in a thermodynam- cillating charges interacting through their instantaneous ical cycle [51, 52], with the frequency changes being the arXiv:2010.08081v3 [quant-ph] 9 May 2021 dipole fields [4]. quantum equivalent to volumetric expansions and con- A more general system may be a time-dependent har- tractions of the medium. In the above topics, adiabatic- monic oscillator (TDHO), i.e., a HO whose parameters: ity is a concept of fundamental importance. mass, frequency or both, are time-dependent. This sys- Adiabaticity in the TDHO, and its relation with tem has been very well studied since the Husimi’s first squeezing, has been extensively studied since its first solution in 1953 [5], and it is a natural scenario to study solution. For instance, it is well-known that any non- the important squeezed states, appearing in a variety adiabatic change in the parameters of a HO, initially in of branches of physics such as quantum optics [6–14], a coherent state, generates squeezing [5, 53]. As an im- gravitational interferometry [15–22], cosmology [23–28], portant case, a sudden change (or, simply a jump) in the metrology [29, 30], telecommunications [31], analogue frequency of a HO, has exact solution and it has been very well studied and characterized [54–63]. On the other hand, adiabatic changes do not produce squeezing and the so-called quantum adiabatic invariants of the system ∗ Correspondence to: [email protected] can be identified and used to solve the dynamics of the 2 system [5, 32, 64–68]. Nevertheless, up to the authors’ wherep ˆ andq ˆ are the usual momentum and position knowledge, in the literature there is not a full charac- operators in the Schr¨odingerpicture. Let us call Hˆ (t) the terization of squeezing for a continuous interpolation be- instantaneous hamiltonian, since it depends explicitly of tween the sudden and adiabatic limits, enabling one to time. Following Rhodes [55], in order to diagonalize Hˆ (t) answer how sudden or adiabatic is a given change in the we introduce the time-dependent operator (we are using parameters of a HO. ~ = 1) In this paper, we study numerically the time evolution r of a HO whose time-dependent frequency can be adjusted ω(t)  pˆ  ˆb(t) ≡ qˆ+ i . (2) in such a way that we can pass from an abrupt change to 2 ω(t) an adiabatic one continuously. Our starting point is the solution for a TDHO developed in Ref. [69], valid for an This operator and its hermitian adjoint satisfy the (equal- arbitrary frequency function and any initial state. This time) commutation relation: solution is extremely convenient to do numerical calcu- h i lations, and given in the initial basis, which diagonalizes ˆb(t), ˆb†(t) = 1 . (3) the hamiltonian at t = 0. We extend this solution to the instantaneous basis description, which diagonalizes the In terms of the above operators, the instantaneous hamil- hamiltonian at any time, considering the initial state as tonian, Eq.(1), reads the fundamental one. The aforementioned bases are re-  1 lated by Bogoliubov transformations, and we do explicit Hˆ (t) = ω(t) ˆb†(t)ˆb(t) + , (4) calculations for the change of basis using BCH-like rela- 2 tions established in Ref. [70]. An important advantage of our results, is that the solution in the instantaneous ba- therefore, it can be diagonalized by the eigenstates of the sis remains suitable for direct numerical implementation. operator ˆb†(t)ˆb(t). Since operators ˆb(t) and ˆb†(t) serve to The frequency transition is modeled with an emblematic diagonalize the instantaneous hamiltonian, let us refer to family of monotonic functions that asymptotically ap- them as the instantaneous representation. proach the initial and final values of the frequency (see Let us now define the usual annihilation operator, Eq. (43) and panel (a) of Fig. 1), containing a rate namely control parameter. Using this model, we initially discuss rω  pˆ  the relation between squeezing and adiabaticity, clarify- aˆ ≡ o qˆ+ i , (5) ing some subtleties and common inaccuracies present in 2 ωo the literature related to the interpretation of the adia- where ωo ≡ ω(t = 0) and which satisfies the commuta- batic theorem. Then, we use numerical analysis to fit a tion relation a,ˆ aˆ† = 1. Consequently,a ˆ† is the creation simple but very accurate analytic expression for the final operator. Notice thata ˆ anda ˆ† coincide with ˆb(t) and squeezing in the HO as a function of the parameter that ˆ† characterizes how fast the frequency changes (once the b (t) at t = 0, respectively, and they are used to diago- ˆ ˆ initial and final frequencies are given). Based on this ex- nalize the hamiltonian at t = 0, namely, H0 ≡ H(t = 0), pression we deduce the adiabatic limit and generate con- as follows tour plots, providing a quick selection of the parameters     ˆ 1 1 to obtain a desired amount of squeezing. The analysis H0 |ni = ωo nˆ + |ni = ωo n + |ni , (6) 2 2 of these results reveals an outstanding difference in the n squeezing production between enhancing or diminishing (aˆ†) withn ˆ =a ˆ†aˆ the number operator, |ni = √ |0i, and, the frequency of a HO. n! This paper is organized as follows. In Section II we consequently, {|ni} is the Fock space at t = 0. Let us present the mathematical survey for the time evolution of refer to Hˆ0 as the initial hamiltonian, {|ni} as the initial the system. In Section III we perform numerical calcula- basis and operatorsa ˆ anda ˆ† as the initial representation. tions to analyze squeezing and adiabaticity in a frequency Using Eq. (5) to write position and momentum operators transition. In Section IV we present the continuous inter- in the initial representation, and then replacing them in polation between sudden and adiabatic transitions. Sec- Eq. (2), we get tion V is left for the conclusions and final remarks. ˆb(t) = cosh(ρ(t))ˆa + sinh(ρ(t))ˆa† , (7)

II. TIME EVOLUTION where 1 ω(t) ρ(t) ≡ ln . (8) Consider a one-dimensional HO of unit mass with ar- 2 ω bitrary time-dependent frequency, ω(t). Its hamiltonian o is given by From Eqs. (3) and (7), we see that operators ˆb(t) and 1 1 ˆb†(t) are given by real Bogoliubov transformations of op- Hˆ (t) = pˆ2 + ω2(t)ˆq2 , (1) † 2 2 eratorsa ˆ anda ˆ . 3

At this time it is convenient to introduce the so-called simplicity of notation. We shall do that along the text vacuum squeezed states (for an introduction on squeezed whenever there is not risk of confusion. Notice that, from states we suggest Ref. [71]). A single-mode vacuum its definition in Eq. (8), ρ(t) is real valued. Nevertheless, squeezed state |ζi (referred to, henceforth, simply by we can identify the SP of Sˆ(ρ) with |ρ(t)|, and the Sph squeezed state) of the initial hamiltonian Hˆ0, can be ob- will be zero, for ω(t) ≥ ωo, and π, for ω(t) < ωo. Using tained by the application of the squeeze operator, Sˆ(ζ), the above results, Eq. (4) can be written as: on the fundamental state of Hˆ , namely, |ζi = Sˆ(ζ)|0i, 0 ω(t) with Hˆ (t) = Sˆ (ρ(t)) Hˆ0Sˆ (−ρ(t)) . (14) ωo ∗  ζ 2 ζ  Sˆ(ζ) ≡ exp − aˆ† + aˆ2 , (9) Accordingly, we can say that the unitary transformation 2 2 carried out by Sˆ (ρ), up to the time-dependent factor where ζ = reiϕ is a complex number and ∗ indicates com- ω(t)/ωo, has the effect to update the initial hamiltonian plex conjugation. Notice that ζ, and hence parameters to the instantaneous one. From the above result and Eq. r and ϕ, determines uniquely the squeezed state. To in- (6), it follows that terpret r and ϕ, let us introduce the initial (since it is    1   given in the initial representation) quadrature operator Hˆ (t) Sˆ(ρ(t)) |ni = ω(t) n + Sˆ(ρ(t)) |ni . (15) 2 Qˆλ, defined by [71] Hence, Sˆ(ρ(t)) |ni is an eigenstate of Hˆ (t). Let us write ˆ 1  iλ † −iλ  Qλ = √ e aˆ + e aˆ , (10) it as 2 ˆ ˆ ˆ |nit ≡ S(ρ(t)) |ni . (16) satisfying the commutation relation [Qλ, Qλ+π/2] = i. It is important to note that, rather than position and A direct consequence is that the fundamental state in momentum operators, the initial quadrature operator de- the instantaneous basis is a squeezed state in the initial pends on the initial frequency. In fact, from its definition basis, explicitly, and Eq. (5), it is related toq ˆ andp ˆ by fixing the value of √ √ the quadrature, λ, as Qˆ = (ˆa† +a ˆ)/ 2 = ω qˆ and ∞ p  n √ λ=0 o p X (2n)! 1 ˆ † √ |0it = sech(ρ) − tanh(ρ) |2ni . (17) Qλ=π/2 = i(ˆa − aˆ)/ 2 =p/ ˆ ωo . The HO is said to n! 2 be squeezed if the variance of one quadrature is smaller n=0 1 than 2 , the so-called coherent limit. It can be shown that The above result will be useful in the next section, when the variance of the initial quadrature operator Qˆλ in a discussing the behavior of the SP in a frequency transi- squeezed state of the initial hamiltonian |ζi, is given by tion. If we definen ˆ(t) = ˆb†(t)ˆb(t) as the instantaneous n [71] (ˆb†(t)) number operator, then |ni = √ |0i , and {|ni } is t n! t t 2r −2r the Fock space at time t. Henceforth, {|ni } will be refer- 2 e e t (∆Q ) = sin2 (λ − ϕ/2) + cos2 (λ − ϕ/2) . 1 λ 2 2 eed to as the instantaneous basis . Notice that the initial (11) and instantaneous bases coincide at t = 0 and whenever Notice the explicit dependence of the variance with r and ω(t) = ωo. ϕ. Even more, from the previous equation we see that At this point, it is convenient to introduce the genera- tors of the su(1,1) Lie algebra in the following form: −2r 2r e 2 e ≤ (∆Q ) ≤ . (12) 2 2 λ 2 aˆ† aˆ2 aˆ†aˆ +a ˆaˆ† Kˆ := , Kˆ := and Kˆ := , (18) + 2 − 2 c 4 From the above relation, parameter r is called squeez- ing parameter (SP), since it determines how much the satisfying the commutation relations variance of the initial quadrature operator is below the 2 h i h i coherent limit. Therefore, r and (∆Qλ) are equivalent Kˆ+, Kˆ− = −2Kˆc and Kˆc, Kˆ± = ±Kˆ± . (19) squeezing quantifiers. Parameter ϕ is referred to as the squeezing phase (Sph). Using Eqs. (7) and (18), the instantaneous hamiltonian Coming back to the squeeze operator, it is unitary, and in Eq. (4) can be written in the initial representation as transforms the annihilation operator as [55, 71] Hˆ (t) = 2ω(t) cosh(2ρ(t))Kˆc + ζ Sˆ(ζ)ˆaSˆ(−ζ) = cosh(|ζ|)ˆa + sinh(|ζ|)ˆa† . (13)  ˆ ˆ  |ζ| +ω(t) sinh(2ρ(t)) K+ + K− . (20)

Comparison of the last equation with Eq. (7), allows us to write ˆb = Sˆ(ρ)ˆaSˆ(−ρ), where we have omitted the temporal dependence in the argument of the operators for 1 This basis is also known as adiabatic basis [49, 50]. 4

Notice that in the above form of the hamiltonian, all where we defined the time-dependence is encoded in the coefficients and not in the operators, as it occurs when the instantaneous Γ1 = cosh(ρ(t)) and Γ2 = sinh(ρ(t)) . (27) representation is used (see Eq. (4)). In the following discussion we shall use the solution Also, note that using the commutation relation, Eq. (3), for a TDHO developed in Ref. [69], for an arbitrary fre- we can define the generators of the su(1,1) Lie algebra quency function and any initial state, but assuming that as: the HO is initially in its fundamental state. This solution ˆ†2 ˆ2 ˆ†ˆ ˆˆ† is given in the initial basis, and is extremely convenient b b b b + bb Tˆ+ := , Tˆ− := and Tˆc := , (28) to do numerical calculations. In this reference the au- 2 2 4 thors considered time splitting in N intervals of equally satisfying commutation relations analogous to those of small enough size τ = t/N, such that the time-dependent Eq. (19), just replacing Kˆ → Tˆ. Now, using Eqs. (26) frequency ω(t) and, consequently, the hamiltonian Eq. and (28) in Eq. (9), a generic squeeze operator in the (20), can be considered constant in each interval. Us- initial representation is transformed to the instantaneous ing algebraic methods the authors showed that the time representation as evolution operator (TEO) of the system can be written as a squeeze operator in the initial representation i.e., ˆ n ˆ ˆ ˆ o Uˆ(t, 0) = Sˆ(z(t)). Consequently, the state of the sys- S (ζ) = exp λ+T+ + λcTc + λ−T− , (29) tem at time t is a squeezed state of the initial hamilto- nian |z(t)i = Uˆ(t, 0)|0i, with z(t) = r(t)eiϕ(t). Defining where we defined the value of the frequency in the j-th time-interval as ∗ 2 2
∗ λ+ = ζ Γ2 − ζ Γ1 = −λ− , ωj := ω(jτ), with j = 1, 2, ..., N, the complex variable z(t) that characterizes the state can be calculated from ∗ the following recurrence relation [69]: λc = 2Γ1Γ2 (ζ − ζ ) . (30) ˆ ˆ χN−1bN Consequently, operators S (z(t)) and S (−ρ(t)) in Eq. χN = aN + (21) (25) are transformed to the instantaneous representation 1 − χN−1aN by substituting in Eq. (29) ζ ≡ z(t) and ζ ≡ ρ(t), respec- where tively. The presence of the term λcTˆc in Eq. (29) implies −i sinh(2ρj) sin(ωjτ) that a generic squeeze operator in the initial representa- aj = , (22) cos(ωjτ) + i cosh(2ρj) sin(ωjτ) tion is not, in general, a squeeze operator in the instan- ∗ ˆ −2 taneous one. However, since ρ = ρ , operator S (−ρ(t)) bj = (cos(ωjτ) + i cosh(2ρj) sin(ωjτ)) (23) in the instantaneous representation is given by

1  ωj  and ρj = ln . Expressing χ in its polar form, χ = n o 2 ωo Sˆ (−ρ) = exp ρTˆ+ − ρTˆ− . (31) |χ| eiθ, the state vector of the system |z(t)i = |r(t)eiϕ(t)i is obtained by the following identifications which turns out to be a squeeze operator. Now, since −1 ˆ ˆ r(t) = tanh |χN | and ϕ(t) = θN ± π , (24) S (z(t)) and S (−ρ(t)) can be identified as elements of the su(1, 1) Lie algebra, therefore we can use new BCH- with N big enough to guarantee convergence in the values like relations [70] to calculate their product as of the SP, r(t), and the Sph, ϕ(t). ˆ ˆ ˆ Recall that Eq.(24) corresponds to the solution using Sˆ (z) Sˆ (−ρ) = eαT+ eln(β)Tc eγT− , (32) the eigenstates of Hˆ0. We shall continue our discussion by describing the system using the eigenstates of Hˆ (t). where To do that, we shall write the state of the system, |z(t)i, in the instantaneous basis. Using the unitary property Γ2Λc Λc α = Λ+ + , β = 2 of the squeeze operator and Eq. (16), the state of the Γ1 − Γ2Λ− (Γ1 − Γ2Λ−) system can be written as Γ Λ − Γ and γ = 1 − 2 , (33) |z(t)i = Sˆ (z(t)) |0i Γ1 − Γ2Λ−

ˆ ˆ with = S (z(t)) S (−ρ(t)) |0it . (25) e−iϕ Γ 2 − eiϕ Γ 2
sinh(r) The solution in this approach depends on rewriting the 2 1 Λ+ = iϕ −iϕ , product of the squeeze operators, Sˆ (z(t)) Sˆ (−ρ(t)), in cosh(r) − Γ1Γ2 (e − e ) sinh(r) the instantaneous representation. To do this, first note e−iϕ Γ 2 − eiϕ Γ 2
sinh(r) that the Bogoliubov transformation given by Eq. (7) 1 2 Λ− = iϕ −iϕ , (and its hermitian adjoint) can be inverted to obtain cosh(r) − Γ1Γ2 (e − e ) sinh(r) ˆ ˆ† † ˆ† ˆ iϕ −iϕ

−2 aˆ = Γ1 b − Γ2 b anda ˆ = Γ1 b − Γ2 b , (26) Λc = cosh(r) − Γ1Γ2 e − e sinh(r) . (34) 5

Finally, using the following results Recall that the variance of the instantaneous quadrature operator could also be calculated using any other basis, ˆ 1p T− |nit = n(n − 1) |n − 2it , including the initial one. 2 We have shown that the state vector of a HO, initially 1p in its fundamental state and whose frequency begins to Tˆ+ |ni = (n + 1)(n + 2) |n + 2i , t 2 t change in time, evolves to a squeezed state of the initial 1  1 hamiltonian, when the initial basis is used to describe Tˆ |ni = n + |ni , (35) c t 2 2 t the system. Even more, when employed the instanta- neous basis to describe the system, the state vector of the Aˆ P∞ 1 ˆn and the well known expansion e = n=0 n! A (valid HO also evolves to a squeezed state, but of the instan- for a general operator Aˆ), substitution of Eq. (32) in Eq. taneous hamiltonian. In the following section we shall (25) yields study the relation between squeezing and adiabaticity in a frequency transition by varying the frequency rate. q ∞ p  n 1/2 X (2n)! 1 |z(t)i = |β| |α| eiϑ |2ni , (36) n! 2 t n=0 III. SQUEEZING AND ADIABATICITY IN A FREQUENCY TRANSITION where the overall phase was removed by the redefinitions

α = |α| eiϑ and β = |β| eiυ. (37) As we have shown in the previous section, the state vector of the HO with a time-dependent frequency can It can be checked that the following relation is satisfied be written in two different, but equivalent ways, namely: as a squeezed state relative to the initial hamiltonian, |α|2 + |β| = 1 , (38) |r(t)eiφ(t)i, or as a squeezed state relative to the instan- iΦ(t) and, accordingly, the state vector of the system in the taneous hamiltonian, |R(t)e i. Here, we shall cal- instantaneous basis can also be identified as a squeezed culate the SP’s, r(t) and R(t), for a given class of fre- state, but now of the instantaneous hamiltonian Hˆ (t). quency transitions, which are monotonic functions con- Let us call the state vector on this basis as |ξ(t)i , with necting the initial and final frequencies of the HO. The t contrast between these two bases descriptions will allow ξ(t) = R(t)eiΦ(t). Therefore, R(t) and Φ(t) are the corre- us to discuss some subtleties and common inaccuracies sponding SP and Sph of |ξ(t)i and they can be calculated t present in the literature related to the interpretation of from the complex variable α in Eq. (33) as the adiabatic theorem. We recall that the numerical im- R(t) = tanh−1 |α| and Φ(t) = ϑ ± π . (39) plementation of our previous results, obtained in section II, enables us to determine the state vector of the system In order to interpret R(t) and Φ(t), let us introduce the as accurate as desired for any frequency function, which ˆ0 instantaneous quadrature operator, Q λ, given by shows the robustness of our method. We model the frequency transition with the following 1 h i Qˆ0 (t) = √ eiλˆb†(t) + e−iλˆb(t) , (40) family of hyperbolic tangent functions [37]: λ 2   ωf + ωo (ωf − ωo) t − t0 ˆ0 ˆ0 ω(t) = + tanh , (43) satisfying the commutation relation [Qλ, Qλ+π/2] = i. 2 2  The above operator is obtained by making the changes aˆ → ˆb anda ˆ† → ˆb† in Eq. (10). Further, the instan- where ωo and ωf are the HO frequencies achieved asymp- taneous and initial quadrature operators are related by totically in the limits of remote past (t → −∞) and dis- an unitary transformation carried out by Sˆ (ρ), namely, tant future (t → +∞), respectively, t0 is the instant for ˆ0 ˆ ˆ ˆ ˆ0 which the HO frequency assumes the halfway value be- Q λ(t) = S(ρ(t))QλS(−ρ(t)). Since Q λ is expressed in the instantaneous representation, it can be shown that tween ωo and ωf , and  is, by assumption, a non-negative its variance is analogous to Eq. (11), but this time with continuous parameter, with dimensions of time, that can the SP and Sph given by Eqs. (39), namely, be adjusted appropriately for our purposes. This kind of time-dependence on the frequency makes it possible e2R e−2R to avoid possible numerical divergences [37], while the (∆Q0 )2 = sin2 (λ − Φ/2) + cos2 (λ − Φ/2) . λ 2 2 values of  control the transition rate. Indeed, from Eq. (41) (43) Moreover, from the previous equation we see that ω − ω  = f o , (44) e−2R e2R 2ω ˙ (t ) ≤ (∆Q0 )2 ≤ , (42) 0 2 λ 2 where the overdot indicates time derivative. Accordingly, which provides an equivalent relation between R(t) and the limiting case  → 0 corresponds to a sudden change 0 2 2 (∆Qλ) as that written in Eq. (12) for r(t) and (∆Qλ) . at t = t0, and the case  → ∞ corresponds to a totally 6 smooth change. Parameter  also allows us to define the transition interval, denoted by I, as the time interval (around t0) in which the frequency changes appreciably. Due to the features of the hyperbolic function and the set of parameters we shall use in the following sections, we choose I ≡ (t0 − 3, t0 + 3). Without any loss of gen- erality, in this section we set for numerical calculations ωo = 1, ωf = 3 and t0 = 10 (in arbitrary units). More- over, although the HO frequency coincides with ωo only in the remote past, we will consider that, for t < t0 − 3, ω(t) can be safely approximated by ωo. Accordingly, the state vector of the HO can be considered as being its fundamental state for t < t0 − 3. Similarly, we assume that ω(t) has already achieved its asymptotic value ωf for t > t0 + 3. In panel (a) of Fig. 1 we plot the frequency given by Eq. (43) as a function of time for different val- ues of , namely,  = 10−3, 0.5, 1.0 and 1.5. Clearly, from this graphic  = 10−3 ∼ 0.0 corresponds to the fastest transition while  = 1.5 corresponds to the slowest one. As we shall see, the previous choices for t0 and  will be quite good for our discussion.

A. SP in the initial basis

Using Eqs. (21)-(24), we calculate the SP of the HO at a generic instant t > 0 in the initial basis, namely, r(t), for the time-dependent frequency given in Eq. (43). We plot r(t) as a function of time in Fig. 1(b), employing the same values for  as those used in Fig. 1(a). From this fig- ure it can be noted that, once the transition has occurred (t > t0 +3), r(t) acquires a periodic behavior. The faster the transition is (the smaller  is), the earlier the peri- odic behavior is achieved. Also, for t > t0 + 3 both, the period and the halfway point between the minimum and maximum values of r(t) (referred to as midpoint of r(t) from now on), are the same for all transitions. With this fact in mind, and recalling that the case of a jump in the frequency ( = 0) has analytical solution, let us FIG. 1. (a) HO frequency given by Eq. (43) as a function of time; (b) time evolution of the SP in the initial basis and infer some relevant characteristics of a general transition (c) time evolution of the SP in the instantaneous basis. In (arbitrary ), from this particular limiting case. For a all panels, different line patterns mean different values of ; frequency jump, the SP can be written as [63]:  ∼ 0.0 (dot-dashed line and sudden limit),  = 0.5 (dotted line),  = 1.0 (long dashed line) and  = 1.5 (short dashed s  2 2 2 line and adiabatic limit). The solid lines ρA(t) (in panel (b)) ωf − ωo 2 r=0(t) = arcosh 1 + sin (ωf t) . (45) and ρB (t) (in panel (c)) are the plots of ρ(t) (Eq. (8)) as a 2ωoωf function of ωot, with ω(t) being the time-dependent frequency considered for  = 1.5 and  ∼ 0.0, respectively. A direct inspection of the above equation shows that the period of r=0(t) is T = π/ωf , since r=0(t + π/ωf ) = oscillations. In summary, the period and the midpoint of r=0(t). Therefore, the period of r(t) for any finite  is the r(t) oscillations depend only on the initial and final also T . Now, from Eq. (45) it can be shown that the frequency values, while the amplitude of oscillations is a (periodic) maximum of r=0(t) is 2ρf (upper horizontal signature of the transition rate. 1  ωf  line) with ρf = ln = ρ(t > t0 + 3) (see Eq. (8)), 2 ωo The behavior of the amplitude can be understood un- while the (periodic) minimum of r=0(t) is zero. Hence, der the light of the adiabatic theorem. Essentially, this once the transition has occurred, the midpoint of r(t) is theorem states that, along an adiabatic transition, the ρf (lower horizontal line) for any finite . It can also instantaneous state of the system remains the same [72]. be noted from Fig. 1(b) that, the slower the transition, Recall that the instantaneous basis is the set of eigen- the smaller the amplitude of r(t) and the smoother the states of Hˆ (t) and the dynamics begins with the sys- 7

tem in its fundamental state. Accordingly, the state of since |ξ(t)it → |0it as  → ∞. This behavior is evidenced the system in the instantaneous basis should approxi- in Fig. 1(c) and, therefore, this basis affords a direct and mate the fundamental of Hˆ (t) as the transition becomes intuitive interpretation of the adiabatic theorem. slower. Also, recall that the fundamental state of Hˆ (t) Notice that, regardless of the basis chosen to represent in the instantaneous basis is a squeezed state in the ini- the state, the behavior of the SP provides a way to inter- ˆ pret the adiabatic theorem for our system. In fact, as we tial basis, namely, |0it = S(ρ(t)) |0i = |ρ(t)i (see Eq. (17)). Hence, when the system approaches the adiabatic have shown, and as it is well known [5], the natural ba- limit, the state vector in the initial basis must satisfy sis for describing this theorem is the instantaneous one. iϕ(t) r(t)e → |ρ(t)i, or equivalently, r(t) → |ρ(t)| and However, in case one chooses another basis, all one has ϕ(t) → 2nπ, with n = 0, 1, 2, 3, .... Furthermore, since to do is to interpret the adiabatic theorem accordingly. ρ(t) can take only positive values because we choose In the following subsection we discuss some common in- ω(t) > ωo, then we should have r(t) → ρ(t). And this is accuracies in the interpretation of the adiabatic theorem exactly what occurs: as we increase the value of , r(t) for this system, and how they are related to the choice of approaches ρ(t). Moreover, as it can be seen from Fig. the basis. 1(b), for  = 1.5 (short dashed line) r(t) almost coincides with ρA(t) (solid line), the latter being ρ(t) with ω(t) given by Eq. (43) considered with  = 1.5. Although C. Discussion on squeezing and the adiabatic this description is consistent with the adiabatic theorem, theorem it may seem a little bit counterintuitive the fact that the SP has a non zero value in an adiabatic transition. How- In section II, we have shown how the squeezing pa- ever, this is a direct consequence of the basis chosen for rameters and the variance of the quadrature operators describing the state of the system and, as we shall see are related (see Eqs. (12) and (42)). Indeed, r(t) and 2 in the following subsection, when employed the instanta- (∆Qλ) are equivalent quantifiers of the squeezing in the neous basis, the SP is zero after an adiabatic transition. 0 2 initial basis, while R(t) and (∆Qλ) fulfill the same role in the instantaneous basis. Nevertheless, in the litera- ture all such quantifiers have been used to describe the B. SP in the instantaneous basis squeezing in the TDHO, and, in some cases, there has been common inaccuracies in the interpretation of the To get a direct and more intuitive interpretation of the adiabatic theorem. Here we want to address such inac- adiabatic theorem, we shall now describe the dynamics curacies, indicate their possible source, and show how of the SP in the instantaneous basis, namely, R(t). Using they can be solved. Eqs. (27), (33), (34) and (39), we calculate R(t) for the Let us initially recall that position and momentum op- frequency given by Eq. (43). Its behavior as a function erators are, up to a constant, proportional to the initial of time is presented in Fig. 1(c), with the same values for quadrature operator, as we showed in Section II. Accord-  as before. As it can be seen in this figure, in contrast to ingly, r(t) describes the squeezing associated to the vari- what happens with r(t), after the transition has occurred ance of all these operators. Also, recall that in this initial R(t) has a constant behavior, instead of a periodic one. basis perspective squeezing will inevitably be generated Notice that the maximum value reached by R(t) is ρf , for any rate of frequency change, sudden or smooth, since occurring for a jump in the frequency, which is the value r(t) 6= 0. Indeed, as we showed in Section III A, for an of r(t) for the adiabatic limit. Naively, one could think adiabatic the transition r(t) → |ρ(t)|. This non-vanishing that, for a frequency jump, the state vector of the system squeezing is a consequence of being out of the instanta- suffers a discontinuity at t = t0. However, this is not the neous basis, and represents a kind of relative squeezing. case, since, as Janszky points out in Ref. [61], even in a On the other hand, in the instantaneous basis adiabatic ˆ changes in the frequency do not produce squeezing, and frequency jump the TEO satisfy limt→t0 U→0(t, t0) = 1. It is worth emphasizing that, although the periodic be- consequently this basis seems very convenient to work havior of the squeezing seems to have been lost, it will with. Nonetheless, some authors, that have also anal- actually be present trough the squeezing phase. Indeed, ysed the TDHO in the adiabatic limit and within the it can be shown that, once the transition has occurred, initial basis perspective, incorrectly have conclude that ˆ0 adiabatic changes in the frequency do not produce this the variance of Q λ(t) for λ = 0 and λ = π/2, will oscil- late around the coherent limit (1/2) with period T . ‘relative’ squeezing (see for instance [55, 58, 73]). We be- Another interesting fact in this description is that R(t) lieve this inaccuracy has its source in the interpretation 2 0 2 approaches ρ(t) in the sudden case, and not in the adia- of r(t) as R(t), or equivalently, of (∆Qλ) as (∆Qλ) . batic one as occur for r(t). In fact, as it can be noted in In order to clarify this issue, let us indicate how 2 0 2 Fig. 1(c), for  ∼ 0 (dot-dashed line and sudden limit), (∆Qλ) and (∆Qλ) are related to each other for the rel- R(t) and ρB(t) are almost superposed, the latter being evant quadratures. First, notice that position and mo- ρ(t) with ω(t) considered with the lowest value of epsilon. mentum operators can also be written as functions of Finally, on this basis the adiabatic theorem indicates that the instantaneous quadrature operator. Explicitly, us- ˆ0 R(t) should approach zero as the transition slows down, ing Eqs. (2) and (40) it can be shown that Q λ=0(t) = 8 p ˆ0 p ω(t)ˆq and Q λ=π/2(t) =p/ ˆ ω(t). Using the simi- lar expressions that we found in Section II for position and momentum operators in terms of the initial quadra- q 0 ω(t) ture operator, we can write Qˆ (t) = Qˆλ=0 and λ=0 ωo q ˆ0 ωo ˆ Q λ=π/2(t) = ω(t) Qλ=π/2 . From the above results we obtain

0 2 ω(t) 2 (∆Qλ=0) (t) = (∆Qλ=0) (t) ωo  2 ω 2 and ∆Q0 (t) = o ∆Q
(t) . (46) λ=π/2 ω(t) λ=π/2 Thus we have found expressions that allows one to calcu- late the variance of the instantaneous quadrature opera- tor from the the variance of the initial one and vice versa. FIG. 2. R(tf ; ) as a function of  for ωo = 1 and different Similarly, expressions for R(t) as a function of r(t) and values of ωf > ωo. The pointed patterns are numerical results ϕ(t), as well as for r(t) as a function of R(t) and Φ(t), and the solid lines correspond to the fitted results using Eq. (47). can be found by substituting Eqs. (11) and (41) in the above equations. Finally, notice from Eqs. (46) that the product of the variances of the instantaneous quadrature increased). In other words, the faster the frequency tran- operator with λ = 0 and λ = π/2 is the same as for the sition occurs, the higher the squeezing. instantaneous quadrature operator with the same values Now, we calculate R(tf ; ), as a function of , for a of lambda. big enough set of different initial and final frequencies within the same order of magnitude. Using the hyper- bolic secant distribution as an Ansatz, we interpolated IV. CONTINUOUS INTERPOLATION: FROM the values of R(t ; ) as a function of the initial and final SUDDEN TO ADIABATIC f values of the frequency (within a ratio of 10), in addition 1  ωf  to . Recalling that ρf = ln = ρ(t > t0 + 3), we In this section we shall study how the SP behaves, 2 ωo once the frequency transition has occurred, for different obtain the following analytical approximate expression initial and final values of the frequency and arbitrary .A comment is in order here: up to the authors knowledge, R(tf ; ) = |ρf | sech (2 [|ρf | + 1] ωmin) , (47) there is not an analytical expression that allows one to where ωmin = min {ωo, ωf }, providing extremely good calculate exactly the correspondent SP for an arbitrary fittings for R(t ; ). Notice that the above equation is time-dependent frequency of the HO. Hence, some efforts f valid for ωf ≤ ωo and ωf ≥ ωo. For instance, for the have been made to accomplish this task. For instance, a choice made previously, namely, ω = 1 and ω = 3ω , successful approximate analytical expression that works o f o we should set ωmin = ωo in Eq. (47) to calculate the well in the sudden limit as well as in the adiabatic one corresponding fitting. Moreover, in Fig. 2 we plot with was obtained in the context of analogue models [37, 74]. solid lines R(tf ; ) as a function of  using our fitting However, this formula is no longer valid in the interme- given by Eq. (47) for the same set of frequencies used diate regime, as correctly stated by these authors. It is in our numerical results (point patterned plots). As it therefore natural to look for an analytical expression for can be seen in that figure, the fitting works very well for the SP as a function of the initial and final frequencies, ωf /ωo ≤ 5 (error smaller than a few percents) and the in addition to . agreement gets better the closer the initial and final fre- Let us denote by R(tf ; ) the value of R(t) for t > quencies are. Notice that Eq. (47) is a precise expression t0 + 3, i.e., after the transition has effectively occurred. from which we can obtain the adiabatic limit for a fre- Initially, we set the initial frequency as ωo = 1, and cal- quency transition. This can be done by taking in such culate R(tf ; ) as a function of , for different final values equation R(tf ; ) → 0 with ωo 6= ωf , resulting in of the frequency with ωf > ωo. These exact (numerical) results correspond to the pointed patterns of Fig. 2. For ω˙ (t0) instance, the (green) square points in this figure corre- << 1 , (48) ωmin(ωf − ωo)[|ρf | + 1] spond to the case ωo = 1 and ωf = 3, and therefore they show how the final height of the SP, plotted in Fig. 1(c), where it has been used Eq. (44). changes from its largest possible value R(tf ;  = 0) = ρf , Now, using Eq. (47), in Fig. 3 we present contour produced by a frequency jump, to its smallest possible plots of R(tf ; ) as a function of ωf /ωo (vertical axis) one R(tf ;  → ∞) = 0, obtained for an adiabatic tran- and ωo (horizontal axis), for the two cases: ωf /ωo > 1 sition. As expected, R(tf ; ) presents monotonically de- (panel (a)) and ωf /ωo < 1 (panel (b)). Notice that ωo, creasing behavior as the transition becomes slower ( is ωf and  are free parameters, and the unique restriction is 9

all. By the same token, we can also see in Fig. 3(a) how R(tf ; ) behaves as we vary the ratio ωf /ωo for fixed val- ues of ωo. It suffices to draw a vertical line and read the intermediate values of R(tf ; ) as we move upwards along this vertical line from the lowest value ωf /ωo = 1. As it can be noted, the maximal values achieved by R(tf ; ) (top of the vertical lines) becomes smaller as we consider higher values of ωo, as expected, since the greater ωo is, the more adiabatic the transition will be. An advantage of the contour plots shown in Fig. 3 is that they provide us with a convenient way to search for the values of the relevant parameters that lead to the desired values of the squeezing parameter. For instance, suppose we want to know the values of ωf /ωo and ωo such that R(tf ; ) < 0.1. Looking at Fig. 3 we see that the answer to this question is the darkest (blue) region limited by the line corresponding to R(tf ; ) = 0.1. For completeness, we show in Fig. 3(b) a contour plot totally analogous to that presented in Fig. 3(a), except for the fact that, now, ωf < ωo, so that ωo is the largest frequency between ωo and ωf . Drawing horizontal and vertical lines, a similar analysis to the previous one can also be made. As it can be checked, for a jump from ωo to ωf (that is, for ωo → 0), the value of R(tf ; ) will be the same if the ratios ωf /ωo in Fig. 3(a) and ωo/ωf in Fig. 3(b), are the same (see the prefactor term in Eq. (47)). However, for finite values of ωo (non-sudden change) this is not true anymore. For instance, take the ratios ωf /ωo = 5 in Fig. 3(a) and ωf /ωo = 0.2 (which means ωo/ωf = 5) and look at the value ωo = 0.4. In the former case, we see that 0.3 < R(tf ; ) < 0.4, while in the latter case 0.8 < R(tf ; ) < 0.9. Therefore, starting with ωo we will get a larger squeezing for finite values of FIG. 3. Contour plots of Eq. (47) for (a) ω /ωo > 1 and (b) f  if condition ω < ω is satisfied instead of ω > ω . ωf /ωo < 1. This graphics enables one to relate the squeezing f o f o parameter R(tf ; ) with relevant parameters of the problem Notice that ω(t) = 0 corresponds to the free particle sit- ωo, ωf and . uation, where the momentum can be infinitely squeezed, while the position is totally spread [75]. Furthermore, in a frequency decreasing (enhancing) the quadrature that that the frequencies do not differ by a factor greater than starts to be squeezed is the one proportional to the mo- 10. These contour plots allow one to explore in a more mentum (position), while the other quadrature spreads convenient way the space of parameters of the problem. [73]. Our results are consistent with the above descrip- To go deeper in the analysis of the above plots, let us tion, since R(tf ; ) increases its value monotonically as the frequency is reduced, independently of . initially consider the case ωf > ωo (Fig. 3(a)). From this plot, it can be seen, for instance, how R(tf ; ) behaves as we vary  with ωo and ωf fixed. In other words, we can recover the information contained in Fig. 2. For our V. CONCLUSIONS AND FINAL REMARKS choice made previously, namely, ωo = 1 and ωf = 3ωo, we can go to Fig. 3(a), draw a straight horizontal line at In this work we considered a harmonic oscillator (HO) ωf /ωo = 3, and read the intermediate values of R(tf ; ), undergoing a frequency transition modeled by a family given by the contour plot, as we move along this hor- of functions containing a continuous parameter, denoted izontal line from ωo = 0 (sudden limit) to the right. by , whose values, when appropriately chosen, can tune It is worth mentioning that, in addition to the choice from a sudden transition ( → 0) to an adiabatic one {ωf = 3, ωo = 1}, there are an infinity of other possibil- ( → ∞). The main purpose of our work was to relate ities with the same ratio ωf /ωo = 3. An inspection of squeezing, described by the so-called squeezing parame- Fig. 3(a) also shows that, as the ratio ωf /ωo tends to ter (SP), to adiabaticity, described by , during the time 1, the squeezing parameter approaches zero as expected, evolution of the transition. Initially, we fixed the initial since in this case there is no change in the frequency at (ωo) and final (ωf ) values of the frequency in the tran- 10 sition and studied how squeezing behaves when the fre- have shown that this result is consistent with obtaining quency rate is changed from abrupt to slow. We did that the free particle situation as the limiting case of a HO using two descriptions, namely, the instantaneous basis, with ω(t) → 0. in which the instantaneous hamiltonian is diagonal, as The model we used for the frequency transition can be well as the initial basis, in which the initial hamiltonian considered as the initial term in a sequence of transitions is diagonal. The contrast between these two bases de- to build an arbitrary ω(t). Therefore, we expect our re- scriptions allowed us to identify that the instantaneous sults to be useful in the construction of analytical expres- basis affords a direct and intuitive interpretation of the sions for the SP in HO’s with arbitrary frequency mod- adiabatic theorem, while for the initial basis the adiabatic ulations. Furthermore, the theoretical framework pre- theorem must be interpreted accordingly. Consequently, sented in this paper can be generalized to study the time we considered a big set of different values of ωo and ωf , evolution of any other initial state, including a thermal and studied how the final value of the SP, denoted by one, interesting in the realms of shortcuts to adiabatic- R(tf ; ), behaves when the frequency rate interpolates ity and quantum thermodynamics, for instance. In such continuously the sudden and adiabatic limits. Using nu- contexts, the clarification of squeezing parameters in dif- merical analysis and the hyperbolic secant distribution ferent basis might be interesting to address important as an Ansatz, we interpolated the values of R(tf ; ) as a questions about the third law of thermodynamics, quan- function of ωo and ωf , in addition to . This is how we tum speed limits, or energetic cost of shortcuts. We also obtained our main result: an analytical expression that think that our results could be useful in the study of har- allows one to calculate, directly, R(tf ; ) as a function monic traps [76–81], characterization of adiabaticity in of the relevant parameters of the problem, namely, , ωo quantum many-body systems [82] as well as in problems and ωf (within the same ratio of 10). We showed that with coupled HO’s [83, 84]. our results provides excellent fitting for R(tf ; ), allowing VI. ACKNOWLEDGMENTS us to answer the question of how sudden or adiabatic is a given change in the frequency of a quantum harmonic oscillator (HO). In order to investigate the values of the The authors acknowledge C. A. D. Zarro, A. L. C. relevant parameters so that R(tf ; ) is restricted to de- Rego, D. Szilard, Reinaldo F. de Melo e Souza, A. Z. sired intervals, we made some contour plots for R(tf ; ) as Khoury and P.A. Maia Neto for enlightening discussions. a function of ωf /ωo and . An interesting result is that, D. M. T. is also grateful to C. M. H. Willach Galliez for a sudden frequency change ( = 0), the same amount and M. Manosalva for their constant support. The au- of squeezing will be produced if we enhance or diminish thors thank the Brazilian agency for scientific and tech- the frequency of the HO by the same factor. However, nological research CAPES for partial financial support. for non-sudden transitions (finite values of  6= 0), the This work was partially supported by Conselho Nacional squeezing will be higher when the frequency is reduced de Desenvolvimento Cient´ıfico e Tecnol´ogico- CNPq, instead of enhanced by the same factor. Additionally, we 310365/2018-0 (C.F.).

[1] Andrew Zangwill. Modern electrodynamics. Cambridge [10] Roy J Glauber and M Lewenstein. Quantum optics of University Press, 2013. dielectric media. Physical Review A, 43(1):467, 1991. [2] John David Jackson. Classical electrodynamics. John [11] Jan Perina. Quantum statistics of linear and nonlinear Wiley & Sons, 2007. optical phenomena. Springer Science & Business Media, [3] Rodney Loudon. The quantum theory of light. OUP 1991. Oxford, 2000. [12] V V Dodonov. Nonclassical’states in quantum optics: [4] Peter W Milonni. The quantum vacuum: an introduction asqueezed’review of the first 75 years. Journal of Optics to quantum electrodynamics. Academic press, 2013. B: Quantum and Semiclassical Optics, 4(1):R1, 2002. [5] KˆodiHusimi. Miscellanea in elementary quantum me- [13] Victor V Dodonov and Vladimir Ivanovich Man’ko. The- chanics, ii. Progress of Theoretical Physics, 9(4):381, ory of nonclassical states of light. CRC Press, 2003. 1953. [14] Roman Schnabel. Squeezed states of light and their appli- [6] Daniel F Walls. Squeezed states of light. nature, cations in laser interferometers. Physics Reports, 684:1, 306(5939):141, 1983. 2017. [7] Rodney Loudon and Peter L Knight. Squeezed light. [15] L P Grishchuk and M V Sazhin. Squeezed quantum states Journal of modern optics, 34(6-7):709, 1987. of a harmonic oscillator in the problem of gravitational- [8] Ling-An Wu, Min Xiao, and H J Kimble. Squeezed states wave detection. Soviet Physics-JETP, 57(6):1128, 1983. of light from an optical parametric oscillator. JOSA B, [16] J Gea-Banacloche and G Leuchs. Squeezed states for 4(10):1465, 1987. interferometric gravitational-wave detectors. Journal of [9] Malvin C Teich and Bahaa E A Saleh. Squeezed state of Modern Optics, 34(6-7):793, 1987. light. Quantum Optics: Journal of the European Optical [17] Jan Harms, Yanbei Chen, Simon Chelkowski, Alexander Society Part B, 1(2):153, 1989. Franzen, Henning Vahlbruch, Karsten Danzmann, and Roman Schnabel. Squeezed-input, optical-spring, signal- 11

recycled gravitational-wave detectors. Physical Review 103(14):147003, 2009. D, 68(4):042001, 2003. [34] J R Johansson, G¨oranJohansson, C M Wilson, and [18] Hartmut Grote, K Danzmann, K L Dooley, R Schnabel, Franco Nori. Dynamical casimir effect in superconduct- J Slutsky, and H Vahlbruch. First long-term application ing microwave circuits. Physical Review A, 82(5):052509, of squeezed states of light in a gravitational-wave obser- 2010. vatory. Physical Review Letters, 110(18):181101, 2013. [35] V V Dodonov. Current status of the dynamical casimir [19] Rafa l Demkowicz-Dobrza´nski, Konrad Banaszek, and effect. Physica Scripta, 82(3):038105, 2010. Roman Schnabel. Fundamental quantum interferome- [36] Christopher M Wilson, G¨oranJohansson, Arsalan Pourk- try bound for the squeezed-light-enhanced gravitational abirian, Michael Simoen, J Robert Johansson, Tim Duty, wave detector geo 600. Physical Review A, 88(4):041802, Franco Nori, and Per Delsing. Observation of the dynam- 2013. ical casimir effect in a superconducting circuit. Nature, [20] S Dwyer, L Barsotti, SSY Chua, M Evans, M Fac- 479(7373):376, 2011. tourovich, D Gustafson, T Isogai, K Kawabe, A Kha- [37] Toshiyuki Fujii, Shigemasa Matsuo, Noriyuki Hatake- laidovski, P K Lam, et al. Squeezed quadrature fluc- naka, Susumu Kurihara, and Anton Zeilinger. Quantum tuations in a gravitational wave detector using squeezed circuit analog of the dynamical casimir effect. Physical light. Optics express, 21(16):19047, 2013. Review B, 84(17):174521, 2011. [21] Junaid Aasi, J Abadie, B P Abbott, Richard Abbott, T D [38] Pasi L¨ahteenm¨aki, GS Paraoanu, Juha Hassel, and Abbott, M R Abernathy, Carl Adams, Thomas Adams, Pertti J Hakonen. Dynamical casimir effect in a joseph- Paolo Addesso, R X Adhikari, et al. Enhanced sensitivity son metamaterial. Proceedings of the National Academy of the ligo gravitational wave detector by using squeezed of Sciences, 110(11):4234, 2013. states of light. Nature Photonics, 7(8):613, 2013. [39] S Felicetti, M Sanz, L Lamata, G Romero, G¨oranJo- [22] Benjamin P Abbott, R Abbott, T D Abbott, M R Aber- hansson, Per Delsing, and E Solano. Dynamical casimir nathy, F Acernese, K Ackley, C Adams, T Adams, P Ad- effect entangles artificial atoms. Physical Review Letters, desso, R X Adhikari, et al. Gw150914: First results 113(9):093602, 2014. from the search for binary black hole coalescence with [40] Masahiro Kitagawa and Masahito Ueda. Squeezed spin advanced ligo. Physical Review D, 93(12):122003, 2016. states. Physical Review A, 47(6):5138, 1993. [23] L P Grishchuk and Yu V Sidorov. Squeezed quantum [41] Ian D Leroux, Monika H Schleier-Smith, and Vladan states of relic gravitons and primordial density fluctua- Vuleti´c. Implementation of cavity squeezing of a collec- tions. Physical Review D, 42(10):3413, 1990. tive atomic spin. Physical Review Letters, 104(7):073602, [24] L P Grishchuk. Quantum effects in cosmology. Classical 2010. and Quantum Gravity, 10(12):2449, 1993. [42] Onur Hosten, Nils J Engelsen, Rajiv Krishnakumar, [25] Andreas Albrecht, Pedro Ferreira, Michael Joyce, and and Mark A Kasevich. Measurement noise 100 times Tomislav Prokopec. Inflation and squeezed quantum lower than the quantum-projection limit using entangled states. Physical Review D, 50(8):4807, 1994. atoms. Nature, 529(7587):505, 2016. [26] Bei Lok Hu, G W Kang, and Andrew Matacz. Squeezed [43] Chin-wen Chou, D B Hume, J C J Koelemeij, David J vacua and the quantum statistics of cosmological particle Wineland, and T Rosenband. Frequency comparison of creation. International Journal of Modern Physics A, two high-accuracy al+ optical clocks. Physical Review 9(07):991, 1994. Letters, 104(7):070802, 2010. [27] Martin B Einhorn and Finn Larsen. Squeezed states in [44] John R Klauder and Bo-Sture Skagerstam. Coherent the de sitter vacuum. Physical Review D, 68(6):064002, states: applications in physics and mathematical physics. 2003. World scientific, 1985. [28] Claus Kiefer, Ingo Lohmar, David Polarski, and Alexei A [45] K Wodkiewicz and J H Eberly. Coherent states, squeezed Starobinsky. Pointer states for primordial fluctuations in fluctuations, and the su (2) am su (1, 1) groups in inflationary cosmology. Classical and Quantum Gravity, quantum-optics applications. JOSA B, 2(3):458, 1985. 24(7):1699, 2007. [46] Jean-Pierre Gazeau. Coherent states in quantum physics. [29] Henning Vahlbruch, Simon Chelkowski, Karsten Danz- Wiley, 2009. mann, and Roman Schnabel. Quantum engineering of [47] Xi Chen, A Ruschhaupt, Sebastian Schmidt, Adolfo squeezed states for quantum communication and metrol- del Campo, David Gu´ery-Odelin,and J Gonzalo Muga. ogy. New Journal of Physics, 9(10):371, 2007. Fast optimal frictionless atom cooling in harmonic traps: [30] Vittorio Giovannetti, Seth Lloyd, and Lorenzo Mac- Shortcut to adiabaticity. Physical Review Letters, cone. Advances in quantum metrology. Nature photonics, 104(6):063002, 2010. 5(4):222, 2011. [48] Jean-Fran¸coisSchaff, Xiao-Li Song, Patrizia Vignolo, and [31] Radan Slav´ık, Francesca Parmigiani, Joseph Kakande, Guillaume Labeyrie. Fast optimal transition between Carl Lundstr¨om, Martin Sj¨odin, Peter A Andrekson, two equilibrium states. Physical Review A, 82(3):033430, Ruwan Weerasuriya, Stylianos Sygletos, Andrew D Ellis, 2010. Lars Gr¨uner-Nielsen, et al. All-optical phase and ampli- [49] A Del Campo. Frictionless quantum quenches in ultra- tude regenerator for next-generation telecommunications cold gases: A quantum-dynamical microscope. Physical systems. Nature Photonics, 4(10):690, 2010. Review A, 84(3):031606, 2011. [32] V V Dodonov and A V Dodonov. Quantum harmonic [50] David Gu´ery-Odelin, Andreas Ruschhaupt, Anthony oscillator and nonstationary casimir effect. Journal of Kiely, Erik Torrontegui, Sofia Mart´ınez-Garaot, and Russian Laser Research, 26(6):445, 2005. Juan Gonzalo Muga. Shortcuts to adiabaticity: Con- [33] J Robert Johansson, G¨oranJohansson, CM Wilson, and cepts, methods, and applications. Reviews of Modern Franco Nori. Dynamical casimir effect in a supercon- Physics, 91(4):045001, 2019. ducting coplanar waveguide. Physical Review Letters, 12

[51] E. Torrontegui and R. Kosloff. Quest for absolute zero in [70] D. Mart´ınez-Tibaduiza, A.H. Arag˜ao, C. Farina, and the presence of external noise. Physical Review E, 88(3), C.A.D. Zarro. New BCH-like relations of the su(1,1), 2013. su(2) and so(2,1) Lie algebras. Physics Letters A, [52] Ronnie Kosloff and Yair Rezek. The quantum harmonic 384(36):126937, 2020. otto cycle. Entropy, 19(4):136, 2017. [71] Stephen Barnett and Paul M Radmore. Methods in the- [53] R Graham. Squeezing and frequency changes in harmonic oretical quantum optics, volume 15. Oxford University oscillations. Journal of Modern Optics, 34(6-7):873, 1987. Press, 2002. [54] J Janszky and Y Y Yushin. Squeezing via frequency [72] M. Born and V. Fock. Beweis des adiabatensatzes. jump. Optics communications, 59(2):151, 1986. Zeitschrift f¨urPhysik, 51(3-4):165–180, mar 1928. [55] Xin Ma and William Rhodes. Squeezing in harmonic os- [73] I Averbukh, B Sherman, and G Kurizki. Enhanced cillators with time-dependent frequencies. Physical Re- squeezing by periodic frequency modulation under para- view A, 39(4):1941, 1989. metric instability conditions. Physical Review A, [56] CF Lo. How does a squeezed state of a general driven 50(6):5301, 1994. time-dependent oscillator evolve? Physica Scripta, [74] Shigemasa Matsuo, Toshiyuki Fujii, and Noriyuki Hatak- 42(4):389, 1990. enaka. Nonadiabatic squeezed-photon generation by a [57] C F Lo. Generating displaced and squeezed number fourier-modified janszky–adam scheme. Physica B: Con- states by a general driven time-dependent oscillator. densed Matter, 468:57, 2015. Physical Review A, 43(1):404, 1991. [75] G. W. Ford and R. F. O’Connell. Wave packet spreading: [58] G S Agarwal and S Arun Kumar. Exact quantum- Temperature and squeezing effects with applications to statistical dynamics of an oscillator with time-dependent quantum measurement and decoherence. American Jour- frequency and generation of nonclassical states. Physical nal of Physics, 70(3):319–324, 2002. Review Letters, 67(26):3665, 1991. [76] Siegfried Grossmann and Martin Holthaus. On bose- [59] J Janszky and P Adam. Strong squeezing by repeated einstein condensation in harmonic traps. Physics Letters frequency jumps. Physical Review A, 46(9):6091, 1992. A, 208(3):188, 1995. [60] T Kiss, P Adam, and J Janszky. Time-evolution of a har- [77] Raam Uzdin, Emanuele G. Dalla Torre, Ronnie Kosloff, monic oscillator: jumps between two frequencies. Physics and Nimrod Moiseyev. Effects of an exceptional point Letters A, 192(5-6):311, 1994. on the dynamics of a single particle in a time-dependent [61] T Kiss, J Janszky, and P Adam. Time evolution of harmonic trap. Physical Review A, 88(2), 2013. harmonic oscillators with time-dependent parameters: [78] K. H. Hoffmann, B. Andresen, and P. Salamon. Optimal A step-function approximation. Physical Review A, control of a collection of parametric oscillators. Physical 49(6):4935, 1994. Review E, 87(6), 2013. [62] H´ectorMoya-Cessa and Manuel Fern´andezGuasti. Co- [79] V.G. Ibarra-Sierra, J.C. Sandoval-Santana, J.L. Cardoso, herent states for the time dependent harmonic oscillator: and A. Kunold. Lie algebraic approach to the time- the step function. Physics Letters A, 311(1):1, 2003. dependent quantum general harmonic oscillator and the [63] D M Tibaduiza, L Pires, D Szilard, C A D Zarro, C Fa- bi-dimensional charged particle in time-dependent elec- rina, and A L C Rego. A time-dependent harmonic os- tromagnetic fields. Annals of Physics, 362:83–117, 2015. cillator with two frequency jumps: an exact algebraic [80] Fabienne Schneiter, Sofia Qvarfort, Alessio Serafini, solution. Brazilian Journal of Physics, 50(5):634, 2020. Andr´e Xuereb, Daniel Braun, Dennis R¨atzel, and [64] H Ralph Lewis Jr and W B Riesenfeld. An exact quantum David Edward Bruschi. Optimal estimation with quan- theory of the time-dependent harmonic oscillator and of tum optomechanical systems in the nonlinear regime. a charged particle in a time-dependent electromagnetic Physical Review A, 101(3):033834, 2020. field. Journal of mathematical physics, 10(8):1458, 1969. [81] Sofia Qvarfort, Alessio Serafini, Andr´eXuereb, Daniel [65] V S Popov and A M Perelomov. Parametric excitation Braun, Dennis R¨atzel,and David Edward Bruschi. Time- of a quantum oscillator. Soviet Physics JETP, 29(4):738, evolution of nonlinear optomechanical systems: Interplay 1969. of mechanical squeezing and non-gaussianity. Journal of [66] I A Malkin and V I Man’Ko. Coherent states and ex- Physics A: Mathematical and Theoretical, 53(7):075304, citation of n-dimensional nonstationary forced oscillator. 2020. Technical report, CM-P00062575, 1970. [82] Amy H. Skelt and Irene D'Amico. Characterizing adia- [67] I A Malkin, V I Man’Ko, and D A Trifonov. Coherent baticity in quantum many-body systems at finite temper- states and transition probabilities in a time-dependent ature. Advanced Quantum Technologies, 3(7):1900139, electromagnetic field. Physical Review D, 2(8):1371, 2020. 1970. [83] Alejandro R Urz´ua, Ir´an Ramos-Prieto, Manuel [68] Frank Boldt, James D. Nulton, Bjarne Andresen, Peter Fern´andez-Guasti, and H´ector M Moya-Cessa. Solu- Salamon, and Karl Heinz Hoffmann. Casimir companion: tion to the time-dependent coupled harmonic oscillators An invariant of motion for hamiltonian systems. Physical hamiltonian with arbitrary interactions. Quantum Re- Review A, 87(2), 2013. ports, 1(1):82–90, 2019. [69] D M Tibaduiza, L Pires, A L C Rego, D Szilard, C Zarro, [84] Alejandro R Urz´ua,Ir´anRamos-Prieto, Francisco Soto- and C Farina. Efficient algebraic solution for a time- Eguibar, V´ıctor Arriz´on, and H´ector M Moya-Cessa. dependent quantum harmonic oscillator. Physica Scripta, Light propagation in inhomogeneous media, coupled 95(10):105102, 2020. quantum harmonic oscillators and phase transitions. Sci- entific Reports, 9(1):1, 2019.