IP Library Granted Patent US 12,260,341
Granted Patent B2
US 12,260,341 · App. 17/933,339 · Granted Mar 25, 2025

Quantum assisted optimization

Inventors: Vasil S. Denchev (West Lafayette, IN); Masoud Mohseni (Redondo Beach, CA); Hartmut Neven (Malibu, CA)
Assignee: Google LLC
G06N3/126G06N5/01G06N10/60
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,260,341
App. No.
17/933,339
Granted
Mar 25, 2025
Kind
B2
Abstract

Methods and apparatus for quantum assisted optimization. In one aspect, a method includes obtaining a set of initial input states, applying one or more of (i) dynamical thermal fluctuations and (ii) cluster update algorithms to the set of input states and subsequent input states when the states evolve within the classical information processors, applying dynamical quantum fluctuations to the set of input states and subsequent states when the states evolve within the quantum systems and repeating the application steps until a desirable output state is obtained.

Claims (56)

1. A method performed by classical information processors, the method comprising:

for a graph partitioned into one or more local regions, the graph representing a quantum system included in a quantum processor, the quantum system comprising a plurality of interacting qubits, wherein the plurality of qubits are represented by respective nodes in the graph, interactions between qubits are represented by respective edges between nodes, a ground state of the quantum system encodes a solution to an optimization task, and the graph is partitioned into the one or more local regions according to a division of the optimization task and embedding capabilities of the quantum processor:

repeatedly, until a target output state that approximates the ground state to within a predetermined threshold is obtained:

applying one or more of (i) dynamical thermal fluctuations or (ii) cluster update algorithms to an initial input state of the quantum system or a subsequent input state to obtain a classically evolved state for the repetition;

for each local region, determining whether application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms are failing to overcome one or more energy barriers in an energy landscape at a given temperature, wherein each energy barrier of the one or more energy barriers is associated with a respective local region of the one or more local regions;

for each local region for which it is determined that the application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms has failed to overcome one or more energy barriers, providing the classically evolved state for the repetition to the quantum processor for application of dynamical quantum fluctuations to the classically evolved state for the repetition, wherein the dynamical quantum fluctuations are only applied, at the given temperature, to qubits of the quantum processor that are represented by graph nodes in the local regions of the graph that are associated with the one or more energy barriers that application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms fail to overcome;

obtaining, from the quantum processor, an output state;

determining whether the output state is the target output state; and

in response to determining that the output state is not the target output state, providing the output state as an input state for a subsequent repetition; or

in response to determining that the output state is the target output state, determining a solution to the optimization task through measurement of the output state.

2. The method of claim 1 , wherein the dynamical thermal fluctuations comprise:

tempered transitions, wherein applying the tempered transitions comprise applying a parallel tempering algorithm; or

weighted dynamical tempered transitions, wherein applying the weighted dynamical tempered transitions comprise applying annealing importance sampling.

3. The method of claim 1 , wherein the cluster update algorithms create:

non-local state transformation in parameter space in various different temperatures; or

non-local isothermal state transformation in parameter space, wherein applying the cluster update algorithms comprise applying Houdayers cluster move algorithm.

4. The method of claim 1 , wherein the dynamical quantum fluctuations comprise increasing and decreasing zero-temperature quantum fluctuations or increasing and decreasing finite-temperature dissipative quantum fluctuations.

5. The method of claim 1 , further comprising partitioning the graph into local regions according to the optimization task.

6. The method of claim 1 , wherein determining whether the output state is the target output state comprises determining whether the output states converge to the predetermined threshold.

7. A system comprising one or more computers and one or more storage devices storing instructions that are operable, when executed by the one or more computers, to cause the one or more computers to perform operations comprising:

for a graph partitioned into one or more local regions, the graph representing a quantum system included in a quantum processor, the quantum system comprising a plurality of interacting qubits, wherein the plurality of qubits are represented by respective nodes in the graph, interactions between qubits are represented by respective edges between nodes, a ground state of the quantum system encodes a solution to an optimization task, and the graph is partitioned into the one or more local regions according to a division of the optimization task and embedding capabilities of the quantum processor:

repeatedly, until a target output state that approximates the ground state to within a predetermined threshold is obtained:

applying one or more of (i) dynamical thermal fluctuations or (ii) cluster update algorithms to an initial input state of the quantum system or a subsequent input state to obtain a classically evolved state for the repetition;

for each local region, determining whether application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms are failing to overcome one or more energy barriers in an energy landscape at a given temperature, wherein each energy barrier of the one or more energy barriers is associated with a respective local region of the one or more local regions;

for each local region for which it is determined that the application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms has failed to overcome one or more energy barriers, providing the classically evolved state for the repetition to the quantum processor for application of dynamical quantum fluctuations to the classically evolved state for the repetition, wherein the dynamical quantum fluctuations are only applied, at the given temperature, to qubits of the quantum processor that are represented by graph nodes in the local regions of the graph that are associated with the one or more energy barriers that application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms fail to overcome;

obtaining, from the quantum processor, an output state;

determining whether the output state is the target output state; and

in response to determining that the output state is not the target output state, providing the output state as an input state for a subsequent repetition; or

in response to determining that the output state is the target output state, determining a solution to the optimization task through measurement of the output state.

8. The system of claim 7 , wherein the dynamical thermal fluctuations comprise:

tempered transitions, wherein applying the tempered transitions comprise applying a parallel tempering algorithm; or

weighted dynamical tempered transitions, wherein applying the weighted dynamical tempered transitions comprise applying annealing importance sampling.

9. The system of claim 7 , wherein the cluster update algorithms create:

non-local state transformation in parameter space in various different temperatures; or non-local isothermal state transformation in parameter space, wherein applying the cluster update algorithms comprise applying Houdayers cluster move algorithm.

10. The system of claim 7 , wherein the dynamical quantum fluctuations comprise increasing and decreasing zero-temperature quantum fluctuations or increasing and decreasing finite-temperature dissipative quantum fluctuations.

11. The system of claim 7 , further comprising partitioning the graph into local regions according to the optimization task.

12. The system of claim 7 , wherein determining whether the output state is the target output state comprises determining whether the output states converge to the predetermined threshold.

13. A non-transitory computer-readable storage medium comprising instructions stored thereon that are executable by a processing device and upon such execution cause the processing device to perform operations for training a machine learning model to predict computer memory failures, the operations comprising:

for a graph partitioned into one or more local regions, the graph representing a quantum system included in a quantum processor, the quantum system comprising a plurality of interacting qubits, wherein the plurality of qubits are represented by respective nodes in the graph, interactions between qubits are represented by respective edges between nodes, a ground state of the quantum system encodes a solution to an optimization task, and the graph is partitioned into the one or more local regions according to a division of the optimization task and embedding capabilities of the quantum processor:

repeatedly, until a target output state that approximates the ground state to within a predetermined threshold is obtained:

applying one or more of (i) dynamical thermal fluctuations or (ii) cluster update algorithms to an initial input state of the quantum system or a subsequent input state to obtain a classically evolved state for the repetition;

for each local region, determining whether application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms are failing to overcome one or more energy barriers in an energy landscape at a given temperature, wherein each energy barrier of the one or more energy barriers is associated with a respective local region of the one or more local regions;

for each local region for which it is determined that the application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms has failed to overcome one or more energy barriers, providing the classically evolved state for the repetition to the quantum processor for application of dynamical quantum fluctuations to the classically evolved state for the repetition, wherein the dynamical quantum fluctuations are only applied, at the given temperature, to qubits of the quantum processor that are represented by graph nodes in the local regions of the graph that are associated with the one or more energy barriers that application of the (i) dynamical thermal fluctuations or (ii) cluster update algorithms fail to overcome;

obtaining, from the quantum processor, an output state;

determining whether the output state is the target output state; and

in response to determining that the output state is not the target output state, providing the output state as an input state for a subsequent repetition; or

in response to determining that the output state is the target output state, determining a solution to the optimization task through measurement of the output state.

14. The non-transitory computer-readable storage medium of claim 13 , wherein the dynamical thermal fluctuations comprise:

tempered transitions, wherein applying the tempered transitions comprise applying a parallel tempering algorithm; or

weighted dynamical tempered transitions, wherein applying the weighted dynamical tempered transitions comprise applying annealing importance sampling.

15. The non-transitory computer-readable storage medium of claim 13 , wherein the cluster update algorithms create:

non-local state transformation in parameter space in various different temperatures; or

non-local isothermal state transformation in parameter space, wherein applying the cluster update algorithms comprise applying Houdayers cluster move algorithm.

16. The non-transitory computer-readable storage medium of claim 13 , wherein the dynamical quantum fluctuations comprise increasing and decreasing zero-temperature quantum fluctuations or increasing and decreasing finite-temperature dissipative quantum fluctuations.

17. The non-transitory computer-readable storage medium of claim 13 , further comprising partitioning the graph into local regions according to the optimization task.

18. The non-transitory computer-readable storage medium of claim 13 , wherein determining whether the output state is the target output state comprises determining whether the output states converge to the predetermined threshold.

Assignments (3)
CORRECTIVE ASSIGNMENT TO CORRECT THE CONVEYING PARTY EXECUTION DATE FROM 9/30/2017 TO 9/29/2017 PREVIOUSLY RECORDED ON REEL 61561 FRAME 918. ASSIGNOR(S) HEREBY CONFIRMS THE CHANGE OF NAME. Recorded Dec 6, 2024
From: GOOGLE INC.
To: GOOGLE LLC
Reel/Frame 069532/0617 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Sep 28, 2022
From: DENCHEV, VASIL S.; MOHSENI, MASOUD; NEVEN, HARTMUT
To: GOOGLE INC.
Reel/Frame 061244/0750 →
CHANGE OF NAME Recorded Sep 28, 2022
From: GOOGLE INC.
To: GOOGLE LLC
Reel/Frame 061561/0918 →
Continuity (3)
Continuation 16096237
Provisional Application 62327384 · Apr 25, 2016
Related Publication 20230008626A1 · Jan 12, 2023
References Cited (29)
US 7008806B1 · Zhao et al. · 2006 [cited by applicant]
US 7875876B1 · Wandzura et al. · 2011 [cited by applicant]
US 9663358B1 · Cory · 2017 [cited by examiner]
US 20060225165A1 · Maassen van den Brink · 2006 [cited by examiner]
US 20070174227A1 · Johnson · 2007 [cited by examiner]
US 20080313114A1 · Rose · 2008 [cited by examiner]
US 20130290919A1 · Narayanaswamy · 2013 [cited by examiner]
US 20140297247A1 · Troyer et al. · 2014 [cited by applicant]
US 20150269124A1 · Hannze et al. · 2015 [cited by applicant]
US 20150363708A1 · Amin et al. · 2015 [cited by applicant]
US 20160267032A1 · Rigetti · 2016 [cited by examiner]
US 20160307649A1 · Yazdanbod · 2016 [cited by examiner]
US 20160335558A1 · Bunyk et al. · 2016 [cited by applicant]
US 20170255629A1 · Thom · 2017 [cited by examiner]
US 20170337155A1 · Novotny · 2017 [cited by examiner]
CN 1672171 · 2005 [cited by applicant]
CN 102082662 · 2011 [cited by applicant]
CN 104063623 · 2014 [cited by applicant]
WO WO2013006836 · 2013 [cited by applicant]
WO WO2015143439 · 2015 [cited by applicant]
MIT Student (and MZB), Lecture 20—Quantum Tunneling of Electrons, Mar. 20, 2009, 10 pages (Year: 2009). [cited by examiner]
Chancellor, “Runback experiments on the D-Wave device” Presentation at Durham University, Jun. 2015, 16 pages. [cited by applicant]
Earl et al., “Parallel Tempering: Theory, Applications, and New Perspectives” Submitted on Aug. 2005, arXiv:physics/0508111v2, 21 pages. [cited by applicant]
Houdayer, “A Cluster Monte Carlo Algorithm for 2-Dimensional Spin Glasses” Submitted on May 2001, arXiv: cond-mat/0101116v3, 6 pages. [cited by applicant]
International Search Report and Written Opinion in International Appln. No. PCT/US2016/069381, mailed on Apr. 20, 2017, 17 pages. [cited by applicant]
International Written Opinion in International Appln. No. PCT/US2016/069381, mailed on Mar. 28, 2018, 8 pages. [cited by applicant]
Office Action in Canadian Appln. No. 3022037, dated Jul. 14, 2020, 5 pages. [cited by applicant]
Office Action in Chinese Appln. No. 201680086099.9, dated May 21, 2021, 35 pages (with English translation). [cited by applicant]
Zhu et al., “Efficient Cluster Algorithm for Spin Glasses in Any Space Dimension” Phys. Rev. Lett. 115, 077201, Aug. 14, 2015, 5 pages. [cited by applicant]