Microscopic wetting of mixed self-assembled monolayers: a molecular dynamics study

Electron microscopy approach to the wetting dynamics of single organosilanized mesopores

Columnar mesoporous silicon (PSi) with hydrophobic vs. hydrophilic chemistries was chosen as a model for the local (pore-by-pore) study of water-pore interactions. Tomographic reconstructions provided a 3D view of the ramified pore structure. An in situ study of PSi wetting was conducted for categorized pore diameters by environmental scanning TEM. An appropriate setting of the contrast allows for the normalization of the gray scale in the images as a function of relative humidity (RH). This allows constructing an isotherm for each single pore and a subsequent averaging provides an isotherm for each pore size range. The isotherms systematically point to an initial adsorption through the formation of water adlayers, followed by a capillary filling process at higher RH. The local isotherms correlate with (global) gravimetric determination of wetting. Our results point at the validation of a technique for the study of aging and stability of single-pore nanoscale devices.

Sequence Modulates Polypeptoid Hydration Water Structure and Dynamics

We use molecular dynamics simulations to investigate the effect of polypeptoid sequence on the structure and dynamics of its hydration waters. Polypeptoids provide an excellent platform to study small-molecule hydration in disordered polymers, as they can be precisely synthesized with a variety of sidechain chemistries. We examine water behavior near a set of peptoid oligomers in which the number and placement of nonpolar versus polar sidechains are systematically varied. To do this, we leverage a new computational workflow enabling accurate sampling of polypeptoid conformations. We find that the hydration waters are less dense, are more tetrahedral, and have slower dynamics compared to bulk water. The magnitude of these shifts increases with the number of nonpolar groups. We also find that shifts in the water structure and dynamics are strongly correlated, suggesting that experimental insight into the dynamics of hydration water obtained by Overhauser dynamic nuclear polarization (ODNP) also contains information about water structural properties. We then demonstrate the ability of ODNP to probe site-specific dynamics of hydration water near these model peptoid systems.

Molecular Dynamics Simulations of Water Condensation on Surfaces with Tunable Wettability

Water condensation plays a major role in a wide range of industrial applications. Over the past few years, many studies have shown interest in designing surfaces with enhanced water condensation and removal properties. It is well known that heterogeneous nucleation outperforms homogeneous nucleation in the condensation process. Because heterogeneous nucleation initiates on a surface at a small scale, it is highly desirable to characterize water-surface interactions at the molecular level. Molecular dynamics (MD) simulations can provide direct insight into heterogeneous nucleation and advance surface designs. Existing MD simulations of water condensation on surfaces were conducted by tuning the solid-water van der Waals interaction energy as a substitute for modeling surfaces with different wettabilities. However, this approach cannot reflect the real intermolecular interactions between the surface and water molecules. Here, we report MD simulations of water condensation on realistic surfaces of alkanethiol self-assembled monolayers with different head group chemistries. We show that decreasing surface hydrophobicity significantly increases the electrostatic forces between water molecules and the surface, thus increasing the water condensation rate. We observe a strong correlation between our rate of condensation results and the results from other surface characterization metrics, such as the interfacial thermal conductance, contact angle, and the molecular-scale wettability metric of Garde and co-workers. This work provides insight into the water condensation process at the molecular scale on surfaces with tunable wettability.

Hydrophilic directional slippery rough surfaces for water harvesting

Multifunctional surfaces that are favorable for both droplet nucleation and removal are highly desirable for water harvesting applications but are rare. Inspired by the unique functions of pitcher plants and rice leaves, we present a hydrophilic directional slippery rough surface (SRS) that is capable of rapidly nucleating and removing water droplets. Our surfaces consist of nanotextured directional microgrooves in which the nanotextures alone are infused with hydrophilic liquid lubricant. We have shown through molecular dynamics simulations that the physical origin of the efficient droplet nucleation is attributed to the hydrophilic surface functional groups, whereas the rapid droplet removal is due to the significantly reduced droplet pinning of the directional surface structures and slippery interface. We have further demonstrated that the SRS, owing to its large surface area, hydrophilic slippery interface, and directional liquid repellency, outperforms conventional liquid-repellent surfaces in water harvesting applications.

Molecular dynamics simulations of nanoscale engravings on an alkanethiol monolayer

Thermal stability of nanoscale engravings on alkanethiol monolayer.

Surface fractal dimension, water adsorption efficiency, and cloud nucleation activity of insoluble aerosol

Surface porosity affects the ability of a substance to adsorb gases. The surface fractal dimension D is a measure that indicates the amount that a surface fills a space, and can thereby be used to characterize the surface porosity. Here we propose a new method for determining D, based on measuring both the water vapour adsorption isotherm of a given substance, and its ability to act as a cloud condensation nucleus when introduced to humidified air in aerosol form. We show that our method agrees well with previous methods based on measurement of nitrogen adsorption. Besides proving the usefulness of the new method for general surface characterization of materials, our results show that the surface fractal dimension is an important determinant in cloud drop formation on water insoluble particles. We suggest that a closure can be obtained between experimental critical supersaturation for cloud drop activation and that calculated based on water adsorption data, if the latter is corrected using the surface fractal dimension of the insoluble cloud nucleus.

Effects of surface pressure on the properties of Langmuir monolayers and interfacial water at the air-water interface

The effects of surface pressure on the physical properties of Langmuir monolayers of palmitic acid (PA) and dipalmitoylphosphatidic acid (DPPA) at the air/water interface are investigated through molecular dynamics simulations with atomistic force fields. The structure and dynamics of both monolayers and interfacial water are compared across the range of surface pressures at which stable monolayers can form. For PA monolayers at T = 300 K, the untilted condensed phase with a hexagonal lattice structure is found at high surface pressure, while the uniformly tilted condensed phase with a centered rectangular lattice structure is observed at low surface pressure, in agreement with the available experimental data. A state with uniform chain tilt but no periodic spatial ordering is observed for DPPA monolayers on a Na(+)/water subphase at both high and low surface pressures. The hydrophobic acyl chains of both monolayers pack efficiently at all surface pressures, resulting in a very small number of gauche defects. The analysis of the hydrogen-bonding structure/dynamics at the monolayer/water interface indicates that water molecules hydrogen-bonded to the DPPA head groups reorient more slowly than those hydrogen-bonded to the PA head groups, with the orientational dynamics becoming significantly slower at high surface pressure. Possible implications for physicochemical processes taking place on marine aerosols in the atmosphere are discussed.

Molecular insight into the dynamical adsorption behavior of nanoscale water droplets on a heterogeneous surface

A heterogeneous surface is constructed by adding one hydrophilic patch at the center of a hydrophobic surface, and the dynamical adsorption process of nanoscale water droplets is investigated adopting molecular dynamics simulations.

Molecular insight into nanoscale water films dewetting on modified silica surfaces

In this work, molecular dynamics simulations are adopted to investigate the microscopic dewetting mechanism of nanoscale water films on methylated silica surfaces.

Molecular mechanism of stabilization of thin films for improved water evaporation protection

All-atom molecular dynamics simulations and experimental characterization have been used to examine the structure and dynamics of novel evaporation-suppressing films where the addition of a water-soluble polymer to an ethylene glycol monooctadecyl ether monolayer leads to improved water evaporation resistance. Simulations and Langmuir trough experiments demonstrate the surface activity of poly(vinyl pyrrolidone) (PVP). Subsequent MD simulations performed on the thin films supported by the PVP sublayer show that, at low surface pressures, the polymer tends to concentrate at the film/water interface. The simulated atomic concentration profiles, hydrogen bonding patterns, and mobility analyses of the water-polymer-monolayer interfaces reveal that the presence of PVP increases the atomic density near the monolayer film, improves the film stability, and reduces the mobility of interfacial waters. These observations explain the molecular basis of the improved efficacy of these monolayer/polymer systems for evaporation protection of water and can be used to guide future development of organic thin films for other applications.

Computer simulation of adhesion between hydrophilic and hydrophobic self-assembled monolayers in water

The grand canonical Monte Carlo technique and atomistic force fields are used to calculate the force-distance relations and free energies of adhesion between carboxyl and methyl terminated alkanethiolate self-assembled monolayers (SAMs) in water. Both symmetric and asymmetric confinements are considered, as formed by like and unlike SAMs, respectively. As the confinement is increased, water confined by the hydrophobic methyl terminated SAMs experiences capillary evaporation. As a consequence, the adhesion energy is determined by the direct interaction between bare SAMs. In the asymmetric system, an incomplete capillary evaporation is observed, with the number of water molecules dropped by more than an order of magnitude. The remaining water molecules are all adsorbed on the hydrophilic SAM, while the hydrophobic SAM is separated from the rest of the system by a thin vapor layer. The calculated free energies of adhesion are in acceptable agreement with experiment.

Water adsorption isotherms on CH3-, OH-, and COOH-terminated organic surfaces at ambient conditions measured with PM-RAIRS

The water adsorption isotherms on methyl (CH(3))-, hydroxyl (OH)-, and carboxylic acid (COOH)-terminated alkylthiol self-assembled monolayers (SAMs) on Au were studied at room temperature and ambient pressure with polarization modulation reflection-absorption infrared spectroscopy (PM-RAIRS). PM-RAIRS analysis showed that water does not adsorb at all on the CH(3)-SAM/Au at subsaturation humidity conditions. In a dry Ar environment, the OH-SAM/Au holds at least 2 layer thick strongly bound water molecules which exhibit a broad O-H stretch vibration peak centered at ∼3360 cm(-1). The peak position implies that the strongly bound water layer on the OH SAM is more like a liquid than an ice. The additional uptake of water in humid environments is relatively weak, and the peak position changes very little. Unlike the OH-SAM/Au, the COOH-SAM/Au does not have strongly bound water layer. This seems to be due to the strong hydrogen bonding between terminal COOH groups in dry conditions. The weak interactions between water and carboxyl groups at low relative humidity (RH) and the solvation of dissociated carboxylic groups in high RH lead to a type III isotherm behavior, based on the BET categories, for water adsorption on the COOH-SAM/Au. The water spectra on the COOH-SAM at RH > 45% are centered at ∼3430 cm(-1) and very broad, indicating that the hydrogen-bonding network of water on the COOH-SAM is much different from that on the OH-SAM.

Molecular dynamics simulations of low-ordered alzheimer β-amyloid oligomers from dimer to hexamer on self-assembled monolayers

Accumulation of small soluble oligomers of amyloid-β (Aβ) in the human brain is thought to play an important pathological role in Alzheimer's disease. The interaction of these Aβ oligomers with cell membrane and other artificial surfaces is important for the understanding of Aβ aggregation and toxicity mechanisms. Here, we present a series of exploratory molecular dynamics (MD) simulations to study the early adsorption and conformational change of Aβ oligomers from dimer to hexamer on three different self-assembled monolayers (SAMs) terminated with CH(3), OH, and COOH groups. Within the time scale of MD simulations, the conformation, orientation, and adsorption of Aβ oligomers on the SAMs is determined by complex interplay among the size of Aβ oligomers, the surface chemistry of the SAMs, and the structure and dynamics of interfacial waters. Energetic analysis of Aβ adsorption on the SAMs reveals that Aβ adsorption on the SAMs is a net outcome of different competitions between dominant hydrophobic Aβ-CH(3)-SAM interactions and weak CH(3)-SAM-water interactions, between dominant electrostatic Aβ-COOH-SAM interactions and strong COOH-SAM-water interactions, and between comparable hydrophobic and electrostatic Aβ-OH-SAM interactions and strong OH-SAM-water interactions. Atomic force microscopy images also confirm that all of three SAMs can induce the adsorption and polymerization of Aβ oligomers. Structural analysis of Aβ oligomers on the SAMs shows a dramatic increase in structural stability and β-sheet content from dimer to trimer, suggesting that Aβ trimer could act as seeds for Aβ polymerization on the SAMs. This work provides atomic-level understanding of Aβ peptides at interface.

Nanoscale water condensation on click-functionalized self-assembled monolayers

We have examined the nanoscale adsorption of molecular water under ambient conditions onto a series of well-characterized functionalized surfaces produced by Cu(I)-catalyzed alkyne-azide cycloaddition (CuAAC or "click") reactions on alkyne-terminated self-assembled monolayers on silicon. Water contact angle (CA) measurements reveal a range of macroscopic hydrophilicity that does not correlate with the tendency of these surfaces to adsorb water at the molecular level. X-ray reflectometry has been used to follow the kinetics of water adsorption on these "click"-functionalized surfaces, and also shows that dense continuous molecular water layers are formed over 30 h. For example, a highly hydrophilic surface, functionalized by an oligo(ethylene glycol) moiety (with a CA = 34°) showed 2.9 Å of adsorbed water after 30 h, while the almost hydrophobic underlying alkyne-terminated monolayer (CA = 84°) showed 5.6 Å of adsorbed water over the same period. While this study highlights the capacity of X-ray reflectometry to study the structure of adsorbed water on these surfaces, it should also serve as a warning for those intending to characterize self-assembled monolayers and functionalized surfaces to avoid contamination by even trace amounts of water vapor. Moreover, contact angle measurements alone cannot be relied upon to predict the likely degree of moisture uptake on such surfaces.

Microscopic wetting of self-assembled monolayers with different surfaces: a combined molecular dynamics and quantum mechanics study

Molecular dynamics simulations are used to study the micronature of the organization of water molecules on the flat surface of well-ordered self-assembled monolayers (SAMs) of 18-carbon alkanethiolate chains bound to a silicon (111) substrate. Six different headgroups (-CH(3), -C═C, -OCH(3), -CN, -NH(2), -COOH) are used to tune the character of the surface from hydrophobic to hydrophilic, while the level of hydration is consistent on all six SAM surfaces. Quantum mechanics calculations are employed to optimize each alkyl chain of the different SAMs with one water molecule and to investigate changes in the configuration of each headgroup under hydration. We report the changes of the structure of the six SAMs with different surfaces in the presence of water, and the area of the wetted surface of each SAM, depending on the terminal group. Our results suggest that a corrugated and hydrophobic surface will be formed if the headgroups of SAM surface are not able to form H-bonds either with water molecules or between adjacent groups. In contrast, the formation of hydrogen bonds not only among polar heads but also between polar heads and water may enhance the SAM surface hydrophilicity and corrugation. We explicitly discuss the micromechanisms for the hydration of three hydrophilic SAM (CN-, NH(2)- and COOH-terminated) surfaces, which is helpful to superhydrophilic surface design of SAM in biomimetic materials.

Amphiphilic silica nanoparticles at the decane-water interface: insights from atomistic simulations

The properties of 3 nm-diameter silica nanoparticles with different surface chemistry were systematically investigated at the decane-water interface using molecular dynamics simulations. Our results show that the decane-water interfacial tension is not much influenced by the presence of the nanoparticles. The three-phase contact angle increases with nanoparticle surface hydrophobicity. Contact angles observed for the nanoparticles at 300 and at 350 K differ very little. The contact angle of the nanoparticle with randomly dispersed hydrophobic groups is smaller than that observed in Janus nanoparticles of equal overall surface chemistry composition. The energy necessary to desorb Janus nanoparticles from the interface is usually higher than that required to desorb the corresponding homogeneous nanoparticles. Desorption from the interface into the aqueous phase is preferred over that into the organic phase for all except one of the nanoparticles considered. Structural and dynamic properties including nanoparticle rotational relaxation, solvent density profiles, and solvent residence autocorrelation functions near the nanoparticles are also presented. The data are useful for designing Pickering emulsions.

Probing surface tension additivity on chemically heterogeneous surfaces by a molecular approach

Surface free energy of a chemically heterogeneous surface is often treated as an approximately additive quantity through the Cassie equation [Cassie ABD (1948) Discuss Faraday Soc 3:11-16]. However, deviations from additivity are common, and molecular interpretations are still lacking. We use molecular simulations to measure the microscopic analogue of contact angle, Θ(c), of aqueous nanodrops on heterogeneous synthetic and natural surfaces as a function of surface composition. The synthetic surfaces are layers of graphene functionalized with prototypical nonpolar and polar head group: methyl, amino, and nitrile. We demonstrate positive as well as negative deviations from the linear additivity. We show the deviations reflect the uneven exposure of mixture components to the solvent and the linear relation is recovered if fractions of solvent-accessible surface are used as the measure of composition. As the spatial variations in polarity become of larger amplitude, the linear relation can no longer be obtained. Protein surfaces represent such natural patterned surfaces, also characterized by larger patches and roughness. Our calculations reveal strong deviations from linear additivity on a prototypical surface comprising surface fragments of melittin dimer. The deviations reflect the disproportionately strong influence of isolated polar patches, preferential wetting, and changes in the position of the liquid interface above hydrophobic patches. Because solvent-induced contribution to the free energy of surface association grows as cos Θ(c), deviations of cos Θ(c) from the linear relation directly reflect nonadditive adhesive energies of biosurfaces.

Reaction of gas phase OH with unsaturated self-assembled monolayers and relevance to atmospheric organic oxidations

The kinetics and mechanisms of the reaction of gas phase OH radicals with organics on surfaces are of fundamental chemical interest, as well as relevant to understanding the degradation of organics on tropospheric surfaces or when they are components of airborne particles. We report here studies of the oxidation of a terminal alkene self-assembled monolayer (7-octenyltrichlorosilane, C8= SAM) on a germanium attenuated total reflectance crystal by OH radicals at a concentration of 2.1 x 10(5) cm(-3) at 1 atm total pressure and 298 K in air. Loss of the reactant SAM and the formation of surface products were followed in real time using infrared spectroscopy. From the rate of loss of the C=C bond, a reaction probability within experimental error of unity was derived. The products formed on the surface include organic nitrates and carbonyl compounds, with yields of 10 +/- 4% and < or = 7 +/- 4%, respectively, and there is evidence for the formation of organic products with C-O bonds such as alcohols, ethers and/or alkyl peroxides and possibly peroxynitrates. The yield of organic nitrates relative to carbonyl compounds is higher than expected based on analogous gas phase mechanisms, suggesting that the branching ratio for the RO(2) + NO reaction is shifted to favor the formation of organic nitrates when the reaction occurs on a surface. Water uptake onto the surface was only slightly enhanced upon oxidation, suggesting that oxidation per se cannot be taken as a predictor of increased hydrophilicity of atmospheric organics. These experiments indicate that the mechanisms for the surface reactions are different from gas phase reactions, but the OH oxidation of surface species will still be a significant contributor to determining their lifetimes in air.