transitionProbability1D_R: Compute transition probability using Boltzmann distribution.

In isingLenzMC: Monte Carlo for Classical Ising Model

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 is pure R implementation.

Usage

  transitionProbability1D_R(ikBT, x, xFlip, J, H)

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

Value

Returns transition probability.

Author(s)

Mehmet Suzen <[e-mail]>

Examples

  n             <- 10 # 10 spin sites
  mySites       <- genConfig1D_R(n) # Generate sites
  mySitesNew    <- flipConfig1D_R(mySites)
  # only short-range part
  transitionProbability1D_R(1.0, mySites, mySitesNew, 1.0, 0.0)

Example output

[1] 0.9820138

isingLenzMC documentation built on May 2, 2019, 2:31 p.m.