Some recent results on the inverse Ising problem. Federico Ricci-Tersenghi Physics Department Sapienza University, Roma
|
|
- Paulina Merritt
- 5 years ago
- Views:
Transcription
1 Some recent results on the inverse Ising problem Federico Ricci-Tersenghi Physics Department Sapienza University, Roma
2 Outline of the talk Bethe approx. for inverse Ising problem Comparison among several mean-field approximations (IP, nmf, TAP, SM, Bethe...) Correlations normalization trick Major limitation in using TAP and Bethe in frustrated models with a field New method: pseudo-likelihood + decimation (Aurelien DECELLE)
3 Inverse Ising problem M configurations of N Ising variables ( s i = ±1) extracted from P (s 1,...,s N )= 1 Z(J, h) exp [alternatively only magnetizations correlations i=j C ij = s i s j m i m j J ij s i s j + i and are given] m i = s i h i s i GOAL: estimate couplings and fields (J, h)
4 Solving inverse Ising pb. Monte Carlo -> unbiased solution...but it is slow! Initial guess for J(0),h(0) Given J(t),h(t) MC compute m(t),c(t) Return J(t),h(t) Yes No dm<eps? dc<eps? Update J(t+1),h(t+1) according to dm(t)=m(t)-m dc(t)=c(t)-c
5 Solving inverse Ising pb. Monte Carlo -> unbiased solution...but it is slow! Initial guess for J(0),h(0) Given J(t),h(t) MC compute m(t),c(t) Return J(t),h(t) Yes No dm<eps? dc<eps? Update J(t+1),h(t+1) according to dm(t)=m(t)-m dc(t)=c(t)-c
6 Solving inverse Ising pb. Mean field approximations (MFA) m MFA i (J, h) =m i Cij MFA (J, h) =C ij
7 Solving inverse Ising pb. Mean field approximations (MFA) m MFA i (J, h) =m i Cij MFA (J, h) =C ij
8 Solving inverse Ising pb. Mean field approximations (MFA) m MFA i (J, h) =m i Cij MFA (J, h) =C ij Log-likelihood L(J, h s) = 1 M log M P (s (k) J, h) = = i h i m i + ij k=1 J ij (C ij + m i m j ) log Z(J, h)
9 Solving inverse Ising pb. Mean field approximations (MFA) m MFA i (J, h) =m i Cij MFA (J, h) =C ij Log-likelihood L(J, h s) = 1 M log M P (s (k) J, h) = = i h i m i + ij k=1 J ij (C ij + m i m j ) log Z(J, h) F MFA (J, h)
10 Solving inverse Ising pb. Mean field approximations (MFA) m MFA i (J, h) =m i Cij MFA (J, h) =C ij Log-likelihood L(J, h s) = 1 M log M = i h i m i + ij k=1 P (s (k) J, h) = J ij (C ij + m i m j ) log Z(J, h) m i = F MFA (J, h)/ h i F MFA (J, h)
11 Mean-field approximations (MFA) to the free-energy naive mean-field (nmf) F nmf = i 1+mi 1 mi H + H i h i m i + i=j J ij m i m j H(x) x ln(x) F nmf m i = j J ij m j + h i atanh(m i )=0 m i = tanh h i + j J ij m j
12 MFA to the free-energy nmf + Onsager reaction term (TAP) F TAP = 1+mi 1 mi H + H i + h i m i + J ij m i m j J ij(1 2 m 2 i )(1 m 2 j) i i=j m i = tanh h i + J ij m j J ij (1 m 2 j)m i j
13 MFA to the free-energy nmf + Onsager reaction term (TAP) F TAP = 1+mi 1 mi H + H i + h i m i + J ij m i m j J ij(1 2 m 2 i )(1 m 2 j) i i=j m i = tanh h i + J ij m j J ij (1 m 2 j)m i j reaction term
14 MFA to the free-energy Plefka expansion in small J F nmf = i H 1+mi + H 2 1 mi 2 + i h i m i + J ij m i m j i=j F TAP = 1+mi 1 mi H + H i + h i m i + J ij m i m j J ij(1 2 m 2 i )(1 m 2 j) i i=j
15 MFA to the free-energy Bethe approximation (BA) F BA = (1 + mi )(1 + m j )+c ij (1 mi )(1 m j )+c ij H + H i=j (1 + mi )(1 m j ) c ij (1 mi )(1 + m j ) c ij + H + H mi 1 mi (1 d i ) H + H + h i m i + J ij (c ij + m i m j ), 2 2 i i i=j
16 MFA to the free-energy Bethe approximation (BA) F BA = (1 + mi )(1 + m j )+c ij (1 mi )(1 m j )+c ij H + H i=j (1 + mi )(1 m j ) c ij (1 mi )(1 + m j ) c ij + H + H mi 1 mi (1 d i ) H + H + h i m i + J ij (c ij + m i m j ), 2 2 i i i=j F BA / c ij = 0 c ij (m i,m j,t ij ) = 1 2t ij 1+t 2 ij t ij = tanh(j ij ) (1 t 2 ij )2 4t ij (m i t ij m j )(m j t ij m i ) m i m j
17 MFA to the free-energy Bethe approximation (BA) and cavity method m i = m(j) i 1+m (j) i m j = t ij m (j) i 1+m (j) i + t ij m (i) j t ij m (i) j + m (i) j t ij m (i) j i m (j) i t ij : magnetization of i in absence of j j
18 MFA to the free-energy Bethe approximation (BA) and cavity method m i = m(j) i 1+m (j) i m j = t ij m (j) i 1+m (j) i + t ij m (i) j t ij m (i) j + m (i) j t ij m (i) j i m (j) i t ij : magnetization of i in absence of j j m (j) i = f(m i,m j,t ij ) m (i) j = f(m j,m i,t ij ) f(m 1,m 2,t)= 1 t2 (1 t 2 ) 2 4t(m 1 m 2 t)(m 2 m 1 t) 2t(m 2 m 1 t)
19 MFA to the free-energy Bethe approximation (BA) and cavity method f(m 1,m 2,t)= 1 t2 (1 t 2 ) 2 4t(m 1 m 2 t)(m 2 m 1 t) 2t(m 2 m 1 t) m i = tanh h i + j ) atanh t ij f(m j,m i,t ij
20 MFA to the free-energy Bethe approximation (BA) and cavity method f(m 1,m 2,t)= 1 t2 (1 t 2 ) 2 4t(m 1 m 2 t)(m 2 m 1 t) 2t(m 2 m 1 t) m i = tanh h i + j ) atanh t ij f(m j,m i,t ij Small J expansion gives nmf, TAP,... h i + j atanh t ij f(m j,m i,t ij ) h i + j J ij m j Jij(1 2 m 2 j)m i +...
21 Computing correlations by linear response Correlations are trivial in MFA C ij =0 in nmf, TAP and BA (between distant spins) Non trivial correlations can be obtained by using the linear response (Kappen Rodriguez, 1998) C ij = m i h j, (C 1 ) ij = h i m j
22 Computing correlations by linear response Correlations are trivial in MFA C ij =0 in nmf, TAP and BA (between distant spins) Non trivial correlations can be obtained by using the linear response (Kappen Rodriguez, 1998) C ij = m i h j, (C 1 ) ij = h i m j (C 1 nmf ) ij = (C 1 TAP ) ij = δ ij 1 m 2 J ij, i 1 1 m 2 + i k J 2 ik (1 m2 k ) δ ij J ij +2J 2 ijm i m j
23 Computing correlations by linear response in BA Analytic expression for the correlations (C 1 BA ) ij = 1 1 m 2 i k t ik f 2 (m k,m i,t ik ) 1 t 2 ik f(m k,m i,t ik ) 2 Coincide with the fixed point of Susceptibility Propagation No need to run any algorithm! δ ij t ijf 1 (m j,m i,t ij ) 1 t 2 ij f(m j,m i,t ij ) 2
24 Numerical results on estimating correlations C 1 N 2 (C ij Cij )2 1e+02 1 i,j 1e+00 1 C 1e-02 1e-04 1e-06 2D ferromagnet diluted (p=0.7) nmf TAP 3 rd 4 th BA 1e β
25 Improving correlations by a normalization trick 2D lattice In ferromagnetic models with loops, inferred correlations are too strong because of loops Usually C MFA ii > 1 Trick: enforce which is unphysical C MFA ii =1 by a normalization C ij C ij Cii C jj
26 2 More numerical results on estimating correlations 1e+02 1e MC data 10 6 MC data 1e-02 C 2D ferromagnet diluted (p=0.7) 1e-04 exact corr 1e-06 1e β TAP TAP norm BA BA norm
27 MFA for the inverse Ising problem
28 MFA for the inverse Ising problem (C 1 nmf ) ij = (C 1 TAP ) ij = δ ij 1 m 2 J ij, i 1 1 m 2 + i k J 2 ik (1 m2 k ) δ ij J ij +2J 2 ijm i m j
29 MFA for the inverse Ising problem (C 1 nmf ) ij = (C 1 TAP ) ij = δ ij 1 m 2 J ij, i 1 1 m 2 + i k J 2 ik (1 m2 k ) δ ij J ij +2J 2 ijm i m j
30 MFA for the inverse Ising problem (C 1 nmf ) ij = (C 1 TAP ) ij = δ ij 1 m 2 J ij, i 1 1 m 2 + i k = J nmf ij J 2 ik (1 m2 k ) = (C 1 ) ij δ ij J ij +2J 2 ijm i m j
31 MFA for the inverse Ising problem (C 1 nmf ) ij = (C 1 TAP ) ij = J TAP ij = δ ij = J nmf = (C 1 ) ij 1 m 2 J ij, ij i 1 1 m 2 + Jik 2 (1 m2 k ) δ ij J ij +2Jijm 2 i m j i k 1 8mi m j (C 1 ) ij 1 4m i m j
32 MFA for the inverse Ising problem J BA ij (C 1 nmf ) ij = (C 1 TAP ) ij = 2(C 1 ) ij J TAP ij = δ ij = J nmf = (C 1 ) ij 1 m 2 J ij, ij i 1 1 m 2 + Jik 2 (1 m2 k ) δ ij J ij +2Jijm 2 i m j i k 1 8mi m j (C 1 ) ij 1 4m i m j 1 = atanh 2(C (1 m 2 i ) )(1 m2 j )(C 1 ) 2 ij m im j ij (1 m 2 i )(1 m2 j )(C 1 ) 2 ij 2m im j (C 1 ) ij 4(C 1 ) 2 ij
33 Normalization trick for the inverse Ising problem TAP with m=0 J TAP ij = (C 1 ) ij
34 Normalization trick for the inverse Ising problem TAP with m=0 (C 1 TAP ) ij = 1+ k J TAP ij = (C 1 ) ij (Jik TAP ) 2 δ ij J TAP ij = C TAP = C
35 Normalization trick for the inverse Ising problem TAP with m=0 (C 1 TAP ) ij = 1+ k J TAP ij = (C 1 ) ij (Jik TAP ) 2 δ ij J TAP ij = C TAP = C ( C TAP ) ij = (C TAP ) ij (CTAP ) ii (C TAP ) jj = λ i λ j (C TAP ) ij λ i are N new parameters to fix s.t. CTAP = C
36 More approximations for the inverse Ising problem Independent pair (IP) approximation J IP ij = 1 4 ln (1 + m i )(1 + m j )+C ij (1 m i )(1 m j )+C ij (1 + m i )(1 m j ) C ij (1 m i )(1 + m j ) C ij Sessak-Monasson (SM) small correlation expansion J SM ij = (C 1 ) ij + J IP ij C ij (1 m 2 i )(1 m2 j ) (C ij) 2 [for m=0 can be derived in 2 lines, ask me if interested]
37 More approximations for the inverse Ising problem Independent pair (IP) approximation J IP ij = 1 4 ln (1 + m i )(1 + m j )+C ij (1 m i )(1 m j )+C ij (1 + m i )(1 m j ) C ij (1 m i )(1 + m j ) C ij Sessak-Monasson (SM) small correlation expansion J SM ij = (C 1 ) ij + J IP ij C ij (1 m 2 i )(1 m2 j ) (C ij) 2 [for m=0 can be derived in 2 lines, ask me if interested]
38 Numerical results for the inverse Ising problem J = i<j (J ij J ij) 2 i<j J 2 ij 10 1 ferromagnet N=20 on Bethe lattice (c=4) 0.1 J M = 10 4 M = 10 5 M = 10 6 exact BA SM β
39 Numerical results for the inverse Ising problem 10 2D ferromagnet N=7 2 diluted (p=0.7) 1 J β IP TAP SM BA BA norm
40 A major problem for frustrated models in a field J TAP ij = 1 8mi m j (C 1 ) ij 1 4m i m j The resulting coupling is real only if TAP 1 8m i m j (C 1 ) ij 0 The same happens also within the BA BA = 1 + 4(1 m 2 i )(1 m2 j )(C 1 ) 2 ij 2m im j (C 1 ) ij 2 4(C 1 ) 2 ij
41 A major problem for frustrated models in a field Three spins interacting with a coupling J in a constant field h P (s 1,s 2,s 3 ) exp[j(s 1 s 2 + s 2 s 3 + s 3 s 1 )+h(s 1 + s 2 + s 3 )] 5 TAP BA 4 < h 1 > J
42 Pseudo-likelihood method (PLM) Approximate P i (s i s \i ) P (s) i Define the pseudo-log-likelihood as PL(h, J s) =log P (s h, J) = i log P i (s h i, J i ) = i f i f i = h i m i + j J ij (C ij + m i m j ) log 2 cosh s i (h i + j J ij s j ) Maximize to estimate and For sparse models use an L1 regularization and maximize f i h i J ij f i λ j J ij
43 PLM vs. MFA 10 2D ferromagnet N=7 2 =49 diluted (p=0.7) M=5000 samples J β IP TAP SM BA PLM BA norm
44 Couplings inferred by PLM ferro 2D N=7 2 diluted (p=0.7) β=0.6 PLM λ=
45 Couplings inferred by PLM ferro 2D N=7 2 diluted (p=0.7) β=0.6 J=0 J=1 PLM λ=
46 Couplings inferred by PLM ferro 2D N=7 2 diluted (p=0.7) β=0.9 PLM λ=
47 Couplings inferred by PLM ferro 2D N=7 2 diluted (p=0.7) β=0.9 J=0 J=1 PLM λ=
48 Decimation procedure Run PLM Set to zero a constant fraction of couplings (those inferred to be the smallest) Re-run PLM only on couplings not set to zero (this is impossible within a MFA) Iterate until...
49 PLM + decimation β = J β = β = 0.6 β = num. iterations
50 PLM + decimation β = J β = β = 0.6 β = num. iterations
51 PLM + decimation PL β = 0.9 β = 0.8 β = 0.6 β = num. iterations
52 PLM + decimation PL log(2) β=0.9 β=0.8 β=0.6 β= #DOF
53 PLM + decimation β=0.9 β=0.8 β=0.6 β= tilted PL log(2) #DOF
54 Summary of recent results about the inverse Ising pb. Analytical expressions for the Bethe approx. Comparison of MFA in a wide temperature range Serious limitation of TAP and Bethe in a field Improvement in the inferred couplings by: normalization trick (in weakly frustrated models) PLM + decimation procedure JSTAT (2012) P08015
A brief introduction to the inverse Ising problem and some algorithms to solve it
A brief introduction to the inverse Ising problem and some algorithms to solve it Federico Ricci-Tersenghi Physics Department Sapienza University, Roma original results in collaboration with Jack Raymond
More informationInference in kinetic Ising models: mean field and Bayes estimators
Inference in kinetic Ising models: mean field and Bayes estimators Ludovica Bachschmid-Romano, Manfred Opper Artificial Intelligence group, Computer Science, TU Berlin, Germany October 23, 2015 L. Bachschmid-Romano,
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 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 informationREVIEW: Derivation of the Mean Field Result
Lecture 18: Mean Field and Metropolis Ising Model Solutions 1 REVIEW: Derivation of the Mean Field Result The Critical Temperature Dependence The Mean Field Approximation assumes that there will be an
More informationBelief Propagation for Traffic forecasting
Belief Propagation for Traffic forecasting Cyril Furtlehner (INRIA Saclay - Tao team) context : Travesti project http ://travesti.gforge.inria.fr/) Anne Auger (INRIA Saclay) Dimo Brockhoff (INRIA Lille)
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 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 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 informationUNDERSTANDING BELIEF PROPOGATION AND ITS GENERALIZATIONS
UNDERSTANDING BELIEF PROPOGATION AND ITS GENERALIZATIONS JONATHAN YEDIDIA, WILLIAM FREEMAN, YAIR WEISS 2001 MERL TECH REPORT Kristin Branson and Ian Fasel June 11, 2003 1. Inference Inference problems
More informationarxiv: v3 [cond-mat.dis-nn] 3 Sep 2012
Replica Cluster Variational Method: the Replica Symmetric solution for the 2D random bond Ising model Alejandro Lage-Castellanos Henri-Poincaré-Group of Complex Systems and Department of Theoretical Physics,
More informationNETADIS Kick off meeting. Torino,3-6 february Payal Tyagi
NETADIS Kick off meeting Torino,3-6 february 2013 Payal Tyagi Outline My background Project: 1) Theoretical part a) Haus master equation and Hamiltonian b) First steps: - Modelling systems with disordered
More informationEfficient inference of interactions from non-equilibrium data and application to multi-electrode neural recordings
Efficient inference of interactions from non-equilibrium data and application to multi-electrode neural recordings ESR: 1 Supervisor: Dr. Yasser Roudi 1, 2 1 Kavli Institute for Systems Neuroscience, Trondheim
More informationStatistical Thermodynamics Solution Exercise 8 HS Solution Exercise 8
Statistical Thermodynamics Solution Exercise 8 HS 05 Solution Exercise 8 Problem : Paramagnetism - Brillouin function a According to the equation for the energy of a magnetic dipole in an external magnetic
More informationQuantum spin systems - models and computational methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Quantum spin systems - models and computational methods Anders W. Sandvik, Boston University Lecture outline Introduction
More informationarxiv: v3 [cond-mat.dis-nn] 18 May 2011
Exact mean field inference in asymmetric kinetic Ising systems M. Mézard and J. Sakellariou Laboratoire de Physique Théorique et Modèles Statistiques, CNRS and Université Paris-Sud, Bât 100, 91405 Orsay
More informationStochastic series expansion (SSE) and ground-state projection
Institute of Physics, Chinese Academy of Sciences, Beijing, October 31, 2014 Stochastic series expansion (SSE) and ground-state projection Anders W Sandvik, Boston University Review article on quantum
More informationGround State Projector QMC in the valence-bond basis
Quantum Monte Carlo Methods at Work for Novel Phases of Matter Trieste, Italy, Jan 23 - Feb 3, 2012 Ground State Projector QMC in the valence-bond basis Anders. Sandvik, Boston University Outline: The
More informationRenormalization Group for the Two-Dimensional Ising Model
Chapter 8 Renormalization Group for the Two-Dimensional Ising Model The two-dimensional (2D) Ising model is arguably the most important in statistical physics. This special status is due to Lars Onsager
More informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
More informationQuantum simulation with string-bond states: Joining PEPS and Monte Carlo
Quantum simulation with string-bond states: Joining PEPS and Monte Carlo N. Schuch 1, A. Sfondrini 1,2, F. Mezzacapo 1, J. Cerrillo 1,3, M. Wolf 1,4, F. Verstraete 5, I. Cirac 1 1 Max-Planck-Institute
More informationHanoi 7/11/2018. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
Hanoi 7/11/2018 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Finite size effects and Reweighting methods 1. Finite size effects 2. Single histogram method 3. Multiple histogram method 4. Wang-Landau
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 information(G). This remains probably true in most random magnets. It raises the question of. On correlation functions in random magnets LETTER TO THE EDITOR
J. Phys. C: Solid State Phys., 14 (1981) L539-L544. Printed in Great Britain LETTER TO THE EDITOR On correlation functions in random magnets Bernard Derridai and Henk Hilhorst$ t Service de Physique ThCorique,
More informationScaling Theory. Roger Herrigel Advisor: Helmut Katzgraber
Scaling Theory Roger Herrigel Advisor: Helmut Katzgraber 7.4.2007 Outline The scaling hypothesis Critical exponents The scaling hypothesis Derivation of the scaling relations Heuristic explanation Kadanoff
More informationThe non-backtracking operator
The non-backtracking operator Florent Krzakala LPS, Ecole Normale Supérieure in collaboration with Paris: L. Zdeborova, A. Saade Rome: A. Decelle Würzburg: J. Reichardt Santa Fe: C. Moore, P. Zhang Berkeley:
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 informationRenormalization group maps for Ising models in lattice gas variables
Renormalization group maps for Ising models in lattice gas variables Department of Mathematics, University of Arizona Supported by NSF grant DMS-0758649 http://www.math.arizona.edu/ e tgk RG in lattice
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationComplex Systems Methods 9. Critical Phenomena: The Renormalization Group
Complex Systems Methods 9. Critical Phenomena: The Renormalization Group Eckehard Olbrich e.olbrich@gmx.de http://personal-homepages.mis.mpg.de/olbrich/complex systems.html Potsdam WS 2007/08 Olbrich (Leipzig)
More informationMonte Carlo study of the Baxter-Wu model
Monte Carlo study of the Baxter-Wu model Nir Schreiber and Dr. Joan Adler Monte Carlo study of the Baxter-Wu model p.1/40 Outline Theory of phase transitions, Monte Carlo simulations and finite size scaling
More information(1) Consider the ferromagnetic XY model, with
PHYSICS 10A : STATISTICAL PHYSICS HW ASSIGNMENT #7 (1 Consider the ferromagnetic XY model, with Ĥ = i
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
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 informationLikelihood Inference for Lattice Spatial Processes
Likelihood Inference for Lattice Spatial Processes Donghoh Kim November 30, 2004 Donghoh Kim 1/24 Go to 1234567891011121314151617 FULL Lattice Processes Model : The Ising Model (1925), The Potts Model
More informationThe Random Matching Problem
The Random Matching Problem Enrico Maria Malatesta Universitá di Milano October 21st, 2016 Enrico Maria Malatesta (UniMi) The Random Matching Problem October 21st, 2016 1 / 15 Outline 1 Disordered Systems
More informationPhase Transitions in the Coloring of Random Graphs
Phase Transitions in the Coloring of Random Graphs Lenka Zdeborová (LPTMS, Orsay) In collaboration with: F. Krząkała (ESPCI Paris) G. Semerjian (ENS Paris) A. Montanari (Standford) F. Ricci-Tersenghi (La
More informationS i J <ij> h mf = h + Jzm (4) and m, the magnetisation per spin, is just the mean value of any given spin. S i = S k k (5) N.
Statistical Physics Section 10: Mean-Field heory of the Ising Model Unfortunately one cannot solve exactly the Ising model or many other interesting models) on a three dimensional lattice. herefore one
More informationExact Solution of Ising Model on a Small-World Network
1/2 Exact Solution of Ising Model on a Small-World Network J. Viana Lopes, Yu. G. Pogorelov, J.M.B. Lopes dos Santos CFP, Departamento de Física, Faculdade de Ciências, Universidade do Porto and Raul Toral
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis
NUMERICAL APPROACHES TO QUANTUM MANY-BODY SYSTEMS, IPAM, January 29, 2009 Quantum Monte Carlo Simulations in the Valence Bond Basis Anders W. Sandvik, Boston University Collaborators Kevin Beach (U. of
More informationZacros. Software Package Development: Pushing the Frontiers of Kinetic Monte Carlo Simulation in Catalysis
Zacros Software Package Development: Pushing the Frontiers of Kinetic Monte Carlo Simulation in Catalysis Jens H Nielsen, Mayeul D'Avezac, James Hetherington & Michail Stamatakis Introduction to Zacros
More informationThe Phase Transition of the 2D-Ising Model
The Phase Transition of the 2D-Ising Model Lilian Witthauer and Manuel Dieterle Summer Term 2007 Contents 1 2D-Ising Model 2 1.1 Calculation of the Physical Quantities............... 2 2 Location of the
More informationReconstruction for Models on Random Graphs
1 Stanford University, 2 ENS Paris October 21, 2007 Outline 1 The Reconstruction Problem 2 Related work and applications 3 Results The Reconstruction Problem: A Story Alice and Bob Alice, Bob and G root
More informationGraphical Representations and Cluster Algorithms
Graphical Representations and Cluster Algorithms Jon Machta University of Massachusetts Amherst Newton Institute, March 27, 2008 Outline Introduction to graphical representations and cluster algorithms
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 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 informationMonte Carlo Simulation of the 2D Ising model
Monte Carlo Simulation of the 2D Ising model Emanuel Schmidt, F44 April 6, 2 Introduction Monte Carlo methods are a powerful tool to solve problems numerically which are dicult to be handled analytically.
More informationSpin glasses and Stein s method
UC Berkeley Spin glasses Magnetic materials with strange behavior. ot explained by ferromagnetic models like the Ising model. Theoretically studied since the 70 s. An important example: Sherrington-Kirkpatrick
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 informationIntro. Each particle has energy that we assume to be an integer E i. Any single-particle energy is equally probable for exchange, except zero, assume
Intro Take N particles 5 5 5 5 5 5 Each particle has energy that we assume to be an integer E i (above, all have 5) Particle pairs can exchange energy E i! E i +1andE j! E j 1 5 4 5 6 5 5 Total number
More informationA variational approach to Ising spin glasses in finite dimensions
. Phys. A: Math. Gen. 31 1998) 4127 4140. Printed in the UK PII: S0305-447098)89176-2 A variational approach to Ising spin glasses in finite dimensions R Baviera, M Pasquini and M Serva Dipartimento di
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 informationPhase transitions and finite-size scaling
Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization
More informationKrammers-Wannier Duality in Lattice Systems
Krammers-Wannier Duality in Lattice Systems Sreekar Voleti 1 1 Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7. (Dated: December 9, 2018) I. INTRODUCTION It was shown by
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 informationMonte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya Interaction
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 663 667 c International Academic Publishers Vol. 46, No. 4, October 15, 2006 Monte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya
More informationPhase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension: Specifically For Square Lattices
IOSR Journal of Applied Physics (IOSR-JAP) ISS: 2278-4861. Volume 2, Issue 3 (ov. - Dec. 2012), PP 31-37 Phase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension:
More informationA Monte Carlo Implementation of the Ising Model in Python
A Monte Carlo Implementation of the Ising Model in Python Alexey Khorev alexey.s.khorev@gmail.com 2017.08.29 Contents 1 Theory 1 1.1 Introduction...................................... 1 1.2 Model.........................................
More information8 Error analysis: jackknife & bootstrap
8 Error analysis: jackknife & bootstrap As discussed before, it is no problem to calculate the expectation values and statistical error estimates of normal observables from Monte Carlo. However, often
More informationReplica cluster variational method: the replica symmetric solution for the 2D random bond Ising model
IOP PUBLISHING OURNAL OF PHYSICS A: MATHEMATICAL AND THEORETICAL. Phys. A: Math. Theor. 46 (2013) 135001 (20pp) doi:10.1088/1751-8113/46/13/135001 Replica cluster variational method: the replica symmetric
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 informationPhase Transitions of Random Binary Magnetic Square Lattice Ising Systems
I. Q. Sikakana Department of Physics and Non-Destructive Testing, Vaal University of Technology, Vanderbijlpark, 1900, South Africa e-mail: ike@vut.ac.za Abstract Binary magnetic square lattice Ising system
More informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ Bayesian paradigm Consistent use of probability theory
More informationSimulation of the two-dimensional square-lattice Lenz-Ising model in Python
Senior Thesis Simulation of the two-dimensional square-lattice Lenz-Ising model in Python Author: Richard Munroe Advisor: Dr. Edward Brash 2012-10-17 Abstract A program to simulate the the two-dimensional
More informationApproximate Inference Part 1 of 2
Approximate Inference Part 1 of 2 Tom Minka Microsoft Research, Cambridge, UK Machine Learning Summer School 2009 http://mlg.eng.cam.ac.uk/mlss09/ 1 Bayesian paradigm Consistent use of probability theory
More informationConcentration of magnetic transitions in dilute magnetic materials
Journal of Physics: Conference Series OPEN ACCESS Concentration of magnetic transitions in dilute magnetic materials To cite this article: V I Beloon et al 04 J. Phys.: Conf. Ser. 490 065 Recent citations
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 informationS j H o = gµ o H o. j=1
LECTURE 17 Ferromagnetism (Refs.: Sections 10.6-10.7 of Reif; Book by J. S. Smart, Effective Field Theories of Magnetism) Consider a solid consisting of N identical atoms arranged in a regular lattice.
More informationMonte Carlo and cold gases. Lode Pollet.
Monte Carlo and cold gases Lode Pollet lpollet@physics.harvard.edu 1 Outline Classical Monte Carlo The Monte Carlo trick Markov chains Metropolis algorithm Ising model critical slowing down Quantum Monte
More informationXVI International Congress on Mathematical Physics
Aug 2009 XVI International Congress on Mathematical Physics Underlying geometry: finite graph G=(V,E ). Set of possible configurations: V (assignments of +/- spins to the vertices). Probability of a configuration
More informationInvaded cluster dynamics for frustrated models
PHYSICAL REVIEW E VOLUME 57, NUMBER 1 JANUARY 1998 Invaded cluster dynamics for frustrated models Giancarlo Franzese, 1, * Vittorio Cataudella, 1, * and Antonio Coniglio 1,2, * 1 INFM, Unità di Napoli,
More informationSpin liquid phases in strongly correlated lattice models
Spin liquid phases in strongly correlated lattice models Sandro Sorella Wenjun Hu, F. Becca SISSA, IOM DEMOCRITOS, Trieste Seiji Yunoki, Y. Otsuka Riken, Kobe, Japan (K-computer) Williamsburg, 14 June
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis. Anders Sandvik, Boston University
Quantum Monte Carlo Simulations in the Valence Bond Basis Anders Sandvik, Boston University Outline The valence bond basis for S=1/2 spins Projector QMC in the valence bond basis Heisenberg model with
More informationPhysics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion
Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to
More informationIsing Model. Ising Lattice. E. J. Maginn, J. K. Shah
Ising Lattice E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop, Brazil Ising Lattice Model Consider a
More informationInverse spin glass and related maximum entropy problems
Inverse spin glass and related maximum entropy problems Michele Castellana and William Bialek,2 Joseph Henry Laboratories of Physics and Lewis Sigler Institute for Integrative Genomics, Princeton University,
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationError Bounds for Sliding Window Classifiers
Error Bounds for Sliding Window Classifiers Yaroslav Bulatov January 8, 00 Abstract Below I provide examples of some results I developed that let me bound error rate of Bayes optimal sliding window classifiers
More informationVSOP19, Quy Nhon 3-18/08/2013. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
VSOP19, Quy Nhon 3-18/08/2013 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Part III. Finite size effects and Reweighting methods III.1. Finite size effects III.2. Single histogram method III.3.
More informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez May 2009, Lecture #2 Statistical Mechanics: Structure S = k B log W! Boltzmann's headstone at the Vienna Zentralfriedhof.!! Photo
More informationThe Missing-Index Problem in Process Control
The Missing-Index Problem in Process Control Joseph G. Voelkel CQAS, KGCOE, RIT October 2009 Voelkel (RIT) Missing-Index Problem 10/09 1 / 35 Topics 1 Examples Introduction to the Problem 2 Models, Inference,
More informationMagnetism at finite temperature: molecular field, phase transitions
Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field
More informationLecture V: Multicanonical Simulations.
Lecture V: Multicanonical Simulations. 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs (2d Ising and Potts models) 4. Re-Weighting to the Canonical Ensemble 5. Energy and Specific
More informationOn the number of circuits in random graphs. Guilhem Semerjian. [ joint work with Enzo Marinari and Rémi Monasson ]
On the number of circuits in random graphs Guilhem Semerjian [ joint work with Enzo Marinari and Rémi Monasson ] [ Europhys. Lett. 73, 8 (2006) ] and [ cond-mat/0603657 ] Orsay 13-04-2006 Outline of the
More informationOptimized statistical ensembles for slowly equilibrating classical and quantum systems
Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,
More informationMatrix-Product states: Properties and Extensions
New Development of Numerical Simulations in Low-Dimensional Quantum Systems: From Density Matrix Renormalization Group to Tensor Network Formulations October 27-29, 2010, Yukawa Institute for Theoretical
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 informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationAdvanced Sampling Algorithms
+ Advanced Sampling Algorithms + Mobashir Mohammad Hirak Sarkar Parvathy Sudhir Yamilet Serrano Llerena Advanced Sampling Algorithms Aditya Kulkarni Tobias Bertelsen Nirandika Wanigasekara Malay Singh
More informationOn Markov chain Monte Carlo methods for tall data
On Markov chain Monte Carlo methods for tall data Remi Bardenet, Arnaud Doucet, Chris Holmes Paper review by: David Carlson October 29, 2016 Introduction Many data sets in machine learning and computational
More informationReal-Space Renormalisation
Real-Space Renormalisation LMU München June 23, 2009 Real-Space Renormalisation of the One-Dimensional Ising Model Ernst Ising Real-Space Renormalisation of the One-Dimensional Ising Model Ernst Ising
More informationKinetic Monte Carlo (KMC) Kinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationMonte Carlo Simulation of the Ising Model. Abstract
Monte Carlo Simulation of the Ising Model Saryu Jindal 1 1 Department of Chemical Engineering and Material Sciences, University of California, Davis, CA 95616 (Dated: June 9, 2007) Abstract This paper
More informationImproving the Convergence of Back-Propogation Learning with Second Order Methods
the of Back-Propogation Learning with Second Order Methods Sue Becker and Yann le Cun, Sept 1988 Kasey Bray, October 2017 Table of Contents 1 with Back-Propagation 2 the of BP 3 A Computationally Feasible
More informationAgeing properties of three-dimensional pure and site-diluted Ising ferromagnets
Journal of Physics: Conference Series OPEN ACCESS Ageing properties of three-dimensional pure and site-diluted Ising ferromagnets To cite this article: Vladimir V Prudnikov et al 2014 J. Phys.: Conf. Ser.
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 information4. Cluster update algorithms
4. Cluster update algorithms Cluster update algorithms are the most succesful global update methods in use. These methods update the variables globally, in one step, whereas the standard local methods
More informationCluster Functional Renormalization Group
Cluster Functional Renormalization Group Johannes Reuther Free University Berlin Helmholtz-Center Berlin California Institute of Technology (Pasadena) Lefkada, September 26, 2014 Johannes Reuther Cluster
More information