The standard Metropolis algorithm [32] was applied with equal
probabilities of move rotation, insertion, and deletion. The cut-off
radius was set at half of the box length and 1.2 108 configurations
were run for both equilibration and sampling.
The average absolute and excess volumetric densities are defined
by: