Formulation and Numerical Solution of Quantum Control Problems

Total Page:16

File Type:pdf, Size:1020Kb

Formulation and Numerical Solution of Quantum Control Problems Index adjoint equation, 86, 113 Coulomb potential, 44, 49, 347 Gâteaux differentiability, 79 Armijo condition, 223, 231 Crank–Nicolson scheme, 182 Galerkin projection, 282 Arzelà–Ascoli theorem, 140 generalized Bloch vectors, 67 azimuthal quantum number, 24 de Broglie’s wave hypothesis, 8 generalized Jacobian, 245 density functional theory, 43 generalized Laguerre BFGS method, 228 density functional theory polynomials, 348 black-box methods, 222 (DFT), 43 global phase, 44 Bloch sphere, 65 density matrix, 27, 29, 31 group velocity, 9 Bochner space, 359 density operator as a gyromagnetic ratio, 18 Bohr radius, 348 rank-one-projector, 29 Born–Oppenheimer DFT’s function sets, 44 Hamilton operator of a approximation, 42 differentiability spherical symmetric Bose–Einstein condensate, 39 real-Fréchet, 133 problem, 345 bosons, 17, 39 real-Gâteaux, 133 Hamiltonian dynamics, 62 bra and ket vectors, 28 dipole approximation, 331, 343 Hartree potential, 50, 128 Dirac formalism, 28 Hermitian operator, 11 canonical momentum, 16 discretize-before-optimize, 236, Hessian, 245, 256 Carathéodory’s theorem, 76 254 heteronuclear spins, 33 cascadic black-box schemes, 230 Duhamel form, 140 hierarchy of objectives, 281 charged particle, 16 Hohenberg–Kohn theorem, 47 coercivity of the Hessian, 98 Ehrenfest’s theorem, 12, 15 homodyne measurements, 64 commutation relation, 15 eigenvalue problem, 14 homotopy continuation commutator, 12, 15 electric scalar potential, 16 scheme, 275 conical eigenvalue intersection, electron density, 43 169, 179 ensemble of quantum states, 27 implicit midpoint scheme, 182, connectivity graph of the exact and approximate 235 control Hamiltonian, controllability, 163 implicit-function theorem, 77 168 exact-controllability problem, infinite-barrier well potential, continuation for initializing 291 155 optimization schemes, exchange potential, 129 infinite-well potential, 14 274 exponential midpoint scheme, intrinsic angular momentum, continuation technique, 292 183 17 continuity equation, 51 control constraints, 231 fermions, 17 KKT system, 111 control Hamiltonian, 164, 354 fidelity, 318 Kohn–Sham (KS) approach, 49 control-to-state map, 77 finite-dimensional Krotov method, 232 controlled Fokker–Planck Schrödinger–Pauli Krylov–Newton methods, 242 equation, 66 equation, 27 KS potential, 50, 128 convolution theorem, 132 Fokker–Planck equation, 65 correlation potential, 129 Fréchet differentiability, 79, Lagrange framework, 88 Coulomb gauge condition, 16 104, 133, 253 Landau levels, 16 389 390 Index Laplacian in spherical Pauli matrices, 20 spin-down, 17 coordinates, 345 Pauli principle, 38, 50 spin quantum number, 18, 25 Larmor frequency, 19, 33 Planck’s constant, 8 spin-spin interaction graph, 172 Legendre polynomials and posterior state vector, 64 spin-up, 17 Rodrigues formula, 347 precession, 19 spinor, 19 Lindblad operator, 63 principal quantum number, 24 stability in the sense of Lindblad–Kossakowski probability density function, 65 Hadamard, 76 equation, 31, 63 projection method, 231 steepest descent method, 223 linearized constraint equation, pure quantum state, 28 stimulated Raman adiabatic 86 passage (STIRAP), 73 Liouville–von Neumann quadratic and superlinear stochastic Belavkin model, 65 equation, 31, 59 convergence, 244 stochastic diffusion process, 64 Lorenz force, 16 quadratic model of the stochastic Schrödinger equation objective, 281 (SE), 64 magnetic moment of a particle, quantum filtering equation, 161 subdifferential, 119, 245 18 quantum Liouville equation, 57 subgradients, 262 magnetic quantum number, 25 quantum nondemolition sufficient optimality magnetic resonance imaging, 19 principle, 161 conditions, 87 mixed quantum state, 28 quantum numbers, 24 symmetric or antisymmetric momentum-stress tensor, 52 quasi-Newton methods, 228 wavefunction, 38 monotonic scheme, 232 symplectic matrix, 183 multigrid optimization radial quantum number, 24 (MGOPT) method, 277 radiofrequency pulse or time reversibility, 182 multiparticle quantum system, control, 19, 37 time-dependent density 37 reachable set, 163, 293 functional theory multiple controls, 164 real-linear, 133 (TDDFT), 50 reduced cost functional, 84 time-dependent Kohn–Sham necessary optimality reduced generalized Jacobian, (TDKS) model, 56 conditions, 87, 119 265 two-level system, 167 nonlinear complementarity reduced gradient, 84 problem (NCP) reduced Hessian, 248 uncertainty principle, 16 functions, 259 reduced Planck’s constant, 8 unitary orbit, 171 nonlinear conjugate gradient rotating-frame transformation, unitary propagator, 182 (NCG) method, 224 33 nonlinear multigrid methods, van Leeuwen, theorem of, 55 276 Schrödinger equation (SE), 9, vector potential, 20 nuclear magnetic resonance 12, 14 (NMR), 19 Schur reduction, 257 wave packet, 9 nuclear magneton, 18 semiclassical limit, 61 weak sequential continuity, 80 open quantum system, 65 Sherman–Morrison–Woodbury Wigner function, 57, 58 operator controllability, 164 formula, 228 Wirtinger calculus, 253 optimality condition, 86 shooting-type method, 299 Wolfe condition, 224, 231 optimality system, 93, 109 Slater determinant, 38, 46 optimize-before-discretize, 236, Slater permanent, 38 Young’s inequality, 132 253 spherical Bessel functions, 349 orbital angular momentum, 17 spherical harmonics, 23, 346, Zeeman effect, 18, 25 350 Zowe–Kurcyusz constraint parity operator, 38 spin, 17 qualification, 90.
Recommended publications
  • Ehrenfest Theorem
    Ehrenfest theorem Consider two cases, A=r and A=p: d 1 1 p2 1 < r > = < [r,H] > = < [r, + V(r)] > = < [r,p2 ] > dt ih ih 2m 2imh 2 [r,p ] = [r.p]p + p[r,p] = 2ihp d 1 d < p > ∴ < r > = < 2ihp > ⇒ < r > = dt 2imh dt m d 1 1 p2 1 < p > = < [p,H] > = < [p, + V(r)] > = < [p,V(r)] > dt ih ih 2m ih [p,V(r)] ψ = − ih∇V(r)ψ − ()− ihV(r)∇ ψ = −ihψ∇V(r) ⇒ [p,V(r)] = −ih∇V(r) d 1 d ∴ < p > = < −ih∇V(r) > ⇒ < p > = < −∇V(r) > dt ih dt d d Compare this with classical result: m vr = pv and pv = − ∇V dt dt We see the correspondence between expectation of an operator and classical variable. Poisson brackets and commutators Poisson brackets in classical mechanics: ⎛ ∂A ∂B ∂A ∂B ⎞ {A, B} = ⎜ − ⎟ ∑ ⎜ ⎟ j ⎝ ∂q j ∂p j ∂p j ∂q j ⎠ dA ⎛ ∂A ∂A ⎞ ∂A ⎛ ∂A ∂H ∂A ∂H ⎞ ∂A = ⎜ q + p ⎟ + = ⎜ − ⎟ + ∑ ⎜ & j & j ⎟ ∑ ⎜ ⎟ dt j ⎝ ∂q j ∂p j ⎠ ∂t j ⎝ ∂q j ∂q j ∂p j ∂p j ⎠ ∂t dA ∂A ∴ = {A,H}+ Classical equation of motion dt ∂ t Compare this with quantum mechanical result: d 1 ∂A < A > = < [A,H] > + dt ih ∂ t We see the correspondence between classical mechanics and quantum mechanics: 1 [A,H] → {A, H} classical ih Quantum mechanics and classical mechanics 1. Quantum mechanics is more general and cover classical mechanics as a limiting case. 2. Quantum mechanics can be approximated as classical mechanics when the action is much larger than h.
    [Show full text]
  • Lecture 8 Symmetries, Conserved Quantities, and the Labeling of States Angular Momentum
    Lecture 8 Symmetries, conserved quantities, and the labeling of states Angular Momentum Today’s Program: 1. Symmetries and conserved quantities – labeling of states 2. Ehrenfest Theorem – the greatest theorem of all times (in Prof. Anikeeva’s opinion) 3. Angular momentum in QM ˆ2 ˆ 4. Finding the eigenfunctions of L and Lz - the spherical harmonics. Questions you will by able to answer by the end of today’s lecture 1. How are constants of motion used to label states? 2. How to use Ehrenfest theorem to determine whether the physical quantity is a constant of motion (i.e. does not change in time)? 3. What is the connection between symmetries and constant of motion? 4. What are the properties of “conservative systems”? 5. What is the dispersion relation for a free particle, why are E and k used as labels? 6. How to derive the orbital angular momentum observable? 7. How to check if a given vector observable is in fact an angular momentum? References: 1. Introductory Quantum and Statistical Mechanics, Hagelstein, Senturia, Orlando. 2. Principles of Quantum Mechanics, Shankar. 3. http://mathworld.wolfram.com/SphericalHarmonic.html 1 Symmetries, conserved quantities and constants of motion – how do we identify and label states (good quantum numbers) The connection between symmetries and conserved quantities: In the previous section we showed that the Hamiltonian function plays a major role in our understanding of quantum mechanics using it we could find both the eigenfunctions of the Hamiltonian and the time evolution of the system. What do we mean by when we say an object is symmetric?�What we mean is that if we take the object perform a particular operation on it and then compare the result to the initial situation they are indistinguishable.�When one speaks of a symmetry it is critical to state symmetric with respect to which operation.
    [Show full text]
  • Schrödinger and Heisenberg Representations
    p. 33 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mathematical formulation of the dynamics of a quantum system is not unique. So far we have described the dynamics by propagating the wavefunction, which encodes probability densities. This is known as the Schrödinger representation of quantum mechanics. Ultimately, since we can’t measure a wavefunction, we are interested in observables (probability amplitudes associated with Hermetian operators). Looking at a time-evolving expectation value suggests an alternate interpretation of the quantum observable: Atˆ () = ψ ()t Aˆ ψ ()t = ψ ()0 U † AˆUψ ()0 = ( ψ ()0 UA† ) (U ψ ()0 ) (4.1) = ψ ()0 (UA † U) ψ ()0 The last two expressions here suggest alternate transformation that can describe the dynamics. These have different physical interpretations: 1) Transform the eigenvectors: ψ (t) →U ψ . Leave operators unchanged. 2) Transform the operators: Atˆ ( ) →U† AˆU. Leave eigenvectors unchanged. (1) Schrödinger Picture: Everything we have done so far. Operators are stationary. Eigenvectors evolve under Ut( ,t0 ) . (2) Heisenberg Picture: Use unitary property of U to transform operators so they evolve in time. The wavefunction is stationary. This is a physically appealing picture, because particles move – there is a time-dependence to position and momentum. Let’s look at time-evolution in these two pictures: Schrödinger Picture We have talked about the time-development of ψ , which is governed by ∂ i ψ = H ψ (4.2) ∂t p. 34 in differential form, or alternatively ψ (t) = U( t, t0 ) ψ (t0 ) in an
    [Show full text]
  • Ph125: Quantum Mechanics
    Section 10 Classical Limit Page 511 Lecture 28: Classical Limit Date Given: 2008/12/05 Date Revised: 2008/12/05 Page 512 Ehrenfest’s Theorem Ehrenfest’s Theorem How do expectation values evolve in time? We expect that, as quantum effects becomes small, the fractional uncertainty in a physical observable Ω, given by q h(∆Ω)2i/hΩi, for a state |ψ i, becomes small and we need only consider the expectation value hΩi, not the full state |ψ i. So, the evolution of expectation values should approach classical equations of motion as ~ → 0. To check this, we must first calculate how expectation values time-evolve, which is easy to do: d „ d « „ d « „ dΩ « hΩi = hψ | Ω |ψ i + hψ |Ω |ψ i + hψ | |ψ i dt dt dt dt i „ dΩ « = − [− (hψ |H)Ω |ψ i + hψ |Ω(H|ψ i)] + hψ | |ψ i ~ dt i „ dΩ « = − hψ |[Ω, H]|ψ i + hψ | |ψ i ~ dt i fi dΩ fl = − h[Ω, H]i + (10.1) ~ dt Section 10.1 Classical Limit: Ehrenfest’s Theorem Page 513 Ehrenfest’s Theorem (cont.) If the observable Ω has no explicit time-dependence, this reduces to d i hΩi = − h[Ω, H]i (10.2) dt ~ Equations 10.1 and 10.2 are known as the Ehrenfest Theorem. Section 10.1 Classical Limit: Ehrenfest’s Theorem Page 514 Ehrenfest’s Theorem (cont.) Applications of the Ehrenfest Theorem To make use of it, let’s consider some examples. For Ω = P, we have d i hPi = − h[P, H]i dt ~ We know P commutes with P2/2 m, so the only interesting term will be [P, V (X )].
    [Show full text]
  • Ehrenfest's Theorem and Nonclassical States of Light
    GENERAL ARTICLE Ehrenfest’s Theorem and Nonclassical States of Light 1. Ehrenfest’s Theorem in Quantum Mechanics Lijo T George, C Sudheesh, S Lakshmibala and V Balakrishnan When it was ¯rst enunciated, Ehrenfest's The- orem provided a necessary and important link between classical mechanics and quantum me- chanics. Today, the content of the theorem is 1 2 understood to be a natural and immediate con- sequence of the equation of motion for opera- tors when quantum mechanics is formulated in 3 4 the Heisenberg picture. Nevertheless, the theo- rem leads to useful approximations when systems 1 Lijo T George got his with Hamiltonians of higher than quadratic or- Masters degree in physics from IIT Madras, and is der in the dynamical variables are considered. presently a doctoral student In this two-part article, we use it to provide a at RRI, Bangalore. His convenient illustration of the di®erences between research interests include so-called classical and nonclassical states of radi- quantum mechanics and ation. cosmology. 2C Sudheesh is with the 1. Paul Ehrenfest and the Quantum{Classical Department of Physics, Connection (IIST), Thiruvanan- thapuram. His research The Austrian-Dutch physicist Paul Ehrenfest (1880{1933) interests are nonlinear stands out among those who worked in the grey area dynamics, quantum optics of semiclassical physics straddling quantum mechanics and quantum information. (QM) and classical mechanics (CM). His work provided 3S Lakshmibala is with the a tangible correspondence between the quantum and Department of Physics, IIT Madras. Her research classical worlds, through the Ehrenfest theorem. This interests are dynamical theorem brought comfort to the practitioners of physics systems, quantum dynamics at a time in history when glimpses of the quantum world and quantum optics.
    [Show full text]
  • Classical from Quantum Arxiv:1803.04700V2 [Quant-Ph] 20 Apr 2018
    Classical from Quantum Timothy J. Hollowood Department of Physics, Swansea University, Swansea, SA2 8PP, United Kingdom E-mail: [email protected] Abstract: We consider the quantum-to-classical transition for macroscopic systems coupled to their environments. By applying Born's rule, we are led to a particular set of quantum trajectories, or an unravelling, that describes the state of the system from the frame of reference of the subsystem. The unravelling involves a branch dependent Schmidt decomposition of the total state vector. The state in the subsystem frame, the conditioned state, is described by a Poisson process that involves a non-linear deterministic effective Schr¨odingerequation interspersed with quantum jumps into orthogonal states. We then consider a system whose classical analogue is a generic chaotic system. Although the state spreads out exponentially over phase space, the state in the frame of the subsystem localizes onto a narrow wave packet that follows the classical trajectory due to Ehrenfest's Theorem. Quantum jumps occur with a rate that is the order of the effective Lyapunov exponent of the classical chaotic system and imply that the wave packet undergoes random kicks described by the classical Langevin equation of Brownian motion. The implication of the analysis is that this theory can explain in detail how classical mechanics arises from quantum mechanics by using only unitary evolution and Born's rule applied to a subsystem. arXiv:1803.04700v2 [quant-ph] 20 Apr 2018 1 Introduction Quantum mechanics is a remarkable fusion of a linear, deterministic theory where evolution involves unitary rotations of a vector in a Hilbert space, with stochastic evolution, in the form of Born's rule.
    [Show full text]
  • Equations of Motion for a Relativistic Wave Packet
    PRAMANA c Indian Academy of Sciences Vol. 78, No. 5 — journal of May 2012 physics pp. 679–685 Equations of motion for a relativistic wave packet LKOCIS1,2 1Julius Kruttschnitt Mineral Research Centre, The University of Queensland, Isles Road, Indooroopilly, Queensland 4068, Australia 2Present address: Peranga Court Unit 4, 43 Fifth Avenue, Sandgate, Queensland 4017, Australia E-mail: [email protected] MS received 15 April 2011; revised 23 November 2011; accepted 23 December 2011 Abstract. The time derivative of the position of a relativistic wave packet is evaluated. It is found that it is equal to the mean value of the momentum of the wave packet divided by the mass of the particle. The equation derived represents a relativistic version of the second Ehrenfest theorem. Keywords. Klein–Gordon equation; classical limit; Ehrenfest theorems. PACS Nos 03.65.−w; 03.65.Pm 1. Introduction In this article we consider a small relativistic wave packet that moves in external static fields and evaluate the time derivative of its position. The classical limit of the Klein–Gordon equation can be achieved by assuming that the Planck’s constant converges to zero [1] or by assuming that the coordinate and momentum dispersions of a relativistic wave packet are small [2]. At this point it is ne- cessary to give a brief review of what was achieved in ref. [2]. In ref. [2] substitution i ψ = R θ , exp (1.1) where R and θ are real functions, was applied in the Klein–Gordon formula for the particle probability density, namely i ∂ψ ∂ψ∗ 2 eφ ρ = ψ∗ − ψ − ψ∗ψ.
    [Show full text]
  • Validity of Ehrenfest's Theorem for Generalized Fields of Dyons
    Turkish Journal of Physics Turk J Phys (2016) 40: 244 { 255 http://journals.tubitak.gov.tr/physics/ ⃝c TUB¨ ITAK_ Research Article doi:10.3906/fiz-1512-11 Validity of Ehrenfest's theorem for generalized fields of dyons Gaurav KARNATAK∗, Praveen Singh BISHT, Om Prakash Singh NEGI Department of Physics, Kumaun University, S. S. J. Campus, Almora, Uttarakhand, India Received: 28.12.2015 • Accepted/Published Online: 18.04.2016 • Final Version: 01.12.2016 Abstract: The validity of Ehrenfest's theorem with its classical correspondence has been justified for the manifestly covariant equations of dyons. We have also developed accordingly the Lagrangian formulation for electromagnetic fields in a minimum coupled source giving rise to conserved current of dyons. Applying the Gupta subsidiary condition we have extended the validity of the Ehrenfest's theorem for U (1) × U (1) Abelian gauge theory of dyons. It is shown that the expectation value of the quantum equation of motion reproduces the classical equation of motion, which is the generalized form of Ehrenfest's theorem in quantum field theory. Key words: Dyons, Ehrenfest's theorem, Lagrangian, quantum field theory, electromagnetic fields 1. Introduction In classical mechanics the dynamical state of each particle is defined by its position and momentum, at a given instant. In quantum mechanics the dynamical state of the system is represented by its wave function [1]. The expectation values of displacement and momentum obey time evolution equations, which are analogous to those of classical mechanics. This result is called Ehrenfest's theorem [2]. The value of the quantum equation of motion taken with the states in the physical subspace reproduces the classical equations of motion.
    [Show full text]
  • Technische Universität München Department of Mathematics
    Technische Universität München Department of Mathematics Evolution of Angular Momentum Expectation in Quantum Mechanics Bachelor’s Thesis by Thomas Wiatowski supervisor : Prof. Dr. Gero Friesecke academic advisor : Dipl.-Math. Mario Koppen submission date : August 9, 2010 Ich erkläre hiermit, daß ich die Bachelorarbeit selbständig und nur mit den angegebe- nen Hilfsmitteln angefertigt habe. ............................................................. (München, der 09. August 2010), (Unterschrift) ii Acknowledgements I would like to thank all people who have helped and inspired me during my bachelor study. Especially I want to thank my academic advisor, Dipl.-Math. Mario Koppen, for his guidance during my bachelor thesis in the sixth semester at the Technische Universität München. He was always accessible and willing to help me with the research on my thesis. As a result, research and writing became smooth and rewarding for me. I was delighted to interact with Prof. Dr. Gero Friesecke by attending his class “Funk- tionalanalysis”, which gave me the mathematical background for this thesis and to have him as my supervisor. I am indebted to many of my colleagues who supported me, especially Friedrich Prade. His perpetual energy and enthusiasm in physics has motivated and helped me throughout my thesis. My deepest gratitude goes to my parents, Anna and Norbert, for their enduring love and support throughout my life, supporting me and my sister and sparing no effort to provide the best possible environment for us to grow up and attend school. Thank you. iii Abstract A mathematically rigorous derivation of the evolution of angular momentum expec- tation value is given, under assumptions which include Hamiltonians with potential wells.
    [Show full text]
  • On the Emergent Aspects of Quantum Mechanics in Relation to the Thermodynamics of Irreversible Processes
    1 ON THE EMERGENT ASPECTS OF QUANTUM MECHANICS IN RELATION TO THE THERMODYNAMICS OF IRREVERSIBLE PROCESSES AND EMERGENT GRAVITY 2 3 Candidato a doctor: Acosta Iglesias, Dagoberto. T´ıtulo: Aspectos emergentes de la mecanica´ cuantica:´ relacion´ con la termodinamica´ de los procesos irreversibles y con la gravedad emergente. Director: Fernandez´ De Cordoba´ Castella,´ Pedro Jose.´ Director: Ferrando Cogollos, Albert. Director: Isidro San Juan, Jose´ Mar´ıa. Entidad Academica´ Responsable: Departamento de Matematica´ Aplicada. Programa: Programa de Doctorado en Matematicas.´ 4 Abstract This PhD thesis elaborates on a proposal made by the Dutch theoretical physicist G. ’t Hooft (1999 Nobel prize in physics), to the effect that quantum mechanics is the emergent theory of some underlying, deterministic theory. According to this proposal, information–loss effects in the underlying deterministic theory lead to the arrangement of states of the latter into equivalence classes, that one identifies as quantum states of the emergent quantum mechanics. In brief, quantisation is dissipation, according to ’t Hooft. In our thesis we present two mechanisms whereby quantum mechanics is explicitly seen to emerge, thus explicitly realising ’t Hooft’s proposal. The first mechanism makes use of Verlinde’s approach to classical mechanics and general relativity via holographic screens. This technique, first presented in 2010 in order to understand the emergent nature of spacetime and gravity, is applied in our thesis to the case of quantum mechanics. The second mechanism presented to support ’t Hooft’s statement is based on a dictionary, also developed by the authors, between semiclassical quantum mechanics, on the one hand, and the classical theory of irreversible thermodynamics, on the other.
    [Show full text]
  • Arxiv:1609.07380V1 [Quant-Ph] 23 Sep 2016 Where Nntl Ag Oeta Nry Uha H Ia Et Ucinpo Function Delta Confi [1]
    Direct manifestation of Ehrenfest’s theorem in the infinite square well model Chyi-Lung Lin∗ Department of Physics, Soochow University, Taipei,Taiwan, R.O.C. Ehrenfest’s theorem in the infinite square well is up to now only manifested indirectly. The manifestation of this theorem is first done in the finite square well, and then consider the infinite square well as the limit of the finite well. For a direct manifestation, we need a more precise formula to describe the degree of infiniteness of the divergent potential energy. We show that the potential energy term term, V (x)Ψn(x), which is the product of the potential energy and the energy eigenfunction, is a well defined function which can be expressed in terms of Dirac delta functions. This means that the infinity in this model is not that vague but has obtained a specification. This results that expectation values can be calculated precisely and Ehrenfest’s thereom can be confirmed directly. PACS numbers: 03.65.-w, Quantum mechanics; 03.65.Ge, solution of wave equations: bound states Keywords: Infinite square well; Ehrenfest theorem; Time independent Schr¨odinger equation; Step function; Dirac delta function; Boundary condition; I. INTRODUCTION The infinite square well is a model using infinitely high potential barrier to confine particles inside a well. This is a basic model in quantum mechanics, and is a standard model for confining particles. We also have other models using infinitely large potential energy, such as the Dirac delta function potential. Belloni and Robinett have given a review on these models [1]. Curiously, for infinite square well, we still do not have a direct manifestation of Ehrenfest’s theorem in this funda- mental model.
    [Show full text]
  • Exact Time Evolution and Master Equations for the Damped Harmonic Oscillator Abstract
    CORE Metadata, citation and similar papers at core.ac.uk Provided by CERN DocumentExact Server Time Evolution and Master Equations for the Damped Harmonic Oscillator Robert Karrlein and Hermann Grabert Fakult¨at f¨ur Physik der Albert–Ludwigs–Universit¨at, Hermann-Herder-Str. 3, D-79104 Freiburg, Germany (August 26, 1996) Abstract Using the exact path integral solution for the damped harmonic oscillator it is shown that in general there does not exist an exact dissipative Liouville operator describing the dynamics of the oscillator for arbitrary initial bath preparations. Exact non-stationary Liouville operators can be found only for particular preparations. Three physically meaningful examples are exam- ined. An exact new master equation is derived for thermal initial conditions. Second, the Liouville operator governing the time-evolution of equilibrium correlations is obtained. Third, factorizing initial conditions are studied. Ad- ditionally, one can show that there are approximate Liouville operators inde- pendent of the initial preparation describing the long time dynamics under appropriate conditions. The general form of these approximate master equa- tions is derived and the coefficients are determined for special cases of the bath spectral density including the Ohmic, Drude and weak coupling cases. The connection with earlier work is discussed. Typeset using REVTEX 1 I. INTRODUCTION Recently the problem of the reduced dynamics of a quantum system in contact with a reservoir has gained renewed interest [1–3]. While the foundations of quantum dissipative processes were layed already in the sixties [4], this early work was mainly concerned with weakly damped systems and has relied on the Born and Markov approximations.
    [Show full text]