Pair force distributions in simple fluids

Abstract

Analytic expressions are derived for the frequency distribution, P(f), of pair forces, f, and those of their α-Cartesian component, f(α), or P(f(α)), for some typical model simple fluids, expressed in terms of the radial distribution function and known constants. For strongly repulsive inverse power (IP), exponential and Yukawa purely repulsive potentials, P(f) diverges at the origin approximately as ∼f(-1), but with different limiting analytic forms. P(f(α)) is also shown to diverge as ∼f(-1) as f → 0 for the IP fluid. For the Lennard-Jones potential fluid, P(f) is finite for all f ≥ 0 but has two singularities for negative f, corresponding to the zero force limit (i.e., f → 0(-)) and the point of inflection in the potential. The corresponding component force distribution is singular as f(α) → 0 from both positive and negative force sides. The large force limit of P(f), which originates from the close neighbor interactions, is nearly exponential for the IP and LJ fluids, as is also found for granular materials. A more complete picture of force distributions in off-lattice particulate systems as a function of force law and state point (particularly the extent of "thermalization" of the particles) is provided.