In this paper we present a new approach to simulation methods for classical statistical mechanics relying on a field-theoretical formalism. It is based on applying the complex Hubbard-Stratonovich transformation to the canonical and grand-canonical partition function, which allows one to reexpress their particle representation interims of a functional integral over a fluctuating auxiliary field. The thermodynamic averages from the resulting field representations can then be calculated with a conventional Monte Carlo algorithm. We explored the applicability of the auxiliary field methodology for both the canonical and grand-canonical ensemble using a system of particles interacting through a purely repulsive Gaussian pair potential in a broad range of external parameters. In the grand-canonical case this technique represents an alternative to standard grand-canonical Monte Carlo methods. Generally providing a framework for simulating classical particle systems within a continuum formalism can be useful for multistage modeling where the field or continuum description naturally appears within quantum mechanics on smaller length scales and within classical mechanics on larger ones.
展开▼