Continuum solvation models enable electronic structure calculations of systems in liquid environments, but because of the large number of empirical parameters, they are limited to the class of systems in their fit set (typically organic molecules). H
ere, we derive a solvation model with no empirical parameters for the dielectric response by taking the linear response limit of a classical density functional for molecular liquids. This model directly incorporates the nonlocal dielectric response of the liquid using an angular momentum expansion, and with a single fit parameter for dispersion contributions it predicts solvation energies of neutral molecules with an RMS error of 1.3 kcal/mol in water and 0.8 kcal/mol in chloroform and carbon tetrachloride. We show that this model is more accurate for strongly polar and charged systems than previous solvation models because of the parameter-free electric response, and demonstrate its suitability for ab initio solvation, including self-consistent solvation in quantum Monte Carlo calculations.
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.
Solid-liquid interfaces are at the heart of many modern-day technologies and provide a challenge to many materials simulation methods. A realistic first-principles computational study of such systems entails the inclusion of solvent effects. In this
work we implement an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code VASP. The implicit solvation model follows the framework of joint density functional theory. We describe the framework, our algorithm and implementation, and benchmarks for small molecular systems. We apply the solvation model to study the surface energies of different facets of semiconducting and metallic nanocrystals and the S$_{text{N}} 2$ reaction pathway. We find that solvation reduces the surface energies of the nanocrystals, especially for the semiconducting ones and increases the energy barrier of the S$_{text{N}} 2$ reaction.
Hybrid density functionals show great promise for chemically-accurate first principles calculations, but their high computational cost limits their application in non-trivial studies, such as exploration of reaction pathways of adsorbents on periodic
surfaces. One factor responsible for their increased cost is the dense Brillouin-zone sampling necessary to accurately resolve an integrable singularity in the exact exchange energy. We analyze this singularity within an intuitive formalism based on Wannier-function localization and analytically prove Wigner-Seitz truncation to be the ideal method for regularizing the Coulomb potential in the exchange kernel. We show that this method is limited only by Brillouin-zone discretization errors in the Kohn-Sham orbitals, and hence converges the exchange energy exponentially with the number of k-points used to sample the Brillouin zone for all but zero-temperature metallic systems. To facilitate the implementation of this method, we develop a general construction for the plane-wave Coulomb kernel truncated on the Wigner-Seitz cell in one, two or three lattice directions. We compare several regularization methods for the exchange kernel in a variety of real systems including low-symmetry crystals and low-dimensional materials. We find that our Wigner-Seitz truncation systematically yields the best k-point convergence for the exchange energy of all these systems and delivers an accuracy to hybrid functionals comparable to semi-local and screened-exchange functionals at identical k-point sets.
Classical density-functional theory provides an efficient alternative to molecular dynamics simulations for understanding the equilibrium properties of inhomogeneous fluids. However, application of density-functional theory to multi-site molecular fl
uids has so far been limited by complications due to the implicit molecular geometry constraints on the site densities, whose resolution typically requires expensive Monte Carlo methods. Here, we present a general scheme of circumventing this so-called inversion problem: compressed representations of the orientation density. This approach allows us to combine the superior iterative convergence properties of multipole representations of the fluid configuration with the improved accuracy of site-density functionals. Next, from a computational perspective, we show how to extend the DFT++ algebraic formulation of electronic density-functional theory to the classical fluid case and present a basis-independent discretization of our formulation for molecular classical density-functional theory. Finally, armed with the above general framework, we construct a simplified free-energy functional for water which captures the radial distributions, cavitation energies, and the linear and non-linear dielectric response of liquid water. The resulting approach will enable efficient and reliable first-principles studies of atomic-scale processes in contact with solution or other liquid environments.
We present an accurate equation of state for water based on a simple microscopic Hamiltonian, with only four parameters that are well-constrained by bulk experimental data. With one additional parameter for the range of interaction, this model yields
a computationally efficient free-energy functional for inhomogeneous water which captures short-ranged correlations, cavitation energies and, with suitable long-range corrections, the non-linear dielectric response of water, making it an excellent candidate for studies of mesoscale water and for use in ab initio solvation methods.
