Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the...

26
Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising model when an initially disordered state is cooled to very low temperature? Sid Redner, Physics Department, Boston University, physics.bu.edu/~redner collaborators: P. L. Krapivsky, V. Spirin, K. Barros, J. Olejarz + NSF support

Transcript of Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the...

Page 1: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Fate of the Kinetic Ising Model

Basic question:

How does magnetic order emerge in a kinetic Ising model when an initially disordered state is cooled to very low temperature?

Sid Redner, Physics Department, Boston University, physics.bu.edu/~rednercollaborators: P. L. Krapivsky, V. Spirin, K. Barros, J. Olejarz + NSF support

Page 2: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

The System

H = !!

!i,j"

!i!j !i = ±1

Ising Hamiltonian

Initial state: disordered, zero magnetizationFinal state: ???

Page 3: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Dynamics of the System

Pick a random spin and compare the outcome after reversing the spin

if ∆E < 0 flip spin

if ∆E > 0 don’t flip

if ∆E = 0 flip with prob. 1/2

Page 4: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

time

+ + + + + + +_ _ _ _ __ space------+++++++---------------++++++++++----------------------++++++++++-----------------+++++++++----------+++++-------------------++++++++-----------+++++

Domain Wall Picture in 1d

Page 5: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

ξ ∼ t1/2

------+++++++---------------++++++++++----------------------++++++++++-----------------+++++++++----------+++++-------------------++++++++-----------+++++

time

+ + + + + + +_ _ _ _ __ space

• no. interfaces: ∝ t−1/2 (Glauber, 1963)• time to ground state: T ∝ L2

Domain Wall Picture in 1d

Page 6: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

t=4

Two Dimensions

coarsening of 256x256 system

Page 7: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

t=4 t=16

Two Dimensions

coarsening of 256x256 system

Page 8: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

t=4 t=16 t=64

Two Dimensions

coarsening of 256x256 system

Page 9: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

t=4 t=16 t=64 t=256

coarsening of 256x256 system

Two Dimensions

Page 10: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Evolution to the Ground State

Page 11: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Evolution to Stripe State

Page 12: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Evolution to Diagonal Stripe State

Page 13: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Two DimensionsQuestion: what is the final state?

0.0 0.1 0.2 0.3 0.41/log2L

0.0

0.1

0.2

0.3

0.4

strip

e pr

obab

ility

Answer from simulations:ground state with probability ≈ 2/3stripe state with probability ≈ 1/3

Page 14: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Basic result: ground state is never reached!typical 20x20x20 system

Features:1. Swiss cheesy2. Zero average curvature

Three Dimensions

Page 15: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Basic result: ground state is never reached!typical 20x20x20 system

Features:1. Swiss cheesy2. Zero average curvature

Three Dimensions

Page 16: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Evolution from Antiferromagnetic State

Page 17: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Evolution from Antiferromagnetic State

Page 18: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

blinker spin

Basic result: ground state is never reached!typical 20x20x20 system

Features:1. Swiss cheesy2. Zero average curvature3. Non-static

Three Dimensions (Olejarz, Krapvisky, & SR)

Page 19: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Blinker Evolution in Three Dimensions

Page 20: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Blinker Evolution in Three Dimensions

Page 21: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Genus Distribution

0 1 2 3 4g/〈g〉

0.0

0.2

0.4

0.6

0.8

1.0

P(g/〈g〉)

L=54

L=76

L=90

g=3

g=16

�g� ∼ L1.7

Page 22: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

Slow Relaxation of Blinkers

deflated inflatedintermediate

synthetic blinker configuration

Page 23: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising
Page 24: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising
Page 25: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

101

ln t10−2

10−1

100

S(t)

100 102 104 106 108t

10−6

10−4

10−2

100

S(t)

Slow Relaxation in 3d

46

8

10

14

20

3040

40

30

20

S(t) ∼ (ln t)−3 ?

slope -3

Page 26: Fate of the Kinetic Ising Modelphysics.bu.edu/~redner/482/13/redner-slides.pdf · Fate of the Kinetic Ising Model Basic question: How does magnetic order emerge in a kinetic Ising

final state: the ground statecompletion time: L²domain length distribution still unsolved

d=1: almost, but not quite, completely soluble

d≥3: topologically complex final statetopological connection between energy & genusperpetually blinking spinsultra-slow relaxation whose functional form is unknownfinite temperaturecorner geometry: partial result

rich state space structure

d=2: final state: usually the ground stateconnection to percolation crossing probabilitiescompletion time: usually L², sometimes L finite temperaturecorner geometry: exactly soluble

ground state/metastable stripe states

3.5

Summary & Open Problems