Arxiv:2105.12767V1 [Quant-Ph] 26 May 2021 III
Total Page:16
File Type:pdf, Size:1020Kb
Fault-Tolerant Quantum Simulations of Chemistry in First Quantization 1, 2, 3, 4, 5, 1, 1 1, Yuan Su, ∗ Dominic W. Berry, y Nathan Wiebe, z Nicholas Rubin, and Ryan Babbush x 1Google Research, Venice, CA 90291, United States 2Institute for Quantum Information and Matter, Caltech, Pasadena, CA 91125, United States 3Department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109, Australia 4Department of Computer Science, University of Toronto, ON M5S 2E4, Canada 5Pacific Northwest National Laboratory, Richland, WA 99354, United States (Dated: May 28, 2021) Quantum simulations of chemistry in first quantization offer some important advantages over approaches in second quantization including faster convergence to the continuum limit and the op- portunity for practical simulations outside of the Born-Oppenheimer approximation. However, since all prior work on quantum simulation of chemistry in first quantization has been limited to asymp- totic analysis, it has been impossible to directly compare the resources required for these approaches to those required for the more commonly studied algorithms in second quantization. Here, we com- pile, optimize and analyze the finite resources required to implement two first quantized quantum algorithms for chemistry from Babbush et al. [1] that realize block encodings for the qubitization and interaction picture frameworks of Low et al. [2,3]. The two algorithms we study enable simulation with gate complexities of Oe(η8=3N 1=3t + η4=3N 2=3t) and Oe(η8=3N 1=3t) where η is the number of electrons, N is the number of plane wave basis functions, and t is the duration of time-evolution (t is linearly inverse to target precision when the goal is to estimate energies). In addition to providing the first explicit circuits and constant factors for any first quantized simulation, and then introducing improvements which reduce circuit complexity by about a thousandfold over naive implementations for modest sized systems, we also describe new algorithms that asymptotically achieve the same scaling in a real space representation. Finally, we assess the resources required to simulate various molecules and materials and conclude that the qubitized algorithm will often be more practical than the interaction picture algorithm. We demonstrate that our qubitization algorithm often requires much less surface code spacetime volume for simulating millions of plane waves than the best second quantized algorithms require for simulating hundreds of Gaussian orbitals. CONTENTS I. Introduction 2 A. Background on representing molecular Hamiltonians in a first quantized plane wave basis3 B. Overview of the Hamiltonian simulation frameworks deployed in this work8 II. The qubitization based algorithm 12 A. Circuit implementation of qubitization 14 B. Block-encoding the kinetic operator 18 C. State preparation for U and V 20 D. The SEL operations for U and V 23 E. Total cost of constructing the qubitization operator 26 arXiv:2105.12767v1 [quant-ph] 26 May 2021 III. The interaction picture based algorithm 31 A. Implementation choices for the interaction picture algorithm 33 B. Exponentiating the kinetic operator 35 C. Implementing SEL and PREP 37 D. Total cost of the interaction picture simulation 38 IV. Results and discussion 43 ∗ corresponding author: [email protected] y corresponding author: [email protected] z corresponding author: [email protected] x corresponding author: [email protected] 2 A. The resources required to simulate real systems 43 B. Conclusion and outlook 49 Acknowledgements 52 References 52 A. Variable and function names 57 B. Background on Galerkin representations 59 1. Galerkin basis representations in first quantization 59 2. Galerkin basis representations in second quantization 60 C. Qubit costings 61 1. Qubit costings for qubitization approach 61 2. Qubit costings for interaction picture approach 63 D. Selection between T , U, and V 66 E. Detailed costing of preparation of k state 69 F. The modified phase gradient state 72 G. Complexity of binary squaring 73 H. Bound on time-discretization error 80 I. Interaction picture algorithm for simulating generic Hamiltonians 84 J. Simulating first quantized real space grid representations 88 I. INTRODUCTION The last several years have seen monumental developments in quantum algorithms, resulting in significant progress towards the goal of showing that a modest sized quantum computer can provide a decisive advantage for scientifically relevant problems. The expectation that quantum simulation might enable an exponential speedup for physical problems stretches back to Feynman's original proposal for quantum computers [4,5] and was first concretely proposed for solving the the central problem of quantum chemistry by Aspuru-Guzik et al.[6]. Since then, progress has accelerated on multiple fronts: improvements in general purpose simulation frameworks [2,3,7{15], more efficient algorithms for quantum chemistry [16{29], more advantageous-to- quantum-simulate representations of molecular Hamiltonians [30{38] and an improved understanding of the problems that genuinely require a quantum computer versus those that can be classically computed [39, 40]. However, a problem lurks within this discussion that has not received as much attention as it deserves. Specifically, most existing methods for simulating quantum chemistry (especially those leveraging simple basis functions like plane waves) become impractical as we scale to the continuum limit. This is because the second quantized simulations that are deployed in the vast majority of quantum computing algorithms for chemistry require a number of qubits that scales linearly with the number of spin-orbital basis functions. Dating back to early work by Kassal et al. [41], first quantized simulations of chemistry have been proposed as a method to ameliorate this problem. The central idea behind such simulations is to track the basis state that each particle is in rather than storing the occupancy of each basis state. This has two significant benefits: the number of qubits needed to represent the state scales logarithmically with the number of basis functions (although linearly with the number of particles) and also, such simulations are easy to adapt to cases where the entanglement between the electronic and nuclear subsystems is non-negligible. Despite the promise of these methods, all prior work on first quantized algorithms has stopped at the level of asymptotics. Our aim in this paper is to address this shortcoming by providing a full costing of the two leading quantum simulation algorithms for chemistry in a first quantized representation, both first described by Babbush et al. [1]. These algorithms comprise asymptotically efficient schemes for realizing the oracles required by the qubitization and interaction picture frameworks of Low et al. [2,3]. Here, we compile and optimize 3 these algorithms within an error-correctable gate set. We develop new techniques to reduce the complexity, resulting in about three orders of magnitude improvement over a naive implementation for even modest sized systems. Our detailed cost analysis reveals that despite the improved scaling of the interaction picture based method, the qubitization based method often proves to be more practical. We conduct numerical experiments to characterize the overhead of both approaches and find that they are highly competitive with second quantized approaches, often requiring significantly fewer resources to reach comparable accuracy. The layout of this paper is as follows. In the remainder of this section, we provide the formalism and context needed to understand the two most efficient schemes for simulating the molecular Hamiltonian in first quantization. Section II develops and analyzes our optimized implementation of the qubitization based algorithm in first quantization. Section III contains a similar analysis for the interaction picture based algorithm. Finally, Section IV numerically compares these algorithms (to each other and to other methods in the literature), estimating the resources required to simulate real systems in a first quantized basis. A. Background on representing molecular Hamiltonians in a first quantized plane wave basis The accurate simulation of quantum systems was Feynman's original vision for quantum computing [4]. There is perhaps no more natural application of this idea than the simulation of atoms and molecules, governed by the interactions of electrons and nuclei. Such systems give rise to the properties of most matter, determining everything from the rates of chemical reactions to the conductivity of materials. In the non- relativistic case, the dynamics of these systems are governed by the Coulomb Hamiltonian: η 2 η L η L 2 L X X X ζ` X 1 X X ζ`ζκ H = ri + r` + : (1) − 2 − R` ri 2 ri rj − 2m` 2 R` Rκ i=1 i=1 `=1 k − k i=j=1 k − k `=1 `=κ=1 k − k | {z } | {z } | 6 {z } | {z } | 6 {z } T U V Tnuclei Vnuclei Here we have used atomic units where ~ as well as the mass and charge of the electron are unity, and denotes the two-norm. In the above expression i; j represent electronic degrees of freedom and `; κ representk · k nuclear degrees of freedom; thus, ri represent the positions of electrons whereas R` represent the positions of nuclei, m` the atomic masses of nuclei, and ζ` the atomic numbers of nuclei. Throughout this work we will use η to denote the number of electrons in our simulation and L to denote the number of nuclei in our simulation. We keep a catalogue of the different symbols used throughout this paper in Appendix A. The methods we analyze and develop in this work are useful for simulations of the dynamics of (1) (treating both nuclear and electronic degrees of freedom as fully quantum) with only minor modifications. But in order to simplify and reduce our analysis, we will primarily focus on simulating (1) under the Born-Oppenheimer approximation: i.e., H(R) = T +U(R)+V +C(R) where C(R) is a constant given by the Coulomb repulsion of the point charges of \classical" nuclei with locations R`. Thus, the Born-Oppenheimer approximation assumes that one can decouple nuclear and electronic degrees of freedom, treating the former as essentially \classical"[42, 43].