Weighted-density functionals for cavity formation and dispersion energies in continuum solvation models

Abstract

Continuum solvation models enable efficient first principles calculations of chemical reactions in solution, but require extensive parametrization and fitting for each solvent and class of solute systems. Here, we examine the assumptions of continuum solvation models in detail and replace empirical terms with physical models in order to construct a minimally-empirical solvation model. Specifically, we derive solvent radii from the nonlocal dielectric response of the solvent from ab initio calculations, construct a closed-form and parameter-free weighted-density approximation for the free energy of the cavity formation, and employ a pair-potential approximation for the dispersion energy. We show that the resulting model with a single solvent-independent parameter: the electron density threshold ($n_c$), and a single solvent-dependent parameter: the dispersion scale factor ($s_6$), reproduces solvation energies of organic molecules in water, chloroform and carbon tetrachloride with RMS errors of 1.1, 0.6 and 0.5 kcal/mol respectively. We additionally show that fitting the solvent-dependent $s_6$ parameter to the solvation energy of a single non-polar molecule does not substantially increase these errors. Parametrization of this model for other solvents, therefore, requires minimal effort and is possible without extensive databases of experimental solvation free energies.

References in corpus (2)

• Joint density-functional theory for electronic structure of solvated systems
• Universal iso-density polarizable continuum model for molecular solvents

Cited by in corpus (5)

• Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry
• JDFTx: software for joint density-functional theory
• The charge-asymmetric nonlocally-determined local-electric (CANDLE) solvation model
• Spicing up continuum solvation models with SaLSA: the spherically-averaged liquid susceptibility ansatz
• Range separated hybrid density functional study of organic dye sensitizers on anatase TiO${}_2$ nanowires