Abstract

Simulating the time-evolution of quantum mechanical systems is BQP-hard and expected to be one of the foremost applications of quantum computers. We consider classical algorithms for the approximation of Hamiltonian dynamics using subsampling methods from randomized numerical linear algebra. We derive a simulation technique whose runtime scales polynomially in the number of qubits and the Frobenius norm of the Hamiltonian. As an immediate application, we show that sample based quantum simulation, a type of evolution where the Hamiltonian is a density matrix, can be efficiently classically simulated under specific structural conditions. Our main technical contribution is a randomized algorithm for approximating Hermitian matrix exponentials. The proof leverages a low-rank, symmetric approximation via the Nyström method. Our results suggest that under strong sampling assumptions there exist classical poly-logarithmic time simulations of quantum computations.

Highlights

  • Special purpose quantum simulators permit the efficient implementation of unitary dynamics governed by physically meaningful families of Hamiltonians, while the general task is BQP-hard – since we can implement any quantum computation by a sequence of Hamiltonian evolutions

  • Recall that in a strong quantum simulation one computes the amplitude of a particular outcome, whereas in a weak simulation one only samples from the output distribution of a quantum circuit

  • Our approach is tackling a problem that is harder than the one solved by quantum computers

Read more

Summary

Introduction

Special purpose quantum simulators permit the efficient implementation of unitary dynamics governed by physically meaningful families of Hamiltonians, while the general task is BQP-hard – since we can implement any quantum computation by a sequence of Hamiltonian evolutions. An important class of methods exploits the entanglement structure of the states to obtain efficient representations, known as tensor networks [VMC08; Orú14], that can be evolved either via Trotterisation [Vid04] or with a time-dependent variational principle [VGC04; Hae+11]. We conclude this introduction by noting that, from a computational perspective, the problem of simulating Hamiltonian dynamics is related, but different, to the problem of simulating quantum circuits. Recent techniques developed for this problem involve the notion of stabilizer rank [Bra+18] or neural network quantum states [JBC18]

Results
Applications and outlook
Classical and quantum approximation of matrix exponentials
Randomized numerical linear algebra and Nyström methods
Organization
Preliminaries
From row-searchability to efficient row-sampling
Algorithm for PSD row-searchable Hermitian matrices
Algorithm for row-searchable Hermitian matrices
F 2 in Eq 25 was bounded using the following inequality log
Application to density matrix simulation
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.