Solute-solvent cavity and bridge functions. I. Varying size of the solute

In this work we present the results of the extensive molecular simulations of solute-solvent cavity and bridge functions. The mixtures of Lennard-Jones solvent with Lennard-Jones solute at infinite dilution are considered for different solute-solvent size ratios-up to 4:1. The Percus-Yevick and hypernetted chain closures deviate substantially from simulation results in the investigated temperature and density ranges. We also find that the behavior of the indirect and cavity correlation functions is non-monotonous within the hard-core region, but the latter can be successfully approximated by mean-field theory if the solute-solvent interaction energy is divided into repulsive and attractive contribution, according to Weeks-Chandler-Andersen theory. Furthermore, in spite of the non-monotonous behavior of logarithm of the cavity function and the indirect correlation function, their difference, i.e., the bridge function remains constant within the hard-core region. Such behavior of the bridge and indirect correlation functions at small distances and for small values of indirect correlation function is well known from the Duh-Haymet plots, where the non-unique relationship results in loops of the bridge function vs. indirect correlation function graphs. We show that the same pathological behavior appears also when distance is small and indirect correlation function is large. We further show that the unique functional behavior of the bridge function can be established when bridge is represented as a function of the renormalized, repulsive indirect correlation function.