IP Library › Granted Patent US 12,282,827
Granted Patent B2
US 12,282,827 · App. 18/036,252 · Granted Apr 22, 2025

Virtual distillation for quantum error mitigation

Inventors: William Huggins (Oakland, CA); Jarrod Ryan McClean (Marina Del Rey, CA)
Assignee: Google LLC
G06N10/70G06N10/20
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,282,827
App. No.
18/036,252
Granted
Apr 22, 2025
Kind
B2
Abstract

Methods, systems and apparatus for determining an error-mitigated expectation value of a target observable with respect to a noisy quantum state. In one aspect a method includes obtaining multiple copies of the noisy quantum state; performing measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of the target observable with respect to an entangled quantum state, wherein M?1 and eigenvalues corresponding to non-dominant eigenvectors of the noisy quantum state in the spectral decomposition of the entangled quantum state are suppressed exponentially in M; and using the computed expectation value of the target observable with respect to an entangled quantum state to determine the error-mitigated expectation value of the target observable with respect to the noisy quantum state.

Claims (61)

1. A method for determining an error-mitigated expectation value of a target observable with respect to a noisy quantum state, the method performed by a system comprising a classical processor and a quantum computer, the method comprising:

obtaining, by the quantum computer, multiple copies of the noisy quantum state, wherein the noisy quantum state comprises stochastic errors resulting from noise on the quantum computer;

performing, by the quantum computer, measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of the target observable with respect to an entangled quantum state, wherein the entangled quantum state is given by ρ M /Tr(ρ M ) where ρ represents the noisy quantum state, M≥1, and eigenvalues corresponding to non-dominant eigenvectors of the noisy quantum state in the spectral decomposition of the entangled quantum state are suppressed exponentially in M;

performing, by the quantum computer, measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of an identity operator with respect to the entangled quantum state; and

dividing, by the classical processor, the computed expectation value of the target observable with respect to the entangled quantum state by the expectation value of the identity operator with respect to the entangled quantum state to determine the error-mitigated expectation value of the target observable with respect to the noisy quantum state, wherein the error-mitigated expectation value of the target observable with respect to the noisy quantum state comprises a reconstructed expectation value of the target observable with respect to a purified version of the noisy quantum state and approximates an expectation value of the target observable computed by the quantum computer in an absence of the stochastic errors and noise on the quantum computer.

2. The method of claim 1 , wherein noise experienced by each copy of the noisy quantum state comprises a same form and strength.

3. The method of claim 1 , wherein performing measurements on tensor products of M copies of the noisy quantum state to compute expectation values of the target observable with respect to an entangled quantum state comprises:

for each of a first number of measurement repetitions, measuring a product of i) the target observable and ii) a cyclic shift operator with respect to a tensor product of M copies of the noisy quantum state to obtain a first measurement outcome for each repetition;

computing, using the first measurement outcomes obtained from the first number of measurement repetitions, an expectation value of the product of i) the target observable and ii) a cyclic shift operator with respect to the tensor product of M copies of the noisy quantum state; and

normalizing the expectation value of the product of i) the target observable and ii) a cyclic shift operator with respect to the tensor product of M copies of the noisy quantum state to obtain the error-mitigated expectation value of the target observable with respect to the noisy quantum state.

4. The method of claim 3 , wherein normalizing the expectation value of the target observable comprises:

for each of a second number of measurement repetitions, measuring the cyclic shift operator with respect to a tensor product of M copies of the noisy quantum state to obtain a second measurement outcome for each repetition;

computing, using the second measurement outcomes obtained from the second number of measurement repetitions, an expectation value of the cyclic shift operator with respect to a tensor product of M copies of the noisy quantum state; and

dividing the expectation value of the product of i) the target observable and ii) a cyclic shift operator with respect to the tensor product of M copies of the noisy quantum state by the expectation value of the cyclic shift operator with respect to the tensor product of M copies of the noisy quantum state.

5. The method of claim 1 , wherein performing measurements on tensor products of M copies of the noisy quantum state comprises performing parallel measurements.

6. The method of claim 1 , wherein M≥2 and the target observable acts on a single qubit.

7. The method of claim 6 , wherein performing measurements on tensor products of M copies of the noisy quantum state comprises, for each of multiple measurement repetitions, wherein the multiple measurement repetitions comprise K measurement repetitions:

applying a diagonalization operator to a tensor product of M copies of the noisy quantum state to obtain an evolved quantum state, wherein the diagonalization operator diagonalizes i) a cyclic shift operator and ii) a product of a symmetrized version of the target observable and the cyclic shift operator, and

measuring a product of the target observable and the cyclic shift operator with respect to the evolved quantum state to obtain a respective measurement outcome for the repetition for each qubit in the evolved quantum state.

8. The method of claim 7 , further comprising:

computing an expectation value of the product of the target observable and the cyclic shift operator with respect to the evolved quantum state using the measurement outcomes obtained from the K measurement repetitions;

computing an expectation value of the cyclic shift operator with respect to the evolved quantum state using the measurement outcomes obtained from the K measurement repetitions; and

dividing the expectation value of the product of the target observable and the cyclic shift operator with respect to the evolved quantum state by the expectation value of the cyclic shift operator with respect to the evolved quantum state.

9. The method of claim 7 , wherein application of the cyclic shift operator couples qubits in a first copy of the noisy quantum state to corresponding qubits in each other copy of the noisy quantum state.

10. The method of claim 1 , wherein performing measurements on tensor products of M copies of the noisy quantum state comprises performing serial measurements.

11. The method of claim 1 , wherein M≥2 and the target observable acts on two or more qubits, the target observable comprising multiple tensor products of one-qubit operators.

12. The method of claim 11 , wherein performing measurements on tensor products of M copies of the noisy quantum state comprises:

for each tensor product of one-qubit operators:

for each of a first number of measurement repetitions:

applying a first diagonalization operator to a tensor product of the M copies of the noisy quantum state to obtain an evolved quantum state, wherein the first diagonalization operator diagonalizes a product of i) the tensor product of one-qubit operators and ii) the cyclic shift operator, and

measuring a product of the tensor product of one-qubit operators and the cyclic shift operator with respect to the evolved quantum state to obtain a respective first measurement outcome for the repetition of the first number of repetitions for each qubit in the evolved quantum state;

for each of a second number of measurement repetitions:

applying a second diagonalization operator to a tensor product of the M copies of the noisy quantum state to obtain a second evolved quantum state, wherein the second diagonalization operator diagonalizes the cyclic shift operator, and

measuring the cyclic shift operator with respect to the second evolved quantum state to obtain a respective second measurement outcome for the repetition of the second number of repetitions for each qubit in the second evolved quantum state.

13. The method of claim 12 , further comprising:

computing an expectation value of a product of the target observable and the cyclic shift operator using the first measurement outcomes obtained from the first number of measurement repetitions,

computing an expectation value of the cyclic shift operator using the second measurement outcomes obtained from the second number of measurement repetitions; and

dividing the expectation value of the product of the target observable and the cyclic shift operator by the expectation value of the cyclic shift operator.

14. The method of claim 1 , wherein performing measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of the target observable with respect to an entangled quantum state comprise performing ancilla-assisted measurements.

15. The method of claim 14 , wherein the target observable comprises a single qubit operator or a multi-qubit operator.

16. The method of claim 14 , wherein M≥2 and performing ancilla assisted measurements comprises, for each of multiple measurement repetitions:

preparing two quantum registers, where each quantum register comprises the noisy quantum state and an ancilla qubit prepared in the 0 state;

performing a non-destructive measurement of a cyclic shift operator, comprising:

applying a Hadamard gate to the ancilla qubit,

applying a cyclic shift operator to the two quantum registers conditioned on the ancilla qubit being in the 1 state, and

measuring the ancilla qubit in the X basis to obtain a first respective measurement outcome; and

measuring qubits in the two quantum registers with respect to a symmetrized version of the target observable to obtain respective second measurement outcomes.

17. The method of claim 16 , further comprising:

computing an expectation value of the target observable using the second measurement outcomes; and

computing an expectation value of the cyclic shift operator using the first measurement outcomes; and

dividing the expectation value of the target observable by the expectation value of the cyclic shift operator.

18. An apparatus comprising:

one or more classical processors; and

one or more quantum computing devices in data communication with the one or more classical processors, wherein the quantum computing hardware comprises:

one or more qubit registers, each qubit register comprising one or more qubits, and

a plurality of control devices configured to operate the one or more qubit registers;

wherein the apparatus is configured to perform the method of determining an error-mitigated expectation value of a target observable with respect to a noisy quantum state, the method comprising:

obtaining, by the one or more quantum computing devices, multiple copies of the noisy quantum state, wherein the noisy quantum state comprises stochastic errors resulting from noise on the quantum computer;

performing, by the one or more quantum computing devices, measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of the target observable with respect to an entangled quantum state, wherein the entangled quantum state is given by ρ M /Tr(ρ M ), where ρ represents the noisy quantum state, M≥1, and eigenvalues corresponding to non-dominant eigenvectors of the noisy quantum state in the spectral decomposition of the entangled quantum state are suppressed exponentially in M;

performing, by the one or more quantum computing devices, measurements on tensor products of M copies of the noisy quantum state to compute an expectation value of an identity operator with respect to the entangled quantum state; and

dividing, by the one or more classical processors, the computed expectation value of the target observable with respect to the entangled quantum state by the expectation value of the identity operator with respect to the entangled quantum state to determine the error-mitigated expectation value of the target observable with respect to the noisy quantum state, wherein the error-mitigated expectation value of the target observable with respect to the noisy quantum state comprises a reconstructed expectation value of the target observable with respect to a purified version of the noisy quantum state and approximates an expectation value of the target observable computed by the quantum computer in an absence of the stochastic errors and noise on the quantum computer.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Aug 28, 2023
From: HUGGINS, WILLIAM; MCCLEAN, JARROD RYAN
To: GOOGLE LLC
Reel/Frame 064725/0210 →
Continuity (2)
Provisional Application 63112593 · Nov 11, 2020
Related Publication 20240013087A1 · Jan 11, 2024
References Cited (74)
US 20200175409A1 · Kandala et al. · 2020 [cited by applicant]
US 20200349457A1 · Low · 2020 [cited by examiner]
JP 2018530931 · 2018 [cited by applicant]
Author: Jordan Cotler et al.; Title: Quantum Virtual Cooling; Publisher: Physical Review X; Edition: Phys. Rev. X 9, 031013—Published Jul. 29, 2019 (Year: 2019). [cited by examiner]
Aharonov et al., “Fault-tolerant quantum computation with constant error” Proceedings of the twenty-ninth annual ACM symposium on Theory of computing, May 1997, 176-188. [cited by applicant]
Arute et al., “Observation of separated dynamics of charge and spin in the Fermi-Hubbard model” CoRR, Submitted on Oct. 2020, arXiv:2010.07965v1, 20 pages. [cited by applicant]
Arute et al., “Quantum supremacy using a programmable superconducting processor,” Nature 574, 2019, 505-510. [cited by applicant]
Bacon et al., “Efficient quantum circuits for schur and Clebsch-Gordan transforms” CoRR, Submitted on Aug. 2004, arXiv:quant-ph/0407082v4, 4 pages. [cited by applicant]
Banchi et al., “Entanglement entropy scaling in solid-state spin arrays via capacitance measurements” CoRR, Submitted on Jan. 2017, arXiv:1608.03970v2, 8 pages. [cited by applicant]
Barenco et al., “Stabilization of quantum computations by symmetrization” CoRR, Submitted on Apr. 1996, arXiv:quant-ph/9604028v1, 20 pages. [cited by applicant]
Berthiaume et al., “The stabilisation of quantum computations” Proceedings Workshop on Physics and Computation. PhysComp'94, Nov. 1994, 60-62. [cited by applicant]
Biamonte et al., “Tensor networks in a nutshell” CoRR, Submitted on Jul. 2017, arXiv:1708.00006v1, 34 pages. [cited by applicant]
Bonet-Monroig et al., “Low-cost error mitigation by symmetry verification” CoRR, Submitted on Jul. 2018, arXiv:1807.10050v1, 11 pages. [cited by applicant]
Bonet-Monroig et al., “Nearly optimal measurement scheduling for partial tomography of quantum states” CoRR, Submitted on Aug. 2019, arXiv:1908.05628v1, 9 pages. [cited by applicant]
Bravyi et al., “Mitigating measurement errors in multi-qubit experiments,” CoRR, Submitted on Jul. 2020, arXiv:2006.14044v2, 13 pages. [cited by applicant]
Bravyi et al., “Universal quantum computation with ideal clifford gates and noisy ancillas” CoRR, Submitted on Mar. 2004, arXiv:quant-ph/0403025v2, 15 pages. [cited by applicant]
Bridgeman et al., “Hand waving and interpretive dance: an introductory course on tensor networks” J. Phys. A: Math. Theor. 50, 223001, May 2017, 61 pages. [cited by applicant]
Brun, “Measuring polynomial functions of states” CoRR, Submitted on Jan. 2004, arXiv:quant-ph/0401067v1, 10 pages. [cited by applicant]
Cai, “Quantum Error Mitigation using Symmetry Expansion” CoRR, Submitted on Jan. 2021, arXiv:2101.03151v1, 15 pages. [cited by applicant]
Campbell, “Random compiler for fast hamiltonian simulation” CoRR, Submitted on Jun. 2019, arXiv:1811.08017v2, 11 pages. [cited by applicant]
Chen et al., “Detector tomography on IBM 5-qubit quantum computers and mitigation of an imperfect measurement” CoRR, Submitted on Apr. 2019, arXiv:1904.11935v1, 16 pages. [cited by applicant]
Childs et al., “Faster Quantum Simulation by Randomization” CoRR, Submitted on May 2018, arXiv:1805.08385v1, 19 pages. [cited by applicant]
Cirac et al., “Optimal purification of single qubits” CoRR, Submitted on Dec. 1998, arXiv:quant-ph/9812075v1, 4 pages. [cited by applicant]
Colless et al., “Computation of molecular spectra on a quantum processor with an Error-Resilient algorithm” Phys. Rev. X 8, 011021, Feb. 2018, 7 pages. [cited by applicant]
Cotler et al., “Quantum virtual cooling.” Physical Review X 9.3, Jul. 2019, 031013:1-10. [cited by applicant]
Cotler et al., “Quantum overlapping tomography” Phys. Rev. Lett. 124, 100401, Mar. 2020, 6 pages. [cited by applicant]
Czarnik et al., “Error mitigation with clifford quantum-circuit data” CoRR, Submitted on May 2020, arXiv:2005.10189v1, 9 pages. [cited by applicant]
Czarnik, “Qubit-Efficient Exponential Suppression of Errors” CoRR, Submitted on Mar. 2021, arXiv:2102.06056v2, 12 pages. [cited by applicant]
Dmitriev et al., “Quasi-one-dimensional anisotropic Heisenberg model in a transverse magnetic field” CoRR, Submitted on Jul. 2004, arXiv:cond-mat/0407203v1, 4 pages. [cited by applicant]
Ekert et al., “Direct estimations of linear and nonlinear functionals of a quantum state” CoRR, Submitted on Mar. 2002, arXiv:quant-ph/0203016v1, 4 pages. [cited by applicant]
Endo et al., “Mitigating algorithmic errors in a hamiltonian simulation” CoRR, Submitted on Aug. 2018, arXiv:1808.03623v1, 8 pages. [cited by applicant]
Endo et al., “Practical quantum error mitigation for Near-Future applications” CoRR, Submitted on Dec. 2017, arXiv:1712.09271v1, 17 pages. [cited by applicant]
Fowler et al., “Surface codes: Towards practical large-scale quantum computation” Phys. Rev. A 86, 032324, Sep. 2012, 54 pages. [cited by applicant]
Garcia-Escartin et al., “The SWAP test and Hong-Ou-Mandel effect are equivalent” CoRR, Submitted on Mar. 2013, arXiv:1303.6814v1, 10 pages. [cited by applicant]
Google AI Quantum and Collaborators, “Hartree-Fock on a superconducting qubit quantum computer” CoRR, Submitted on Sep. 2020, arXiv:2004.04174v4, 30 pages. [cited by applicant]
Haah et al., “Codes and protocols for distilling T, controlled-S, and toffoli gates” Quantum 2, 71, Jun. 2018, 29 pages. [cited by applicant]
Hastings et al., “Measuring renyi entanglement entropy with quantum monte carlo” CoRR, Submitted on Mar. 2010, arXiv:1001.2335v2, 4 pages. [cited by applicant]
Ho et al., “Enhancing Quantum Models of Stochastic Processes with Error Mitigation” CoRR, Submitted on May 2021, arXiv:2105.06448v1, 20 pages. [cited by applicant]
Horodecki et al., “Direct detection of quantum entanglement” CoRR, Submitted on Nov. 2001, arXiv:quant-ph/0111064v1, 5 pages. [cited by applicant]
Horsman et al., “Surface code quantum computing by lattice surgery” New J. Phys. 14, 123011, Dec. 2012, 27 pages. [cited by applicant]
Huggins et al., “Virtual Distillation for Quantum Error Mitigation” CoRR, Submitted on Aug. 2021, arXiv:2011.07064v3, 26 pages. [cited by applicant]
Huggins et al., “Efficient and noise resilient measurements for quantum chemistry on Near-Term quantum computers” CoRR, Submitted on Jul. 2019, arXiv:1907.13117v1, 8 pages. [cited by applicant]
International Search Report and Written Opinion in International Appln. No. PCT/US2021/058958, mailed on Feb. 23, 2023, 21 pages. [cited by applicant]
Islam et al., “Measuring entanglement entropy in a quantum many-body system,” Nature 528, Dec. 2015, 77-83. [cited by applicant]
Johri et al., “Entanglement spectroscopy on a quantum computer” CoRR, Submitted on Jul. 2017, arXiv:1707.07658v1, 7 pages. [cited by applicant]
Kandala et al., “Error mitigation extends the computational reach of a noisy quantum processor” Nature 567, Mar. 2019, 491-495. [cited by applicant]
Kitaev, “Quantum measurements and the abelian stabilizer problem,” CoRR, Submitted on Nov. 1995, arXiv:quant-ph/9511026v1, 22 pages. [cited by applicant]
Kivlichan et al., “Quantum simulation of electronic structure with linear depth and connectivity” CoRR, Submitted on Feb. 2018, arXiv:1711.04789v2, 8 pages. [cited by applicant]
Knill, “Quantum Computing with Realistically Noisy Devices” Nature 434, Mar. 2005, 39-44. [cited by applicant]
Koczor, “Exponential error suppression for Near-Term quantum devices” CoRR, Submitted on Nov. 2020, arXiv:2011.05942v1, 21 pages. [cited by applicant]
Krantz et al., “A quantum engineer's guide to superconducting qubits” Applied Physics Reviews 6, Jun. 2019, 021318-1-57. [cited by applicant]
Litinski, “A Game of Surface Codes: Large-Scale Quantum Computing with Lattice Surgery,” Quantum 3, 128, Mar. 2019, 37 pages. [cited by applicant]
Maciejewski et al., “Mitigation of readout noise in near-term quantum devices by classical post-processing based on detector tomography” Quantum 4, 257, Apr. 2020, 23 pages. [cited by applicant]
McArdle et al., “Error Mitigated digital quantum simulation” Phys. Rev. Lett. 122, May 2019, 180501-6. [cited by applicant]
McClean et al., “Decoding quantum errors with subspace expansions” Nat. Commun. 11, 636, Jan. 2020, 9 pages. [cited by applicant]
McClean et al., “Hybrid quantum-classical hierarchy for mitigation of decoherence and determination of excited states” Phys. Rev. A 95, 042308, Apr. 2017, 11 pages. [cited by applicant]
O'Brien et al., “Error mitigation via verified phase estimation” CoRR, Submitted on Oct. 2020, arXiv:2010.02538v1, 25 pages. [cited by applicant]
O'Gorman et al., “Generalized swap networks for near-term quantum computing.” CoRR, Submitted on May 2019, arXiv:1905.05118v1, 16 pages. [cited by applicant]
Orus, “A practical introduction to tensor networks: Matrix product states and projected entangled pair states” CoRR, Submitted on Jun. 2013, arXiv:1306.2164v1, 50 pages. [cited by applicant]
Ouyang et al., “Compilation by stochastic Hamiltonian sparsification” Quantum 4, 235, Feb. 2020, 17 pages. [cited by applicant]
Peres, “Error symmetrization in quantum computers” CoRR, Submitted on May 1996, arXiv:quant-ph/9605009v1, 3 pages. [cited by applicant]
Perturbation theory for linear operators, Springer Science & Business Media, Jun. 2013, 643 pages. [cited by applicant]
Preskill, “Lecture notes for physics 229: Quantum information and computation” California Insitute of Technology, Sep. 1998, 321 pages. [cited by applicant]
Preskill, “Quantum computing in the NISQ era and beyond” Quantum 2, 79, Aug. 2018, 20 pages. [cited by applicant]
Sagastizabal et al., Experimental error mitigation via symmetry verification in a variational quantum eigensolver Phys. Rev. A 100, Jul. 2019, 010302:1-6. [cited by applicant]
Strikis et al., “Learning-based quantum error mitigation” CoRR, Submitted on Mar. 2021, arXiv:2005.07601v2, 28 pages. [cited by applicant]
Subasi et al., “Entanglement spectroscopy with a depth-two quantum circuit” CoRR, Submitted on Jun. 2018, arXiv:1806.08863v1, 7 pages. [cited by applicant]
Temme et al., “Error mitigation for Short-Depth quantum circuits” CoRR, Submitted on Nov. 2017, arXiv:1612.02058v3, 15 pages. [cited by applicant]
Tran et al., “Faster Digital Quantum Simulation by Symmetry Protection” CoRR, Submitted on Jun. 2020, arXiv:2006.16248v1, 25 pages. [cited by applicant]
Wang et al., “Threshold error rates for the toric and surface codes” CoRR, Submitted on May 2009, arXiv:0905.0531v1, 14 pages. [cited by applicant]
Office Action in Australian Appln. No. 2021379675, mailed on Oct. 18, 2023, 3 pages. [cited by applicant]
Office Action in Australian Appln. No. 2021379675, mailed on Sep. 26, 2024, 3 pages. [cited by applicant]
Office Action in Canada Appln. No. 3,198,490, mailed on Oct. 8, 2024, 4 pages. [cited by applicant]
Office Action in Japanese Appln. No. 2023-528316, mailed on Jul. 16, 2024, 13 pages (with English translation). [cited by applicant]
Cited By (2)
US 12,682,272 US 12,743,646