Unraveling Quantum Annealers using Classical Hardness

Recent advances in quantum technology have led to the development and manufacturing of experimental programmable quantum annealing optimizers that contain hundreds of quantum bits. These optimizers, named `D-Wave' chips, promise to solve practical optimization problems potentially faster than conventional `classical' computers. Attempts to quantify the quantum nature of these chips have been met with both excitement and skepticism but have also brought up numerous fundamental questions pertaining to the distinguishability of quantum annealers from their classical thermal counterparts. Here, we propose a general method aimed at answering these, and apply it to experimentally study the D-Wave chip. Inspired by spin-glass theory, we generate optimization problems with a wide spectrum of `classical hardness', which we also define. By investigating the chip's response to classical hardness, we surprisingly find that the chip's performance scales unfavorably as compared to several analogous classical algorithms. We detect, quantify and discuss purely classical effects that possibly mask the quantum behavior of the chip.Comment: 12 pages, 9 figure

Testing statics-dynamics equivalence at the spin-glass transition in three dimensions

The statics-dynamics correspondence in spin glasses relate non-equilibrium results on large samples (the experimental realm) with equilibrium quantities computed on small systems (the typical arena for theoretical computations). Here we employ statics-dynamics equivalence to study the Ising spin-glass critical behavior in three dimensions. By means of Monte Carlo simulation, we follow the growth of the coherence length (the size of the glassy domains), on lattices too large to be thermalized. Thanks to the large coherence lengths we reach, we are able to obtain accurate results in excellent agreement with the best available equilibrium computations. To do so, we need to clarify the several physical meanings of the dynamic exponent close to the critical temperature.Comment: Version to appear in Physical Review

An Ising Model for Metal-Organic Frameworks

We present a three-dimensional Ising model where lines of equal spins are frozen in such that they form an ordered framework structure. The frame spins impose an external field on the rest of the spins (active spins). We demonstrate that this "porous Ising model" can be seen as a minimal model for condensation transitions of gas molecules in metal-organic frameworks. Using Monte Carlo simulation techniques, we compare the phase behavior of a porous Ising model with that of a particle-based model for the condensation of methane (CH$_4$) in the isoreticular metal-organic framework IRMOF-16. For both models, we find a line of first-order phase transitions that end in a critical point. We show that the critical behavior in both cases belongs to the 3D Ising universality class, in contrast to other phase transitions in confinement such as capillary condensation.Comment: 11 pages, 9 figure

Advantages of Unfair Quantum Ground-State Sampling

The debate around the potential superiority of quantum annealers over their classical counterparts has been ongoing since the inception of the field by Kadowaki and Nishimori close to two decades ago. Recent technological breakthroughs in the field, which have led to the manufacture of experimental prototypes of quantum annealing optimizers with sizes approaching the practical regime, have reignited this discussion. However, the demonstration of quantum annealing speedups remains to this day an elusive albeit coveted goal. Here, we examine the power of quantum annealers to provide a different type of quantum enhancement of practical relevance, namely, their ability to serve as useful samplers from the ground-state manifolds of combinatorial optimization problems. We study, both numerically by simulating ideal stoquastic and non-stoquastic quantum annealing processes, and experimentally, using a commercially available quantum annealing processor, the ability of quantum annealers to sample the ground-states of spin glasses differently than classical thermal samplers. We demonstrate that i) quantum annealers in general sample the ground-state manifolds of spin glasses very differently than thermal optimizers, ii) the nature of the quantum fluctuations driving the annealing process has a decisive effect on the final distribution over ground-states, and iii) the experimental quantum annealer samples ground-state manifolds significantly differently than thermal and ideal quantum annealers. We illustrate how quantum annealers may serve as powerful tools when complementing standard sampling algorithms.Comment: 13 pages, 11 figure

Precursors of the Spin Glass Transition in Three Dimensions

We study energy landscape and dynamics of the three-dimensional Heisenberg Spin Glass model in the paramagnetic phase, i.e. for temperature $T$ larger than the critical temperature $T_\mathrm{c}$. The landscape is non-trivially related to the equilibrium states even in the high-temperature phase, and reveals an onset of non-trivial behavior at a temperature $T_\mathrm{o}$, which is also seen through the behavior of the thermoremanent magnetization. We also find a power-law growth of the relaxation times far from the spin-glass transition, indicating a dynamical crossover at a temperature $T_\mathrm{d}$, $T_\mathrm{c}<T_\mathrm{d}<T_\mathrm{o}$. The arising picture is reminiscent of the phenomenology of supercooled liquids, and poses questions on which mean-field models can describe qualitatively well the phenomenology in three dimensions. On the technical side, local energy minima are found with the Successive Overrelaxation algorithm, which reveals very efficient for energy minimization in this kind of models.Comment: 16 pages, 6 figure

Analog Errors in Ising Machines

Recent technological breakthroughs have precipitated the availability of specialized devices that promise to solve NP-Hard problems faster than standard computers. These `Ising Machines' are however analog in nature and as such inevitably have implementation errors. We find that their success probability decays exponentially with problem size for a fixed error level, and we derive a sufficient scaling law for the error in order to maintain a fixed success probability. We corroborate our results with experiment and numerical simulations and discuss the practical implications of our findings.Comment: 14 pages, 15 figures. v2: Updated to published versio

Specific-heat exponent and modified hyperscaling in the 4D random-field Ising model

We report a high-precision numerical estimation of the critical exponent $\alpha$ of the specific heat of the random-field Ising model in four dimensions. Our result $\alpha = 0.12(1)$ indicates a diverging specific-heat behavior and is consistent with the estimation coming from the modified hyperscaling relation using our estimate of $\theta$ via the anomalous dimensions $\eta$ and $\bar{\eta}$. Our analysis benefited form a high-statistics zero-temperature numerical simulation of the model for two distributions of the random fields, namely a Gaussian and Poissonian distribution, as well as recent advances in finite-size scaling and reweighting methods for disordered systems. An original estimate of the critical slowing down exponent $z$ of the maximum-flow algorithm used is also provided.Comment: 11 pages, 3 figures, 1 table. arXiv admin note: text overlap with arXiv:1512.0657

Finite-size scaling of the random-field Ising model above the upper critical dimension

Finite-size scaling above the upper critical dimension is a long-standing puzzle in the field of Statistical Physics. Even for pure systems various scaling theories have been suggested, partially corroborated by numerical simulations. In the present manuscript we address this problem in the even more complicated case of disordered systems. In particular, we investigate the scaling behavior of the random-field Ising model at dimension $D = 7$, i.e., above its upper critical dimension $D_{\rm u} = 6$, by employing extensive ground-state numerical simulations. Our results confirm the hypothesis that at dimensions $D > D_{\rm u}$, linear length scale $L$ should be replaced in finite-size scaling expressions by the effective scale $L_{\rm eff} = L^{D / D_{\rm u}}$. Via a fitted version of the quotients method that takes this modification, but also subleading scaling corrections into account, we compute the critical point of the transition for Gaussian random fields and provide estimates for the full set of critical exponents. Thus, our analysis indicates that this modified version of finite-size scaling is successful also in the context of the random-field problem.Comment: 19 pages preprint style, 5 figures, Appendix include

On the high-density expansion for Euclidean random matrices

Diagrammatic techniques to compute perturbatively the spectral properties of Euclidean random matrices (ERM) in the high-density regime are introduced and discussed in detail. Such techniques are developed in two alternative and very different formulations of the mathematical problem and are shown to give identical results up to second order in the perturbative expansion. One method, based on writing the so-called resolvent function as a Taylor series, allows us to group the diagrams into a small number of topological classes, providing a simple way to determine the infrared (small momenta) behaviour of the theory up to third order, which is of interest for the comparison with experiments. The other method, which reformulates the problem as a field theory, can instead be used to study the infrared behaviour at any perturbative order.Facultad de Ciencias ExactasInstituto de Investigaciones Fisicoquímicas Teóricas y Aplicada