BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Helmut G. Katzgraber

Helmut G. Katzgraber appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.dis-nnquant-phcond-mat.stat-mechNeural and Evolutionary Computingphysics.comp-phcond-mat.mes-hallmath.OCArtificial Intelligencecond-mat.mtrl-sciphysics.chem-phcond-mat.str-elData Structures and AlgorithmsMachine LearningMathematical SoftwareNumerical Analysisphysics.data-anphysics.med-phphysics.soc-ph

Catalog footprint

What is connected

66works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2025arXiv

A Random-Key Optimizer for Combinatorial Optimization

This paper introduces the Random-Key Optimizer (RKO), a versatile and efficient stochastic local search method tailored for combinatorial optimization problems. Using the random-key concept, RKO encodes solutions as vectors of random keys that are subsequently decoded into feasible solutions via problem-specific decoders. The RKO framework is able to combine a plethora of classic metaheuristics, each capable of operating independently or in parallel, with solution sharing facilitated through an elite solution pool. This modular approach allows for the adaptation of various metaheuristics, including simulated annealing, iterated local search, and greedy randomized adaptive search procedures, among others. The efficacy of the RKO framework, implemented in C++ and publicly available (Github public repository: github.com/RKO-solver), is demonstrated through its application to three NP-hard combinatorial optimization problems: the alpha-neighborhood p-median problem, the tree of hubs location problem, and the node-capacitated graph partitioning problem. The results highlight the framework's ability to produce high-quality solutions across diverse problem domains, underscoring its potential as a robust tool for combinatorial optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Automated Design of Pulse Sequences for Magnetic Resonance Fingerprinting using Physics-Inspired Optimization

Magnetic Resonance Fingerprinting (MRF) is a method to extract quantitative tissue properties such as T1 and T2 relaxation rates from arbitrary pulse sequences using conventional magnetic resonance imaging hardware. MRF pulse sequences have thousands of tunable parameters which can be chosen to maximize precision and minimize scan time. Here we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggests systematic errors dominate over random errors in MRF scans under clinically-relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in T1 and T2. Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Combinatorial Optimization with Physics-Inspired Graph Neural Networks

Combinatorial optimization problems are pervasive across science and industry. Modern deep learning tools are poised to solve these problems at unprecedented scales, but a unifying framework that incorporates insights from statistical physics is still outstanding. Here we demonstrate how graph neural networks can be used to solve combinatorial optimization problems. Our approach is broadly applicable to canonical NP-hard problems in the form of quadratic unconstrained binary optimization problems, such as maximum cut, minimum vertex cover, maximum independent set, as well as Ising spin glasses and higher-order generalizations thereof in the form of polynomial unconstrained binary optimization problems. We apply a relaxation strategy to the problem Hamiltonian to generate a differentiable loss function with which we train the graph neural network and apply a simple projection to integer variables once the unsupervised training process has completed. We showcase our approach with numerical results for the canonical maximum cut and maximum independent set problems. We find that the graph neural network optimizer performs on par or outperforms existing solvers, with the ability to scale beyond the state of the art to problems with millions of variables.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adding color: Visualization of energy landscapes in spin glasses

Disconnectivity graphs are used to visualize the minima and the lowest energy barriers between the minima of complex systems. They give an easy and intuitive understanding of the underlying energy landscape and, as such, are excellent tools for understanding the complexity involved in finding low-lying or global minima of such systems. We have developed a classification scheme that categorizes highly-degenerate minima of spin glasses based on similarity and accessibility of the individual states. This classification allows us to condense the information pertained in different dales of the energy landscape to a single representation using color to distinguish its type and a bar chart to indicate the average size of the dales at their respective energy levels. We use this classification to visualize disconnectivity graphs of small representations of different tile-planted models of spin glasses. An analysis of the results shows that different models have distinctly different features in the total number of minima, the distribution of the minima with respect to the ground state, the barrier height and in the occurrence of the different types of minimum energy dales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computational hardness of spin-glass problems with tile-planted solutions

We investigate the computational hardness of spin-glass instances on a square lattice, generated via a recently introduced tunable and scalable approach for planting solutions. The method relies on partitioning the problem graph into edge-disjoint subgraphs, and planting frustrated, elementary subproblems that share a common local ground state, which guarantees that the ground state of the entire problem is known a priori. Using population annealing Monte Carlo, we compare the typical hardness of problem classes over a large region of the multi-dimensional tuning parameter space. Our results show that the problems have a wide range of tunable hardness. Moreover, we observe multiple transitions in the hardness phase space, which we further corroborate using simulated annealing and simulated quantum annealing. By investigating thermodynamic properties of these planted systems, we demonstrate that the harder samples undergo magnetic ordering transitions which are also ultimately responsible for the observed hardness transitions on changing the sample composition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Learning to find order in disorder

We introduce the use of neural networks as classifiers on classical disordered systems with no spatial ordering. In this study, we implement a convolutional neural network trained to identify the spin-glass state in the three-dimensional Edwards-Anderson Ising spin-glass model from an input of Monte Carlo sampled configurations at a given temperature. The neural network is designed to be flexible with the input size and can accurately perform inference over a small sample of the instances in the test set. Using the neural network to classify instances of the three-dimensional Edwards-Anderson Ising spin-glass in a (random) field we show that the inferred phase boundary is consistent with the absence of an Almeida-Thouless line.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Machine learning in physics: The pitfalls of poisoned training sets

Known for their ability to identify hidden patterns in data, artificial neural networks are among the most powerful machine learning tools. Most notably, neural networks have played a central role in identifying states of matter and phase transitions across condensed matter physics. To date, most studies have focused on systems where different phases of matter and their phase transitions are known, and thus the performance of neural networks is well controlled. While neural networks present an exciting new tool to detect new phases of matter, here we demonstrate that when the training sets are poisoned (i.e., poor training data or mislabeled data) it is easy for neural networks to make misleading predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Wishart planted ensemble: A tunably-rugged pairwise Ising model with a first-order phase transition

We propose the Wishart planted ensemble, a class of zero-field Ising models with tunable algorithmic hardness and specifiable (or planted) ground state. The problem class arises from a simple procedure for generating a family of random integer programming problems with specific statistical symmetry properties, but turns out to have intimate connections to a sign-inverted variant of the Hopfield model. The Hamiltonian contains only 2-spin interactions, with the coupler matrix following a type of Wishart distribution. The class exhibits a classical first-order phase transition in temperature. For some parameter settings the model has a locally-stable paramagnetic state, a feature which correlates strongly with difficulty in finding the ground state and suggests an extremely rugged energy landscape. We analytically probe the ensemble thermodynamic properties by deriving the Thouless-Anderson-Palmer equations and free energy and corroborate the results with a replica and annealed approximation analysis; extensive Monte Carlo simulations confirm our predictions of the first-order transition temperature. The class exhibits a wide variation in algorithmic hardness as a generation parameter is varied, with a pronounced easy-hard-easy profile and peak in solution time towering many orders of magnitude over that of the easy regimes. By deriving the ensemble-averaged energy distribution and taking into account finite-precision representation, we propose an analytical expression for the location of the hardness peak and show that at fixed precision, the number of constraints in the integer program must increase with system size to yield truly hard problems. The Wishart planted ensemble is interesting for its peculiar physical properties and provides a useful and analytically-transparent set of problems for benchmarking optimization algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Perspectives of quantum annealing: Methods and implementations

Quantum annealing is a computing paradigm that has the ambitious goal of efficiently solving large-scale combinatorial optimization problems of practical importance. However, many challenges have yet to be overcome before this goal can be reached. This perspectives article first gives a brief introduction to the concept of quantum annealing, and then highlights new pathways that may clear the way towards feasible and large scale quantum annealing. Moreover, since this field of research is to a strong degree driven by a synergy between experiment and theory, we discuss both in this work. An important focus in this article is on future perspectives, which complements other review articles, and which we hope will motivate further research.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Best-case performance of quantum annealers on native spin-glass benchmarks: How chaos can affect success probabilities

Recent tests performed on the D-Wave Two quantum annealer have revealed no clear evidence of speedup over conventional silicon-based technologies. Here, we present results from classical parallel-tempering Monte Carlo simulations combined with isoenergetic cluster moves of the archetypal benchmark problem-an Ising spin glass-on the native chip topology. Using realistic uncorrelated noise models for the D-Wave Two quantum annealer, we study the best-case resilience, i.e., the probability that the ground-state configuration is not affected by random fields and random-bond fluctuations found on the chip. We thus compute classical upper-bound success probabilities for different types of disorder used in the benchmarks and predict that an increase in the number of qubits will require either error correction schemes or a drastic reduction of the intrinsic noise found in these devices. We outline strategies to develop robust, as well as hard benchmarks for quantum annealing devices, as well as any other computing paradigm affected by noise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bond chaos in spin glasses revealed through thermal boundary conditions

Spin glasses have competing interactions that lead to a rough energy landscape which is highly susceptible to small perturbations. These chaotic effects strongly affect numerical simulations and, as such, gaining a deeper understanding of chaos in spin glasses is of much importance. The use of thermal boundary conditions is an effective approach to study chaotic phenomena. Here, we generalize population annealing Monte Carlo, combined with thermal boundary conditions, to study bond chaos due to small perturbations in the spin-spin couplings of the three-dimensional Edwards-Anderson Ising spin glass. We show that bond and temperature-induced chaos share the same scaling exponents and that bond chaos is stronger than temperature chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

borealis - A generalized global update algorithm for Boolean optimization problems

Optimization problems with Boolean variables that fall into the nondeterministic polynomial (NP) class are of fundamental importance in computer science, mathematics, physics and industrial applications. Most notably, solving constraint-satisfaction problems, which are related to spin-glass-like Hamiltonians in physics, remains a difficult numerical task. As such, there has been great interest in designing efficient heuristics to solve these computationally difficult problems. Inspired by parallel tempering Monte Carlo in conjunction with the rejection-free isoenergetic cluster algorithm developed for Ising spin glasses, we present a generalized global update optimization heuristic that can be applied to different NP-complete problems with Boolean variables. The global cluster updates allow for a wide-spread sampling of phase space, thus considerably speeding up optimization. By carefully tuning the pseudo-temperature (needed to randomize the configurations) of the problem, we show that the method can efficiently tackle optimization problems with over-constraints or on topologies with a large site-percolation threshold. We illustrate the efficiency of the heuristic on paradigmatic optimization problems, such as the maximum satisfiability problem and the vertex cover problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Charge avalanches and depinning in the Coulomb glass: The role of long-range interactions

We explore the stability of far-from-equilibrium metastable states of a three-dimensional Coulomb glass at zero temperature by studying charge avalanches triggered by a slowly varying external electric field. Surprisingly, we identify a sharply defined dynamical ("depinning") phase transition from stationary to nonstationary charge displacement at a critical value of the external electric field. Using particle-conserving dynamics, scale-free system-spanning avalanches are observed only at the critical field. We show that the qualitative features of this depinning transition are completely different for an equivalent short-range model, highlighting the key importance of long-range interactions for nonequilibrium dynamics of Coulomb glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Error tolerance of topological codes with independent bit-flip and measurement errors

Topological quantum error correction codes are currently among the most promising candidates for efficiently dealing with the decoherence effects inherently present in quantum devices. Numerically, their theoretical error threshold can be calculated by mapping the underlying quantum problem to a related classical statistical-mechanical spin system with quenched disorder. Here, we present results for the general fault-tolerant regime, where we consider both qubit and measurement errors. However, unlike in previous studies, here we vary the strength of the different error sources independently. Our results highlight peculiar differences between toric and color codes. This study complements previous results published in New J. Phys. 13, 083006 (2011).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evolutionary Approaches to Optimization Problems in Chimera Topologies

Chimera graphs define the topology of one of the first commercially available quantum computers. A variety of optimization problems have been mapped to this topology to evaluate the behavior of quantum enhanced optimization heuristics in relation to other optimizers, being able to efficiently solve problems classically to use them as benchmarks for quantum machines. In this paper we investigate for the first time the use of Evolutionary Algorithms (EAs) on Ising spin glass instances defined on the Chimera topology. Three genetic algorithms (GAs) and three estimation of distribution algorithms (EDAs) are evaluated over $1000$ hard instances of the Ising spin glass constructed from Sidon sets. We focus on determining whether the information about the topology of the graph can be used to improve the results of EAs and on identifying the characteristics of the Ising instances that influence the success rate of GAs and EDAs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite-size critical scaling in Ising spin glasses in the mean-field regime

We study in Ising spin glasses the finite-size effects near the spin-glass transition in zero field and at the de Almeida-Thouless transition in a field by Monte Carlo methods and by analytical approximations. In zero field, the finite-size scaling function associated with the spin-glass susceptibility of the Sherrington-Kirkpatrick mean-field spin-glass model is of the same form as that of one-dimensional spin-glass models with power-law long-range interactions in the regime where they can be a proxy for the Edwards-Anderson short-range spin-glass model above the upper critical dimension. We also calculate a simple analytical approximation for the spin-glass susceptibility crossover function. The behavior of the spin-glass susceptibility near the de Almeida-Thouless transition line has also been studied, but here we have only been able to obtain analytically its behavior in the asymptotic limit above and below the transition. We have also simulated the one-dimensional system in a field in the non-mean-field regime to illustrate that when the Imry-Ma droplet length scale exceeds the system size one can then be erroneously lead to conclude that there is a de Almeida-Thouless transition even though it is absent.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interface free-energy exponent in the one-dimensional Ising spin glass with long-range interactions in both the droplet and broken replica symmetry regions

The one-dimensional Ising spin-glass model with power-law long-range interactions is a useful proxy model for studying spin glasses in higher space dimensions and for finding the dimension at which the spin-glass state changes from having broken replica symmetry to that of droplet behavior. To this end we have calculated the exponent that describes the difference in free energy between periodic and antiperiodic boundary conditions. Numerical work is done to support some of the assumptions made in the calculations and to determine the behavior of the interface free-energy exponent of the power law of the interactions. Our numerical results for the interface free-energy exponent are badly affected by finite-size problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Retrieving the ground state of spin glasses using thermal noise: Performance of quantum annealing at finite temperatures

We study the problem to infer the ground state of a spin-glass Hamiltonian using data from another Hamiltonian with interactions disturbed by noise from the original Hamiltonian, motivated by the ground-state inference in quantum annealing on a noisy device. It is shown that the average Hamming distance between the inferred spin configuration and the true ground state is minimized when the temperature of the noisy system is kept at a finite value, and not at zero temperature. We present a spin-glass generalization of a well-established result that the ground state of a purely ferromagnetic Hamiltonian is best inferred at a finite temperature in the sense of smallest Hamming distance when the original ferromagnetic interactions are disturbed by noise. We use the numerical transfer-matrix method to establish the existence of an optimal finite temperature in one- and two-dimensional systems. Our numerical results are supported by mean-field calculations, which give an explicit expression of the optimal temperature to infer the spin-glass ground state as a function of variances of the distributions of the original interactions and the noise. The mean-field prediction is in qualitative agreement with numerical data. Implications on postprocessing of quantum annealing on a noisy device are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Site and bond percolation thresholds in $K_{n,n}$-based lattices: Vulnerability of quantum annealers to random qubit and coupler failures on chimera topologies

We estimate the critical thresholds of bond and site percolation on nonplanar, effectively two-dimensional graphs with chimera like topology. The building blocks of these graphs are complete and symmetric bipartite subgraphs of size $2n$, referred to as $K_{n,n}$ graphs. For the numerical simulations we use an efficient union-find based algorithm and employ a finite-size scaling analysis to obtain the critical properties for both bond and site percolation. We report the respective percolation thresholds for different sizes of the bipartite subgraph and verify that the associated universality class is that of standard two-dimensional percolation. For the canonical chimera graph used in the D-Wave Systems Inc.~quantum annealer ($n = 4$), we discuss device failure in terms of network vulnerability, i.e., we determine the critical fraction of qubits and couplers that can be absent due to random failures prior to losing large-scale connectivity throughout the device.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Boltzmann Equation for Magnetic Relaxation in High-Spin Molecules

We introduce the stochastic Boltzmann equation (SBE) as an approach for exploring the spin dynamics of magnetic molecules coupled to a stochastic environment. The SBE is a time-evolution equation for the probability density of the spin density matrix of the system. This probability density is relevant to experiments which take measurements on single molecules, in which probabilities of observing particular spin states (rather than ensemble averages) are of interest. By analogy with standard treatments of the regular Boltzmann equation, we propose a relaxation-time approximation for the SBE, and show that solutions to the SBE under the relaxation-time approximation can be obtained by performing simple trajectory simulations for the case of a boson gas environment. Cases where the relaxation-time approximation are satisfied can therefore be investigated by careful choice of the parameters for the boson gas environment, even if the actual environment is quite different from a boson gas. The application of the SBE approach is demonstrated through an illustrative example.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strengths and weaknesses of weak-strong cluster problems: A detailed overview of state-of-the-art classical heuristics vs quantum approaches

To date, a conclusive detection of quantum speedup remains elusive. Recently, a team by Google Inc.~[V.~S.~Denchev {\em et al}., Phys.~Rev.~X {\bf 6}, 031015 (2016)] proposed a weak-strong cluster model tailored to have tall and narrow energy barriers separating local minima, with the aim to highlight the value of finite-range tunneling. More precisely, results from quantum Monte Carlo simulations, as well as the D-Wave 2X quantum annealer scale considerably better than state-of-the-art simulated annealing simulations. Moreover, the D-Wave 2X quantum annealer is $\sim 10^8$ times faster than simulated annealing on conventional computer hardware for problems with approximately $10^3$ variables. Here, an overview of different sequential, nontailored, as well as specialized tailored algorithms on the Google instances is given. We show that the quantum speedup is limited to sequential approaches and study the typical complexity of the benchmark problems using insights from the study of spin glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chaos in spin glasses revealed through thermal boundary conditions

We study the fragility of spin glasses to small temperature perturbations numerically using population annealing Monte Carlo. We apply thermal boundary conditions to a three-dimensional Edwards-Anderson Ising spin glass. In thermal boundary conditions all eight combinations of periodic versus antiperiodic boundary conditions in the three spatial directions are present, each appearing in the ensemble with its respective statistical weight determined by its free energy. We show that temperature chaos is revealed in the statistics of crossings in the free energy for different boundary conditions. By studying the energy difference between boundary conditions at free-energy crossings, we determine the domain-wall fractal dimension. Similarly, by studying the number of crossings, we determine the chaos exponent. Our results also show that computational hardness in spin glasses and the presence of chaos are closely related.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparing Monte Carlo methods for finding ground states of Ising spin glasses: population annealing, simulated annealing and parallel tempering

Population annealing is a Monte Carlo algorithm that marries features from simulated annealing and parallel tempering Monte Carlo. As such, it is ideal to overcome large energy barriers in the free-energy landscape while minimizing a Hamiltonian. Thus, population annealing Monte Carlo can be used as a heuristic to solve combinatorial optimization problems. We illustrate the capabilities of population annealing Monte Carlo by computing ground states of the three-dimensional Ising spin glass with Gaussian disorder, whilst comparing to simulated annealing and parallel tempering Monte Carlo. Our results suggest that population annealing Monte Carlo is significantly more efficient than simulated annealing but comparable to parallel tempering Monte Carlo for finding spin-glass ground states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient Cluster Algorithm for Spin Glasses in Any Space Dimension

Spin systems with frustration and disorder are notoriously difficult to study both analytically and numerically. While the simulation of ferromagnetic statistical mechanical models benefits greatly from cluster algorithms, these accelerated dynamics methods remain elusive for generic spin-glass-like systems. Here we present a cluster algorithm for Ising spin glasses that works in any space dimension and speeds up thermalization by at least one order of magnitude at temperatures where thermalization is typically difficult. Our isoenergetic cluster moves are based on the Houdayer cluster algorithm for two-dimensional spin glasses and lead to a speedup over conventional state-of-the-art methods that increases with the system size. We illustrate the benefits of the isoenergetic cluster moves in two and three space dimensions, as well as the nonplanar chimera topology found in the D-Wave Inc.~quantum annealing machine.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Population annealing: Theory and application in spin glasses

Population annealing is an efficient sequential Monte Carlo algorithm for simulating equilibrium states of systems with rough free energy landscapes. The theory of population annealing is presented, and systematic and statistical errors are discussed. The behavior of the algorithm is studied in the context of large-scale simulations of the three-dimensional Ising spin glass and the performance of the algorithm is compared to parallel tempering. It is found that the two algorithms are similar in efficiency though with different strengths and weaknesses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Seeking Quantum Speedup Through Spin Glasses: The Good, the Bad, and the Ugly

There has been considerable progress in the design and construction of quantum annealing devices. However, a conclusive detection of quantum speedup over traditional silicon-based machines remains elusive, despite multiple careful studies. In this work we outline strategies to design hard tunable benchmark instances based on insights from the study of spin glasses - the archetypal random benchmark problem for novel algorithms and optimization devices. We propose to complement head-to-head scaling studies that compare quantum annealing machines to state-of-the-art classical codes with an approach that compares the performance of different algorithms and/or computing architectures on different classes of computationally hard tunable spin-glass instances. The advantage of such an approach lies in having to only compare the performance hit felt by a given algorithm and/or architecture when the instance complexity is increased. Furthermore, we propose a methodology that might not directly translate into the detection of quantum speedup, but might elucidate whether quantum annealing has a "`quantum advantage" over corresponding classical algorithms like simulated annealing. Our results on a 496 qubit D-Wave Two quantum annealing device are compared to recently-used state-of-the-art thermal simulated annealing codes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of quenched bond disorder on first-order phase transitions

We investigate the effect of quenched bond disorder on the two-dimensional three-color Ashkin-Teller model, which undergoes a first-order phase transition in the absence of impurities. This is one of the simplest and striking models in which quantitative numerical simulations can be carried out to investigate emergent criticality due to disorder rounding of first-order transition. Utilizing extensive cluster Monte Carlo simulations on large lattice sizes of up to $128 \times 128$ spins, each of which is represented by three colors taking values $\pm 1$, we show that the rounding of the first-order phase transition is an emergent criticality. We further calculate the correlation length critical exponent, $ν$, and the magnetization critical exponent, $β$, from finite size scaling analysis. We find that the critical exponents, $ν$ and $β$, change as the strength of disorder or the four-spin coupling varies, and we show that the critical exponents appear not to be in the Ising universality class. We know of no analytical approaches that can explain our non-perturbative results. However our results should inspire further work on this important problem, either numerical or analytical.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A probabilistic evolutionary optimization approach to compute quasiparticle braids

Topological quantum computing is an alternative framework for avoiding the quantum decoherence problem in quantum computation. The problem of executing a gate in this framework can be posed as the problem of braiding quasiparticles. Because these are not Abelian, the problem can be reduced to finding an optimal product of braid generators where the optimality is defined in terms of the gate approximation and the braid's length. In this paper we propose the use of different variants of estimation of distribution algorithms to deal with the problem. Furthermore, we investigate how the regularities of the braid optimization problem can be translated into statistical regularities by means of the Boltzmann distribution. We show that our best algorithm is able to produce many solutions that approximates the target gate with an accuracy in the order of $10^{-6}$, and have lengths up to 9 times shorter than those expected from braids of the same accuracy obtained with other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Boolean decision problems with competing interactions on scale-free networks: Equilibrium and nonequilibrium behavior in an external bias

We study the equilibrium and nonequilibrium properties of Boolean decision problems with competing interactions on scale-free networks in an external bias (magnetic field). Previous studies at zero field have shown a remarkable equilibrium stability of Boolean variables (Ising spins) with competing interactions (spin glasses) on scale-free networks. When the exponent that describes the power-law decay of the connectivity of the network is strictly larger than 3, the system undergoes a spin-glass transition. However, when the exponent is equal to or less than 3, the glass phase is stable for all temperatures. First, we perform finite-temperature Monte Carlo simulations in a field to test the robustness of the spin-glass phase and show that the system has a spin-glass phase in a field, i.e., exhibits a de Almeida-Thouless line. Furthermore, we study avalanche distributions when the system is driven by a field at zero temperature to test if the system displays self-organized criticality. Numerical results suggest that avalanches (damage) can spread across the whole system with nonzero probability when the decay exponent of the interaction degree is less than or equal to 2, i.e., that Boolean decision problems on scale-free networks with competing interactions can be fragile when not in thermal equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dealing with correlated choices: How a spin-glass model can help political parties select their policies

Starting from preferences on N proposed policies obtained via questionnaires from a sample of the electorate, an Ising spin-glass model in a field can be constructed from which a political party could find the subset of the proposed policies which would maximize its appeal, form a coherent choice in the eyes of the electorate, and have maximum overlap with the party's existing policies. We illustrate the application of the procedure by simulations of a spin glass in a random field on scale-free networks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error Thresholds for Abelian Quantum Double Models: Increasing the bit-flip Stability of Topological Quantum Memory

Current approaches for building quantum computing devices focus on two-level quantum systems which nicely mimic the concept of a classical bit, albeit enhanced with additional quantum properties. However, rather than artificially limiting the number of states to two, the use of d-level quantum systems (qudits) could provide advantages for quantum information processing. Among other merits, it has recently been shown that multi-level quantum systems can offer increased stability to external disturbances - a key problem in current technologies. In this study we demonstrate that topological quantum memories built from qudits, also known as abelian quantum double models, exhibit a substantially increased resilience to noise. That is, even when taking into account the multitude of errors possible for multi-level quantum systems, topological quantum error correction codes employing qudits can sustain a larger error rate than their two-level counterparts. In particular, we find strong numerical evidence that the thresholds of these error-correction codes are given by the hashing bound. Considering the significantly increased error thresholds attained, this might well outweigh the added complexity of engineering and controlling higher dimensional quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evidence against a mean field description of short-range spin glasses revealed through thermal boundary conditions

A theoretical description of the low-temperature phase of short-range spin glasses has remained elusive for decades. In particular, it is unclear if theories that assert a single pair of pure states, or theories that are based infinitely many pure states-such as replica symmetry breaking-best describe realistic short-range systems. To resolve this controversy, the three-dimensional Edwards-Anderson Ising spin glass in thermal boundary conditions is studied numerically using population annealing Monte Carlo. In thermal boundary conditions all eight combinations of periodic vs antiperiodic boundary conditions in the three spatial directions appear in the ensemble with their respective Boltzmann weights, thus minimizing finite-size corrections due to domain walls. From the relative weighting of the eight boundary conditions for each disorder instance a sample stiffness is defined, and its typical value is shown to grow with system size according to a stiffness exponent. An extrapolation to the large-system-size limit is in agreement with a description that supports the droplet picture and other theories that assert a single pair of pure states. The results are, however, incompatible with the mean-field replica symmetry breaking picture, thus highlighting the need to go beyond mean-field descriptions to accurately describe short-range spin-glass systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Existence of a Thermodynamic Spin-Glass Phase in the Zero-Concentration Limit of Anisotropic Dipolar Systems

The nature of ordering in dilute dipolar interacting systems dates back to the work of Debye and is one of the most basic, oldest and as-of-yet unsettled problems in magnetism. While spin-glass order is readily observed in several RKKY-interacting systems, dipolar spin-glasses are subject of controversy and ongoing scrutiny, e.g., in ${{\rm LiHo_xY_{1-x}F_4}}$, a rare-earth randomly diluted uniaxial (Ising) dipolar system. In particular, it is unclear if the spin-glass phase in these paradigmatic materials persists in the limit of zero concentration or not. We study an effective model of ${{\rm LiHo_xY_{1-x}F_4}}$ using large-scale Monte Carlo simulations that combine parallel tempering with a special cluster algorithm tailored to overcome the numerical difficulties that occur at extreme dilutions. We find a paramagnetic to spin-glass phase transition for all Ho ion concentrations down to the smallest concentration numerically accessible of 0.1%, and including Ho ion concentrations which coincide with those studied experimentally up to 16.7%. Our results suggest that randomly-diluted dipolar Ising systems have a spin-glass phase in the limit of vanishing dipole concentration, with a critical temperature vanishing linearly with concentration, in agreement with mean field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Glassy Chimeras could be blind to quantum speedup: Designing better benchmarks for quantum annealing machines

Recently, a programmable quantum annealing machine has been built that minimizes the cost function of hard optimization problems by adiabatically quenching quantum fluctuations. Tests performed by different research teams have shown that, indeed, the machine seems to exploit quantum effects. However experiments on a class of random-bond instances have not yet demonstrated an advantage over classical optimization algorithms on traditional computer hardware. Here we present evidence as to why this might be the case. These engineered quantum annealing machines effectively operate coupled to a decohering thermal bath. Therefore, we study the finite-temperature critical behavior of the standard benchmark problem used to assess the computational capabilities of these complex machines. We simulate both random-bond Ising models and spin glasses with bimodal and Gaussian disorder on the D-Wave Chimera topology. Our results show that while the worst-case complexity of finding a ground state of an Ising spin glass on the Chimera graph is not polynomial, the finite-temperature phase space is likely rather simple: Spin glasses on Chimera have only a zero-temperature transition. This means that benchmarking optimization methods using spin glasses on the Chimera graph might not be the best benchmark problems to test quantum speedup. We propose alternative benchmarks by embedding potentially harder problems on the Chimera topology. Finally, we also study the (reentrant) disorder-temperature phase diagram of the random-bond Ising model on the Chimera graph and show that a finite-temperature ferromagnetic phase is stable up to 19.85(15)% antiferromagnetic bonds. Beyond this threshold the system only displays a zero-temperature spin-glass phase. Our results therefore show that a careful design of the hardware architecture and benchmark problems is key when building quantum annealing machines.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low-temperature behavior of the statistics of the overlap distribution in Ising spin-glass models

Using Monte Carlo simulations, we study in detail the overlap distribution for individual samples for several spin-glass models including the infinite-range Sherrington-Kirkpatrick model, short-range Edwards-Anderson models in three and four space dimensions, and one-dimensional long-range models with diluted power-law interactions. We study three long-range models with different powers as follows: the first is approximately equivalent to a short-range model in three dimensions, the second to a short-range model in four dimensions, and the third to a short-range model in the mean-field regime. We study an observable proposed earlier by some of us which aims to distinguish the "replica symmetry breaking" picture of the spin-glass phase from the "droplet picture," finding that larger system sizes would be needed to unambiguously determine which of these pictures describes the low-temperature state of spin glasses best, except for the Sherrington-Kirkpatrick model which is unambiguously described by replica symmetry breaking. Finally, we also study the median integrated overlap probability distribution and a typical overlap distribution, finding that these observables are not particularly helpful in distinguishing the replica symmetry breaking and the droplet pictures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

2D molecular magnets with weak topological invariant magnetic moments: Mathematical prediction of targets for chemical synthesis

An open problem in applied mathematics is to predict interesting molecules which are realistic targets for chemical synthesis. In this paper, we use a spin Hamiltonian-type model to predict molecular magnets (MMs) with magnetic moments that are intrinsically robust under random shape deformations to the molecule. Using the concept of convergence in probability, we show that for MMs in which all spin centers lie in-plane and all spin center interactions are ferromagnetic, the total spin of the molecule is a `weak topological invariant' when the number of spin centers is sufficiently large. By weak topological invariant, we mean that the total spin of the molecule only depends upon the arrangement of spin centers in the molecule, and is unlikely to change under shape deformations to the molecule. Our calculations show that only between 20 and 50 spin centers are necessary for the total spin of these MMs to be a weak topological invariant. The robustness effect is particularly enhanced for 2D ferromagnetic MMs that possess a small number of spin rings in the structure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Correlations between the dynamics of parallel tempering and the free-energy landscape in spin glasses

We present the results of a large-scale numerical study of the equilibrium three-dimensional Edwards-Anderson Ising spin glass with Gaussian disorder. Using parallel tempering (replica exchange) Monte Carlo we measure various static, as well as dynamical quantities, such as the autocorrelation times and round-trip times for the parallel tempering Monte Carlo method. The correlation between static and dynamic observables for 5000 disorder realizations and up to 1000 spins down to temperatures at 20% of the critical temperature is examined. Our results show that autocorrelation times are directly correlated with the roughness of the free-energy landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Diluted antiferromagnets in a field seem to be in a different universality class than the random-field Ising model

We perform large-scale Monte Carlo simulations using the Machta-Newman-Chayes algorithms to study the critical behavior of both the diluted antiferromagnet in a field with 30% dilution and the random-field Ising model with Gaussian random fields for different field strengths. Analytical calculations by Cardy [Phys. Rev. B 29, 505 (1984)] predict that both models map onto each other and share the same universality class in the limit of vanishing fields. However, a detailed finite-size scaling analysis of both the Binder cumulant and the two-point finite-size correlation length suggests that even in the limit of small fields, where the mapping is expected to work, both models are not in the same universality class. Therefore, care should be taken when interpreting (experimental) data for diluted antiferromagnets in a field using the random-field Ising model. Based on our numerical data, we present analytical expressions for the phase boundaries of both models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Genetic braid optimization: A heuristic approach to compute quasiparticle braids

In topologically-protected quantum computation, quantum gates can be carried out by adiabatically braiding two-dimensional quasiparticles, reminiscent of entangled world lines. Bonesteel et al. [Phys. Rev. Lett. 95, 140503 (2005)], as well as Leijnse and Flensberg [Phys. Rev. B 86, 104511 (2012)] recently provided schemes for computing quantum gates from quasiparticle braids. Mathematically, the problem of executing a gate becomes that of finding a product of the generators (matrices) in that set that approximates the gate best, up to an error. To date, efficient methods to compute these gates only strive to optimize for accuracy. We explore the possibility of using a generic approach applicable to a variety of braiding problems based on evolutionary (genetic) algorithms. The method efficiently finds optimal braids while allowing the user to optimize for the relative utilities of accuracy and/or length. Furthermore, when optimizing for error only, the method can quickly produce efficient braids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Novel disordering mechanism in ferromagnetic systems with competing interactions

Ferromagnetic Ising systems with competing interactions are considered in the presence of a random field. We find that in three space dimensions the ferromagnetic phase is disordered by a random field which is considerably smaller than the typical interaction strength between the spins. This is the result of a novel disordering mechanism triggered by an underlying spin-glass phase. Calculations for the specific case of the long-range dipolar LiHo_xY_{1-x}F_4 compound suggest that the above mechanism is responsible for the peculiar dependence of the critical temperature on the strength of the random field and the broadening of the susceptibility peaks as temperature is decreased, as found in recent experiments by Silevitch et al. [Nature (London) 448, 567 (2007)]. Our results thus emphasize the need to go beyond the standard Imry-Ma argument when studying general random-field systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reply to Comment on "Evidence of Non-Mean-Field-Like Low-Temperature Behavior in the Edwards-Anderson Spin-Glass Model"

In Phys. Rev. Lett. 110, 219701 (2013) [arXiv:1211.0843] Billoire et al. criticize the conclusions of our Letter [Phys. Rev. Lett. 109, 177204 (2012), arxiv:1206.0783]. They argue that considering the Edwards-Anderson and Sherrington-Kirkpatrick models at the same temperature is inappropriate and propose an interpretation based on the replica symmetry breaking theory. Here we show that the theory presented in the Comment does not explain our data on the Edwards-Anderson spin glass and we stand by our assertion that the low-temperature behavior of the Edwards-Anderson spin glass model does not appear to be mean-field like.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-Organized Criticality in Glassy Spin Systems Requires a Diverging Number of Neighbors

We investigate the conditions required for general spin systems with frustration and disorder to display self-organized criticality, a property which so far has been established only for the fully-connected infinite-range Sherrington-Kirkpatrick Ising spin-glass model [Phys. Rev. Lett. 83, 1034 (1999)]. Here we study both avalanche and magnetization jump distributions triggered by an external magnetic field, as well as internal field distributions in the short-range Edwards-Anderson Ising spin glass for various space dimensions between 2 and 8, as well as the fixed-connectivity mean-field Viana-Bray model. Our numerical results, obtained on systems of unprecedented size, demonstrate that self-organized criticality is recovered only in the strict limit of a diverging number of neighbors, and is not a generic property of spin-glass models in finite space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spin glasses in a field: Three and four dimensions as seen from one space dimension

We study the existence of a line of transitions of an Ising spin glass in a magnetic field-known as the de Almeida-Thouless line-using one-dimensional power-law diluted Ising spin-glass models. We choose the power-law exponent to have values that approximately correspond to three- and four-dimensional nearest-neighbor systems and perform a detailed finite-size scaling analysis of the data for large linear system sizes, both using a new approach proposed recently [Phys. Rev. Lett. 103, 267201 (2009)], as well as traditional approaches. Our results for the model corresponding to a three-dimensional system are consistent with there being no de Almeida-Thouless line, although the new finite-size scaling approach does not rule one out. For the model corresponding to four space dimensions, the new and traditional finite-size scaling analyses give conflicting results, indicating the need for a better understanding of finite-size scaling of spin glasses in a magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stability of topologically-protected quantum computing proposals as seen through spin glasses

Sensitivity to noise makes most of the current quantum computing schemes prone to error and nonscalable, allowing only for small proof-of-principle devices. Topologically-protected quantum computing aims at solving this problem by encoding quantum bits and gates in topological properties of the hardware medium that are immune to noise that does not impact the entire system at once. There are different approaches to achieve topological stability or active error correction, ranging from quasiparticle braidings to spin models and topological color codes. The stability of these proposals against noise can be quantified by their error threshold. This figure of merit can be computed by mapping the problem onto complex statistical-mechanical spin-glass models with local disorder on nontrival lattices that can have many-body interactions and are sometimes described by lattice gauge theories. The error threshold for a given source of error then represents the point in the temperature-disorder phase diagram where a stable symmetry-broken phase vanishes. An overview of the techniques used to estimate the error thresholds is given, as well as a summary of recent results on the stability of different topologically-protected quantum computing schemes to different error sources.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bond disorder induced criticality of the three-color Ashkin-Teller model

An intriguing result of statistical mechanics is that a first-order phase transition can be rounded by disorder coupled to energy-like variables. In fact, even more intriguing is that the rounding may manifest itself as a critical point, quantum or classical. In general, it is not known, however, what universality classes, if any, such criticalities belong to. In order to shed light on this question we examine in detail the disordered three-color Ashkin-Teller model by Monte Carlo methods. Extensive analyses indicate that the critical exponents define a new universality class. We show that the rounding of the first-order transition of the pure model due to the impurities is manifested as criticality. However, the magnetization critical exponent, (β), and the correlation length critical exponent, (ν), are found to vary with disorder and the four-spin coupling strength, and we conclusively rule out that the model belongs to the universality class of the two-dimensional Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Boolean decision problems with competing interactions on scale-free networks: Critical thermodynamics

We study the critical behavior of Boolean variables on scale-free networks with competing interactions (Ising spin glasses). Our analytical results for the disorder-network-decay-exponent phase diagram are verified using Monte Carlo simulations. When the probability of positive (ferromagnetic) and negative (antiferromagnetic) interactions is the same, the system undergoes a finite-temperature spin-glass transition if the exponent that describes the decay of the interaction degree in the scale-free graph is strictly larger than 3. However, when the exponent is equal to or less than 3, a spin-glass phase is stable for all temperatures. The robustness of both the ferromagnetic and spin-glass phases suggests that Boolean decision problems on scale-free networks are quite stable to local perturbations. Finally, we show that for a given decay exponent spin glasses on scale-free networks seem to obey universality. Furthermore, when the decay exponent of the interaction degree is larger than 4 in the spin-glass sector, the universality class is the same as for the mean-field Sherrington-Kirkpatrick Ising spin glass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence of non-mean-field-like low-temperature behavior in the Edwards-Anderson spin-glass model

The three-dimensional Edwards-Anderson and mean-field Sherrington-Kirkpatrick Ising spin glasses are studied via large-scale Monte Carlo simulations at low temperatures, deep within the spin-glass phase. Performing a careful statistical analysis of several thousand independent disorder realizations and using an observable that detects peaks in the overlap distribution, we show that the Sherrington-Kirkpatrick and Edwards-Anderson models have a distinctly different low-temperature behavior. The structure of the spin-glass overlap distribution for the Edwards-Anderson model suggests that its low-temperature phase has only a single pair of pure states.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal error correction in topological subsystem codes

A promising approach to overcome decoherence in quantum computing schemes is to perform active quantum error correction using topology. Topological subsystem codes incorporate both the benefits of topological and subsystem codes, allowing for error syndrome recovery with only 2-local measurements in a two-dimensional array of qubits. We study the error threshold for topological subsystem color codes under very general external noise conditions. By transforming the problem into a classical disordered spin model, we estimate using Monte Carlo simulations that topological subsystem codes have an optimal error tolerance of 5.5(2)%. This means there is ample space for improvement in existing error-correcting algorithms that typically find a threshold of approximately 2%.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

STM Studies of Isolated Mn12-Ph Single Molecule Magnets

We study Mn12O12(C6H5COO)16(H2O)4 (Mn12-Ph) single-molecule magnets on highly ordered pyrolytic graphite (HOPG) using low temperature scanning tunneling microscopy (LT-STM) experiments. We report Mn12-Ph in isolation, resembling single molecules with metallic core atoms and organic outer ligands. The local tunneling current observed within the molecular structure shows a strong bias voltage dependency, which is distinct from that of the HOPG surface. Further, evidence of internal inhomogeneity in the local density of states has been observed with high spatial resolution, and this inhomogeneity appears to be due to localized metallic behavior. These results facilitate magneto-metric studies of single molecule magnets in isolation. As compared to bulk crystal studies, our experiments allow the specific investigation of atomic sites in the molecule.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong Resilience of Topological Codes to Depolarization

The inevitable presence of decoherence effects in systems suitable for quantum computation necessitates effective error-correction schemes to protect information from noise. We compute the stability of the toric code to depolarization by mapping the quantum problem onto a classical disordered eight-vertex Ising model. By studying the stability of the related ferromagnetic phase both via large-scale Monte Carlo simulations and via the duality method, we are able to demonstrate an increased error threshold of 18.9(3)% when noise correlations are taken into account. Remarkably, this agrees within error bars with the result for a different class of codes-topological color codes-where the mapping yields interesting new types of interacting eight-vertex models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ultrametric probe of the spin-glass state in a field

We study the ultrametric structure of phase space of one-dimensional Ising spin glasses with random power-law interaction in an external random field. Although in zero field the model in both the mean-field and non-mean-field universality classes shows an ultrametric signature [Phys. Rev. Lett. 102, 037207 (2009)], when a field is applied ultrametricity seems only present in the mean-field regime. The results for the non-mean field case in an external field agree with data for spin glasses studied within the Migdal-Kadanoff approximation. Our results therefore suggest that the spin-glass state might be fragile to external fields below the upper critical dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Critical behavior and universality in Levy spin glasses

Using large-scale Monte Carlo simulations that combine parallel tempering with specialized cluster updates, we show that Ising spin glasses with Levy-distributed interactions share the same universality class as Ising spin glasses with Gaussian or bimodal-distributed interactions. Corrections to scaling are large for Levy spin glasses. In order to overcome these and show that the critical exponents agree with the Gaussian case, we perform an extended scaling of the two-point finite size correlation length and the spin glass susceptibility. Furthermore, we compute the critical temperature and compare its dependence on the disorder distribution width to recent analytical predictions [J. Stat. Mech. (2008) P04006].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evidence of a glass transition in a 10-state non-mean-field Potts glass

Potts glasses are prototype models that have been used to understand the structural glass transition. However, in finite space dimensions a glass transition remains to be detected in the 10-state Potts glass. Using a one-dimensional model with long-range power-law interactions we present evidence that a glass transition below the upper critical dimension can exist for short-range systems at low enough temperatures. Gaining insights into the structural glass transition for short-range systems using spin models is thus potentially possible, yet difficult.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Introduction to Monte Carlo Methods

Monte Carlo methods play an important role in scientific computation, especially when problems have a vast phase space. In this lecture an introduction to the Monte Carlo method is given. Concepts such as Markov chains, detailed balance, critical slowing down, and ergodicity, as well as the Metropolis algorithm are explained. The Monte Carlo method is illustrated by numerically studying the critical behavior of the two-dimensional Ising ferromagnet using finite-size scaling methods. In addition, advanced Monte Carlo methods are described (e.g., the Wolff cluster algorithm and parallel tempering Monte Carlo) and illustrated with nontrivial models from the physics of glassy systems. Finally, we outline an approach to study rare events using a Monte Carlo sampling with a guiding function.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimizing glassy p-spin models

Computing the ground state of Ising spin-glass models with p-spin interactions is, in general, an NP-hard problem. In this work we show that unlike in the case of the standard Ising spin glass with two-spin interactions, computing ground states with p=3 is an NP-hard problem even in two space dimensions. Furthermore, we present generic exact and heuristic algorithms for finding ground states of p-spin models with high confidence for systems of up to several thousand spins.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sampling the ground-state magnetization of d-dimensional p-body Ising models

We demonstrate that a recently introduced heuristic optimization algorithm [Phys. Rev. E 83, 046709 (2011)] that combines a local search with triadic crossover genetic updates is capable of sampling nearly uniformly among ground-state configurations in spin-glass-like Hamiltonians with p-spin interactions in d space dimensions that have highly degenerate ground states. Using this algorithm we probe the zero-temperature ferromagnet to spin-glass transition point q_c of two example models, the disordered version of the two-dimensional three-spin Baxter-Wu model [q_c = 0.1072(1)] and the three-dimensional Edwards-Anderson model [q_c = 0.2253(7)], by computing the Binder ratio of the ground-state magnetization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simplest model to study reentrance in physical systems

We numerically investigate the necessary ingredients for reentrant behavior in the phase diagram of physical systems. Studies on the possibly simplest model that exhibits reentrance, the two-dimensional random bond Ising model, show that reentrant behavior is generic whenever frustration is present in the model. For both discrete and continuous disorder distributions, the phase diagram in the disorder-temperature plane is found to be reentrant, where for some disorder strengths a paramagnetic phase exists at both high and low temperatures, but an ordered ferromagnetic phase exists for intermediate temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tricolored Lattice Gauge Theory with Randomness: Fault-Tolerance in Topological Color Codes

We compute the error threshold of color codes, a class of topological quantum codes that allow a direct implementation of quantum Clifford gates, when both qubit and measurement errors are present. By mapping the problem onto a statistical-mechanical three-dimensional disordered Ising lattice gauge theory, we estimate via large-scale Monte Carlo simulations that color codes are stable against 4.5(2)% errors. Furthermore, by evaluating the skewness of the Wilson loop distributions, we introduce a very sensitive probe to locate first-order phase transitions in lattice gauge theories.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universality in phase boundary slopes for spin glasses on self dual lattices

We study the effects of disorder on the slope of the disorder--temperature phase boundary near the Onsager point (Tc = 2.269...) in spin-glass models. So far, studies have focused on marginal or irrelevant cases of disorder. Using duality arguments, as well as exact Pfaffian techniques we reproduce these analytical estimates. In addition, we obtain different estimates for spin-glass models on hierarchical lattices where the effects of disorder are relevant. We show that the phase-boundary slope near the Onsager point can be used to probe for the relevance of disorder effects.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical studies of a one-dimensional 3-spin spin-glass model with long-range interactions

We study a p-spin spin-glass model to understand if the finite-temperature glass transition found in the mean-field regime of p-spin models, and used to model the behavior of structural glasses, persists in the non-mean-field regime. By using a 3-spin spin-glass model with long-range power-law diluted interactions we are able to continuously tune the (effective) space dimension via the exponent of the interactions. Monte Carlo simulations of the spin-glass susceptibility and the two-point finite-size correlation length show that deep in the non-mean-field regime the finite-temperature transition is lost, whereas this is not the case in the mean-field regime, in agreement with the prediction of Moore and Drossel [Phys. Rev. Lett. 89, 217202 (2002)] that 3-spin models are in the same universality class as an Ising spin glass in a magnetic field. However, slightly in the non-mean-field region, we find an apparent transition in the 3-spin model, in contrast to results for the Ising spin glass in a field. This may indicate that even larger sizes are needed to probe the asymptotic behavior in this region.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random Numbers in Scientific Computing: An Introduction

Random numbers play a crucial role in science and industry. Many numerical methods require the use of random numbers, in particular the Monte Carlo method. Therefore it is of paramount importance to have efficient random number generators. The differences, advantages and disadvantages of true and pseudo random number generators are discussed with an emphasis on the intrinsic details of modern and fast pseudo random number generators. Furthermore, standard tests to verify the quality of the random numbers produced by a given generator are outlined. Finally, standard scientific libraries with built-in generators are presented, as well as different approaches to generate nonuniform random numbers. Potential problems that one might encounter when using large parallel machines are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reply to Comment "Density of states and critical behavior of the three-dimensional Coulomb Glass"

We reply to the Comment by Mobius and Richter [arXiv:0908.3092, Phys. Rev. Lett. 105, 039701 (2010)] on "Density of States and Critical Behavior of the Coulomb Glass" [arXiv:0805.4640, Phys. Rev. Lett. 102, 067205 (2009)] and address the issues raised with our results on the density of states. In addition, we correct our statements about the random displacement version of the Coulomb glass model where the Wigner crystal is not as robust to disorder as stated. Still, our main result of a lack of a finite-temperature transition in the Coulomb glass remains unchallenged.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Scientific Software Engineering in a Nutshell

Writing complex computer programs to study scientific problems requires careful planning and an in-depth knowledge of programming languages and tools. In this chapter the importance of using the right tool for the right problem is emphasized. Common tools to organize computer programs, as well as to debug and improve them are discussed, followed by simple data reduction strategies and visualization tools. Furthermore, some useful scientific libraries such as boost, GSL, LEDA and numerical recipes are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Topological color codes on Union Jack lattices: A stable implementation of the whole Clifford group

We study the error threshold of topological color codes on Union Jack lattices that allow for the full implementation of the whole Clifford group of quantum gates. After mapping the error-correction process onto a statistical mechanical random 3-body Ising model on a Union Jack lattice, we compute its phase diagram in the temperature-disorder plane using Monte Carlo simulations. Surprisingly, topological color codes on Union Jack lattices have similar error stability than color codes on triangular lattices, as well as the Kitaev toric code. The enhanced computational capabilities of the topological color codes on Union Jack lattices with respect to triangular lattices and the toric code demonstrate the inherent robustness of this implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Absence of a structural glass phase in a monoatomic model liquid predicted to undergo an ideal glass transition

We study numerically a monodisperse model of interacting classical particles predicted to exhibit a static liquid-glass transition. Using a dynamical Monte Carlo method we show that the model does not freeze into a glassy phase at low temperatures. Instead, depending on the choice of the hard-core radius for the particles the system either collapses trivially or a polycrystalline hexagonal structure emerges.

0 citations0 reviewsNo code linkedOpen access