Shawn M. Kathmann

Understanding the Surface Potential of Water

We have resolved the inconsistency in quantifying the surface potential at the liquid-vapor interface when using explicit ab initio electronic charge density and effective atomic partial charge models of liquid water. This is related, in part, to the fact that the resulting electric potentials from partial-charge models and ab initio charge distributions are quite different except for those regions of space between the molecules. We show that the electrostatic surface potential from a quantum mechanical charge distribution compares well to high-energy electron diffraction and electron holography measurements, as opposed to the comparison with electrochemical measurements. We suggest that certain regions of space be excluded when comparing computed surface potentials with electrochemical measurements. This work describes a novel interpretation of ab initio computed surface potentials through high-energy electron holography measurements as useful benchmarks toward a better understanding of electrochemistry.

Spectroscopic Studies of the Phase Transition in Ammonia Borane: Raman spectroscopy of single crystal NH3BH3 as a function of temperature from 88 to 330 K

Raman spectra of single crystal ammonia borane, NH3BH3, were recorded as a function of temperature from 88 to 300 K using Raman microscopy and a variable temperature stage. The orthorhombic to orientationally disordered tetragonal phase transition at 225 K was clearly evident from the decrease in the number of vibrational modes. However, some of the modes in the orthorhombic phase appeared to merge 10-12 K below the phase transition perhaps suggesting the presence of an intermediate phase. Factor group analysis of vibrational spectra for both orthorhombic and tetragonal phase is provided. In addition, electronic structure calculations are used to assist in the interpretation and assignment of the normal modes.

Hybrid approach for free energy calculations with high-level methods: Application to the SN2 reaction of CHCl3 and OH− in water

We present an approach to calculate the free energy profile along a condensed-phase reaction path based on high-level electronic structure methods for the reactive region. The bulk of statistical averaging is shifted toward less expensive descriptions by using a hierarchy of representations that includes molecular mechanics, density functional theory, and coupled cluster theories. As an application of this approach we study the reaction of CHCl3 with OH- in aqueous solution.

Molecular simulations of heterogeneous ice nucleation. I. Controlling ice nucleation through surface hydrophilicity

Ice formation is one of the most common and important processes on earth and almost always occurs at the surface of a material. A basic understanding of how the physicochemical properties of a materials surface affect its ability to form ice has remained elusive. Here, we use molecular dynamics simulations to directly probe heterogeneous ice nucleation at a hexagonal surface of a nanoparticle of varying hydrophilicity. Surprisingly, we find that structurally identical surfaces can both inhibit and promote ice formation and analogous to a chemical catalyst, it is found that an optimal interaction between the surface and the water exists for promoting ice nucleation. We use our microscopic understanding of the mechanism to design a modified surface in silico with enhanced ice nucleating ability.

Snapshots of Proton Accommodation at a Microscopic Water Surface: Understanding the Vibrational Spectral Signatures of the Charge Defect in Cryogenically Cooled H(+)(H2O)(n=2-28) Clusters.

We review the role that gas-phase, size-selected protonated water clusters, H(+)(H2O)n, have played in unraveling the microscopic mechanics responsible for the spectroscopic behavior of the excess proton in bulk water. Because the larger (n ≥ 10) assemblies are formed with three-dimensional cage morphologies that more closely mimic the bulk environment, we report the spectra of cryogenically cooled (10 K) clusters over the size range 2 ≤ n ≤ 28, over which the structures evolve from two-dimensional arrangements to cages at around n = 10. The clusters that feature a complete second solvation shell around a surface-embedded hydronium ion yield spectral signatures of the proton defect similar to those observed in dilute acids. The origins of the large observed shifts in the proton vibrational signature upon cluster growth were explored with two types of theoretical analyses. First, we calculate the cubic and semidiagonal quartic force constants and use these in vibrational perturbation theory calculations to establish the couplings responsible for the large anharmonic red shifts. We then investigate how the extended electronic wave functions that are responsible for the shapes of the potential surfaces depend on the nature of the H-bonded networks surrounding the charge defect. These considerations indicate that, in addition to the sizable anharmonic couplings, the position of the OH stretch most associated with the excess proton can be traced to large increases in the electric fields exerted on the embedded hydronium ion upon formation of the first and second solvation shells. The correlation between the underlying local structure and the observed spectral features is quantified using a model based on Badgers rule as well as via the examination of the electric fields obtained from electronic structure calculations.

Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers

Coarse grained molecular dynamics simulations are presented in which the sensitivity of the ice nucleation rate to the hydrophilicity of a graphene nanoflake is investigated. We find that an optimal interaction strength for promoting ice nucleation exists, which coincides with that found previously for a face centered cubic (111) surface. We further investigate the role that the layering of interfacial water plays in heterogeneous ice nucleation and demonstrate that the extent of layering is not a good indicator of ice nucleating ability for all surfaces. Our results suggest that to be an efficient ice nucleating agent, a surface should not bind water too strongly if it is able to accommodate high coverages of water.

Neutron powder diffraction and molecular simulation study of the structural evolution of ammonia borane from 15 to 340 K.

The structural behavior of (11)B-, (2)H-enriched ammonia borane, ND(3)(11)BD(3), over the temperature range from 15 to 340 K was investigated using a combination of neutron powder diffraction and ab initio molecular dynamics simulations. In the low temperature orthorhombic phase, the progressive displacement of the borane group under the amine group was observed leading to the alignment of the B-N bond near parallel to the c-axis. The orthorhombic to tetragonal structural phase transition at 225 K is marked by dramatic change in the dynamics of both the amine and borane group. The resulting hydrogen disorder is problematic to extract from the metrics provided by Rietveld refinement but is readily apparent in molecular dynamics simulation and in difference Fourier transform maps. At the phase transition, Rietveld refinement does indicate a disruption of one of two dihydrogen bonds that link adjacent ammonia borane molecules. Metrics determined by Rietveld refinement are in excellent agreement with those determined from molecular simulation. This study highlights the valuable insights added by coupled experimental and computational studies.

Analysis of the activation and heterolytic dissociation of H2 by frustrated Lewis pairs: NH3/BX3 (X = H, F, and Cl).

We performed a computational study of H(2) activation and heterolytic dissociation promoted by prototype Lewis acid/base pairs NH(3)/BX(3) (X = H, F, and Cl) to understand the mechanism in frustrated Lewis pairs (FLPs). Although the NH(3)/BX(3) pairs form strong dative bonds, electronic structure theories make it possible to explore the potential energy surface away from the dative complex, in regions relevant to H(2) activation in FLPs. A weakly bound precursor complex, H(3)N·H(2)·BX(3), was found in which the H(2) molecule interacts side-on with B and end-on with N. The BX(3) group is pyramidal in the case of X = H, similar to the geometry of BH(5), but planar in the complexes with X = F and Cl. The latter complexes convert to ion pairs, [NH(4)(+)][BHX(3)(-)] with enthalpy changes of 7.3 and -9.4 kcal/mol, respectively. The minimum energy paths between the FLP and the product ion pair of the chloro and fluoro complexes were calculated and analyzed in great detail. At the transition state (TS), the H(2) bond is weakened and the BX(3) moiety has undergone significant pyramidal distortion. As such, the FLP is prepared to accept the incipient proton and hydride ion on the product-side. The interaction energy of the H(2) with the acid/base pair and the different contributions for the precursor and TS complex from an energy decomposition analysis expose the dominant factors affecting the reactivity. We find that structural reorganization of the precursor complex plays a significant role in the activation and that charge-transfer interactions are the dominant stabilizing force in the activated complex. The electric field clearly has a role in polarizing H(2), but its contribution to the overall interaction energy is small compared to that from the overlap of the p(N), σ(H-H), σ*(H-H), and p(B) orbitals at the TS. Our detailed analysis of the interaction of H(2) with the FLP provides insight into the important components that should be taken into account when designing related systems to activate H(2).

Hydrated Structure of Ag(I) Ion from Symmetry-Dependent, K- and L-Edge XAFS Multiple Scattering and Molecular Dynamics Simulations

Details of the first-shell water structure about Ag(+) are reported from a corefinement of the K- and L(2)-edge multiple scattering signal in the X-ray absorption fine structure (XAFS) spectra. Detailed fits of the Ag K-edge data that include the contributions from multiple scattering processes in the hydrated ion structure cannot distinguish between models containing tetrahedral symmetry versus those containing collinear O-Ag-O bonds. However, we show that the multiple scattering oscillations at the L(2)-edges have distinctly different phase and amplitude functions than at the K-edge. These phase and amplitude functions depend not only on the symmetry of the multiple scattering paths but also on the nature of the final state electronic wave function probed by the dipole-allowed transition. Hence the multiple scattering portions of K- and L(2)-edge spectra provide independent measurements of the local symmetry--not a redundant measurement as is commonly believed. On the basis of the enhanced information content obtained by the simultaneous assessment of both the K- and L(2)-edges, we report that the hydrated Ag(+) structure contains five or six water molecules in the first shell with a significant number of nearly collinear and 90 degrees O-Ag-O bond angles. Finally, the K- and L(2)-edge spectra are used to benchmark the hydration structure that is generated from both DFT-based and classical molecular dynamics simulations. Simulated first-shell structures are compared to the experimental structures.

Non-hexagonal ice at hexagonal surfaces: the role of lattice mismatch

It has long been known that ice nucleation usually proceeds heterogeneously on the surface of a foreign body. However, little is known at the microscopic level about which properties of a material determine its effectiveness at nucleating ice. This work focuses on the long standing, conceptually simple, view on the role of a good crystallographic match between bulk ice and the underlying substrate. We use grand canonical Monte Carlo to generate the first overlayer of water at the surface and find that the traditional view of heterogeneous nucleation does not adequately account for the array of structures that water may form at the surface. We find that, in order to describe the structures formed, a good match between the substrate and the nearest neighbour oxygen-oxygen distance is a better descriptor than a good match to the bulk ice lattice constant.