Numerical simulations
James Jeffrey Binney, N J Dowrick, Andrew J. Fisher, Mark E. J. Newman
Abstract
Abstract We have now seen how powerful the tools of statistical mechanics can be in the study of complex systems in thermal equilibrium and how they can be used to determine exactly the behaviour of some simplified, but nevertheless interesting, models of phase transitions. However, these exact solutions are sometimes lengthy and difficult and there is no guarantee that any particular model will yield exact results. For example, we were able to solve the two-dimensional, but not the three-dimensional, Ising model, and that only in zero magnetic field. So, it is natural to consider whether we can use a computer to simulate these systems. In this chapter we shall see that this is indeed possible; a great deal of our understanding of critical phenomena is based wholly or partly on the results of such numerical work. However, as we shall also see, such simulations are not entirely straightforward and a number of tricks are needed to obtain accurate results. This is partly a consequence of the very physics that makes the phase transitions interesting to start with. One should probably think of numerical simulation as a whole subject, almost an art, of its own; we shall have time only to scratch at the surface and discuss a few of the most important ideas.