A rejection-free Monte Carlo method for the hard-disk system

TL;DR

This paper introduces a rejection-free Monte Carlo method tailored for the hard-disk system, enhancing computational efficiency at high densities while accurately preserving the system's dynamic behavior.

Contribution

The paper presents a novel rejection-free Monte Carlo algorithm for hard-disk systems, improving efficiency and maintaining accurate dynamic evolution compared to standard methods.

Findings

Enhanced efficiency at high densities

Preserves nonequilibrium relaxation behavior

Applicable to two-dimensional solid dynamics

Abstract

We construct a rejection-free Monte Carlo method for the hard-disk system. Rejection-free Monte Carlo methods preserve the time-evolution behavior of the standard Monte Carlo method, and this relationship is confirmed for our method by observing nonequilibrium relaxation of a bond-orientational order parameter. The rejection-free method gives a greater computational efficiency than the standard method at high densities. The rejection free method is implemented in a shrewd manner using optimization methods to calculate a rejection probability and to update the system. This method should allow an efficient study of the dynamics of two-dimensional solids at high density.