TY - JOUR

T1 - Reduced density profile of small particles near a large particle

T2 - Results of an integral equation theory with an accurate bridge function and a Monte Carlo simulation

AU - Nakamura, Yuka

AU - Arai, Shota

AU - Kinoshita, Masahiro

AU - Yoshimori, Akira

AU - Akiyama, Ryo

N1 - Funding Information:
We thank Dr. A. Suematsu for the software program used in MC simulations to obtain the reduced density profiles. This work was supported by JSPS KAKENHI Grant Nos. JP19H01863, JP18H03673, JP18K03555, JP17H03663, JP16K05512, and JP15K05249.
Publisher Copyright:
© 2019 Author(s).

PY - 2019/7/28

Y1 - 2019/7/28

N2 - Solute-solvent reduced density profiles of hard-sphere fluids were calculated by using several integral equation theories for liquids. The traditional closures, Percus-Yevick (PY) and the hypernetted-chain (HNC) closures, as well as the theories with bridge functions, Verlet, Duh-Henderson, and Kinoshita (named MHNC), were used for the calculation. In this paper, a one-solute hard-sphere was immersed in a one-component hard-sphere solvent and various size ratios were examined. The profiles between the solute and solvent particles were compared with those calculated by Monte Carlo simulations. The profiles given by the integral equations with the bridge functions were much more accurate than those calculated by conventional integral equation theories, such as the Ornstein-Zernike (OZ) equation with the PY closure. The accuracy of the MHNC-OZ theory was maintained even when the particle size ratio of solute to solvent was 50. For example, the contact values were 5.7 (Monte Carlo), 5.6 (MHNC), 7.8 (HNC), and 4.5 (PY), and the first minimum values were 0.48 (Monte Carlo), 0.46 (MHNC), 0.54 (HNC), and 0.40 (PY) when the packing fraction of the hard-sphere solvent was 0.38 and the size ratio was 50. The asymptotic decay and the oscillation period for MHNC-OZ were also very accurate, although those given by the HNC-OZ theory were somewhat faster than those obtained by Monte Carlo simulations.

AB - Solute-solvent reduced density profiles of hard-sphere fluids were calculated by using several integral equation theories for liquids. The traditional closures, Percus-Yevick (PY) and the hypernetted-chain (HNC) closures, as well as the theories with bridge functions, Verlet, Duh-Henderson, and Kinoshita (named MHNC), were used for the calculation. In this paper, a one-solute hard-sphere was immersed in a one-component hard-sphere solvent and various size ratios were examined. The profiles between the solute and solvent particles were compared with those calculated by Monte Carlo simulations. The profiles given by the integral equations with the bridge functions were much more accurate than those calculated by conventional integral equation theories, such as the Ornstein-Zernike (OZ) equation with the PY closure. The accuracy of the MHNC-OZ theory was maintained even when the particle size ratio of solute to solvent was 50. For example, the contact values were 5.7 (Monte Carlo), 5.6 (MHNC), 7.8 (HNC), and 4.5 (PY), and the first minimum values were 0.48 (Monte Carlo), 0.46 (MHNC), 0.54 (HNC), and 0.40 (PY) when the packing fraction of the hard-sphere solvent was 0.38 and the size ratio was 50. The asymptotic decay and the oscillation period for MHNC-OZ were also very accurate, although those given by the HNC-OZ theory were somewhat faster than those obtained by Monte Carlo simulations.

UR - http://www.scopus.com/inward/record.url?scp=85070102868&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85070102868&partnerID=8YFLogxK

U2 - 10.1063/1.5100040

DO - 10.1063/1.5100040

M3 - Article

C2 - 31370562

AN - SCOPUS:85070102868

VL - 151

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 4

M1 - 044506

ER -