Understanding the physical nature of the D-Wave annealers remains a subject of active investigation. In this study, we analyze the sampling behavior of these systems and explore whether their results can be replicated using quantum and Markovian models. Employing the standard and the fast annealing protocols, we observe that the D-Wave annealers sample states with frequencies matching the Gibbs distribution for sufficiently long annealing times. Using Bloch equation simulations for single-qubit problems and Lindblad and Markovian master equations for two-qubit systems, we compare experimental data with theoretical predictions. Our results provide insights into the role of quantum mechanics in these devices.
Understanding the physical nature of the D-Wave annealers remains a subject of active investigation. In this study, we analyze the sampling behavior of these systems and explore whether their results can be replicated using quantum and Markovian models. Employing the standard and the fast annealing protocols, we observe that the D-Wave annealers sample states with frequencies matching the Gibbs distribution for sufficiently long annealing times. Using Bloch equation simulations for single-qubit problems and Lindblad and Markovian master equations for two-qubit systems, we compare experimental data with theoretical predictions. Our results provide insights into the role of quantum mechanics in these devices.
We evaluate using programmable superconducting flux qubit D-Wave quantum annealers to approximate the partition function of Ising models. We propose the use of two distinct quantum annealer sampling methods: chains of Monte Carlo-like reverse quantum anneals, and standard linear-ramp quantum annealing. The control parameters used to attenuate the quality of the simulations are the effective analog energy scale of the J coupling, the total annealing time, and for the case of reverse annealing the anneal-pause. The core estimation technique is to sample across the energy spectrum of the classical Hamiltonian of interest, and therefore obtain a density of states estimate for each energy level, which in turn can be used to compute an estimate of the partition function with some sampling error. This estimation technique is powerful because once the distribution is sampled it allows thermodynamic quantity computation at arbitrary temperatures. On a
25 spin
±J hardware graph native Ising model we find parameter regimes of the D-Wave processors that provide comparable result quality to two standard classical Monte Carlo methods, Multiple Histogram Reweighting and Wang-Landau. Remarkably, we find that fast quench-like anneals can quickly generate ensemble distributions that are very good estimates of the true partition function of the classical Ising model; on a Pegasus graph-structured QPU we report a logarithmic relative error of
7.6×10−6, from
171,000 samples generated using
0.2 seconds of QPU time with an anneal time of
8 nanoseconds per sample which is interestingly within the closed system dynamics timescale of the superconducting qubits.