Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Temperature-dependent slip length for water and electrolyte solution

HYPOTHESIS

The temperature dependence of boundary slip at liquid-solid interface is critical both for the fundamental theory and applications of fluid mechanics on micro and nanoscale, such as sustainable cooling of electronic devices. However, there is a controversy on the temperature dependence of boundary slip which lacks experimental evidence, we aim to resolve it by hypothesizing that the temperature dependent slip length depends on the variation in the interfacial energy barrier.

EXPERIMENTS

Here, we measured l_s - T relation of water and NaCl solution on self-assembled FDTS (Perfluorodecyltrichlorosilane) surface using colloidal probe AFM. The transition of l_s - T monotonicity is found. For water and 0.1 M NaCl solution, l_s is negatively correlated with T, while for 1 M NaCl solution, l_s is positively correlated with T.

FINDINGS

Together with molecular dynamics simulations, such observation is quantitatively explained with an analytical model based on rate theory, where the l_s - T monotonicity depends on the difference between liquid-solid interfacial energy barrier and liquid internal energy barrier. Our results provide not only solid experimental evidence for the boundary slip being a rate process, but also a basis for the thermal-hydrodynamic design of microfluidic and nanofluidic devices.

Accurate estimation of dynamical quantities for nonequilibrium nanoscale systems

Fluctuations of dynamical quantities are fundamental and inevitable. For the booming research in nanotechnology, huge relative fluctuation comes with the reduction of system size, leading to large uncertainty for the estimates of dynamical quantities. Thus, increasing statistical efficiency, i.e., reducing the number of samples required to achieve a given accuracy, is of great significance for accurate estimation. Here we propose a theory as a fundamental solution for such problem by constructing auxiliary path for each real path. The states on auxiliary paths constitute canonical ensemble and share the same macroscopic properties (NVT) with the initial states of the real path. By implementing the theory in molecular dynamics simulations, we obtain a nanoscale Couette flow field with an accuracy of 0.2μm/s with relative standard error <0.1. The required number of samples is reduced by 12 orders compared to conventional method. The predicted thermolubric behavior of water sliding on a self-assembled surface is directly validated by experiment under the same velocity. This theory only assumes the system is initially in thermal equilibrium, then driven from that equilibrium by an external perturbation. It could serve as a general approach for extracting the accurate estimate of dynamical quantities from large fluctuations to provide insights on atomic level under experimental conditions, and benefit the studies on mass transport through (biological) nanochannels and fluid film lubrication of nanometer thickness.

Hydration Structures on γ-Alumina Surfaces With and Without Electrolytes Probed by Atomistic Molecular Dynamics Simulations

A wide range of systems, both engineered and natural, feature aqueous electrolyte solutions at interfaces. In this study, the structure and dynamics of water at the two prevalent crystallographic terminations of gamma-alumina, [110] and [100], and the influence of salts─sodium chloride, ammonium acetate, barium acetate, and barium nitrate on such properties─were investigated using equilibrium molecular dynamics simulations. The resulting interfacial phenomena were quantified from simulation trajectories via atomic density profiles, angle probability distributions, residence times, 2-D density distributions within the hydration layers, and hydrogen bond density profiles. Analysis and interpretation of the results are supported by simulation snapshots. Taken together, our results show stronger interaction and closer association of water with the [110] surface, compared to [100], while ion-induced disruption of interfacial water structure was more prevalent at the [100] surface. For the latter, a stronger association of cations is observed, namely sodium and ammonium, and ion adsorption appears determined by their size. The differences in surface-water interactions between the two terminations are linked to their respective surface features and distributions of surface groups, with atomistic-scale roughness of the [110] surface promoting closer association of interfacial water. The results highlight the fundamental role of surface characteristics in determining surface-water interactions, and the resulting effects on ion-surface and ion-water interactions. Since the two terminations of gamma-alumina considered represent interfaces of significance to numerous industrial applications, the results provide insights relevant for catalyst preparation and adsorption-based water treatment, among other applications.

Breakthrough pressure of oil displacement by water through the ultra-narrow kerogen pore throat from the Young-Laplace equation and molecular dynamic simulations

As one common unconventional reservoir, shale plays a pivotal role to compensate the depletion of conventional oil resources. There are numerous nanoscale pores and ultra-narrow pore throats (sub 2 nm) in shale media. To displace oil through ultra-narrow pore throats by water, one needs to overcome excessively-high capillary pressure. Understanding the water-oil two-phase displacement process through pore throats is critical to numerical simulation on tight/shale oil exploitation and ultimate oil recovery estimation. In this work, we use molecular dynamics simulations to investigate oil (represented by n-octane) displacement by water through a ~2 nm kerogen (represented by Type II-C kerogen) pore throat. Besides, the applicability of the Young-Laplace equation to the ultra-narrow kerogen pore throat has been assessed. We find that although the Type II-C kerogen is generally oil-wet, water has an excellent displacement efficiency without the oil film on the substrate, thanks to the hydrogen bonding formed between water and heteroatoms (such as O, N, and S) on the kerogen surface. Unlike previous studies, the capillary pressure obtained from the widely used Young-Laplace equation shows a good agreement with the breakthrough pressure obtained from MD simulations for the ∼2 nm kerogen pore throat. Our work indicates that explicitly considering intermolecular interactions as well as atomistic and molecular level characteristics is imperative to study the two-phase displacement process through ultra-narrow pore throats.

Fluorination to enhance superlubricity performance between self-assembled monolayer and graphite in water

HYPOTHESIS

Achievement of superlubricity is an effective method to reduce friction and wear, which has a prominent influence on the operational efficiency and lifetime of a device. However, some burning issues still remain to be solved for the practical applications of superlubricity, such as the poor load-bearing capacity, especially in liquid superlubricity. Therefore, exploring an effective method to enhance the superlubricity performance is essential to accelerate the application of superlubricity.

EXPERIMENTS

The friction properties between two different self-assembled monolayers (SAMs)-a perfluorocarbon SAM and a hydrocarbon SAM-and graphite in water were explored and compared by atomic force microscopy (AFM).

FINDINGS

Enhanced superlubricity performance due to the fluorination was observed. Specifically, we observed an approximately 85% reduction of the friction coefficient after fluorination, and superlubricity was achieved with extremely low friction coefficient of 0.0003. Moreover, 2.4-fold greater load-bearing capacity of the superlubricity was obtained after fluorination. The molecular origin of the superlubricity enhancement by fluorination was revealed by molecular dynamics (MD) simulations, indicating that the greater load-bearing capacity of the perfluorocarbon SAM was ascribed to the enhanced interaction between the water and SAM by fluorination to form a more robust layered water structure confined in the contact zone, which played a pivotal role in the superlubricity.

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D₂O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

PDADMAC/PSS Oligoelectrolyte Multilayers: Internal Structure and Hydration Properties at Early Growth Stages from Atomistic Simulations

We analyze the internal structure and hydration properties of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) oligoelectrolyte multilayers at early stages of their layer-by-layer growth process. Our study is based on large-scale molecular dynamics simulations with atomistic resolution that we presented recently [Sánchez et al., Soft Matter2019, 15, 9437], in which we produced the first four deposition cycles of a multilayer obtained by alternate exposure of a flat silica substrate to aqueous electrolyte solutions of such polymers at 0.1M of NaCl. In contrast to any previous work, here we perform a local structural analysis that allows us to determine the dependence of the multilayer properties on the distance to the substrate. We prove that the large accumulation of water and ions next to the substrate observed in previous overall measurements actually decreases the degree of intrinsic charge compensation, but this remains as the main mechanism within the interface region. We show that the range of influence of the substrate reaches approximately 3 nm, whereas the structure of the outer region is rather independent from the position. This detailed characterization is essential for the development of accurate mesoscale models able to reach length and time scales of technological interest.

Structure and sum-frequency generation spectra of water on uncharged Q4 silica surfaces: a molecular dynamics study

The structural characteristics and sum-frequency generation (SFG) spectra of water near neutral Q₄ silica surfaces were investigated by molecular dynamics simulations. The interactions of water molecules with atoms of the solid were described by different potential models, in particular by the CLAYFF [Cygan et al., J. Phys. Chem. B, 2004, 108, 1255] and INTERFACE [Heinz et al. Langmuir, 2013, 29, 1754] force fields. The calculations of the contact angle of water have shown that the silica surface modeled with CLAYFF behaves as macroscopically hydrophilic, in contrast to the surface described with the INTERFACE model. The hydrophilicity of CLAYFF stems from too attractive electrostatic surface-water interactions. Regardless of the surface's affinity for water, the aqueous phase has a layered structure in the direction perpendicular to the surface with density fluctuations decaying within a distance of 10 Å from the surface. The orientational ordering of H₂O molecules was found to be more short-range than the density fluctuations, especially for the hydrophobic surfaces. Modeling the SFG spectra has shown that the spectra of all studied hydrophobic silica-water interfaces are similar and have features in common with the spectrum of the water-vapor interface. The spectra fairly agree with experimental results obtained for the silica-water interface at low pH conditions [Myalitsin et al., J. Phys. Chem. C, 2016, 120, 9357]. The spectral response for the hydrophobic interface was computed to primarily arise from the topmost molecules of the first layer of interfacial water. In contrast, the SFG signal from the hydrophilic silica-water interface is accumulated over a greater distance extending for several water layers due to more long-range perturbation of the structure by the surface.

Influence of surface hydrophilicity and hydration on the rotational relaxation of supercooled water on graphene oxide surfaces

The presence of a bulk water film influences the dynamical transitions of supercooled water on graphene oxide surfaces.

Atomistic simulation of PDADMAC/PSS oligoelectrolyte multilayers: overall comparison of tri- and tetra-layer systems

By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1 M NaCl electrolyte solutions and the swollen structures are then subsequently exposed to varying added salt concentration. We investigated the microscopic properties of the films, analyzing in detail the differences between three- and four-layer systems. Our simulations provide insights into the early stages of growth of a multilayer, which are particularly challenging for experimental observations. We found rather strong complexation of the oligoelectrolytes, with fuzzy layering of the film structure. The main charge compensation mechanism is for all cases intrinsic, whereas extrinsic compensation is relatively enhanced for the layer of the last deposition cycle. In addition, we quantified other fundamental observables of these systems, such as the film thickness, water uptake, and overcharge fractions for each deposition layer.

Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

Origin of hydration lubrication of zwitterions on graphene

Formation of a hydration layer on charge sites can support normal pressure, and meanwhile it retains excellent fluidity to provide efficient boundary lubrication; however, it is limited to the sliding system between two similarly charged surfaces. In the present study, we report extremely low friction as the zwitterions in a lipid bilayer slide on the topmost graphene layer of graphite across pure water, with the friction coefficient falling to the level of 0.001, which provides direct evidence that hydration lubrication is effective even between such dissimilar surfaces. The origin of hydration lubrication on graphene was studied by atomic force microscopy and molecular dynamics simulation simultaneously. It reveals that a subnanometer hydration layer is confined between zwitterions and graphene, which remains as a liquid phase under normal pressure. The shear occurs between water molecules and graphene because of the extremely low shear strength of the water/graphene interface, which contributes to extremely low friction. Our finding demonstrates that the formation of a hydration layer is possible to lubricate layered materials efficiently, which has potential implications for designing efficient boundary lubrication with layered materials.

Enhanced Ion Adsorption on Mineral Nanoparticles

Classical molecular dynamics simulation was used to study the adsorption of Na⁺, Ca²⁺, Ba²⁺, and Cl^- ions on gibbsite edge (1 0 0), basal (0 0 1), and nanoparticle (NP) surfaces. The gibbsite NP consists of both basal and edge surfaces. Simulation results indicate that Na⁺ and Cl^- ions adsorb on both (1 0 0) and (0 0 1) surfaces as inner-sphere species (i.e., no water molecules between an ion and the surface). Outer-sphere Cl^- ions (i.e., one water molecule between an ion and the surface) were also found on these surfaces. On the (1 0 0) edge, Ca²⁺ ions adsorb as inner-sphere and outer-sphere complexes, whereas on the (0 0 1) surface, outer-sphere Ca²⁺ ions are the dominant species. Ba²⁺ ions were found as inner-sphere and outer-sphere complexes on both surfaces. Calculated ion surface coverages indicate that, for all ions, surface coverages are always higher on the basal surface compared to those on the edge surface. More importantly, surface coverages for cations on the gibbsite NP are always higher than those calculated for the (1 0 0) and (0 0 1) surfaces. This enhanced ion adsorption behavior for the NP is due to the significant number of inner-sphere cations found at NP corners. Outer-sphere cations do not contribute to the enhanced surface coverage. In addition, there is no ion adsorption enhancement observed for the Cl^- ion. Our work provides a molecular-scale understanding of the relative significance of ion adsorption onto gibbsite basal versus edge surfaces and demonstrates the corner effect on ion adsorption on NPs.

Propane-Water Mixtures Confined within Cylindrical Silica Nanopores: Structural and Dynamical Properties Probed by Molecular Dynamics

Despite the multiple length and time scales over which fluid-mineral interactions occur, interfacial phenomena control the exchange of matter and impact the nature of multiphase flow, as well as the reactivity of C-O-H fluids in geologic systems. In general, the properties of confined fluids, and their influence on porous geologic phenomena are much less well understood compared to those of bulk fluids. We used equilibrium molecular dynamics simulations to study fluid systems composed of propane and water, at different compositions, confined within cylindrical pores of diameter ∼16 Å carved out of amorphous silica. The simulations are conducted within a single cylindrical pore. In the simulated system all the dangling silicon and oxygen atoms were saturated with hydroxyl groups and hydrogen atoms, respectively, yielding a total surface density of 3.8 -OH/nm2. Simulations were performed at 300 K, at different bulk propane pressures, and varying the composition of the system. The structure of the confined fluids was quantified in terms of the molecular distribution of the various molecules within the pore as well as their orientation. This allowed us to quantify the hydrogen bond network and to observe the segregation of propane near the pore center. Transport properties were quantified in terms of the mean square displacement in the direction parallel to the pore axis, which allows us to extract self-diffusion coefficients. The diffusivity of propane in the cylindrical pore was found to depend on pressure, as well as on the amount of water present. It was found that the propane self-diffusion coefficient decreases with increasing water loading because of the formation of water bridges across the silica pores, at sufficiently high water content, which hinder propane transport. The rotational diffusion, the lifespan of hydrogen bonds, and the residence time of water molecules at contact with the silica substrate were quantified from the simulated trajectories using the appropriate autocorrelation functions. The simulations contribute to a better understanding of the molecular phenomena relevant to the behavior of fluids in the subsurface.

Confinement Effects on Carbon Dioxide Methanation: A Novel Mechanism for Abiotic Methane Formation

An important scientific debate focuses on the possibility of abiotic synthesis of hydrocarbons during oceanic crust-seawater interactions. While on-site measurements near hydrothermal vents support this possibility, laboratory studies have provided data that are in some cases contradictory. At conditions relevant for sub-surface environments it has been shown that classic thermodynamics favour the production of CO2 from CH4, while abiotic methane synthesis would require the opposite. However, confinement effects are known to alter reaction equilibria. This report shows that indeed thermodynamic equilibrium can be shifted towards methane production, suggesting that thermal hydrocarbon synthesis near hydrothermal vents and deeper in the magma-hydrothermal system is possible. We report reactive ensemble Monte Carlo simulations for the CO2 methanation reaction. We compare the predicted equilibrium composition in the bulk gaseous phase to that expected in the presence of confinement. In the bulk phase we obtain excellent agreement with classic thermodynamic expectations. When the reactants can exchange between bulk and a confined phase our results show strong dependency of the reaction equilibrium conversions, [Formula: see text], on nanopore size, nanopore chemistry, and nanopore morphology. Some physical conditions that could shift significantly the equilibrium composition of the reactive system with respect to bulk observations are discussed.

Surface Structure and Stability of Partially Hydroxylated Silica Surfaces

Surface energies of silicates influence crack propagation during brittle fracture and decrease with surface relaxation caused by annealing and hydroxylation. Molecular-level simulations are particularly suited for the investigation of surface processes. In this work, classical MD simulations of silica surfaces are performed with two force fields (ClayFF and ReaxFF) to investigate the effect of force field reactivity on surface structure and energy as a function of surface hydroxylation. An unhydroxylated fracture surface energy of 5.1 J/m² is calculated with the ClayFF force field, and 2.0 J/m² is calculated for the ReaxFF force field. The ClayFF surface energies are consistent with the experimental results from double cantilever beam fracture tests (4.5 J/m²), whereas ReaxFF underestimated these surface energies. Surface relaxation via annealing and hydroxylation was performed by creating a low-energy equilibrium surface. Annealing condensed neighboring siloxane bonds increased the surface connectivity, and decreased the surface energies by 0.2 J/m² for ClayFF and 0.8 J/m² for ReaxFF. Posthydroxylation surface energies decreased further to 4.6 J/m² with the ClayFF force field and to 0.2 J/m² with the ReaxFF force field. Experimental equilibrium surface energies are ∼0.35 J/m², consistent with the ReaxFF force field. Although neither force field was capable of replicating both the fracture and equilibrium surface energies reported from experiment, each was consistent with one of these conditions. Therefore, future computational investigations that rely on accurate surface energy values should consider the surface state of the system and select the appropriate force field.

Molecular dynamics simulations of propane in slit shaped silica nano-pores: direct comparison with quasielastic neutron scattering experiments

MD simulations reveal the origin of anomalous pressure dependence of propane diffusion in silica mesopores.

Effects of Aspect Ratio on Water Immersion into Deep Silica Nanoholes

Understanding the influence of aspect ratio on water immersion into silica nanoholes is of significant importance to the etching process of semiconductor fabrication and other water immersion-related physical and biological processes. In this work, the processes of water immersion into silica nanoholes with different height/width aspect ratios (ϕ = 0.87, 1.92, 2.97, 4.01, 5.06) and different numbers of water molecules (N = 9986, 19972, 29958, 39944) were studied by molecular dynamics simulations. A comprehensive analysis has been conducted about the detailed process of water immersion and the influence of aspect ratios on water immersion rates. Five distinguishable stages were identified for the immersion process with all studied models. The results reveal that water can easily immerse into the silica nanoholes with larger ϕ and smaller N. The calculation also suggests that aspect ratios have a greater effect on water immersion rates for larger N numbers. The mechanism of the water immersion process is discussed in this work. We also propose a mathematical model to correlate the complete water immersion process for different aspect ratios.

Confined Water Determines Transport Properties of Guest Molecules in Narrow Pores

We computed the transport of methane through 1 nm wide slit-shaped pores carved out of solid substrates. Models for silica, magnesium oxide, and alumina were used as solid substrates. The pores were filled with water. The results show that the methane permeability through the hydrated pores is strongly dependent on the solid substrate. Detailed analysis of the simulated systems reveals that local properties of confined water, including its structure, and more importantly, evolution of solvation free energy and hydrogen bond structure are responsible for the pronounced differences observed. The simulations are extended to multicomponent systems representative of natural gas, containing methane, ethane, and H2S. The results show that all pores considered have high affinity for H2S, moderate affinity for methane, and low affinity for ethane. The H2S/methane transport selectivity through the hydrated alumina pore is comparable, or superior, to that reported for existing commercial membranes. A multiscale approach was then implemented to demonstrate that a Smoluchowski one-dimensional model is able to reproduce the molecular-level results for short pores when appropriate values for the local self-diffusion coefficients are used as input parameters. We propose that the model can be extended to predict methane transport through uniform hydrated pores of macroscopic length. When verified by experiments, our simulation results could have important implications in applications such as natural gas sweetening and predictions of methane migration through hydraulically fractured shale formations.

Factors governing the behaviour of aqueous methane in narrow pores

All-atom equilibrium molecular dynamics simulations were employed to investigate the behaviour of aqueous methane confined in 1-nm-wide pores obtained from different materials. Models for silica, alumina and magnesium oxide were used to construct the slit-shaped pores. The results show that methane solubility in confined water strongly depends on the confining material, with silica yielding the highest solubility in the systems considered here. The molecular structure of confined water differs within the three pores, and density fluctuations reveal that the silica pore is effectively less 'hydrophilic' than the other two pores considered. Comparing the water fluctuation autocorrelation function with local diffusion coefficients of methane across the hydrated pores we observed a direct proportional coupling between methane and water dynamics. These simulation results help to understand the behaviour of gas in water confined within narrow subsurface formations, with possible implications for fluid transport.

The importance of chemical potential in the determination of water slip in nanochannels

We investigate the slip properties of water confined in graphite-like nanochannels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle--but a different microscopic surface potential--do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.

Interfacial water properties in the presence of surfactants

Water, because of its fundamental role in biology, geology, and many industrial applications and its anomalous behavior compared to that of simple fluids, continues to fascinate and attract extensive scientific interest. Building on previous studies of water in contact with different surfaces, in this study, we report results obtained from molecular dynamics simulations of water near hydrophilic and hydrophobic interfaces in the presence of nonionic and ionic amphiphilic molecules, hexaethylene glycol monododecyl ether (C12E6) and sodium dodecyl sulfate (SDS). We elucidate how these surfactants affect the packing (i.e., density profiles) and orientation of interfacial water. The results highlight the interplay of both surfactant charges and the substrate charge distribution predominantly with respect to the orientation of water molecules, up to distances longer than those expected based on simulation results on flat solid surfaces. We also quantify the dynamics of interfacial water molecules by computing the residence probability for water in contact with various substrates. We compare our results to those previously obtained for interfacial water on silica and graphite and also with experimental sum-frequency vibrational spectroscopy results at the air-water interface in the presence of surfactants. Our analysis could be useful for a better understanding of interfacial water not only near solid substrates but also near self-assembled/aggregated molecules at a variety of interfaces.

Properties of water in the interfacial region of a polyelectrolyte bilayer adsorbed onto a substrate studied by computer simulations

We study the static and dynamic properties of water near a poly(styrene sulfonate)/poly(diallyldimethylammonium) (PSS/PDADMA) bilayer adsorbed onto a substrate by atomistic molecular dynamics simulations. Qualitative changes in the dynamics of water in the proximity of the adsorbed bilayer are observed - such as in the lateral diffusion, residence time and hydrogen-bonding lifetime - as compared with water in the presence of the bare substrate. Static properties of water are similarly influenced, and a high polarization of water molecules is found to be present surprisingly far from the adsorbed bilayer.

Aqueous NaCl and CsCl solutions confined in crystalline slit-shaped silica nanopores of varying degree of protonation

All-atom molecular dynamics simulations were conducted to study the dynamics of aqueous electrolyte solutions confined in slit-shaped silica nanopores of various degrees of protonation. Five degrees of protonation were prepared by randomly removing surface hydrogen atoms from fully protonated crystalline silica surfaces. Aqueous electrolyte solutions containing NaCl or CsCl salt were simulated at ambient conditions. In all cases, the ionic concentration was 1 M. The results were quantified in terms of atomic density distributions within the pores, and the self-diffusion coefficient along the direction parallel to the pore surface. We found evidence for ion-specific properties that depend on ion-surface, water-ion, and only in some cases ion-ion correlations. The degree of protonation strongly affects the structure, distribution, and the dynamic behavior of confined water and electrolytes. Cl(-) ions adsorb on the surface at large degrees of protonation, and their behavior does not depend significantly on the cation type (either Na(+) or Cs(+) ions are present in the systems considered). The cations show significant ion-specific behavior. Na(+) ions occupy different positions within the pore as the degree of protonation changes, while Cs(+) ions mainly remain near the pore center at all conditions considered. For a given degree of protonation, the planar self-diffusion coefficient of Cs(+) is always greater than that of Na(+) ions. The results are useful for better understanding transport under confinement, including brine behavior in the subsurface, with important applications such as environmental remediation.

Investigating the quartz (1010)/water interface using classical and ab initio molecular dynamics

Two different terminations of the (1010) surface of quartz (α and β) interacting with water are simulated by classical (CMD) (using two different force fields) and ab initio molecular dynamics (AIMD) and compared with previously published X-ray reflectivity (XR) experiments. Radial distribution functions between hydroxyl and water show good agreement between AIMD and CMD using the ClayFF force field for both terminations. The Lopes et al. (Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D. J. Phys. Chem. B2006, 110, 2782-2792) force field (LFF), however, underestimates the extent of hydroxyl-water hydrogen bonding. The β termination is found to contain hydroxyl-hydroxyl hydrogen bonds; the quartz surface hydroxyl hydrogens and oxygens that hydrogen bond with each other exhibit greatly reduced hydrogen bonding to water. Conversely, the hydroxyl hydrogen and oxygens that are not hydrogen bonded to other surface hydroxyls but are connected to those that are show a considerable amount of hydrogen bonding to water. The electron density distribution of an annealed surface of quartz (1010) obtained by XR is in qualitative agreement with electron densities calculated by CMD and AIMD. In all simulation methods, the interfacial water peak appears farther from the surface than observed by XR. Agreement among AIMD, LFF, and XR is observed for the relaxation of the near-surface atoms; however, ClayFF shows a larger discrepancy. Overall, results show that for both terminations of (1010), LFF treats the near-surface structure more accurately whereas ClayFF treats the interfacial water structure more accurately. It is shown that the number of hydroxyl and water hydrogen bonds to the bridging Si-O-Si oxygens connecting the surface silica groups to the rest of the crystal is much greater for the α than the β termination. It is suggested that this may play a role in the greater resistance to dissolution of the β termination than that of the α termination.