REVIEWS OF MODERN PHYSICS, VOLUME 91, OCTOBER–DECEMBER 2019 Shortcuts to adiabaticity: Concepts, methods, and applications
D. Gu´ery-Odelin Laboratoire de Collisions Agr´egats R´eactivit´e, CNRS UMR 5589, IRSAMC, Universit´e de Toulouse (UPS), 118 Route de Narbonne, 31062 Toulouse CEDEX 4, France
A. Ruschhaupt and A. Kiely Department of Physics, University College Cork, Cork, Ireland
E. Torrontegui Instituto de Física Fundamental IFF-CSIC, Calle Serrano 113b, 28006 Madrid, Spain
S. Martínez-Garaot and J. G. Muga Departamento de Química Física, UPV/EHU, Apdo. 644, 48080 Bilbao, Spain
(published 24 October 2019) Shortcuts to adiabaticity (STA) are fast routes to the final results of slow, adiabatic changes of the controlling parameters of a system. The shortcuts are designed by a set of analytical and numerical methods suitable for different systems and conditions. A motivation to apply STA methods to quantum systems is to manipulate them on timescales shorter than decoherence times. Thus shortcuts to adiabaticity have become instrumental in preparing and driving internal and motional states in atomic, molecular, and solid-state physics. Applications range from information transfer and processing based on gates or analog paradigms to interferometry and metrology. The multiplicity of STA paths for the controlling parameters may be used to enhance robustness versus noise and perturbations or to optimize relevant variables. Since adiabaticity is a widespread phenomenon, STA methods also extended beyond the quantum world to optical devices, classical mechanical systems, and statistical physics. Shortcuts to adiabaticity combine well with other concepts and techniques, in particular, with optimal control theory, and pose fundamental scientific and engineering questions such as finding speed limits, quantifying the third law, or determining process energy costs and efficiencies. Concepts, methods, and applications of shortcuts to adiabaticity are reviewed and promising prospects are outlined, as well as open questions and challenges ahead.
DOI: 10.1103/RevModPhys.91.045001
CONTENTS 2. Examples of invariant-based inverse engineering 13 a. Two-level system 13 I. Introduction 2 b. Lewis-Leach family 14 A. Overview of shortcuts to adiabaticity 2 3. Scaling laws 15 B. Motivation and scope of this review 3 4. Connection with Lax pairs 16 II. Methods 4 D. Variational methods 16 A. Overview of inverse engineering approaches 4 E. Fast forward 17 1. Quantum transport 5 1. The original formalism 17 2. Spin manipulation 5 2. Streamlined fast-forward approach 18 3. Beyond mean values 5 3. Generalizations and terminology 19 B. Counterdiabatic driving 6 F. FAQUAD and related approaches 19 1. Superadiabatic iterations 8 G. Optimal control and shortcuts to adiabaticity 20 2. Beyond the basic formalism 9 H. Robustness 21 “ ” a. Physical unitary transformations 9 1. Error sensitivity and its optimization using b. Schemes that focus on one state 10 perturbation theory 21 c. Effective counterdiabatic field 10 a. Illustrative example: Control of a two-level d. Dressed-states approach 10 system 21 e. Variational approach 11 b. Optimization using perturbation theory f. Counterdiabatic Born-Oppenheimer dynamics 12 in other settings 22 g. Constant CD-term approximation 12 2. Other approaches 23 C. Invariants and scaling laws 12 I. Three-level systems 23 1. Lewis-Riesenfeld invariants 12
0034-6861=2019=91(4)=045001(54) 045001-1 © 2019 American Physical Society D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
1. Apply counterdiabatic shortcuts to the full three-level Λ system 23 2. Applying invariant-based inverse engineering shortcuts to the full three-level Λ system 24 3. Using STA techniques after mapping to a two-level system 24 J. Motional states mapped into a discrete system 24 III. Applications in Quantum Science and Technology 24 A. Trapped ions 25 1. Dynamical normal modes 25 2. Ion transport 27 3. Other operations 28 B. Double wells 29 FIG. 1. A turtle on wheels is a good metaphor for shortcuts to C. Cavity quantum electrodynamics 29 adiabaticity. The image is inspired by the artist work of Andree D. Superconducting circuits 30 Richmond. E. Spin-orbit coupling 31 F. Nitrogen-vacancy centers 31 G. Many-body and spin-chain models 32 illustrates the concept pictorially by comparing diabatic, H. Metrology 34 adiabatic, and STA paths in a discrete system. IV. The Energy Cost of STA, Engines, and the Third Law 34 The shortcuts rely on specific time dependences of the A. Energy costs 34 control parameters and/or on the addition of auxiliary time- B. Engines and refrigerators 37 dependent couplings or interactions with respect to some C. Third law 38 reference Hamiltonian or, more generally, a Liouvillian or V. Open Quantum Systems 38 transition-rate matrix. STA methods were first applied in A. Concept of adiabaticity for open systems 38 simple quantum systems: two- and three-level systems or a B. Engineering the environment 39 particle in a time-dependent harmonic oscillator. They have C. Mitigating the effect of environment 39 since come to encompass a much broader domain since slow D. Non-Hermitian Hamiltonians 40 processes are quite common as a simple way to prepare the VI. Optical Devices 40 state of a system or to change conditions avoiding excitations VII. Extension to Classical and Statistical Physics 42 A. Counterdiabatic methods in classical mechanics 42 in a wide spectrum of areas from atomic, molecular, and B. Mechanical engineering 43 optical physics, solid state, or chemistry to classical mechani- 1. Cranes 43 cal systems and engineering. In parallel to such a large scope, 2. Robustness issues 44 different methods have been developed and applied. C. STA for isolated dilute gases 44 This description of shortcuts needs some caveats and 1. Boltzmann equation 44 clarifications. In many quantum mechanical applications, 2. Extension to Navier-Stokes equation 45 the final state in an STA process reproduces the set of D. Shortcuts for systems in contact with a thermostat 45 adiabatic probabilities to find the system in the eigenstates 1. The overdamped regime 45 of the final Hamiltonian. The mapping of probabilities from 2. Connection with free energy and irreversible work 46 3. Extensions 46 VIII. Outlook and Open Questions 46 List of Symbols and Abbreviations 47 Acknowledgments 47 Appendix A: Example of Lie Transform 47 Appendix B: Counterdiabatic Hamiltonian for a Two-level Hamiltonian with Complex Coupling 48 References 48
I. INTRODUCTION
A. Overview of shortcuts to adiabaticity “ ” Shortcuts to adiabaticity (STA) are fast routes to the final FIG. 2. Schematic example of adiabatic, diabatic, and STA results of slow, adiabatic changes of the controlling param- processes. The system is initially in the third level. In the eters of a system. Adiabatic processes are here broadly defined adiabatic path (short-dashed red) the system evolves always as those for which the slow changes of the controls leave some along the third level. In a diabatic evolution (long-dashed blue) dynamical properties invariant, the “adiabatic invariants,” the system gets excited by jumping across avoided crossings. such as the quantum number in quantum systems or phase- Along an STA path (green dots) the system does not always travel space areas in classical systems. Figure 1 informally portrays along the third level but it arrives at the third level in a shorter the central idea of shortcuts to adiabaticity. Figure 2 also time (the time dimension is not explicitly shown in the figure).
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-2 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … initial to final settings, without final excitation but allowing the design of operations explicitly intended to achieve it (see for excitations en route, holds for any possible initial state in an example in Sec. III.A), as well as indirectly, since a some state-independent STA protocols and simple enough mitigated decoherence reduces the need for highly demand- systems. These state-independent shortcuts are ideal for ing, qubit-consuming error-correction codes. maximal robustness and to reduce dependence on temper- • The second main paradigm of quantum information ature, so as to avoid costly and time-consuming cooling to the processing is adiabatic computing or quantum annealing, ground state. However, except for idealized models, state which may be accelerated or made possible by changing independence holds only approximately and for a domain of initial and/or final Hamiltonians, modifying the interpolation parameters, for example, within the harmonic approximation path between initial and final Hamiltonians, or adding in a trapped ion. In other applications state independence is auxiliary terms to the transient Hamiltonian. For adiabatic not really necessary or it may be difficult or impossible to quantum computing and other applications with complex achieve as in chaotic systems. Yet shortcuts may still be found systems, progress is being made to find effective STA without for a chosen specific state, typically the ground state or for a using difficult-to-find information on spectra and eigenstates. state subspace. The more recent paradigms of topological quantum informa- A further caveat is that the scope of STA methods has in tion processing or measured based quantum computation may practice outstripped the original, pure aim in many applica- also benefit from STA. tions. For example, these methods can be extended with Shortcuts to adiabaticity are also having an impact or are minimal changes to drive general transitions, regardless of expected to be useful in quantum technologies other than whether the initial and final states can be connected adiabati- quantum information processing such as interferometry and cally, such as transitions where the initial state is an eigenstate metrology, communications, and microengines or refrigera- of the initial Hamiltonian, whereas the final state is not an tors in quantum thermodynamics. eigenstate of the final Hamiltonian. This broader perspective Fundamental questions.—Finally, apart from being a prac- merges with inverse engineering methods of the Hamiltonian tical aid to process design, STA also played a role and will to achieve arbitrary transitions or unitary transformations. continue to be instrumental in clarifying fundamental con- Motivations.—There are different motivations for the cepts such as quantum-classical relations, quantum speed speedup that depend on the setting. In optics, time is often limits and trade-off relations between timing, energy, robust- substituted by length to quantify the rate of change so the ness, entropy or information, the third law of thermodynam- shortcuts imply shorter, more compact optical devices. In ics, and the proper characterization of energy costs of mechanical engineering, we look for fast and safe protocols, processes. say of robotic cranes, to enhance productivity. In microscopic quantum systems, slowness often implies decoherence, the B. Motivation and scope of this review accumulation of errors and perturbations, or even the escape of the system from its confinement. The shortcuts provide a This review provides a broad overview of concepts, methods, useful toolbox to avoid or mitigate these problems and thus to and applications of shortcuts to adiabaticity. A strong motivation develop quantum technologies. Moreover, with shorter proc- for writing it is the breadth of systems and areas where STA are ess times experiments can be repeated more often to increase used. The hope is that this work will assist in stimulating signal-to-noise ratios. discussions and information transfer between different domains. A generic and important feature of STA apart from the There are already examples that demonstrate the benefits of such speed achieved is that there are typically many alternative an interaction. Analogies found between systems as disparate as routes for the control parameters, and this flexibility can be individual ions in Paul traps and mechanical cranes or optical used to optimize physically relevant variables, for example, waveguide devices and atomic internal states have proved to minimize transient energy excitations and/or energy con- fruitful. The link between Lewis-Riesenfeld invariants and sumptions or to maximize robustness against perturbations. Lax pairs is another unexpected example. Shortcuts to adiabaticity in quantum technologies.— The development of STA methods and applications has Indeed, shortcuts have been mostly developed and applied been quite rapid since Chen, Ruschhaupt et al. (2010), where for quantum systems, and much of the current interest in STA the expression shortcuts to adiabaticity was coined. is rooted in the quest for quantum technologies. STA combine Precedents before 2010 exist where STA concepts had been well with the two main paradigms of quantum information applied (Unanyan et al., 1997; Emmanouilidou et al., 2000; processing. Couvert et al., 2008; Motzoi et al., 2009; Muga et al., 2009; • In the gate-based paradigm, shortcuts to adiabaticity Rezek et al., 2009; Salamon et al., 2009; Schmiedl et al., contribute to improve and speed up gates or state preparations 2009; Masuda and Nakamura, 2010). In the post-2010 period, and to perform elementary operations like moving atoms the early work of Demirplak and Rice (2003, 2005, 2008), and leaving them unexcited. Shortcuts to adiabaticity may be Berry (2009) on “counterdiabatic (CD) driving” has been applied to all physical platforms, such as trapped ions, cold particularly influential, as well as methodologies such as atoms, nitrogen-vacancy (NV) centers, superconducting cir- “inverse engineering,”“invariants,”“scaling laws,”“fast- cuits, quantum dots, or atoms in cavities. forward,” (FF) or “local adiabatic” methods. To add further Apart from fighting decoherence, implementing “scalable” flexibility to this already rich scenario of approaches, some architectures toward larger quantum systems is a second methods provide in general a multiplicity of control protocols. major challenge to develop quantum information processing. Moreover, STA methods relate synergistically to or overlap Scalability benefits from STA in two ways: directly, through partially with other control methods. In particular STA blend
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-3 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … well with optimal control theory (OCT), decoherence-free experiments, and theoretical results since its publication subspaces (DFSs) (S. L. Wu et al., 2017), linear response clearly surpasses the work reviewed there. We mostly pay theory (Acconcia, Bonança, and Deffner, 2015), or perturba- attention to the work done after Torrontegui, Ibáñez et al. tive and variational schemes (see Secs. II.H and II.D). There is (2013), but for presenting the key concepts some overlap is in summary a dense network of different approaches and STA allowed. Work on the counterdiabatic method was reviewed protocols available, frequently hybridized. recently by Kolodrubetz et al. (2017) with emphasis on Some of the basic STA techniques rely on specific geometric aspects and classical-quantum relations. For the formalisms (invariants and scaling, CD driving, or fast- approaches to speed up adiabatic computing, see Albash and forward) that, within specific domains, can be related to each Lidar (2018) and Takahashi (2019). For a review on STA to other and be made potentially equivalent because of under- control quantum critical dynamics, see del Campo and lying common structures. For example, a particular protocol to Sengupta (2015). See also Menchon-Enrich et al. (2016) speed up the transport of a particle by defining a trap path may about the spatial adiabatic passage and Vitanov et al. (2017) be found as the result of “invariant-based engineering,” fast- about the stimulated Raman adiabatic passage (STIRAP), a forward, or “local-CD” approaches. This convergence for technique that has often been sped up with STA. Finally, specific operations and methods might be misleading because Masuda and Rice (2016) reviewed the counterdiabatic and it does not extend to all systems and circumstances. No single fast-forward methods focusing on applications to polyatomic all-inclusive theory exists for all STA, and each of the major molecules and Bose-Einstein condensates (BEC). Also, a methodologies carries different limitations, different construc- Focus issue of New Journal of Physics on Shortcuts to tion recipes, and different natural application domains. Adiabaticity provides a landscape of current tendencies (del “ ” A widespread misconception is to identify STA with one Campo and Kim, 2019). particular approach, often CD driving, and even with one Structure.—The review is organized into “Methods” of particular protocol within one approach. We shall indeed pay STA in Sec. II and “Applications” from Sec. III to VII.We due attention to unifying concepts and connections among address applications of STA to quantum science and tech- the STA formalisms, which are of course useful and worth nology in Sec. III. This includes the main physical platforms stressing, but at the same time it is hard to overemphasize that and the quantum system is mainly considered as a closed the diversity of approaches is a powerful asset of the shortcuts quantum system. The connections between STA and quantum to adiabaticity that explains their versatility. thermodynamic concepts, a topic between closed and open New hybrid or approximate methods are being created as quantum systems, are reviewed in Sec. IV. STA and open we write and more will be devised to adapt to diverse needs quantum systems is the topic of Sec. V. Sections VI and VII and systems. This review is also intended to map and deal with results and peculiarities of two fields that also offer a characterize the options, help users to navigate among them, promising arena for practical applications: optics in Sec. VI and encourage the invention of novel, more efficient or goal- and classical systems in Sec. VII. adapted approaches. A number of techniques that fit naturally We tried to keep acronyms to a minimum but some terms into the definition of STA but have not been usually tagged as [among the most prominent: shortcuts to adiabaticity (STA) STA will also be mentioned to promote transfer of ideas. and counterdiabatic (CD)] appear so many times that the use Examples of this are the derivative removal of adiabatic gate of the acronym is well justified. Other acronyms are already (DRAG) and weak anharmonicity with average Hamiltonian of widespread use, such as OCT for optical control theory. To (WAHWAH) approaches to implement fast pulses free from facilitate reading, a full list of acronyms is provided. As for spurious transitions in superconducting qubits (Motzoi et al., notation, an effort toward consistency has been made, so the 2009; Schutjens et al., 2013). symbols often differ from the ones in the original papers. Terminology.—Not surprisingly, the various approaches Beware of possible multiple uses of some symbols (constants, to STA have been described and used with inconsistent in particular). The context and text should help to avoid any terminologies by different groups and communities. Indeed confusion. the rapid growth of STA-related work and its extension across many disciplines has brought up different uses for the same “ ” “ ” words ( superadiabaticity and fast forward are clear II. METHODS examples of polysemic terms), and different expressions for the same concept or method (for example, “transitionless A. Overview of inverse engineering approaches quantum driving” and “counterdiabatic approach”). This review is also intended to clarify some commonly found uses We begin with an overview of inverse engineering. This or expressions, making explicit our preference when the expression refers to inferring the time variation of the control polysemy may lead to confusion. parameters from a chosen evolution of the physical system of Scope and related reviews.—We intend the review to be interest. The notion is broadly applicable for quantum, didactical for the noninitiated but also comprehensive so that classical, or stochastic dynamics and embraces many STA the experts in the different subfields can find a good starting techniques. In this section we deal first with two simple point for exploring other STA-related areas. Several recent examples for motional and internal degrees of freedom, where reviews on partial aspects or overlapping topics are useful more sophisticated or auxiliary concepts and formalisms, e.g., companions: Torrontegui, Ibáñez et al. (2013) is the previous related to invariants, counterdiabatic driving, or fast forward, most comprehensive review on the subject, but the number are not explicitly used. Instead, the connection between of new applications in previously unexplored fields, dynamics and control is performed by solving for the control
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-4 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … function(s) in the effective equation of motion for a classical Section II.H provides a deeper view on the techniques variable or a quantum mean value. developed to enhance the robustness of a given protocol. It is often necessary or useful to go beyond this simple inversion approach for several reasons. For example: the 2. Spin manipulation dynamical path of the system is not easy to design, the inversion is nontrivial, we are interested in finding state- A similar approach can be used for designing the magnetic independent protocols, the direct inversion leads to unrealiz- field components to induce a given trajectory of the mean able control functions, or we look for stable controls, robust value of a spin 1=2 SðtÞ on the Bloch sphere (Berry, 2009), with respect to specific perturbations. To address these 1 problems, in the following sections several auxiliary concepts BðtÞ¼B0ðtÞSðtÞþ SðtÞ × ∂tSðtÞ; ð3Þ and formal superstructures will be added to the simple inverse γ engineering idea in the different STA approaches. The first γ hint of a family of inversion methods dealing with detailed where is the gyromagnetic ratio and B0ðtÞ is any time- dynamics beyond mean values is given in Sec. II.A.3 after the dependent function. This solution is found by inverse engi- two examples. neering of the precession equation for the mean value of the spin, 1. Quantum transport ∂tSðtÞ¼γBðtÞ × SðtÞ: ð4Þ Finding the motion of a harmonic trap to transport a quantum particle so that it starts in the ground state and ends In practice, we usually define the initial and target state and up in the ground state of the displaced potential amounts to build up the solution by interpolation as before. Such an inversely solving a classical Newton equation (Schmiedl et al., approach has been used for instance to manipulate a spin by 2009; Torrontegui et al., 2011). The equation of motion of the defining the time evolution of the spherical angle of the spin particle position x inside the potential reads on the Bloch sphere (Vitanov and Shore, 2015; Zhang, Chen, and Gu´ery-Odelin, 2017). The very same problem can be ̈ ω2 ω2 x þ 0x ¼ 0x0ðtÞ; ð1Þ reformulated using other formalisms such as the Madelung representation (Zhang, Chen, and Gu´ery-Odelin, 2017) that ω where 0 is the angular frequency of the harmonic trap, x0ðtÞ yields an equation of motion that can be readily reversed. is the instantaneous position of its center, and the dots Other formulations of the inverse engineering approach can be represent, here and in the following, time derivatives. made by a proper shaping of the evolution operator (Jing et al., Equation (1) describes a forced oscillator driven by the 2013; Kang et al., 2016c) or by time rescaling (Bernardo, ω2 time-dependent force FðtÞ¼m 0x0ðtÞ. We can interpolate 2019). This method has been applied to two- (Zhang, Chen, the trajectory xðtÞ between the initial position of the particle and Gu´ery-Odelin, 2017), three- (Kang et al., 2016c; Kang, and the desired final position d. In addition, to ensure that the Huang et al., 2017), and four-level systems (Li, Martínez- transport ends up at the lowest energy state, one needs to Cercós et al., 2018). Inverse engineering was also used for 0 cancel out the first and second derivatives at the initial t ¼ open quantum systems (Jing et al., 2013; Impens and Gu´ery- and final time tf (Torrontegui et al., 2011). A simple Odelin, 2017). polynomial interpolation of fifth degree can account for such boundary conditions (Torrontegui et al., 2011), 3. Beyond mean values t 3 t 4 t 5 Suppose now that a more detailed specification of the xðtÞ¼d 10 − 15 þ 6 : ð2Þ tf tf tf dynamics is needed and we focus on closed, linear quantum systems. We shall design the unitary evolution operator UðtÞ ψ Once xðtÞ is defined, Eq. (1) can be easily inverted to give the by specifying a complete basis of dynamical states j jðtÞi corresponding expression for the driving term x0ðtÞ. Note that assumed to satisfy a time-dependent Schrödinger equation there are infinitely many interpolating functions consistent driven by a, to be determined, Hamiltonian, with the boundary conditions at initial and final times. This X freedom is quite typical of different STA methods and can be UðtÞ¼ jψ jðtÞihψ jð0Þj ð5Þ exploited to satisfy other conditions, e.g., minimizing average j energy of a particle during displacement. More parameters can be added in the interpolation functions to minimize the [see other proposals for the form of U by Kang et al., 2016c quantity of interest (Torrontegui et al., 2011). For instance, the and Santos, 2018]. The corresponding Hamiltonian is given robustness against errors in the value of the angular frequency from the assumed dynamics by of the trap ω0 can be enhanced using a Fourier reformulation of the transport problem (Gu´ery-Odelin and Muga, 2014); HðtÞ¼iℏUU_ †: ð6Þ see an experimental application by An et al. (2016) and Sec. VII.B.2. The simple example provided here was gener- As in the previous examples, a typical scenario is that the alized to take into account anharmonicities (Zhang, Chen, and initial and final Hamiltonians are fixed by the experiment or Gu´ery-Odelin, 2015) and in 3D to manipulate Bose-Einstein the intended operation. With these boundary conditions the condensates with an atom chip (Corgier et al., 2018). functions jψ jðtÞi at initial and final times are usually chosen as
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-5 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … the eigenstates of the corresponding Hamiltonians, but there is initially an eigenstate of H0ð0Þ will continue to be so under freedom to interpolate them, and therefore HðtÞ, in between. slow enough driving with the form Several approaches depend on different ways to choose the ξ ψ i nðtÞ orthogonal basis functions j jðtÞi: (a) in the counterdiabatic jψ nðtÞi ¼ e jnðtÞi; ð8Þ driving approach they are instantaneous eigenstates of a reference Hamiltonian H0ðtÞ; (b) in invariant-based engineer- where the adiabatic phases ξnðtÞ are found by substituting ing, they are eigenstates of the invariant of an assumed Eq. (8) as an Ansatz into the time-dependent Schrödinger Hamiltonian form; and (c) more generally they can be just equation driven by H0ðtÞ, convenient functions. In particular, they can be parametrized Z Z 1 t t so that HðtÞ obeys certain constraints, such as making zero 0 0 0 0 0 ξnðtÞ¼− dt Enðt Þþi dt hnðt Þj∂t0 nðt Þi: ð9Þ undesired terms or matrix elements. Note that (a), (b), and (c) ℏ 0 0 are in a certain sense equivalent as they can be reformulated in each other’s language, and they all rely on Eq. (6). For We now seek a Hamiltonian HðtÞ for which the approximate ψ example, (b) or (c) do not explicitly need an H0ðtÞ, and (a) or states j nðtÞi become the exact evolving states, (c) do not explicitly need the invariants, but H0ðtÞ or the ℏ∂ ψ ψ invariants couldP be found if needed. In particular, any linear i tj nðtÞi ¼ HðtÞj nðtÞi: ð10Þ combination cjjψ jðtÞihψ jðtÞj with constant coefficients cj is by construction an invariant of motion of HðtÞ. Inverse HðtÞ is constructed using Eq. (6) from the unitary evolution engineering may also be based on partial information, such as, operator e.g., using a single function of the set and imposing some X ξ additional condition on the Hamiltonian, for example, that the UðtÞ¼ ei nðtÞjnðtÞihnð0Þj; ð11Þ potential is local (i.e., diagonal) and real in coordinate space. n These conditions are the essence of the streamlined version ℏ∂ (Torrontegui, Martínez-Garaot et al., 2012) of the fast-forward which obeys i tUðtÞ¼HðtÞUðtÞ, so that an arbitrary state approach (Masuda and Nakamura, 2008, 2010). The follow- evolves as X ing explores all these approaches in more detail. ξ jψðtÞi ¼ jnðtÞiei nðtÞhnð0Þjψð0Þi: ð12Þ B. Counterdiabatic driving n
The basic idea of counterdiabatic driving is to add auxiliary After substituting Eq. (11) into Eq. (6), or alternatively interactions to some reference Hamiltonian H0ðtÞ so that the differentiating Eq. (12), the Hamiltonian becomes dynamics follows exactly the approximate adiabatic evolution driven by H0ðtÞ. An illustrative analogy is a flat, horizontal HðtÞ¼H0ðtÞþHCDðtÞ; ð13Þ road turn (the reference) that is modified by inclining the X ℏ ∂ roadway surface about its longitudinal axis with a bank angle HCDðtÞ¼i ½j tnðtÞihnðtÞj so that the vehicles can go faster without sliding off the road. n After some precedents (Unanyan et al., 1997; Emmanouilidou − hnðtÞj∂tnðtÞijnðtÞihnðtÞj ; ð14Þ et al., 2000), the CD driving paradigm was worked out and developed systematically by Demirplak and Rice (2003, 2005, where HCD is Hermitian and nondiagonal in the jnðtÞi basis. 2008) for internal state transfer using control fields, then As a simple example of H0 and HCD, consider a two-level rediscovered in a different but equivalent way as “transition- system with reference Hamiltonian ” less tracking by Berry (2009), and used to design many control schemes after Chen, Lizuain et al. (2010), unaware of ℏ −ΔðtÞ ΩRðtÞ ’ H0ðtÞ¼ ; ð15Þ the work of Demirplak and Rice, employed Berry s method to 2 Ω ðtÞ ΔðtÞ control two- and three-level systems. R Berry’s formulation.—We start with Berry’s formulation where ΔðtÞ is the detuning and ΩRðtÞ is the real Rabi because it is somewhat simpler. In Berry (2009), the starting frequency. The counterdiabatic Hamiltonian has the form point is a reference Hamiltonian X ℏ 0 −iΩaðtÞ H0ðtÞ¼ jnðtÞiE ðtÞhnðtÞj: ð7Þ HCDðtÞ¼ ; ð16Þ n 2 iΩ ðtÞ 0 n a Ω Ω Δ_ − Ω_ Δ Ω2 Ω We adopt for simplicity a notation appropriate for a discrete withpffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiaðtÞ¼½ RðtÞ ðtÞ RðtÞ ðtÞ = ðtÞ and ðtÞ¼ 1 0 Δ2 Ω2 (real) spectrum and no degeneracies. A state jnð Þi that is ðtÞþ RðtÞ. See further simple examples in Appendix B. Returning to the general expression (13), we note that HCD † 1Generalizations for degenerate levels, relevant, for example, to is orthogonal to H0 considering the scalar product trðA BÞ 0 speed up holonomic quantum gates, may be found by Takahashi of two operators. Using dhnðtÞjmðtÞi=dt ¼ it can be seen _ (2013b) and Zhang et al. (2015) or Karzig et al. (2015). General- that HCD is also orthogonal to H0 (Petiziol et al., 2018). izations for non-Hermitian Hamiltonians H0ðtÞ are discussed in An alternative form for HCD is found by differentiating Sec. V.D. H0ðtÞjnðtÞi ¼ EnðtÞjnðtÞi,
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-6 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
X X _ ℏ jmðtÞihmðtÞjH0jnðtÞihnðtÞj HCD ¼ i − ; ð17Þ m≠n n En Em which, using the scaled time s ¼ t=tf, gives the scaling HCD ∼ 1=tf. In general HCDðtÞ vanishes for t<0 and t>tf, either suddenly or continuously at the extreme times, so that the jnðtÞi become eigenstates of the full Hamiltonian at 0− þ the boundary times t ¼ and tf . In terms of the phases the total Hamiltonian (13) can be written as (Chen, Torrontegui, and Muga, 2011) X HðtÞ¼JðtÞþiℏ j∂tnðtÞihnðtÞj; ð18Þ n with X _ JðtÞ¼−ℏ jnðtÞiξnðtÞhnðtÞj: ð19Þ n
Subtracting H − HCD we get an alternative form of H0 consistent with Eqs. (7) and (9), X _ H0 ¼ ℏ jnðtÞi½−ξn þ ihnðtÞjn_ðtÞi hnðtÞj: ð20Þ n
If we write the same H0ðtÞ [Eq. (7)] as beforeP in an alternative 0 0 basis with different phases H0ðtÞ¼ njn ðtÞiEnðtÞhn ðtÞj, ϕ where jn0ðtÞi ¼ ei n jnðtÞi, the formalism goes through using 0 primed functions ξn and corresponding primed operators. ξ 0 nðtÞ 0 FIG. 3. Scheme for superadiabatic iterative interaction pictures. [One example is jn ðtÞi ¼ e jnðtÞi so that ξnðtÞ¼0.] It is easy to check though that the terms compensate so that At the end of each iteration it is possible, apart from starting a 0 new one, to either neglect the nondiagonal coupling in H 1 H ðtÞ¼H ðtÞ, i.e., a change of representation for H0ðtÞ jþ CD CD 0 does not change HðtÞ, nor the CD-driving term, nor the (adiabatic approximation if j ¼ or superadiabatic approxima- j ≥ 1 j 0 physics. tion for ) or add a term (CD for ¼ or super-CD otherwise) to the Hamiltonian H 1 to exactly cancel the Quite a different issue is to change the physics when jþ coupling. imposing a set of basis functions fjnðtÞig and a set of phases fξnðtÞg which are not a priori regarded as adiabatic (Chen, Torrontegui, and Muga, 2011). This procedure is essentially The formulation of Demirplak and Rice.—To follow some invariant-based inverse engineering since one is imposing past and recent developments and generalizations it is worth some specific dynamics (thus some invariants) without pre- finding the auxiliary Hamiltonian HCD by the equivalent supposing a given H0ðtÞ. UðtÞ in Eq. (11) becomes the formulation of Demirplak and Rice (2003, 2005, 2008), which primary object and the driving Hamiltonian is given by is the zeroth iteration j ¼ 0 in Fig. 3. Eq. (18) with diagonal part (20), which defines H0ðtÞ, and In the following discussion we concentrate on j ¼ 0. a coupling, nondiagonal part (14). Now, changing the phases Among the possible phase choices for the eigenstates of ξnðtÞ modifies H0ðtÞ and has an impact on the physical H0ðtÞ, we now use the one that satisfies the “parallel evolution of a general wave function jψðtÞi. Of course the transport” condition hn0ðtÞjn_0ð0Þi ¼ 0. Regardless of the 2 “ ” populations PnðtÞ ≡ jhnðtÞjψðtÞij ¼ Pnð0Þ driven by working basis of eigenvectors jnðtÞi one starts with, the _ Eq. (18) are not affected by phase shifts. The ξnðtÞ can be parallel transported basis is found as optimized to minimize energy costs (Hu et al., 2018), robust- ness against decoherence (Santos and Sarandy, 2018), or R − thnðt0Þjn_idt0 intensity of the extra CD controls (Santos et al., 2019). jn0ðtÞi ¼ e 0 jnðtÞi: ð21Þ A related procedure is to set the jnðtÞi but consider H0ðtÞ and the eigenvalues EnðtÞ controllable elements rather than Later on we shall see the consequences of applying or not given (Berry, 2009). For example, setting all EnðtÞ¼0 cancels the dynamical phase. That means that for a given applying parallel transport. A dynamical solution driven by H0ðtÞ is written as jψ 0ðtÞi and the instantaneous eigenvalues set of states jnðtÞi, HCDðtÞ alone (without an H0) drives the ð0Þ same populations than HðtÞ¼H0ðtÞþHCDðtÞ for any choice as En ðtÞ. The extra notational burden of the index j ¼ 0 may of EnðtÞ (Chen, Lizuain et al., 2010). be ignored in many applications, in particular, for ordinary CD
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-7 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … driving, but it will be useful to cope with higher order ½n ℏ HCD ¼ i ½dPn=dt; Pn −; ð25Þ suparadiabatic iterations later on. The way to find the auxiliary driving term is to first express ½n which fulfils PmHCDPn ¼ PmHCDPn and PnHCDPm ¼ the dynamics in an adiabatic frame. For that we apply the ½n ½n † 0 0 ψ ψ PnHCDPm for all m, as well as PmHCDPm ¼ for m, transformation j 1ðtÞi ¼ A0ðtÞj 0ðPtÞi to the dynamical states 0 ½n 0 m ≠ n, uncouples the dynamics of level n. For H0 H , driven by H0ðtÞ, with A0ðtÞ¼ njn0ðtÞihn0ð Þj. In the þ CD ξ resulting adiabatic frame the coupling, diabatic terms ei nðtÞjnðtÞi is an exact solution of the dynamics. This is an − † ℏ _ † interesting simplification as we are often interested in only one Pare made obvious, A0K0A0, where K0 ¼ i A0A0 ¼ _ state, typically the ground state. The last condition imposed by njn0ðtÞihn0ðtÞj. In this adiabatic frame we may (a) neglect ½n 0 0 0 ≠ the coupling (adiabatic approximation), or (b) cancel the Demirplak and Rice (2008), PmHCDPm ¼ for m, m n,is † coupling by adding A0K0A to the Hamiltonian; see again not really necessary for the uncoupling of the nth level. Fig. 3. In the original, Schrödinger picture, also referred to Without it, a broad set of “state-dependent” CD operators ½n 1 − as the laboratory frame hereafter, this addition amounts to HCD þ QnBQn, where Qn ¼ Pn and B is any Hermitian using H ¼ H0 þ K0. operator, can be generated. This multiplicity may be useful An alternative useful form for K0ðtÞ is (Messiah, 1962) and explains why different auxiliary state-dependent CD terms have been proposed (Patra and Jarzynski, 2017b; X dPnðtÞ Setiawan et al., 2017). K0ðtÞ¼iℏ PnðtÞ; ð22Þ n dt 1. Superadiabatic iterations where PnðtÞ¼jn0ðtÞihn0ðtÞj ¼ jnðtÞihnðtÞj. K0 may Let us recap before moving ahead. The adiabatic interaction as well be written in an arbitrary basis of eigenvectors picture (IP) corresponds to expressing the quantum dynamics jnðPtÞi, i.e., not necessarily parallel transported, K0ðtÞ¼ in the adiabatic basis of instantaneous eigenstates of H0ðtÞ. ℏ _ − _ i njnðtÞihnðtÞj jnðtÞihnðtÞjnihnðtÞj ¼ HCDðtÞ. This The dynamical equation in the adiabatic IP includes an result proves that K0 is identical to HCD in Eq. (14). effective Hamiltonian H1ðtÞ with a diagonal (adiabatic) term Since the projectors PnðtÞ¼jn0ðtÞihn0ðtÞj ¼ jnðtÞihnðtÞj and a coupling term; see Fig. 3. are invariant with respect to phase choices of the eigenvectors, We can repeat the sequence iteratively. In the first “super- ” the form in which HCD is usually written in Eq. (14) is adiabatic iteration, H1 is diagonalized to find its instanta- invariant under different choices of phases for the eigenstates. neous (parallel-transported) eigenstates jn1ðtÞi. With the new “ ” In particular, this implies that K0 is purely nondiagonal superadiabatic basis a new IP is generated driven by an [Kato’s condition (Kato, 1950; Demirplak and Rice, 2008)] effective Hamiltonian H2 with a diagonal part in the super- in the arbitrarily chosen basis of eigenvectors of H0ðtÞ, jnðtÞi. adiabatic basis and a coupling term. This new coupling term † If we now return to the parallel transported basis in the −A1K1A1 may be (a) neglected (superadiabatic approxima- † interaction (adiabatic) picture, the addition of A0K0A0 cancels tion) or (b) canceled by adding its negative, and so on. The H0 the couplings so the dynamics is trivially solved with canceling term to be added to in the SchrödingerQ picture ðjÞ † 1 j−1 ≥ 1 dynamical phase factors. In the Schrödinger picture, is Hcd ¼ BjKjBj , with B0 ¼ and Bj ¼ k¼0 Ak for j . R Note that in general only the one generated in the zeroth − ℏ t 0 0 0 ði= Þ Enðt Þdt iteration agrees with the standard CD term H ¼ H . jψ 0ðtÞi ¼ e 0 jn0ðtÞihn0ð0Þjψ 0ð0Þi; ð23Þ cd CD The recursive iterations were worked out by Garrido (1964), without considering the cancellations, to find out which is exactly Eq. (12) as can be seen by using Eq. (21). ˜ generalizations of the adiabatic approximation. Berry (1987) PIf, instead of A0, a more general transformation A ¼ 0 also used them to calculate a sequence of corrections to njnðtÞihnð Þj is used, the coupling term in the interaction Berry’s phase and introduced the name “superadiabatic trans- − ˜ † ˜ ˜ picture becomes A0K0A0, where formations.” Demirplak and Rice (2008) proposed to apply the superadiabatic iterative frame to generate alternative (to the ˜ ˜_ ˜ † K0 ¼ iℏA0A0 simple CD approach) higher order coupling-canceling terms. X Later Ibáñez, Chen, and Muga (2013) made explicit the _ ℏ _ ¼ jnðtÞihnðtÞj ¼ K0 þ i jnðtÞihnðtÞjnðtÞihnðtÞj conditions that the derivatives of H0ðtÞ must satisfy at the n time boundaries in order to really generate a shortcut to adiabaticity (rather than just a shortcut to superadiabaticity), and its cancellation would lead in the laboratory picture to a i.e., a protocol that takes instantaneous eigenstates of different state, H0ðt ¼ 0Þ to corresponding eigenstates of H0ðtfÞ. X R The naive expectation that each iteration will produce −ði=ℏÞ t E ðt0Þdt0 e 0 n jnðtÞihnð0Þjψ 0ð0Þi; ð24Þ smaller and smaller couplings does not hold in general. n They decrease up to an optimal iteration and then grow (Berry, 1987). Working with the optimal frame may or may although the probabilities are not affected. not be worthwhile depending on whether the boundary Demirplak and Rice (2008) also pointed out that the conditions for derivatives of H0ðtÞ are fulfilled (Ibáñez, operator Chen, and Muga, 2013).
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-8 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
An interesting feature of the superadiabatic sequence of a. “Physical” unitary transformations coupling terms is that their operator structure changes with A useful method to generate alternative, physically feasible the iteration. For example, for a two-level system with shortcuts from an existing shortcut generated by counter- σ σ σ Hamiltonian X x þ ZðtÞ z, X is constant, where the x;y;z diabatic driving or otherwise, is to perform physical, rather are Pauli matrices, the first (adiabatic) CD term K0 is of the than formal, unitary transformations (Ibáñez et al., 2012), σ form YðtÞ y, see Eq. (16), whereas the second (first-order or corresponding canonical transformations in classical sys- superadiabatic) coupling term in the Schrödinger frame tems (Deffner, Jarzynski, and del Campo, 2014). Given a reproduces the structure of H0 with x and z components Hamiltonian HðtÞ that drives the wave function jψðtÞi, the but not a y component. [For three-level systems see Huang unitarily transformed state jψ 0ðtÞi ¼ U†ðtÞjψðtÞi is driven by et al. (2016), Kang et al. (2016b), Song, Ai et al. (2016), the Hamiltonian (the primes here distinguish the picture, they and Wu, Su, Ji, and Zhang (2017).] Unitarily transforming K0 do not represent derivatives) (see Sec. II.B.2) also provides a Hamiltonian without a 1 ð Þ 0 † y component, which is different from the term Hcd ¼ H ðtÞ¼U ðH − KÞU; ð26Þ † A0K1A0 generated from the first superadiabatic iteration, see an explicit comparison by Ibáñez et al. (2012), where K ¼ iℏUU_ †: ð27Þ ð1Þ Hcd was shown to have smaller intensity than HCD. To avoid confusion we discourage the use of the expression If we set Uð0Þ¼UðtfÞ¼1, then the wave functions coincide 0 superadiabatic to refer to the regular CD approach (in fact a at the boundary times jψð0Þi¼jψ ð0Þi and jψðtfÞi ¼ zeroth order in the superadiabatic iterative frame) or its unitarily jψ 0ðt Þi. If, in addition, U_ð0Þ¼U_ðt Þ¼0, then also the transformed versions discussed in Sec. II.B.2. This use of the f f Hamiltonians coincide at boundary times Hð0Þ¼H0ð0Þ word superadiabatic, as being equivalent to CD driving or even and Hðt Þ¼H0ðt Þ. generically to all shortcuts is, however, somewhat extended. f f Note that here the alternative Hamiltonian form (26) and the DRAG controls.—The DRAG framework was developed to state jψ 0ðtÞi are not just convenient mathematical transforms avoid diabatic transitions to undesired levels in the context of of HðtÞ or jψðtÞi representing the same physics as in superconducting quantum devices. It was recently reviewed conventional interaction or Heisenberg picture transforma- by Theis et al. (2018) but we sketch the main ideas here since tions. Instead, at intermediary times, H0ðtÞ and HðtÞ represent it is related to the superadiabatic scheme as formulated, e.g., indeed different (laboratory) drivings, and ψ t and ψ 0 t by Ibáñez, Chen, and Muga (2013). A typical scenario is that j ð Þi j ð Þi ðjÞ † different dynamical states. This was emphasized by S. Ibáñez the (super)-CD term Hcd ¼ BjKjBj is not physically feasible et al. (2012) by calling the different alternatives “multiple and does not match the controls in the lab. The DRAG Schrödinger pictures.” approach addresses this problem by decomposing the con- The art is to find a useful UðtÞ to make H0ðtÞ feasible. This Ptrollable Hamiltonian that “corrects” the dynamics as Hctrl ¼ approach was applied in many works; see, e.g., Hollenberg kukðtÞhk þ H:c: with control fields ukðtÞ and coupling (2012), Takahashi (2015), Agundez et al. (2017), and Sels and terms hk. In a given superadiabatic frame partial contributions Polkovnikov (2017), and several experiments (Bason et al., ðjÞ 2012; Zhang et al., 2013; An et al., 2016; Du et al., 2016). to ukðtÞ are found by projecting Hcd into the assumed Hamiltonian structure. This generates an approximation to Deffner, Jarzynski, and del Campo (2014) applied it to the exact uncoupling term so the diabatic coupling, even if not generate feasible (i.e., involving local potentials in coordinate canceled exactly, is reduced. An important point is that further space, independent of momentum) Hamiltonians for scale- iterations, in contrast with the bare superadiabatic iterations, invariant dynamical processes. typically converge, so that the coupling eventually vanishes. When HðtÞ is a linear combination of generators Ga of A variant of this approach applies when the error terms are not some Lie algebra, independently controlled, because they all depend on some common control, e.g., a single laser field. Different perturba- XN α tive approximations using a power series in the inverse gap ½Gb;Gc ¼ abcGa; ð28Þ 1 energies were worked out systematically (Theis et al., 2018). a¼
where the αabc are the structure constants, U may be con- 2. Beyond the basic formalism structed by exponentiating elements of the algebra and imposing the vanishing of the unwanted terms (Martínez- The CD Hamiltonian HCD often implies different operators from those in H0 that typically may be hard or even impossible Garaot et al., 2014). To carry out the transformation an to generate in the laboratory. Moreover, a lot of spectral element G of the Lie algebra of the Hamiltonian is chosen, information is in principle used to build HCD,specifically −igðtÞG the eigenvectors of H0, so a number of strategies, reviewed UðtÞ¼e ; ð29Þ hereafter, are put forward to avoid some terms in the auxiliary Hamiltonian, and/or the spectral information needed. Changing where gðtÞ is a real function to be set. This type of unitary the phases ξnðtÞ without modifying the jnðtÞi changes H0 but operator UðtÞ constitutes a “Lie transform.” Note that K in _ 0 not HCD, so it is not enough for these purposes (Ibáñez, Eq. (27) becomes −ℏgðtÞG and commutes with G so H is Martínez-Garaot et al.,2011). given by
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-9 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
U†ðH − KÞU ¼ eigGðH − KÞe−igG b. Schemes that focus on one state 2 A simplifying assumption that helps one to find simpler − ℏ _ − g ¼ H gG þ ig½G; H 2! ½G; ½G; H decoupling terms is to focus on only one state, typically the 3 ground state j0ðtÞi or a subset of states. Explicit exact forms of g − the driving that uncouples that state were worked out quite early i 3! ½G; ½G; ½G; H þ ; ð30Þ (Demirplak and Rice, 2008), see Eq. (25), and more recently for systems described in coordinate space in Patra and Jarzynski which depends only on G, H, and its nested commutators with (2017b) or in a discrete basis (Setiawan et al., 2017). As for G, so it stays in the algebra. If we can choose G and gðtÞ so approximate schemes, Opatrný and MølmerP (2014) proposed K that the undesired generator components in HðtÞ cancel out to add to H0ðtÞ the Hamiltonian HBðtÞ¼ k fkðtÞTk using and the boundary conditions for U are satisfied, the method only feasible interactions (i.e., available experimentally) Tk, provides a feasible, alternative shortcut. A simple example for where the fkðtÞ are amplitudes found by minimizing the norm − 0 the two-level Hamiltonian is given in Appendix A. Martínez- of ðHB HCDÞj ðtÞi. A related approach uses the Lyapunov Garaot et al. (2014) and Kang et al. (2018) provided more control theory (Ran et al., 2017). Similarly, Chen, Xia et al. examples. (2016) and Chen et al. (2017) achieved feasible auxiliary Interaction pictures.—Frequently the CD terms are Hamiltonians by adding adjustable Hamiltonians that nullify found formally in a transformed picture, I in Fig. 4, which unwanted nonadiabatic couplings for specific transitions. is intended only as a mathematical aid. The effective c. Effective counterdiabatic field Hamiltonian HI in this picture, without the CD term, P represents the same physics as some original Hamiltonian In Petiziol et al. (2018, 2019) H0ðtÞ¼ kukðtÞHk is H in the Schrödinger picture S in Fig. 4. Once the CD term is S written as before in terms of control functions ukðtÞ and added we get a new Hamiltonian H0 that when transformed I available time-independent control Hamiltonians Hk. Using back gives H 0 . It happens often that simple-looking auxiliary S control theory arguments it is found that HCD necessarily terms in the transformed frames become difficult to implement belongs to the corresponding dynamical Lie algebra, i.e., the in the laboratory frame. In Sec. III.G we comment on some N- smallest algebra that contains the −iHk and the nested body models with easily solvable CD terms in a convenient commutators. As well, the action of HCD can be approximated transformed frame but hard to realize in the lab frame. by using only the initially available HamiltoniansP in an Alternative unitary transformations Uφ˜ , different from the “effective counterdiabatic” Hamiltonian HE ¼ ekðtÞHk. one used to go between the I and S pictures Uφ, may help to To implement HE in Petiziol et al. (2018 , 2019), the solve the problem, giving from HI0 feasible shortcuts driven control functions ekðtÞ are chosen as periodic functions of by a new lab frame Hamiltonian HS00 . Ibáñez et al. (2015) Pperiod T with the form of a truncated Fourier expansion A ω B ω ω 2π worked out this alternative route for two-level systems when k k sinðk tÞþ k cosðk tÞ, where ¼ =T. the rotating wave approximation is not applicable; see also The coefficients are determined by setting the first terms Ibáñez, Peralta Conde et al. (2011), Chen and Wei (2015), of the Magnus expansion generated by HE to match those of and H. Li et al. (2017). Physics beyond the rotating wave the desired evolution stroboscopically at multiples of T and approximation is of much current interest due to the increasing interpolating smoothly in between. Avoided crossing prob- use of strong fields and microwave frequencies, for example, lems and entanglement creation are addressed with this in NV centers. technique. Note alternative uses of the Magnus expansion in the WAHWAH technique of Schutjens et al. (2013), aimed at producing fast pulses to operate on a qubit without interfering with other qubits in frequency-crowded systems, and in Claeys et al. (2019).
d. Dressed-states approach CD driving is generalized by Baksic, Ribeiro, and Clerk (2016) by considering a dressed-states approach which uses three different dynamical pictures. In this summary the notation and even the terminology differs from Baksic, Ribeiro, and Clerk (2016); see Table I. First consider a Schrödinger picture description where the driving Hamiltonian is H0ðtÞþHcðtÞ with reference Hamiltonian H0ðtÞ as in Eq. (7), and driven wave function ψ t . H t will be found later so that the total Hamiltonian FIG. 4. Schematic relation between different Schrödinger and j ð Þi cð Þ interaction pictures. Each node corresponds also to different satisfies some conditions. Then a first rotating picture defined ψ † ψ Hamiltonians. The rectangular boxes enclose nodes that represent byP j IðtÞi ¼ UðtÞ j ðtÞi is introduced, where UðtÞ¼ the same underlying physics. The solid lines represent unitary njnðtÞihnj and the jni are time independent. A simple 0 0 1 relations for the linked states and the dashed line represents a choice is jni¼jnð Þi that insures Uð Þ¼UðtfÞ¼ . Since nonunitary addition of an auxiliary term to the Hamiltonian. jnðtÞi are instantaneous (adiabatic) eigenstates of H0ðtÞ
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-10 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
TABLE I. Scheme for “dressed-state” driving (Baksic, Ribeiro, and Clerk, 2016).
Picture Wave function Unitary transformation Hamiltonian
Schrödinger ψðtÞ P H ¼ H0 þ Hc “ ” ψ †ψ _ † † First rotating picture ( adiabatic ) IðtÞ¼U ðtÞ U ¼ PnjnðtÞihnj HI ¼ iℏU U þ U HU “ ” ψ †ψ ˜ _ † † Second rotating picture ( dressed ) IIðtÞ¼V IðtÞ V ¼ njnðtÞihnj HII ¼ iℏV V þ V HIV
“ ” _ we may naturally call this picture the adiabatic frame with HðtÞ¼H0 þ λAλ; ð33Þ driving Hamiltonian HIðtÞ; see Table I. Then a second rotating ψ † ψ ˜ † † frameP is introduced as j IIðtÞi ¼ VðtÞ j IðtÞi, with VðtÞ¼ where Aλ ¼ UAλU ¼ iℏð∂λUÞU such that in the moving ˜ ˜ eff ˜ njnðtÞihnj, driven by HIIðtÞ. fjnðtÞig are called dressed frame H ¼ H0 is diagonal and no transitions are allowed. 0 1 0 states and it is assumed that Vð Þ¼VðtfÞ¼ . The method to Up to now we have only introduced a new notation and ˜ generate shortcuts is to choose VðtÞ, i.e., the functions jnðtÞi terminology in Eqs. (31)–(33).2 In the adiabatic limit jλ_j → 0 and H so that H ðtÞ is diagonal in the basis fjnig, c P II the Hamiltonian (33) reduces to the original one HðtÞ¼H0ðtÞ. II HIIðtÞ¼ njniEn ðtÞhnj. Back to the Schrödinger picture After differentiating the first equation in Eq. (32), the gauge this means that no transitions occur among states potential satisfies (Jarzynski, 2013) UðtÞVðtÞjni¼UðtÞjn˜ðtÞi. VðtÞ is chosen to ensure that HcðtÞ is feasible. The two iℏð∂λH0 þ FadÞ¼½Aλ;H0 ; ð34Þ nested transformations make the method more involved and P less intuitive than standard CD driving. To gain insight note − ∂ λ λ λ where Fad ¼ n λEnð Þjnð Þihnð Þj is the adiabatic force that if VðtÞ¼1 for all t the method reduces to standard CD operator. Since, by construction, Fad commutes with H0, driving. One can think of V as a way to add flexibility to the Eq. (34) implies that inverse engineering so as to drive states in the Schrödinger picture along uncoupled UðtÞVðtÞjni vectors rather than along ½iℏ∂λH0 − ½Aλ;H0 ;H0 ¼0; ð35Þ vectors UðtÞjni¼jnðtÞi. The latter are given, up to phases, once H0ðtÞ is specified, whereas UðtÞVðtÞjni may still be where the difficult-to-calculate force has been eliminated. This manipulated to find a convenient HcðtÞ and possibly minimize equation can be used to find the adiabatic gauge potentials the occupancy of some state to be avoided, e.g., because of directly without diagonalizing the Hamiltonian. Moreover, spontaneous decay (Baksic, Ribeiro, and Clerk, 2016). solving this equation is analogous to minimizing the Hilbert- For applications see Baksic et al. (2017), Coto et al. (2017), Schmidt norm of the operator Liu et al. (2017), Wu, Ji, and Zhang (2017b), B. B. Zhou et al. (2017), and X. Zhou et al. (2017). i ∂ Gλ ¼ λH0 þ ℏ ½Aλ;H0 ð36Þ e. Variational approach with respect to Aλ, where Aλ is a trial gauge potential. For an Motivated by difficulties to diagonalize H0ðtÞ and the application to prepare ground states in a lattice gauge model, nonlocalities in the exact counterdiabatic Hamiltonian in see Hartmann and Lechner (2018). many-body systems, Kolodrubetz et al. (2017) and Sels Potential problems with this scheme are that it may be and Polkovnikov (2017) developed a variational method to difficult to know what local operators should be included in construct approximate counterdiabatic Hamiltonians without the variational basis, and moreover their practical realizability using spectral information. The starting point in Sels and is not guaranteed. To solve these two difficulties, Claeys et al. λ PPolkovnikov (2017) is a unitary transformation U½ ðtÞ ¼ (2019) expanded the gauge potential in terms of nested λ ψ njn½ ðtÞ ihnj to rotate the state j ðtÞi that evolves under a commutators of the Hamiltonian H0 and the driving term λ time-dependent Hamiltonian H0½ ðtÞ , to the moving frame ∂λH0, even if these do not close a Lie algebra. Since the † state jψ˜ ðtÞi ¼ U ½λðtÞ jψðtÞi, which satisfies the effective commutators also arise in the Magnus expansion in Floquet Schrödinger equation systems, they can be realized up to arbitrary order using Floquet engineering, i.e., periodical driving (Boyers et al., ℏ∂ ψ˜ ˜ λ − λ_ ˜ ψ˜ i tj i¼fH0½ ðtÞ Aλgj i; ð31Þ 2018). The expansion coefficients can be calculated either ˜ ˜ analytically or variationally. The method can be adapted where H0 is diagonal in the jni basis and Aλ is the adiabatic to suppress excitations in a known frequency window. “gauge potential” in the moving frame, X Application examples were provided for a two- and a ˜ † three-level system, where one, respectively, two terms H0½λðtÞ ¼ U H0½λðtÞ U ¼ EnðλÞjnihnj; n ˜ † 2 Aλ ¼ iℏU ∂λU: ð32Þ We assumed for simplicity a dependence on time via a single parameter λ. A more general dependence on a vector λ is worked out All nonadiabatic transitions are produced by the gauge by Deffner, Jarzynski, and del Campo (2014) and Nishimura and potential. In the counterdiabatic approach the system is driven Takahashi (2018) which leads to a decomposition of the CD term and by the Hamiltonian a “zero curvature condition” among the CD components.
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-11 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … returned the exact gauge potential, and to a many-body spin general solution to the Schrödinger equation can be chain, where a limited number of terms and an approximate written as gauge potential resulted in a drastic increase in fidelity. X Further useful properties of the expansion are that it relates jψðtÞi ¼ cnjψ nðtÞi; ð38Þ the locality of the gauge potential to the order in the n expansion, and that it remains well defined in the thermody- namic and classical limits. where the cn are independent of time. These invariants were originally used to solve for the state f. Counterdiabatic Born-Oppenheimer dynamics driven by a known time-dependent Hamiltonian (Lewis and Riesenfeld, 1969). In shortcuts to adiabaticity this idea is Duncan and del Campo (2018) proposed to exploit the reversed (Chen, Ruschhaupt et al., 2010) and the Hamiltonian separation of fast and slow variables using the Born- is found from a prescribed state evolution. Formally, a time Oppenheimer approximation. Simpler CD terms (compared evolution operator of the form to the exact CD term) for the fast and slow variables can be X found in two steps avoiding the diagonalization of the full α U ¼ ei nðtÞjϕ ðtÞihϕ ð0Þj ð39Þ Hamiltonian. The method is tested for two coupled harmonic n n n oscillators and a system of two charged particles. implies a Hamiltonian HðtÞ¼iℏUU_ †. This is the essence of g. Constant CD-term approximation invariant-based inverse engineering. Oh and Kais (2014) proposed, in the context of the adiabatic For a given Hamiltonian there are many possible invariants. Grover’s search algorithm implemented by a two-level system, For example, the density operator describing the evolution of to substitute the exact CD term by a constant term. This method a system is a dynamical invariant. The choice of which improves the scaling of the nonadiabatic transitions with particular invariant to use is made on the basis of mathematical respect to the running time. J. Zhang et al. (2018) performed convenience. Invariants have also been generalized to non- an experiment with a single trapped ion choosing the constant Hermitian invariants and Hamiltonians (Gao, Xu, and Qian, term with the aid of a numerical simulation. 1992; Lohe, 2009; Ibáñez, Martínez-Garaot et al., 2011), as In a different vein Santos and Sarandy (2018) discussed the well as open systems (see Sec. V). conditions for H0ðtÞ and phase choice necessary in the unitary A connection to adiabaticity is asymptotic. In the limit of evolution operator (11) to implement an exact constant long operation times (i.e., adiabatic), Eq. (37) becomes driving. ½H; I ≈ 0 and the dynamics prescribed by invariant-based engineering will approach adiabatic dynamics driven by H. C. Invariants and scaling laws Hence for long times, HðtÞ and IðtÞ have approximately a common eigenbasis for all times. Dynamical invariants and invariant-based engineering con- Other relations to adiabaticity do not need long times nor stitute a major route to design STA protocols. The basic common bases for HðtÞ and IðtÞ at all times. In particular, reason, already sketched in Eqs. (5) and (6), is that, in linear demanding only that the invariant and the Hamiltonian systems, determining the desired dynamics amounts to setting commute at the start and the end of the process, i.e., dynamical invariants of motion, and the Hamiltonian may in ½Ið0Þ;Hð0Þ ¼ ½IðtfÞ;HðtfÞ ¼ 0, the eigenstates of the principle be found from them. This connection is quite general invariant and the Hamiltonian coincide at initial and final and is also valid classically (Lewis and Leach, 1982), but it is times, but may differ with each other at intermediate times mostly applied for systems in which the Hamiltonian form and because the commutativity is not imposed. If no level cross- corresponding invariants are known explicitly as functions of ings take place, the final state will keep the initial populations auxiliary parameters that satisfy auxiliary equations consistent for each nth level, as in an adiabatic process, but in a finite with the dynamical equation and the invariant. (short, faster-than-adiabatic) time. This leaves freedom to choose how the state evolves in the intermediate time and then 1. Lewis-Riesenfeld invariants use Eq. (37) to find the Hamiltonian that drives such a state Originally proposed in 1969, a Lewis-Riesenfeld invariant evolution. The flexibility can be exploited to improve the (Lewis and Riesenfeld, 1969) for a Hamiltonian HðtÞ is a stability of the schemes against noise and systematic errors; Hermitian operator IðtÞ which satisfies see Sec. II.H. The relation of the invariant-based approach to CD driving dI ∂I i (Chen, Torrontegui, and Muga, 2011) implies a different ¼ þ ½H; I ¼0; ð37Þ dt ∂t ℏ connection to adiabaticity. So far we have not mentioned in this section any “reference” H0ðtÞ, but by reinterpreting so that the expectation values for states driven by HðtÞ jϕnðtÞi as eigenvectors of H0ðtÞ, and the Lewis-Riesenfeld are constant in time. Since IðtÞ is a constant of motion it phases as adiabatic phases, an implicit H0ðtÞ operator may be ϕ has time-independent eigenvalues. If j nðtÞi is an instanta- written down using Eq. (20). Hence an implicit HCDðtÞ neous eigenstate of IðtÞ, a solution of the Schrödinger follows from subtraction, H ðtÞ¼HðtÞ−H0ðtÞ. In invariant- ℏ∂ ψ ψ CD equation i tj nðtÞi ¼ HðtÞj nðtÞi can beR constructed as based engineering, however, these implicit operators are not iα ðtÞ t jψ nðtÞi ¼ e n jϕnðtÞi. Here αnðtÞ¼ð1=ℏÞ 0 hϕnðsÞj½iℏ∂s− really used and do not play any role in practice. As a HðsÞ jϕnðsÞids is the Lewis-Riesenfeld phase. Hence a consequence of the different emphases and construction
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-12 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … recipes for HðtÞ, given some initial and final Hamiltonians, 2. Examples of invariant-based inverse engineering STA protocols designed via CD or invariant approaches are Lewis-Riesenfeld invariants have been used to design state often very different. transfer schemes in two-level (Ruschhaupt et al., 2012; — Lie algebras. As first noted by Sarandy, Duzzioni, and Martínez-Garaot et al., 2013; Kiely and Ruschhaupt, 2014), Serra (2011), invariant-based inverse engineering can also three-level (Chen and Muga, 2012a; Kiely and Ruschhaupt, be formulated in terms of dynamical Lie algebras. The 2014; Benseny et al., 2017), and four-level systems assumption of a Lie structure is also used to classify and (G¨ungörd¨u et al., 2012; Herrera et al., 2014; Kiely et al., construct dynamical invariants for four-level or smaller non- 2016), or Hamiltonians quadratic in creation and annihilation trivial systems for specific applications (G¨ungörd¨u et al., operators (Stefanatos and Paspalakis, 2018a). They are also 2012; Herrera et al., 2014; Kiely et al., 2016). The approaches useful to design motional dynamics in harmonic traps or by Martínez-Garaot et al. (2014) and Petiziol et al. (2018), otherwise. Many more applications specific for different already discussed, also make use of the Lie algebraic structure. system types may be found in Sec. III. Here we provide Torrontegui, Martínez-Garaot, and Muga (2014) provided two basic examples, the two-level model and the Lewis-Leach a bottom-up construction procedure of the Hamiltonian family of potentials for a particle of mass m moving in 1D. using invariants and Lie algebras. Let us assume that the Hamiltonian of a system HðtÞ and the invariant IðtÞ can be a. Two-level system written as linear combinations of Hermitian operators Ga (generators), We now present an example of a two-level system with a Hamiltonian given by XN XN HðtÞ¼ haðtÞGa;IðtÞ¼ faðtÞGa; ð40Þ ℏ −ΔðtÞ ΩRðtÞ − iΩIðtÞ a¼1 a¼1 HðtÞ¼2 ; ð44Þ ΩRðtÞþiΩIðtÞ ΔðtÞ that form a Lie algebra closed under commutation; see and an invariant of the form Eq. (28). Inserting these forms into Eq. (37), we get that ℏ cos ½θðtÞ sin ½θðtÞ e−iαðtÞ XN IðtÞ¼ : ð45Þ _ 2 sin ½θðtÞ eiαðtÞ − cos ½θðtÞ faðtÞ − AabðtÞhbðtÞ¼0; ð41Þ b¼1 From the equation that defines the invariant, Eq. (37), θðtÞ and α t must satisfy where the N × N matrix A is defined by ð Þ
θ_ ¼ Ω cos α − Ω sin α; ð46Þ 1 XN I R A ≡ α f ðtÞ; ð42Þ ab iℏ abc c c¼1 α_ ¼ −Δ − cot θðΩR cos α þ ΩI sin αÞ: ð47Þ
The eigenvectors of IðtÞ are and the αabc are the structure constants. In vector form, cos θ=2 e−iα=2 ⃗ ⃗ ϕ ð Þ ∂tfðtÞ¼AhðtÞ; ð43Þ j þðtÞi ¼ ; ð48Þ sin ðθ=2Þeiα=2 ⃗ ⃗ where the vectors fðtÞ and hðtÞ represent the invariant and sin ðθ=2Þe−iα=2 Hamiltonian, respectively. The inversion trick is to first choose jϕ−ðtÞi ¼ ; ð49Þ − cos ðθ=2Þeiα=2 the auxiliary functions f⃗ðtÞ (and therefore the state evolution) ⃗ and infer hðtÞ from this. Technically the inversion requires with eigenvalues ℏ=2. The angles α and θ can be thought introducing a projector Q for the null subspace of A and the of as spherical coordinates on the Bloch sphere. The complementary projector P.InP subspace a pseudoinverse general solution of the Schrödinger equation is then a linear matrix can be defined, and the Q component of h⃗is chosen to combination of the eigenvectors of IðtÞ, i.e., jΨðtÞi ¼ κ κ i þðtÞ ϕ i −ðtÞ ϕ ∈ C make the resulting Hamiltonian realizable (Torrontegui, cþe j þðtÞi þ c−e j −ðtÞi, where c and Martínez-Garaot, and Muga, 2014). Examples for two- [SU(2) algebra] and three-level systems (with the four- 1 κ_ ðtÞ¼ hϕ ðtÞj½iℏ∂ − HðtÞ jϕ ðtÞi. dimensional algebra U3S3) were provided. Levy et al. ℏ t (2018) reformulated this formalism in terms of the density operator and used it to design robust control protocols against Therefore, it is possible to construct a particular solution the influence of different types of noise. In particular, they ψ ϕ −iγðtÞ=2 developed a method to construct a control protocol which is j ðtÞi ¼ j þðtÞie ; ð50Þ robust against dissipation of the population and minimizes the γ 2κ effect of dephasing. where ¼ and
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-13 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, …
1 ϕ γ_ Ω α Ω α nðq; tÞ ¼ θ ð R cos þ I sin Þ: ð51Þ sin im q − q ρ_ 2 2ρ _ ρ − ρ_ ρ ρ−1=2Φ c ¼ exp ℏ ½ q = þðqc qc Þq= n ρ ; Using Eqs. (46), (47), and (51) we can retrieve the physical quantities in terms of the auxiliary functions, ð59Þ
_ Φ σ ΩR ¼ cos α sin θγ_ − sin αθ; ð52Þ where nð Þ is a solution of the stationary Schrödinger equation Ω ¼ sin α sin θγ_ þ cos αθ_; ð53Þ I ℏ2 ∂2 1 − mω2σ2 U σ Φ λ Φ ; 2 ∂σ2 þ 2 0 þ ð Þ n ¼ n n ð60Þ Δ ¼ − cos θγ_ − α_: ð54Þ m with σ ¼ðq − q Þ=ρ. This quadratic invariant has been If the functions α, γ, and θ are chosen with the appropriate c instrumental in designing many of the schemes which boundary conditions, different state manipulations are pos- manipulate trapping potentials for expansions, compressions, sible. For example, θð0Þ¼0 and θðt Þ¼π imply perfect f transport, launching, stopping, or combined processes. It is population inversion at a time t . Note the freedom to f key, for example, to manipulate the motion of trapped ions, interpolate along different paths. see Sec. III.A, and in proposals to implement STA-based interferometry (Dupont-Nivet, Westbrook, and Schwartz, b. Lewis-Leach family 2016; Martínez-Garaot, Rodriguez-Prieto, and Muga, 2018). 2 Designing first the function ρðtÞ which determines the wave Consider a one-dimensional Hamiltonian H ¼ p =2m þ function width, and q ðtÞ (a classical particle trajectory), the Vðq; tÞ with potential (Lewis and Leach, 1982) c force FðtÞ and ωðtÞ can be determined using Eq. (57). The boundary values of the auxiliary functions at the time limits − m ω2 2 Vðq; tÞ¼ FðtÞq þ 2 ðtÞq are fixed to satisfy physical conditions and commutativity of H and the invariant. The nonunique interpolation is usually 1 q − q ðtÞ U c g t : done with polynomials or trigonometric functions. þ ρ 2 ρ þ ð Þ ð55Þ ðtÞ ðtÞ Expansions.—Invariant-based inverse engineering for Lewis-Leach Hamiltonians was first implemented by Chen, These Hamiltonians have a quadratic in momentum invariant Ruschhaupt et al. (2010) to cool down a trapped atom by 1 expanding the trap. Expansions using invariant-based engi- ρ − _ − ρ_ − 2 neering were first implemented with ultracold atoms in a I ¼ 2 ½ ðp mqcÞ m ðq qcÞ m pioneering experiment on STA techniques by Schaff et al. 1 − 2 − 2 q qc q qc (2010) and Schaff, Capuzzi et al. (2011). Such cooling þ mω0 þ U ; ð56Þ 2 ρ ρ protocols have also been envisioned to optimize sympathetic cooling (Choi, Onofrio, and Sundaram, 2011; Onofrio, 2016). ρ ω provided the functions , qc, , and F satisfy the auxiliary For very short processes, tf < 1=ð2ωfÞ, where ωf is the trap equations frequency at the final time, ω2ðtÞ may become negative during some time interval. While this implies a transient repulsive 2 2 ω0 potential, the atoms always remain confined (Chen, ρ̈þ ω ðtÞρ ¼ ; ρ3 Ruschhaupt et al., 2010). A repulsive potential may or may q̈þ ω2ðtÞq ¼ FðtÞ=m; ð57Þ not be difficult to implement depending on the physical c c setting. For example, the analysis by Torrontegui et al. (2018) suggested that it is viable for trapped ions. with ω0 a constant. The first equation is known as the Compared to the simplicity of invariant-based engineering, Ermakov equation (Ermakov, 1880), while the second is the CD approach for expansions or compressions provides, for the Newton equation of motion for a forced harmonic an H0 characterized by some predetermined time-dependent oscillator. They can be found by inserting the quadratic-in- frequency ω t , a nonlocal, cumbersome counterdiabatic term p invariant, Eq. (56), into Eq. (37). The properties of such ð Þ H − pq qp ω_= 4ω (Muga et al., 2010). However, a invariants have also been formulated in terms of Feynman CD ¼ ð þ Þ ð Þ unitary transformation produces a new shortcut with local propagators (Dhara and Lawande, 1984). potential and modified frequency (Ibáñez et al., 2012) For this family of Hamiltonians we can explicitly calculate the Lewis-Riesenfeld phase, 3ω_ 2 ω̈ 1=2 ω0 ω2 − Z ¼ 4ω2 þ 2ω ; ð61Þ 1 t λ m q_ ρ − q ρ_ 2 − ω2q2=ρ2 α − 0 n ½ð c c Þ 0 c nðtÞ¼ ℏ dt ρ2 þ 2ρ2 þ g ; 0 see further connections among the two Hamiltonians in ð58Þ del Campo (2013) and Mishima and Izumida (2017). So far we have considered only one-dimensional motion. and the eigenvectors in coordinate representation, Formally, the three coordinates in an ideal harmonic trap are
Rev. Mod. Phys., Vol. 91, No. 4, October–December 2019 045001-14 D. Gu´ery-Odelin et al.: Shortcuts to adiabaticity: Concepts, methods, … uncoupled so expansion or transport processes can be treated avoided this difficulty for an infinite-range Ising model in a independently. However, in many cold atom experiments transverse field by constructing an invariant using a mean- changing the intensity of a laser beam simultaneously affects field Ansatz. However, another approach is to exploit scaling the longitudinal and transversal frequencies (Torrontegui, laws. If the Hamiltonian fulfills certain scaling laws one can Chen, Modugno, Ruschhaput et al., 2012). determine the invariant in a similar manner as in Sec. II.C.2.b. Transport.—Designing particle transport has been another Specifically the Hamiltonian major application of the invariant (56) (Torrontegui et al., 2011; Ness et al., 2018; Tobalina, Alonso, and Muga, 2018); XN p2 X see also Secs. II.G and III.A. Note that U in Eq. (55) is HðtÞ¼ i þ U½q ; λðtÞ þ ϵðtÞ Vðq − q Þ 2m i i j arbitrary. Setting ω ¼ ω0 ¼ 0, and ρ ¼ 1, F ¼ mq̈ plays the i¼1 i