DOCSLIB.ORG

Efficient Estimation of Pauli Observables by Derandomization

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Efficient Estimation of Pauli Observables by Derandomization PHYSICAL REVIEW LETTERS 127, 030503 (2021) Efficient Estimation of Pauli Observables by Derandomization Hsin-Yuan Huang ,1,2,* Richard Kueng,3 and John Preskill1,2,4,5 1Institute for Quantum Information and Matter, Caltech, Pasadena, California 91125, USA 2Department of Computing and Mathematical Sciences, Caltech, Pasadena, California 91125, USA 3Institute for Integrated Circuits, Johannes Kepler University Linz, A-4040, Austria 4Walter Burke Institute for Theoretical Physics, Caltech, Pasadena, California 91125, USA 5AWS Center for Quantum Computing, Pasadena, California 91125, USA (Received 19 March 2021; accepted 14 June 2021; published 16 July 2021) We consider the problem of jointly estimating expectation values of many Pauli observables, a crucial subroutine in variational quantum algorithms. Starting with randomized measurements, we propose an efficient derandomization procedure that iteratively replaces random single-qubit measurements by fixed Pauli measurements; the resulting deterministic measurement procedure is guaranteed to perform at least as well as the randomized one. In particular, for estimating any L low-weight Pauli observables, a deterministic measurement on only of order logðLÞ copies of a quantum state suffices. In some cases, for example, when some of the Pauli observables have high weight, the derandomized procedure is substantially better than the randomized one. Specifically, numerical experiments highlight the advantages of our derandomized protocol over various previous methods for estimating the ground-state energies of small molecules. DOI: 10.1103/PhysRevLett.127.030503 Introduction.—Noisy intermediate-scale quantum (NISQ) randomized protocol that achieves our goal quite efficiently: devices are becoming available [1].Thoughlesspowerful For each of M copies of ρ, and for each of the n qubits, than fully error-corrected quantum computers, NISQ devices we chose uniformly at random to measure X, Y,orZ.Then used as coprocessors might have advantages over classical we can achieve the desired prediction accuracy with high computers for solving some problems of practical interest. success probability if M ¼ Oð3w log L=ε2Þ, assuming that For example, variational algorithms using NISQ hardware all L operators on our list have weight no larger than w have potential applications to chemistry, materials science, [11,12]. If the list contains high-weight operators, however, and optimization [2–10]. this randomized method is not likely to succeed unless M is In a typical NISQ variational algorithm, we need to very large. estimate expectation values for a specified set of operators In this Letter, we describe a deterministic protocol for fO1;O2; …;OLg in a quantum state ρ that can be prepared estimating Pauli-operator expectation values that always repeatedly using a programmable quantum system. To performs at least as well as the randomized protocol and obtain precise estimates, each operator must be measured performs much better in some cases. This deterministic many times, and finding a reasonably efficient procedure protocol is constructed by derandomizing the randomized for extracting the desired information is not easy in general. protocol. The key observation is that we can compute a In this Letter, we consider the special case where each Oj lower bound on the probability that randomized measure- is a Pauli operator; this case is of particular interest for ments on M copies successfully achieve the desired error ε near-term applications. for every one of our L target Pauli operators. Furthermore, Suppose we have quantum hardware that produces we can compute this lower bound even when the meas- multiple copies of the n-qubit state ρ. Furthermore, for urement protocol is partially deterministic and partially every copy, we can measure all the qubits independently, randomized; that is, when some of the measured single- choosing at our discretion to measure each qubit in the qubit Pauli operators are fixed, and others are still sampled X, Y,orZ basis. We are given a list of Ln-qubit Pauli uniformly from fX; Y; Zg. operators (each one a tensor product of n Pauli matrices), Hence, starting with the fully randomized protocol, we and our task is to estimate the expectation values of all L can proceed step by step to replace each randomized single- operators in the state ρ, with an additive error no larger than qubit measurement by a deterministic one, taking care in ε for each operator. We would like to perform this task each step to ensure that the new partially randomized using as few copies of ρ as possible. protocol, with one additional fixed measurement, has suc- If all L Pauli operators have relatively low weight (act cess probability at least as high as the preceding protocol. nontrivially on only a few qubits), there is a simple When all measurements have been fixed, we have a fully 0031-9007=21=127(3)=030503(7) 030503-1 © 2021 American Physical Society PHYSICAL REVIEW LETTERS 127, 030503 (2021) deterministic protocol. In numerical experiments, we find Lemma 1. (Confidence bound). Fix ε ∈ ð0; 1Þ (accu- that this deterministic protocol substantially outperforms racy) and 1 − δ ∈ ð0; 1Þ (confidence). Suppose that Pauli randomized protocols [12–16]. The improvement is espe- observables O ¼½o1; …; oL and Pauli measurements cially significant when the list of target observables includes P ¼½p1; …; pM are such that operators with relatively high weight. Further performance gains are possible by executing (at least) linear-depth circuits XL ε2 δ – O; P ≔ − h o ; P ≤ : before measurements [17 20]. Such procedures do, how- Confεð Þ exp 2 ð l Þ 2 ð3Þ ever, require deep quantum circuits. In contrast, our protocol l¼1 only requires single-qubit Pauli measurements, which are more amenable to execution on near-term devices. Then, the associated empirical averages (2) obey The manuscript is organized as follows. We first provide some statistical background, explain the randomized meas- jωˆ l − ωlðρÞj ≤ ε for all 1 ≤ l ≤ L ð4Þ urement protocol, then analyze the derandomization pro- cedure. We then provide numerical results showing how the with probability (at least) 1 − δ. derandomized protocol improves on previous methods. We See Supplemental Material Sec. B.1 for a detailed conclude with remarks and outlooks. Further examples and derivation [21]. We call the function defined in Eq. (3) details of proofs are in the Supplemental Material [21]. the “confidence bound.” It is a statistically sound summary Statistical background.—Let ρ be a fixed, but unknown, parameter that checks whether a set of Pauli measurements quantum state on n qubits. We want to accurately predict L (P) allows for confidently predicting a collection of Pauli expectation values observables (O) up to accuracy ε each. Randomized Pauli measurements.—Intuitively speaking, ωl ρ tr Oo ρ for 1 ≤ l ≤ L; 1 ð Þ¼ ð l Þ ð Þ a small confidence bound (3) implies a good Pauli estimation protocol. But how should we choose our M Pauli measure- O σ ⊗ ⊗ σ where each ol ¼ ol½1 ÁÁÁ ol½n is a tensor product ments (P) in order to achieve ConfεðO; PÞ ≤ δ=2?The of single-qubit Pauli matrices, i.e., ol ¼ [ol½1; …; ol½n] randomized measurement toolbox [12,13,16,22,23] provides with ol½k ∈ fI;X;Y;Zg. To extract meaningful informa- a perhaps surprising answer to this question. Let wðolÞ tion, we perform M (single-shot) Pauli measurements on denote the weight of Pauli observable ol, i.e., the number ρ 3n independent copies of . There are possible measure- Pof qubits on which the observable acts nontrivially: wðolÞ¼ n ment choices. Each of them is characterized by a full- k 1 1fol½k ≠ Ig. These weights capture the probability n ¼ weight Pauli string pm ∈ fX; Y; Zg and produces a of hitting ol with a completely random measurement string: n n q ∈ 1 wðolÞ random string of outcome signs m f g . Probp½ol▹p¼1=3 . In turn, a total of M randomly p 1 ≤ m ≤ M Not every Pauli measurement m ( )pro- selected Pauli measurements will, on average, achieve o vides actionable advice about every target observable l wðolÞ EP½hðol;PÞ¼M=3 hits, regardless of the actual Pauli (1 ≤ l ≤ L). The two must be compatible in the sense observable ol in question. This insight allows us to compute that the latter corresponds to a marginal of the former; expectation values of the confidence bound (3) i.e., it is possible to obtain ol from pm by replacing I some local nonidentity Pauli matrices with .Ifthisis L o ▹p p X the case, we write l m and say that measurement m wðolÞ M EP½ConfεðO; PÞ ¼ ð1 − ν=3 Þ ; ð5Þ “hits” target observable ol. For instance, X; I ; I;X ; ½ ½ l¼1 ½X; X▹½X; X,but½Z; I; ½I;Z; ½Z; Z ▹ ½X; X.Wecan approximate each ωlðρÞ by empirically averaging (appro- 2 where ν ¼ 1 − expð−ε =2Þ ∈ ð0; 1Þ. Each of the L terms is priately marginalized) measurement outcomes that exponentially suppressed in ε2M=3wðolÞ. Concrete realiza- belong to Pauli measurements that hit ol, tions of a randomized measurement protocol are extremely 1 X Y unlikely to deviate substantially from this expected behavior ωˆ q j ; (see, e.g., [11]). Combined with Lemma 1, this observation l ¼ h o ; p ; …; p m½ ð2Þ ð l ½ 1 MÞ m∶o ▹p l m j∶ol½j≠I implies a powerful error bound. Theorem 1. — P (Theorem 3 in Ref. [11].) Empirical h o ; p ;…;p M 1 o ▹p ∈ 0;1;…;M averages (2) obtained from M randomized Pauli measure- where ð l ½ 1 MÞ¼ m¼1 f l mg f g counts how many Pauli measurements hit target observ- ments allow for ε-accurately predicting L Pauli expectation O ρ ; …; O ρ ε able ol. values trð o1 Þ trð oL Þ up to additive error given wðolÞ 2 It is easy to check that each ωˆ l exactly reproduces ωlðρÞ that M ∝ logðLÞ maxl 3 =ε .
Load more

Recommended publications
  • Computational Complexity: a Modern Approach
    i Computational Complexity: A Modern Approach Draft of a book: Dated January 2007 Comments welcome! Sanjeev Arora and Boaz Barak Princeton University [email protected] Not to be reproduced or distributed without the authors’ permission This is an Internet draft. Some chapters are more finished than others. References and attributions are very preliminary and we apologize in advance for any omissions (but hope you will nevertheless point them out to us). Please send us bugs, typos, missing references or general comments to [email protected] — Thank You!! DRAFT ii DRAFT Chapter 9 Complexity of counting “It is an empirical fact that for many combinatorial problems the detection of the existence of a solution is easy, yet no computationally efficient method is known for counting their number.... for a variety of problems this phenomenon can be explained.” L. Valiant 1979 The class NP captures the difficulty of finding certificates. However, in many contexts, one is interested not just in a single certificate, but actually counting the number of certificates. This chapter studies #P, (pronounced “sharp p”), a complexity class that captures this notion. Counting problems arise in diverse fields, often in situations having to do with estimations of probability. Examples include statistical estimation, statistical physics, network design, and more. Counting problems are also studied in a field of mathematics called enumerative combinatorics, which tries to obtain closed-form mathematical expressions for counting problems. To give an example, in the 19th century Kirchoff showed how to count the number of spanning trees in a graph using a simple determinant computation. Results in this chapter will show that for many natural counting problems, such efficiently computable expressions are unlikely to exist.
    [Show full text]
  • The Complexity Zoo
    The Complexity Zoo Scott Aaronson www.ScottAaronson.com LATEX Translation by Chris Bourke [email protected] 417 classes and counting 1 Contents 1 About This Document 3 2 Introductory Essay 4 2.1 Recommended Further Reading ......................... 4 2.2 Other Theory Compendia ............................ 5 2.3 Errors? ....................................... 5 3 Pronunciation Guide 6 4 Complexity Classes 10 5 Special Zoo Exhibit: Classes of Quantum States and Probability Distribu- tions 110 6 Acknowledgements 116 7 Bibliography 117 2 1 About This Document What is this? Well its a PDF version of the website www.ComplexityZoo.com typeset in LATEX using the complexity package. Well, what’s that? The original Complexity Zoo is a website created by Scott Aaronson which contains a (more or less) comprehensive list of Complexity Classes studied in the area of theoretical computer science known as Computa- tional Complexity. I took on the (mostly painless, thank god for regular expressions) task of translating the Zoo’s HTML code to LATEX for two reasons. First, as a regular Zoo patron, I thought, “what better way to honor such an endeavor than to spruce up the cages a bit and typeset them all in beautiful LATEX.” Second, I thought it would be a perfect project to develop complexity, a LATEX pack- age I’ve created that defines commands to typeset (almost) all of the complexity classes you’ll find here (along with some handy options that allow you to conveniently change the fonts with a single option parameters). To get the package, visit my own home page at http://www.cse.unl.edu/~cbourke/.
    [Show full text]
  • The Computational Complexity Column by V
    The Computational Complexity Column by V. Arvind Institute of Mathematical Sciences, CIT Campus, Taramani Chennai 600113, India [email protected] http://www.imsc.res.in/~arvind The mid-1960’s to the early 1980’s marked the early epoch in the field of com- putational complexity. Several ideas and research directions which shaped the field at that time originated in recursive function theory. Some of these made a lasting impact and still have interesting open questions. The notion of lowness for complexity classes is a prominent example. Uwe Schöning was the first to study lowness for complexity classes and proved some striking results. The present article by Johannes Köbler and Jacobo Torán, written on the occasion of Uwe Schöning’s soon-to-come 60th birthday, is a nice introductory survey on this intriguing topic. It explains how lowness gives a unifying perspective on many complexity results, especially about problems that are of intermediate complexity. The survey also points to the several questions that still remain open. Lowness results: the next generation Johannes Köbler Humboldt Universität zu Berlin [email protected] Jacobo Torán Universität Ulm [email protected] Abstract Our colleague and friend Uwe Schöning, who has helped to shape the area of Complexity Theory in many decisive ways is turning 60 this year. As a little birthday present we survey in this column some of the newer results related with the concept of lowness, an idea that Uwe translated from the area of Recursion Theory in the early eighties. Originally this concept was applied to the complexity classes in the polynomial time hierarchy.
    [Show full text]
  • Physics 212B, Ken Intriligator Lecture 6, Jan 29, 2018 • Parity P Acts As P
    Physics 212b, Ken Intriligator lecture 6, Jan 29, 2018 Parity P acts as P = P † = P −1 with P~vP = ~v and P~aP = ~a for vectors (e.g. ~x • − and ~p) and axial pseudo-vectors (e.g. L~ and B~ ) respectively. Note that it is not a rotation. Parity preserves all the basic commutation relations, but it is not necessarily a symmetry of the Hamiltonian. The names scalar and pseudo scalar are often used to denote parity even vs odd: P P = for “scalar” operators (like ~x2 or L~ 2, and P ′P = ′ for pseudo O O O −O scalar operators (like S~ ~x). · Nature respects parity aside from some tiny effects at the fundamental level. E.g. • in electricity and magnetism, we can preserve Maxwell’s equations via V~ ~ V and → − V~ +V~ . axial → axial In QM, parity is represented by a unitary operator such that P ~r = ~r and • | i | − i P 2 = 1. So the eigenvectors of P have eigenvalue 1. A state is even or odd under parity ± if P ψ = ψ ; in position space, ψ( ~x) = ψ(~x). In radial coordinates, parity takes | i ±| i − ± cos θ cos θ and φ φ + π, so P α,ℓ,m =( 1)ℓ α,ℓ,m . This is clear also from the → − → | i − | i ℓ = 1 case, any getting general ℓ from tensor products. If [H,P ] = 0, then the non-degenerate energy eigenstates can be written as P eigen- states. Parity in 1d: PxP = x and P pP = p. The Hamiltonian respects parity if V ( x) = − − − V (x).
    [Show full text]
  • Constructive Galois Connections
    ZU064-05-FPR main 23 July 2018 12:36 1 Constructive Galois Connections DAVID DARAIS University of Vermont, USA (e-mail: [email protected]) DAVID VAN HORN University of Maryland, College Park, USA (e-mail: [email protected]) Abstract Galois connections are a foundational tool for structuring abstraction in semantics and their use lies at the heart of the theory of abstract interpretation. Yet, mechanization of Galois connections using proof assistants remains limited to restricted modes of use, preventing their general application in mechanized metatheory and certified programming. This paper presents constructive Galois connections, a variant of Galois connections that is effective both on paper and in proof assistants; is complete with respect toa large subset of classical Galois connections; and enables more general reasoning principles, including the “calculational” style advocated by Cousot. To design constructive Galois connections we identify a restricted mode of use of classical ones which is both general and amenable to mechanization in dependently-typed functional programming languages. Crucial to our metatheory is the addition of monadic structure to Galois connections to control a “specification effect.” Effectful calculations may reason classically, while pure calculations have extractable computational content. Explicitly moving between the worlds of specification and implementation is enabled by our metatheory. To validate our approach, we provide two case studies in mechanizing existing proofs from the literature: the first uses calculational abstract interpretation to design a static analyzer; the second forms a semantic basis for gradual typing. Both mechanized proofs closely follow their original paper-and-pencil counterparts, employ reasoning principles arXiv:1807.08711v1 [cs.PL] 23 Jul 2018 not captured by previous mechanization approaches, support the extraction of verified algorithms, and are novel.
    [Show full text]
  • Parity Detection of Propagating Microwave Fields
    Parity Detection of Propagating Microwave Fields Jean-Claude Besse,1, ∗ Simone Gasparinetti,1, y Michele C. Collodo,1 Theo Walter,1 Ants Remm,1 Jonas Krause,1 Christopher Eichler,1 and Andreas Wallraff1 1Department of Physics, ETH Zurich, CH-8093 Zurich, Switzerland (Dated: December 23, 2019) The parity of the number of elementary excitations present in a quantum system provides impor- tant insights into its physical properties. Parity measurements are used, for example, to tomographi- cally reconstruct quantum states or to determine if a decay of an excitation has occurred, information which can be used for quantum error correction in computation or communication protocols. Here we demonstrate a versatile parity detector for propagating microwaves, which distinguishes between radiation fields containing an even or odd number n of photons, both in a single-shot measurement and without perturbing the parity of the detected field. We showcase applications of the detector for direct Wigner tomography of propagating microwaves and heralded generation of Schr¨odinger cat states. This parity detection scheme is applicable over a broad frequency range and may prove useful, for example, for heralded or fault-tolerant quantum communication protocols. I. INTRODUCTION qubits for measurement based entanglement generation [13], for elements of error correction [14{16], and entan- In quantum physics, the parity P of a wavefunction glement stabilization [17], an experiment which was also governs whether a system has an even or odd number of performed with ions [18, 19]. excitations n. The parity affects, for example, the sys- tem's statistical properties such as the likelihood of tran- sitions occurring between distinct quantum states [1,2].
    [Show full text]
  • Lecture 24 1 the Complexity of Counting
    Notes on Complexity Theory Last updated: November, 2011 Lecture 24 Jonathan Katz 1 The Complexity of Counting We explore three results related to hardness of counting. Interestingly, at their core each of these results relies on a simple | yet powerful | technique due to Valiant and Vazirani. 1.1 Hardness of Unique-SAT Does SAT become any easier if we are guaranteed that the formula we are given has at most one solution? Alternately, if we are guaranteed that a given boolean formula has a unique solution does it become any easier to ¯nd it? We show here that this is not likely to be the case. De¯ne the following promise problem: USAT def= fÁ : Á has exactly one satisfying assignmentg USAT def= fÁ : Á is unsatis¯ableg: Clearly, this problem is in promise-NP. We show that if it is in promise-P, then NP = RP. We begin with a lemma about pairwise-independent hashing. n m m+1 Lemma 1 Let S ⊆ f0; 1g be an arbitrary set with 2 · jSj · 2 , and let Hn;m+2 be a family of pairwise-independent hash functions mapping f0; 1gn to f0; 1gm+2. Then Pr [there is a unique x 2 S with h(x) = 0m+2] ¸ 1=8: h2Hn;m+2 Proof Let 0 def= 0m+2, and let p def= 2¡(m+2). Let N be the random variable (over choice of random h 2 Hn;m+2) denoting the number of x 2 S for which h(x) = 0. Using the inclusion/exclusion principle, we have X 1 X Pr[N ¸ 1] ¸ Pr[h(x) = 0] ¡ ¢ Pr[h(x) = h(x0) = 0] 2 x2S x6=x02S µ ¶ jSj = jSj ¢ p ¡ p2; 2 ¡ ¢ P 0 jSj 2 while Pr[N ¸ 2] · x6=x02S Pr[h(x) = h(x ) = 0] = 2 p .
    [Show full text]
  • A Delayed Promotion Policy for Parity Games$
    A Delayed Promotion Policy for Parity GamesI Massimo Benerecettia, Daniele Dell’Erbaa, Fabio Mogaverob aUniversita` degli Studi di Napoli Federico II, Napoli, Italy bUniversity of Oxford, Oxford, UK Abstract Parity games are two-player infinite-duration games on graphs that play a crucial role in various fields of theoretical computer science. Finding efficient algorithms to solve these games in practice is widely acknowledged as a core problem in formal verification, as it leads to efficient solutions of the model-checking and satisfiability problems of expressive temporal logics, e.g., the modal µCALCULUS. Their solution can be reduced to the problem of identifying sets of positions of the game, called dominions, in each of which a player can force a win by remaining in the set forever. Recently, a novel technique to compute dominions, called priority promotion, has been proposed, which is based on the notions of quasi dominion, a relaxed form of dominion, and dominion space. The underlying framework is general enough to accommodate different instantiations of the solution procedure, whose correctness is ensured by the nature of the space itself. In this paper we propose a new such instantiation, called delayed promotion, that tries to reduce the possible exponential behaviors exhibited by the original method in the worst case. The resulting procedure often outperforms both the state-of-the-art solvers and the original priority promotion approach. Keywords: Parity games; Infinite-duration games on graphs; Formal methods. 1. Introduction The abstract concept of game has proved to be a fruitful metaphor in theoretical computer science [2]. Several decision problems can, indeed, be encoded as path- forming games on graphs, where a player willing to achieve a certain goal, usually the verification of some property on the plays derived from the original problem, has to face an opponent whose aim is to pursue the exact opposite task.
    [Show full text]
  • The Complexity of Partial-Observation Parity Games
    The Complexity of Partial-Observation Parity Games Krishnendu Chatterjee1 and Laurent Doyen2 1 IST Austria (Institute of Science and Technology Austria) 2 LSV, ENS Cachan & CNRS, France Abstract. We consider two-player zero-sum games on graphs. On the basis of the information available to the players these games can be classified as follows: (a) partial-observation (both players have partial view of the game); (b) one-sided partial-observation (one player has partial-observation and the other player has complete-observation); and (c) complete-observation (both players have com- plete view of the game). We survey the complexity results for the problem of de- ciding the winner in various classes of partial-observation games with ω-regular winning conditions specified as parity objectives. We present a reduction from the class of parity objectives that depend on sequence of states of the game to the sub-class of parity objectives that only depend on the sequence of observations. We also establish that partial-observation acyclic games are PSPACE-complete. 1 Introduction Games on graphs. Games played on graphs provide the mathematical framework to analyze several important problems in computer science as well as mathematics. In par- ticular, when the vertices and edges of a graph represent the states and transitions of a reactive system, then the synthesis problem (Church’s problem) asks for the construc- tion of a winning strategy in a game played on a graph [5,21,20, 18]. Game-theoretic formulations have also proved useful for the verification [1], refinement [13], and com- patibility checking [9] of reactive systems.
    [Show full text]
  • Lecture 1 1 Turing Machines
    Notes on Complexity Theory Last updated: August, 2011 Lecture 1 Jonathan Katz 1 Turing Machines I assume that most students have encountered Turing machines before. (Students who have not may want to look at Sipser's book [3].) A Turing machine is de¯ned by an integer k ¸ 1, a ¯nite set of states Q, an alphabet ¡, and a transition function ± : Q £ ¡k ! Q £ ¡k¡1 £ fL; S; Rgk where: ² k is the number of (in¯nite, one-dimensional) tapes used by the machine. In the general case we have k ¸ 3 and the ¯rst tape is a read-only input tape, the last is a write-once output tape, and the remaining k ¡2 tapes are work tapes. For Turing machines with boolean output (which is what we will mostly be concerned with in this course), an output tape is unnecessary since the output can be encoded into the ¯nal state of the Turing machine when it halts. ² Q is assumed to contain a designated start state qs and a designated halt state qh. (In the case where there is no output tape, there are two halting states qh;0 and qh;1.) ² We assume that ¡ contains f0; 1g, a \blank symbol", and a \start symbol". ² There are several possible conventions for what happens when a head on some tape tries to move left when it is already in the left-most position, and we are agnostic on this point. (Anyway, by our convention, below, that the left-most cell of each tape is \marked" there is really no reason for this to ever occur.
    [Show full text]
  • Arxiv:2103.07510V1 [Quant-Ph] 12 Mar 2021
    Efficient estimation of Pauli observables by derandomization Hsin-Yuan Huang,1, 2, ∗ Richard Kueng,3 and John Preskill1, 2, 4, 5 1Institute for Quantum Information and Matter, Caltech, Pasadena, CA, USA 2Department of Computing and Mathematical Sciences, Caltech, Pasadena, CA, USA 3Institute for Integrated Circuits, Johannes Kepler University Linz, Austria 4Walter Burke Institute for Theoretical Physics, Caltech, Pasadena, CA, USA 5AWS Center for Quantum Computing, Pasadena, CA, USA (Dated: March 16, 2021) We consider the problem of jointly estimating expectation values of many Pauli observables, a crucial subroutine in variational quantum algorithms. Starting with randomized measurements, we propose an efficient derandomization procedure that iteratively replaces random single-qubit measurements with fixed Pauli measurements; the resulting deterministic measurement procedure is guaranteed to perform at least as well as the randomized one. In particular, for estimating any L low-weight Pauli observables, a deterministic measurement on only of order log(L) copies of a quantum state suffices. In some cases, for example when some of the Pauli observables have a high weight, the derandomized procedure is substantially better than the randomized one. Specifically, numerical experiments highlight the advantages of our derandomized protocol over various previous methods for estimating the ground-state energies of small molecules. I. INTRODUCTION efficiently: For each of M copies of ρ, and for each of the n qubits, we chose uniformly at random to mea- Noisy Intermediate-Scale Quantum (NISQ) devices sure X, Y , or Z. Then we can achieve the desired prediction accuracy with high success probability if are becoming available [39]. Though less powerful w 2 than fully error-corrected quantum computers, NISQ M = O(3 log L/ ), assuming that all L operators devices used as coprocessors might have advantages on our list have weight no larger than w [15, 21].
    [Show full text]
  • Energy Parity Games✩
    Theoretical Computer Science 458 (2012) 49–60 Contents lists available at SciVerse ScienceDirect Theoretical Computer Science journal homepage: www.elsevier.com/locate/tcs Energy parity gamesI Krishnendu Chatterjee a, Laurent Doyen b,∗ a IST Austria (Institute of Science and Technology Austria), Austria b LSV, ENS Cachan & CNRS, France article info a b s t r a c t Article history: Energy parity games are infinite two-player turn-based games played on weighted graphs. Received 19 April 2011 The objective of the game combines a (qualitative) parity condition with the (quantitative) Received in revised form 14 May 2012 requirement that the sum of the weights (i.e., the level of energy in the game) must remain Accepted 30 July 2012 positive. Beside their own interest in the design and synthesis of resource-constrained Communicated by M. Mavronicolas omega-regular specifications, energy parity games provide one of the simplest model of games with combined qualitative and quantitative objectives. Our main results are as Keywords: follows: (a) exponential memory is sufficient and may be necessary for winning strategies Games on graphs Parity objectives in energy parity games; (b) the problem of deciding the winner in energy parity games Quantitative objectives can be solved in NP \ coNP; and (c) we give an algorithm to solve energy parity by reduction to energy games. We also show that the problem of deciding the winner in energy parity games is logspace-equivalent to the problem of deciding the winner in mean-payoff parity games, which can thus be solved in NP \ coNP. As a consequence we also obtain a conceptually simple algorithm to solve mean-payoff parity games.
    [Show full text]