Numerical Calculation of the Bridge Function for Soft-Sphere Supercooled Fluids via Molecular Dynamics Simulations

Abstract

Isokinetic molecular dynamics simulations have been performed for 13,500 soft-spheres interacting through the inverse-power potential, ε(sgrave]/r)ⁿ , near and below the freezing temperature. The bridge function for the integral equation of the theory of liquids is extracted from the pair distribution function (PDF) obtained by the computer simulations for n = 6 and 12. The result is compared with that of approximate theories, i.e., the Rogers-Young (RY) approximation and a modified hypernetted-chain approximation for supercooled soft-sphere fluids (MHNCS approximation). Below the freezing temperature, the bridge function obtained by the computer simulation begins to oscillate around zero at intermediate distances where the second peak of the PDF appears. Such oscillatory behavior of the bridge function is well reproduced by the MHNCS approximation which includes correlations given by the leading elementary diagram, in remarkable contrast to that of the RY approximation. The present result suggests that the split second peak of the PDF for highly supercooled liquids is essentially dominated by the intermediate-distance-range correlation of the leading elementary diagram.