Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2605.19381

A calibrated diagnostic for reverse-anneal sampling in programmable quantum annealers

Luis Lozano

quant-phcond-mat.dis-nnSubmitted May 19, 2026• Announced May 20, 2026• 0 citations• Revised Aug 9, 2026 (v2)
v1 · May 19, 2026v2 · Aug 9, 2026
arXivHTMLPDF

Abstract

Reverse annealing is widely used as a heuristic sampler on programmable quantum annealers, but it is often unclear whether the readout has erased its initialization memory or whether the resulting samples represent any calibrated target distribution. We introduce a subsystem-level validation protocol that pairs a memory order parameter $\mathcal{M}$ with the total-variation distance $D_{\mathrm{TV}}$ between the measured subsystem readout and a fixed, independently calibrated conditional-Boltzmann reference, and apply it on two D-Wave QPU generations. The protocol delivers three things: a diagnostic that separates relaxed, memory-retaining and "wrong-basin" readouts; the isolation of relaxed-but-non-thermal trapping, in which a readout passes the initial-state-independence test while concentrating on the wrong basin; and an engineering control showing that SDK auto-scale, embedding choice and environment preparation shift apparent thresholds, so calibration metadata must be logged to interpret them. Relaxed ferromagnetic readouts are near-deterministic, so the small distances seen there are a consistency check rather than a thermometric measurement; the diagnostic earns its value in the opposite regime. Across 494 conditions on one processor, 113 pass the memory criterion and 110 of those agree with the reference at $D_{\mathrm{TV}}<0.05$, while a wrong-basin instance, reproduced at two Hamiltonian scales, passes the memory test yet reaches $D_{\mathrm{TV}}>0.94$; a 60-condition replication confirms that the diagnostic transfers to a second processor. We provide the protocol as a transferable, reproducible validation workflow for annealer-as-sampler applications.

Comments
v2: substantially revised: title and framing updated, benchmarking section removed following a reproducibility audit, claims narrowed; reformatted to revtex4-2 (11 pages, 5 figures + supplemental material as ancillary file). All remaining results regenerate from the deposited raw data

Related papers

Schedule-dependent basin occupation and instance-specific effective scales in cycled reverse annealing on programmable quantum annealers

Luis Lozano

Reverse annealing initializes a programmable quantum annealer in a chosen classical state, reintroduces quantum fluctuations down to a pause point, and returns to readout; repeated cycles form a configurable dissipative dynamics whose schedule knobs are routinely treated as a temperature. On mixed-frustration 12-qubit Ising instances on two D-Wave generations (Advantage2 and Advantage_system6.4), we measure schedule-driven redistribution of the subsystem readout and calibrate each schedule's effective scales, using exhaustive enumeration of the $2^{12}$ configurations, a locked analysis plan with a prior-knowledge disclosure, and regenerated classical reference families. First, schedule shape redistributes basin occupation: 11 of 20 random instances respond to a single-pause-to-two-pause change (15 of 20 for long-hold; block permutation, FDR-controlled), with shifts up to 38 percentage points and changes of dominant-configuration identity. Second, the single-qubit probe temperature does not transfer to the many-body problem: moment-matched inverse temperatures span 0.63-9.1 against a probe value of 7.22, and pause depth and duration each shift the scale on most instances (19/20 and 15/20). Third, no observation requires more than the tested classical families: with update budgets extended below one sweep per cycle, all 33 (device, instance, schedule) cells are reproduced within the tested family (minimum intersection-union $p=0.76$), and schedule-aware spin-vector Monte Carlo brings the late-window statistics within semiclassical reach. A pre-registered linear predictor fails on held-out instances, and an earlier pilot's two-pause enhancement claim is contradicted (Fisher $p=0.0013$). Reverse-anneal schedules are instance-specific basin-occupation controls whose effective scales must be calibrated per instance and schedule; nothing here licenses a quantum-mechanistic reading.

quant-pharXiv:2605.17604•May 19, 2026•0 citations

Quantum annealers as programmable thermal machines

Jakub Pawłowski, Tomasz Śmierzchalski, Fengping Jin +3 more

Programmable quantum annealers are used for optimization, probabilistic sampling, and simulation, but their performance is commonly reported without the energy exchanged during computation. Here we characterize the D-Wave quantum annealer as a closed thermodynamic cycle. From initial and final Ising energies and an effective temperature fitted to the output distribution, we obtain lower bounds on entropy production, environment energy exchange, work, and power. By varying the prepared distribution and the reverse annealing turning point, we map heater-, accelerator-, refrigerator-, and engine-compatible regimes in one dimensional chains and higher connectivity instances, and apply the same analysis to Advantage and Advantage2 hardware. For an encoded optimization problem, the measured processor energy change states whether final candidates improve or worsen the programmed objective on average. For sampling, the fitted temperature provides an operational measure of how strongly probability is concentrated among low energy configurations. The thermodynamic mode therefore adds information absent from solution quality or runtime alone: it distinguishes driven refinement, net heating, and heat pumping while quantifying their energetic consequences. This framework connects quantum optimization, probabilistic computing, statistical physics simulation, hardware diagnostics, and energy-aware assessment without assuming that a thermodynamic label alone determines computational performance.

quant-pharXiv:2608.04564•Aug 6, 2026•0 citations

Partition Function Estimation Using Analog Quantum Processors

Thinh Le, Elijah Pelofske

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 $\pm 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 \times 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.

quant-pharXiv:2512.19685•Dec 23, 2025•0 citations

Testing a Quantum Annealer as a Quantum Thermal Sampler

Zoe Gonzalez Izquierdo, Itay Hen, Tameem Albash

Motivated by recent experiments in which specific thermal properties of complex many-body systems were successfully reproduced on a commercially available quantum annealer, we examine the extent to which quantum annealing hardware can reliably sample from the thermal state in a specific basis associated with a target quantum Hamiltonian. We address this question by studying the diagonal thermal properties of the canonical one-dimensional transverse-field Ising model on a D-Wave 2000Q quantum annealing processor. We find that the quantum processor fails to produce the correct expectation values predicted by Quantum Monte Carlo. Comparing to master equation simulations, we find that this discrepancy is best explained by how the measurements at finite transverse fields are enacted on the device. Specifically, measurements at finite transverse field require the system to be quenched from the target Hamiltonian to a Hamiltonian with negligible transverse field, and this quench is too slow. The limitations imposed by such hardware make it an unlikely candidate for thermal sampling, and it remains an open question what thermal expectation values can be robustly estimated in general for arbitrary quantum many-body systems.

arXiv:2003.00361•ACM Transactions on Quantum Computing•Jun 30, 2021•18 citations

Analysis of a Programmable Quantum Annealer as a Random Number Generator

Elijah Pelofske

Quantum devices offer a highly useful function - that is generating random numbers in a non-deterministic way since the measurement of a quantum state is not deterministic. This means that quantum devices can be constructed that generate qubits in a uniform superposition and then measure the state of those qubits. If the preparation of the qubits in a uniform superposition is unbiased, then quantum computers can be used to create high entropy, secure random numbers. Typically, preparing and measuring such quantum systems requires more time compared to classical pseudo random number generators (PRNGs) which are inherently deterministic algorithms. Therefore, the typical use of quantum random number generators (QRNGs) is to provide high entropy secure seeds for PRNGs. Quantum annealing (QA) is a type of analog quantum computation that is a relaxed form of adiabatic quantum computation and uses quantum fluctuations in order to search for ground state solutions of a programmable Ising model. Here we present extensive experimental random number results from a D-Wave 2000Q quantum annealer, totaling over 20 billion bits of QA measurements, which is significantly larger than previous D-Wave QA random number generator studies. Current quantum annealers are susceptible to noise from environmental sources and calibration errors, and are not in general unbiased samplers. Therefore, it is of interest to quantify whether noisy quantum annealers can effectively function as an unbiased QRNG. The amount of data that was collected from the quantum annealer allows a comprehensive analysis of the random bits to be performed using the NIST SP 800-22 Rev 1a testsuite, as well as min-entropy estimates from NIST SP 800-90B. The randomness tests show that the generated random bits from the D-Wave 2000Q are biased, and not unpredictable random bit sequences. With no server-side sampling post-processing, the 1 microsecond annealing time measurements had a min-entropy of 0.824.

arXiv:2307.02573•IEEE Transactions on Information Forensics and Security•Jan 1, 2024•4 citations

Optimization versus fair sampling in quantum annealing: An exactly characterized frustrated Ising benchmark

Masaru Hitomi, Masayuki Ohzeki

Reaching a degenerate ground-state manifold and sampling it uniformly are distinct properties of a quantum annealer. We quantify this distinction using an Ising ring with additional couplings between maximally separated sites. For antiferromagnetic benchmark couplings and $N \geq 5$, we derive the exact ground-state energy, degeneracy, and defect structure for each congruence class of $N$ modulo four. This classification separates ground-state hit probability from conditional sampling bias. We embed instances with up to 162 logical spins on a D-Wave Advantage processor using length-two chains and compare its output with simulated annealing on the embedded and logical problems. For a frustrated $N=83$ instance, the processor visits 164 of the 166 exact ground states, yet the conditional distributions over states and defect sectors remain nonuniform. For $N=160$, the processor has a larger ground-state hit probability than logical simulated annealing (0.869 versus 0.519), but also a larger total-variation distance from uniformity (0.317 versus 0.136). A pause during the anneal raises the ground-state hit probability from 0.42 to 0.86 without restoring uniform sampling. These results identify a relaxation window and demonstrate that optimization success and sampling quality need not improve together, providing an exactly characterized benchmark for open-system quantum annealing.

quant-pharXiv:2610.29035•Oct 22, 2026•0 citations
Reverse annealing is widely used as a heuristic sampler on programmable quantum annealers, but it is often unclear whether the readout has erased its initialization memory or whether the resulting samples represent any calibrated target distribution. We introduce a subsystem-level validation protocol that pairs a memory order parameter M\mathcal{M}M with the total-variation distance DTVD_{\mathrm{TV}}DTV​ between the measured subsystem readout and a fixed, independently calibrated conditional-Boltzmann reference, and apply it on two D-Wave QPU generations. The protocol delivers three things: a diagnostic that separates relaxed, memory-retaining and "wrong-basin" readouts; the isolation of relaxed-but-non-thermal trapping, in which a readout passes the initial-state-independence test while concentrating on the wrong basin; and an engineering control showing that SDK auto-scale, embedding choice and environment preparation shift apparent thresholds, so calibration metadata must be logged to interpret them. Relaxed ferromagnetic readouts are near-deterministic, so the small distances seen there are a consistency check rather than a thermometric measurement; the diagnostic earns its value in the opposite regime. Across 494 conditions on one processor, 113 pass the memory criterion and 110 of those agree with the reference at DTV<0.05D_{\mathrm{TV}}<0.05DTV​<0.05, while a wrong-basin instance, reproduced at two Hamiltonian scales, passes the memory test yet reaches DTV>0.94D_{\mathrm{TV}}>0.94DTV​>0.94; a 60-condition replication confirms that the diagnostic transfers to a second processor. We provide the protocol as a transferable, reproducible validation workflow for annealer-as-sampler applications.
Reverse annealing initializes a programmable quantum annealer in a chosen classical state, reintroduces quantum fluctuations down to a pause point, and returns to readout; repeated cycles form a configurable dissipative dynamics whose schedule knobs are routinely treated as a temperature. On mixed-frustration 12-qubit Ising instances on two D-Wave generations (Advantage2 and Advantage_system6.4), we measure schedule-driven redistribution of the subsystem readout and calibrate each schedule's effective scales, using exhaustive enumeration of the 2122^{12}212 configurations, a locked analysis plan with a prior-knowledge disclosure, and regenerated classical reference families. First, schedule shape redistributes basin occupation: 11 of 20 random instances respond to a single-pause-to-two-pause change (15 of 20 for long-hold; block permutation, FDR-controlled), with shifts up to 38 percentage points and changes of dominant-configuration identity. Second, the single-qubit probe temperature does not transfer to the many-body problem: moment-matched inverse temperatures span 0.63-9.1 against a probe value of 7.22, and pause depth and duration each shift the scale on most instances (19/20 and 15/20). Third, no observation requires more than the tested classical families: with update budgets extended below one sweep per cycle, all 33 (device, instance, schedule) cells are reproduced within the tested family (minimum intersection-union p=0.76p=0.76p=0.76), and schedule-aware spin-vector Monte Carlo brings the late-window statistics within semiclassical reach. A pre-registered linear predictor fails on held-out instances, and an earlier pilot's two-pause enhancement claim is contradicted (Fisher p=0.0013p=0.0013p=0.0013). Reverse-anneal schedules are instance-specific basin-occupation controls whose effective scales must be calibrated per instance and schedule; nothing here licenses a quantum-mechanistic reading.
Programmable quantum annealers are used for optimization, probabilistic sampling, and simulation, but their performance is commonly reported without the energy exchanged during computation. Here we characterize the D-Wave quantum annealer as a closed thermodynamic cycle. From initial and final Ising energies and an effective temperature fitted to the output distribution, we obtain lower bounds on entropy production, environment energy exchange, work, and power. By varying the prepared distribution and the reverse annealing turning point, we map heater-, accelerator-, refrigerator-, and engine-compatible regimes in one dimensional chains and higher connectivity instances, and apply the same analysis to Advantage and Advantage2 hardware. For an encoded optimization problem, the measured processor energy change states whether final candidates improve or worsen the programmed objective on average. For sampling, the fitted temperature provides an operational measure of how strongly probability is concentrated among low energy configurations. The thermodynamic mode therefore adds information absent from solution quality or runtime alone: it distinguishes driven refinement, net heating, and heat pumping while quantifying their energetic consequences. This framework connects quantum optimization, probabilistic computing, statistical physics simulation, hardware diagnostics, and energy-aware assessment without assuming that a thermodynamic label alone determines computational performance.
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 252525 spin ±J\pm J±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−67.6 \times 10^{-6}7.6×10−6, from 171,000171,000171,000 samples generated using 0.20.20.2 seconds of QPU time with an anneal time of 888 nanoseconds per sample which is interestingly within the closed system dynamics timescale of the superconducting qubits.
Motivated by recent experiments in which specific thermal properties of complex many-body systems were successfully reproduced on a commercially available quantum annealer, we examine the extent to which quantum annealing hardware can reliably sample from the thermal state in a specific basis associated with a target quantum Hamiltonian. We address this question by studying the diagonal thermal properties of the canonical one-dimensional transverse-field Ising model on a D-Wave 2000Q quantum annealing processor. We find that the quantum processor fails to produce the correct expectation values predicted by Quantum Monte Carlo. Comparing to master equation simulations, we find that this discrepancy is best explained by how the measurements at finite transverse fields are enacted on the device. Specifically, measurements at finite transverse field require the system to be quenched from the target Hamiltonian to a Hamiltonian with negligible transverse field, and this quench is too slow. The limitations imposed by such hardware make it an unlikely candidate for thermal sampling, and it remains an open question what thermal expectation values can be robustly estimated in general for arbitrary quantum many-body systems.
Quantum devices offer a highly useful function - that is generating random numbers in a non-deterministic way since the measurement of a quantum state is not deterministic. This means that quantum devices can be constructed that generate qubits in a uniform superposition and then measure the state of those qubits. If the preparation of the qubits in a uniform superposition is unbiased, then quantum computers can be used to create high entropy, secure random numbers. Typically, preparing and measuring such quantum systems requires more time compared to classical pseudo random number generators (PRNGs) which are inherently deterministic algorithms. Therefore, the typical use of quantum random number generators (QRNGs) is to provide high entropy secure seeds for PRNGs. Quantum annealing (QA) is a type of analog quantum computation that is a relaxed form of adiabatic quantum computation and uses quantum fluctuations in order to search for ground state solutions of a programmable Ising model. Here we present extensive experimental random number results from a D-Wave 2000Q quantum annealer, totaling over 20 billion bits of QA measurements, which is significantly larger than previous D-Wave QA random number generator studies. Current quantum annealers are susceptible to noise from environmental sources and calibration errors, and are not in general unbiased samplers. Therefore, it is of interest to quantify whether noisy quantum annealers can effectively function as an unbiased QRNG. The amount of data that was collected from the quantum annealer allows a comprehensive analysis of the random bits to be performed using the NIST SP 800-22 Rev 1a testsuite, as well as min-entropy estimates from NIST SP 800-90B. The randomness tests show that the generated random bits from the D-Wave 2000Q are biased, and not unpredictable random bit sequences. With no server-side sampling post-processing, the 1 microsecond annealing time measurements had a min-entropy of 0.824.
Reaching a degenerate ground-state manifold and sampling it uniformly are distinct properties of a quantum annealer. We quantify this distinction using an Ising ring with additional couplings between maximally separated sites. For antiferromagnetic benchmark couplings and N≥5N \geq 5N≥5, we derive the exact ground-state energy, degeneracy, and defect structure for each congruence class of NNN modulo four. This classification separates ground-state hit probability from conditional sampling bias. We embed instances with up to 162 logical spins on a D-Wave Advantage processor using length-two chains and compare its output with simulated annealing on the embedded and logical problems. For a frustrated N=83N=83N=83 instance, the processor visits 164 of the 166 exact ground states, yet the conditional distributions over states and defect sectors remain nonuniform. For N=160N=160N=160, the processor has a larger ground-state hit probability than logical simulated annealing (0.869 versus 0.519), but also a larger total-variation distance from uniformity (0.317 versus 0.136). A pause during the anneal raises the ground-state hit probability from 0.42 to 0.86 without restoring uniform sampling. These results identify a relaxation window and demonstrate that optimization success and sampling quality need not improve together, providing an exactly characterized benchmark for open-system quantum annealing.