Equilibrium, out of equilibrium and consequences
|
|
- Myrtle Beasley
- 5 years ago
- Views:
Transcription
1 Equilibrium, out of equilibrium and consequences Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland SF-MTPT Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 1 / 63
2 Agenda 1 Equilibrium and H-Boltzmann theorem 2 Phase transitions in 1D 3 Example: Zero temperature Glauber dynamics in 1D 4 Markov chains and Monte Carlo simulations 5 Spin dynamics Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 2 / 63
3 H-Boltzmann theorem Microscopic reversibility P r (t) - probability of finding a system in the state r (one of possible states) at time t A system is certainly in one of possible states, i.e. r P r(t) = 1 W r s W rs - transition probability per unit time from state r to s Evolution equation (balance): dp r dt = P s W s r P r W r s s r s r }{{}}{{} s r r s If there is a symmetry W sr = W rs (microscopic reversibility): (1) dp r dt = s r W rs (P s P r ) (2) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 3 / 63
4 What is H? Let s define H as: H-Boltzmann theorem H < ln P r > r P r ln P r, (3) therefore: dh dt = d dt = r = r ( ) P r ln P r r ( ln P r dp r ) dt + P 1 dp r r P r dt ) ( dp r ln P r dt + dp r dt = r dp r dt (ln P r + 1) (4) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 4 / 63
5 Evolution of H H-Boltzmann theorem dh dt = r dp r dt (ln P r + 1), (5) but: dp r dt = s W rs (P s P r ), (6) and finally: dh dt = r W rs (P s P r ) (ln P r + 1), (7) s Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 5 / 63
6 H-Boltzmann theorem On the other hand we can replace r by s: H = P s ln P s (8) s and we obtain: dh dt = s Recalling that also: dh dt = r we can write: dh dt W sr (P r P s ) (ln P s + 1). (9) r W rs (P s P r ) (ln P r + 1), (10) s = 1 W rs (P s P r ) (ln P r + 1) 2 r s + 1 W sr (P r P s ) (ln P s + 1) (11) 2 s r Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 6 / 63
7 H-Boltzmann theorem Use the symmetry... dh dt = 1 W rs (P s P r ) (ln P r + 1) 2 r s + 1 W sr (P r P s ) (ln P s + 1) (12) 2 s r We use the microscopic reversibility (symmetry) condition W rs = W sr : dh = 1 W rs (P r P s ) (ln P r + 1) dt 2 r s + 1 W rs (P r P s ) (ln P s + 1) 2 s r = 1 W rs (P r P s ) (ln P r ln P s ). (13) 2 r s Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 7 / 63
8 H-Boltzmann theorem Evolution of H and symmetry dh dt = 1 W rs (P r P s ) (ln P r ln P s ). (14) 2 r s Note that (P r P s ) (ln P r ln P s ) 0, because: P r > P s then ln P r > ln P s (P r P s ) (ln P r ln P s ) > 0, P r < P s then ln P r < ln P s (P r P s ) (ln P r ln P s ) > 0, P r = P s then ln P r = ln P s (P r P s ) (ln P r ln P s ) = 0. Moreover W rs 0 (probability), therefore: dh dt 0 (= for P s = P r ) (15) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 8 / 63
9 H-Boltzmann theorem H-Boltzmann theorem dh dt 0 (= for P s = P r ) (16) H < ln P r > decreases in time until P s = P r (equal a priori probabilities - equilibrium) What is H? hence: S k B < ln P r > S = k B H, (17) dh dt 0 (= for P s = P r ) ds dt 0 (= for P s = P r ) (18) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 9 / 63
10 H-Boltzmann theorem H-Boltzmann theorem dh ds 0 (= for Ps = Pr ) 0 (= for Ps = Pr ) dt dt Microscopic interpretation of a second law of thermodynamics Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT (19) 10 / 63
11 H-Boltzmann theorem Comments on H-Boltzmann theorem Valid for isolated systems: dp r dt = P s W sr P r W rs +? (20) s s }{{}}{{} s r r s Microscopic reversibility is needed: W rs = W sr Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 11 / 63
12 Detailed balance Master equation and detailed balance P r (t) - probability of finding a system in the state r (one of possible states) at time t A system is certainly in one of possible states, i.e. r P r(t) = 1 W r s W rs - transition probability per unit time from state r to s Master equation: Detailed balance: dp r dt = P s W s r P r W r s s r s r }{{}}{{} s r r s (21) P s W sr = P r W rs dp r dt = 0 (22) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 12 / 63
13 Master equation and detailed balance Steady state and stationarity It is possible that detailed balance is not fulfilled but: dp r = 0 (23) dt Condition for a stationarity (steady state): P s W s r = P r W r s (24) s r s r }{{}}{{} s r r s Steady state in and out of equilibrium: Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 13 / 63
14 Questions: Master equation and detailed balance Consider a system with 3 possible states A, B, C. Transition probabilities w i j 0 for any i j (where i and j denotes one of possible states A, B, C). What are conditions for w i j, to fulfill detailed balance: What are stationarity conditions? P s W sr = P r W rs? (25) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 14 / 63
15 Phase transitions in 1D Diffusion on a circle 1D lattice with periodic boundary conditions Each site is occupied by a particle (A, B, C) C B A B B Random sequential updating: B A AB q BA, BC q CB, CA q AC, BA 1 AB, CB 1 BC, AC 1 CA. A C C B C The number of particles (N A, N B, N C ) is conserved [1] M. R. Evans, Y. Kafri, H. M. Koduvely and D. Mukamel, Phys. Rev. Lett. 80, 425 (1998) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 15 / 63
16 Phase transitions in 1D Phase separation For q = 1 disorder For q < 1 walls AB, BC, CA are stable Domains are created...aabbccaaab... B C B A B B A Domains grow in time ln t/ ln q Thermodynamic limit separtion into 3 domains A C C B C Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 16 / 63
17 Phase transitions in 1D Phase transitions in 1D? - Landau argument Equilibrium minimum of a free energy: F = E TS Change of a free energy (extensive): F = E T S = ɛ T S Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 17 / 63
18 What is entropy? Phase transitions in 1D S = k B ln W L {}}{ AAAABBBB }{{} L 1 Change of a free energy W = L 1 S = k B ln(l 1) F = ɛ T S = ɛ Tk B ln(l 1) In thermodynamic limit, i.e. for L F < 0 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 18 / 63
19 Phase transitions in 1D Landau assumptions Energy of a domain wall ɛ is finite F = ɛ Tk B ln(l 1) Short-range interactions no interactions between domain walls In 1D Ising model with J(r) r 1 σ there is an order in low temperatures T > 0 for σ 1 T > 0 otherwise entropy argument is not important F = E TS Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 19 / 63
20 Phase transitions in 1D Other argument - dynamics of a system AAAABBBBBAAA Transition disorder order annihilation of domain walls Domain ended with two walls costs 2ɛ In 1D there is no energy profit related to shrinking of a domain There is no force that eliminates domain walls (bonds) in 1D Assumption: no long-range interactions between bonds Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 20 / 63
21 Phase transitions in 1D Other suggestions: 1 Exact result for 1D Ising model 2 Van Hove theorem: No phase transitions for in 1D system of hard cores with short-range interactions 3 Mermin-Wagner theorem: No phase transitions in 1D and 2D systems with continuous symmetry for short-range interactions There is no general theorem for non-existence of a phase transitions in 1D! Literature: M. R. Evans, Braz. J. Phys. 30, 42 (2000) J. A. Cuesta, A. Sanchez, J. Stat. Phys. 115, 869 (2004) N. Theodorakopoulos, Physica D 216, 185 (2006) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 21 / 63
22 Phase transitions in 1D L. van Hove, Physica 16 (1950) 137 Free energy has been calculated for a hard-core model with short-range interactions Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 22 / 63
23 Phase transitions in 1D Limitations Identical particles were considered (hard cores of radius d 0 ) Homogeneous system identical elements Integrations on a distance d d 1 Short-range interactions Interaction energy depended only on a distance between particles (and not positions) No heterogeneous field or interactions Transfer matrix method have been used only maximum eigenvalue in thermodynamic limit Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 23 / 63
24 Phase transitions in 1D Equilibrium phase transitions in 1D Heterogeneous systems Long-range interactions Infinite energy of interactions At temperature T = 0 Out of equilibrium? Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 24 / 63
25 Out of equilibrium Phase transitions in 1D H-Boltzmann theorem Isolated system W sr = W rs (micro rev) H < ln P r >= S/k B dh dt 0 (= dla P s = P r ) Detailed balance: P eq s W sr = Pr eq W rs Type I Hermitian Hamiltonian Steady states given by the Gibbs-Boltzmann dist Initial state far from equilibrium In thermodynamic limit equilibrium is never reached Type II No Hamiltonian Dynamics given by transition probabilities No detailed balance Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 25 / 63
26 Example: Zero temperature Glauber dynamics in 1D Ising model in 1D 1D lattice with periodic boundary conditions Each site occupied with σ i = ±1 Interactions described by: H = J L σ i σ i+1 i=1 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 26 / 63
27 Example: Zero temperature Glauber dynamics in 1D Generalized Glauber dynamics at T = 0 Roy J. Glauber, Journal of Mathematical Physics (1963) σ i (t) σ i (t + dt) with prob. 1 ifδe < 0, W (δe) = W 0 if δe = 0, 0 if δe > 0, σ i 1 σ i σ i+1 W W 0 W 0 W 0 W 0 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 27 / 63
28 Example: Zero temperature Glauber dynamics in 1D Updating 1 Sequential consider single spin in an elementary time step dt (Glauber, Metropolis) 2 Synchronous all spins are considered at the same time (as in cellular automata); for W 0 = 1 absorbing antiferromagnetic state: 3 c-synchronous partially synchronous K. Sznajd-Weron and S. Krupa, Phys. Rev. E 74, (2006) sequential }{{} c=1/l synchronous }{{} c=1 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 28 / 63
29 Example: Zero temperature Glauber dynamics in 1D Order parameter and phase diagram ρ = 1 L (1 σ i σ i+1 ) 2L i=1 <ρ st > W c Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 29 / 63
30 Example: Zero temperature Glauber dynamics in 1D Discontinuous phase transitions for c = 1 1 <ρ> W 0 =0.8 W =0.7 0 W 0 =0.6 W =0.5 0 W 0 =0.4 W 0 =0.3 W = t Monte Carlo simulations: time evolution of bonds for the systems size L = 160, averaging over 10 4 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 30 / 63
31 Example: Zero temperature Glauber dynamics in 1D Scaling for c = 1: β = 0, ν = 1 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 31 / 63
32 Example: Zero temperature Glauber dynamics in 1D Discontinuous phase transitions Jump of an order parameter ρ Trivial critical exponents β = 0, ν = 1 Phase coexistence for W 0 = 1 2 Hysteresis cycle metastable states, dependence on the initial conditions 1 1 <ρ st > 0.5 <ρ st > W 0 W 0 Figure : left: c = 1, right: c = 0.95 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 32 / 63
33 Katarzyna Sznajd-Weron Figure (WUT) : c = Equilibrium, 0.99, Wout c of = equilibrium 0.53, and β = , ν = 1 SF-MTPT 33 / 63 Example: Zero temperature Glauber dynamics in 1D Continuous phase transitions for c < 1 <ρ st > W 0 L β/ν <ρ st > 5 4 L=32 3 L=48 2 L=64 L=100 1 L=123 L= (W 0 W c )L 1/ν Figure : c = 0.9, W c = 0.6, β = 0.4, ν = <ρ st > W 0 L β/ν <ρ st > L=32 2 L=48 L=64 1 L=100 L=123 L= (W 0 W c )L 1/ν
34 Example: Zero temperature Glauber dynamics in 1D Critical exponents and the type of transition c W c β ν Discontinuous c = 1, Continuous(W 0 = 1) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 34 / 63
35 Monte Carlo and dynamics Monte Carlo and dynamics Monte Carlo is a method of calculating quantities that can be presented as the expected value of certain probabilistic distributions Let s m is our quantity and we know that m = EX of a certain random variable X If we are able to generate S 1, S 2,... that are iid (independent identically distributed) from distribution of X than the law of large numbers: 1 lim n n (S S n ) = m. (26) Monte Carlo estimate a quantity m by a properly chosen sample of n elements Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 35 / 63
36 Monte Carlo and dynamics How to choose a sample properly in Statistical Physics? System in a given external conditions equilibrium (or approaches equilibrium) Probability density function in equilibrium: P eq (α) p α = 1 Z exp( βe α), (27) where Z = α exp( βe α) partition function Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 36 / 63
37 Monte Carlo and dynamics Example: evolution of a system Consider Ising model Initial condition: m(0) = m 0 = 1 (low temperature state T 0, in general m = ±1). Initial distribution of our random variable m is given by P m (0) = δ m,m0 = δ m,1. What will distribution of m at temperature T = T 1 >> 0? From statistical physics system approaches equilibrium given by the Boltzmann-Gibbs distribution Evolution of m(t) and its distribution In the probability theory such an evolution stochastic processes Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 37 / 63
38 Monte Carlo and dynamics Stochastic processes Stochastic processes evolution of some random value X t, or system, over time t (function) From a statistical physics point of view ensemble of trajectories (realizations of a given process) What can be our random variable? Particularly convenient is a microstate (configuration) of a system. My recomendation: Ch. M. Grinstead and J. L. Snell Introduction to Probability Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 38 / 63
39 Monte Carlo and dynamics Markov process We have a set of states S = s 1, s 2,..., s r The process starts in one of these states and moves successively from one state to another If the chain is currently in state s i, then it moves to state s j at the next step with a probability W ij The probability W ij does not depend upon which states the chain was in before the current state R. A. Howard a frog jumping on a set of lily pads. The frog starts on one of the pads and then jumps from lily pad to lily pad with the appropriate transition probabilities Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 39 / 63
40 Monte Carlo and dynamics Example: Weather in the Land of Oz The Land of Oz is blessed by many things, but not by good weather. They never have two nice days in a row. If they have a nice day, they are just as likely to have snow as rain the next day. If they have snow or rain, they have an even chance of having the same the next day. If there is change from snow or rain, only half of the time is this a change to a nice day. Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 40 / 63
41 Monte Carlo and dynamics Example: Weather in the Land of Oz We form a Markov chain We take as states the kinds of weather: R (rain), N (nice), and S (snow). From the information we determine the transition probabilities (transition matrix): Check it! W = R N S R 1/2 1/4 1/4 N 1/2 0 1/2 S 1/4 1/4 1/2 (28) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 41 / 63
42 Monte Carlo and dynamics Markov chain Let s define Markov process with discrete time t 1 < t 2 <... < t n, for a finite number of states σ 1, σ 2, σ 3,..., σ N : Random variable X t describes a state of a system at t. Conditional probability that X tn = σ in : P(X tn = σ in X tn 1 = σ in 1, X tn 2 = σ in 2,..., X t1 = σ i1 ), (29) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 42 / 63
43 Markov chain Monte Carlo and dynamics Markov property: P(X tn = σ in X tn 1 = σ in 1,..., X t1 = σ i1 ) = P(X tn = σ in X tn 1 = σ in 1 ). (30) Above conditional probability transition probability from i to j therefore: W ij = W (σ i σ j ) = P(X tn = σ j X tn 1 = σ i ). (31) W ij 0, W ij = 1, (32) j Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 43 / 63
44 Monte Carlo and dynamics Master equation Probability that at time t n system will be in a state σ j P(X tn = σ j ) = P(X tn = σ j X tn 1 = σ i )P(X tn 1 = σ i ) (33) = W ij P(X tn 1 = σ i ). (34) Evolution described by the Master Equation: dp(σ j, t) dt = i W ji P(σ j, t) + i W ij P(σ i, t). (35) W ji P eq (σ j ) = W ij P eq (σ i ). (36) P eq (σ j ) = Z 1 e E(σ j )/T. (37) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 44 / 63
45 Spin dynamics Voter model Toy model Each site of a graph (or lattice) is occupied by a particle in one of q states In each elementary time step we choose one spin It takes a state of a randomly chosen nn Dynamics each finite system reaches of q absorbing states (all the same) Exit time (time to reach an absorbing state) depends on the system size, dimension and q Most common q = 2 (spin s = 1/2) s( x) = ±1 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 45 / 63
46 Spin dynamics Voter model - algorithm Set an initial conditions e.g. for each x = 1,..., N with probability p set s(x, 0) = 1, and with 1 p set s(x, 0) = 1 1 Choose randomly (uniform distribution) a spin s(x, t) 2 Choose randomly one of its nn s(y, t) 3 If s(x, t) s(y, t) then s(x, t + t) = s(y, t) 4 Increase time t t + t 5 Check if all spins are the same. If not goto 1 N elementary steps increase time by 1, i.e. t = 1/N Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 46 / 63
47 Spin dynamics Voter model - transition probability Transition rate w x ( s) spin s(x) flips: [ w x ( s) = 1 1 s(x) 2 z hence: s(x, t + t) = y <x> s(y) { s(x, t) with probability 1 wx ( s) t s(x, t) with probability w x ( s) t ] (38) (39) prawdopodobieństwo = 0 prawdopodobieństwo = 1/4 prawdopodobieństwo = 1/2 prawdopodobieństwo = 3/4 prawdopodobieństwo = 1 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 47 / 63
48 Spin dynamics Voter model - master equation (ME) State is fully described by P( s, t) probability that at time t a system has a state s = [s 1, s 2,..., s N ]. Number of states 2 N Probability distribution evolves due to ME dp( s) dt = x w x ( s)p( s) + x w x ( s x )P( s x ), (40) where s x differs from s by spin at site x. Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 48 / 63
49 Spin dynamics Voter model magnetization Spin at site x changes by 2s(x) with a rate w x ( s) w x, therefore average spin at site denoted by S(x) < s(x) >: Recalling that: ds(x) dt w x ( s) = 1 2 = 2 < s(x)w x > (41) [ 1 s(x) z y <x> s(y) ] We obtain (calculate this!): ds(x) dt = S(x) + 1 z y <x> S(y) (42) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 49 / 63
50 Spin dynamics Voter model magnetization ds(x) dt Average magnetization: = S(x) + 1 z y <x> S(y) (43) m x S(x) N (44) From (43) we obtain that on average: What are consequences? dm dt = 0 (45) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 50 / 63
51 Spin dynamics Voter model exit Probability Consider a finite system with the initial probability of +1 equal p, and 1 equal 1 p Therefore initial magnetization m 0 = 2p 1 Eventually system reaches absorbing state m = 1 with probability E(p) or m = 1 z 1 E(p) Final magnetization: m = E(p) 1 + [1 E(p)] ( 1) = 2E(p) 1 (46) Previously, we have obtained that m = m 0 = 2p 1 E(p) = p Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 51 / 63
52 Spin dynamics Glauber dynamics for an Ising model H = <i,j> J ij S i S j S i = ±1 (47) We choose in each time step single spin and S i S i with probability that depends on the energy change energia maleje energia rośnie brak zmiany energii Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 52 / 63
53 Spin dynamics Glauber dynamics - transition rates (TR) Detailed balance has to be fulfilled: P eq ( s)w i ( s) = P eq ( s i )w i ( s i ), (48) where s denotes one of 2 N possible states s i is created from s by flipping i-th spin Probability density at equilibrium: hence: P eq ( s) = exp[ βh( s)]. (49) Z w i ( s) w i ( s i ) = P eq( s i ) P eq ( s) = exp[ βs i j <i> J ijs j ] exp[βs i j <i> J ijs j ] (50) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 53 / 63
54 Spin dynamics Glauber dynamics - transition rates (TR) w i ( s) w i ( s i ) = P eq( s i ) P eq ( s) = exp[ βs i j <i> J ijs j ] exp[βs i j <i> J ijs j ] Use the relation for s ± 1 (show it!): =? (51) e As = cosh(a) + s sinh(a) = cosh(a)[1 + s tanh(a)]. (52) then: w i ( s) 1 s i tanh (β ) w i ( s i ) = j <i> J ijs j 1 + s i tanh (β ) (53) j <i> J ijs j Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 54 / 63
55 Spin dynamics Glauber dynamics - transition rates (TR) The simplest choice w i ( s) 1 s i tanh (β ) w i ( s i ) = j <i> J ijs j 1 + s i tanh (β ) (54) j <i> J ijs j w i ( s) = W 0 [1 s i tanh ( β j <i> J ij s j )], (55) where W 0 is a constant chosen by Glauber as 1/2. What is a physical meaning of W 0? Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 55 / 63
56 Spin dynamics Ising model - MFA Let s write the Hamiltonia H = <i,j> JS i S j H = i h i S i, h i = 1 2 Js j. (56) j <i> Local field (from other spins) acts on every spin MFA: h i = h = Jzm/2. As we have seen it can be obtained from < δs i δs j >= 0. (57) Because s i =< s i > +δs i we obtain: < s i s j >=< s i >< s j >, (58) from mathematical point of view s i and s j are independent. Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 56 / 63
57 Spin dynamics Ising model on a complete graph Every spin interact with any other: H = 1 S i S j. (59) N i<j then: w i ( s) = W 0 [1 s i tanh = W 0 [1 s i tanh ( β ( β N j <i> J ij s j )] )] s j j i = W 0 [1 s i tanh (βm)]. (60) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 57 / 63
58 Spin dynamics Glauber dynamics - transition rates (TR) We obtain: therefore: s i (t + t) = w i ( s) = W 0 [1 s i tanh (βm)] (61) { si (t) with probability 1 w i ( s) t s i (t) with probability w i ( s) t Which means that spin at site i changes by 2s i with rate w i ( s) w i, therefore on average S i < s i >: (62) ds i dt = 2 < s iw i >= S i + tanh(βm). (63) Sum over all i: dm = m + tanh(βm). (64) dt Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 58 / 63
59 Spin dynamics Critical point and scaling At equilibrium: dm dt = 0 m = tanh(βm). (65) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 59 / 63
60 Spin dynamics Critical point and scaling At equilibrium: dm dt = 0 m = tanh(βm). (66) For β β c = 1 the only solution is m = 0. For β > β c = 1 there are two nontrivial solutions m ± Expand the equation above a critical point in βm: dm dt = tanh(βm) = (β c β)m 1 3 (βm) (67) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 60 / 63
61 Dynamics Spin dynamics dm dt = tanh(βm) = (β c β)m 1 3 (βm) (68) for T we have dm dt = (β c β)m m e t/τ, τ = (β c β) 1 (69) for T = T c : dm dt = 1 3 (βm) m t 1/2 (70) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 61 / 63
62 Spin dynamics Glauber dynamics in 1D We know that: in 1D simply: w i ( s) = W 0 [1 s i tanh w i ( s) = W 0 [1 s i tanh But foe ɛ = 0, ±1 we have (show it!): Hence: ( β j <i> J ij s j )], (71) ( βj s )] i 1 + s i+1. (72) 2 tanh(ɛx) = ɛ tanh(x). (73) [ ] s i 1 + s i+1 w i ( s) = W 0 1 γs i, γ = tanh(βj). (74) 2 Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 62 / 63
63 Spin dynamics Relaxation in 1D Glauber-Ising model As previously: ds i dt = 2 < s iw i >= S i + γ 2 (S i 1 + S i+1 ). (75) This equation can be solved assuming certain initial conditions (give an example) m(t) = m(0)e (1 γ)t γ = tanh(βj) (76) Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium and... SF-MTPT 63 / 63
Physics 115/242 Monte Carlo simulations in Statistical Physics
Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly
More information3.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 informationKinetic 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 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 informationThe Ising model Summary of L12
The Ising model Summary of L2 Aim: Study connections between macroscopic phenomena and the underlying microscopic world for a ferromagnet. How: Study the simplest possible model of a ferromagnet containing
More informationCoarsening process in the 2d voter model
Alessandro Tartaglia (LPTHE) Coarsening in the 2d voter model May 8, 2015 1 / 34 Coarsening process in the 2d voter model Alessandro Tartaglia LPTHE, Université Pierre et Marie Curie alessandro.tartaglia91@gmail.com
More informationJ ij S i S j B i S i (1)
LECTURE 18 The Ising Model (References: Kerson Huang, Statistical Mechanics, Wiley and Sons (1963) and Colin Thompson, Mathematical Statistical Mechanics, Princeton Univ. Press (1972)). One of the simplest
More informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More informationAny live cell with less than 2 live neighbours dies. Any live cell with 2 or 3 live neighbours lives on to the next step.
2. Cellular automata, and the SIRS model In this Section we consider an important set of models used in computer simulations, which are called cellular automata (these are very similar to the so-called
More informationarxiv:cond-mat/ v1 [cond-mat.other] 4 Aug 2004
Conservation laws for the voter model in complex networks arxiv:cond-mat/0408101v1 [cond-mat.other] 4 Aug 2004 Krzysztof Suchecki, 1,2 Víctor M. Eguíluz, 1 and Maxi San Miguel 1 1 Instituto Mediterráneo
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 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 informationarxiv: v1 [cond-mat.stat-mech] 18 Dec 2008
preprintaps/13-qed Updating schemes in zero-temperature single-spin flip dynamics Sylwia Krupa and Katarzyna Sznajd Weron Institute of Theoretical Physics, University of Wroc law, pl. Maxa Borna 9, 50-04
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 information7.1 Coupling from the Past
Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov
More 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 informationLecture 11: Non-Equilibrium Statistical Physics
Massachusetts Institute of Technology MITES 2018 Physics III Lecture 11: Non-Equilibrium Statistical Physics In the previous notes, we studied how systems in thermal equilibrium can change in time even
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationThermodynamics of phase transitions
Thermodynamics of phase transitions Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland 7 Oct 2013, SF-MTPT Katarzyna Sznajd-Weron (WUT) Thermodynamics of phase transitions
More informationLECTURE 10: Monte Carlo Methods II
1 LECTURE 10: Monte Carlo Methods II In this chapter, we discuss more advanced Monte Carlo techniques, starting with the topics of percolation and random walks, and then continuing to equilibrium statistical
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 informationDirectional Ordering in the Classical Compass Model in Two and Three Dimensions
Institut für Theoretische Physik Fakultät für Physik und Geowissenschaften Universität Leipzig Diplomarbeit Directional Ordering in the Classical Compass Model in Two and Three Dimensions vorgelegt von
More information8.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 informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationThermodynamics of phase transitions
Thermodynamics of phase transitions Katarzyna Sznajd-Weron Department of Theoretical of Physics Wroc law University of Science and Technology, Poland March 12, 2017 Katarzyna Sznajd-Weron (WUST) Thermodynamics
More informationLecture 3: September 10
CS294 Markov Chain Monte Carlo: Foundations & Applications Fall 2009 Lecture 3: September 10 Lecturer: Prof. Alistair Sinclair Scribes: Andrew H. Chan, Piyush Srivastava Disclaimer: These notes have not
More informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
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 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 informationAssignment 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 informationSpontaneous Symmetry Breaking
Spontaneous Symmetry Breaking Second order phase transitions are generally associated with spontaneous symmetry breaking associated with an appropriate order parameter. Identifying symmetry of the order
More informationPhase Transitions. µ a (P c (T ), T ) µ b (P c (T ), T ), (3) µ a (P, T c (P )) µ b (P, T c (P )). (4)
Phase Transitions A homogeneous equilibrium state of matter is the most natural one, given the fact that the interparticle interactions are translationally invariant. Nevertheless there is no contradiction
More informationPhysics 127b: Statistical Mechanics. Second Order Phase Transitions. The Ising Ferromagnet
Physics 127b: Statistical Mechanics Second Order Phase ransitions he Ising Ferromagnet Consider a simple d-dimensional lattice of N classical spins that can point up or down, s i =±1. We suppose there
More informationKinetics in a non-glauberian Ising model: some exact results
Kinetics in a non-glauberian Ising model: some exact results Malte Henkel Laboratoire de Physique de Matériaux Université Henri Poincaré Nancy I France collaborators: S.B. Dutta, H. Park (KIAS Seoul, Corea)
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 informationBasic Concepts and Tools in Statistical Physics
Chapter 1 Basic Concepts and Tools in Statistical Physics 1.1 Introduction Statistical mechanics provides general methods to study properties of systems composed of a large number of particles. It establishes
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 2 Apr 1998
Nonequilibrium Phase Transition in the Kinetic Ising model: Dynamical symmetry breaking by randomly varying magnetic field arxiv:cond-mat/979v2 [cond-mat.stat-mech] 2 Apr 998 Muktish Acharyya Institute
More information3. General properties of phase transitions and the Landau theory
3. General properties of phase transitions and the Landau theory In this Section we review the general properties and the terminology used to characterise phase transitions, which you will have already
More informationMarkov Chain Monte Carlo Method
Markov Chain Monte Carlo Method Macoto Kikuchi Cybermedia Center, Osaka University 6th July 2017 Thermal Simulations 1 Why temperature 2 Statistical mechanics in a nutshell 3 Temperature in computers 4
More 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 informationMicroscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/ v1 [cond-mat.stat-mech] 22 Sep 1999
Microscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/9909323v1 [cond-mat.stat-mech] 22 Sep 1999 B. Zheng FB Physik, Universität Halle, 06099 Halle, Germany Abstract Numerically we
More informationthe renormalization group (RG) idea
the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed
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 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 informationAdvanced Monte Carlo Methods Problems
Advanced Monte Carlo Methods Problems September-November, 2012 Contents 1 Integration with the Monte Carlo method 2 1.1 Non-uniform random numbers.......................... 2 1.2 Gaussian RNG..................................
More information3. The Voter Model. David Aldous. June 20, 2012
3. The Voter Model David Aldous June 20, 2012 We now move on to the voter model, which (compared to the averaging model) has a more substantial literature in the finite setting, so what s written here
More informationarxiv: v1 [physics.soc-ph] 9 Jan 2010
arxiv:1001.1441v1 [physics.soc-ph] 9 Jan 2010 Spontaneous reorientations in a model of opinion dynamics with anticonformists Grzegorz Kondrat and Katarzyna Sznajd-Weron Institute of Theoretical Physics,
More informationClusters and Percolation
Chapter 6 Clusters and Percolation c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 5 November 2012 6.1 Introduction In this chapter we continue our investigation of nucleation near the spinodal. We
More informationMarkov Chains and MCMC
Markov Chains and MCMC Markov chains Let S = {1, 2,..., N} be a finite set consisting of N states. A Markov chain Y 0, Y 1, Y 2,... is a sequence of random variables, with Y t S for all points in time
More informationFluctuation relations for equilibrium states with broken discrete symmetries
Journal of Statistical Mechanics: Theory and Experiment (0) P080 (9 pages) Fluctuation relations for equilibrium states with broken discrete symmetries Pierre Gaspard Center for Nonlinear Phenomena and
More informationSTOCHASTIC PROCESSES Basic notions
J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving
More informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
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 informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
More informationMathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge
Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking
More informationIntroduction 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 informationExam 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 informationThe Dynamics of Consensus and Clash
The Dynamics of Consensus and Clash Annual CNLS Conference 2006 Question: How do generic opinion dynamics models with (quasi)-ferromagnetic interactions evolve? Models: Voter model on heterogeneous graphs
More informationStochastic Processes
Stochastic Processes 8.445 MIT, fall 20 Mid Term Exam Solutions October 27, 20 Your Name: Alberto De Sole Exercise Max Grade Grade 5 5 2 5 5 3 5 5 4 5 5 5 5 5 6 5 5 Total 30 30 Problem :. True / False
More informationarxiv:cond-mat/ v1 10 Aug 2002
Model-free derivations of the Tsallis factor: constant heat capacity derivation arxiv:cond-mat/0208205 v1 10 Aug 2002 Abstract Wada Tatsuaki Department of Electrical and Electronic Engineering, Ibaraki
More informationDiscrete Markov Chain. Theory and use
Discrete Markov Chain. Theory and use Andres Vallone PhD Student andres.vallone@predoc.uam.es 2016 Contents 1 Introduction 2 Concept and definition Examples Transitions Matrix Chains Classification 3 Empirical
More informationChapter 2 Ensemble Theory in Statistical Physics: Free Energy Potential
Chapter Ensemble Theory in Statistical Physics: Free Energy Potential Abstract In this chapter, we discuss the basic formalism of statistical physics Also, we consider in detail the concept of the free
More informationQuasi-Stationary Simulation: the Subcritical Contact Process
Brazilian Journal of Physics, vol. 36, no. 3A, September, 6 685 Quasi-Stationary Simulation: the Subcritical Contact Process Marcelo Martins de Oliveira and Ronald Dickman Departamento de Física, ICEx,
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationMarkov Processes. Stochastic process. Markov process
Markov Processes Stochastic process movement through a series of well-defined states in a way that involves some element of randomness for our purposes, states are microstates in the governing ensemble
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 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 informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationSPIN DYNAMICS. Chapter The Voter Model
Chapter 7 SPIN DYNAMICS Kinetic spin systems play a crucial role in the development of non-equilibrium statistical physics. The prototypical example is the appealing-simple kinetic Ising model, in which
More informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
More informationVI. 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 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 informationCritical 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 informationChapter 4: Going from microcanonical to canonical ensemble, from energy to temperature.
Chapter 4: Going from microcanonical to canonical ensemble, from energy to temperature. All calculations in statistical mechanics can be done in the microcanonical ensemble, where all copies of the system
More informationIntroduction to Phase Transitions in Statistical Physics and Field Theory
Introduction to Phase Transitions in Statistical Physics and Field Theory Motivation Basic Concepts and Facts about Phase Transitions: Phase Transitions in Fluids and Magnets Thermodynamics and Statistical
More informationRensselaer Polytechnic Institute, 110 8th Street, Troy, NY , USA. Florida State University, Tallahassee, Florida , USA
5 Dynamic Phase Diagram for a Periodically Driven Kinetic Square-lattice Ising Ferromagnet: Finite-size Scaling Evidence for the Absence of a Tri-critical Point G. Korniss 1, P.A. Rikvold 2,3, and M.A.
More informationClassical and quantum simulation of dissipative quantum many-body systems
0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 Classical and quantum simulation of dissipative quantum many-body systems
More informationDecay of correlations in 2d quantum systems
Decay of correlations in 2d quantum systems Costanza Benassi University of Warwick Quantissima in the Serenissima II, 25th August 2017 Costanza Benassi (University of Warwick) Decay of correlations in
More informationWang-Landau sampling for Quantum Monte Carlo. Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart
Wang-Landau sampling for Quantum Monte Carlo Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart Overview Classical Monte Carlo First order phase transitions Classical Wang-Landau
More informationA 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 informationPhase transitions in Nowak-Sznajd opinion dynamics
arxiv:physics/69178v1 [physics.soc-ph] 21 Sep 26 Phase transitions in Nowak-Sznajd opinion dynamics Maciej Wo loszyn, Dietrich Stauffer* and Krzysztof Ku lakowski Faculty of Physics and Applied Computer
More informationAppendix 4. Appendix 4A Heat Capacity of Ideal Gases
Appendix 4 W-143 Appendix 4A Heat Capacity of Ideal Gases We can determine the heat capacity from the energy content of materials as a function of temperature. The simplest material to model is an ideal
More informationThermal and Statistical Physics Department Exam Last updated November 4, L π
Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =
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 informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationarxiv: v3 [cond-mat.stat-mech] 20 Nov 2014
Effective ergodicity in single-spin-flip dynamics Mehmet Süzen Applied Mathematical Physiology Lab Bonn University, Sigmund-Freud-Str.25, 53127 Bonn, Germany (Dated: March 9, 2018) arxiv:1405.4497v3 [cond-mat.stat-mech]
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 information2D 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 informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More informationMetastability for the Ising model on a random graph
Metastability for the Ising model on a random graph Joint project with Sander Dommers, Frank den Hollander, Francesca Nardi Outline A little about metastability Some key questions A little about the, and
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 informationA 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 informationCH 240 Chemical Engineering Thermodynamics Spring 2007
CH 240 Chemical Engineering Thermodynamics Spring 2007 Instructor: Nitash P. Balsara, nbalsara@berkeley.edu Graduate Assistant: Paul Albertus, albertus@berkeley.edu Course Description Covers classical
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationPhase transitions for particles in R 3
Phase transitions for particles in R 3 Sabine Jansen LMU Munich Konstanz, 29 May 208 Overview. Introduction to statistical mechanics: Partition functions and statistical ensembles 2. Phase transitions:
More informationV. SCALING. A. Length scale - correlation length. h f 1 L a. (50)
V. SCALING A. Length scale - correlation length This broken symmetry business gives rise to a concept pretty much missed until late into the 20th century: length dependence. How big does the system need
More informationMetastability and finite-size effects in magnetization switching
Metastability and finite-size effects in magnetization switching Per Arne Rikvold Department of Physics and MARTECH Florida State University, Tallahassee, FL Supported by NSF and FSU http://www.physics.fsu.edu/users/rikvold/info/rikvold.htm
More informationComputational Statistical Physics
Computational Statistical Physics Lecture Notes After a lecture by H.J. Herrmann written by Nicola Offeddu Marcel Thielmann Madis Ollikainen Summer Semester 2014 Preface This is the first part of a provvisory
More informationThe general reason for approach to thermal equilibrium of macroscopic quantu
The general reason for approach to thermal equilibrium of macroscopic quantum systems 10 May 2011 Joint work with S. Goldstein, J. L. Lebowitz, C. Mastrodonato, and N. Zanghì. Claim macroscopic quantum
More information