Variational Ansatz-Based Quantum Simulation of Imaginary Time Evolution
Total Page:16
File Type:pdf, Size:1020Kb
Variational ansatz-based quantum simulation of imaginary time evolution Sam McArdle,1 Tyson Jones,1 Suguru Endo,1 Ying Li,2 Simon Benjamin,1 and Xiao Yuan1, ∗ 1Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom 2Graduate School of China Academy of Engineering Physics, Beijing 100193, China (Dated: September 17, 2019) Imaginary time evolution is a powerful tool for studying quantum systems. While it is possible to simulate with a classical computer, the time and memory requirements generally scale exponentially with the system size. Conversely, quantum computers can efficiently simulate quantum systems, but not non-unitary imaginary time evolution. We propose a variational algorithm for simulating imaginary time evolution on a hybrid quantum computer. We use this algorithm to find the ground state energy of many-particle systems; specifically molecular Hydrogen and Lithium Hydride, finding the ground state with high probability. Our method can also be applied to general optimisation problems and quantum machine learning. As our algorithm is hybrid, suitable for error mitigation, and can exploit shallow quantum circuits, it can be implemented with current quantum computers. INTRODUCTION to represent the real time propagator with a sequence of Imaginary time is an unphysical, yet powerful, mathe- gates stems from its unitarity. In contrast, because the matical concept. It has been utilised in numerous phys- imaginary time operator is non-unitary, it is not straight- ical domains including: quantum mechanics, statistical forward to decompose it into a sequence of unitary gates mechanics, and cosmology. Often referred to as perform- using Trotterization, and thus directly realise it with a ing a `Wick rotation' [1], replacing real time with imag- quantum circuit. As a result, alternative methods are inary time connects Euclidean and Minkowski space [2], required to implement imaginary time evolution using a quantum and statistical mechanics [3], and static prob- quantum computer. lems to problems of dynamics [4]. In quantum mechan- Classically, we can simulate real (imaginary) time evo- ics, propagating a wavefunction in imaginary time en- lution of parametrised trial states by repeatedly solving ables: the study of finite temperature properties [5{7], the (Wick-rotated) Schr¨odingerequation over a small finding the ground state wavefunction and energy (such timestep, and updating the parameters for the next as in density matrix renormalisation group) [8{11], and timestep [8, 9, 11, 24{27]. This method has recently simulating real time dynamics (such as time dependent been extended to quantum computing, where it was used Hartree) [12, 13]. For a system with Hamiltonian, H, to simulate real time dynamics [28]. Closely related evolving in real time, t, the propagator is given by e−iHt. are the variational quantum eigensolver (VQE) [29{35] The corresponding propagator in imaginary time, τ = it, and the quantum approximate optimisation algorithm is given by e−Hτ ; a non-unitary operator. (QAOA) [36], which update the parameters using a clas- Using a classical computer, we can simulate imaginary sical optimisation routine, to find the minimum energy time evolution by evaluating the propagator and applying eigenvalue of a given Hamiltonian. As `hybrid quantum- it to the system wavefunction. There also exist various classical methods', these algorithms use a small quantum related classical methods, such as quantum Monte Carlo computer to carry out a classically intractable subrou- [14, 15] and density matrix renormalization group [16, 17] tine, and a classical computer to solve the higher level for solving different problems. However, because the di- problem. The quantum subroutine may only require a mension of the wavefunction grows exponentially with small number of qubits and a low depth circuit, present- the number of particles, classical simulation of many- ing a potential use for noisy intermediate-scale quantum body quantum systems is generally hard [18]. While ef- hardware [37]. arXiv:1804.03023v4 [quant-ph] 16 Sep 2019 ficient variational trial states have been developed for a In this paper, we propose a method to simulate number of applications [19], powerful trial wavefunctions imaginary time evolution on a quantum computer, using typically require classical computational resources which a hybrid quantum-classical variational algorithm. The scale exponentially with the system size [11]. proposed method thus combines the power of quantum Quantum computing can naturally and efficiently store computers to efficiently represent many-body quantum many-body quantum states, and hence is suitable for sim- states, with classical computers' ability to simulate ulating quantum systems [20]. We can map the system arbitrary (including unphysical) processes. We discuss Hamiltonian to a qubit Hamiltonian, and simulate real using this method to find the ground state energy of time evolution (as described by the Schr¨odingerequa- many-body quantum systems, and to solve optimisation tion) by realising the corresponding unitary evolution problems. We then numerically test the performance with a quantum circuit [21]. Using Trotterization [22], of our algorithm at finding the ground state energy of the real time propagator can be decomposed into a se- both the Hydrogen molecule (H2) and Lithium Hydride quence of single and two qubit gates [23]. The ability (LiH). We compare our results for LiH to those obtained 2 using the VQE with gradient descent. As our algorithm where kρk = Tr[pρρy] denotes the trace norm of a state. only requires a low depth circuit, it can be realised with By replacing j (τ)i with jφ(τ)i = jφ(θ~(τ))i, we effec- current and near-term quantum processors. tively project the desired imaginary time evolution onto the manifold of the ansatz space. The evolution of the pa- RESULTS rameters is obtained from the resulting differential equa- Variational imaginary time evolution. We focus on tion many-body systems that are described by Hamiltonians P X _ H = i λihi, with real coefficients, λi, and observ- Aijθj = Ci; (4) ables, hi, that are tensor products of Pauli matrices. j We assume that the number of terms in this Hamilto- nian scales polynomially with the system size, which is where true for many physical systems, such as molecules or the @ hφ(τ)j @ jφ(τ)i Fermi-Hubbard model. Given an initial state j i, the Aij = < ; normalised imaginary time evolution is defined by @θi @θj ! (5) X @ hφ(τ)j −Hτ C = < − λ h jφ(τ)i ; j (τ)i = A(τ)e j (0)i; (1) i α @θ α α i where A(τ) = 1=ph (0)j e−2Hτ j (0)i is a normalisation and hα and λα are the Pauli terms and coefficients of factor. In the instance that the initial state is a maxi- the Hamiltonian, as described above. The derivation of mally mixed state, the state at time τ is a thermal or Eq. (4) can be found in the Supplementary Materials. −Hτ −Hτ Gibbs state ρT =1/τ = e =Tr[e ], with temperature _ As both Aij and Ci are real, the derivative θj is also T = 1/τ. When the initial state has a non-zero overlap real, as required for parametrising a quantum circuit. with the ground state, the state at τ ! 1 is the ground Interestingly, although the average energy term Eτ state of H. Equivalently, the Wick rotated Schr¨odinger appears in Eq. (2), it does not appear in Eq. (4). This is equation is, because the ansatz applied maintains normalisation, as it is composed of unitary operators. @ j (τ)i = −(H − E ) j (τ)i ; (2) @τ τ Imaginary time evolution with quantum circuits. By where the term Eτ = h (τ)jHj (τ)i results from enforc- following a similar method to that introduced in ing normalisation. Even if j (τ)i can be represented by a Ref. [28], we can efficiently measure Aij and Ci us- quantum computer, the non-unitary imaginary time evo- ing a quantum computer. We assume that the deriva- lution cannot be naively mapped to a quantum circuit. tive of a unitary gate Ui(θi) can be expressed as P In our variational method, instead of directly encod- @Ui(θi)=@θi = k fk;iUi(θi)σk;i, with unitary op- ing the quantum state j (τ)i at time τ, we approxi- erator σk;i. The derivative of the trial state is P ~ ¯ ~ mate it using a parametrised trial state jφ(θ~(τ))i, with given by @ jφ(τ)i=@θi = k fk;iVk;i j0i, with Vk;i = ~ UN (θN ) :::Ui+1(θi+1)Ui(θi)σk;i :::U1(θ1). There are θ(τ) = (θ1(τ); θ2(τ); : : : ; θN (τ)). This stems from the in- tuition that the physically relevant states are contained typically only one or two terms resulting from each in a small subspace of the full Hilbert space [38]. The trial derivative. As an example, when Ui(θi) is a single qubit −iθiσz =2 state is referred to as the ansatz. In condensed matter rotation Rz(θi) = e , the derivative @Ui(θi)=@θi = −iθiσz =2 physics and computational chemistry, a wide variety of −i=2×σze . The coefficients Aij and Ci are given ans¨atzehave been proposed for both classical and quan- by tum variational methods [11, 20, 39, 40]. 0 1 Using a quantum circuit, we prepare the trial state, X A = < f ∗ f h0¯j V~ y V~ j0¯i ; jφ(θ~)i, by applying a sequence of parametrised unitary ij @ k;i l;j k;i l;j A k;l gates, V (θ~) = U (θ ) :::U (θ ) :::U (θ ) to our initial N N k k 1 1 0 1 (6) ~ ~ state, j0¯i. We express this as jφ(θ)i = V (θ) j0¯i and re- X C = < f ∗ λ h0¯j V~ y h V j0¯i : mark that V (θ~) is also referred to as the ansatz.