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 -