Quantum Simulation: Solving Schrödinger Equation on a Quantum Computer

Similar documents
The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.

Quantum Annealing for Heisenberg Spin Chains

5.74 TIME-DEPENDENT QUANTUM MECHANICS

C/CS/Phys C191 Deutsch and Deutsch-Josza algorithms 10/20/07 Fall 2007 Lecture 17

Classification of problem & problem solving strategies. classification of time complexities (linear, logarithmic etc)

Quantum Information & Quantum Computation

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Quantum Computing Lecture 7. Quantum Factoring

C191 - Lecture 2 - Quantum states and observables

Notes on Quantum Computing

1 Adiabatic and diabatic representations

Lanczos-Haydock Recursion

l -State Solutions of a New Four-Parameter 1/r^2 Singular Radial Non-Conventional Potential via Asymptotic Iteration Method

Lecture 11: Pseudorandom functions

True Nature of Potential Energy of a Hydrogen Atom

Analysis of Deutsch-Jozsa Quantum Algorithm

Last time, we talked about how Equation (1) can simulate Equation (2). We asserted that Equation (2) can also simulate Equation (1).

Chapter 4. Fourier Series

Lecture 25 (Dec. 6, 2017)

1. Szabo & Ostlund: 2.1, 2.2, 2.4, 2.5, 2.7. These problems are fairly straightforward and I will not discuss them here.

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

arxiv: v1 [math.co] 23 Mar 2016

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Chapter 5 Vibrational Motion

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.

Lecture #3. Math tools covered today

Quantum Signal Processing and Non-Orthogonal State Detection

Kinetics of Complex Reactions

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields

Multiple Groenewold Products: from path integrals to semiclassical correlations

3/21/2017. Commuting and Non-commuting Operators Chapter 17. A a

NUMERICAL METHODS FOR SOLVING EQUATIONS

Chapter 8 Approximation Methods, Hueckel Theory

Quantum Computing - A new Implementation of Simon Algorithm for 3-Dimensional Registers

) +m 0 c2 β K Ψ k (4)

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions

Algorithm of Superposition of Boolean Functions Given with Truth Vectors

DEGENERACY AND ALL THAT

NEW IDENTIFICATION AND CONTROL METHODS OF SINE-FUNCTION JULIA SETS

Signal Processing in Mechatronics

Math 155 (Lecture 3)

Matsubara-Green s Functions

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Exercises and Problems

PHY4905: Nearly-Free Electron Model (NFE)

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Office: JILA A709; Phone ;

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

TIME-CORRELATION FUNCTIONS

REPRESENTING MARKOV CHAINS WITH TRANSITION DIAGRAMS

010 Zeg Ke et al Vol. 14 level j+i ad a lower level j i. The upper level j+i is coupled to the lower level j i by the sigle-mode vacuum field. I the i

Some New Iterative Methods for Solving Nonlinear Equations

BIOSTATISTICS. Lecture 5 Interval Estimations for Mean and Proportion. dr. Petr Nazarov

IP Reference guide for integer programming formulations.

Addition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

PROPERTIES OF AN EULER SQUARE

EE260: Digital Design, Spring n MUX Gate n Rudimentary functions n Binary Decoders. n Binary Encoders n Priority Encoders

A Block Cipher Using Linear Congruences

Chimica Inorganica 3

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

CSE 1400 Applied Discrete Mathematics Number Theory and Proofs

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

Stability Analysis of the Euler Discretization for SIR Epidemic Model

ADVANCED DIGITAL SIGNAL PROCESSING

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan


Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.

1 Introduction: within and beyond the normal approximation

Modified Logistic Maps for Cryptographic Application

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods

Computability and computational complexity

Essential Question How can you recognize an arithmetic sequence from its graph?

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

A Note on the Symmetric Powers of the Standard Representation of S n

Quantum computers can search rapidly by using almost any transformation

The Wave Function and Quantum Reality

Building an NMR Quantum Computer

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

Large Deviations in Quantum Information Theory

Matrix Operators and Functions Thereof

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

On Extracting Properties of Lie Groups from Their Lie Algebras

arxiv:quant-ph/ v1 21 Jan 1998

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

Machine Learning for Data Science (CS 4786)

FFTs in Graphics and Vision. The Fast Fourier Transform

Stochastic Matrices in a Finite Field

The Pendulum. Purpose

Hilbert Space Methods Used in a First Course in Quantum Mechanics

Disjoint set (Union-Find)

multiplies all measures of center and the standard deviation and range by k, while the variance is multiplied by k 2.

Chapter 7: The z-transform. Chih-Wei Liu

Research Article A Unified Weight Formula for Calculating the Sample Variance from Weighted Successive Differences

Physics 2D Lecture Slides Lecture 25: Mar 2 nd

PAPER : IIT-JAM 2010

Notes on iteration and Newton s method. Iteration

A new iterative algorithm for reconstructing a signal from its dyadic wavelet transform modulus maxima

Transcription:

Purdue Uiversity Purdue e-pubs Birc Poster Sessios Birc Naotechology Ceter 4-14-008 Quatum Simulatio: Solvig Schrödiger Equatio o a Quatum Computer Hefeg Wag Purdue Uiversity, wag10@purdue.edu Sabre Kais Birc Naotechology Ceter ad Departmet of Chemistry, Purdue Uiversity, ais@purdue.edu Ala Aspuru-Guzi Harvard Uiversity Follow this ad additioal wors at: http://docs.lib.purdue.edu/aoposter Wag, Hefeg; Kais, Sabre; ad Aspuru-Guzi, Ala, "Quatum Simulatio: Solvig Schrödiger Equatio o a Quatum Computer" (008). Birc Poster Sessios. Paper 18. http://docs.lib.purdue.edu/aoposter/18 This documet has bee made available through Purdue e-pubs, a service of the Purdue Uiversity Libraries. Please cotact epubs@purdue.edu for additioal iformatio.

Quatum Simulatio: Solvig Schrödiger Equatio o a Quatum Computer Hefeg Wag ad Sabre Kais Departmet of Chemistry ad Birc Naotechology Ceter Ala Aspuru-Guzi Departmet of Chemistry ad Chemical Biology, Harvard

Why Quatum Computig? Come forth ito the light of thigs, let ature be your teacher --W. Wordsworth o Why ot? If ature is quatum mechaical o Rolf Ladauer: Iformatio is physical o Quatum computatio, quatum commuicatio, quatum teleportatio, superdese codig, quatum cryptography, etc.

Bit vs. Qubit Classical Bit State 0 or 1 Measuremets does t chage the state Determiistic result Ca mae a copy of bit Qubit State 0>, 1> or Superpositio Measuremets chage the system Probability results Ca ot cloe the qubit

Power of Quatum Computig Quatum parallelism: superpositio priciple 1 i 0 c i i A 1 i 0 c i A i Quatum etaglemet: olocal correlatio c 1, 1, 1 Bell state: 00 00 + 11 Jozef Grusa, Quatum computig, 1999, McGraw-Hill Publishig Compay

Shor s Algorithm The difficulty of factorizatio uderpis the security of may commo methods of ecryptio

Quatum Simulatio Problems of simulatig quatum systems o a classical computer Expoetial complexity growth of quatum systems. The Schrödiger equatio, which i most cases is too complicated to be exactly solvable. As the Hilbert space icreases expoetially with the icrease of the size of the system, umerical methods lie Cofiguratio Iteractio (CI) are too expesive.

Cotiue 198, Feyma: The expoetial complexity of quatum systems might be put to good use to simulate dyamics of aother quatum system. Seth Lloyd: the coecture is correct, it is possible to build uiversal quatum simulators. For a fairly geeral class of quatum systems, especially discrete systems, it is possible to achieve a expoetial speedup. Richard P. Feyma, Simulatig physics with computers. Iter. J. Theo. Phys., 198, 1, 467-488. Seth Lloyd, Uiversal Quatum Simulators. Sciece, 1996, 73, 1073-1078.

Solvig Schrödiger Equatio The Schrödiger equatio The solutio i H ˆ t The eergy spectrum P P ( r, t) A u ( r) exp( ie t) () E A ( E ) () t ( r,0) ( r, t) A exp( ie t) E

Quatum Phase Estimatio [0,1) To estimate a umber give uitary operator U with eigevector, eigevalue Prepare two qubit registers: idex register ad target register U Perform the cotrolled- o the target register Iverse QFT o the first register, the mae a measuremet. e i

Quatum Algorithm to Obtai Eigeeergy through Phase Estimatio Costruct two registers: The idex register ad the target register, a approximate eigevector. 0 0 V a H 1 1 M M 1 0 V a CotrolU 1 M M 1 0 U V a U exp ( iht) V a c

Meas QFT M i M M e c M c M c U M 3., 1 0 1 0 1 0 1 1 ) ( 1 ) ( 1 i e c p V a c t E e e U i iht / ;

Basis Set Method HF wave fuctio: oe determiat Mappig the Foc space of the FCI wave fuctio to the Hilbert space of the qubits. Direct mappig: each qubit represets the fermioic occupatio state of a atomic orbital. Occupied or ot. Compact mappig: oly a subspace of Foc space with fixed electros is mapped oto the qubits. Compact (spi) mappig Ala Aspuru-Guzi, et. al Sciece, 005, 309, 1704-1707.

Problems of Usig HF Wave Fuctio as Iitial Guess Excited states have differet leadig cofiguratios For some complicated curves, several states have to be cosidered simultaeously i order to obtai a accurate descriptio. The success probability is small at regios far from the equilibrium geometry ad avoided crossig.

Geeral CI Scheme FCI is too expesive eve o a quatum computer for moderate size molecules. For most chemical problems, the iformatio of groud state ad a few low-lyig excited states is sufficiet. MCSCF provides good iitial guess wave fuctios, which ca be evolved efficietly to the CI wave fuctio. The success probability ca be icreased dramatically at critical regios by usig appropriate MCSCF wave fuctio as iitial guess.

MCSCF Wave Fuctio MCSCF wave fuctio MCSCF A K K K K A i ik i μ μ C μ i Averaged MCSCF wave fuctio E H E i w i E i

Implemet CI Based o MCSCF Wave Fuctio Mappig the Foc space of the averaged MCSCF wave fuctio to the Hilbert space of the qubits. More compact mappig, cosider the symmetry of the molecule, states belog to differet irreducible represetatio do ot iteract. Usig the MCSCF wave fuctio as the iitial guess, for either the groud or excited state. Probability of gettig the correspodig CI wave fuctio i i is MCSCF

Applicatio to Water Molecule Groud state: cc-pvdz basis set, siglet state, Cv symmetry. MCSCF method: freeze first two a1 orbitals, active space, 15 CSFs: MRCI usig the same active space, but the sigle ad double excitatios to the exteral space, 1387 CSFs.

Results Hefeg Wag, Sabre Kais, Ala Aspuru-Guzi ad Mar R. Hoffma, Quatum Algorithm for Solvig Schrödiger Equatio to obtai the eergy spectrum. submitted.

Cotiue

Cotiue

Summary Usig the MCSCF wave fuctio as the iitial guess ca improve the success probability dramatically, eve ust a few CSFs. The method ca be geeralized to geeral MRCI scheme ad the etire potetial eergy surface ca be explored for groud state ad excited states. Istead of startig from a sigle elemet as i the HF wave fuctio, MCSCF method starts from a small matrix. This maes the evolutio safer ad faster. The idea ca be geeralized to Fiite elemet method ad ay umerical method.