Distance-based optimal sampling in a hypercube: Analogies to N-body systems

A method is proposed for the construction of uniformly distributed point sets within a design domain using an analogy to a dynamical system of interacting particles. The possibility of viewing various distance-based optimality criteria as formulas representing the potential energy of a system of charged particles is discussed. The potential energy is employed in deriving the equations of motion of the particles. The particles are either attracted or repelled and dissipative dynamical systems can be simulated to achieve optimal and near-optimal arrangements of points. The design domain is set up as an N-var dimensional unit hypercube, with N-var being the number of variables (factors). The number of points is equal to the number of simulations (levels). The periodicity assumption of the design domain is shown to be an elegant way to obtain statistically uniform coverage of the design domain. The phi(p) criterion, which is a generalization of the Maximin criterion, is selected in order to demonstrate its analogy with an N-body system. This criterion guarantees that the points are spread uniformly within the design domain. The solution to such an N-body system is presented. The obtained designs are shown to outperform the existing optimal designs in various types of applications: multidimensional numerical integration, statistical exploration of computer models, reliability analyses of engineering systems, and screenings or exploratory designs for the global optimization/minimization of functions.