TY - JOUR

T1 - Solvation effects on diffusion processes of a macromolecule

T2 - Accuracy required for radial distribution function to calculate diffusion coefficient

AU - Nakamura, Yuka

AU - Yoshimori, Akira

AU - Akiyama, Ryo

N1 - Funding Information:
We thank Professor Masahiro Kinoshita of Kyoto University for the software program used in the numerical preparation to obtain the radial distribution functions. This work was supported by the JSPS KAKENHI (Grant Nos. JP19K03772, JP19H01863, JP18H03673, JP18K03555, and JP16K05512). R.A. thanks the Research Center for Computational Science, Okazaki, Japan, and the Research Institute for Information Technology, Kyushu University.
Publisher Copyright:
© 2021 Author(s).

PY - 2021/2/28

Y1 - 2021/2/28

N2 - We investigate the dependence of the diffusion coefficient of a large solute particle on the solvation structure around a solute. The diffusion coefficient of a hard-sphere system is calculated by using a perturbation theory of large-particle diffusion with radial distribution functions around the solute. To obtain the radial distribution function, some integral equation theories are examined, such as the Percus-Yevick (PY), hypernetted-chain (HNC), and modified HNC theories using a bridge function proposed by Kinoshita (MHNC) closures. In one-component solvent systems, the diffusion coefficient depends on the first-minimum value of the radial distribution function. The results of the MHNC closure are in good agreement with those of calculation using the radial distribution functions of Monte Carlo simulations since the MHNC closure very closely reproduces the radial distribution function of Monte Carlo simulations. In binary-solvent mixtures, the diffusion coefficient is affected by the larger solvent density distribution in the short-range part, particularly the height and sharpness of the first peak and the depth of the first minimum. Since the HNC closure gives the first peak that is higher and sharper than that of the MHNC closure, the calculated diffusion coefficient is smaller than the MHNC closure result. In contrast, the results of the PY closure are qualitatively and quantitatively different from those of the MHNC and HNC closures.

AB - We investigate the dependence of the diffusion coefficient of a large solute particle on the solvation structure around a solute. The diffusion coefficient of a hard-sphere system is calculated by using a perturbation theory of large-particle diffusion with radial distribution functions around the solute. To obtain the radial distribution function, some integral equation theories are examined, such as the Percus-Yevick (PY), hypernetted-chain (HNC), and modified HNC theories using a bridge function proposed by Kinoshita (MHNC) closures. In one-component solvent systems, the diffusion coefficient depends on the first-minimum value of the radial distribution function. The results of the MHNC closure are in good agreement with those of calculation using the radial distribution functions of Monte Carlo simulations since the MHNC closure very closely reproduces the radial distribution function of Monte Carlo simulations. In binary-solvent mixtures, the diffusion coefficient is affected by the larger solvent density distribution in the short-range part, particularly the height and sharpness of the first peak and the depth of the first minimum. Since the HNC closure gives the first peak that is higher and sharper than that of the MHNC closure, the calculated diffusion coefficient is smaller than the MHNC closure result. In contrast, the results of the PY closure are qualitatively and quantitatively different from those of the MHNC and HNC closures.

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

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

U2 - 10.1063/5.0038894

DO - 10.1063/5.0038894

M3 - Article

C2 - 33639733

AN - SCOPUS:85101346877

VL - 154

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 8

M1 - 38894

ER -