Arxiv:2010.08081V3 [Quant-Ph] 9 May 2021 Dipole fields [4]
Total Page:16
File Type:pdf, Size:1020Kb
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); ^by(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) ^by(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 ^by(t)^b(t). Since operators ^b(t) and ^by(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^y = 1. Consequently,a ^y is the creation simple but very accurate analytic expression for the final operator. Notice thata ^ anda ^y coincide with ^b(t) and squeezing in the HO as a function of the parameter that ^y 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 jni = !o n^ + jni = !o n + jni ; (6) 2 2 of these results reveals an outstanding difference in the n squeezing production between enhancing or diminishing (a^y) withn ^ =a ^ya^ the number operator, jni = p j0i, and, the frequency of a HO. n! This paper is organized as follows. In Section II we consequently, fjnig 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, fjnig as the initial the system. In Section III we perform numerical calcula- basis and operatorsa ^ anda ^y 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))^ay ; (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 ^by(t) are given by real Bogoliubov transformations of op- H^ (t) = p^2 + !2(t)^q2 ; (1) y 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.