RC Expansion/Edwards Sokal Measure

Size: px
Start display at page:

Download "RC Expansion/Edwards Sokal Measure"

Transcription

1 RC Expansion/Edwards Sokal Measure April 11, 2007 Let Λ d be finite and consider the q state Potts model with couplings {J x, } and free boundar condition (no boundar condition). The Hamiltonian (q = 2 is Ising) is H() = J x, (δ x, 1), where x {1,..., q}. The partition function is [] = e βh() = e β[jx,(δx, 1)] (1) Now consider the random cluster model with parameter q and bond weights {p x, } with partition function given b [] = q c() (p x, ) (1 p x, ), (2) : x,=1 : x,=0 where c() denotes the number of connected components of the bond configuration. Note that if we set e βjx,(δx, 1) = δ x, + (1 δ x, )e βjx, = δ x, p x, + (1 p x, ), p x, = 1 e βjx,. Plugging this into (1), we see that e βh() = (δ x, p x, + (1 p x, )). (3) 1

2 If we expand out the product, we get (δ x, p x, ) Ω Ω(Λ) Ω / Ω (1 p x, ), where Ω(Λ) denotes the set of bonds connecting points in Λ. Each term above can be identified as a bond configuration where each bond is occupied if and onl if it lies in Ω. The above sum can now be rewritten as (δ x, p x, ) (1 p x, ). : x,=1 The partition function (1) now becomes (δ x, p x, ) : x,=1 : x,=0 : x,=0 (1 p x, ). (4) Observe that the above sum is finite and so far is simpl summed over all possible bond and spin configurations, and hence we ma interchange the order of summation to obtain (δ x, p x, ) (1 p x, ). (5) : x,=1 : x,=0 From (4) we see that if is such that a bond is present when x, then the contribution to the sum is zero. Similarl, in (5), we see that if is such that x when there is a bond present between x and, then the contribution to the sum is zero. So to get rid of the δ x,, we introduce the constraint function (, ) = () = (), which is the indicator function of the fact that and are compatible in the sense that x = whenever x, = 1 and x, = 0 whenever x. We can now rewrite the Potts partition function as [] = () (p x, ) (1 p x, ). : x,=1 : x,=0 Now we observe that given, () = q c(), (6) 2

3 since the functions forces each connected cluster of to have the same spin and there are q possible spins. We can now perform the inner sum in the penultimate displa to see that in fact (with p x, = 1 e βjx, ) [] = [] [,]. So this quantit is the normalization constant for three (probabilit) measures: Potts, random cluster, and Edwards Sokal, and from now on will simpl be denoted. Explicitl, the Potts measure assigns a spin configuration the probabilit µ P () = 1 (δ x, p x, + (1 p x, )), the random cluster measure assigns a bond configurations the probabilit µ RC () = 1 (δ x,,1p x, + δ x,,0(1 p x, )) and the Edwards Sokal measure assigns a spin bond configuration (, ) the probabilit µ ES (, ) = 1 (, ) (δ x,,1p x, + δ x,,0(1 p x, )). (7) B (6), if we fix and sum (7) over all spins, then we get exactl µ RC () (so the marginal distribution of the bond variables is the RC model). Similarl, if we fix and sum (7) over all bonds (putting in δ x, and erasing (, )), then we get 1 (1 p x,) + p x, δ x, ), which b (3) is exactl µ P () (so the marginal distribution of the spin variables is the Potts model). Now let s find the conditional distributions. Let g() be a (dumm) spin observable. We compute the expected value of g(), g() Λ f, with free boundar conditions in Λ. We would like to write it as a nested expectation E [E (g() )]. Writing out g() λ f and expanding as before using (3), 3

4 we have g() Λ f = 1 e βh() g() = 1 [δ x, p x, + (1 p x, ] g() = 1 () [δ x,,1 p x, + δ x,,0 (1 p x, )] g() = 1 [δ x,,1 p x, + δ x,,0 (1 p x, )] (, )g() = 1 q c() [δ x,,1 p x, + δ x,,0 (1 p x, )] (, )g() (, ) = ( ) (, ) µ RC (), (, )g() where we have used (6) ( (, ) = qc() ). So we see that µ ES ( ) = 1 (, ), qc() that is, given, the conditional measure concentrates on s which are compatible with, assigning each such probabilit 1/q c(). Equivalentl, given, the clusters of are labeled 1,..., q with uniform probabilit. Here is a quick application of this: consider q = 2 (the Ising case). Then x Λ f = ( ) (, ) µ RC () 2 c() (2δ x, 1). If x and and in the same connected component of, then x = and the contribution of the inner sum is 1. If x and are in different connected components, then with probabilit 1/2 (so half the contributing configurations) the will be labeled the same spin, in which case the inner sum contributes -1/2. Therefore the net contribution from bond configurations with x and in different clusters is 0, and we conclude x Λ f = µ RC(x ). 4

5 This formula can be generalized to integer values of q, if we take the tetrahedron representation, i.e. we represent the spin variables as unit vectors pointing to vertices of a (q 1) dimensional tetrahedron, so that we have x = (qδ x, 1)/(q 1). Now let g() be a (dumm) bond observable. Then we ma reverse expand using (3) and (6) and get g() = 1 q c() [δ x,,1 p x, + δ x,,0 (1 p x, )] g() = 1 (, ) [δ x, δ x,,1 p x, + δ x,,0 (1 p x, )] g() = 1 (, ) δ x, p x, (1 p x, ) g() : x,=1 : x,=0 = 1 x, [p x,δ x, + (1 p x, )] (, ) x, [δ x,,1 p x, + δ x,,0 (1 p x, )] x, [δ x,,1 p x, + δ x,,0 (1 p x, )] (, )g() = e βh() ( (, ) x, [δ ) x,,1 p x, + δ x,,0 (1 p x, )] (, ) x, [δ x,,1 p x, + δ x,,0 (1 p x, )] g(). So we see that µ ES ( ) = (, ) x, [δ x,,1 p x, + δ x,,0 (1 p x, )] (, ) x, [δ x,,1 p x, + δ x,,0 (1 p x, )]. Again the conditional measure is concentrated on bond configurations which are compatible with : Given, bonds are placed between x and with probabilit p x, if the have the same spin. References [1] M. Aizenman, J. T. Chaes, L. Chaes, and C. M. Newman. Discontinuit of the Magnetization in One-Dimensional 1/ x 2 Ising and Potts Models. J. Stat. Phs. 77, (1994). 5

Graphical Representations and Cluster Algorithms

Graphical 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 information

(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?

(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T? 1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction

More information

Advanced Topics in Probability

Advanced Topics in Probability Advanced Topics in Probability Conformal Methods in 2D Statistical Mechanics Pierre Nolin Different lattices discrete models on lattices, in two dimensions: square lattice Z 2 : simplest one triangular

More information

Critical Dynamics of Two-Replica Cluster Algorithms

Critical Dynamics of Two-Replica Cluster Algorithms University of Massachusetts Amherst From the SelectedWorks of Jonathan Machta 2001 Critical Dynamics of Two-Replica Cluster Algorithms X. N. Li Jonathan Machta, University of Massachusetts Amherst Available

More information

MATH LECTURE NOTES FIRST ORDER SEPARABLE DIFFERENTIAL EQUATIONS OVERVIEW

MATH LECTURE NOTES FIRST ORDER SEPARABLE DIFFERENTIAL EQUATIONS OVERVIEW MATH 234 - LECTURE NOTES FIRST ORDER SEPARABLE DIFFERENTIAL EQUATIONS OVERVIEW Now will will begin with the process of learning how to solve differential equations. We will learn different techniques for

More information

A = {(x, u) : 0 u f(x)},

A = {(x, u) : 0 u f(x)}, Draw x uniformly from the region {x : f(x) u }. Markov Chain Monte Carlo Lecture 5 Slice sampler: Suppose that one is interested in sampling from a density f(x), x X. Recall that sampling x f(x) is equivalent

More information

Statistical Thermodynamics Solution Exercise 8 HS Solution Exercise 8

Statistical 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 information

The 1+1-dimensional Ising model

The 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 information

8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization

8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization 8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization This problem set is partly intended to introduce the transfer matrix method, which is used

More information

Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution

Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution Journal of Physics: Conference Series OPEN ACCESS Ising Model with Competing Interactions on Cayley Tree of Order 4: An Analytic Solution To cite this article: Rukiah bin Ali et al 2013 J. Phys.: Conf.

More information

Lecture 7. The Cluster Expansion Lemma

Lecture 7. The Cluster Expansion Lemma Stanford Universit Spring 206 Math 233: Non-constructive methods in combinatorics Instructor: Jan Vondrák Lecture date: Apr 3, 206 Scribe: Erik Bates Lecture 7. The Cluster Expansion Lemma We have seen

More information

Antiferromagnetic Potts models and random colorings

Antiferromagnetic Potts models and random colorings Antiferromagnetic Potts models and random colorings of planar graphs. joint with A.D. Sokal (New York) and R. Kotecký (Prague) October 9, 0 Gibbs measures Let G = (V, E) be a finite graph and let S be

More information

The near-critical planar Ising Random Cluster model

The near-critical planar Ising Random Cluster model The near-critical planar Ising Random Cluster model Gábor Pete http://www.math.bme.hu/ gabor Joint work with and Hugo Duminil-Copin (Université de Genève) Christophe Garban (ENS Lyon, CNRS) arxiv:1111.0144

More information

Kyle Reing University of Southern California April 18, 2018

Kyle Reing University of Southern California April 18, 2018 Renormalization Group and Information Theory Kyle Reing University of Southern California April 18, 2018 Overview Renormalization Group Overview Information Theoretic Preliminaries Real Space Mutual Information

More information

The Statistical Mechanics of Belief Sharing in Multi Agent Systems

The Statistical Mechanics of Belief Sharing in Multi Agent Systems The Statistical Mechanics of Belief Sharing in Multi Agent Systems Robin Glinton, Katia Sycara, David Scerri, Paul Scerri Robotics Institute Carnegie Mellon University Pittsburgh, PA, U.S.A. (rglinton,

More information

Introduction to the Renormalization Group

Introduction to the Renormalization Group Introduction to the Renormalization Group Gregory Petropoulos University of Colorado Boulder March 4, 2015 1 / 17 Summary Flavor of Statistical Physics Universality / Critical Exponents Ising Model Renormalization

More information

3.320 Lecture 18 (4/12/05)

3.320 Lecture 18 (4/12/05) 3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,

More information

Critical exponents of two-dimensional Potts and bond percolation models

Critical exponents of two-dimensional Potts and bond percolation models J. Phys. A: Math. Gen. 14 (1981) L45-L49. Printed in Great Britain LETTER TO THE EDITOR Critical exponents of two-dimensional Potts and bond percolation models H W J Blotet, M P Nightingale? and B DerridaS

More information

3 Statistical physics models

3 Statistical physics models 3 Statistical physics models 3.1 Percolation In Werner s St. Flour lectures he discusses percolation in the first section and then in more detail in section 10. In the article by Kager and Nienhuis percolation

More information

Decimation Technique on Sierpinski Gasket in External Magnetic Field

Decimation Technique on Sierpinski Gasket in External Magnetic Field Egypt.. Solids, Vol. (), No. (), (009 ) 5 Decimation Technique on Sierpinski Gasket in External Magnetic Field Khalid Bannora, G. Ismail and M. Abu Zeid ) Mathematics Department, Faculty of Science, Zagazig

More information

Intro. 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. 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 information

8.334: Statistical Mechanics II Spring 2014 Test 3 Review Problems

8.334: Statistical Mechanics II Spring 2014 Test 3 Review Problems 8.334: Statistical Mechanics II Spring 014 Test 3 Review Problems The test is closed book, but if you wish you may bring a one-sided sheet of formulas. The intent of this sheet is as a reminder of important

More information

Classical behavior of magnetic dipole vector. P. J. Grandinetti

Classical behavior of magnetic dipole vector. P. J. Grandinetti Classical behavior of magnetic dipole vector Z μ Y X Z μ Y X Quantum behavior of magnetic dipole vector Random sample of spin 1/2 nuclei measure μ z μ z = + γ h/2 group μ z = γ h/2 group Quantum behavior

More information

On Markov Chain Monte Carlo

On Markov Chain Monte Carlo MCMC 0 On Markov Chain Monte Carlo Yevgeniy Kovchegov Oregon State University MCMC 1 Metropolis-Hastings algorithm. Goal: simulating an Ω-valued random variable distributed according to a given probability

More information

VI. Series Expansions

VI. Series Expansions VI. Series Expansions VI.A Low-temperature expansions Lattice models can also be studied by series expansions. Such expansions start with certain exactly solvable limits, and typically represent perturbations

More information

Quadratic Forms. Ferdinand Vanmaele. July 5, 2017

Quadratic Forms. Ferdinand Vanmaele. July 5, 2017 Quadratic Forms Ferdinand Vanmaele July 5, 2017 1 Introduction Question. Consider Pythagorean triples (, y, z), that is integers for which 2 + y 2 = z 2. More generally, we wish to find all integers n

More information

Corrections. ir j r mod n 2 ). (i r j r ) mod n should be d. p. 83, l. 2 [p. 80, l. 5]: d. p. 92, l. 1 [p. 87, l. 2]: Λ Z d should be Λ Z d.

Corrections. ir j r mod n 2 ). (i r j r ) mod n should be d. p. 83, l. 2 [p. 80, l. 5]: d. p. 92, l. 1 [p. 87, l. 2]: Λ Z d should be Λ Z d. Corrections We list here typos and errors that have been found in the published version of the book. We welcome any additional corrections you may find! We indicate page and line numbers for the published

More information

Graphical Representations for Ising Systems in External Fields

Graphical Representations for Ising Systems in External Fields Journal of Statistical Physics. Vol. 93, Nos. 1/2, 1998 Graphical Representations for Ising Systems in External Fields L. Chayes, 1 J. Machta, 2 and O. Redner 3 Received May 8. 1998; final June 26. 1998

More information

Meron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations

Meron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations Meron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations Uwe-Jens Wiese Bern University IPAM Workshop QS2009, January 26, 2009 Collaborators: B. B. Beard

More information

Quantum mechanics in one hour

Quantum mechanics in one hour Chapter 2 Quantum mechanics in one hour 2.1 Introduction The purpose of this chapter is to refresh your knowledge of quantum mechanics and to establish notation. Depending on your background you might

More information

CHAPTER 4. Cluster expansions

CHAPTER 4. Cluster expansions CHAPTER 4 Cluster expansions The method of cluster expansions allows to write the grand-canonical thermodynamic potential as a convergent perturbation series, where the small parameter is related to the

More information

Exam TFY4230 Statistical Physics kl Wednesday 01. June 2016

Exam TFY4230 Statistical Physics kl Wednesday 01. June 2016 TFY423 1. June 216 Side 1 av 5 Exam TFY423 Statistical Physics l 9. - 13. Wednesday 1. June 216 Problem 1. Ising ring (Points: 1+1+1 = 3) A system of Ising spins σ i = ±1 on a ring with periodic boundary

More information

Percolation and the Potts Model

Percolation and the Potts Model Journal of Statistical Physics, Vol. 18, No. 2, 1978 Percolation and the Potts Model F. Y. Wu 1 Received July 13, 1977 The Kasteleyn-Fortuin formulation of bond percolation as a lattice statistical model

More information

Critical Behavior for 2D Uniform and Disordered Ferromagnets at Self-Dual Points

Critical Behavior for 2D Uniform and Disordered Ferromagnets at Self-Dual Points Commun. Math. Phys. 204, 353 366 (1999) Communications in Mathematical Physics Springer-Verlag 1999 Critical Behavior for 2D Uniform and Disordered Ferromagnets at Self-Dual Points L. Chayes 1, K. Shtengel

More information

Spin liquids in frustrated magnets

Spin liquids in frustrated magnets May 20, 2010 Contents 1 Frustration 2 3 4 Exotic excitations 5 Frustration The presence of competing forces that cannot be simultaneously satisfied. Heisenberg-Hamiltonian H = 1 J ij S i S j 2 ij The ground

More information

Basic Physical Chemistry Lecture 2. Keisuke Goda Summer Semester 2015

Basic Physical Chemistry Lecture 2. Keisuke Goda Summer Semester 2015 Basic Physical Chemistry Lecture 2 Keisuke Goda Summer Semester 2015 Lecture schedule Since we only have three lectures, let s focus on a few important topics of quantum chemistry and structural chemistry

More information

arxiv: Statistical mechanical models for quantum codes subject to correlated noise

arxiv: Statistical mechanical models for quantum codes subject to correlated noise Statistical mechanical models for quantum codes subject to correlated noise Christopher Chubb (EQUS USyd), Steve Flammia (EQUS USyd/YQI Yale) (me@)christopherchubb.com arxiv:1809.10704 UNM 2018-10 Quantum

More information

Potts And XY, Together At Last

Potts And XY, Together At Last Potts And XY, Together At Last Daniel Kolodrubetz Massachusetts Institute of Technology, Center for Theoretical Physics (Dated: May 16, 212) We investigate the behavior of an XY model coupled multiplicatively

More information

MEMORIAL UNIVERSITY OF NEWFOUNDLAND

MEMORIAL UNIVERSITY OF NEWFOUNDLAND MEMORIAL UNIVERSITY OF NEWFOUNDLAND DEPARTMENT OF MATHEMATICS AND STATISTICS Volumes Math 11 Winter 17 SOLUTIONS 1. (a) i. The axis of smmetr is a horizontal line, so we integrate with respect to x. The

More information

Non-standard finite-size scaling at first-order phase transitions

Non-standard finite-size scaling at first-order phase transitions Non-standard finite-size scaling at first-order phase transitions Marco Mueller, Wolfhard Janke, Des Johnston Coventry MECO39, April 2014 Mueller, Janke, Johnston Non-Standard First Order 1/24 Plan of

More information

A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ M. H. KALOS J. L. LEBOWITZ

A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ M. H. KALOS J. L. LEBOWITZ JOURNAL OF COMPUTATIONAL PHYSICS 17, lo-18 (1975) A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ Berfer Graduate School of Science, Yeshiva University, New York, New York

More information

POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS

POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS POSITIVE AND NEGATIVE CORRELATIONS FOR CONDITIONAL ISING DISTRIBUTIONS CAMILLO CAMMAROTA Abstract. In the Ising model at zero external field with ferromagnetic first neighbors interaction the Gibbs measure

More information

Introduction to Functions

Introduction to Functions Introduction to Functions A relation is a correspondence between two sets that assigns to each element of the first set an element of the second set. For us, we will consider a relation to be a set of

More information

Random geometric analysis of the 2d Ising model

Random geometric analysis of the 2d Ising model Random geometric analysis of the 2d Ising model Hans-Otto Georgii Bologna, February 2001 Plan: Foundations 1. Gibbs measures 2. Stochastic order 3. Percolation The Ising model 4. Random clusters and phase

More information

The Quantum Harmonic Oscillator

The Quantum Harmonic Oscillator The Classical Analysis Recall the mass-spring system where we first introduced unforced harmonic motion. The DE that describes the system is: where: Note that throughout this discussion the variables =

More information

Electronic Structure Models

Electronic Structure Models Electronic Structure Models Hückel Model (1933) Basic Assumptions: (a) One orbital per atom contributes to the basis set; all orbitals "equal" (b) The relevant integrals involving the Hamiltonian are α

More information

11.1 Double Riemann Sums and Double Integrals over Rectangles

11.1 Double Riemann Sums and Double Integrals over Rectangles Chapter 11 Multiple Integrals 11.1 ouble Riemann Sums and ouble Integrals over Rectangles Motivating Questions In this section, we strive to understand the ideas generated b the following important questions:

More information

The Random-Cluster Model

The Random-Cluster Model Geoffrey Grimmett The Random-Cluster Model With 37 Figures Springer Geoffrey Grimmett Statistical Laboratory Centre for Mathematical Sciences University of Cambridge Wilberforce Road Cambridge CB3 0WB

More information

Universal phase transitions in Topological lattice models

Universal phase transitions in Topological lattice models Universal phase transitions in Topological lattice models F. J. Burnell Collaborators: J. Slingerland S. H. Simon September 2, 2010 Overview Matter: classified by orders Symmetry Breaking (Ferromagnet)

More information

Brazilian Journal of Physics ISSN: Sociedade Brasileira de Física Brasil

Brazilian Journal of Physics ISSN: Sociedade Brasileira de Física Brasil Brazilian Journal of Physics ISSN: 0103-9733 luizno.bjp@gmail.com Sociedade Brasileira de Física Brasil Braga, Gastão A.; Sá Barreto, F. C.; Santos, Jander P. Upper Bounds on the Critical Temperature of

More information

Simulation of the two-dimensional square-lattice Lenz-Ising model in Python

Simulation 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 information

Assignment 6: MCMC Simulation and Review Problems

Assignment 6: MCMC Simulation and Review Problems Massachusetts Institute of Technology MITES 2018 Physics III Assignment 6: MCMC Simulation and Review Problems Due Monday July 30 at 11:59PM in Instructor s Inbox Preface: In this assignment, you will

More information

Degeneracy & in particular to Hydrogen atom

Degeneracy & in particular to Hydrogen atom Degeneracy & in particular to Hydrogen atom In quantum mechanics, an energy level is said to be degenerate if it corresponds to two or more different measurable states of a quantum system. Conversely,

More information

ON THE ISING MODEL WITH RANDOM BOUNDARY CONDITION

ON THE ISING MODEL WITH RANDOM BOUNDARY CONDITION ON THE ISING MODEL WITH RANDOM BOUNDARY CONDITION A. C. D. VAN ENTER, K. NETOČNÝ, AND H. G. SCHAAP ABSTRACT. The infinite-volume limit behavior of the 2d Ising model under possibly strong random boundary

More information

0. Construction of d-dimensional Isotropic and Anisotropic Hierarchical Lattices 1. Frustrated Systems and Chaotic Renormalization- Group

0. Construction of d-dimensional Isotropic and Anisotropic Hierarchical Lattices 1. Frustrated Systems and Chaotic Renormalization- Group Hierarchical Lattices: Renormalization-Group Solutions of Plain, Anisotropic,Chaotic, Heterogeneous, and Clustered Systems Collaborators: D. Andelman, A. Erbaş, A. Falicov, K. Hui, M. Hinczewski, A. Kabakçıoğlu,

More information

Metropolis Monte Carlo simulation of the Ising Model

Metropolis 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 information

8 Potts model. 2. These two actions commute and each action generates the full centraliser

8 Potts model. 2. These two actions commute and each action generates the full centraliser 8 Potts model When writing the integrable R-matrix of the six-vertex model, in (6.18), we have briefly come across a new type of algebraic structure: the Temperley- Lieb algebra. This is an example of

More information

Nematic phase of spin 1 quantum systems

Nematic phase of spin 1 quantum systems Nematic phase of spin 1 quantum systems Daniel Ueltschi Department of Mathematics, University of Warwick International Conference on Applied Mathematics, Hράκλειo, 17 September 2013 D. Ueltschi (Univ.

More information

S j H o = gµ o H o. j=1

S 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 information

Translation invariant Gibbs states for the Ising model

Translation invariant Gibbs states for the Ising model Translation invariant Gibbs states for the Ising model Thierry Bodineau To cite this version: Thierry Bodineau. Translation invariant Gibbs states for the Ising model. 2004. HAL Id: hal-00003038

More information

2D Critical Systems, Fractals and SLE

2D Critical Systems, Fractals and SLE 2D Critical Systems, Fractals and SLE Meik Hellmund Leipzig University, Institute of Mathematics Statistical models, clusters, loops Fractal dimensions Stochastic/Schramm Loewner evolution (SLE) Outlook

More information

INTRODUCTION TO STATISTICAL MECHANICS Exercises

INTRODUCTION TO STATISTICAL MECHANICS Exercises École d hiver de Probabilités Semaine 1 : Mécanique statistique de l équilibre CIRM - Marseille 4-8 Février 2012 INTRODUCTION TO STATISTICAL MECHANICS Exercises Course : Yvan Velenik, Exercises : Loren

More information

Kempe reconfiguration and Potts antiferromagnets

Kempe reconfiguration and Potts antiferromagnets Kempe reconfiguration and Potts antiferromagnets Jesús Salas Grupo de Teorías de Campos y Física Estadística Instituto Gregorio Millán Barbany, UC3M Unidad Asociada al Instituto de Estructura de la Materia,

More information

3 Undirected Graphical Models

3 Undirected Graphical Models Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.438 Algorithms For Inference Fall 2014 3 Undirected Graphical Models In this lecture, we discuss undirected

More information

Name:... Section:... Physics 208 Quiz 8. April 11, 2008; due April 18, 2008

Name:... Section:... Physics 208 Quiz 8. April 11, 2008; due April 18, 2008 Name:... Section:... Problem 1 (6 Points) Physics 8 Quiz 8 April 11, 8; due April 18, 8 Consider the AC circuit consisting of an AC voltage in series with a coil of self-inductance,, and a capacitor of

More information

Complex Systems Methods 9. Critical Phenomena: The Renormalization Group

Complex 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 information

6. This sum can be rewritten as 4( ). We then recall the formula n =

6. This sum can be rewritten as 4( ). We then recall the formula n = . c = 9b = 3 b = 3 a 3 = a = = 6.. (3,, ) = 3 + + 3 = 9 + + 3 = 6 6. 3. We see that this is equal to. 3 = ( +.) 3. Using the fact that (x + ) 3 = x 3 + 3x + 3x + and replacing x with., we find that. 3

More information

. (70.1) r r. / r. Substituting, we have the following equation for f:

. (70.1) r r. / r. Substituting, we have the following equation for f: 7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical

More information

Local Maximums and Local Minimums of Functions. f(x, y) has a local minimum at the point (x 0, y 0 ) if f(x 0, y 0 ) f(x, y) for

Local Maximums and Local Minimums of Functions. f(x, y) has a local minimum at the point (x 0, y 0 ) if f(x 0, y 0 ) f(x, y) for Local Extrema Previousl we have taken the partial derivative of a function f(x, ). But those partial derivatives were themselves functions and so we can take their partial derivatives. Local Maximums and

More information

VI.D Self Duality in the Two Dimensional Ising Model

VI.D Self Duality in the Two Dimensional Ising Model VI.D Self Duality in the Two Dimensional Ising Model Kramers and Wannier discovered a hidden symmetry that relates the properties of the Ising model on the square lattice at low and high temperatures.

More information

INTRODUCTION TO SIGMA NOTATION

INTRODUCTION TO SIGMA NOTATION INTRODUCTION TO SIGMA NOTATION The notation itself Sigma notation is a way of writing a sum of many terms, in a concise form A sum in sigma notation looks something like this: 3k The Σ sigma) indicates

More information

Random even graphs and the Ising model

Random even graphs and the Ising model Random even graphs and the Ising model Geoffrey Grimmett, Svante Janson September 19, 007 Abstract We explore the relationship between the Ising model with inverse temperature β, the q = random-cluster

More information

INF Introduction to classifiction Anne Solberg

INF Introduction to classifiction Anne Solberg INF 4300 8.09.17 Introduction to classifiction Anne Solberg anne@ifi.uio.no Introduction to classification Based on handout from Pattern Recognition b Theodoridis, available after the lecture INF 4300

More information

A THREE-COMPONENT MOLECULAR MODEL WITH BONDING THREE-BODY INTERACTIONS

A THREE-COMPONENT MOLECULAR MODEL WITH BONDING THREE-BODY INTERACTIONS STATISTICAL PHYSICS, SOLID STATE PHYSICS A THREE-COMPONENT MOLECULAR MODEL WITH BONDING THREE-BODY INTERACTIONS FLORIN D. BUZATU Department of Theoretical Physics, National Institute of Physics and Nuclear

More information

Macroscopic Degeneracy and FSS at 1st Order Phase Transitions

Macroscopic Degeneracy and FSS at 1st Order Phase Transitions Macroscopic Degeneracy and FSS at 1st Order Phase Transitions Marco Mueller (hero), Wolfhard Janke (good guy), Des Johnston (villain) Krakow, Oct 2015 Mueller, Janke, Johnston Degeneracy/FSS 1/22 Plan

More information

Reteaching (continued)

Reteaching (continued) Quadratic Functions and Transformations If a, the graph is a stretch or compression of the parent function b a factor of 0 a 0. 0 0 0 0 0 a a 7 The graph is a vertical The graph is a vertical compression

More information

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation. Review 1 1) State the largest possible domain o deinition or the unction (, ) = 3 - ) Determine the largest set o points in the -plane on which (, ) = sin-1( - ) deines a continuous unction 3) Find the

More information

UNIT #8 QUADRATIC FUNCTIONS AND THEIR ALGEBRA REVIEW QUESTIONS

UNIT #8 QUADRATIC FUNCTIONS AND THEIR ALGEBRA REVIEW QUESTIONS Answer Ke Name: Date: UNIT #8 QUADRATIC FUNCTIONS AND THEIR ALGEBRA REVIEW QUESTIONS Part I Questions. For the quadratic function shown below, the coordinates of its verte are, (), 7 6,, 6 The verte is

More information

Read carefully the instructions on the answer book and make sure that the particulars required are entered on each answer book.

Read carefully the instructions on the answer book and make sure that the particulars required are entered on each answer book. THE UNIVERSITY OF WARWICK FOURTH YEAR EXAMINATION: APRIL 2007 INTRODUCTION TO STATISTICAL MECHANICS Time Allowed: 3 hours Read carefully the instructions on the answer book and make sure that the particulars

More information

arxiv: v2 [math.pr] 26 Aug 2017

arxiv: v2 [math.pr] 26 Aug 2017 CONSTRAINED PERCOLATION, ISING MODEL AND XOR ISING MODEL ON PLANAR LATTICES ZHONGYANG LI arxiv:1707.04183v2 [math.pr] 26 Aug 2017 Abstract. We study constrained percolation models on planar lattices including

More information

Decision tables and decision spaces

Decision tables and decision spaces Abstract Decision tables and decision spaces Z. Pawlak 1 Abstract. In this paper an Euclidean space, called a decision space is associated with ever decision table. This can be viewed as a generalization

More information

Lesson 29. Exact Differential Equations: Integrating Factors

Lesson 29. Exact Differential Equations: Integrating Factors Module 3: Ordinar Differential Equations Lesson 9 Eact Differential Equations: Integrating Factors In general, equations of the tpe M,)d + N,)d = 0 are not eact. However, it is sometimes possible to transform

More information

Bond Dilution Effects on Bethe Lattice the Spin-1 Blume Capel Model

Bond Dilution Effects on Bethe Lattice the Spin-1 Blume Capel Model Commun. Theor. Phys. 68 (2017) 361 365 Vol. 68, No. 3, September 1, 2017 Bond Dilution Effects on Bethe Lattice the Spin-1 Blume Capel Model Erhan Albayrak Erciyes University, Department of Physics, 38039,

More information

Layering in the SOS Model Without External Fields

Layering in the SOS Model Without External Fields Layering in the SOS Model Without External Fields Ken Alexander Univ. of Southern California Coauthors: François Dunlop (Université de Cergy-Pontoise) Salvador Miracle-Solé (CPT, Marseille) March 2011

More information

Magnetoelastics in the frustrated spinel ZnCr 2 O 4. Roderich Moessner CNRS and ENS Paris

Magnetoelastics in the frustrated spinel ZnCr 2 O 4. Roderich Moessner CNRS and ENS Paris Magnetoelastics in the frustrated spinel ZnCr 2 O 4 Roderich Moessner CNRS and ENS Paris CEA Saclay, June 2005 Overview Neutron scattering experiments on ZnCr 2 O 4 Modelling magnetism on the B sublattice

More information

Cluster Expansion in the Canonical Ensemble

Cluster Expansion in the Canonical Ensemble Università di Roma Tre September 2, 2011 (joint work with Dimitrios Tsagkarogiannis) Background Motivation: phase transitions for a system of particles in the continuum interacting via a Kac potential,

More information

PH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5

PH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 Spring 2009 POP QUIZ

More information

1 2n. i=1. 2m+1 i = m, while. 2m + 1 = 1. is the best possible for any odd n.

1 2n. i=1. 2m+1 i = m, while. 2m + 1 = 1. is the best possible for any odd n. Problem 1. Let n be a positive integer. Find the largest nonnegative real number f(n) (depending on n) with the following property: whenever a 1, a,..., a n are real numbers such that a 1 + a + + a n is

More information

i=1 h n (ˆθ n ) = 0. (2)

i=1 h n (ˆθ n ) = 0. (2) Stat 8112 Lecture Notes Unbiased Estimating Equations Charles J. Geyer April 29, 2012 1 Introduction In this handout we generalize the notion of maximum likelihood estimation to solution of unbiased estimating

More information

Martin Schnabl. Institute of Physics AS CR. Collaborators: T. Kojita, M. Kudrna, C. Maccaferri, T. Masuda and M. Rapčák

Martin Schnabl. Institute of Physics AS CR. Collaborators: T. Kojita, M. Kudrna, C. Maccaferri, T. Masuda and M. Rapčák Martin Schnabl Collaborators: T. Kojita, M. Kudrna, C. Maccaferri, T. Masuda and M. Rapčák Institute of Physics AS CR 36th Winter School Geometry and Physics, Srní, January 22nd, 2016 2d Conformal Field

More information

Continuity, End Behavior, and Limits. Unit 1 Lesson 3

Continuity, End Behavior, and Limits. Unit 1 Lesson 3 Unit Lesson 3 Students will be able to: Interpret ke features of graphs and tables in terms of the quantities, and sketch graphs showing ke features given a verbal description of the relationship. Ke Vocabular:

More information

A lower bound on the two arms exponent for critical percolation on the lattice

A lower bound on the two arms exponent for critical percolation on the lattice A lower bound on the two arms exponent for critical percolation on the lattice Raphaël Cerf Université Paris Sud and IUF June 13, 2013 Abstract We consider the standard site percolation model on the d

More information

NAME: KEY CHEMISTRY 419, Spring 2017(17S) Section Number: 10 Midterm Examination #2, May 4, 2017

NAME: KEY CHEMISTRY 419, Spring 2017(17S) Section Number: 10 Midterm Examination #2, May 4, 2017 NAME: KEY CHEMISTRY 419, Spring 017(17S) Section Number: 10 Midterm Examination #, Ma 4, 017 Answer each question in the space provided; use back of page if extra space is needed. Answer questions so the

More information

Stochastic Interpretation for the Arimoto Algorithm

Stochastic Interpretation for the Arimoto Algorithm Stochastic Interpretation for the Arimoto Algorithm Serge Tridenski EE - Sstems Department Tel Aviv Universit Tel Aviv, Israel Email: sergetr@post.tau.ac.il Ram Zamir EE - Sstems Department Tel Aviv Universit

More information

2.8 Implicit Differentiation

2.8 Implicit Differentiation .8 Implicit Differentiation Section.8 Notes Page 1 Before I tell ou what implicit differentiation is, let s start with an example: EXAMPLE: Find if x. This question is asking us to find the derivative

More information

Phase transitions beyond the Landau-Ginzburg theory

Phase transitions beyond the Landau-Ginzburg theory Phase transitions beyond the Landau-Ginzburg theory Yifei Shi 21 October 2014 1 Phase transitions and critical points 2 Laudau-Ginzburg theory 3 KT transition and vortices 4 Phase transitions beyond Laudau-Ginzburg

More information

CHAPTER EIGHT: SOLVING QUADRATIC EQUATIONS Review April 9 Test April 17 The most important equations at this level of mathematics are quadratic

CHAPTER EIGHT: SOLVING QUADRATIC EQUATIONS Review April 9 Test April 17 The most important equations at this level of mathematics are quadratic CHAPTER EIGHT: SOLVING QUADRATIC EQUATIONS Review April 9 Test April 17 The most important equations at this level of mathematics are quadratic equations. They can be solved using a graph, a perfect square,

More information

Comparing and Improving Quark Models for the Triply Bottom Baryon Spectrum

Comparing and Improving Quark Models for the Triply Bottom Baryon Spectrum Comparing and Improving Quark Models for the Triply Bottom Baryon Spectrum A thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Science degree in Physics from the

More information

X-ray Spectroscopy Theory Lectures

X-ray Spectroscopy Theory Lectures TIMES Lecture Series SIMES-SLAC-Stanford Winter, 2017 X-ray Spectroscopy Theory Lectures J. J. Rehr I. Introduction to the Theory of X-ray spectra II. Real-space Green's function Theory and FEFF III. Inelastic

More information

Series expansions from the corner transfer matrix renormalization group method

Series expansions from the corner transfer matrix renormalization group method Series expansions from the corner transfer matrix renormalization group method 1 Andrew Rechnitzer 2 1 LaBRI/The University of Melbourne 2 University of British Columbia January 27, 2011 What is the CTMRG

More information