期刊
JOURNAL OF CHEMICAL PHYSICS
卷 134, 期 19, 页码 -出版社
AMER INST PHYSICS
DOI: 10.1063/1.3589142
关键词
-
资金
- Agence Nationale pour la Recherche [ANR-06-BLAN-0147]
- Agence Nationale de la Recherche (ANR) [ANR-06-BLAN-0147] Funding Source: Agence Nationale de la Recherche (ANR)
A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3589142]
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据