Ice tessellation on a hydroxylated silica surface

Molecular Dynamics Study on the Thermal Aspects of the Effect of Water Molecules at the Adhesive Interface on an Adhesive Structure

The presence of adsorbed water on hydrophilic solid surfaces should be taken into account, especially in humid environments. It significantly reduces the adhesive strength between the epoxy resin and the adherend surface. Here, the adhesion structure of interfacial water sandwiched between bisphenol A epoxy resin and a hydroxylated silica (001) surface is investigated with microsecond molecular dynamics simulations. Specifically, interfacial water layers with initial thicknesses of 7.5, 10, and 20 Å are modeled. The density curves of water and the diglycidyl ether of bisphenol A show that at room temperature, the surface of the silica with hydroxyl groups is completely covered with a thick layer of water. For water layers thinner than 10 Å, the density of epoxy resin on the silica surface increases when the system is heated and does not return to the original density when the system is cooled. Furthermore, calculation of the interaction energy revealed that the exclusion of water from the hydroxylated surface by epoxy resin during heating can contribute to the increase in the adhesive interaction between the epoxy resin and the silica surface with hydroxyl groups.

Structure of the Water Molecule Layer between Ice and Amorphous/Crystalline Surfaces Based on Molecular Dynamics Simulations

The structure of the water layer between the ice interface and the hydroxylated amorphous/crystalline silica surfaces was investigated using molecular dynamics simulations. The results indicate that the density profile in the direction perpendicular to the surface has two density peaks in the water layer at the ice-silica interface, which are affected by the silanol group density on the wall and the degree of supercooling in the system. In the two density peaks, the one facing the ice interface side has the same structure as the ice crystal, while the other density peak facing the silica surface has an icelike structure. In the solidification process, the ice and icelike structures in the layer progress more on the amorphous silica surface where the density of the silanol groups is low. The relationship between the ice crystallization and the thickness of the layer has been studied in detail; the lower the temperature, the more the ice crystallization progresses and the thinner the layer becomes.

Modification of adsorption, aggregation and wetting properties of surfactants by short chain alcohols

The adsorption of methanol, ethanol and propan-1-ol at the solution-air and solid-solution interfaces, their aggregation in the aqueous media as well as wetting properties regarding their applications as additives or co-surfactants in the surfactants aqueous solution were discussed based on the literature data. Mutual influence of alcohols and surfactants on the solution-air and solid-solution interface tension was considered. For this purpose there were used different methods allowing to describe or predict changes of water surface tension as a function of alcohols concentration. These, in turn, as a function of alcohol and/or surfactant concentration were also analyzed by means of the methods applied for prediction of surface tension of aqueous solution of the classical surfactants mixture. The same considerations related to the behaviour of alcohol and surfactant at the solid-solution and solution-air interfaces were made. To explain the behaviour of alcohols and surfactants mixture at the solution-air and solid-solution interfaces the components and parameters of water, alcohols, surfactants and solids surface tension as well as the Gibbs free energy changes during the adsorption process were taken into account. It was proved that wettability of some solids can be predicted based on alcohol and surfactants adsorption as well as surface tension components and parameters. As follows the mutual influence of alcohol and surfactant on their adsorption at the solution-air and solid-solution interfaces as well as on the wetting properties at the alcohol concentration from zero to its critical aggregation concentration (CAC) is different from that at its concentration higher than CAC.

Exchange Reactions at Mineral Interfaces

Exchange reactions are a family of chemical reactions that appear when mineral surfaces come into contact with protic solvents. Exchange reactions can also be understood as a unique interaction at mineral interfaces. Particularly significant interactions occurring at mineral surfaces are those with water and CO₂. The rather complex process occurring when minerals such as calcium silicate hydrate (C-S-H) phases come into contact with aqueous environments is referred to as a metal-proton exchange reaction (MPER). This process leads to the leaching of calcium ions from the near-surface region, the first step in the corrosion of cement-bound materials. Among the various corrosion reactions of C-S-H phases, the MPER appears to be the most important one. A promising approach to bridging certain problems caused by MPER and carbonation is the passivation of C-S-H surfaces. Today, such passivation is reached, for instance, by the functionalization of C-S-H surfaces with water-repelling organic films. Unfortunately, these organic films are weak against temperature and especially weak against abrasion. Exchange reactions at mineral interfaces allow the preparation of intrinsic, hydrophobic surfaces of C-S-H phases just at room temperature via a metal-metal exchange reaction.

Rich topologies of monolayer ices via unconventional electrowetting

Accurate manipulation of a substance on the nanoscale and ultimately down to the level of a single atom or molecule is an ongoing subject of frontier research. Herein, we show that topologies of water monolayers on substrates, in the complete wetting condition, can be manipulated into rich forms of ordered structures via electrowetting. Notably, two new topologies of monolayer ices were identified from our molecular dynamics simulations: one stable below room temperature and the other one having the ability to be stable at room temperature. Moreover, the wettability of the substrate can be tuned from superhydrophobic to superhydrophilic by uniformly changing the charge of each atomic site of the dipole or quadrupole distributed in an orderly manner on the model substrate. At a certain threshold value of the atomic charge, water droplets on the substrate can spread out spontaneously, achieving a complete electrowetting. Importantly, unlike the conventional electrowetting, which involves application of a uniform external electric field, we proposed non-conventional electrowetting, for the first time, by invoking the electric field of dipoles and quadrupoles embedded in the substrate. Moreover, different topologies of water monolayers can be achieved by using the non-conventional electrowetting. A major advantage of the non-conventional electrowetting is that the contact-angle saturation, a long-standing and known limitation in the field of electrowetting, can be overcome by tuning uniformly the lattice atomic charge at the surface, thereby offering a new way to mitigate the contact-angle saturation for various electrowetting applications.

Closest-Packing Water Monolayer Stably Intercalated in Phyllosilicate Minerals under High Pressure

The directional hydrogen-bond (HB) network and nondirectional van der Waals (vdW) interactions make up the specificity of water. Directional HBs could construct an ice-like monolayer in hydrophobic confinement even in the ambient regime. Here, we report a water monolayer dominated by vdW interactions confined in a phyllosilicate interlayer under high pressure. Surprisingly, it was in a thermodynamically stable state coupled with bulk water at the same pressure (P) and temperature (T), as revealed by the thermodynamic integration approach on the basis of molecular dynamics (MD) simulations. Both classical and ab initio MD simulations showed water O atoms were stably trapped and exhibited an ordered hexagonal closest-packing arrangement, but OH bonds of water reoriented frequently and exhibited a specific two-stage reorientation relaxation. Strikingly, hydration in the interlayer under high pressure had no relevance with surface hydrophilicity rationalized by the HB forming ability, which, however, determines wetting in the ambient regime. Intercalated water molecules were trapped by vdW interactions, which shaped the closest-packing arrangement and made hydration energetically available. The high pressure-volume term largely drives hydration, as it compensates the entropy penalty which is restricted by a relatively lower temperature. This vdW water monolayer should be ubiquitous in the high pressure but low-temperature regime.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

Thickness and Structure of Adsorbed Water Layer and Effects on Adhesion and Friction at Nanoasperity Contact

Most inorganic material surfaces exposed to ambient air can adsorb water, and hydrogen bonding interactions among adsorbed water molecules vary depending on, not only intrinsic properties of material surfaces, but also extrinsic working conditions. When dimensions of solid objects shrink to micro- and nano-scales, the ratio of surface area to volume increases greatly and the contribution of water condensation on interfacial forces, such as adhesion (Fa) and friction (Ft), becomes significant. This paper reviews the structural evolution of the adsorbed water layer on solid surfaces and its effect on Fa and Ft at nanoasperity contact for sphere-on-flat geometry. The details of the underlying mechanisms governing water adsorption behaviors vary depending on the atomic structure of the substrate, surface hydrophilicity and atmospheric conditions. The solid surfaces reviewed in this paper include metal/metallic oxides, silicon/silicon oxides, fluorides, and two-dimensional materials. The mechanism by which water condensation influences Fa is discussed based on the competition among capillary force, van der Waals force and the rupture force of solid-like water bridge. The condensed meniscus and the molecular configuration of the water bridge are influenced by surface roughness, surface hydrophilicity, temperature, sliding velocity, which in turn affect the kinetics of water condensation and interfacial Ft. Taking the effects of the thickness and structure of adsorbed water into account is important to obtain a full understanding of the interfacial forces at nanoasperity contact under ambient conditions.

Phase transition-like behavior of the water monolayer close to the polarized surface of a nanotube

By molecular dynamics simulations, we have investigated effects of temperature on the dynamical behavior of water layers at the charged surface of a nanotube. The behavior of the first water monolayer at the charged surface is very different from that of bulk water. There are three different temperature regions for the axial diffusion coefficient and they increase in different ways (linearly or exponentially) with temperature. The dipole distribution of water molecules was chosen as the order parameter to analyze the phase transition-like behavior. The simulation results indicate that the transition from ordered water to disordered water is continuous, which has not been found in the bulk counterpart. The mechanism behind the unexpected phenomenon was also investigated.

Phase behaviors of deeply supercooled bilayer water unseen in bulk water

Akin to bulk water, water confined to an isolated nanoslit can show a wealth of new 2D phases of ice and amorphous ice, as well as unusual phase behavior. Indeed, 2D water phases, such as bilayer hexagonal ice and monolayer square ice, have been detected in the laboratory, confirming earlier computational predictions. Herein, we report theoretical evidence of a hitherto unreported state, namely, bilayer very low density amorphous ice (BL-VLDA), as well as evidence of a strong first-order transition between BL-VLDA and the BL amorphous ice (BL-A), and a weak first-order transition between BL-VLDA and the BL very low density liquid (BL-VLDL) water. The diffusivity of BL-VLDA is typically in the range of 10-9 cm2/s to 10-10 cm2/s. Similar to bulk (3D) water, 2D water can exhibit two forms of liquid in the deeply supercooled state. However, unlike supercooled bulk water, for which the two forms of liquid can coexist and merge into one at a critical point, the 2D BL-VLDL and BL high-density liquid (BL-HDL) phases are separated by the highly stable solid phase of BL-A whose melting line exhibits the isochore end point (IEP) near 220 K in the temperature-pressure diagram. Above the IEP temperature, BL-VLDL and BL-HDL are indistinguishable. At negative pressures, the metastable BL-VLDL exhibits a spatially and temporally heterogeneous structure induced by dynamic changes in the nanodomains, a feature much less pronounced in the BL-HDL.

Thermodynamic, structural, and dynamical properties of nano-confined water using SPC/E and TIP4P models by molecular dynamics simulations

Different morphologies of water molecules are confined between two parallel graphene surfaces.

Self-assembly of glycine on Cu(001): the effect of temperature and polarity

Glycine on Cu(001) is studied as an example to illustrate the critical role of finite temperature and molecular polarity in the self-assembly of biomolecules at a metal surface.

Two Dimensional Ice from First Principles: Structures and Phase Transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2  GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Structural and configurational properties of nanoconfined monolayer ice from first principles

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

Computational Free Energy Studies of a New Ice Polymorph Which Exhibits Greater Stability than Ice Ih

The absolute free energies of several ice polymorphs were calculated using thermodynamic integration. These polymorphs are predicted by computer simulations using a variety of common water models to be stable at low pressures. A recently discovered ice polymorph that has as yet only been observed in computer simulations (Ice-i) was determined to be the stable crystalline state for all the water models investigated. Phase diagrams were generated, and phase coexistence lines were determined for all of the known low-pressure ice structures. Additionally, potential truncation was shown to play a role in the resulting shape of the free energy landscape.

Liquid-solid and solid-solid phase transition of monolayer water: high-density rhombic monolayer ice

Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.

Highly confined water: two-dimensional ice, amorphous ice, and clathrate hydrates

Understanding phase behavior of highly confined water, ice, amorphous ice, and clathrate hydrates (or gas hydrates), not only enriches our view of phase transitions and structures of quasi-two-dimensional (Q2D) solids not seen in the bulk phases but also has important implications for diverse phenomena at the intersection between physical chemistry, cell biology, chemical engineering, and nanoscience. Relevant examples include, among others, boundary lubrication in nanofluidic and lab-on-a-chip devices, synthesis of antifreeze proteins for ice-growth inhibition, rapid cooling of biological suspensions or quenching emulsified water under high pressure, and storage of H2 and CO2 in gas hydrates. Classical molecular simulation (MD) is an indispensable tool to explore states and properties of highly confined water and ice. It also has the advantage of precisely monitoring the time and spatial domains in the sub-picosecond and sub-nanometer scales, which are difficult to control in laboratory experiments, and yet allows relatively long simulation at the 10(2) ns time scale that is impractical with ab initio molecular dynamics simulations. In this Account, we present an overview of our MD simulation studies of the structures and phase behaviors of highly confined water, ice, amorphous ice, and clathrate, in slit graphene nanopores. We survey six crystalline phases of monolayer (ML) ice revealed from MD simulations, including one low-density, one mid-density, and four high-density ML ices. We show additional supporting evidence on the structural stabilities of the four high-density ML ices in the vacuum (without the graphene confinement), for the first time, through quantum density-functional theory optimization of their free-standing structures at zero temperature. In addition, we summarize various low-density, high-density, and very-high-density Q2D bilayer (BL) ice and amorphous ice structures revealed from MD simulations. These simulations reinforce the notion that the nanoscale confinement not only can disrupt the hydrogen bonding network in bulk water but also can allow satisfaction of the ice rule for low-density and high-density Q2D crystalline structures. Highly confined water can serve as a generic model system for understanding a variety of Q2D materials science phenomena, for example, liquid-solid, solid-solid, solid-amorphous, and amorphous-amorphous transitions in real time, as well as the Ostwald staging during these transitions. Our simulations also bring new molecular insights into the formation of gas hydrate from a gas and water mixture at low temperature.

Two-dimensional water diffusion at a graphene-silica interface

Because of the dominant role of the surface of molecules and their individuality, molecules behave distinctively in a confined space, which has far-reaching implications in many physical, chemical, and biological systems. Here we demonstrate that graphene forms a unique atom-thick interstitial space that enables the study of molecular diffusion in two dimensions with underlying silica substrates. Raman spectroscopy visualized intercalation of water from the edge to the center underneath graphene in real time, which was dictated by the hydrophilicity of the substrates. In addition, graphene undergoes reversible deformation to conform to intercalating water clusters or islands. Atomic force microscopy confirmed that the interfacial water layer is only ca. 3.5 Å thick, corresponding to one bilayer unit of normal ice. This study also demonstrates that oxygen species responsible for the ubiquitous hole doping are located below graphene. In addition to serving as a transparent confining wall, graphene and possibly other two-dimensional materials can be used as an optical indicator sensitive to interfacial mass transport and charge transfer.

Ferroelectric hexagonal and rhombic monolayer ice phases

Two new monolayer ice phases are predicted from molecular dynamics simulations, both proven to be ferroelectric.

Absorption cross sections of surface-adsorbed H2O in the 295-370 nm region and heterogeneous nucleation of H2O on fused silica surfaces

We have determined absorption cross sections of a monolayer of H2O adsorbed on the fused silica surfaces in the 295-370 nm region at 293 ± 1 K by using Brewster angle cavity ring-down spectroscopy. Absorption cross sections of surface-adsorbed H2O vary between (4.66 ± 0.83) × 10(-20) and (1.73 ± 0.52) × 10(-21) cm(2)/molecule over this wavelength range, where errors quoted represent experimental scatter (1σ). Our experimental study provides direct evidence that surface-adsorbed H2O is an absorber of the near UV solar radiation. We also varied the H2O pressure in the surface study cell over the 0.01-17 Torr range and obtained probe laser absorptions at 295, 340, and 350 nm by multilayer of adsorbed H2O molecules until the heterogeneous nucleation of water occurred on fused silica surfaces. The average absorption cross sections of multilayer adsorbed H2O are (2.17 ± 0.53) × 10(-20), (2.48 ± 0.67) × 10(-21), and (2.34 ± 0.59) × 10(-21) cm(2)/molecule at 295, 340, and 350 nm. The average absorption cross sections of transitional H2O layer are (6.06 ± 2.73) × 10(-20), (6.48 ± 3.85) × 10(-21), and (8.04 ± 4.92) × 10(-21) cm(2)/molecule at 295, 340, and 350 nm. The average thin water film absorption cross sections are (2.39 ± 0.50) × 10(-19), (3.21 ± 0.81) × 10(-20), and (3.37 ± 0.94) × 10(-20) cm(2)/molecule at 295 nm, 340 nm, and 350 nm. Atmospheric implications of the results are discussed.

Polymorphism and polyamorphism in bilayer water confined to slit nanopore under high pressure

A distinctive physical property of bulk water is its rich solid-state phase behavior, which includes 15 crystalline (ice I-ice XIV) and at least 3 glassy forms of water, namely, low-density amorphous, high-density amorphous, and very-high-density amorphous (VHDA). Nanoscale confinement adds a new physical variable that can result in a wealth of new quasi-2D phases of ice and amorphous ice. Previous computer simulations have revealed that when water is confined between two flat hydrophobic plates about 7-9 Å apart, numerous bilayer (BL) ices (or polymorphs) can arise [e.g., BL-hexagonal ice (BL-ice I)]. Indeed, growth of the BL-ice I through vapor deposition on graphene/Pt(111) substrate has been achieved experimentally. Herein, we report computer simulation evidence of pressure-induced amorphization from BL-ice I to BL-amorphous and then to BL-VHDA(2) at 250 K and 3 GPa. In particular, BL-VHDA(2) can transform into BL-VHDA(1) via decompression from 3 to 1.5 GPa at 250 K. This phenomenon of 2D polyamorphic transition is akin to the pressure-induced amorphization in 3D ice (e.g., from hexagonal ice to HDA and then to VHDA via isobaric annealing). Moreover, when the BL-ice I is compressed instantly to 6 GPa, a new very-high-density BL ice is formed. This new phase of BL ice can be viewed as an array of square ice nanotubes. Insights obtained from pressure-induced amorphization and crystallization of confined water offer a guide with which to seek a thermodynamic path to grow a new form of methane clathrate whose BL ice framework exhibits the Archimedean 4⋅8(2) (square-octagon) pattern.

Simulation of Forces between Humid Amorphous Silica Surfaces: A Comparison of Empirical Atomistic Force Fields

Atmospheric humidity strongly influences the interactions between dry granular particles in process containers. To reduce the energy loss in industrial production processes caused by particle agglomeration, a basic understanding of the dependence of particle interactions on humidity is necessary. Hence, in this study, molecular dynamic simulations were carried out to calculate the adhesion between silica surfaces in the presence of adsorbed water. For a realistic description, the choice of force field is crucial. Because of their frequent use and transferability to biochemical systems, the Clay and CWCA force fields were investigated with respect to their ability to describe the water-silica interface in comparison to the more advanced Reax force field, ab initio calculations, and experiments.

Constrained evolutionary algorithm for structure prediction of molecular crystals: methodology and applications

Evolutionary crystal structure prediction proved to be a powerful approach for studying a wide range of materials. Here we present a specifically designed algorithm for the prediction of the structure of complex crystals consisting of well defined molecular units. The main feature of this new approach is that each unit is treated as a whole body, which drastically reduces the search space and improves the efficiency, but necessitates the introduction of new variation operators described here. To increase the diversity of the population of structures, the initial population and part (~20%) of the new generations are produced using space-group symmetry combined with random cell parameters, and random positions and orientations of molecular units. We illustrate the efficiency and reliability of this approach by a number of tests (ice, ammonia, carbon dioxide, methane, benzene, glycine and butane-1,4-diammonium dibromide). This approach easily predicts the crystal structure of methane A containing 21 methane molecules (105 atoms) per unit cell. We demonstrate that this new approach also has a high potential for the study of complex inorganic crystals as shown on examples of a complex hydrogen storage material Mg(BH(4))(2) and elemental boron.

Effect of anionic surfactant and short-chain alcohol mixtures on adsorption at quartz/water and water/air interfaces and the wettability of quartz

Measurements of the advancing contact angles for aqueous solutions of sodium dodecyl sulfate (SDDS) or sodium hexadecyl sulfonate (SHS) in mixtures with methanol, ethanol, or propanol on a quartz surface were carried out. On the basis of the obtained results and Young and Gibbs equations the critical surface tension of quartz wetting, the composition of the surface layer at the quartz-water interface, and the activity coefficients of the anionic surfactants and alcohols in this layer as well as the work of adhesion of aqueous solutions of anionic surfactant and alcohol mixtures to the quartz surface were determined. The analysis of the contact angle data showed that the wettability of quartz changed visibly only in the range of alcohol and anionic surfactant concentration at which these surface-active agents were present in the solution in the monomeric form. The analysis also showed that there was a linear dependence between the adhesion and the surface tension of aqueous solutions of anionic surfactant and alcohol mixtures. This dependence can be described by linear equations for which the constants depend on the anionic surfactant and alcohol concentrations. The slope of all linear dependence between adhesion and surface tension was positive. The critical surface tension of quartz wetting determined from this dependence by extrapolating the adhesion tension to the value equal to the surface tension (for contact angle equal zero) depends on the assumption whether the concentration of anionic surfactant or alcohol was constant. Its average value is equal to 29.95mN/m and it is considerably lower than the quartz surface tension. The positive slope of the adhesion-surface tension curves was explained by the possibility of the presence of liquid vapor film beyond the solution drop which settled on the quartz surface and the adsorption of surface-active agents at the quartz/monolayer water film-water interface. This conclusion was confirmed by the work of adhesion of aqueous solutions of anionic surfactants and short-chain alcohol mixtures to the quartz surface determined on the basis of the contact angle data and molar fraction of anionic surfactants and alcohols and their activity coefficient in the surface layer.

Experimental and theoretical characterization of adsorbed water on self-assembled monolayers: understanding the interaction of water with atmospherically relevant surfaces

A combination of experiments and molecular dynamic (MD) simulations has been applied to elucidate the nature of water on organic self-assembled monolayers (SAMs) before and after oxidation. SAMs mimic organics adsorbed on environmental urban surfaces. Water on clean or SAM-coated borosilicate glass surfaces was measured at equilibrium as a function of relative humidity (RH), using transmission Fourier transform infrared (FTIR) spectroscopy at 1 atm and 22 +/- 1 degrees C. The SAMs included C18 and C8 alkanes, as well as the C8 terminal alkene. Oxidation of the terminal alkene SAM was carried out with either KMnO(4) solution or gaseous O(3). The FTIR data showed at least two distinct peaks due to water on these surfaces, one at approximately 3200 cm(-1), which dominates at low RH (20%), and one at approximately 3400 cm(-1) at high RH (80%), which is similar to that in bulk liquid water. Temperature-programmed desorption (TPD) experiments showed that oxidation leads to more strongly adsorbed water. However, the amount of water in equilibrium with water vapor on the oxidized alkene was not significantly different from that on the unoxidized SAM, although there was a change in the relative intensities of the two contributing infrared peaks at 80% RH. MD simulations with hydrogen bond analysis suggest that molecules on the surface of small water clusters that dominate on SAM surfaces at low RH have fewer hydrogen bonds, while those in the interior of the clusters have three and four hydrogen bonds similar to bulk liquid water. Taken together, the experimental infrared data and MD simulations suggest a correlation between the relative intensities of the 3200 cm(-1)/3400 cm(-1) bands and the hydrogen-bonding patterns of the water on the surface and in the interior of clusters on the SAM surfaces. These studies suggest that water clusters will be present even on hydrophobic surfaces in the atmosphere and hence are available to participate in heterogeneous chemistry. In addition, oxidation of organic coatings on atmospheric particles or surfaces in the boundary layer may not lead to enhanced water uptake as is often assumed.

Body-centered tetragonal C4: a viable sp3 carbon allotrope

We have investigated by first principles the electronic, vibrational, and structural properties of bct C4, a new form of crystalline sp{3} carbon recently found in molecular dynamics simulations of carbon nanotubes under pressure. This phase is transparent, dynamically stable at zero pressure, and more stable than graphite beyond 18.6 GPa. Coexistence of bct C4 with M carbon can explain better the x-ray diffraction pattern of a transparent and hard phase of carbon produced by the cold compression of graphite. Its structure appears to be intermediate between that of graphite and hexagonal diamond. These facts suggest that bct C4 is an accessible form of sp{3} carbon along the graphite-to-hexagonal diamond transformation path.

Guest-free monolayer clathrate and its coexistence with two-dimensional high-density ice

Three-dimensional (3D) gas clathrates are ice-like but distinguished from bulk ices by containing polyhedral nano-cages to accommodate small gas molecules. Without space filling by gas molecules, standalone 3D clathrates have not been observed to form in the laboratory, and they appear to be unstable except at negative pressure. Thus far, experimental evidence for guest-free clathrates has only been found in germanium and silicon, although guest-free hydrate clathrates have been found, in recent simulations, able to grow from cold stretched water, if first nucleated. Herein, we report simulation evidence of spontaneous formation of monolayer clathrate ice, with or without gas molecules, within hydrophobic nano-slit at low temperatures. The guest-free monolayer clathrate ice is a low-density ice (LDI) whose geometric pattern is identical to Archimedean 4.8(2)-truncated square tiling, i.e. a mosaic of tetragons and octagons. At large positive pressure, a second phase of 2D monolayer ice, i.e. the puckered square high-density ice (HDI) can form. The triple point of the LDI/liquid/HDI three-phase coexistence resembles that of the ice-I(h)/water/ice-III three-phase coexistence. More interestingly, when the LDI is under a strong compression at 200 K, it transforms into the HDI via a liquid intermediate state, the first direct evidence of Ostwald's rule of stages at 2D. The tensile limit of the 2D LDI and water are close to that of bulk ice-I(h) and laboratory water.

Molecular dynamics studies of the interactions of water and amino acid analogues with quartz surfaces

The interactions of silica surfaces with water and biomolecules are of considerable significance in bio- and nanotechnology and in geochemistry. An important goal in the fields of biomineralization and biomimetics is to fine-tune these interactions for the control, e.g., of assembly of materials at the nanoscale. Here we report molecular dynamics simulations of fully hydroxylated alpha-quartz (1010), (0001), and (0111) surfaces in explicit water. We also present free energy estimates of adsorbing water and analogues of amino acid side chains onto the quartz (1010) surface. We find that at least two layers of structured water form on the surface, which is driven by the formation of a strong hydrogen bond network at the interface. Interestingly, we find that the free energy change to move methane (analogue of the side chain of alanine) from bulk water to the (1010) interface is favorable. We ascribe this to the presence of microscopic voids on the surface, which can accommodate small hydrophobic moieties and shield them from the solvent. These observations draw some useful insights into the possible mechanisms by which biomolecules, in particular peptides and proteins, bind to quartz and other silica surfaces.

Growth and structure of water on SiO2 films on Si Investigated by Kelvin probe microscopy and in Situ X-ray spectroscopies

The growth of water on thin SiO2 films on Si wafers at vapor pressures between 1.5 and 4 Torr and temperatures between -10 and 21 degrees C has been studied in situ using Kelvin probe microscopy and X-ray photoemission and absorption spectroscopies. From 0 to 75% relative humidity (RH), water adsorbs forming a uniform film 4-5 layers thick. The surface potential increases in that RH range by about 400 mV and remains constant upon further increase of the RH. Above 75% RH, the water film grows rapidly, reaching 6-7 monolayers at around 90% RH and forming a macroscopic drop near 100%. The O K-edge near-edge X-ray absorption spectrum around 75% RH is similar to that of liquid water (imperfect H-bonding coordination) at temperatures above 0 degrees C and is ice-like below 0 degrees C.

Stability of two-dimensional tessellation ice on the hydroxylated beta-cristobalite (100) surface

Monolayer adsorbed water on the beta-cristobalite (100) surface is studied via classical molecular dynamics simulations. The ordered two-dimensional (2D) tessellation ice structure (i.e., the four-membered and the eight-membered rings appear alternatively) is justified at low temperatures in the simulations. The stability of this possible new ice phase is further investigated by heating the system from 5 to 300 K. An order-disorder structural transition is observed between 100 and 200 K, featuring the melting process of the tessellation ice. This process is characterized by the water oxygen-oxygen radial distribution function, the coordination number, the distance vector between the center of mass of the oxygen and the hydrogen atoms in water, the mean square displacement of oxygen in water, and the vibrational density of state. The above techniques show consistency on that the order-disorder transition temperature of the 2D tessellation ice is far below 300 K. The 2D tessellation ice structure is also obtained via density functional calculations with different generalized gradient approximations. By comparing the calculated adsorption and the lateral energies between different methods, we find that the melting temperature of the specific 2D ice structure is strongly method dependent. Therefore, further experimental works are urged to justify this possible new ice phase and probe its stability.