Title: WORM ALGORITHM FOR CLASSICAL AND QUANTUM
1WORM ALGORITHM FOR CLASSICAL AND QUANTUM
STATISTICAL MODELS
Nikolay Prokofiev, Umass, Amherst
Boris Svistunov, Umass, Amherst
Many thanks to collaborators on major algorithm
developments
Igor Tupitsyn, PITP
Vladimir Kashurnikov, MEPI, Moscow
Evgeni Burovski, Umass, Amherst
Massimo Boninsegni, UAlberta, Edmonton
NASA
Les Houches, June 2006
2Worm algorithm idea
Consider
- configuration space arbitrary closed loops
- each cnf. has a weight factor
- quantity of interest
3conventional sampling scheme
local shape change
Add/delete small loops
No sampling of topological classes
can not evolve to
dynamical critical exponent in many
cases
Critical slowing down
4Worm algorithm idea
draw and erase
Masha
Ira
Ira
Masha
Masha
Masha
keep drawing
or
Topological classes are (whatever you can draw!)
No critical slowing down in most cases
Disconnected loops relate to important physics
(correlation functions) and are not merely an
algorithm trick!
5High-T expansion for the Ising model
where
3
4
2
1
4
4
2
number of lines enter/exit rule
6Spin-spin correlation function
Worm algorithm cnf. space
Same as for generalized partition
1
I
4
3
M
4
2
7Getting more practical since
Complete algorithm
- If , select a new site for
at random
- select direction to move , let it be bond
- If accept
with prob.
8IM
I
M
M
M
M
Correlation function
Magnetization fluctuations
Energy either
or
9Ising lattice field theory
expand
if
closed oriented loops
where
tabulated numbers
10Flux in Flux out closed oriented loops
of integer
N-currents
I
(one open loop)
M
Z-configurations have
11Same algorithm
IM
sectors, prob. to accept
M
M
draw
M
M
erase
Keep drawing/erasing
12Multi-component gauge field-theory (deconfined
criticality, XY-VBS and Neel-VBS quantum phase
transitions
XY-VBS transition understood (?) no DCP, always
first-order
Neel-VBS transition, unknown !
13Winding numbers
Homogeneous gauge in x-direction
Ceperley Pollock 86