Lukas Vlcek

Surface protonation at the rutile (110) interface: explicit incorporation of solvation structure within the refined MUSIC model framework.

The detailed solvation structure at the (110) surface of rutile (alpha-TiO2) in contact with bulk liquid water has been obtained primarily from experimentally verified classical molecular dynamics (CMD) simulations of the ab initio-optimized surface in contact with SPC/E water. The results are used to explicitly quantify H-bonding interactions, which are then used within the refined MUSIC model framework to predict surface oxygen protonation constants. Quantum mechanical molecular dynamics (QMD) simulations in the presence of freely dissociable water molecules produced H-bond distributions around deprotonated surface oxygens very similar to those obtained by CMD with nondissociable SPC/E water, thereby confirming that the less computationally intensive CMD simulations provide accurate H-bond information. Utilizing this H-bond information within the refined MUSIC model, along with manually adjusted Ti-O surface bond lengths that are nonetheless within 0.05 A of those obtained from static density functional theory (DFT) calculations and measured in X-ray reflectivity experiments (as well as bulk crystal values), give surface protonation constants that result in a calculated zero net proton charge pH value (pHznpc) at 25 degrees C that agrees quantitatively with the experimentally determined value (5.4+/-0.2) for a specific rutile powder dominated by the (110) crystal face. Moreover, the predicted pHznpc values agree to within 0.1 pH unit with those measured at all temperatures between 10 and 250 degrees C. A slightly smaller manual adjustment of the DFT-derived Ti-O surface bond lengths was sufficient to bring the predicted pHznpcvalue of the rutile (110) surface at 25 degrees C into quantitative agreement with the experimental value (4.8+/-0.3) obtained from a polished and annealed rutile (110) single crystal surface in contact with dilute sodium nitrate solutions using second harmonic generation (SHG) intensity measurements as a function of ionic strength. Additionally, the H-bond interactions between protolyzable surface oxygen groups and water were found to be stronger than those between bulk water molecules at all temperatures investigated in our CMD simulations (25, 150 and 250 degrees C). Comparison with the protonation scheme previously determined for the (110) surface of isostructural cassiterite (alpha-SnO2) reveals that the greater extent of H-bonding on the latter surface, and in particular between water and the terminal hydroxyl group (Sn-OH) results in the predicted protonation constant for that group being lower than for the bridged oxygen (Sn-O-Sn), while the reverse is true for the rutile (110) surface. These results demonstrate the importance of H-bond structure in dictating surface protonation behavior, and that explicit use of this solvation structure within the refined MUSIC model framework results in predicted surface protonation constants that are also consistent with a variety of other experimental and computational data.

Optimized Unlike-Pair Interactions for Water—Carbon Dioxide Mixtures Described by the SPC/E and EPM2 Models

The unlike-pair interaction parameters for the SPC/E-EPM2 models have been optimized to reproduce the mutual solubility of water and carbon dioxide at the conditions of liquid-supercritical fluid phase equilibria. An efficient global optimization of the parameters is achieved through an implementation of the coupling parameter approach, adapted to phase equilibria calculations in the Gibbs ensemble, that explicitly corrects for the overpolarization of the SPC/E water molecule in the nonpolar CO(2) environments. The resulting H(2)O-CO(2) force field accurately reproduces the available experimental solubilities at the two fluid phases in equilibria as well as the corresponding species tracer diffusion coefficients.

Dynamic and structural properties of room-temperature ionic liquids near silica and carbon surfaces

The dynamic and structural properties of a room-temperature ionic liquid (RTIL) 1-butyl-3-methyl-imidazolium(trifluoromethanesulfonimide) ([C4mim][Tf2N]) confined in silica and carbon mesopores were investigated by molecular dynamics (MD) simulations and nuclear magnetic resonance (NMR) experiments. The complex interfacial microstructures of confined [C4mim][Tf2N] are attributed to the distinctive surface features of the silica mesopore. The temperature-dependent diffusion coefficients of [C4mim][Tf2N] confined in the silica or carbon mesopore exhibit divergent behavior. The loading fraction (f = 1.0, 0.5, and 0.25) has a large effect on the magnitude of the diffusion coefficient in the silica pore and displays weaker temperature dependence as the loading fraction decreases. The diffusion coefficients of mesoporous carbon-confined [C4mim][Tf2N] are relatively insensitive to the loading faction and exhibit a temperature dependence that is similar to the bulk dependence at all loading levels. Such phenomena can be attributed to the unique surface heterogeneity, dissimilar interfacial microstructures, and interaction potential profile of RTILs near silica and carbon walls.

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.

Correspondence between Cluster-Ion and Bulk Solution Thermodynamic Properties: On the Validity of the Cluster-Pair-Based Approximation

Since the single-ion thermodynamic properties of bulk solutions are not directly accessible from experiments, extrapolations have been devised to estimate them from experimental measurements on small-clusters. Extrapolations based on the cluster-pair-based approximation (CPA) technique (Tissandier et al. J. Phys. Chem. A 1998, 102, 7787-7794) and its variants are currently considered one of the most reliable source of single-ion hydration thermodynamic data and have been used as a benchmark for the development of molecular and continuum solvation models. Despite its importance, the CPA has not been thoroughly tested and recent studies have indicated inconsistencies with molecular simulations. The present work challenges the key CPA assumptions that the hydration properties of single cations and anions in growing clusters rapidly converge to each other following a monotonous trend. Using a combination of simulation techniques to study the transition between alkali halide ions in small clusters and bulk solution, we show that this convergence is rather slow and involves a surprising change in trends, which can result in significant errors in the original estimated single-ion properties. When these cluster-size-dependent effects are taken into account, the inconsistencies between molecular models and experimental predictions disappear, and the value of the proton hydration enthalpy based on the CPA aligns with estimates based on other principles.

Comparison of cation adsorption by isostructural rutile and cassiterite.

Macroscopic net proton charging curves for powdered rutile and cassiterite specimens with the (110) crystal face predominant, as a function of pH in RbCl and NaCl solutions, trace SrCl(2) in NaCl, and trace ZnCl(2) in NaCl and Na Triflate solutions, are compared to corresponding molecular-level information obtained from static DFT optimizations and classical MD simulations, as well as synchrotron X-ray methods. The similarities and differences in the macroscopic charging behavior of rutile and cassiterite largely reflect the cation binding modes observed at the molecular level. Cation adsorption is primarily inner-sphere on both isostructural (110) surfaces, despite predictions that outer-sphere binding should predominate on low bulk dielectric constant oxides such as cassiterite (ε(bulk) ≈ 11). Inner-sphere adsorption is also significant for Rb(+) and Na(+) on neutral surfaces, whereas Cl(-) binding is predominately outer-sphere. As negative surface charge increases, relatively more Rb(+), Na(+), and especially Sr(2+) are bound in highly desolvated tetradentate fashion on the rutile (110) surface, largely accounting for enhanced negative charge development relative to cassiterite. Charging curves in the presence of Zn(2+) are very steep but similar for both oxides, reflective of Zn(2+) hydrolysis (and accompanying proton release) during the adsorption process, and the similar binding modes for ZnOH(+) on both surfaces. These results suggest that differences in cation adsorption between high and low bulk dielectric constant oxides are more subtly related to the relative degree of cation desolvation accompanying inner-sphere binding (i.e., more tetradentate binding on rutile), rather than distinct inner- and outer-sphere adsorption modes. Cation desolvation may be favored at the rutile (110) surface in part because inner-sphere water molecules are bound further from and less tightly than on the cassiterite (110) surface. Hence, their removal upon inner-sphere cation binding is relatively more favorable.

Charging Properties of Cassiterite (α-SnO2) 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.

Molecular dynamics simulation of ss-DNA translocation between copper nanoelectrodes incorporating electrode charge dynamics.

Molecular dynamics simulations have been performed to study the translocation of single-stranded (ss)-DNA through the nanoscale gap between the nanoscale electrodes of a proposed genomic sequencing device. Using a fixed gap width between the electrodes and a small sample segment of ss-DNA as initial starting points in this project, the effect of applied electric fields on translocation velocity was studied. To describe the electrostatic interactions of the water, ions, and ss-DNA with the nanoscale electrodes, we applied the electrode charge dynamics (ECD) method. Through the density profile and comparison of translocation velocities to extrapolated experimental data, we found the ECD potential to be a better descriptor of the metal/nonmetal electrostatic interactions compared to the commonly used universal force field (UFF). Translocation velocities obtained using the ECD potential were consistent with simulated bulk data.

Vapor–Liquid Equilibrium and Polarization Behavior of the GCP Water Model: Gaussian Charge-on-Spring versus Dipole Self-Consistent Field Approaches to Induced Polarization

We developed the Gaussian charge-on-spring (GCOS) version of the original self-consistent field implementation of the Gaussian Charge Polarizable water model and test its accuracy to represent the polarization behavior of the original model involving smeared charges and induced dipole moments. For that purpose we adapted the recently proposed multiple-particle-move (MPM) within the Gibbs and isochoric-isothermal ensembles Monte Carlo methods for the efficient simulation of polarizable fluids. We assessed the accuracy of the GCOS representation by a direct comparison of the resulting vapor-liquid phase envelope, microstructure, and relevant microscopic descriptors of water polarization along the orthobaric curve against the corresponding quantities from the actual GCP water model.

Thermodynamic perturbation theory of the second order: Implementation for models with double-bonded sites

Contributions of three classes of graphs contributing to the second order of the thermodynamic perturbation theory have been evaluated for two types of qualitatively different primitive models: (i) MeOH3 model of methanol, and (ii) EPM3 model of water. It is shown that the contributions of linear chains, that is the graphs satisfying the condition of steric incompatibility, bring only a marginal improvement over the first-order theory. The most significant contribution comes from the graph accounting for double bonding of the oxygen site of the models. Neglecting the linear chain diagrams and retaining only this graph, we derive general analytic expressions for the thermodynamic properties of the considered models and find that the theory within this approximation is in agreement with simulation data.