A reference-modified density functional theory: An application to solvation free-energy calculations for a Lennard-Jones solution.

In the conventional classical density functional theory (DFT) for simple fluids, an ideal gas is usually chosen as the reference system because there is a one-to-one correspondence between the external field and the density distribution function, and the exact intrinsic free-energy functional is available for the ideal gas. In this case, the second-order density functional Taylor series expansion of the excess intrinsic free-energy functional provides the hypernetted-chain (HNC) approximation. Recently, it has been shown that the HNC approximation significantly overestimates the solvation free energy (SFE) for an infinitely dilute Lennard-Jones (LJ) solution, especially when the solute particles are several times larger than the solvent particles [T. Miyata and J. Thapa, Chem. Phys. Lett. 604, 122 (2014)]. In the present study, we propose a reference-modified density functional theory as a systematic approach to improve the SFE functional as well as the pair distribution functions. The second-order density functional Taylor series expansion for the excess part of the intrinsic free-energy functional in which a hard-sphere fluid is introduced as the reference system instead of an ideal gas is applied to the LJ pure and infinitely dilute solution systems and is proved to remarkably improve the drawbacks of the HNC approximation. Furthermore, the third-order density functional expansion approximation in which a factorization approximation is applied to the triplet direct correlation function is examined for the LJ systems. We also show that the third-order contribution can yield further refinements for both the pair distribution function and the excess chemical potential for the pure LJ liquids.

[1]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energy representation. III. Treatment of the molecular flexibility , 2003 .

[2]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energy representation. II. Functional for the chemical potential , 2002 .

[3]  David Chandler,et al.  Free energy functions in the extended RISM approximation , 1985 .

[4]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[5]  Jianzhong Wu,et al.  Structures of hard-sphere fluids from a modified fundamental-measure theory , 2002 .

[6]  Seiichiro Ten-no,et al.  Free energy of solvation for the reference interaction site model: Critical comparison of expressions , 2001 .

[7]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[8]  Daniel Borgis,et al.  Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections. , 2014, The journal of physical chemistry letters.

[9]  B. Montgomery Pettitt,et al.  Application of an extended RISM equation to dipolar and quadrupolar fluids , 1982 .

[10]  Tatsuhiko Miyata,et al.  Accuracy of solvation free energy calculated by hypernetted chain and Kovalenko–Hirata approximations for two-component system of Lennard-Jones liquid , 2014 .

[11]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[12]  J. A. Barker,et al.  Fourth Virial Coefficients for the 12–6 Potential , 1962 .

[13]  S. Ichimaru.  Strongly coupled plasmas: high-density classical plasmas and degenerate electron liquids , 1982 .

[14]  W. Curtin.  Freezing in the density functional approach: Effect of third‐order contributions , 1988 .

[15]  D. Henderson.  Fundamentals of Inhomogeneous Fluids , 1992 .

[16]  Jason Crain,et al.  Improved estimates for hydration free energy obtained by the reference interaction site model , 2007 .

[17]  Arieh Ben-Naim,et al.  Solvation: from small to macro molecules , 1994 .

[18]  Yutaka Maruyama,et al.  A solvation‐free‐energy functional: A reference‐modified density functional formulation , 2015, J. Comput. Chem..

[19]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[20]  Hansen,et al.  Factorization of the triplet direct correlation function in dense fluids. , 1987, Physical review letters.

[21]  Lloyd L. Lee,et al.  Constructing a new closure theory based on the third-order Ornstein-Zernike equation and a study of the adsorption of simple fluids. , 2011, The Journal of chemical physics.

[22]  Seiichiro Ten-no,et al.  Comparative study on solvation free energy expressions in reference interaction site model integral equation theory. , 2005, The journal of physical chemistry. B.

[23]  H. Iyetomi,et al.  Density-functional theory of correlations in dense plasmas: Improvement on the hypernetted-chain scheme , 1983 .

[24]  Rosenfeld,et al.  Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing. , 1989, Physical review letters.

[25]  Kenichiro Koga,et al.  Solvation of hydrophobes in water and simple liquids. , 2011, Physical chemistry chemical physics : PCCP.

[26]  Lloyd L. Lee,et al.  Correlation functions of classical fluids. III. The method of partition function variation applied to the chemical potential: Cases of PY and HNC2 , 1974 .

[27]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energetic representation. I. Formulation , 2000 .

[28]  B. Montgomery Pettitt,et al.  Integral equation predictions of liquid state structure for waterlike intermolecular potentials , 1982 .

[29]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[30]  Jianzhong Wu,et al.  A Site Density Functional Theory for Water: Application to Solvation of Amino Acid Side Chains. , 2013, Journal of chemical theory and computation.

[31]  H. Iyetomi,et al.  Improvement on the hypernetted-chain equations for dense plasmas , 1982 .

[32]  Gerhard Kahl,et al.  Fundamental measure theory for hard-sphere mixtures revisited: the White Bear version , 2002 .

[33]  Fumio Hirata,et al.  Hydration free energy of hydrophobic solutes studied by a reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method , 2000 .

[34]  Paul Labute,et al.  A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies , 2014, Journal of chemical theory and computation.

[35]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[36]  Neil W. Ashcroft,et al.  Theory of simple classical fluids: Universality in the short-range structure , 1979 .

[37]  Maxim V Fedorov,et al.  An accurate prediction of hydration free energies by combination of molecular integral equations theory with structural descriptors. , 2010, The journal of physical chemistry. B.

[38]  J. Lebowitz,et al.  Statistical Thermodynamics of Nonuniform Fluids , 1963 .

[39]  Weis,et al.  Iterative predictor-corrector method for extraction of the pair interaction from structural data for dense classical liquids. , 1986, Physical review. A, General physics.

[40]  H. Totsuji,et al.  Dielectric Response Function of Electron Liquids. II , 1970 .

[41]  Hideo Sekino,et al.  A Self-Consistent Density-Functional Approach for Homogeneous and Inhomogeneous Classical Fluids(Condensed matter: structure and mechanical and thermal properties) , 2008 .

[42]  D. Lévesque,et al.  Free-energy model for the inhomogeneous hard-sphere fluid mixture: Triplet and higher-order direct correlation functions in dense fluids , 1990 .

[43]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[44]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[45]  Y. Rosenfeld,et al.  Free energy model for inhomogeneous fluid mixtures: Yukawa‐charged hard spheres, general interactions, and plasmas , 1993 .

[46]  Jerome Percus,et al.  Approximation Methods in Classical Statistical Mechanics , 1962 .

[47]  K. E. Starling,et al.  Equation of State for Nonattracting Rigid Spheres , 1969 .

[48]  J. Hansen,et al.  On the equilibrium structure of dense fluids , 1988 .

[49]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. I. General Theory and Variational Formulation of the Mean Spherical Model and Hard Sphere Percus‐Yevick Equations , 1972 .