Compute transition probability using Boltzmann distribution.

Share:

Description

Given a vector of flip sites, 1s or -1s, representing up and down spins respectively, and an other flip sites, return the transition probability, applying periodic boundary conditions, i.e., cyclic. This function calls the C function 'transitionProbability1D'.

Usage

1
   transitionProbability1D(ikBT, x, xflip, J, H, probSel)

Arguments

ikBT

1/kB*T (Boltzmann factor)

x

1D Spin sites on the lattice.

xflip

1D Spin sites on the lattice: after a flip.

J

Interaction strength

H

External field

probSel

Which transition probability to use. 1 for Metropolis 2 for Glauber

Value

Returns transition probability.

Author(s)

Mehmet Suzen <mehmet.suzen@physics.org>

Examples

1
2
3
4
5
6
  n             <- 10             # 10 spin sites
  mySites       <- genConfig1D(n) # Generate sites
  mySitesNew    <- flipConfig1D(mySites)
  # only short-range part
  transitionProbability1D(1.0, mySites, mySitesNew, 1.0, 0.0, 1) # Metropolis
  transitionProbability1D(1.0, mySites, mySitesNew, 1.0, 0.0, 2) # Glauber