A brief introduction to the inverse Ising problem and some algorithms to solve it
|
|
- Beverly Williamson
- 5 years ago
- Views:
Transcription
1 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 and Aurelien Decelle
2 Simplifications for this talk Ising variables s i = ±1 At most pairwise (two body) interactions. The most general Hamiltonian is H(s J, h) = X J ij s i s j + X i h i s i with corresponding measure P (s 1,...,s N )= 1 Z(J, h) exp i=j J ij s i s j + i h i s i
3 The direct problem (main pb. in stat. mech.) Given the Hamiltonian, compute the free energy F (J, h) = log Z(J, h) = log X exp X J ij s i s j + X {s i } i h i s i! and average values hoi = X s O(s)P (s) The sum is over exponentially many terms
4 The inverse Ising problem Given data generated according to which may be P (s 1,...,s N )= 1 Z(J, h) exp i=j J ij s i s j + i h i s i either M configurations of spins either magnetizations and correlations m i = hs i i C ij = hs i s j i m i m j GOAL: estimate couplings and fields (J, h)
5 The exact solution Maximize the log-likelihood = X i L(J, h s) = 1 M log M Y h i hs i i + X ij k=1 P (s (k) J, h) = J ij hs i s j i log Z(J, h) = X i h i m i + X ij J ij (C ij + m i m j )+F (J, h) Taking the derivatives m i hi F (J, h) =0 =) m i (DATA) = m i (J, h) C ij + m i m j Jij F (J, h) =0 =) C ij (DATA) = C ij (J, h)
6 Input data: Magnetizations and Correlations Less information than having configurations If only m i = hs i i C ij = hs i s j i m i m j are given then maximum entropy principle implies Hamiltonian contains only single-body and two-body interactions H(s J, h) = X J ij s i s j + X i h i s i
7 Brute force solution by Monte Carlo 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
8 Mean field approximations (MFA) Log-likelihood L(J, h s) = 1 M log M Y P (s (k) J, h) = = X i h i m i + X ij k=1 J ij (C ij + m i m j ) log Z(J, h) F MFA (J, h) m MFA i hi F MFA (J, h) m MFA i (J, h) =m i (DATA) Cij MFA (J, h) =C ij (DATA)
9 MFA to the free-energy naive mean-field (nmf) F nmf = i 1+mi 1 H + H 2 P (s) = Q i P i(s i ) mi 2 + 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
10 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
11 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 = i + i H 1+mi + H 2 h i m i + i=j 1 mi 2 + J ij m i m j J 2 ij(1 m 2 i )(1 m 2 j)
12 MFA to the free-energy Bethe approximation (BA) tries to include any correlations between n.n. spins in principle no small J required (but beware to phase transitions) factorization over links exact on trees P (s) = Q i P i(s i ) Q ij P ij (s i,s j ) P i (s i )P j (s j )
13 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 t ij = tanh(j ij ) F BA / c ij = 0 + c ij (m i,m j,t ij ) = 1 1+t 2 ij (1 t 2 ij 2t )2 4t ij (m i t ij m j )(m j t ij m i ) ij m i m j
14 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 m (j) i i : magnetization of i in absence of j t ij 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)
15 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 ij m j Jij(1 2 m 2 j)m i +... j
16 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) ij j ( 1 ) ij j ( (CnMF 1 ) ij = ( 1 1 TAP ) ij (C 1 TAP ) ij = ij 1 m 2 i 1 1 m 2 i J ij, + k J 2 ik (1 m2 k ) ij J ij +2J 2 ijm i m j
17 Computing correlations by linear response in BA Analytic expression for the linear responses in BA (C 1 BA ) ij = ( BA ) ij 1 1 m 2 i Coincide with the fixed point of Susceptibility Propagation No need to run any algorithm! k t ik f 2 (m k,m i,t ik ) 1 t 2 ik f(m k,m i,t ik ) 2 ij t ij f 1 (m j,m i,t ij ) 1 t 2 ij f(m j,m i,t ij ) 2
18 Solving the inverse problem by MFA Match measured magnetizations and correlations with MF approximated magnetizations and linear responses m MFA i (J, h) =m i (DATA) MFA ij (J, h) =C ij (DATA) Under the Bethe approx. one could use either c BA ij or BA ij for n.n. correlations. Which one is better?
19 Zero field case is simpler If all field are zero, then magnetizations are null by symmetry, and expressions simplify to naive MF (C 1 nmf ) ij = ij J ij, " # TAP (C 1 TAP ) ij = Bethe (C 1 BA ) ij = " 1+ X k 1+ X k J 2 ik t 2 ik 1 t 2 ik ij J ij, # ij t ij, 1 t 2 ij
20 Exactly solvable case for the inverse Ising problem? Curie-Weiss model, fully connected nmf approximation J ij = /(N 1) ij with i 6= j N=10,20,40 exact correlations ii
21 Exactly solvable case for the inverse Ising problem? Curie-Weiss model, fully connected more MFA (N=20) ij with i 6= j ii J ij = /(N 1) TAP ~ BA nmf 3rd th
22 Exactly solvable case for the inverse Ising problem? Bethe approximation on trees is ok Bethe approximation on random graph is ok only far from the critical point (as nmf for the Curie-Weiss model) How much the paramagnetic properties of a model on a finite size random graph are different from those of the same model defined on a tree? see recent works on finite size corrections for models defined on random graphs (Lucibello and Morone)
23 Numerical results on estimating correlations 1e+02 C 1 N 2 i,j (C ij C ij ) 2 1 1e+00 1 C 1e-02 1e-04 1e-06 2D ferromagnet diluted (p=0.7) nmf TAP 3 rd 4 th BA 1e
24 Matching data and MF predictions 0 1 C ij... 1? C A = 0 11 ij... 1 C A C ij 1 ij NN Usually only off-diagonal elements are used J ij = (C 1 ) ij and diagonal elements are ignored...
25 MFA for the inverse Ising problem J BA ij = atanh (C 1 nmf ) ij = (C 1 TAP ) ij = J TAP ij = + ij =) J nmf 1 m 2 J ij, i 1 1 m 2 + Jik 2 (1 m2 k ) i k 1 8mi m j (C 1 ) ij 1 4m i m j ij = (C 1 ) ij 1 2(C (1 m 2 i ) )(1 m2 j )(C 1 ) 2 ij m i m j ij ij J ij +2J 2 ijm i m j 1 2(C 1 ) ij 1 + 4(1 m 2 i )(1 m2 j )(C 1 ) 2 ij 2m i m j (C 1 ) ij 2 4(C 1 ) 2 ij equal for m=0
26 More MFA 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 Jij SM = (C 1 ) ij + Jij IP C ij (1 m 2 i )(1 m2 j ) (C ij) 2
27 Numerical results for the inverse Ising problem 10 2D ferromagnet N=7 2 diluted (p=0.7) J IP TAP SM BA BA norm
28 Improving correlations by a normalization trick In ferromagnetic models with loops, linear response correlations in BA are too strong because of loops, which are unexpected in SuscProp 2D lattice Leading to ii > 1 which is unphysical Trick: enforce ii =1 by a normalization b ij = p ij ii jj
29 Make MFA & LR consistent Consistency is more important than truth (S. Ting) Add Lagrange multipliers to your referred MF free-energy F MFA ({m i }, {C ij },...) to enforce consistency with linear response estimates ii =1 m 2 i ij = C ij free energy minimum curvature free energy minimum location
30 General framework (MFA + LR) F = F MFA ({m i }, {C ij },...)+ X i im 2 i + X i<j ijc ij Your preferred MFA can be set to zero to recover known approx. or used to satisfy ii =1 m 2 i ij = C ij
31 Nearest-neighbor correlation (2D square lattice) 0.8 NN correlations Bethe approx. LR ME exact ME+LR
32 Off-diagonal constraint only 0.01 NMF/TAP (λ=0) Bethe (λ=0) Bethe (λ=0) [KR] Bethe Plaquette J β Bethe + off-diagonal constraints = SM
33 Random field Ising model 2D square lattice D RFIM <h> = 0.0 σ h = 0.2 β = Bethe, no constraint Bethe, diag. constraint Bethe, both constraints m infer - m exact m exact
34 Random field Ising model 2D square lattice 1e-02 2D RFIM <h> = 0.0 σ h = 0.2 β = e-03 m infer - m exact 1e-04 1e-05 Bethe, no constraint Bethe, diag. constraint Bethe, both constraints 1e m exact
35 Random field Ising model 2D square lattice 1e-01 2D RFIM <h> = 0.0 σ h = 0.2 β = e-02 self correlations 1e-03 C infer - C exact 1e-04 1e-05 1e-06 Bethe, no constraint (LR) Bethe, no constraint (ME) Bethe, diag. constraint (LR) Bethe, both constraints (ME=LR) 1e C exact
36 Random field Ising model 2D square lattice 1e-01 2D RFIM <h> = 0.0 σ h = 0.2 β = 0.25 NNN correlations NN correlations 1e-02 C infer - C exact 1e-03 1e-04 1e-05 Bethe, no constraint (LR) Bethe, diag. constraint (LR) Bethe, both constraints (ME=LR) C exact
37 Input data: Configurations More information than knowing only m i = hs i i C ij = hs i s j i m i m j In principle one can access to all higher order correlations (but these are much more noisy) Many different inference algorithms. Among these: Adaptive cluster expansion cluster configurations & apply MFA within any state
38 Pseudo-likelihood method (PLM) For each variable define a conditional probability P i (s i s \i )= exp[s i(h i + P j J ijs j )] 2 cosh(h i + P j J ijs j ) Maximize the local log-likelihood L i = hlog P i (s i s \i )i = J ij (C ij + m i m j ) hlog 2 cosh(h i + X j J ij s j )i h i m i + X j to estimate h i and J ij Note that for each coupling J ij PLM returns 2 estimates Better maximizing PL(h, J) = X i L i
39 PLM vs. MFA 10 2D ferromagnet N=7 2 =49 diluted (p=0.7) 1 M=5000 samples J IP TAP SM BA PLM BA norm
40 Inferring topology in sparse models For simplicity let s assume J ij 2 {0, } non-zero couplings are sparse Maximize L1-regularized pseudo-likelihood PL (h, J) = X i L1-regularization gives a bias to the estimates! L i X J ij ij
41 Couplings inferred by PLM β=0.5 β=0.9 β=0.9 L1 regularized 1e-4 1e-2 1 1e-4 1e-2 1 log(j inferred) 1e-12 1e-8 1e-4 1 2D Ising model (30% dilution) M=4500
42 Decimation procedure No L1-regularizer -> no bias Maximize PL(h, J) Set to zero a constant fraction of couplings (those inferred to be the smallest) Maximize again PL(h, J) to zero (this is impossible within a MFA) only on couplings still not set Iterate until maximum of sensibly PL(h, J) starts decreasing
43 PLM + decimation tplf Error tplf Error Number of non-decimated couplings
44 PLM + decimation True Positive Rate PLM+L1 δ=1e-8 PLM+L1 δ=1e-6 PLM+L1 δ=1e-4 PLM+L1 δ=1e-2 PLM+L1 δ=1e-1 PLM+Decimation STOP point True Negative Rate
45 Some conclusions Mean field approximations Inverse problem harder than direct problem Requires (at least) improvement in the direct problem Fundamental problem of going beyond Bethe and trees... Pseudo-likelihood method Better performances in general Specially well suited for inferring topology in sparse models via L1-regularization, thresholding or decimation
Some recent results on the inverse Ising problem. Federico Ricci-Tersenghi Physics Department Sapienza University, Roma
Some recent results on the inverse Ising problem Federico Ricci-Tersenghi Physics Department Sapienza University, Roma Outline of the talk Bethe approx. for inverse Ising problem Comparison among several
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 information2 The Curie Weiss Model
2 The Curie Weiss Model In statistical mechanics, a mean-field approximation is often used to approximate a model by a simpler one, whose global behavior can be studied with the help of explicit computations.
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 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 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 informationMagnetic ordering of local moments
Magnetic ordering Types of magnetic structure Ground state of the Heisenberg ferromagnet and antiferromagnet Spin wave High temperature susceptibility Mean field theory Magnetic ordering of local moments
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 informationPhase transition and spontaneous symmetry breaking
Phys60.nb 111 8 Phase transition and spontaneous symmetry breaking 8.1. Questions: Q1: Symmetry: if a the Hamiltonian of a system has certain symmetry, can the system have a lower symmetry? Q: Analyticity:
More informationBayesian Learning in Undirected Graphical Models
Bayesian Learning in Undirected Graphical Models Zoubin Ghahramani Gatsby Computational Neuroscience Unit University College London, UK http://www.gatsby.ucl.ac.uk/ Work with: Iain Murray and Hyun-Chul
More informationStructured Variational Inference
Structured Variational Inference Sargur srihari@cedar.buffalo.edu 1 Topics 1. Structured Variational Approximations 1. The Mean Field Approximation 1. The Mean Field Energy 2. Maximizing the energy functional:
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 information(1) Consider the ferromagnetic XY model, with
PHYSICS 10A : STATISTICAL PHYSICS HW ASSIGNMENT #7 (1 Consider the ferromagnetic XY model, with Ĥ = i
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 informationLinear Response for Approximate Inference
Linear Response for Approximate Inference Max Welling Department of Computer Science University of Toronto Toronto M5S 3G4 Canada welling@cs.utoronto.ca Yee Whye Teh Computer Science Division University
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 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 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 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 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 informationDoes Better Inference mean Better Learning?
Does Better Inference mean Better Learning? Andrew E. Gelfand, Rina Dechter & Alexander Ihler Department of Computer Science University of California, Irvine {agelfand,dechter,ihler}@ics.uci.edu Abstract
More information12 : Variational Inference I
10-708: Probabilistic Graphical Models, Spring 2015 12 : Variational Inference I Lecturer: Eric P. Xing Scribes: Fattaneh Jabbari, Eric Lei, Evan Shapiro 1 Introduction Probabilistic inference is one of
More informationFractional Belief Propagation
Fractional Belief Propagation im iegerinck and Tom Heskes S, niversity of ijmegen Geert Grooteplein 21, 6525 EZ, ijmegen, the etherlands wimw,tom @snn.kun.nl Abstract e consider loopy belief propagation
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 informationPhase Transitions in Physics and Computer Science. Cristopher Moore University of New Mexico and the Santa Fe Institute
Phase Transitions in Physics and Computer Science Cristopher Moore University of New Mexico and the Santa Fe Institute Magnetism When cold enough, Iron will stay magnetized, and even magnetize spontaneously
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 information13 : Variational Inference: Loopy Belief Propagation and Mean Field
10-708: Probabilistic Graphical Models 10-708, Spring 2012 13 : Variational Inference: Loopy Belief Propagation and Mean Field Lecturer: Eric P. Xing Scribes: Peter Schulam and William Wang 1 Introduction
More informationThe mixed-spins 1/2 and 3/2 Blume Capel model with a random crystal field
The mixed-spins 1/2 and 3/2 Blume Capel model with a random crystal field Erhan Albayrak Erciyes University, Department of Physics, 38039, Kayseri, Turkey (Received 25 August 2011; revised manuscript received
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 informationThe 1+1-dimensional Ising model
Chapter 4 The 1+1-dimensional Ising model The 1+1-dimensional Ising model is one of the most important models in statistical mechanics. It is an interacting system, and behaves accordingly. Yet for a variety
More information13 : Variational Inference: Loopy Belief Propagation
10-708: Probabilistic Graphical Models 10-708, Spring 2014 13 : Variational Inference: Loopy Belief Propagation Lecturer: Eric P. Xing Scribes: Rajarshi Das, Zhengzhong Liu, Dishan Gupta 1 Introduction
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
1 Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 2, 24 March 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
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 informationPHYSICS 219 Homework 2 Due in class, Wednesday May 3. Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235.
PHYSICS 219 Homework 2 Due in class, Wednesday May 3 Note: Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235. No lecture: May 8 (I m away at a meeting) and May 29 (holiday).
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 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 informationQuantum Integrability and Algebraic Geometry
Quantum Integrability and Algebraic Geometry Marcio J. Martins Universidade Federal de São Carlos Departamento de Física July 218 References R.J. Baxter, Exactly Solved Models in Statistical Mechanics,
More informationThe XY-Model. David-Alexander Robinson Sch th January 2012
The XY-Model David-Alexander Robinson Sch. 08332461 17th January 2012 Contents 1 Introduction & Theory 2 1.1 The XY-Model............................... 2 1.2 Markov Chains...............................
More informationLearning MN Parameters with Approximation. Sargur Srihari
Learning MN Parameters with Approximation Sargur srihari@cedar.buffalo.edu 1 Topics Iterative exact learning of MN parameters Difficulty with exact methods Approximate methods Approximate Inference Belief
More informationTHE PHYSICS OF COUNTING AND SAMPLING ON RANDOM INSTANCES. Lenka Zdeborová
THE PHYSICS OF COUNTING AND SAMPLING ON RANDOM INSTANCES Lenka Zdeborová (CEA Saclay and CNRS, France) MAIN CONTRIBUTORS TO THE PHYSICS UNDERSTANDING OF RANDOM INSTANCES Braunstein, Franz, Kabashima, Kirkpatrick,
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 informationSpin models on random graphs with controlled topologies beyond degree constraints
Spin models on random graphs with controlled topologies beyond degree constraints ACC Coolen (London) and CJ Pérez-Vicente (Barcelona) (work in progress) Motivation - processes on random graphs Deformed
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More informationVariational algorithms for marginal MAP
Variational algorithms for marginal MAP Alexander Ihler UC Irvine CIOG Workshop November 2011 Variational algorithms for marginal MAP Alexander Ihler UC Irvine CIOG Workshop November 2011 Work with Qiang
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 informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
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 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 informationQuantum annealing by ferromagnetic interaction with the mean-field scheme
Quantum annealing by ferromagnetic interaction with the mean-field scheme Sei Suzuki and Hidetoshi Nishimori Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan
More informationarxiv: v1 [cond-mat.dis-nn] 16 Feb 2018
Parallel Tempering for the planted clique problem Angelini Maria Chiara 1 1 Dipartimento di Fisica, Ed. Marconi, Sapienza Università di Roma, P.le A. Moro 2, 00185 Roma Italy arxiv:1802.05903v1 [cond-mat.dis-nn]
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 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 informationPhase transitions in discrete structures
Phase transitions in discrete structures Amin Coja-Oghlan Goethe University Frankfurt Overview 1 The physics approach. [following Mézard, Montanari 09] Basics. Replica symmetry ( Belief Propagation ).
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 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 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 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 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 information17 Variational Inference
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.438 Algorithms for Inference Fall 2014 17 Variational Inference Prompted by loopy graphs for which exact
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 informationProbabilistic Graphical Models
School of Computer Science Probabilistic Graphical Models Variational Inference IV: Variational Principle II Junming Yin Lecture 17, March 21, 2012 X 1 X 1 X 1 X 1 X 2 X 3 X 2 X 2 X 3 X 3 Reading: X 4
More information9 Classification. 9.1 Linear Classifiers
9 Classification This topic returns to prediction. Unlike linear regression where we were predicting a numeric value, in this case we are predicting a class: winner or loser, yes or no, rich or poor, positive
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationA Monte Carlo Study of the Specific Heat of Diluted Antiferromagnets
A Monte Carlo Study of the Specific Heat of Diluted Antiferromagnets M. Staats and K. D. Usadel email: michael@thp.uni-duisburg.de Theoretische Physik and SFB 166 Gerhard-Mercator-Universität Gesamthochschule
More informationProbabilistic Graphical Models Lecture Notes Fall 2009
Probabilistic Graphical Models Lecture Notes Fall 2009 October 28, 2009 Byoung-Tak Zhang School of omputer Science and Engineering & ognitive Science, Brain Science, and Bioinformatics Seoul National University
More informationProbabilistic and Bayesian Machine Learning
Probabilistic and Bayesian Machine Learning Day 4: Expectation and Belief Propagation Yee Whye Teh ywteh@gatsby.ucl.ac.uk Gatsby Computational Neuroscience Unit University College London http://www.gatsby.ucl.ac.uk/
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 informationBayesian Learning in Undirected Graphical Models
Bayesian Learning in Undirected Graphical Models Zoubin Ghahramani Gatsby Computational Neuroscience Unit University College London, UK http://www.gatsby.ucl.ac.uk/ and Center for Automated Learning and
More informationStatistical description of magnetic domains in the Ising model
arxiv:0804.3522v1 [cond-mat.stat-mech] 22 Apr 2008 Statistical description of magnetic domains in the Ising model K. Lukierska-Walasek Institute of Physics University of Zielona Góra ul. Z. Szafrana 4a,
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 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 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 informationMagnets, 1D quantum system, and quantum Phase transitions
134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider
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 informationThe Generalized Distributive Law and Free Energy Minimization
The Generalized Distributive Law and Free Energy Minimization Srinivas M. Aji Robert J. McEliece Rainfinity, Inc. Department of Electrical Engineering 87 N. Raymond Ave. Suite 200 California Institute
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 informationRenormalization Group analysis of 2D Ising model
Renormalization Group analysis of D Ising model Amir Bar January 7, 013 1 Introduction In this tutorial we will see explicitly how RG can be used to probe the phase diagram of d > 1 systems, focusing as
More informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More information