Resolving local configurational contributions to X-ray and neutron radial distribution functions within solutions of concentrated electrolytes - a case study of concentrated NaOH.

Extreme conditions of complex materials often lead to a manifold of local environments that challenge characterization and require new advances at the intersection of modern experimental and theoretical techniques. In this contribution, highly caustic and viscous aqueous NaOD solutions were characterized with a combination of X-ray and neutron radial distribution function (RDF) analyses, molecular dynamics simulations and sub-ensemble analysis. While this system has been the topic of some study, the current work expands upon the state of knowledge regarding the extent to which water is perturbed within this chemically extreme solution. Further, we introduce analyses that goes beyond merely identifying the different local environments (ion solvation and coordination environments) that are present, but toward understanding their relative contributions to the ensemble solution RDF. This integrated approach yields unique insight into the experimental sensitivity of RDFs to changes in local geometries, the composition of solvation environments about ions, and the challenge of experimentally differentiating the ensemble of all superimposed local environments-a feature of increasing importance within the extreme condition of high ionic strength.

Computational Insights into Materials and Interfaces for Capacitive Energy Storage.

Supercapacitors such as electric double-layer capacitors (EDLCs) and pseudocapacitors are becoming increasingly important in the field of electrical energy storage. Theoretical study of energy storage in EDLCs focuses on solving for the electric double-layer structure in different electrode geometries and electrolyte components, which can be achieved by molecular simulations such as classical molecular dynamics (MD), classical density functional theory (classical DFT), and Monte-Carlo (MC) methods. In recent years, combining first-principles and classical simulations to investigate the carbon-based EDLCs has shed light on the importance of quantum capacitance in graphene-like 2D systems. More recently, the development of joint density functional theory (JDFT) enables self-consistent electronic-structure calculation for an electrode being solvated by an electrolyte. In contrast with the large amount of theoretical and computational effort on EDLCs, theoretical understanding of pseudocapacitance is very limited. In this review, we first introduce popular modeling methods and then focus on several important aspects of EDLCs including nanoconfinement, quantum capacitance, dielectric screening, and novel 2D electrode design; we also briefly touch upon pseudocapactive mechanism in RuO2. We summarize and conclude with an outlook for the future of materials simulation and design for capacitive energy storage.

Calorimetric study of alkali and alkaline-earth cation adsorption and exchange at the quartz-solution interface.

Cations in natural solutions significantly impact interfacial processes, particularly dissolution and surface charge measurements for quartz and silica, which are amongst the most naturally abundant and technologically important solids. Thermodynamic parameters for cation-specific interfacial reactions have heretofore been mostly derived instead of directly measured experimentally. This work investigates the energetics of adsorption and exchange reactions of alkali metal (M+) and alkaline earth (M2+) cations with the quartz surface by flow adsorption microcalorimetry, in tandem with in-situ pH measurements. The magnitudes of the heats of adsorption and exchange were found to increase along the Hofmeister series i.e., Li+

Ion Exchange Thermodynamics at the Rutile-Water Interface: Flow Microcalorimetric Measurements and Surface Complexation Modeling of Na-K-Rb-Cl-NO3 Adsorption.

Flow microcalorimetry was used to investigate the energetics associated with Rb+, K+, Na+, Cl-, and NO3- exchange at the rutile-water interface. Heats of exchange reflected differences in bulk hydration/dehydration enthalpies (Na+ > K+ > Rb+, and Cl- > NO3-) such that exchanging Na+ or Cl- from the surface was exothermic, reflecting their greater bulk hydration enthalpies. Exchange heats were measured at pH 2, 3.25, 5.8, and 11 and exhibited considerable differences as well as pH dependence. These trends were rationalized with the aid of a molecularly constrained surface complexation model (SCM) that incorporated the inner-sphere binding observed for the cations on the rutile (110) surface. Explicitly accounting for the inner-sphere binding configuration differences between Rb+, K+, and Na+, as well as accompanying differences in negative surface charge development, resulted in much better agreement with measured exchange ratios than by considering bulk hydration enthalpies alone. The observation that calculated exchange ratios agreed with those measured experimentally lends additional credence to the SCM. Consequently, flow microcalorimetry and surface complexation modeling are a useful complement of techniques for probing the energetics associated with ion exchange and adsorption processes and should also serve to help validate molecular simulations of interfacial energetics.

Ion distribution and selectivity of ionic liquids in microporous electrodes.

The energy density of an electric double layer capacitor, also known as supercapacitor, depends on ion distributions in the micropores of its electrodes. Herein we study ion selectivity and partitioning of symmetric, asymmetric, and mixed ionic liquids among different pores using the classical density functional theory. We find that a charged micropore in contact with mixed ions of the same valence is always selective to the smaller ions, and the ion selectivity, which is strongest when the pore size is comparable to the ion diameters, drastically falls as the pore size increases. The partitioning behavior in ionic liquids is fundamentally different from those corresponding to ion distributions in aqueous systems whereby the ion selectivity is dominated by the surface energy and entropic effects insensitive to the degree of confinement.

Precise determination of water exchanges on a mineral surface.

Solvent exchanges on solid surfaces and dissolved ions are a fundamental property important for understanding chemical reactions, but the rates of fast exchanges are poorly constrained. We probed the diffusional motions of water adsorbed onto nanoparticles of the mineral barite (BaSO4) using quasi-elastic neutron scattering (QENS) and classical molecular dynamics (MD) to reveal the complex dynamics of water exchange along mineral surfaces. QENS data as a function of temperature and momentum transfer (Q) were fit using scattering functions derived from MD trajectories. The simulations reproduce the dynamics measured in the experiments at ambient temperatures, but as temperature is lowered the simulations overestimate slower motions. Decomposition of the MD-computed QENS intensity into contributions from adsorbed and unbound water shows that the majority of the signal arises from adsorbed species, although the dynamics of unbound water cannot be dismissed. The mean residence times of water on each of the four surface sites present on the barite {001} were calculated using MD: at room temperature the low barium site is 194 ps, whereas the high barium site contains two distributions of motions at 84 and 2.5 ps. These contrast to 13 ps residence time on both sulfate sites, with an additional surface diffusion exchange of 66 ps. Surface exchanges are similar to those of the aqueous ions calculated using the same force field: Baaq2+ is 208 ps and SO4aq2- is 5.8 ps. This work demonstrates how MD can be a reliable method to deconvolute solvent exchange reactions when quantitatively validated by QENS measurements.

Molecular Origins of the Zeta Potential.

The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO2) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to -0.4 C/m2), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl2 aqueous solutions, were modeled and compared to experimental ZPs for TiO2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that "anomalous" ZP values and inequalities between the point of zero charge derived from electrokinetic versus pH titration measurements both arise from differing co- and counterion sorption affinities. We show that at the molecular level the ZP arises from the delicate interplay of spatially varying dynamics, structure, and electrostatics in a narrow interfacial region within about 15 Å of the surface, even in dilute salt solutions. This contrasts fundamentally with continuum descriptions of such interfaces, which predict the ZP response region to be inversely related to ionic strength. In reality the properties of this interfacial region are dominated by relatively immobile and structured water. Consequently, viscosity values are substantially greater than in the bulk, and electrostatic potential profiles are oscillatory in nature.

Quantum Tunneling of Water in Beryl: A New State of the Water Molecule.

Using neutron scattering and ab initio simulations, we document the discovery of a new "quantum tunneling state" of the water molecule confined in 5 Å channels in the mineral beryl, characterized by extended proton and electron delocalization. We observed a number of peaks in the inelastic neutron scattering spectra that were uniquely assigned to water quantum tunneling. In addition, the water proton momentum distribution was measured with deep inelastic neutron scattering, which directly revealed coherent delocalization of the protons in the ground state.

Aqueous proton transfer across single-layer graphene.

Proton transfer across single-layer graphene proceeds with large computed energy barriers and is therefore thought to be unfavourable at room temperature unless nanoscale holes or dopants are introduced, or a potential bias is applied. Here we subject single-layer graphene supported on fused silica to cycles of high and low pH, and show that protons transfer reversibly from the aqueous phase through the graphene to the other side where they undergo acid-base chemistry with the silica hydroxyl groups. After ruling out diffusion through macroscopic pinholes, the protons are found to transfer through rare, naturally occurring atomic defects. Computer simulations reveal low energy barriers of 0.61-0.75 eV for aqueous proton transfer across hydroxyl-terminated atomic defects that participate in a Grotthuss-type relay, while pyrylium-like ether terminations shut down proton exchange. Unfavourable energy barriers to helium and hydrogen transfer indicate the process is selective for aqueous protons.

Structure and stability of SnO2 nanocrystals and surface-bound water species.

The structure of SnO2 nanoparticles (avg. 5 nm) with a few layers of water on the surface has been elucidated by atomic pair distribution function (PDF) methods using in situ neutron total scattering data and molecular dynamics (MD) simulations. Analysis of PDF, neutron prompt gamma, and thermogravimetric data, coupled with MD-generated surface D2O/OD configurations demonstrates that the minimum concentration of OD groups required to prevent rapid growth of nanoparticles during thermal dehydration corresponds to ~0.7 monolayer coverage. Surface hydration layers not only stabilize the SnO2 nanoparticles but also induce particle-size-dependent structural modifications and are likely to promote interfacial reactions through hydrogen bonds between adjacent particles. Upon heating/dehydration under vacuum above 250 °C, nanoparticles start to grow with low activation energies, rapid increase of nanoparticle size, and a reduction in the a lattice dimension. This study underscores the value of neutron diffraction and prompt-gamma analysis, coupled with molecular modeling, in elucidating the influence of surface hydration on the structure and metastable persistence of oxide nanomaterials.

High-pressure cell for neutron reflectometry of supercritical and subcritical fluids at solid interfaces.

A new high-pressure cell design for use in neutron reflectometry (NR) for pressures up to 50 MPa and a temperature range of 300-473 K is described. The cell design guides the neutron beam through the working crystal without passing through additional windows or the bulk fluid, which provides for a high neutron transmission, low scattering background, and low beam distortion. The o-ring seal is suitable for a wide range of subcritical and supercritical fluids and ensures high chemical and pressure stability. Wafers with a diameter of 5.08 cm (2 in.) and 5 mm or 10 mm thickness can be used with the cells, depending on the required pressure and momentum transfer range. The fluid volume in the sample cell is very small at about 0.1 ml, which minimizes scattering background and stored energy. The cell design and pressure setup for measurements with supercritical fluids are described. NR data are shown for silicon/silicon oxide and quartz wafers measured against air and subsequently within the high-pressure cell to demonstrate the neutron characteristics of the high-pressure cell. Neutron reflectivity data for supercritical CO(2) in contact with quartz and Si/SiO(2) wafers are also shown.

Direct measurements of pore fluid density by vibrating tube densimetry.

The densities of pore-confined fluids were measured for the first time by means of vibrating tube densimetry (VTD). A custom-built high-pressure, high-temperature vibrating tube densimeter was used to measure the densities of propane at subcritical and supercritical temperatures (between 35 and 97 °C) and carbon dioxide at supercritical temperatures (between 32 and 50 °C) saturating hydrophobic silica aerogel (0.2 g/cm(3), 90% porosity) synthesized inside Hastelloy U-tubes. Additionally, supercritical isotherms of excess adsorption for CO(2) and the same porous solid were measured gravimetrically using a precise magnetically coupled microbalance. Pore fluid densities and total adsorption isotherms increased monotonically with increasing density of the bulk fluid, in contrast to excess adsorption isotherms, which reached a maximum and then decreased toward zero or negative values above the critical density of the bulk fluid. The isotherms of confined fluid density and excess adsorption obtained by VTD contain additional information. For instance, the maxima of excess adsorption occur below the critical density of the bulk fluid at the beginning of the plateau region in the total adsorption, marking the end of the transition of pore fluid to a denser, liquidlike pore phase. Compression of the confined fluid significantly beyond the density of the bulk fluid at the same temperature was observed even at subcritical temperatures. The effect of pore confinement on the liquid-vapor critical temperature of propane was less than ~1.7 K. The results for propane and carbon dioxide showed similarity in the sense of the principle of corresponding states. Good quantitative agreement was obtained between excess adsorption isotherms determined from VTD total adsorption results and those measured gravimetrically at the same temperature, confirming the validity of the vibrating tube measurements. Thus, it is demonstrated that vibrating tube densimetry is a novel experimental approach capable of providing directly the average density of pore-confined fluids, and hence complementary to the conventional gravimetric or volumetric/piezometric adsorption techniques, which yield the excess adsorption (the Gibbsian surface excess).

Neutron reflectivity study of substrate surface chemistry effects on supported phospholipid bilayer formation on (11 Ì 20) sapphire.

Oxide-supported phospholipid bilayers (SPBs) used as biomimetic membranes are significant for a broad range of applications including improvement of biomedical devices and biosensors, and in understanding biomineralization processes and the possible role of mineral surfaces in the evolution of pre-biotic membranes. Continuous-coverage and/or stacked SPBs retain properties (e.g., fluidity) more similar to native biological membranes, which is desirable for most applications. Using neutron reflectivity, we examined the role of oxide surface charge (by varying pH and ionic strength) and of divalent Ca(2+) in controlling surface coverage and potential stacking of dipalmitoylphosphatidylcholine (DPPC) bilayers on the (11 ̅20) face of sapphire (α-Al(2)O(3)). Nearly full bilayers were formed at low to neutral pH, when the sapphire surface is positively charged, and at low ionic strength (I=15 mM NaCl). Coverage decreased at higher pH, close to the isoelectric point of sapphire, and also at high I≥210 mM, or with addition of 2mM Ca(2+). The latter two effects are not additive, suggesting that Ca(2+) mitigates the effect of higher I. These trends agree with previous results for phospholipid adsorption on α-Al(2)O(3) particles determined by adsorption isotherms and on single-crystal (10 ̅10) sapphire by atomic force microscopy, suggesting consistency of oxide surface chemistry-dependent effects across experimental techniques.

Diffusion processes in water on oxide surfaces: quasielastic neutron scattering study of hydration water in rutile nanopowder.

Quasielastic neutron scattering (QENS) was used to investigate the diffusion dynamics of hydration water on the surface of rutile (TiO(2)) nanopowder. The dynamics measurements utilizing two inelastic instruments, a backscattering spectrometer and a disk chopper spectrometer, probed the fast, intermediate, and slow motions of the water molecules on the time scale of picoseconds to more than a nanosecond. We employed a model-independent analysis of the data collected at each value of the scattering momentum transfer to investigate the temperature dependence of several diffusion components. All of the probed components were present in the studied temperature range of 230-320 K, providing, at a first sight, no evidence of discontinuity in the hydration water dynamics. However, a qualitative change in the elastic scattering between 240 and 250 K suggested a surface freezing-melting transition, when the motions that were localized at lower temperatures became delocalized at higher temperatures. On the basis of our previous molecular dynamics simulations of this system, we argue that interpretation of QENS data from such a complex interfacial system requires at least qualitative input from simulations, particularly when comparing results from spectrometers with very different energy resolutions and dynamic ranges.

Faster proton transfer dynamics of water on SnO2 compared to TiO2.

Proton jump processes in the hydration layer on the iso-structural TiO(2) rutile (110) and SnO(2) cassiterite (110) surfaces were studied with density functional theory molecular dynamics. We find that the proton jump rate is more than three times faster on cassiterite compared with rutile. A local analysis based on the correlation between the stretching band of the O-H vibrations and the strength of H-bonds indicates that the faster proton jump activity on cassiterite is produced by a stronger H-bond formation between the surface and the hydration layer above the surface. The origin of the increased H-bond strength on cassiterite is a combined effect of stronger covalent bonding and stronger electrostatic interactions due to differences of its electronic structure. The bridging oxygens form the strongest H-bonds between the surface and the hydration layer. This higher proton jump rate is likely to affect reactivity and catalytic activity on the surface. A better understanding of its origins will enable methods to control these rates.

Charging properties of cassiterite (alpha-SnO(2)) surfaces in NaCl and RbCl ionic media.

The acid-base properties of cassiterite (alpha-SnO2) surfaces at 10-50 degrees C were studied using potentiometric titrations of powder suspensions in aqueous NaCl and RbCl media. The proton sorption isotherms exhibited common intersection points in the pH range of 4.0-4.5 under all conditions, and the magnitude of charging was similar but not identical in NaCl and RbCl. The hydrogen bonding configuration at the oxide-water interface, obtained from classical molecular dynamics (MD) simulations, was analyzed in detail, and the results were explicitly incorporated in calculations of protonation constants for the reactive surface sites using the revised MUSIC model. The calculations indicated that the terminal SnOH2 group is more acidic than the bridging Sn2OH group, with protonation constants (log KH) of 3.60 and 5.13 at 25 degrees C, respectively. This is contrary to the situation on the isostructural alpha-TiO2 (rutile), apparently because of the difference in electronegativity between Ti and Sn. MD simulations and speciation calculations indicated considerable differences in the speciation of Na+ and Rb+, despite the similarities in overall charging. Adsorbed sodium ions are almost exclusively found in bidentate surface complexes, whereas adsorbed rubidium ions form comparable numbers of bidentate and tetradentate complexes. Also, the distribution of adsorbed Na+ between the different complexes shows a considerable dependence on the surface charge density (pH), whereas the distribution of adsorbed Rb+ is almost independent of pH. A surface complexation model (SCM) capable of accurately describing both the measured surface charge and the MD-predicted speciation of adsorbed Na+/Rb+ was formulated. According to the SCM, the deprotonated terminal group (SnOH(-0.40)) and the protonated bridging group (Sn2OH+0.36) dominate the surface speciation over the entire pH range of this study (2.7-10). The complexation of medium cations increases significantly with increasing negative surface charge, and at pH 10, roughly 40% of the terminal sites are predicted to form cation complexes, whereas anion complexation is minor throughout the studied pH range.

Magnesium silicate dissolution investigated by 29Si MAS, 1H-29Si CPMAS, 25Mg QCPMG, and 1H-25Mg CP QCPMG NMR.

Olivine-(Mg,Fe)(2)SiO(4)-has been the subject of frequent investigation in the earth sciences because of its simple structure and rapid dissolution kinetics. Several studies have observed a preferential release of the divalent cation with respect to silicon during weathering under acidic conditions, which has been correlated to the formation of a silicon-rich leached layer. While leached layer formation has been inferred through the changing solution chemistry, a thorough spectroscopic investigation of olivine reacted under acidic conditions has not been conducted. The pure magnesium end member of the olivine series (forsterite-Mg(2)SiO(4)) was chosen for detailed investigations in this study because paramagnetic iron hinders NMR investigations by providing an extra mode of relaxation for neighboring nuclei, causing lineshapes to become significantly broadened and unobservable in the NMR spectrum. For reacting forsterite, spectroscopic interrogations using nuclear magnetic resonance (NMR) can elucidate the changing magnesium coordination and bonding environment. In this study, we combine analysis of the changing solution chemistry with advanced NMR techniques ((29)Si MAS, (1)H-(29)Si CP MAS, (25)Mg QCPMG, and (1)H-(25)Mg CP QCPMG NMR) to probe leached layer formation and secondary phase precipitation during the dissolution of forsterite at 150 degrees C.

Electrophoretic study of the SnO2/aqueous solution interface up to 260 degrees C.

An electrophoresis cell developed in our laboratory was utilized to determine the zeta potential at the SnO(2) (cassiterite)/aqueous solution (10(-3) mol kg(-1) NaCl) interface over the temperature range from 25 to 260 degrees C. Experimental techniques and methods for the calculation of zeta potential at elevated temperature are described. From the obtained zeta potential data as a function of pH, the isoelectric points (IEPs) of SnO(2) were obtained for the first time. From these IEP values, the standard thermodynamic functions were calculated for the protonation-deprotonation equilibrium at the SnO(2) surface, using the 1-pK surface complexation model. It was found that the IEP values for SnO(2) decrease with increasing temperature, and this behavior is compared to the predicted values by the multisite complexation (MUSIC) model and other semitheoretical treatments, and were found to be in excellent agreement.

Suppression of the dynamic transition in surface water at low hydration levels: a study of water on rutile.

Our quasielastic neutron-scattering experiments and molecular-dynamics simulations probing surface water on rutile (TiO2) have demonstrated that a sufficiently high hydration level is a prerequisite for the temperature-dependent crossover in the nanosecond dynamics of hydration water. Below the monolayer coverage of mobile surface water, a weak temperature dependence of the relaxation times with no apparent crossover is observed. We associate the dynamic crossover with interlayer jumps of the mobile water molecules, which become possible only at a sufficiently high hydration level.