Experimental demonstration of perturbative anticrossing mitigation using non-uniform driver Hamiltonians

Experimental demonstration of perturbative anticrossing mitigation using non-uniform driver Hamiltonians

Trevor Lanting tlanting@dwavesys.com D-Wave Systems Inc., Burnaby B.C.    Andrew D. King aking@dwavesys.com D-Wave Systems Inc., Burnaby B.C.    Bram Evert bevert@dwavesys.com D-Wave Systems Inc., Burnaby B.C.    Emile Hoskinson ehoskinson@dwavesys.com D-Wave Systems Inc., Burnaby B.C.
August 3, 2019

Perturbative anticrossings have long been identified as a potential computational bottleneck for quantum annealing. This bottleneck can appear, for example, when a uniform transverse driver Hamiltonian is applied to each qubit. Previous theoretical research sought to alleviate such anticrossings by adjusting the transverse driver Hamiltonians on individual qubits according to a perturbative approximation. Here we apply this principle to a physical implementation of quantum annealing in a D-Wave 2000Q system. We use samples from the quantum annealing hardware and per-qubit anneal offsets to produce nonuniform driver Hamiltonians. On small instances with severe perturbative anticrossings, our algorithm yields an increase in minimum eigengaps, ground state success probabilities, and escape rates from metastable valleys. We also demonstrate that the same approach can mitigate biased sampling of degenerate ground states.

I Introduction

Quantum annealing (QA) holds the potential to confer a computational advantage over classical methods via multiqubit tunneling and entanglement Lanting et al. (2014); Denchev et al. (2016). In the near term it offers a scalable alternative to circuit-model quantum computing Andriyash and Amin (2017); Denchev et al. (2016). QA is based on physical evolution of a time-dependent quantum Hamiltonian towards a classical Hamiltonian; the fact that the mid-anneal Hamiltonian may not resemble the final Hamiltonian provides both an advantage and a challenge. In some cases, QA fails with high probability when the instantaneous quantum ground state has overwhelming support from excited states in the computational basis. As the anneal progresses this support can vanish at a small-gap perturbative anticrossing, where QA fails via Landau-Zener transition.

Previously available quantum annealing systems have allowed only uniform QA, in which each qubit is initialized with the same transverse-field driver Hamiltonian and all qubits are annealed in unison. In this case small-gap perturbative anticrossings can arise in unfavorably structured inputs, even those that are not particularly hard for classical solvers Boixo et al. (2013); Steiger et al. (2015); King et al. (2016).

Just as this computational bottleneck can be identified via a perturbative expansion, it can be mitigated with perturbative expansion. Such approaches, in particular those that modify the transverse field on a per-qubit basis, have been proposed and simulated Dickson and Amin (2012). Here we experimentally demonstrate a similar approach in which qubit dynamics are tuned via anneal offsets, which allow for suppression or enhancement of a qubit’s dynamics by annealing the qubit slightly in advance of or behind other qubits, respectively. We find that both the longitudinal and the transverse perturbative corrections contribute to improved performance across a set of tailored inputs. Measurements of tunneling dynamics and spectral analysis support these findings.

Ii Mitigating perturbative anticrossings with non-uniform driver Hamiltonians

We consider a system of qubits connected to one another with tunable longitudinal spin-spin interactions . We can write the uniform QA Hamiltonian as a time-dependent linear combination of an initial driver Hamiltonian and a final problem Hamiltonian


where is the transverse driver energy, is the problem energy, is a normalized annealing parameter, are Pauli matrices operating on the th qubit. At the beginning of the QA algorithm, and at the end of the annealing algorithm . Fig. 1 shows the annealing schedule for the system used herein. All instances studied have for each qubit.

Figure 1: and versus for the QA system under the uniform QA algorithm. A global time-dependent annealing bias tunes both and simultaneously for all qubits throughout the QA algorithm.

ii.1 Generalizing the uniform QA algorithm

The terms and represent the transverse and longitudinal energy scales, which in the uniform QA algorithm do not differ from qubit to qubit. Allowing non-uniform energy scales and on qubit , we can rewrite :


ii.2 Connecting perturbative anticrossings and degeneracy

A spin in a classical state is said to be free or floppy if flipping the spin results in another classical state that is degenerate. Floppy qubits, and associated degeneracies, have been implicated in perturbative anticrossings King et al. (2016); Boixo et al. (2014, 2013). Here we provide an argument, based on first-order perturbation calculations, to connect observed probabilities of qubit floppiness to the change in energy of degenerate excited states when a finite transverse field is present. This analysis, from which we derive our mitigation strategy, closely mirrors the approach of Dickson and Amin in Ref. Dickson and Amin (2012), where the authors use a second-order perturbative expansion because the first order is zero in the problem studied.

We first consider a system with qubits with two-fold degeneracy in the second classical energy level. We can label these excited classical states and , each having energy . We also suppose that and differ by a single bit flip. At , a small transverse field is present which lifts the degeneracy between these states. The new eigenenergies are given by


where are the two eigenvalues of the matrix


In this simple example, and . To first order, the presence of lifts the degeneracy, lowering the energy of the symmetric superposition of and by . Since grows as , this produces a perturbative anticrossing between the ground state and the first excited state.

We now consider a system with qubits and degenerate excited states with energy . When these states are connected via single bit flips, the presence of a transverse field again lifts the degenerate of these states. To first order, the lowest energy




where is the set of pairs of states connected by a single bit flip (see Appendix B and Amin and Choi (2009)). Defining as the qubit index associated with the single bit flip difference, we can write




where denotes the Dirac delta function. We define a per-qubit floppiness metric


which can be interpreted as the fraction of the states in which the th qubit is floppy. To first order, then, there is a relationship between qubit floppiness and the reduction in excited state energy:


Again, grows as , producing a perturbative anticrossing between the ground state and the first excited state.

This suggests a simple mitigation algorithm based on producing a non-uniform per-qubit driver Hamiltonian : If an initial set of samples obtained via the “uniform” annealing algorithm is dominated by states from a highly degenerate cluster of first excited states, we can reduce by decreasing for qubits for which is high relative to other qubits, based on the empirical samples returned by the hardware. If there are multiple clusters of first excited states, this algorithm based on should also lift the lowest energy of each of these clusters. More generally, the goal is to suppress the appearance of states that are disproportionately favored by a perturbative transverse field, and simply mitigating based on these observed states will suppress oversampled states disproportionately.

Before describing our algorithmic implementation of this approach, we need to lay out the limitations of our in situ per-qubit anneal control.

ii.3 Per-qubit anneal offsets

Currently available QA systems do not facilitate independent control of and . Instead we control and together on a per-qubit basis using anneal offsets.

The D-Wave 2000Q system by default implements Hamiltonian (1) with networks of coupled rf-SQUID flux qubits Lanting et al. (2014). See Appendix A for more details. The QA algorithm is run by adjusting a global time-dependent bias, , simultaneously for all of the qubits. We define a normalized global bias, , as


where and for the particular processor used in this study. Like the annealing parameter , ranges from at the beginning of the anneal to at the end of the anneal. Increasing over the course of the anneal both decreases and increases for all qubits; the ratio is a fixed function of determined by the macroscopic rf-SQUID flux qubit device parameters (see Appendix A). Fig. 1 shows and for the experiments discussed herein, and Fig. 2 shows . Note that we choose this particular global trajectory to produce a quadratic growth of .

Figure 2: (Top) relationship between and the tunable parameber . (Bottom) and versus for the QA hardware and for , , and .

We can modify the annealing trajectory for an individual qubit through the use of a tunable in situ static bias that adds or subtracts from the global time-dependent annealing bias signal , giving each qubit a bias signal . We then define as and as . Fig. 2 shows and for . On current hardware, the tunable static bias allows an adjustment of up to . This allows one to advance or delay the annealing signal locally for each qubit.

ii.4 Algorithmic mitigation approach

Following the perturbative analysis in Section II.2, we employ an iterative method similar to the one described in Ref. Dickson and Amin (2012). At each iteration , per-qubit anneal offsets are adjusted to slow the dynamics of qubits that are frequently floppy in the observed states. Adjustments are scaled down for each successive iteration as the iterative search is refined by combining the new and current offsets in an ansatz ratio of .

  1. Choose a static offset magnitude .

  2. Initialize each anneal offset to zero.

  3. For iteration

    1. Anneal times, saving each result.

    2. Compute , the probability of floppiness, for each qubit based on the saved results of the current annealing run.

    3. Adjust each anneal offset according to :

As a simple demonstration of the method we apply one iteration of the method to the 16-qubit system studied in Ref. Dickson et al. (2013), similar to those studied elsewhere Boixo et al. (2013); Albash et al. (2015), shown in Fig. 3. This instance has a unique ground state and a metastable valley of first excited states in which all eight outer qubits (those coupled to only one other qubit) are floppy, leading to a small-gap anticrossing and low ground state probability when run without mitigation. In the first iteration if mitigation, the outer qubits have and the others have . Applying mitigation with increases both ground state probability and minimum eigengap, while the opposite is true when .

Figure 3: (Top) 16-qubit instance studied in Ref. Dickson et al. (2013); circles represent qubits with fields of value (red), (gray), and (blue), and lines represent FM couplings with value . (Bottom) A single iteration of our mitigation strategy advances the outer qubits relative to the inner qubits. Observed ground state probability (bottom left) increases with mitigation magnitude , as does the minimum eigengap (bottom right).

Iii Improving optimization with anneal offsets

iii.1 Optimization testbed

To study the effect of mitigating perturbative crossings with nonuniform driver Hamiltonians, we need small instances—small enough to exactly diagonalize Hamiltonian (1)—in which the uniform QA algorithm will be drawn into a large valley of metastable states by perturbative anticrossings as described in Section II.2. We begin with a qubit connectivity graph with 24 qubits, with each qubit coupled to four others, and construct many thousands of random Ising instances by randomly assigning each a value of or . Fig. 4 shows an example.

Figure 4: A randomly-generated 24-qubit instance shown on a grid of unit cells in the D-Wave 2000Q qubit topology. Qubits are represented by circles; couplers are represented by line segments. Each nonzero coupler is (red) or (blue). This instance has two classical ground states and , and first excited states.

The median ground state probability observed on these instances using the uniform QA algorithm is over . To identify a set of instances for which uniform QA fails due to perturbative anticrossings, we first discard all instances with more than three ground states or fewer than 50 first excited states. We then run all remaining instances with the uniform QA algorithm, performing anneals and collecting all spin configurations. Our testbed consists of those 100 instances—unique up to graph isomorphism and spin reversal—for which the processor yielded the lowest observed ground state probability. Each instance has exactly two classical ground states forming a symmetric pair: and . Each has large valleys of first excited states connected by floppy qubits.

iii.2 Ground state probability

Fig. 5

Figure 5: Algorithmic mitigation improves QA success probability. Shown are ground state probabilities after zero, one, and four iterations of mitigation. Median increases from 62% to 85% over four iterations.

shows measurements of ground state probability for the 100 testbed instances after zero, one, and four iterations of the algorithm outlined in Section II.4 for all 100 instances. We run with and . The median improves from to over four iterations of mitigation.

We examine the minimum eigengaps calculated for these instances—see Appendix C—and compare to the minimum eigengaps calculated after a single iteration of the mitigation algorithm applying offsets for both and , which we term mitigation and antimitigation respectively. Fig. 6 shows the mitigated and antimitigated eigengaps plotted versus the baseline eigengap. The algorithm systematically increases the minimum eigengap for and systematically decreases the eigengap for .

Figure 6: Algorithmic mitigation increases the spectral gap. Shown are the minimum gap between the instantaneous ground state and third excited state.

iii.3 Approximating orthogonal control of and

Anneal offsets do not offer control of a qubit’s transverse field independent of classical energy scale. Advancing a qubit slows its dynamics by both reducing its transverse field and increasing its Ising energy scale. To determine whether or not the positive effects of mitigation can be explained by variation of classical energy scale alone, we repeat the experiment from Section III.2, now compensating for variation in classical energy scale.

To balance classical energy scales at a particular point , we find a new set of couplings that equalize effective energy scales for each coupler at :

Following (4), applying QA to in the place of gives us


approximating independent control of at . We fix , where and differ the most. Results are shown in Fig. 7 using .

Mitigation continues to provide a systematic performance improvement, but less so than what was achieved without balancing classical energy scales. This indicates that both factors contribute to improved performance in the metric. When running these experiments we rescaled all Ising Hamiltonians by a factor of to keep all couplings within the available interval . The increase in minimum baseline when downscaling the Hamiltonian is consistent with previous work Albash et al. (2015); King et al. (2016).

Figure 7: Application of anneal offsets causes imbalance in longitudinal fields between qubits. We approximate orthogonal transverse field control by compensating for this imbalance in the Ising Hamiltonian sent to the QA system. Once this compensation is applied, mitigation still achieves a systematic improvement in success probability, albeit smaller than the improvement seen in Fig. 5.
Figure 8: Measurements of escape rates out of the first excited state. For each instance, we prepared the first excited state and measured the escape rate out of this state at the annealing parameter that corresponded to the location of the anticrossing. We show data for baseline and mitigated .

iii.4 Metastable valley escape rates

To further illustrate the effect of the mitigaiton algorithm, we study the dynamics of relaxation from an excited state during the QA process with and without mitigation. We prepare the QA hardware in a particular classical state at and run the QA algorithm backwards to an intermediate , where we remain for some dwell time before completing the QA algorithm back to . For a given classical state preparation, we can measure the initial escape rate out of this state at .

This measurement protocol is a generalization of the macroscopic resonant tunneling (MRT) protocol described elsewhere Harris et al. (2008). We measure for excited-state preparations for the uniform QA Hamiltonian and for the mitigated and antimitigated case. For each measurement, we choose an that corresponds to the anticrossing identified by exact diagonalization of Hamitonian 4. Fig. 8 shows the rate of escape from the metastable valley for the 100 instances. Mitigation systematically enhances this rate and antimitigation suppresses it.

Iv Sampling degenerate ground states

The initial testbed of 100 instances was designed to exhibit perturbative anticrossings late in the anneal. Now we consider a second testbed of 100 24-qubit instances with highly degenerate ground states in which some ground states are observed much more frequently than others. In these instances deflection of eigenvalues late in the anneal leads to nonuniform sampling of classical states by QA Matsuda et al. (2009); Mandrà et al. (2017); Zhang et al. (2017)—this deflection does not lead to perturbative anticrossing because the eigenvalues converge at . Again the same mitigation approach—which does not consider whether or not observed states are excited—gives both empirical and spectral improvement.

With the optimization testbed, our performance metric was ground state probability . Here, with the sampling testbed, our performance metric is the expected number of samples required to see every ground state (up to antipodal symmetry) given by the empirical probabilities of each ground state, and normalized on a per-instance basis by the same expected value given perfectly uniform sampling of all ground states. This metric is based on the coupon collector problem, and we denote it Flajolet et al. (1992); King et al. (2017). Fig. 9 shows the results analogous to Fig. 5. As before we see a systematic improvement in performance.

Figure 9: Algorithmic mitigation improves uniformity of QA ground state sampling. Shown are the expected number of samples required to observe every ground state, divided by the expected number for equal sampling so that each instance has minimum value 1. Expectation values for 100 instances are shown after zero, one, and four iterations of mitigation.

To measure spectral bias among classical ground states, we look at the support of the mid-anneal instantaneous ground state in the computational basis, specifically those basis vectors corresponding to classical ground states. Where from an empirical perspective we can view classical ground state probabilities output by QA, from a spectral perspective we can view probabilities of classical ground states in the zero-temperature quantum Boltzmann distribution at some fixed point in the anneal. In Fig. 10 we show this data for , contrasting the choice of in Ref. Zhang et al. (2017); if we increase to the results look qualitatively similar but with larger .

Figure 10: Algorithmic mitigation systematically improves uniformity of zero-temperature quantum Boltzmann distribution over classical ground states at .

V Outstanding questions

This work demonstrates a simple approach to a computational bottleneck in quantum annealing. To move towards a more general demonstration there are some outstanding questions that need to be addressed:

  • Orthogonal control: Currently available QA systems do not allow independent adjustment of and ; there is distortion of the classical Hamiltonian when . This distortion may be helpful or harmful; we attempt to compensate for it in Section III.3 and show that our approximation to transverse-only mitigation is successful for the instances studied. The success of this initial demonstration highlights the need for future QA hardware to allow for at least partially independent adjustment of and .

  • Higher-order expansion: We have used a first-order perturbative expansion to estimate the contribution of a small transverse field; this estimate is zero in the absence of floppy qubits. In other problem classes such as Maximum Stable Set Dickson and Amin (2012), disordered Ising spin glasses Katzgraber et al. (2015), or binary Ising spin glasses on graphs with odd connectivity Matsuda et al. (2009); King et al. (2016), higher-order perturbative expansion is required because the first order term is zero. The corresponding second-order algorithm has been shown to work in simulation Dickson and Amin (2012) but remains to be studied in greater detail.

  • Multiple avoided crossings: The instances studied in this work have a relatively simple spectrum with a single avoided crossing. The mitigation algorithm effectively lifts this avoided crossing. In more complicated instances there may be multiple avoided crossings at different points during the QA algorithm Knysh (2016). A single iteration of the mitigation algorithm may alleviate multiple avoided crossings but if the number of avoided crossings increases quickly with problem size or first-order perturbation is insufficient to resolve all (or any) bottlenecks, the algorithm may fail to boost ground state probability.

  • Offset granularity in large systems: For this work we focused on small instances with clear avoided crossings and eigenspectra that are tractable to calculate throughout the QA algorithm. At this problem scale the granularity of available anneal offsets— values are quantized at steps of approximately —is not a barrier to success. Larger systems, in addition to facing the attendant issues mentioned already, may require finer control of anneal offsets.

  • Computational advantage: This mitigation technique must ultimately confer a computational advantage in an algorithmic application in order to motivate practical use. Our choice of 24-qubit instances allowed us to study the eigenspectra of the input instances, but strong evidence for a computational advantage would require the study of larger and more difficult inputs.

Vi Conclusions

We have experimentally demonstrated an algorithm for mitigating perturbative anticrossings in quantum annealing. The algorithm uses anneal offsets to slow the dynamics of floppy qubits, reducing the associated splitting of energy levels that causes perturbative anticrossings. This results in an improvement of optimization performance in 100 small spin-glass instances designed to exhibit a response to this first-order perturbative mitigation. Exact calculations of eigenspectra throughout the anneal confirm that this algorithm increases the eigengap and antimitigating reduces the eigengap. Dynamics measurements confirm that under mitigation, relaxation out of excited states is enhanced; under antimitigation the opposite holds.

This work is the first experimental demonstration that a key computational bottleneck for quantum annealing—perturbative anticrossings—can be removed or mitigated with a targeted change in driving transverse Hamiltonian. Extending these ideas to more complex input instances and more sophisticated mitigation algorithms will be the subject of future research.

Vii Acknowledgments

We thank Richard Harris, Evgeny Andriyash, Mohammad Amin, Jack Raymond, and others at D-Wave for helpful conversations on this work and related topics.


Appendix A QA apparatus and methods

We acquired the data reported in this manuscript with a D-Wave 2000Q system in Burnaby, British Columbia. Each qubit is a compound-compound Josephson-junction rf SQUID Harris et al. (2010). The total critical current of the rf-SQUID junctions is A, the rf-SQUID body inductance is pH, and the rf-SQUID capacitance is fF. The qubits are tunably coupled with inter-qubit magnetic coupling elements that produce a maximum anti-ferromagnetic, pH Harris et al. (2009).

We run the QA algorithm by adjusting the external bias on the compound-compound Josephson-junction loop from at to at according to the curve shown in Fig. 1. We ran the QA algorithm with at a processor temperature of . Of the qubits physically present, were operational. This allowed us to embed disjoint copies of a -qubit system on the processor, each subject to an independent spin reversal transformation King and McGeoch (2014).

Figure 11: First, second, and third eigengaps of three exemplary instances using antimitigation (left), no mitigation (middle), and mitigation (right). First eigengap always goes to zero late in the anneal; the third eigengap gives a clean representation of the desired tunneling event between superpositions of classical ground and excited states.

Appendix B Perturbation calculation details

Consider a system with degenerate eigenstates with energy connected to one another via single bit flips. The presence of a transverse field lifts this degeneracy by an amount . To first order in degenerate perturbation theory, is the smallest eigenvalue of the matrix where and denote specific degenerate eigenstates in the set of .

is the negative adjacency matrix for states in this excited state manifold. For simplicity, we first assume that each state is connected to the same number of other states via single bit flips, i.e. the connectivity of these states is a regular graph. It is trivial, then, to check that the element vector


is an eigenvector of with eigenvalue


where is the set of pairs of states connected by a single bit flip.

Because all elements , this eigenvalue corresponds to the minimum eigenvalue Perron (1907). Thus, and


For more complicated connection topologies of degenerate first excited states, that is, non-regular graphs, Eq. 17 is still correct to first order and a good approximation as long as the connectivity graph is approximately regular.

Appendix C Using third eigengap

Spectral analysis of quantum annealing typically considers the first eigengap Farhi et al. (2000, 2001); Smelyanskiy et al. (2001), since this value governs guarantees of adiabaticity in a closed system. Since the 100 instances considered here have two antipodal classical ground states, this gap goes to zero at the end of the anneal, so analyzing its minimum is meaningless. Instead we consider the third eigengap . The minimum third eigengap gives a consistent representation of the tunneling between ground-state and excited-state manifolds when the instantaneous ground state transitions from being mainly supported by classical excited states to being mainly supported by classical ground states. Fig. 11 shows three exemplary instances under antimitigation, no mitigation, and mitigation, where the first and second gaps do not necessarily represent the desired tunneling event.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description