The Measurement of Free Energy by
Monte Carlo Computer Simulation
Graham R Smith
A thesis submitted in ful lment of the requirements
for the degree of Do ctor of Philosophy
to the
University of Edinburgh
Abstract
One of the most imp ortant problems in statistical mechanics is the measurement of free energies
these b eing the quantities that determine the direction of chemical reactions and the concern
of this thesis the lo cation of phase transitions While Monte Carlo MC computer simulation
is a well established and invaluable aid in statistical mechanical calculations it is well known
that in its most commonly practised form where samples are generated from the Boltzmann
distribution it fails if applied directly to the free energy problem This failure o ccurs b ecause
the measurement of free energies requires a much more extensive exploration of the system s
con guration space than do most statistical mechanical calculations con gurations which have
a very low Boltzmann probability make a substantial contribution to the free energy and the
imp ortant regions of con guration space may b e separated by p otential barriers
We b egin the thesis with an introduction and then give a review of the very substantial
literature that the problem of the MC measurement of free energy has pro duced explaining
and classifying the various di erent approaches that have b een adopted We then pro ceed to
present the results of our own investigations
First we investigate metho ds in which the con gurations of the system are sampled from a
distribution other than the Boltzmann distribution concentrating in particular on a recently
developed technique known as the multicanonical ensemble The principal di culty in using
the multicanonical ensemble is the di culty of constructing it implicit in it is at least partial
knowledge of the very free energy that we are trying to measure and so to pro duce it requires an
iterative pro cess Therefore we study this iterative pro cess using Bayesian inference to extend
the usual metho d of MC data analysis and introducing a new MC metho d in which inferences
are made based not on the macrostates visited by the simulation but on the transitions made
b etween them We present a detailed comparison b etween the multicanonical ensemble and
the traditional metho d of free energy measurement thermo dynamic integration and use the i
former to make a high accuracy investigation of the critical magnetisation distribution of the