A Monte Carlo method to sample the classical configurational canonical ensemble is introduced. In contrast to the Metropolis algorithm, where trial moves can be rejected, in this approach collisions take place. The implementation is event-driven; i.e., at scheduled times the collisions occur. A unique feature of the new method is that smooth potentials (instead of only step-wise changing ones) can be used. In addition to an event-driven approach, where all particles move simultaneously, we introduce a straight event-chain implementation. As proof of principle, a system of Lennard-Jones particles is simulated.
|Number of pages
|Physical Review E - Statistical, Nonlinear, and Soft Matter Physics
|Published - 2012