Quantum versus Thermal annealing, the role of Temperature Chaos
|
|
- Britney Gray
- 6 years ago
- Views:
Transcription
1 Quantum versus Thermal annealing, the role of Temperature Chaos Víctor Martín-Mayor Dep. Física Teórica I, Universidad Complutense de Madrid Janus Collaboration In collaboration with Itay Hen (Information Sciences Institute, USC). CPRCS 2014, Capri, September V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
2 Foreword Monte Carlo is an extremely bad method; it should be used only when all alternative methods are worse. (Alan Sokal, 1989). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
3 Foreword Monte Carlo is an extremely bad method; it should be used only when all alternative methods are worse. (Alan Sokal, 1989). Such desperate problems are common in Theoretical Physics. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
4 Foreword Monte Carlo is an extremely bad method; it should be used only when all alternative methods are worse. (Alan Sokal, 1989). Such desperate problems are common in Theoretical Physics. SG simulations are extremely computer intensive but simple: the Janus collaboration has dared to produce dedicated hardware. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
5 Foreword Monte Carlo is an extremely bad method; it should be used only when all alternative methods are worse. (Alan Sokal, 1989). Such desperate problems are common in Theoretical Physics. SG simulations are extremely computer intensive but simple: the Janus collaboration has dared to produce dedicated hardware. Janus is a great success, but classical Monte Carlo is hitting an algorithmic wall (temperature chaos). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
6 Foreword Monte Carlo is an extremely bad method; it should be used only when all alternative methods are worse. (Alan Sokal, 1989). Such desperate problems are common in Theoretical Physics. SG simulations are extremely computer intensive but simple: the Janus collaboration has dared to produce dedicated hardware. Janus is a great success, but classical Monte Carlo is hitting an algorithmic wall (temperature chaos). Is quantum computing our breakthrough? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
7 Plan of the presentation 1 Desperate problem, desperate solutions: the Janus computer. 2 The temperature chaos algorithmic wall. 3 A more conventional approach to temperature chaos. 4 D-wave, the chimera lattice and temperature chaos. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
8 The Janus Collaboration Team from 5 universities in Spain and Italy: Universidad Complutense de Madrid: M. Baity-Jesi, L.A. Fernandez, V. Martin-Mayor, A. Muñoz Sudupe Universidad de Extremadura: A. Gordillo-Guerrero, J.J. Ruiz-Lorenzo Università di Ferrara: M. Pivanti, S.F. Schifano, R. Tripiccione La Sapienza Università di Roma: A. Maiorano, E. Marinari, G. Parisi, F. Ricci-Tersenghi, D. Yllanes, B. Seoane Universidad de Zaragoza: R.A. Baños, A. Cruz, J.M. Gil-Narvión, M. Guidetti, D. Iñiguez, J. Monforte-Garcia, D. Navarro, S. Perez-Gaviro, A. Tarancon, P. Tellez. Physicists and engineers dedicated to the design and exploitation of special-purpose computers, optimised for Monte Carlo simulations in condensed matter physics. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
9 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
10 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. 2 Single instance simulation very long. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
11 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. 2 Single instance simulation very long. For modest system sizes (i.e. N = 32 3 = spins): Typical instance: 4.5 standard-cpu years (i.e updates). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
12 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. 2 Single instance simulation very long. For modest system sizes (i.e. N = 32 3 = spins): Typical instance: 4.5 standard-cpu years (i.e updates). Worst in 10 3 instances: 800 standard-cpu years (i.e updates). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
13 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. 2 Single instance simulation very long. For modest system sizes (i.e. N = 32 3 = spins): Typical instance: 4.5 standard-cpu years (i.e updates). Worst in 10 3 instances: 800 standard-cpu years (i.e updates). Fortunately, the spin update (the core algorithm) is very simple and (in principle) trivial to parallelize. But... V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
14 Desperate problems, desperate solutions: Janus (I) Even with binary spins, simulation of spin-glasses is heavy in two respects: 1 Many ( 10 3 ) problem instances embarrassingly parallel. 2 Single instance simulation very long. For modest system sizes (i.e. N = 32 3 = spins): Typical instance: 4.5 standard-cpu years (i.e updates). Worst in 10 3 instances: 800 standard-cpu years (i.e updates). Fortunately, the spin update (the core algorithm) is very simple and (in principle) trivial to parallelize. But... Modern architectures (GPU, Xeon, Xeon-φ) efficient only for larger N astronomical number of updates ( e cn, probably). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
15 Desperate problems, desperate solutions: Janus (II) Parallelizable problem Parallelise within each instance We divide the lattice in a checkerboard scheme, all sites of the same colour can be updated simultaneously Memory bandwith: 13 bits to update one bit! Only solution: Memory local to the processor. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
16 Desperate problems, desperate solutions: Janus (II) Parallelizable problem FPGA opportunity window: Large on-chip memory (several Mbits). Huge bandwidth on-chip distributed memory ( bits in and out per clock cycle). Large amount of logic 1024 Spin-Update Engines. Janus 1 (2008): 1000 boost in spin-glasses simulations. Green computer: energy consumption per update. Janus 2: Fall 2014 V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
17 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
18 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
19 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
20 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) Clustering correlation functions in equilibrium. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
21 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) Clustering correlation functions in equilibrium. Quantitative correspondence: equilibrium, finite size non-equilibrium, finite time. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
22 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) Clustering correlation functions in equilibrium. Quantitative correspondence: equilibrium, finite size non-equilibrium, finite time. Precision estimates of critical exponents. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
23 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) Clustering correlation functions in equilibrium. Quantitative correspondence: equilibrium, finite size non-equilibrium, finite time. Precision estimates of critical exponents. Ultrametricity in the spin-glass phase.... V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
24 Desperate problems, desperate solutions: Janus (III) Main problems tackled with Janus ( ) Spin-glasses in a magnetic field (Juan s talk, tomorrow). Isothermal aging (from 1 picosecond to 0.1 seconds) Clustering correlation functions in equilibrium. Quantitative correspondence: equilibrium, finite size non-equilibrium, finite time. Precision estimates of critical exponents. Ultrametricity in the spin-glass phase.... Temperature chaos V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
25 Temperature chaos: the showstopper (I) Increasing computing speed x1000, not such a big deal Pre-Janus era: up to N = 16 3 spins. Janus era: up to N = 32 3 spins. Why? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
26 Temperature chaos: the showstopper (I) Increasing computing speed x1000, not such a big deal Pre-Janus era: up to N = 16 3 spins. Janus era: up to N = 32 3 spins. Why? We need to talk a bit about algorithms: Simulating at fixed temperature, simply not enough. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
27 Temperature chaos: the showstopper (I) Increasing computing speed x1000, not such a big deal Pre-Janus era: up to N = 16 3 spins. Janus era: up to N = 32 3 spins. Why? We need to talk a bit about algorithms: Simulating at fixed temperature, simply not enough. Temperature needs to become dynamic. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
28 Temperature chaos: the showstopper (II) Simulated Annealing Simplest protocol: thermal energy 1 High T : easy exploration V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
29 Temperature chaos: the showstopper (II) Simulated Annealing Simplest protocol: thermal energy 1 High T : easy exploration 2 T -lowering protocol: Trapped at nearby local minimum. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
30 Temperature chaos: the showstopper (II) Simulated Annealing Simplest protocol: thermal energy 1 High T : easy exploration 2 T -lowering protocol: Trapped at nearby local minimum. Outdated algorithm. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
31 Temperature chaos: the showstopper (III) Parallel Tempering T raised or lowered: thermal energy 1 Low T : local exploration V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
32 Temperature chaos: the showstopper (III) Parallel Tempering T raised or lowered: thermal energy 1 Low T : local exploration 2 High T : global exploration V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
33 Temperature chaos: the showstopper (III) Parallel Tempering T raised or lowered: thermal energy 1 Low T : local exploration 2 High T : global exploration 3 No trapping better solution. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
34 Temperature chaos: the showstopper (III) Parallel Tempering T raised or lowered: thermal energy 1 Low T : local exploration 2 High T : global exploration 3 No trapping better solution. N T temperatures: simultaneous simulation of N T clones (one at each temperature). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
35 Temperature chaos: the showstopper (III) Parallel Tempering T raised or lowered: thermal energy 1 Low T : local exploration 2 High T : global exploration 3 No trapping better solution. N T temperatures: simultaneous simulation of N T clones (one at each temperature). Periodically, clones attempt to exchange their temperature. The rule preserves detailed balance. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
36 Temperature chaos: the showstopper (IV) It looks perfect! What can go wrong? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
37 Temperature chaos: the showstopper (IV) It looks perfect! What can go wrong? Each clone performs a temperature Random Walk. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
38 Temperature chaos: the showstopper (IV) It looks perfect! What can go wrong? Each clone performs a temperature Random Walk. The simulation is long enough if all the clones visited all the temperatures several times. Mixing time: τ. Random Walk in temperatures of a clone A mixing Random Walk A stuck Random Walk i T i T Wallclock time (hours) Wallclock time (hours) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
39 Temperature chaos: the showstopper (V) Temperature chaos: Relevant minima, completely different at nearby temperatures. T -random walk refuses to go across. T T - dt V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
40 Temperature chaos: the showstopper (VI) τ : Operational definition of Temperature chaos. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
41 Temperature chaos: the showstopper (VI) τ : Operational definition of Temperature chaos. τ-pdf L = 32, T c 1.1 Extreme sample-to-sample fluctuations. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
42 Temperature chaos: the showstopper (VI) τ : Operational definition of Temperature chaos. τ-pdf L = 32, T c 1.1 Extreme sample-to-sample fluctuations. L and T sensitivity. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
43 Temperature chaos: the showstopper (VI) τ : Operational definition of Temperature chaos. τ-pdf L = 32, T c 1.1 Extreme sample-to-sample fluctuations. L and T sensitivity. At variance with standard T -chaos studies, it is easy to observe the effect. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
44 Contacting with conventional approach (I) Defining T -chaos through a Markov MC dynamics is: Mathematically clear. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
45 Contacting with conventional approach (I) Defining T -chaos through a Markov MC dynamics is: Mathematically clear. Cumbersome (analytical computations??). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
46 Contacting with conventional approach (I) Defining T -chaos through a Markov MC dynamics is: Mathematically clear. Cumbersome (analytical computations??). Unsatisfying: T -chaos is supposed to be a static effect! V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
47 Contacting with conventional approach (I) Defining T -chaos through a Markov MC dynamics is: Mathematically clear. Cumbersome (analytical computations??). Unsatisfying: T -chaos is supposed to be a static effect! However, it provides a useful definition. Instead, the static approach: Hard for some analytically tractable models. Sherrington-Kirkpatrick: Rizzo-Crisanti (2003) Migdal-Kadanoff: McKay, Nihat-Berker, Kirkpatrick, (1982). Numerically, very hard to identify. Scaling laws barely known (Katzgraber& Krzakala, 2007). We still lack predictions relevant for experiments. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
48 Contacting with conventional approach (II) Our main ingredients (Fernandez, V.M.-M,Parisi,Seoane, 2013) Janus data base (2010): O(10 3 ) samples, L 32, well thermalized at low temperatures. Wash-out thermal fluctuations (Ney-Nifle and Young, 1997) Look at whole distribution (not only average!) Large-deviation functional (the successful analytical approach for SK) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
49 Contacting with conventional approach (II) Our main ingredients (Fernandez, V.M.-M,Parisi,Seoane, 2013) Janus data base (2010): O(10 3 ) samples, L 32, well thermalized at low temperatures. Wash-out thermal fluctuations (Ney-Nifle and Young, 1997) Look at whole distribution (not only average!) Large-deviation functional (the successful analytical approach for SK) Consistency checks: Must correlate with dynamic approach Previously subtle effects should become visible. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
50 Contacting with conventional approach (III) Washing-out thermal fluctuations Useful technicality: the chaotic parameter X J T 1,T 2 = qt 2 1,T 2 J qt 2 1,T 1 J qt 2 2,T 2 J X J = 1 no chaos; X J = 0 strong chaos. Mind that, for T 1, T 2 < T c and large L, one expects XT J 1,T 2 qt 2 1,T 2 J qt 2 1,T 1 J qt 2 2,T 2 J 1. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
51 Contacting with conventional approach (IV) Average or full distribution? X T1 =0.7026,T 2 X J T 1 =0.7026,T L = 32 L = 24 L = 16 L = 12 L = 8 average S 050 S 150 S 250 S 350 S 450 S 550 S 650 S 750 S 850 S 950 Expectations: E J (X J T 1,T 2 ) = 0 if T 1 T 2 We are far from that. T 2 > T c : nothing happens. 0.2 σ (X T1 =0.7026,T 2 ) q XT J T 2 1,T 2 = 1,T 2 J qt 2 1,T 1 J qt 2 2,T 2 J T 2 V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
52 Contacting with conventional approach (IV) Average or full distribution? 1.0 X T1 =0.7026,T 2 X J T 1 =0.7026,T L = 32 L = 24 L = 16 L = 12 L = 8 average S 050 S 150 S 250 S 350 S 450 S 550 S 650 S 750 S 850 S 950 Expectations: E J (X J T 1,T 2 ) = 0 if T 1 T 2 We are far from that. T 2 > T c : nothing happens. In the slowest samples we identify chaotic events ( level-crossings in Quantum Mechanics). σ (X T1 =0.7026,T 2 ) I J = Tmax T min dt 2 X J T min,t T 2 V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
53 Contacting with conventional approach (V) I J correlates with τ! We are on the right track... I (X T1 =0.7026,T 2 ) r = 0.766(13) averaged median 2 25 log τ exp V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
54 Contacting with conventional approach (VI) Large deviations functional 3D, T 1 = 0.7, T 2 = 0.84, (T c = 1.1). Ω L (ɛ, T 1, T 2 ) Probability[X J T 1,T 2 > ɛ] = e LD Ω L (ɛ) Ω > 0 chaos! p(x T1,T 2 ε) Ω L (ε) L = 32 L = 24 L = 16 L = 12 L = ε L = 32 L = L = 16 L = 12 L = ε 2 V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
55 Contacting with conventional approach (VI) Large deviations functional 3D, T 1 = 0.7, T 2 = 0.84, (T c = 1.1). p(x T1,T 2 ε) Ω L (ε) L = 32 L = 24 L = 16 L = 12 L = ε L = 32 L = L = 16 L = 12 L = ε 2 Ω L (ɛ, T 1, T 2 ) Probability[X J T 1,T 2 > ɛ] = e LD Ω L (ɛ) Ω > 0 chaos! Ω L approaches convergence within Janus L-range. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
56 Contacting with conventional approach (VI) Large deviations functional 3D, T 1 = 0.7, T 2 = 0.84, (T c = 1.1). p(x T1,T 2 ε) Ω L (ε) L = 32 L = 24 L = 16 L = 12 L = ε L = 32 L = L = 16 L = 12 L = ε 2 Ω L (ɛ, T 1, T 2 ) Probability[X J T 1,T 2 > ɛ] = e LD Ω L (ɛ) Ω > 0 chaos! Ω L approaches convergence within Janus L-range. L = 8, 12 chaos is extremely rare. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
57 Contacting with conventional approach (VI) Large deviations functional 3D, T 1 = 0.7, T 2 = 0.84, (T c = 1.1). p(x T1,T 2 ε) Ω L (ε) L = 32 L = 24 L = 16 L = 12 L = ε L = 32 L = L = 16 L = 12 L = ε 2 Ω L (ɛ, T 1, T 2 ) Probability[X J T 1,T 2 > ɛ] = e LD Ω L (ɛ) Ω > 0 chaos! Ω L approaches convergence within Janus L-range. L = 8, 12 chaos is extremely rare. Crude parameterization Ω T1,T 2 (ɛ) T 1 T 2 b ɛ β V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
58 Contacting with conventional approach (VI) Large deviations functional 3D, T 1 = 0.7, T 2 = 0.84, (T c = 1.1). p(x T1,T 2 ε) Ω L (ε) L = 32 L = 24 L = 16 L = 12 L = ε L = 32 L = L = 16 L = 12 L = ε 2 Ω L (ɛ, T 1, T 2 ) Probability[X J T 1,T 2 > ɛ] = e LD Ω L (ɛ) Ω > 0 chaos! Ω L approaches convergence within Janus L-range. L = 8, 12 chaos is extremely rare. Crude parameterization Ω T1,T 2 (ɛ) T 1 T 2 b ɛ β Weak-chaos scaling (Katzgraber& Krzakala, 2007) explained: ζ = D/b = 1.07(2) Chaos length: ξ C = L a unless β = 1. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
59 Contacting with conventional approach (VII) SK finally yields to numerics (Billoire 2014) Rizzo-Crisanti and Rizzo-Parisi compute Ω(q T 1,T 2 )... SK model 256 sites SK model 4096 sites V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
60 Contacting with conventional approach (VII) SK finally yields to numerics (Billoire 2014) Rizzo-Crisanti and Rizzo-Parisi compute Ω(q T 1,T 2 )... But Ω(X J T 1,T 2 > ɛ) Ω(q T 1,T 2 ) (see also Rizzo 2014) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
61 Contacting with conventional approach (VIII) Take-home messages Dynamic methods might be preferable in real calculations. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
62 Contacting with conventional approach (VIII) Take-home messages Dynamic methods might be preferable in real calculations. Three (rather than 2) scaling variables: N, T and ɛ. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
63 Contacting with conventional approach (VIII) Take-home messages Dynamic methods might be preferable in real calculations. Three (rather than 2) scaling variables: N, T and ɛ. Temperature chaos is generic for large problem size N. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
64 Contacting with conventional approach (VIII) Take-home messages Dynamic methods might be preferable in real calculations. Three (rather than 2) scaling variables: N, T and ɛ. Temperature chaos is generic for large problem size N. In practice, specially for small N: 1 The large majority of problem instances are easy (small τ). 2 For some of them, though, τ inordinately large. 3 The larger is N, the more frequently missbehaving instances appear difficult to assess algorithmic scaling with N. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
65 Is quantum-computing our breakthrough? (I) From an impressive insight (Richard P. Feynman, 1982) NP-problems, specially simulation of quantum systems: best solved on quantum computers... V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
66 Is quantum-computing our breakthrough? (I)... to (possibly) quantum-computing objects (2014). D-wave Two V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
67 Is quantum-computing our breakthrough? (I)... to (possibly) quantum-computing objects (2014). A quantum annealer should: 1 Read accurately an instance. D-wave Two 2 Add a strong transverse magnetic field. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
68 Is quantum-computing our breakthrough? (I)... to (possibly) quantum-computing objects (2014). A quantum annealer should: 1 Read accurately an instance. D-wave Two 2 Add a strong transverse magnetic field. 3 At low enough T... 4 With low noise... 5 Slowly take field 0. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
69 Is quantum-computing our breakthrough? (I)... to (possibly) quantum-computing objects (2014). A quantum annealer should: 1 Read accurately an instance. D-wave Two 2 Add a strong transverse magnetic field. 3 At low enough T... 4 With low noise... 5 Slowly take field 0. All requirements met? global minimum. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
70 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
71 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: Chimera Each blob: 8 q-bits V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
72 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: No SG phase for T > 0 T c = 0 easier problems. Chimera Each blob: 8 q-bits V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
73 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: Chimera No SG phase for T > 0 T c = 0 easier problems. Small decycling set Main graph in Selby heuristics: 78% in a single tree (no loops!) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
74 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: Chimera No SG phase for T > 0 T c = 0 easier problems. Small decycling set T = 0 heuristics better than thermal methods (i.e. PT). Main graph in Selby heuristics: 78% in a single tree (no loops!) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
75 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: Chimera No SG phase for T > 0 T c = 0 easier problems. Small decycling set T = 0 heuristics better than thermal methods (i.e. PT). T > 0 2D-methods better than Monte Carlo (Middleton 2011, square lattice: N = ). Main graph in Selby heuristics: 78% in a single tree (no loops!) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
76 Is quantum-computing our breakthrough? (II) D-wave solves a toy problem: Small problems N = 512 (actually, N = 503 in USC). Chimera graph: non-planar but 2D-like. Two-dimensional penalties: Chimera No SG phase for T > 0 T c = 0 easier problems. Small decycling set T = 0 heuristics better than thermal methods (i.e. PT). T > 0 2D-methods better than Monte Carlo (Middleton 2011, square lattice: N = ). Are we learning something? Main graph in Selby heuristics: 78% in a single tree (no loops!) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
77 Is quantum-computing our breakthrough? (III) In three spatial dimensions only thermal annealing works. Real question: Is there chaos in the chimera lattice? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
78 Is quantum-computing our breakthrough? (III) In three spatial dimensions only thermal annealing works. Real question: Is there chaos in the chimera lattice? Middleton et al.: chaos in square lattice, but N = Chaos with only N = 503 q-bits? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
79 Is quantum-computing our breakthrough? (III) In three spatial dimensions only thermal annealing works. Real question: Is there chaos in the chimera lattice? Middleton et al.: chaos in square lattice, but N = Chaos with only N = 503 q-bits? Not at first sight Probability density (instances) e+06 τ (Metropolis sweeps) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
80 Is quantum-computing our breakthrough? (III) In three spatial dimensions only thermal annealing works. Real question: Is there chaos in the chimera lattice? Middleton et al.: chaos in square lattice, but N = Chaos with only N = 503 q-bits? Not at first sight... But look at that fat tail! 2 in 10 4 instances: τ τ -0.8 Probability density (instances) e+06 τ (Metropolis sweeps) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
81 Is quantum-computing our breakthrough? (IV) As usual, τ pinpoints peculiar samples... Overlap betwen Ground-State and 1 st Excited-State... 1e GS-GS, easy sample GS-ES 1e+07 1e+06 GS-GS, hard sample GS-ES q q V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
82 Is quantum-computing our breakthrough? (IV) As usual, τ pinpoints peculiar samples or energy-dependence with T. e = e T - e GS Same easy sample Same hard sample T e = e T - e GS Same easy sample Same hard sample exp[-2/t] V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
83 Is quantum-computing our breakthrough? (V) Meaningful algorithmic classification at fixed N: τ-scaling. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
84 Is quantum-computing our breakthrough? (V) Meaningful algorithmic classification at fixed N: τ-scaling. Parallel-Tempering: τ 1, Selby heuristics (2D!): τ b 0.3, D-wave: τ a Time to solution (seconds in Xeon) e-05 τ 0.3 PT (x 0.002) Selby s average Selby s median D-wave median,τ e+06 τ (Metropolis sweeps) V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
85 Conclusions The Janus computer: new window for spin-glasses. A new generation will soon enter into operation. Not everything in Spin-Glasses Physics is self-averaging. Temperature chaos is a clear example. Sometimes, studying dynamics (τ) might be the easiest way to learn about statics. Temperature chaos is a major obstacle. Is quantum annealing an alternative? V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
86 Conclusions The Janus computer: new window for spin-glasses. A new generation will soon enter into operation. Not everything in Spin-Glasses Physics is self-averaging. Temperature chaos is a clear example. Sometimes, studying dynamics (τ) might be the easiest way to learn about statics. Temperature chaos is a major obstacle. Is quantum annealing an alternative? D-wave is a candidate quantum-annealer. Object to be experimentally investigated, rather than a finished product. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
87 Conclusions The Janus computer: new window for spin-glasses. A new generation will soon enter into operation. Not everything in Spin-Glasses Physics is self-averaging. Temperature chaos is a clear example. Sometimes, studying dynamics (τ) might be the easiest way to learn about statics. Temperature chaos is a major obstacle. Is quantum annealing an alternative? D-wave is a candidate quantum-annealer. Object to be experimentally investigated, rather than a finished product. Currently, performance not competitive with Parallel Tempering. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
88 Conclusions The Janus computer: new window for spin-glasses. A new generation will soon enter into operation. Not everything in Spin-Glasses Physics is self-averaging. Temperature chaos is a clear example. Sometimes, studying dynamics (τ) might be the easiest way to learn about statics. Temperature chaos is a major obstacle. Is quantum annealing an alternative? D-wave is a candidate quantum-annealer. Object to be experimentally investigated, rather than a finished product. Currently, performance not competitive with Parallel Tempering. Reasons for failure intrinsic? Current investigation. V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
89 Many thanks to... The Janus collaboration, Alain Billoire, Itay Hen, The meeting organizers,... and to you (the audience), for your attention! V. Martin-Mayor (Física Teórica I, UCM) Quantum vs. Classical annealing Capri, September / 28
Quantum versus Thermal annealing (or D-wave versus Janus): seeking a fair comparison
Quantum versus Thermal annealing (or D-wave versus Janus): seeking a fair comparison Víctor Martín-Mayor Dep. Física Teórica I, Universidad Complutense de Madrid Janus Collaboration In collaboration with
More informationConnected spatial correlation functions and the order parameter of the 3D Ising spin glass
Connected spatial correlation functions and the order parameter of the 3D Ising spin glass David Yllanes for the Janus Collaboration 1 Dep. Física Teórica, Universidad Complutense de Madrid http://teorica.fis.ucm.es/grupos/grupo-tec.html
More informationarxiv: v2 [cond-mat.dis-nn] 11 Oct 2016
arxiv:165.35v [cond-mat.dis-nn] 11 Oct 16 Temperature chaos is a non-local effect L. A. Fernandez 1,, E. Marinari 3,4, V. Martin-Mayor 1,, G. Parisi 3,4 and D. Yllanes 5, 1 Depto. de Física Teórica I.
More informationSpin glasses and Adiabatic Quantum Computing
Spin glasses and Adiabatic Quantum Computing A.P. Young alk at the Workshop on heory and Practice of Adiabatic Quantum Computers and Quantum Simulation, ICP, rieste, August 22-26, 2016 Spin Glasses he
More informationJanus: FPGA Based System for Scientific Computing Filippo Mantovani
Janus: FPGA Based System for Scientific Computing Filippo Mantovani Physics Department Università degli Studi di Ferrara Ferrara, 28/09/2009 Overview: 1. The physical problem: - Ising model and Spin Glass
More informationAn FPGA-based supercomputer for statistical physics: the weird case of Janus
An FPGA-based supercomputer for statistical physics: the weird case of Janus M. Baity-Jesi, R. A. Baños, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, M. Guidetti, D. Iñiguez, A. Maiorano,
More informationarxiv: v2 [cond-mat.dis-nn] 2 Nov 2011
Sample-to-sample fluctuations of the overlap distributions in the three-dimensional Edwards-Anderson spin glass arxiv:7.577v [cond-mat.dis-nn] Nov R. A. Baños, A. Cruz,, L.A. Fernandez, 3, J.M. Gil-Narvion,
More informationUniversal critical behavior of the two-dimensional Ising spin glass
PHYSICAL REVIEW B 94, 024402 (2016) Universal critical behavior of the two-dimensional Ising spin glass L. A. Fernandez, 1,2 E. Marinari, 3,4 V. Martin-Mayor, 1,2 G. Parisi, 3,4 and J. J. Ruiz-Lorenzo
More informationSpin glass phase in the four-state three-dimensional Potts model
Spin glass phase in the four-state three-dimensional Potts model A. Cruz, 1,2 L. A. Fernandez, 2,3 A. Gordillo-Guerrero, 2,4 M. Guidetti, 5 A. Maiorano, 2,5 F. Mantovani, 5 E. Marinari, 6 V. Martin-Mayor,
More informationPhase Transitions in Spin Glasses
Phase Transitions in Spin Glasses Peter Young Talk available at http://physics.ucsc.edu/ peter/talks/sinica.pdf e-mail:peter@physics.ucsc.edu Supported by the Hierarchical Systems Research Foundation.
More informationThe Last Survivor: a Spin Glass Phase in an External Magnetic Field.
The Last Survivor: a Spin Glass Phase in an External Magnetic Field. J. J. Ruiz-Lorenzo Dep. Física, Universidad de Extremadura Instituto de Biocomputación y Física de los Sistemas Complejos (UZ) http://www.eweb.unex.es/eweb/fisteor/juan
More informationNumerical Studies of the Quantum Adiabatic Algorithm
Numerical Studies of the Quantum Adiabatic Algorithm A.P. Young Work supported by Colloquium at Universität Leipzig, November 4, 2014 Collaborators: I. Hen, M. Wittmann, E. Farhi, P. Shor, D. Gosset, A.
More informationLearning About Spin Glasses Enzo Marinari (Cagliari, Italy) I thank the organizers... I am thrilled since... (Realistic) Spin Glasses: a debated, inte
Learning About Spin Glasses Enzo Marinari (Cagliari, Italy) I thank the organizers... I am thrilled since... (Realistic) Spin Glasses: a debated, interesting issue. Mainly work with: Giorgio Parisi, Federico
More informationEquilibrium Study of Spin-Glass Models: Monte Carlo Simulation and Multi-variate Analysis. Koji Hukushima
Equilibrium Study of Spin-Glass Models: Monte Carlo Simulation and Multi-variate Analysis Koji Hukushima Department of Basic Science, University of Tokyo mailto:hukusima@phys.c.u-tokyo.ac.jp July 11 15,
More informationSpin glasses, where do we stand?
Spin glasses, where do we stand? Giorgio Parisi Many progresses have recently done in spin glasses: theory, experiments, simulations and theorems! In this talk I will present: A very brief introduction
More informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationThermodynamic States in Finite Dimensional Spin Glasses
Thermodynamic States in Finite Dimensional Spin Glasses J. J. Ruiz-Lorenzo Dep. Física & ICCAEx (Univ. de Extremadura) & BIFI (Zaragoza) http://www.eweb.unex.es/eweb/fisteor/juan/juan_talks.html Instituto
More informationarxiv: v2 [cond-mat.dis-nn] 13 Dec 2017
arxiv:709.09829v2 [cond-mat.dis-nn] 3 Dec 207 Dynamic Variational Study of Chaos: Spin Glasses in Three Dimensions A. Billoire, L. A. Fernandez 2,3, A. Maiorano 4,3, E. Marinari 4,5,6, V. Martin-Mayor
More informationMind the gap Solving optimization problems with a quantum computer
Mind the gap Solving optimization problems with a quantum computer A.P. Young http://physics.ucsc.edu/~peter Work supported by NASA future technologies conference, January 17-212, 2012 Collaborators: Itay
More informationSpin glass simulations on Janus
Spin glass simulations on Janus R. (lele) Tripiccione Dipartimento di Fisica, Universita' di Ferrara raffaele.tripiccione@unife.it UCHPC, Rodos (Greece) Aug. 27 th, 2012 Warning / Disclaimer / Fineprints
More informationIs there a de Almeida-Thouless line in finite-dimensional spin glasses? (and why you should care)
Is there a de Almeida-Thouless line in finite-dimensional spin glasses? (and why you should care) Peter Young Talk at MPIPKS, September 12, 2013 Available on the web at http://physics.ucsc.edu/~peter/talks/mpipks.pdf
More informationPhase Transitions in Spin Glasses
p.1 Phase Transitions in Spin Glasses Peter Young http://physics.ucsc.edu/ peter/talks/bifi2008.pdf e-mail:peter@physics.ucsc.edu Work supported by the and the Hierarchical Systems Research Foundation.
More informationMind the gap Solving optimization problems with a quantum computer
Mind the gap Solving optimization problems with a quantum computer A.P. Young http://physics.ucsc.edu/~peter Work supported by Talk at Saarbrücken University, November 5, 2012 Collaborators: I. Hen, E.
More informationMind the gap Solving optimization problems with a quantum computer
Mind the gap Solving optimization problems with a quantum computer A.P. Young http://physics.ucsc.edu/~peter Work supported by Talk at the London Centre for Nanotechnology, October 17, 2012 Collaborators:
More informationIs the Sherrington-Kirkpatrick Model relevant for real spin glasses?
Is the Sherrington-Kirkpatrick Model relevant for real spin glasses? A. P. Young Department of Physics, University of California, Santa Cruz, California 95064 E-mail: peter@physics.ucsc.edu Abstract. I
More informationExploring reverse annealing as a tool for hybrid quantum/classical computing
Exploring reverse annealing as a tool for hybrid quantum/classical computing University of Zagreb QuantiXLie Seminar Nicholas Chancellor October 12, 2018 Talk structure 1. Background Quantum computing:
More informationThe Quantum Adiabatic Algorithm
The Quantum Adiabatic Algorithm A.P. Young http://physics.ucsc.edu/~peter Work supported by Talk at SMQS-IP2011, Jülich, October 18, 2011 The Quantum Adiabatic Algorithm A.P. Young http://physics.ucsc.edu/~peter
More informationComplexity of the quantum adiabatic algorithm
Complexity of the quantum adiabatic algorithm Peter Young e-mail:peter@physics.ucsc.edu Collaborators: S. Knysh and V. N. Smelyanskiy Colloquium at Princeton, September 24, 2009 p.1 Introduction What is
More informationSlightly off-equilibrium dynamics
Slightly off-equilibrium dynamics Giorgio Parisi Many progresses have recently done in understanding system who are slightly off-equilibrium because their approach to equilibrium is quite slow. In this
More informationarxiv: v1 [cond-mat.dis-nn] 29 Oct 2015
Composite operators in Cubic Field Theories and Link Overlap Fluctuations in Spin-Glass Models arxiv:50.08774v [cond-mat.dis-nn] 29 Oct 205 Ada Altieri,2, Giorgio Parisi,3,4 and Tommaso Rizzo 5, Dip. Fisica,
More informationNumerical Studies of Adiabatic Quantum Computation applied to Optimization and Graph Isomorphism
Numerical Studies of Adiabatic Quantum Computation applied to Optimization and Graph Isomorphism A.P. Young http://physics.ucsc.edu/~peter Work supported by Talk at AQC 2013, March 8, 2013 Collaborators:
More informationSIMULATED TEMPERING: A NEW MONTECARLO SCHEME
arxiv:hep-lat/9205018v1 22 May 1992 SIMULATED TEMPERING: A NEW MONTECARLO SCHEME Enzo MARINARI (a),(b) and Giorgio PARISI (c) Dipartimento di Fisica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,
More informationCURRICULUM VITAE ET STUDIORUM
CURRICULUM VITAE ET STUDIORUM Surname and Name: Ricci Tersenghi, Federico Place and Date of Birth: Roma (Italy). July 9, 1972 Nationality: Italian Current Position: Associate Professor in Theoretical Physics
More informationMonte Carlo Methods for Rough Free Energy Landscapes: Population Annealing and Parallel Tempering. Abstract
Monte Carlo Methods for Rough Free Energy Landscapes: Population Annealing and Parallel Tempering J. Machta Physics Department, University of Massachusetts, Amherst, MA 01003 USA R. S. Ellis Department
More informationLow-temperature behavior of the statistics of the overlap distribution in Ising spin-glass models
University of Massachusetts Amherst From the SelectedWorks of Jonathan Machta 2014 Low-temperature behavior of the statistics of the overlap distribution in Ising spin-glass models Matthew Wittmann B.
More information[7] M. Falcioni, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti, Phys. Lett. B 108
[5] G. Parisi, Statistical Field Theory, Addisson Wesley 1992. [6] U. Wol, Phys. Lett. 228B 3(1989) [7] M. Falcioni, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti, Phys. Lett. B 108 (1982) 331.
More informationNeural Networks for Machine Learning. Lecture 11a Hopfield Nets
Neural Networks for Machine Learning Lecture 11a Hopfield Nets Geoffrey Hinton Nitish Srivastava, Kevin Swersky Tijmen Tieleman Abdel-rahman Mohamed Hopfield Nets A Hopfield net is composed of binary threshold
More informationarxiv: v1 [cond-mat.dis-nn] 15 Dec 2018
arxiv:1812.6334v1 [cond-mat.dis-nn] 15 Dec 218 Study of longitudinal fluctuations of the Sherrington-Kirkpatrick model 1. Introduction Giorgio Parisi 1, Leopoldo Sarra 2, Lorenzo Talamanca 2 1 Dipartimento
More informationWeak first-order transition in the three-dimensional site-diluted Ising antiferromagnet in a magnetic field
Weak first-order transition in the three-dimensional site-diluted Ising antiferromagnet in a magnetic field A. Maiorano,,2 V. Martín-Mayor, 3,2 J. J. Ruiz-Lorenzo,,2 and A. Tarancón 4,2 Departamento de
More informationLista delle Pubblicazioni di Vincenzo Marinari
Lista delle Pubblicazioni di Vincenzo Marinari 1. Computation of the Critical Exponents of Percolation, F. Fucito and E. Marinari, J. Phys. A: Math. Gen. 14 (1981) L85; 2. A Proposal for Monte Carlo Simulations
More informationLow T scaling behavior of 2D disordered and frustrated models
Low T scaling behavior of 2D disordered and frustrated models Enzo Marinari (Roma La Sapienza, Italy) A. Galluccio, J. Lukic, E.M., O. C. Martin and G. Rinaldi, Phys. Rev. Lett. 92 (2004) 117202. Additional
More informationPhase Transitions (and their meaning) in Random Constraint Satisfaction Problems
International Workshop on Statistical-Mechanical Informatics 2007 Kyoto, September 17 Phase Transitions (and their meaning) in Random Constraint Satisfaction Problems Florent Krzakala In collaboration
More informationarxiv: v1 [cond-mat.dis-nn] 23 Apr 2008
lack of monotonicity in spin glass correlation functions arxiv:0804.3691v1 [cond-mat.dis-nn] 23 Apr 2008 Pierluigi Contucci, Francesco Unguendoli, Cecilia Vernia, Dipartimento di Matematica, Università
More informationQuantum and classical annealing in spin glasses and quantum computing. Anders W Sandvik, Boston University
NATIONAL TAIWAN UNIVERSITY, COLLOQUIUM, MARCH 10, 2015 Quantum and classical annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Cheng-Wei Liu (BU) Anatoli Polkovnikov (BU)
More informationQuantum vs Classical Optimization: A status update on the arms race. Helmut G. Katzgraber
Quantum vs Classical Optimization: A status update on the arms race Helmut G. Katzgraber https://intractable.lol Quantum vs Classical Optimization: A status update on the arms race Helmut G. Katzgraber
More informationQuantum annealing. Matthias Troyer (ETH Zürich) John Martinis (UCSB) Dave Wecker (Microsoft)
Quantum annealing (ETH Zürich) John Martinis (UCSB) Dave Wecker (Microsoft) Troels Rønnow (ETH) Sergei Isakov (ETH Google) Lei Wang (ETH) Sergio Boixo (USC Google) Daniel Lidar (USC) Zhihui Wang (USC)
More informationarxiv:cond-mat/ v2 [cond-mat.dis-nn] 14 May 1998
1 arxiv:cond-mat/9805125v2 [cond-mat.dis-nn] 14 May 1998 Scaling corrections: site-percolation and Ising model in three dimensions H. G. Ballesteros a, L. A. Fernández a, V. Martín-Mayor a, A. Muñoz Sudupe
More informationMarkov Chain Monte Carlo The Metropolis-Hastings Algorithm
Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Anthony Trubiano April 11th, 2018 1 Introduction Markov Chain Monte Carlo (MCMC) methods are a class of algorithms for sampling from a probability
More information27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
More informationPopulation Annealing: An Effective Monte Carlo Method for Rough Free Energy Landscapes
Population Annealing: An Effective Monte Carlo Method for Rough Free Energy Landscapes Jon Machta SFI WORKING PAPER: 21-7-12 SFI Working Papers contain accounts of scientific work of the author(s) and
More informationOptimization Methods via Simulation
Optimization Methods via Simulation Optimization problems are very important in science, engineering, industry,. Examples: Traveling salesman problem Circuit-board design Car-Parrinello ab initio MD Protein
More informationThere are self-avoiding walks of steps on Z 3
There are 7 10 26 018 276 self-avoiding walks of 38 797 311 steps on Z 3 Nathan Clisby MASCOS, The University of Melbourne Institut für Theoretische Physik Universität Leipzig November 9, 2012 1 / 37 Outline
More informationPhysics 115/242 Monte Carlo simulations in Statistical Physics
Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly
More informationCorrelations Between the Dynamics of Parallel Tempering and the Free Energy Landscape in Spin Glasses
University of Massachusetts Amherst From the SelectedWorks of Jonathan Machta January, 2013 Correlations Between the Dynamics of Parallel Tempering and the Free Energy Landscape in Spin Glasses Burcu Yucesoy,
More informationMarkov Chain Monte Carlo Method
Markov Chain Monte Carlo Method Macoto Kikuchi Cybermedia Center, Osaka University 6th July 2017 Thermal Simulations 1 Why temperature 2 Statistical mechanics in a nutshell 3 Temperature in computers 4
More informationProblem Set 2. MAS 622J/1.126J: Pattern Recognition and Analysis. Due: 5:00 p.m. on September 30
Problem Set 2 MAS 622J/1.126J: Pattern Recognition and Analysis Due: 5:00 p.m. on September 30 [Note: All instructions to plot data or write a program should be carried out using Matlab. In order to maintain
More informationEnergy-Decreasing Dynamics in Mean-Field Spin Models
arxiv:cond-mat/0210545 v1 24 Oct 2002 Energy-Decreasing Dynamics in Mean-Field Spin Models L. Bussolari, P. Contucci, M. Degli Esposti, C. Giardinà Dipartimento di Matematica dell Università di Bologna,
More informationPopulation annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice
Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice Michal Borovský Department of Theoretical Physics and Astrophysics, University of P. J. Šafárik in Košice,
More informationA Tour of Spin Glasses and Their Geometry
A Tour of Spin Glasses and Their Geometry P. Le Doussal, D. Bernard, LPTENS A. Alan Middleton, Syracuse University Support from NSF, ANR IPAM: Random Shapes - Workshop I 26 March, 2007 Goals Real experiment
More information12. LOCAL SEARCH. gradient descent Metropolis algorithm Hopfield neural networks maximum cut Nash equilibria
12. LOCAL SEARCH gradient descent Metropolis algorithm Hopfield neural networks maximum cut Nash equilibria Lecture slides by Kevin Wayne Copyright 2005 Pearson-Addison Wesley h ttp://www.cs.princeton.edu/~wayne/kleinberg-tardos
More informationTheir Statistical Analvsis. With Web-Based Fortran Code. Berg
Markov Chain Monter rlo Simulations and Their Statistical Analvsis With Web-Based Fortran Code Bernd A. Berg Florida State Univeisitfi USA World Scientific NEW JERSEY + LONDON " SINGAPORE " BEIJING " SHANGHAI
More informationand B. Taglienti (b) (a): Dipartimento di Fisica and Infn, Universita di Cagliari (c): Dipartimento di Fisica and Infn, Universita di Roma La Sapienza
Glue Ball Masses and the Chameleon Gauge E. Marinari (a),m.l.paciello (b),g.parisi (c) and B. Taglienti (b) (a): Dipartimento di Fisica and Infn, Universita di Cagliari Via Ospedale 72, 09100 Cagliari
More informationA Search and Jump Algorithm for Markov Chain Monte Carlo Sampling. Christopher Jennison. Adriana Ibrahim. Seminar at University of Kuwait
A Search and Jump Algorithm for Markov Chain Monte Carlo Sampling Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Adriana Ibrahim Institute
More informationModern WalkSAT algorithms
Modern WalkSAT algorithms J. Rosti Summary of workshop presentation Special research course T-79.7003, Autumn 2007 Phase Transitions in Optimization Problems ABSTRACT In this workshop presentation summary
More informationSimulated Quantum Annealing For General Ising Models
Simulated Quantum Annealing For General Ising Models Thomas Neuhaus Jülich Supercomputing Centre, JSC Forschungszentrum Jülich Jülich, Germany e-mail : t.neuhaus@fz-juelich.de November 23 On the Talk Quantum
More informationResults on top physics by CMS
EPJ Web of Conferences 95, 04069 (2015) DOI: 10.1051/ epjconf/ 20159504069 C Owned by the authors, published by EDP Sciences, 2015 Results on top physics by CMS Silvano Tosi 1,2,a, on behalf of the CMS
More informationQuantum Computing. Separating the 'hope' from the 'hype' Suzanne Gildert (D-Wave Systems, Inc) 4th September :00am PST, Teleplace
Quantum Computing Separating the 'hope' from the 'hype' Suzanne Gildert (D-Wave Systems, Inc) 4th September 2010 10:00am PST, Teleplace The Hope All computing is constrained by the laws of Physics and
More informationSimulated Tempering: a New Monte Carlo Scheme.
EUROPHYSICS LETTERS Europhys. Lett., 19 (6), pp. 451458 (1992) 15 July 1992 Simulated Tempering: a New Monte Carlo Scheme. E. MARINARI (*) (**)(9 and G. PARISI (*)(@) (*> Dipartimento di Fisica, Universitd
More informationTrivially parallel computing
Parallel Computing After briefly discussing the often neglected, but in praxis frequently encountered, issue of trivially parallel computing, we turn to parallel computing with information exchange. Our
More informationParallel tempering and 3D spin glass models
Journal of Physics: Conference Series OPEN ACCESS Parallel tempering and 3D spin glass models To cite this article: T Papakonstantinou and A Malakis 2014 J. Phys.: Conf. Ser. 487 012010 View the article
More informationMulticanonical methods
Multicanonical methods Normal Monte Carlo algorithms sample configurations with the Boltzmann weight p exp( βe). Sometimes this is not desirable. Example: if a system has a first order phase transitions
More informationImprovement of Monte Carlo estimates with covariance-optimized finite-size scaling at fixed phenomenological coupling
Improvement of Monte Carlo estimates with covariance-optimized finite-size scaling at fixed phenomenological coupling Francesco Parisen Toldin Max Planck Institute for Physics of Complex Systems Dresden
More informationarxiv:hep-lat/ v2 8 Jan 2003
Dynamical generation of a gauge symmetry in the Double-Exchange model arxiv:hep-lat/0301004v2 8 Jan 2003 J.M. Carmona a, A. Cruz a,e, L. A. Fernández b,e, S. Jiménez a,e, V. Martín-Mayor b,e, A. Muñoz-Sudupe
More informationMarkov Chains for Tensor Network States
Sofyan Iblisdir University of Barcelona September 2013 arxiv 1309.4880 Outline Main issue discussed here: how to find good tensor network states? The optimal state problem Markov chains Multiplicative
More informationThe cavity method. Vingt ans après
The cavity method Vingt ans après Les Houches lectures 1982 Early days with Giorgio SK model E = J ij s i s j i
More informationPhase Transition in Vector Spin Glasses. Abstract
Phase Transition in Vector Spin Glasses A. P. Young Department of Physics, University of California, Santa Cruz, California 95064 (Dated: June 3, 2004) Abstract We first give an experimental and theoretical
More informationQuadratic replica coupling in the Sherrington-Kirkpatrick mean field spin glass model
arxiv:cond-mat/0201091v2 [cond-mat.dis-nn] 5 Mar 2002 Quadratic replica coupling in the Sherrington-Kirkpatrick mean field spin glass model Francesco Guerra Dipartimento di Fisica, Università di Roma La
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationExperiments with and Applications of the D-Wave Machine
Experiments with and Applications of the D-Wave Machine Paul Warburton, University College London p.warburton@ucl.ac.uk 1. Brief introduction to the D-Wave machine 2. Black box experiments to test quantumness
More informationQuantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC
Quantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC Itay Hen Information Sciences Institute, USC Workshop on Theory and Practice of AQC and Quantum Simulation Trieste, Italy August
More informationThe ultrametric tree of states and computation of correlation functions in spin glasses. Andrea Lucarelli
Università degli studi di Roma La Sapienza Facoltà di Scienze Matematiche, Fisiche e Naturali Scuola di Dottorato Vito Volterra Prof. Giorgio Parisi The ultrametric tree of states and computation of correlation
More informationA. Time dependence from separation of timescales
Lecture 4 Adiabatic Theorem So far we have considered time independent semiclassical problems. What makes these semiclassical is that the gradient term (which is multiplied by 2 ) was small. In other words,
More informationHill climbing: Simulated annealing and Tabu search
Hill climbing: Simulated annealing and Tabu search Heuristic algorithms Giovanni Righini University of Milan Department of Computer Science (Crema) Hill climbing Instead of repeating local search, it is
More informationFilippo Tramonto. Miniworkshop talk: Quantum Monte Carlo simula9ons of low temperature many- body systems
Miniworkshop talk: Quantum Monte Carlo simulations of low temperature many-body systems Physics, Astrophysics and Applied Physics Phd school Supervisor: Dott. Davide E. Galli Outline Interests in quantum
More informationStatistics and Quantum Computing
Statistics and Quantum Computing Yazhen Wang Department of Statistics University of Wisconsin-Madison http://www.stat.wisc.edu/ yzwang Workshop on Quantum Computing and Its Application George Washington
More informationExtremal Optimization for Low-Energy Excitations of very large Spin-Glasses
Etremal Optimization for Low-Energy Ecitations of very large Spin-Glasses Stefan Boettcher Dept of Physics Collaborators: Senior: Allon Percus (Los Alamos) Paolo Sibani (Odense-DK) Michelangelo Grigni
More informationPhase Transition in Vector Spin Glasses. Abstract
Phase Transition in Vector Spin Glasses A. P. Young Department of Physics, University of California, Santa Cruz, California 95064 (Dated: June 2, 2004) Abstract We first give an experimental and theoretical
More informationPhysics Sep Example A Spin System
Physics 30 7-Sep-004 4- Example A Spin System In the last lecture, we discussed the binomial distribution. Now, I would like to add a little physical content by considering a spin system. Actually this
More informationCritical properties of disordered XY model on sparse random graphs
Manuale di Identità Visiva Sapienza Università di Roma Sapienza Università di Roma Scuola di Dottorato Vito Volterra Dipartimento di Fisica Critical properties of disordered XY model on sparse random graphs
More informationThe Markov Chain Monte Carlo Method
The Markov Chain Monte Carlo Method Idea: define an ergodic Markov chain whose stationary distribution is the desired probability distribution. Let X 0, X 1, X 2,..., X n be the run of the chain. The Markov
More informationNew insights from one-dimensional spin glasses.
New insights from one-dimensional spin glasses http://katzgraber.org/stuff/leipzig 1 New insights from one-dimensional spin glasses Helmut G. Katzgraber http://katzgraber.org 2 New insights from one-dimensional
More information7.1 Coupling from the Past
Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov
More informationMonte Carlo methods for sampling-based Stochastic Optimization
Monte Carlo methods for sampling-based Stochastic Optimization Gersende FORT LTCI CNRS & Telecom ParisTech Paris, France Joint works with B. Jourdain, T. Lelièvre, G. Stoltz from ENPC and E. Kuhn from
More informationGPU-based computation of the Monte Carlo simulation of classical spin systems
Perspectives of GPU Computing in Physics and Astrophysics, Sapienza University of Rome, Rome, Italy, September 15-17, 2014 GPU-based computation of the Monte Carlo simulation of classical spin systems
More informationMonte Carlo Methods. PHY 688: Numerical Methods for (Astro)Physics
Monte Carlo Methods Random Numbers How random is random? From Pang: we want Long period before sequence repeats Little correlation between numbers (plot ri+1 vs ri should fill the plane) Fast Typical random
More informationLOCKHEED MARTIN SITE UPDATE
LOCKHEED MARTIN SITE UPDATE 25 SEPTEMBER 2018 Julia Kwok Software Engineer Quantum Applications THE USC-LM QUANTUM COMPUTING CENTER Dr. Edward H. Ned Allen Chief Scientist and LM Senior Fellow Lockheed
More informationBranislav K. Nikolić
Interdisciplinary Topics in Complex Systems: Cellular Automata, Self-Organized Criticality, Neural Networks and Spin Glasses Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware,
More informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 3 More Markov Chain Monte Carlo Methods The Metropolis algorithm isn t the only way to do MCMC. We ll
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More information