How heteroepitaxy occurs on strontium titanate.

In traditional models of heteroepitaxy, the substrate serves mainly as a crystalline template for the thin-film lattice, dictating the initial roughness of the film and the degree of coherent strain. Here, performing in situ surface x-ray diffraction during the heteroepitaxial growth of LaTiO3 on SrTiO3 (001), we find that a TiO2 adlayer composed of the ( 13 × 13 ) R33.7° and ( 2 × 2 ) R45.0° reconstructions is a highly active participant in the growth process, continually diffusing to the surface throughout deposition. The effects of the TiO2 adlayer on layer-by-layer growth are investigated using different deposition sequences and anomalous x-ray scattering, both of which permit detailed insight into the dynamic layer rearrangements that take place. Our work challenges commonly held assumptions regarding growth on TiO2-terminated SrTiO3 (001) and demonstrates the critical role of excess TiO2 surface stoichiometry on the initial stages of heteroepitaxial growth on this important perovskite oxide substrate material.

Improving accuracy of electrochemical capacitance and solvation energetics in first-principles calculations.

Reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model.

JDFTx: software for joint density-functional theory.

Density-functional theory (DFT) has revolutionized computational prediction of atomic-scale properties from first principles in physics, chemistry and materials science. Continuing development of new methods is necessary for accurate predictions of new classes of materials and properties, and for connecting to nano- and mesoscale properties using coarse-grained theories. JDFTx is a fully-featured open-source electronic DFT software designed specifically to facilitate rapid development of new theories, models and algorithms. Using an algebraic formulation as an abstraction layer, compact C++11 code automatically performs well on diverse hardware including GPUs (Graphics Processing Units). This code hosts the development of joint density-functional theory (JDFT) that combines electronic DFT with classical DFT and continuum models of liquids for first-principles calculations of solvated and electrochemical systems. In addition, the modular nature of the code makes it easy to extend and interface with, facilitating the development of multi-scale toolkits that connect to ab initio calculations, e.g. photo-excited carrier dynamics combining electron and phonon calculations with electromagnetic simulations.

Structure of the Photo-catalytically Active Surface of SrTiO3.

A major goal of energy research is to use visible light to cleave water directly, without an applied voltage, into hydrogen and oxygen. Although SrTiO3 requires ultraviolet light, after four decades, it is still the "gold standard" for the photo-catalytic splitting of water. It is chemically robust and can carry out both hydrogen and oxygen evolution reactions without an applied bias. While ultrahigh vacuum surface science techniques have provided useful insights, we still know relatively little about the structure of these electrodes in contact with electrolytes under operating conditions. Here, we report the surface structure evolution of a n-SrTiO3 electrode during water splitting, before and after "training" with an applied positive bias. Operando high-energy X-ray reflectivity measurements demonstrate that training the electrode irreversibly reorders the surface. Scanning electrochemical microscopy at open circuit correlates this training with a 3-fold increase of the activity toward the photo-induced water splitting. A novel first-principles joint density functional theory simulation, constrained to the X-ray data via a generalized penalty function, identifies an anatase-like structure as the more active, trained surface.

Spicing up continuum solvation models with SaLSA: the spherically averaged liquid susceptibility ansatz.

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). Here, 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 a 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.

A recipe for free-energy functionals of polarizable molecular fluids.

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. This recipe combines the exact free energy of the non-interacting system with fundamental measure theory for the repulsive contribution and a weighted density functional for the short-ranged attractive interactions. We add to these ingredients a weighted polarization functional for the long-range correlations in both the rotational and molecular-polarizability contributions to the dielectric response. We also perform molecular dynamics calculations for the free energy of cavity formation and the high-field dielectric response, and show that our free-energy functional adequately describes these properties (which are key for accurate solvation calculations) for all three solvents in our study: water, chloroform, and carbon tetrachloride.

Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways.

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. 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 SN2 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 SN2 reaction.

A computationally efficacious free-energy functional for studies of inhomogeneous liquid water.

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.