Introduction


So far, in these lectures, we have concentrated on equilibrium statistical mechanics and related computational-physics approaches, notably the equilibrium Monte external image arrow-10x10.png method. These and other approaches allowed us to external image arrow-10x10.png partition functions, energies, superfluid densities, etc. Physical time played a minor role, as the observables were generally time-independent. Likewise, Monte Carlo time was
treated as of secondary interest, if not a nuisance: we strove only to make things happen as quickly as possible, that is, to have algorithms converge rapidly.

In this lecture, we reach beyond equilibrium statistical mechanics, and explore time-dependent phenomena such as the crystallization of hard spheres after a sudden increase in pressure or the external image arrow-10x10.png response of Ising external image arrow-10x10.png to an external field switched on at some initial time. The local Monte Carlo algorithm often provides an excellent framework for studying dynamical phenomena.

We must be sure to understand the difference in paradigm between the role of the Monte Carlo time in equilibrium methods and in external image arrow-10x10.png algorithms: In the former, it is often unphysical, in the latter, the Monte Carlo time is taken as the model for the physical time, often the time-scale for diffusion.

The metropolis dynamics for a spin embedded in its molecular field



NE_ising_spin_field.jpg
A single spin in a field
We consider a single external image arrow-10x10.png σ in an external field h, and use the Metropolis algorithm as a dynamic model:



This transition probability satisfies detailed balance and ensures that at large times, the two external image arrow-10x10.png configurations
appear with their Boltzmann weights. Note that, if at time t the external image arrow-10x10.png is opposite to the field, it will be aligned with it at time t+1. The sampling problem is non-trivial only if the external image arrow-10x10.png σ = + 1.

Slow Metropolis dynamics and Faster than a clock algorithms


Ising model and the (BKL) algorithm

For the single-external image arrow-10x10.png model, the "faster-than-the-clock" algorithm, besides having a nice name, is actually quite efficient. It is therefore tempting to apply it to a non-trivial case, such as the Ising model.

In the one-external image arrow-10x10.png model the relevant parameter was λ , the probability "to do nothing". Let us therefore consider a external image arrow-10x10.png external image arrow-10x10.png, σ, and the same configuration, with spin k flipped, σk. The Metropolis probability to flip spin k is



The first term, 1/N, gives the probability of selecting spin k, followed by the Metropolis probability of accepting a flip of that spin. In the Ising model and its variants at low external image arrow-10x10.png, most spin-flips are rejected. It can then be interesting to implement a faster-than-the-clock algorithm which first samples the time of the next spin-flip, and then the spin to be flipped, just as in the earlier deposition problem. Rejections are avoided altogether, although the method is not unproblematic.


The probability to do nothing is



This external image arrow-10x10.png expresses that to determine the probability of doing nothing, we must know all the probabilities for N flipping spins. Naively, we
can recalculate λ after each step, and sample the waiting time as. After finding out when to flip the next spin, we must decide on which of the N spins to flip. This problem is solved through a second application of tower sampling. However, each spin-flip requires of the order of N operations.

We can actually deal with this cumbersome constructions using an algorithm due to Bortz, Kalos, and Lebowitz [1] . The algorithm is based on classes and is not in the course external image arrow-10x10.png.



References


  1. ^ Bortz A. B., Kalos M. H., Lebowitz J. L. (1975) A new algorithm for Monte Carlo simulations of Ising spin systems Journal of Chemical Physics 17, 10