Lyapunov-based control of quantum systems
|
|
- Kathlyn Lawrence
- 5 years ago
- Views:
Transcription
1 Lyapunov-based control of quantum systems Symeon Grivopoulos Bassam Bamieh Department of Mechanical and Environmental Engineering University of California, Santa Barbara, CA Abstract We propose feedback control laws for quantum systems, based on Lyapunov design. These feedback laws locally asymptotically stabilize any desired eigenstate of the system Hamiltonian. The target eigenstate and the size of its region of attraction can be tailored by the choice of design parameters. These feedback control laws are used to design open loop controls for state transfer problems. We illustrate the method with a transfer problem between bound states of a molecular Morse oscillator. I. INTRODUCTION The model of a controlled quantum system (Schroedinger equation) is i ψ = ( H + α H α u α (t) ) ψ, () where ψ(t) is the state vector of the quantum system that belongs to an appropriate (finite or infinite-dimensional) complex Hilbert space, and H, H α are self-adjoint operators in this space, referred to as the system and control Hamiltonians respectively. H describes the natural dynamics of the system, and H α the coupling of the system to external classical control fields u α (t). The dynamics generated by () is unitary (norm preserving). By convention, we normalize ψ =. Perhaps the most significant question about controlled quantum systems concerns their controllability properties, such as controllability conditions on () (see [], [2]) and the design of controls for steering (see e.g. [3] for constructive control, [4], [5] for optimal control and [6] for time-optimal control). This work also concerns itself with a steering problem inspired by Quantum Chemistry: Design of controls that drive the system to an eigenstate of H. These states are very important in their own right: Each of them corresponds to a definite measurable energy of the system and being on one or the other energy eigenstate determines the chemical behavior of the system. The class of controls we propose are in state feedback form and are based on Lyapunov design. We prove that they asymptotically drive any initial state to an eigenstate of H. Although the controls are in feedback form, the collapse of the state associated with measurements in Quantum Mechanics would prevent us from using them as such in an experimental setting. Rather, we only use our method to generate open loop controls for specific transfer problems. For related approaches, see [7], [8]. In the following, we consider only finite-dimensional quantum systems, use units where = and consider () with only one control: ı ψ = ( H + H u(t) ) ψ. (2) II. CONTROL DESIGN We use V (ψ) = ψ P ψ, P = P, as a Lyapunov function. Our first task is to analyze the structure of its critical points on the sphere ψ =. Lemma : The set of critical points of V (ψ) = ψ P ψ on the sphere ψ ψ = is given by the normalized eigenvectors of P. The eigenvectors with the largest eigenvalue are the maxima, the eigenvectors with the smallest eigenvalue are the minima and all others are saddle points. Proof: To find the critical points of V (ψ) under the constraint ψ ψ =, consider the augmented function Ṽ (ψ, λ) = ψ P ψ λ (ψ ψ ), where λ is a Lagrange multiplier. Varying Ṽ with respect to ψ we have, ( ) Ṽ = = P ψ λψ, ψ so, ψ is an eigenvector of P and λ the corresponding eigenvalue. Varying Ṽ with respect to λ enforces the normalization of ψ. Denote the normalized eigenvectors of P by q i, i =,..., N, and the corresponding eigenvalues by p i. Then, P = N p iq i q i. To determine the structure of V around one of its critical points, q N for example, consider a finite variation ψ such that q N + ψ 2 =. Express ψ in the basis of the eigenvectors of P, ψ = N ψ i q i. The normalization condition q N + ψ 2 = implies that ( ψ N + ψn + ψ N ψn) + ψ i ψi =,
2 that is, not all components of the variation can be chosen arbitrarily. Consider the difference V (q N + ψ) V (q N ) = p N ( ψ N + ψn + ψ N ψn) + p i ψ i ψi = (p i p N ) ψ i ψi. Taking the ψ i, i =,..., N and arg( ψ N ), as the independent parameters of the variation, we see that the structure of V around q N depends on the ordering of the eigenvalues: q N is a local maximum iff p N is the largest eigenvalue, a local minimum iff p N is the smallest eigenvalue and a saddle point otherwise. We calculate the rate of change of V along a trajectory of (2) : dv dt = ıψ P (H + H u)ψ + ıψ (H + H u)p ψ = ıψ [H + H u, P ]ψ = ıψ [H, P ]ψ + ıψ [H, P ]ψ u, where [A, B] := A B B A, denotes the commutator of two matrices. Since tr([a, B]) = tr(ab) tr(ba) =, a commutator can never be sign definite. For this reason, we choose P to commute with H : [H, P ] =. This implies that the two matrices must have the same eigenvectors and hence, by Lemma, the critical points of V are the eigenvectors of H. Now, dv dt = ıψ [H, P ]ψ u. To make V, we choose the following (feedback) control law: u = f(ıψ [H, P ]ψ), (3) where f is a continuous function with f() = and xf(x) >, x. With this choice of u, (2) becomes ı ψ = ( H f(ıψ [H, P ]ψ) H ) ψ. (4) We show now that the set of eigenvectors of H is the largest invariant set of (4). For ψ to be in the invariant set of (4), V = u 2 =, which means that must be satisfied along with ı ψ (t) [H, P ] ψ(t) = (5) ı ψ(t) = ( H f( ı ψ (t)[h, P ] ψ(t) ) H ) ψ(t), with initial condition ψ. The solution is given by ψ(t) = e ıht ψ N = c i e ıeit q i, where ψ = N c i q i and the E i s are the eigenvalues of H. Substituting this in (5), we obtain N ( ı) (p i p j ) c i c j e ıωjit (q jh q i ) =, (6) i,j where ω ji := E j E i is the transition frequency between levels i and j. We make the following assumptions: () ω ji ω lk, (i, j) (k, l). (2) q j H q i, i j. (3) p i p j, i j. Assumption () guarantees that the exponential functions in (6) are linearly independent, while (2) and (3) ensure that the resulting equations, c i c j =, hold i j. The only non-trivial solutions of these have only one c i non-zero. This means that ψ is one of the q i s. Hence, the largest invariant set of (4) is the set of eigenvectors of H. LaSalle s invariance principle [9] guarantees its asymptotic stability. We summarize our result in Theorem : Consider system (2) with the feedback control law (3), where H and P commute. Given the assumptions () (3), this control law asymptotically drives any trajectory of (2) to an eigenspace of H. This analysis can be generalized: (i) One may also consider a P with multiple eigenvalues. Suppose that P has only M < N distinct eigenvalues. Group the eigenvectors of P in M corresponding eigenspaces. One can see that the conclusions of the theorem continue to hold if we substitute eigenspaces for eigenvectors. One still needs assumptions ()-(3), but the indices need only belong to levels in different eigenspaces of P. (ii) For systems with more controls, instead of (2) one would need the weaker condition (2 ) q j H αq i, i j, for some H α. This theorem only guarantees that all trajectories of (2) will approach asymptotically one of the eigenstates of H. For it to have any practical value, one should be able to specify the target energy eigenstate. One should also know the region of attraction of this eigenstate, that is the set of initial conditions that are asymptotically driven to it. Suppose that the desired eigenstate is q N. Make the corresponding eigenvalue, p N, strictly smaller than the others and then, by the Lemma, q N will be the unique minimum of V. Then, every state in the set A N = {ψ S 2 V (ψ) < min(p,..., p )} will asymptotically end up to q N, because energy decreases and the only critical point inside A N is q N. Incidentally, the region of attraction of q N may be bigger than A N but we can only prove convergence from A N. This is an open set whose volume depends on the sizes of the p i s. If one lumps p,..., p close together and makes p N much smaller, one can make it be
3 almost the whole sphere. As a limiting case, one can take p =... = p > p N. Then, V becomes a true Lyapunov function and A N (which will now contain the whole sphere except for the equator {ψ S 2 ψ q N = } = {ψ S 2 ψ span{q,..., q }}) is the whole region of attraction of q N, see also [8]. Nevertheless, the case of more than two distinct p i s is of interest in applications. This occurs, for example, when one wants to realize a transition between two levels and at the same time avoid populating a third level. Finally, 3 remarks: (i)numerical experimentation with these feedback laws lead us to suspect that the minimum of V is almost-globally asymptotically attractive, that is, with the exception of a set of initial states of measure zero, every trajectory of (2) will asymptotically approach the minimum of V. This is only an intuition which we cannot prove at the moment. (ii)assumptions () and (2) concern only the system and imply controllability of (2). As a matter of fact, they are somewhat stronger than controllability as given by the Lie Algebra Rank Condition [], [2]. Nevertheless, they are met by many realistic systems or, in some cases, can be satisfied by turning on a small constant control u and considering H = H + H u. Satisfying assumption (3) and, more generally, the ordering of the p i s, is part of the design. (iii) The controls generated by (3) are bounded, u < max f(x). x 2 H P III. APPLICATION TO A MOLECULAR SYSTEM We apply our method to transfer problems in a model that describes the vibrational dynamics of a two atom molecule, the Morse oscillator. The relevant Hamiltonian is with the Morse potential d 2 H = 2 2m dr 2 + V (r), V (r) = V {exp[ α(r r )] } 2 V, This is an infinite dimensional quantum system. r is the intramolecular distance, m is the reduced mass of the system, r the equilibrium distance of the molecule, α determines the width of the molecular well and V its depth. For these parameters, we take the values that correspond to the OH bond. The theory of the Morse oscillator predicts that the OH system supports 22 bound states as well as continuum states. The control is the time-varying electric field E(t) of a laser that enters through the semi-classical interaction Hamiltonian (in the electric dipole approximation) H int (t) = µ(r)e(t) and µ(r) is the molecular dipole function µ(r) = µ r exp( r/r ) (µ and r are parameters of the interaction). For the application of our technique, we truncate the (infinite dimensional) state space down to the subspace of bound states. Ignoring the continuous part of the spectrum, introduces a small error to our numerical results, but these are mostly meant to be an illustration. We consider a system initially in its ground state and construct a control that excites it to the th level and then to its last (22nd) level. Figures () and (2) display the control and the populations of some bound states during the transition. The population of state increases to 9% in about picoseconds. The final two figures display the control and the populations of some bound states for the transition 22. This transition is 9% complete also in picoseconds. IV. CONCLUSIONS We have presented a family of state feedback laws for the control of quantum systems. These feedback laws generate arbitrarily bounded controls that asymptotically drive the system to a desired eigenstate of H. The region of attraction of the target eigenspace can be specified at will. An application to a model molecular system was given. V. REFERENCES [] R. Brockett, Lie Theory and Control Systems defined on Spheres, SIAM Journal of Applied Mathematics, vol. 25, no. 2, p. 23, 973. [2] C. Altafini, Controllability of quantum mechanical systems by root space decomposition of su(n), Journal of Mathematical Physics, vol. 43, no. 5, p. 25, 983. [3] V. Ramakrishna, K. Flores, H. Rabitz, and R. Ober, Quantum control by decompositions of SU(2), Physical Review A, vol. 62, no. 5349, 2. [4] W. Zhu and H. Rabitz, A rapid monotonically convergent iteration algorithm for quantum optimal control over the expectation value of a positive definite operator, Journal of Chemical Physics, vol. 9, no. 2, p. 385, 998. [5] S. Grivopoulos and B. Bamieh, Iterative algorithms for optimal control of quantum systems, in Proceedings of the 4st IEEE Conference on Decision and Control, Dec 22. [6] N. Khaneja, R. Brockett, and S. Glaser, Time optimal control in spin systems, Physical Review A, vol. 63, no. 3238, 2. [7] A. Ferrante, M. Pavon, and G. Raccanelli, Driving the propagator of a spin system: a feedback approach, in Proceedings of the 4st IEEE Conference on Decision and Control, Dec 22. [8] P. Vettori, On the convergence of a feedback control strategy for multilevel quantum systems, in Proceedings of the MTNS Conference, 22. [9] H. Khalil, Nonlinear Systems. Prentice Hall, second ed., 996.
4 control (a.u) Fig.. transition : electric field.9.8 populations Fig. 2. transition : state populations
5 control (a.u) Fig. 3. transition 22: electric field populations Fig. 4. transition 22: state populations
On the convergence of a feedback control strategy for multilevel quantum systems
On the convergence of a feedback control strategy for multilevel quantum systems Paolo Vettori Departamento de Matemática Universidade de Aveiro Campus de Santiago, 3810-193 Aveiro, Portugal email: pvettori@mat.ua.pt.
More informationOptimal Control of Quantum Systems
University of California, Santa Barbara Optimal Control of Quantum Systems PhD Thesis Symeon Grivopoulos 1 Average populations.8.6.4.2 1 3 4 6 9 1.1.2.3.4.5.6.7.8.9 1 Scaled time Profile intensities.1.2.3.4.5.6.7.8.9
More informationImplicit Lyapunov control of finite dimensional Schrödinger equations
Implicit Lyapunov control of finite dimensional Schrödinger equations Karine Beauchard a, Jean Michel Coron b, Mazyar Mirrahimi c, Pierre Rouchon c, a CMLA, ENS Cachan, avenue du president Wilson, 94230
More informationCONTROLLABILITY OF QUANTUM SYSTEMS. Sonia G. Schirmer
CONTROLLABILITY OF QUANTUM SYSTEMS Sonia G. Schirmer Dept of Applied Mathematics + Theoretical Physics and Dept of Engineering, University of Cambridge, Cambridge, CB2 1PZ, United Kingdom Ivan C. H. Pullen
More informationGeometric control for atomic systems
Geometric control for atomic systems S. G. Schirmer Quantum Processes Group, The Open University Milton Keynes, MK7 6AA, United Kingdom S.G.Schirmer@open.ac.uk Abstract The problem of explicit generation
More informationControl of quantum systems using model-based feedback strategies
Control of quantum systems using model-based feedback strategies Augusto Ferrante Dipartimento di Elettronica e Informatica Università di Padova via Gradenigo 6/A 35131 Padova, Italy. E-mail:augusto@dei.unipd.it
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationOptimal Control of Quantum Systems
University of California Santa Barbara Optimal Control of Quantum Systems A Dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Mechanical Engineering
More informationExplicit approximate controllability of the Schrödinger equation with a polarizability term.
Explicit approximate controllability of the Schrödinger equation with a polarizability term. Morgan MORANCEY CMLA, ENS Cachan Sept. 2011 Control of dispersive equations, Benasque. Morgan MORANCEY (CMLA,
More informationEnergy Level Sets for the Morse Potential
Energy Level Sets for the Morse Potential Fariel Shafee Department of Physics Princeton University Princeton, NJ 08540 Abstract: In continuation of our previous work investigating the possibility of the
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationLecture 2: Linear operators
Lecture 2: Linear operators Rajat Mittal IIT Kanpur The mathematical formulation of Quantum computing requires vector spaces and linear operators So, we need to be comfortable with linear algebra to study
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationState and Parameter Estimation Based on Filtered Transformation for a Class of Second-Order Systems
State and Parameter Estimation Based on Filtered Transformation for a Class of Second-Order Systems Mehdi Tavan, Kamel Sabahi, and Saeid Hoseinzadeh Abstract This paper addresses the problem of state and
More informationSymmetries, Groups, and Conservation Laws
Chapter Symmetries, Groups, and Conservation Laws The dynamical properties and interactions of a system of particles and fields are derived from the principle of least action, where the action is a 4-dimensional
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationAutomorphic Equivalence Within Gapped Phases
1 Harvard University May 18, 2011 Automorphic Equivalence Within Gapped Phases Robert Sims University of Arizona based on joint work with Sven Bachmann, Spyridon Michalakis, and Bruno Nachtergaele 2 Outline:
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationA Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,
More informationOn the Stabilization of Neutrally Stable Linear Discrete Time Systems
TWCCC Texas Wisconsin California Control Consortium Technical report number 2017 01 On the Stabilization of Neutrally Stable Linear Discrete Time Systems Travis J. Arnold and James B. Rawlings Department
More informationControl of bilinear systems: multiple systems and perturbations
Control of bilinear systems: multiple systems and perturbations Gabriel Turinici CEREMADE, Université Paris Dauphine ESF OPTPDE Workshop InterDyn2013, Modeling and Control of Large Interacting Dynamical
More informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More information1 Fundamental physical postulates. C/CS/Phys C191 Quantum Mechanics in a Nutshell I 10/04/07 Fall 2007 Lecture 12
C/CS/Phys C191 Quantum Mechanics in a Nutshell I 10/04/07 Fall 2007 Lecture 12 In this and the next lecture we summarize the essential physical and mathematical aspects of quantum mechanics relevant to
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationINTRODUCTION TO LIE ALGEBRAS. LECTURE 2.
INTRODUCTION TO LIE ALGEBRAS. LECTURE 2. 2. More examples. Ideals. Direct products. 2.1. More examples. 2.1.1. Let k = R, L = R 3. Define [x, y] = x y the cross-product. Recall that the latter is defined
More informationPath Integral for Spin
Path Integral for Spin Altland-Simons have a good discussion in 3.3 Applications of the Feynman Path Integral to the quantization of spin, which is defined by the commutation relations [Ŝj, Ŝk = iɛ jk
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationHybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching
Hybrid Control and Switched Systems Lecture #11 Stability of switched system: Arbitrary switching João P. Hespanha University of California at Santa Barbara Stability under arbitrary switching Instability
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
More informationA gradient-based approach to feedback control of quantum systems
A gradient-based approach to feedback control of quantum systems Gerasimos G. Rigatos Unit of Industrial Automation Industrial Systems Institute 26504, Rion Patras, Greece Email: grigat@isi.gr Abstract
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationL 2 -induced Gains of Switched Systems and Classes of Switching Signals
L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit
More informationThe Postulates of Quantum Mechanics
p. 1/23 The Postulates of Quantum Mechanics We have reviewed the mathematics (complex linear algebra) necessary to understand quantum mechanics. We will now see how the physics of quantum mechanics fits
More information4.3 Lecture 18: Quantum Mechanics
CHAPTER 4. QUANTUM SYSTEMS 73 4.3 Lecture 18: Quantum Mechanics 4.3.1 Basics Now that we have mathematical tools of linear algebra we are ready to develop a framework of quantum mechanics. The framework
More informationA Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets
A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets George A. Hagedorn Happy 60 th birthday, Mr. Fritz! Abstract. Although real, normalized Gaussian wave packets minimize the product
More information0.1 Schrödinger Equation in 2-dimensional system
0.1 Schrödinger Equation in -dimensional system In HW problem set 5, we introduced a simpleminded system describing the ammonia (NH 3 ) molecule, consisting of a plane spanned by the 3 hydrogen atoms and
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem
More informationLecture Notes 2: Review of Quantum Mechanics
Quantum Field Theory for Leg Spinners 18/10/10 Lecture Notes 2: Review of Quantum Mechanics Lecturer: Prakash Panangaden Scribe: Jakub Závodný This lecture will briefly review some of the basic concepts
More informationQuantum Computing Lecture 3. Principles of Quantum Mechanics. Anuj Dawar
Quantum Computing Lecture 3 Principles of Quantum Mechanics Anuj Dawar What is Quantum Mechanics? Quantum Mechanics is a framework for the development of physical theories. It is not itself a physical
More informationPart III Symmetries, Fields and Particles
Part III Symmetries, Fields and Particles Theorems Based on lectures by N. Dorey Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often
More informationQUANTUM FIELD THEORY: THE WIGHTMAN AXIOMS AND THE HAAG-KASTLER AXIOMS
QUANTUM FIELD THEORY: THE WIGHTMAN AXIOMS AND THE HAAG-KASTLER AXIOMS LOUIS-HADRIEN ROBERT The aim of the Quantum field theory is to offer a compromise between quantum mechanics and relativity. The fact
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationThe Framework of Quantum Mechanics
The Framework of Quantum Mechanics We now use the mathematical formalism covered in the last lecture to describe the theory of quantum mechanics. In the first section we outline four axioms that lie at
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationPhase control in the vibrational qubit
Phase control in the vibrational qubit THE JOURNAL OF CHEMICAL PHYSICS 5, 0405 006 Meiyu Zhao and Dmitri Babikov a Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin
More informationDifferential Topology Final Exam With Solutions
Differential Topology Final Exam With Solutions Instructor: W. D. Gillam Date: Friday, May 20, 2016, 13:00 (1) Let X be a subset of R n, Y a subset of R m. Give the definitions of... (a) smooth function
More informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
More informationC/CS/Phys 191 Quantum Mechanics in a Nutshell I 10/04/05 Fall 2005 Lecture 11
C/CS/Phys 191 Quantum Mechanics in a Nutshell I 10/04/05 Fall 2005 Lecture 11 In this and the next lecture we summarize the essential physical and mathematical aspects of quantum mechanics relevant to
More informationThe Principles of Quantum Mechanics: Pt. 1
The Principles of Quantum Mechanics: Pt. 1 PHYS 476Q - Southern Illinois University February 15, 2018 PHYS 476Q - Southern Illinois University The Principles of Quantum Mechanics: Pt. 1 February 15, 2018
More informationIntroduction to Quantum Spin Systems
1 Introduction to Quantum Spin Systems Lecture 2 Sven Bachmann (standing in for Bruno Nachtergaele) Mathematics, UC Davis MAT290-25, CRN 30216, Winter 2011, 01/10/11 2 Basic Setup For concreteness, consider
More informationLecture 4: Postulates of quantum mechanics
Lecture 4: Postulates of quantum mechanics Rajat Mittal IIT Kanpur The postulates of quantum mechanics provide us the mathematical formalism over which the physical theory is developed. For people studying
More informationConsensus Protocols for Networks of Dynamic Agents
Consensus Protocols for Networks of Dynamic Agents Reza Olfati Saber Richard M. Murray Control and Dynamical Systems California Institute of Technology Pasadena, CA 91125 e-mail: {olfati,murray}@cds.caltech.edu
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationFEEDBACK DESIGN OF CONTROL ALGORITHMS FOR DISSOCIATION OF DIATOMIC MOLECULES
FEEDBACK DESIGN OF CONTROL ALGORITHMS FOR DISSOCIATION OF DIATOMIC MOLECULES Alexander Efimov, Alexander Fradkov, Anton Krivtsov Institute for Problems of Mechanical Engineering, Russian Academy of Sciences,
More informationLecture 6 Positive Definite Matrices
Linear Algebra Lecture 6 Positive Definite Matrices Prof. Chun-Hung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Spring 2017 2017/6/8 Lecture 6: Positive Definite Matrices
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More informationOn complexified quantum mechanics and space-time
On complexified quantum mechanics and space-time Dorje C. Brody Mathematical Sciences Brunel University, Uxbridge UB8 3PH dorje.brody@brunel.ac.uk Quantum Physics with Non-Hermitian Operators Dresden:
More informationIntroduction to Quantum Computing
Introduction to Quantum Computing Petros Wallden Lecture 3: Basic Quantum Mechanics 26th September 2016 School of Informatics, University of Edinburgh Resources 1. Quantum Computation and Quantum Information
More informationContents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31
Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More information1 Time-Dependent Two-State Systems: Rabi Oscillations
Advanced kinetics Solution 7 April, 16 1 Time-Dependent Two-State Systems: Rabi Oscillations a In order to show how Ĥintt affects a bound state system in first-order time-dependent perturbation theory
More informationNumerical Linear and Multilinear Algebra in Quantum Tensor Networks
Numerical Linear and Multilinear Algebra in Quantum Tensor Networks Konrad Waldherr October 20, 2013 Joint work with Thomas Huckle QCCC 2013, Prien, October 20, 2013 1 Outline Numerical (Multi-) Linear
More informationHarmonic Oscillator I
Physics 34 Lecture 7 Harmonic Oscillator I Lecture 7 Physics 34 Quantum Mechanics I Monday, February th, 008 We can manipulate operators, to a certain extent, as we would algebraic expressions. By considering
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationL = e i `J` i `K` D (1/2,0) (, )=e z /2 (10.253)
44 Group Theory The matrix D (/,) that represents the Lorentz transformation (.4) L = e i `J` i `K` (.5) is D (/,) (, )=exp( i / /). (.5) And so the generic D (/,) matrix is D (/,) (, )=e z / (.53) with
More informationwhere P a is a projector to the eigenspace of A corresponding to a. 4. Time evolution of states is governed by the Schrödinger equation
1 Content of the course Quantum Field Theory by M. Srednicki, Part 1. Combining QM and relativity We are going to keep all axioms of QM: 1. states are vectors (or rather rays) in Hilbert space.. observables
More informationMic ael Flohr Representation theory of semi-simple Lie algebras: Example su(3) 6. and 20. June 2003
Handout V for the course GROUP THEORY IN PHYSICS Mic ael Flohr Representation theory of semi-simple Lie algebras: Example su(3) 6. and 20. June 2003 GENERALIZING THE HIGHEST WEIGHT PROCEDURE FROM su(2)
More informationMATRICES ARE SIMILAR TO TRIANGULAR MATRICES
MATRICES ARE SIMILAR TO TRIANGULAR MATRICES 1 Complex matrices Recall that the complex numbers are given by a + ib where a and b are real and i is the imaginary unity, ie, i 2 = 1 In what we describe below,
More informationVideo 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 8.1 Vijay Kumar 1 Definitions State State equations Equilibrium 2 Stability Stable Unstable Neutrally (Critically) Stable 3 Stability Translate the origin to x e x(t) =0 is stable (Lyapunov stable)
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationControl of the Inertia Wheel Pendulum by Bounded Torques
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 5 Seville, Spain, December -5, 5 ThC6.5 Control of the Inertia Wheel Pendulum by Bounded Torques Victor
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationarxiv: v1 [quant-ph] 4 Sep 2008
The Schroedinger equation as the time evolution of a deterministic and reversible process Gabriele Carcassi Brookhaven National Laboratory, Upton, NY 11973 (Dated: September 4, 008) In this paper we derive
More informationLinear Algebra 2 Spectral Notes
Linear Algebra 2 Spectral Notes In what follows, V is an inner product vector space over F, where F = R or C. We will use results seen so far; in particular that every linear operator T L(V ) has a complex
More informationA Glimpse of Quantum Computation
A Glimpse of Quantum Computation Zhengfeng Ji (UTS:QSI) QCSS 2018, UTS 1. 1 Introduction What is quantum computation? Where does the power come from? Superposition Incompatible states can coexist Transformation
More informationPassivity-based Stabilization of Non-Compact Sets
Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained
More informationGeorgia Institute of Technology Nonlinear Controls Theory Primer ME 6402
Georgia Institute of Technology Nonlinear Controls Theory Primer ME 640 Ajeya Karajgikar April 6, 011 Definition Stability (Lyapunov): The equilibrium state x = 0 is said to be stable if, for any R > 0,
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationDissipativity. Outline. Motivation. Dissipative Systems. M. Sami Fadali EBME Dept., UNR
Dissipativity M. Sami Fadali EBME Dept., UNR 1 Outline Differential storage functions. QSR Dissipativity. Algebraic conditions for dissipativity. Stability of dissipative systems. Feedback Interconnections
More informationDuke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014
Duke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014 Linear Algebra A Brief Reminder Purpose. The purpose of this document
More information04. Five Principles of Quantum Mechanics
04. Five Principles of Quantum Mechanics () States are represented by vectors of length. A physical system is represented by a linear vector space (the space of all its possible states). () Properties
More informationDOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS. Graziano Chesi
DOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS Graziano Chesi Dipartimento di Ingegneria dell Informazione Università di Siena Email: chesi@dii.unisi.it Abstract: Estimating the Domain
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationDifferential of the Exponential Map
Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements
More informationAppendix: SU(2) spin angular momentum and single spin dynamics
Phys 7 Topics in Particles & Fields Spring 03 Lecture v0 Appendix: SU spin angular momentum and single spin dynamics Jeffrey Yepez Department of Physics and Astronomy University of Hawai i at Manoa Watanabe
More informationMath 1553, Introduction to Linear Algebra
Learning goals articulate what students are expected to be able to do in a course that can be measured. This course has course-level learning goals that pertain to the entire course, and section-level
More informationPHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep
Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V
More informationC/CS/Phys C191 Quantum Mechanics in a Nutshell 10/06/07 Fall 2009 Lecture 12
C/CS/Phys C191 Quantum Mechanics in a Nutshell 10/06/07 Fall 2009 Lecture 12 In this lecture we summarize the essential physical and mathematical aspects of quantum mechanics relevant to this course. Topics
More informationVibrations and Rotations of Diatomic Molecules
Chapter 6 Vibrations and Rotations of Diatomic Molecules With the electronic part of the problem treated in the previous chapter, the nuclear motion shall occupy our attention in this one. In many ways
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More information