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.

Classical density-functional theory is the most direct approach to equilibrium structures and free energies of inhomogeneous liquids, but requires the construction of an approximate free-energy functional for each liquid of interest. We present a general recipe for constructing functionals for small-molecular liquids based only on bulk experimental properties and ab initio calculations of a single solvent molecule.

A major challenge in the development of new battery materials is understanding their fundamental mechanisms of operation and degradation. Their microscopically inhomogeneous nature calls for characterization tools that provide operando and localized information from individual grains and particles. Here, we describe an approach that enables imaging the nanoscale distribution of ions during electrochemical charging of a battery in a transmission electron microscope liquid flow cell.

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 fluids 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.

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 Vienna ab initio Software Package. The implicit solvation model follows the framework of joint density functional theory.

Delivering the full benefits of first-principles calculations to battery materials demands the development of accurate and computationally efficient electronic structure methods that incorporate the effects of the electrolyte environment and electrode potential. Realistic electrochemical interfaces containing polar surfaces are beyond the regime of validity of existing continuum solvation theories developed for molecules, due to the presence of significantly stronger electric fields.

Hybrid density functionals show great promise for chemically accurate first-principles calculations, but their high computational cost limits their application in nontrivial 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.

This work explores the use of joint density functional theory, an extension of density functional theory for the ab initio description of electronic systems in thermodynamic equilibrium with a liquid environment, to describe electrochemical systems. After reviewing the physics of the underlying fundamental electrochemical concepts, we identify the mapping between commonly measured electrochemical observables and microscopically computable quantities within an, in principle, exact theoretical framework.

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 nonlinear dielectric response of water, making it an excellent candidate for the studies of mesoscale water and for use in ab initio solvation methods.

Employing a classical density-functional description of liquid environments, we introduce a rigorous method for the diffusion quantum Monte Carlo calculation of free energies and thermodynamic averages of solvated systems that requires neither thermodynamic sampling nor explicit solvent electrons. We find that this method yields promising results and small convergence errors for a set of test molecules. It is implemented readily and is applicable to a range of challenges in condensed matter, including the study of transition states of molecular and surface reactions in liquid environments.