Hydration and dewetting near fluorinated superhydrophobic plates

The water dynamics near nanoscale fluorinated (CF(3)(CF(2))(7)(CH(2))(2)SiH(3)) monolayers (plates) as well as possible dewetting transitions in-between two such plates have been studied with molecular dynamics simulations in this paper. A "weak water depletion" is found near the single fluorinated surface, with an average water density in the first solvation shells 6-8% lower than its hydrogenated counterpart. The fluorinated molecules are also found to be water impermeable, consistent with experimental findings. More surprisingly, a dewetting transition is found in the interplate region with a critical distance D(c) of 10 A (3-4 water diameters) for double plates with 8 x 8 molecules each (plate size approximately 4 nm x 4 nm). This transition, although occurring on a microscopic length scale, is reminiscent of a first-order phase transition from liquid to vapor. The unusual superhydrophobicity of fluorocarbons is found to be related to their larger size (or surface area) as compared to hydrocarbons, which "dilutes" their physical interactions with water. The water-plate interaction profile shows that the fluorinated carbons have a 10-12% weaker water-plate interaction than their hydrogenated counterparts in the nearest solvation shell, even though the fluorocarbons do have a stronger electrostatic interaction with water due to their larger partial charges. However, the van der Waals interactions dominate the water-plate interaction within the nearest shell, with up to 90% contributions to the total interaction energy, and fluorocarbons have a noticeably weaker (by 10-15%) van der Waals interaction with water in the nearest shell than do hydrocarbons. Both the slightly weaker water-plate interaction and larger surface area contribute to the stronger dewetting transition in the current fluorinated carbon plates.     

Hydration and dewetting near graphite-CH(3) and graphite-COOH plates

The dynamics of water near the nanoscale hydrophobic (graphite-CH(3)) and hydrophilic (graphite-COOH) plates has been studied in detail with molecular dynamics simulations in this paper. It is shown that these designed surfaces (by growing a layer of methyl or carboxyl groups on top of graphite) can have a significant impact on the neighboring water dynamics, with the hydrophilic carboxyl surface having even more profound effects. The water hydrogen bond lifetime is much longer near both types of surfaces than that in the bulk, while on the other hand the water diffusion constant is much smaller than that in the bulk. The difference in the diffusion constant can be as large as a factor of 8 and the difference in the hydrogen bond lifetime can be as large as a factor of 2, depending on the distance from the surface. Furthermore, the water molecules in the first solvation shell of surface atoms show a strong bias in hydroxyl group orientation near the surface, confirming some of the previous findings. Finally, the possible water dewetting transition between two graphite-CH(3) plates and the effect of the strength of the solute-solvent attractions on the water drying transition are investigated. The relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on the graphite-CH(3) surface is also analyzed on the basis of a simple macroscopic theory, which can be used to predict the dewetting transition critical distance.     

水与乙醇分子在金管内吸附作用的分子动力学研究

使用分子动力学研究了乙醇与水分子在纳米金管内按照不同比例混合时的吸附现象,并利用径向密度分布函数及水和乙醇分子所形成的平均氢键数来探讨纳米限制效应.结果表明,径向密度分布函数和氢键数目受纳米金管影响较大.另外,水与金管之间的作用力比乙醇与金管之间的大,导致水分子形成的平均氢键数不同于乙醇分子的.     

Water inside a hydrophobic cavitand molecule

The structure and dynamics of water inside a water-soluble, bowl-shaped cavitand molecule with a hydrophobic interior are studied using molecular dynamics computer simulations. The simulations find that the number of inside water molecules is about 4.5, but it fluctuates from being completely empty to full on a time scale of tens of nanoseconds. The transition from empty to full is energetically favorable and entropically unfavorable. The water molecules inside have fewer hydrogen bonds than the bulk and in general weaker interactions; the lower energy results from the nearest-neighbor interactions with the cavitand atoms and the water molecules at the entrance of the cavitand, interactions that are lost upon dewetting. An analysis of translational and rotational motion suggests that the lower entropy of the inside water molecules is due to decreased translational entropy, which outweighs an increased orientational entropy. The cavitand molecule acts as a host binding hydrophobic guests, and dewetting can be induced by the presence of a hydrophobic guest molecule about 3 A above the entrance. At this position, the guest displaces the water molecules which stabilize the inside water molecules and the empty cavitand becomes more stable than the full.     

Theoretical investigation of gas phase ethanol–(water)n (n = 1–5) clusters and comparison with gas phase pure water clusters (water)n (n = 2–6)

Various properties (such as optimal structures, structural parameters, hydrogen bonds, natural bond orbital charge distributions, binding energies, electron densities at hydrogen bond critical points, cooperative effects, and so on) of gas phase ethanol–(water)_n (n = 1–5) clusters with the change in the number of water molecules have been systematically explored at the MP2/aug‐cc‐pVTZ//MP2/6‐311++G(d,p) computational level. The study of optimal structures shows that the most stable ethanol‐water heterodimer is the one where exists one primary hydrogen bond (O? H…O) and one secondary hydrogen bond (C? H …O) simultaneously. The cyclic geometric pattern formed by the primary hydrogen bonds, where all the molecules are proton acceptor and proton donor simultaneously, is the most stable configuration for ethanol–(water)_n (n = 2–4) clusters, and a transition from two‐dimensional cyclic to three‐dimensional structures occurs at n = 5. At the same time, the cluster stability seems to correlate with the number of primary hydrogen bonds, because the secondary hydrogen bond was extremely weaker than the primary hydrogen bond. Furthermore, the comparison of cooperative effects between ethanol–water clusters and gas phase pure water clusters has been analyzed from two aspects. First of all, for the cyclic structure, the cooperative effect in the former is slightly stronger than that of the latter with the increasing of water molecules. Second, for the ethanol–(water)₅ and (water)₆ structure, the cooperative effect in the former is also correspondingly stronger than that of the latter except for the ethanol–(water)₅ book structure. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulations of the hydration of poly(vinyl methyl ether):Hydrogen bonds and quasi-hydrogen bonds

Atomistic detailed hydration structures of poly(vinyl methyl ether)(PVME) have been investigated by molecular dynamics simulations under 300 K at various concentrations. Both radial distribution functions and the distance distributions between donors and acceptors in hydrogen bonds show that the hydrogen bonds between the polymer and water are shorter by 0.005 nm than those between water molecules. The Quasi-hydrogen bonds take only 7.2% of the van der Waals interaction pairs. It was found the hydrogen bonds are not evenly distributed along the polymer chain,and there still exists a significant amount(10%) of ether oxygen atoms that are not hydrogen bonded to water at a concentration as low as 3.3%. This shows that in polymer solutions close contacts occur not only between polymer chains but also between chain segments within the polymer,which leads to inefficient contacts between ether oxygen atoms and water molecules. Variation of the quasi-hydrogen bonds with the concentration is similar to that of hydrogen bonds,but the ratio of the repeat units forming quasi-hydrogen bonds to those forming hydrogen bonds approaches 0.2. A transition was found in the demixing enthalpy at around 30% measured by dynamic testing differential scanning calorimetry(DTDSC) for aqueous solutions of a mono-dispersed low molecular weight PVME,which can be related to the transition of the fractions of hydrogen bonds and quasi-hydrogen bonds at ～27%. The transition of the fractions of hydrogen bonds and quasi-hydrogen bonds at ～27% can be used to explain the demixing enthalpy transition at 30% at a molecular scale. In addition,at the concentration of 86%,each ether oxygen atom bonded with water is assigned 1.56 water molecules on average,and 'free' water molecules emerge at the concentration of around 54%.     

The effects of confinement on the behavior of water molecules between parallel Au plates of (001) planes

Molecular dynamics simulation is utilized to investigate the behavior of water molecules confined between two Au plates of (001) planes separated by gaps of 24.48, 16.32, 12.24, 11.22, and 10.20 A. The simulation results indicate that the arrangements of the water molecules are dependent on the gap size. For the largest gap size, adsorption of the Au surface creates two permanent water layers in the vicinity of each Au plate. Furthermore, in this case, the gap size is sufficiently large to permit the formation of a central region within which the water molecules are randomly oriented in a similar manner to bulk water molecules. The results indicate that the orientation of the first water layer directly absorbed by the plate surface does not change as the gap size between the two Au plates is reduced. However, the orientations of the O-H bonds in the second water layer parallel to the surface rearrange to form hydrogen bonds between the water layers as the separation between the plates is decreased. Finally, an inspection of the variation of the self-diffusion coefficients with the gap size suggests that the difference between the dynamic properties of the water molecules in the z direction and the x-y plane decreases as the distance between the two Au plates increases.     

Sensitivity of hydrogen bond lifetime dynamics to the presence of ethanol at the interface of a phospholipid bilayer

Atomistic molecular dynamics simulations of a fully hydrated liquid crystalline lamellar phase of a dimyrystoylphosphatidylcholine lipid bilayer containing ethanol at 1:1 composition as well as of the pure lamellar phase of the bilayer have been performed. Detailed analyses have been carried out to investigate the effects of ethanol, if any, on the lifetime dynamics of lipid-water and water-water hydrogen bonds in the hydration layer of the lipid headgroups. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed the identification of the bound states at the bilayer interface and the quantification of the dynamic equilibrium between the bound and the free water molecules, in terms of time-dependent relaxation rates. The calculations show that the overall relaxation of phosphate-water hydrogen bonds is faster in the presence of ethanol. Studies of the residence time and the number fluctuation of the hydration layer water molecules reveal that the presence of ethanol molecules decreases the rigidity of the lipid hydration layer.     

Perturbation of phospholipid bilayer properties by ethanol at a high concentration

Atomistic molecular dynamics (MD) simulations have been carried out at 30 degrees C on a fully hydrated liquid crystalline lamellar phase of dimyrystoylphosphatidylcholine (DMPC) lipid bilayer with embedded ethanol molecules at 1:1 composition, as well as on the pure bilayer phase. The ethanol molecules are found to exhibit a preference to occupy regions near the upper part of the lipid acyl chains and the phosphocholine headgroups. The calculations revealed that the phosphocholine headgroup dipoles (P- --> N+) of the lipids prefer to orient more toward the aqueous layer in the presence of ethanol. It is noticed that the ethanol molecules modify the dynamic properties of both lipids as well as the water molecules in the hydration layer of the lipid headgroups. Both the in-plane "rattling" and out-of-plane "protrusion" motions of the lipids have been found to increase in the presence of ethanol. Most importantly, it is observed that the water molecules within the hydration layer of the lipid headgroups exhibit faster translational and rotational motions in the presence of ethanol. This arises due to faster dynamics of hydrogen bonds between lipid headgroups and water in the presence of ethanol.     

Study of water clusters in the n = 2–34 size regime, based on the ABEEM/MM model

Various properties of water clusters in the n = 2–34 size regime with the change of cluster size have been systemically explored based on the newly developed flexible-body and charge-fluctuating ABEEM/MM water potential model. The ABEEM/MM water model is to take ABEEM charges of all atoms, bonds, and lone-pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The computed correlating properties characterizing water clusters (H₂O) _n (n = 2–34) include optimal structures, structural parameters, ABEEM charge distributions, binding energies, hydrogen bonds, dipole moments, and so on. The study of optimal structures shows that the ABEEM/MM model can correctly predict the following important structural features, such as the transition from two-dimensional (from dimer to pentamer) to three-dimensional (for clusters larger than the hexamer) structures at hexamer region, the transition from cubes to cages at dodecamer (H₂O)₁₂, the transition from all-surface (all water molecules on the surface of the cluster) to one water-centered (one water molecule at the center of the cluster, fully solvated) structures at (H₂O)₁₇, the transition from one to two internal molecules in the cage at (H₂O)₃₃, and so on. The first three structural transitions are in good agreement with those obtained from previous work, while the fourth transition is different from that identified by Hartke. Subsequently, a systematic investigation of structural parameters, ABEEM charges, energetic properties, and dipole moments of water clusters with increasing cluster size can provide important reference for describing the objective trait of hydrogen bonds in water cluster system, and also provide a strong impetus toward understanding how the water clusters approach the bulk limit.     

Variations in the concentration of microwave dielectric spectra in a water-acetonitrile system

Relaxation models for describing the microwave dielectric spectra of acetonitrile (AN) solutions are presented. It is noted that the models, which contain two fundamental Debye spectral components, can be used for most concentrations. The separate influence of polar and nonpolar molecules in groups of AN molecules on the initial structure of water in the solution is established and found to have different signs. A schematic of the concentration’s structural transition from water to pure AN is given. It is concluded that integrated clusters consisting of water and AN molecules are found in the system’s region of medium concentrations (0.4–0.6 m.f. AN), where the motions of water molecules are slowed down and it is stabilized due to the hydrophobic hydration of nonpolar AN groups. It is shown that the AN molecules are more tightly bound compared than pure AN due to the formation of AN-H₂O bonds over the range of concentrations.     

Effects of Hydrogen Bonding on the Transition Properties of Ethanol–Water Clusters: A TD-DFT Study

In this work, time-dependent density functional theory method was used to study the electronic transitions of hydrogen-bonded ethanol–water complexes Dimer-I, Dimer-II and Trimer. The intermolecular hydrogen bonds H₁···O₁ and O···H₂ were demonstrated by the optimized geometric structures of the three hydrogen-bonded ethanol–water complexes. It is demonstrated that the S₁-state electronic transitions for ethanol monomer and the hydrogen-bonded complex Dimer-I (through HB-I) should be of LE nature on the ethanol molecule, while those of complexes Dimer-II and Trimer should be of CT character from the hydrogen-bonded water molecule (through HB-II) to the ethanol moiety. The different electronic transition types should be the reasons for the tiny redshift of the S₁-state electronic energy for Dimer-I and the large blueshifts for Dimer-II and the Trimer compared with that of the ethanol monomer.     

Study on the thermoresponsive two phase transition processes of hydroxypropyl cellulose concentrated aqueous solution: from a microscopic perspective

In this paper, it was discovered that during the heating process from 35 to 63 °C, hydroxypropyl cellulose (HPC) concentrated aqueous solution (20 wt%) would first go through coil-to-globule transition and then sol–gel transition with temperature elevation. The microdynamic mechanisms of the two phase transitions were thoroughly illustrated using mid and near infrared spectroscopy in combination with two-dimensional correlation spectroscopy (2Dcos) and perturbation correlation moving window (PCMW) technique. Mid infrared spectroscopy is an effective way to study the hydrophobic interactions in HPC molecules. And near infrared spectroscopy is a potent method to study hydrogen bonds between HPC molecules and water molecules. Boltzmann fitting and PCMW could help determine the exact transition temperatures of each involving functional groups in the two processes. Moreover, 2Dcos was used to discern the sequential moving orders of the functional groups during the two phase transitions. Depending on the structure of HPC and the thermodynamic conditions, the dominating associative elements in either process might vary. During the coil-to-globule transition, HPC molecules precipitated to form an opaque system with mobility.It was discovered that the driving force of the coil-to-globule transition process in microdynamics could only be the dehydration and hydrophobic interactions of C–H groups. However, in the sol–gel transition, the system crosslinked to form a physical network with no mobility. The driving force of this process in microdynamics was primarily the self-assembly behavior of O–H groups in HPC “active molecules”.     

Properties of a water layer on hydrophilic and hydrophobic self-assembled monolayer surfaces: A molecular dynamics study

The microscopic behaviors of a water layer on different hydrophilic and hydrophobic surfaces of well ordered self-assembled monolayers (SAMs) are studied by molecular dynamics simulations. The SAMs consist of 18-carbon alkyl chains bound to a silicon(111) substrate, and the characteristic of its surface is tuned from hydrophobic to hydrophilic by using different terminal functional groups ( CH 3 , COOH). In the simulation, the properties of water membranes adjacent to the surfaces of SAMs were reported by comparing pure water in mobility, structure, and orientational ordering of water molecules. The results suggest that the mobility of water molecules adjacent to hydrophilic surface becomes weaker and the molecules have a better ordering. The distribution of hydrogen bonds indicates that the number of water-water hydrogen bonds per water molecule tends to be lower. However, the mobility of water molecules and distribution of hydrogen bonds of a water membrane in hydropho- bic system are nearly the same as those in pure water system. In addition, hydrogen bonds are mainly formed between the hydroxyl of the COOH group and water molecules in a hydrophilic system, which is helpful in understanding the structure of interfacial water.     

Hydration state of nonionic surfactant monolayers at the liquid/vapor interface: structure determination by vibrational sum frequency spectroscopy

The OH stretching region of water molecules in the vicinity of nonionic surfactant monolayers has been investigated using vibrational sum frequency spectroscopy (VSFS) under the polarization combinations ssp, ppp, and sps. The surface sensitivity of the VSFS technique has allowed targeting the few water molecules present at the surface with a net orientation and, in particular, the hydration shell around alcohol, sugar, and poly(ethylene oxide) headgroups. Dramatic differences in the hydration shell of the uncharged headgroups were observed, both in comparison to each another and in comparison to the pure water surface. The water molecules around the rigid glucoside and maltoside sugar rings were found to form strong hydrogen bonds, similar to those observed in tetrahedrally coordinated water in ice. In the case of the poly(ethylene oxide) surfactant monolayer a significant ordering of both strongly and weakly hydrogen bonded water was observed. Moreover, a band common to all the surfactants studied, clearly detected at relatively high frequencies in the polarization combinations ppp and sps, was assigned to water species located in proximity to the surfactant hydrocarbon tail phase, with both hydrogen atoms free from hydrogen bonds. An orientational analysis provided additional information on the water species responsible for this band.     

Structural studies on the hydration of L-glutamic acid in solution

A combination of neutron diffraction augmented with isotopic substitution and computer modeling using empirical potential structure refinement has been used to extract detailed structural information for L-glutamic acid dissolved in 2 M NaOH solution. This work shows that the tetrahedral hydrogen bonding network in water is severely disrupted by the addition of glutamic acid and NaOH, with the number of water-water hydrogen bonds being reduced from 1.8 bonds per water molecule in pure water to 1.4 bonds per water molecule in the present solution. In the glutamic acid molecule, each carboxylate oxygen atom forms an average of three hydrogen bonds with the surrounding water solvent with one of these hydrogens being shared between the two oxygen atoms on each carboxylate group, while each amine hydrogen forms a single hydrogen bond with the surrounding water solvent. Additionally, the average conformation of the glutamic acid molecules in these solutions is extracted.     

Local heterogeneous dynamics of water around lysozyme: a computer simulation study

Water present near the surface of a protein exhibits dynamic properties different from that of water in the pure bulk state. In this work, we have carried out atomistic molecular dynamics simulation of an aqueous solution of hen egg-white lysozyme. Attempts have been made to explore the correlation between the local heterogeneous mobility of water around the protein segments and the rigidity of the hydration layers with the microscopic dynamics of hydrogen bonds formed by water molecules with the protein residues. The kinetics of breaking and reformation of hydrogen bonds involving the surface water molecules have been calculated. It is found that the reformations of broken hydrogen bonds are more frequent for the hydration layers of those segments of the protein that are more rigid. The calculation of the low-frequency vibrational modes of hydration layer water molecules reveals that the protein influences the transverse and longitudinal degrees of freedom of water around it in a differential manner. These findings can be verified by appropriate experimental studies.     

Spectroscopic characterization of microscopic hydrogen-bonding disparities in supercritical water

The local hydrogen-bonding environment in supercritical water (380 degrees C, 300 bars, density 0.54 gcm3) was studied by x-ray Raman scattering at the oxygen K edge. The spectra are compared to those of the gas phase, liquid surface, bulk liquid, and bulk ice, as well as to calculated spectra. The experimental model systems are used to assign spectral features and to quantify specific local hydrogen-bonding situations in supercritical water. The first coordination shell of the molecules is characterized in more detail with the aid of the calculations. Our analysis suggests that approximately 65% of the molecules in supercritical water are hydrogen bonded in configurations that are distinctly different from those in liquid water and ice. In contrast to liquid water the bonded molecules in supercritical water have four intact hydrogen bonds and in contrast to ice large variations of bond angles and distances are observed. The remaining approximately 35% of the molecules exhibit two free O-H bonds and are thus either not involved in hydrogen bonding at all or have one or two hydrogen bonds on the oxygen side. We determine an average O-O distance of 3.1+/-0.1 A in supercritical water for the H bonded molecules at the conditions studied here. This and the corresponding hydrogen bond lengths are shown to agree with neutron- and x-ray-diffraction data at similar conditions. Our results on the local hydrogen-bonding environment with mainly two disparate hydrogen-bonding configurations are consistent with an extended structural model of supercritical water as a heterogeneous system with small patches of bonded molecules in various tetrahedral configurations and surrounding nonbonded gas-phase-like molecules.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.