Kathleen A. Schwarz

Nanoscale imaging of lithium ion distribution during in situ operation of battery electrode and electrolyte.

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. We use valence energy-loss spectroscopy to track both solvated and intercalated ions, with electronic structure fingerprints of the solvated ions identified using an ab initio nonlinear response theory. Equipped with the new electrochemical cell holder, nanoscale spectroscopy and theory, we have been able to determine the lithiation state of a LiFePO4 electrode and surrounding aqueous electrolyte in real time with nanoscale resolution during electrochemical charge and discharge. We follow lithium transfer between electrode and electrolyte and image charging dynamics in the cathode. We observe competing delithiation mechanisms such as core-shell and anisotropic growth occurring in parallel for different particles under the same conditions. This technique represents a general approach for the operando nanoscale imaging of electrochemically active ions in the electrode and electrolyte in a wide range of electrical energy storage systems.

The importance of nonlinear fluid response in joint density-functional theory studies of battery systems

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. We present an ab initio theory of the nonlinear dielectric and ionic response of solvent environments within the framework of joint density-functional theory, with precisely the same optimizable parameters as conventional polarizable continuum models. We demonstrate that the resulting nonlinear theory agrees with the standard linear models for organic molecules and metallic surfaces under typical operating conditions. However, we find that the saturation effects in the rotational response of polar solvent molecules, inherent to our nonlinear theory, are crucial for a qualitatively correct description of the ionic surfaces typical of the solid electrolyte interface.

Spectroscopic Characterization of Charged Defects in Polycrystalline Pentacene by Time‐ and Wavelength‐Resolved Electric Force Microscopy

Spatial maps of topography and trapped charge are acquired for polycrystalline pentacene thin-film transistors using electric and atomic force microscopy. In regions of trapped charge, the rate of trap clearing is studied as a function of the wavelength of incident radiation.

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.

Evaluating continuum solvation models for the electrode-electrolyte interface: Challenges and strategies for improvement

Ab initio modeling of electrochemical systems is becoming a key tool for understanding and predicting electrochemical behavior. Development and careful benchmarking of computational electrochemical methods are essential to ensure their accuracy. Here, using charging curves for an electrode in the presence of an inert aqueous electrolyte, we demonstrate that most continuum models, which are parameterized and benchmarked for molecules, anions, and cations in solution, undersolvate metal surfaces, and underestimate the surface charge as a function of applied potential. We examine features of the electrolyte and interface that are captured by these models and identify improvements necessary for realistic electrochemical calculations of metal surfaces. Finally, we reparameterize popular solvation models using the surface charge of Ag(100) as a function of voltage to find improved accuracy for metal surfaces without significant change in utility for molecular and ionic solvation.

Partial oxidation of step-bound water leads to anomalous pH effects on metal electrode step-edges.

The design of better heterogeneous catalysts for applications such as fuel cells and electrolyzers requires a mechanistic understanding of electrocatalytic reactions and the dependence of their activity on operating conditions such as pH. A satisfactory explanation for the unexpected pH dependence of electrochemical properties of platinum surfaces has so far remained elusive, with previous explanations resorting to complex co-adsorption of multiple species and resulting in limited predictive power. This knowledge gap suggests that the fundamental properties of these catalysts are not yet understood, limiting systematic improvement. Here, we analyze the change in charge and free energies upon adsorption using density-functional theory (DFT) to establish that water adsorbs on platinum step edges across a wide voltage range, including the double-layer region, with a loss of approximately 0.2 electrons upon adsorption. We show how this as-yet unreported change in net surface charge due to this water explains the anomalous pH variations of the hydrogen underpotential deposition (Hupd) and the potentials of zero total charge (PZTC) observed in published experimental data. This partial oxidation of water is not limited to platinum metal step edges, and we report the charge of the water on metal step edges of commonly used catalytic metals, including copper, silver, iridium, and palladium, illustrating that this partial oxidation of water broadly influences the reactivity of metal electrodes.

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.

Framework for solvation in quantum Monte Carlo

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.

Direct visualization of sulfur cathodes: new insights into Li–S batteries via operando X-ray based methods

As the need for the development of “beyond lithium” ion battery technologies continues unabated, lithium sulfur batteries have attracted widespread attention due to their very high theoretical energy density of 2600 W h kg−1. However, despite much effort, the detailed reaction mechanism remains poorly understood. In this study, we have combined operando X-ray diffraction and X-ray microscopy along with X-ray tomography, to visualize the evolution of both the morphology and crystal structure of the materials during the entire battery cycling (discharging/charging) process. The dissolution and reformation of sulfur clusters is clearly observed during cycling. In addition, we demonstrate, for the first time, the critical role of current density and temperature in determining the size of both the resulting sulfur clusters and Li2S particles. This study provides new insights about promising avenues for the continued development of lithium sulfur batteries, which we believe may lead to their broad deployment and application.

Computationally efficient dielectric calculations of molecular crystals

The microscopic dielectric response is a key quantity for electronic materials such as organic semiconductors. Calculations of this response for molecular crystals are currently either expensive or rely on extreme simplifications such as multipole expansions which lack microscopic detail. We present an alternate approach using a microscopic analogue of the Clausius-Mossotti equation, which constructs the dielectric response of a crystal from an eigenvalue decomposition of the dielectric response of individual molecules. This method can potentially be used to examine the effects of defects, disorder, and surfaces on the dielectric properties of molecular solids.