Overcoming the minimum image constraint using the closest point search

Abstract

Finding the set of nearest images of a point in a simulation cell with periodic (torus) boundary conditions is of central importance for molecular dynamics algorithms. To compute all pairwise distances closer than a given cutoff in linear time requires region-based neighbor-listing algorithms. Available algorithms encounter increasing difficulties when the cutoff distance exceeds half the shortest cell length. This work provides details on two ways to directly and efficiently generate region–region interaction lists in n-dimensional space, free from the minimum image restriction. The solution is based on a refined version of existing algorithms solving the closest vector problem. A self-contained discussion of lattice reduction methods for efficient higher-dimensional searches is also provided. In the MD setting, these reduction criteria provide useful guidelines for lattice compaction.