Optimization with a Rydberg atom-based quantum processor

2 months ago 59

Amazon Braket recently launched the Aquila quantum processing unit (QPU) based on Rydberg atoms by QuEra Computing. In a previous post, we explored how researchers can use Rydberg devices to study problems and quantum phenomena in fundamental physics, for instance the emergence of a spin liquid phase [Semeghini et al., 2021]. But QuEra’s device also enables researchers to study quantum algorithms in new regimes. In this post, we review a recent publication, “Quantum optimization of maximum independent set using Rydberg atom arrays” [Ebadi et al., Science (2022)], in which the authors report the largest deployment of a variational algorithm to date, used to solve combinatorial optimization problems on a Rydberg device similar to the one launched on Nov 1 on Amazon Braket.

In the publication [Ebadi et al., Science (2022)], researchers from Harvard University, in collaboration with scientists from QuEra Computing, MIT, and the University of Innsbruck, reported on using Rydberg atom-based quantum processors to compute solutions for a particular graph coloring problem: finding the maximum independent set (MIS), which is relevant for optimal resource allocation, for example. By running a quantum-classical hybrid algorithm on the device, researchers were able to find the MIS for a class of graphs with a probability that scales favorably with the difficulty of the problem when compared to a classical simulated annealing approach. This work not only highlights the potential of Rydberg-atom devices for solving hard optimization problems, but also opens new research directions.

Maximum independent set

First, let us look at what an MIS problem really is. Edges between nodes of a graph may represent a wide variety of relationships, e.g., social ties between people, railroad connections between cities, or wiring between electrical components. Many optimization problems can be simplified to a task of coloring (or labeling) graph nodes while respecting certain rules.

Figure 1. Left: In a unit-disk graph, edges are created between nodes that are closer to each other than a certain fixed distance on the 2-d plane. Right: Red nodes mark a MIS of the graph corresponding to the largest number of nodes (six) that are not directly connected. [source: Ebadi et al., Science (2022)]

One such task is finding the largest set of nodes (colored red in Figure 1 above) such that no two nodes are directly connected. Such a set of red nodes is called an MIS of the graph. Computing an MIS has application, for instance, in planning the placement of cellphone towers, where no two towers may be closer than a given minimum distance due to signal interference. This particular example corresponds to a sub-class of MIS problems, the so-called unit-disk MIS, where the graph in question is formed by a set of nodes on a 2-dimensional plane and two nodes are directly connected if they are closer than a predefined distance (see Figure 1 above).

Finding an MIS is related to other computational tasks on graphs (e.g., finding cliques and dominant sets), and can be extended to sampling or enumerating all MISs of a graph. Many optimization problems in finance, bioinformatics, logistics, and automation can be mapped to such tasks [Wurtz et al., (2022)].

Discovering an MIS is a computationally hard problem even for the class of 2-d unit-disk graphs (Pichler et al., (2018)) and deciding whether an MIS of a certain size exists for such a graph is NP-complete. The running time of an exhaustive search increases exponentially with the size of the graph. Promising approaches to this include the “branch and bound” algorithm, which guarantees the best solutions but has exponential complexity in the worst case. This means that while these approaches may be fast for easy problem instances, for the hardest graphs (where the number of almost maximal independent sets is much larger than the number true MISs) they still take exponentially long to arrive at the solution. Another approach is to use heuristic algorithms, such as simulated annealing, which can compute independent sets quickly even for large graphs (with a few thousand nodes), and these provide good-enough practical solutions to many use cases, but these methods provide no guarantees that the result is actually an MIS. In particular, the chance of arriving at an MIS also decreases with the difficulty of the problem. Quantifying this difficulty for heuristic algorithms is one of the main achievements of this work. Another achievement of this work is showing that quantum effects can improve the algorithms’ likelihood of success.

Computing the maximum independent set with Rydberg atoms

Rydberg-atom quantum computers, like QuEra’s Aquila device on Amazon Braket, are well suited to study quantum algorithms for MIS problems on unit disk graphs. By arranging the individual atom using optical tweezers, the graph can be directly encoded in the atom positions with each graph node represented by an atom.

Figure 2. Left: Atoms (green spots) are arranged in a 2-d plane to realize a unit disk graph, where nodes are connected by an edge if they are closer than the Rydberg blockade radius (yellow circle). Middle: A coherent driving field with time-dependent amplitude and detuning  is used to steer the evolution of the system from a non-excited state to a maximally excited state constrained by Rydberg blockade () between nearby atoms. Right: An image of the final state is captured with an optical microscope. Missing spots (red circles) mark the locations of excited atoms, the collection of which corresponds to the discovered MIS of the graph. A classical optimizer is used to enhance the probability of finding an MIS. [source: Ebadi et al., Science (2022)]

There are two states of interest for each atom: atoms in the ground state represent white nodes, while those in the highly excited Rydberg state stand for red nodes. For a single shot of the MIS search, all atoms are initialized to be in the ground state (see left panel of Figure 2) and the system undergoes a coherent evolution governed by the time-dependent Hamiltonian  (see Eq. 1).

Equation 1. The driving Hamiltonian  induces coherent mixing between ground  and Rydberg  states of all atoms . The time-dependent driving amplitude  is ramped up from zero to high values at the beginning of the evolution and ramped down to zero at the end (See Figure 3, red curve). This mixing ensures that the cost of the joint state of all atoms stays low even as the cost Hamiltonian  is gradually changed by slowly increasing the detuning parameter  from negative to positive values (See Figure 3, blue curve).

Here the driving Hamiltonian  transitions the atoms between ground states and excited states for searching an MIS. The pairwise Rydberg interactions , however, prevents neighboring atoms from both being excited if their distance is within a certain range. Hence, this Rydberg “blockade” interaction naturally implements the rule that nodes connected by an edge cannot both be part of an independent set, and its range is illustrated with a yellow circle on Figure 2. To find an MIS, a high reward for Rydberg excitations is added to the objective Hamiltonian , which favors solutions with the highest number of red nodes. Throughout the evolution, the system remains in the approximate instantaneous ground state of the Hamiltonian. At the end of the evolution, the driving Hamiltonian is turned off, and an image of the final state is taken with an optical microscope, from which the value of  is computed.

An MIS of the graph corresponds to an atomic configuration with the lowest . To enhance the probability of finding an MIS, the evolution is repeated multiple times, to estimate the average  across shots, followed by using a classical optimizer to update the parameters of the quantum evolution. This variational optimization method continually updates the parameters of driving amplitude the detuning, until convergence where  is minimized.

Figure 3. The schedule of the amplitude and the detuning  of the driving field for the variational quantum adiabatic algorithm. The amplitude is a constant, except the ramping periods, and the smoothened piecewise-linear time series of the detuning parameter  are optimized with a closed-loop optimization. [source: Ebadi et al., Science (2022)]

Comparison with simulated annealing

In Ebadi et al., Science (2022), the authors tested different parametrizations of the schedules and found that the variational quantum adiabatic algorithm (VQAA, shown in Figure 3) can find the MIS with higher probability than the quantum approximate optimization algorithm (QAOA, where the detuning is kept at zero and the driving phase is changed in a stepwise fashion). They also compared VQAA to its classical counterpart, namely the simulated annealing algorithm. This algorithm works by simulating the thermal jitter of a fictitious mechanical system whose energy landscape matches the cost function in question. Throughout the algorithm, the temperature parameter is gradually decreased, causing the system to relax to low-energy states, thereby discovering the local optima of the cost function.

The authors found that the probability of finding an MIS for either VQAA or simulated annealing depends on the “hardness” of the graph. Intuitively, if the graph supports many MISs, it is easy for both algorithms to find an MIS; conversely, if the graph has only few MISs compared to a large number of sub-optimal independent sets with one less node , then both algorithms struggle. The authors compared the probability  with which the MIS is returned by the two methods for a range of graphs and for a fixed number of steps that both algorithms are allowed to perform. For some of the hardest graphs, with up to 80 nodes, they found that VQAA was able to approach the MIS in fewer steps, while simulated annealing remained stuck in local minima with fewer red nodes. For a particularly hard case, the quantum algorithm produced the correct MIS with 10 times higher probability than simulated annealing. To quantify the expected advantage of the VQAA on the quantum device over simulated annealing, the authors compared how  is decreasing with increasing hardness for the two algorithms: VQAA exhibited a more favorable scaling (with an exponent of -0.63) compared to simulated annealing (which scaled with an exponent of -1.03). The difference amounts to a superlinear improvement of success probability PMIS.

The probability of finding an MIS with the quantum algorithm was found to crucially depend on the minimal spectral gap: the smaller the energy difference between the two lowest energy states, the slower the adiabatic transition needs to be to guarantee a high probability of finding an MIS. Yet, if time evolution is too slow, decoherence and atom loss limits the performance. Specifically, the paper finds that the required time scales polynomially with the inverse of the minimum gap; the dependence of the gap size with problem size and specific instance is unknown. Still, this is an optimistic result; as technology evolves and coherence times improves, the ensemble of MIS problems for which these quantum processors offer an improved performance compared to simulated annealing is expected to increase.

Solving MIS with Rydberg simulator on Amazon Braket

Let us now take a look how we can put these ideas into action. In this section we will demonstrate how to create the graph from Figure 2, and solve the corresponding unit disc problem with Aquila on Braket. Using the Braket python SDK, the following code can be used to define the analog Hamiltonian simulation program that realizes an MIS search on a fixed unit disk graph with a fixed schedule and run it on QuEra’s Aquila QPU.

First, we define the atom arrangement that realizes the unit disk graph defined by our optimization problem. Later, we are going to program the time-dependent parameters so that they realize a blockade radius of about 6 micrometers, which means placing the atoms in a square grid with lattice spacing of 4 microns realizes the unit disk graph shown in Figure 2.

from braket.ahs.atom_arrangement import AtomArrangement separation = 4e-6  # in meters row_max = 14 col_max = 14 vacancies = { # (row_idx, col_idx) (0, 1), (0, 8), (0, 11), (0, 12), (1, 8), (1, 9), (2, 10), (3, 12), (3, 13), (3, 14), (4, 3), (4, 4), (4, 7), (4, 13), (4, 14), (5, 7), (5, 13), (6, 6), (6, 8), (6, 11), (6, 14), (7, 0), (7, 4), (7, 5), (7, 8), (7, 11), (8, 1), (8, 6), (8, 13), (9, 6), (9, 9), (10, 1), (10, 2), (10, 5), (10, 6), (10, 7), (11, 4), (11, 10), (12, 8), (13, 5), (14, 2), (14, 4), (14, 8), (14, 9), (14, 11) } register = AtomArrangement() for row in range(row_max + 1): for col in range(col_max + 1): if (row, col) not in vacancies: # x coordinate labels columns x = col * separation # y coordinate labels rows (in reverse order) y = (row_max - row) * separation site = (x, y) register.add(site)

To visually confirm the correctness of the register, we plot its atomic positions with the following code:

import numpy as np import matplotlib.pyplot as plt plt.figure(figsize=(7,7)) plt.plot( register.coordinate_list(0), register.coordinate_list(1), '.', ms=10 ) plt.xticks(np.linspace(0, col_max * separation, col_max + 1)) plt.yticks(np.linspace(0, row_max * separation, row_max + 1)) plt.grid(color='gray', alpha=0.5, lw=1, ls=':') plt.show()

Next, we define the time-dependent parameters of the Hamiltonian, the amplitude, detuning, and phase of the driving field, to realize a simple adiabatic transition sequence.

from braket.timings.time_series import TimeSeries from braket.ahs.driving_field import DrivingField amplitude_min = 0  # rad / s amplitude_max = 15800000   # rad / s detuning_min = -125000000  # rad / s detuning_max = 125000000   # rad / s time_max = 4e-6  # s time_ramp = 1e-7  # s time_points = [0, time_ramp, time_max - time_ramp, time_max] amplitude_values = [amplitude_min, amplitude_max, amplitude_max, amplitude_min] detuning_values = [detuning_min, detuning_min, detuning_max, detuning_max] phase_values = [0, 0, 0, 0] phase = TimeSeries() amplitude = TimeSeries() for t, v in zip(time_points, amplitude_values): amplitude.put(t, v) detuning = TimeSeries() for t, v in zip(time_points, detuning_values): detuning.put(t, v) for t, v in zip(time_points, phase_values): phase.put(t, v) hamiltonian = DrivingField( amplitude=amplitude, detuning=detuning, phase=phase )

Finally, we combine the atom array object and the Hamiltonian into an analog Hamiltonian simulation program and run it on QuEra’s Aquila QPU for 100 shots.

from braket.ahs.analog_hamiltonian_simulation import ( AnalogHamiltonianSimulation ) from braket.aws import AwsDevice aquila_qpu = AwsDevice( "arn:aws:braket:us-east-1::device/qpu/quera/Aquila" ) ahs_program = AnalogHamiltonianSimulation( register=register, hamiltonian=Hamiltonian ).discretize(aquila_qpu) task = aquila_qpu.run(ahs_program, shots=300)

Once the status of the task is COMPLETED (which can be checked with task.metadata()["status"], or on the tasks page of the Braket console), we can query the results withresult = task.result()

The result.measurements object contains a list of shot results. The list, measurements[shot_idx].pre_sequence may contain “0”s which mark sites where atom preparation failed (hence unintentionally missing atoms), and “1”s everywhere else. The measurements[shot_idx].post_sequence measurement results indicate the states of the atoms (“0” for Rydberg state, “1” for ground state) in each of the shot_idx in range(300) shots, where the sequence of atoms is identical to the sequence in which they were added to register, earlier. Shot results containing the most Rydberg stats are the QPU’s best estimates for solving the MIS problem, as they represent the Rydberg atoms, which correspond to the red nodes.

Note, we only run one half iteration of the VQAA optimization algorithm, as we did not adjust the time series parameters of the AHS program. To go beyond this, you can use Amazon Braket Hybrid jobs feature to evolve the system iteratively to minimize the cost function and find the optimal MIS.


The work by Ebadi et al. (Ebadi et al., Science (2022)) realizes a quantum variational optimization algorithm that achieves superlinear improvement compared to classically simulated annealing. Despite the exciting results, many open questions remain. It is unclear for even larger graphs how well VQAA can converge and if the superlinear scaling improvement is retained. Further investigations are needed to compare VQAA’s performance with other state-of-the-art classical approaches, and study how to find the MIS for non-unit-disc graphs on Rydberg systems [Nguyen et al., (2022)].

Our goal is to enable researchers in industry and academia to investigate these questions. With the launch of QuEra’s Aquila QPU to Amazon Braket, a global community of researchers is now able to access state-of-the-art Rydberg computing technology on demand and over the cloud, to advance the understanding of quantum algorithms.

To get started with using the QuEra QPU on Amazon Braket, refer to our AHS example notebooks. For more information about pricing, refer to our Pricing Page. Researchers from accredited institutions can apply for credits by submitting a research proposal here.

Read Entire Article