Hydrogen-bond dynamics of interior and surface molecules in a water nanocluster: temperature and size effects

Molecular dynamics simulations are employed to investigate the effects of temperature and size on the hydrogen-bond dynamics of interior molecules and surface molecules in a water nanocluster. The flexible three-centred (F3C) water model is invoked in the simulations. To inspect the dynamics of the interior hydrogen bonds and the surface hydrogen bonds, a spherical water nanocluster is modelled and then divided into interior molecules and surface molecules according to the density profile of the water nanocluster. It is observed that at higher temperatures the average number of hydrogen bonds decreases and yields faster hydrogen-bond relaxation for both interior molecules and surface molecules of the water nanocluster. Furthermore, the surface molecules have a lower average number of hydrogen bonds than the interior molecules. The lifetime of the surface hydrogen bonds is slightly longer than that of the interior hydrogen bonds, whereas the hydrogen-bond structural relaxation time of the surface molecules is more obviously slower than that of the interior molecules. Regarding the size effect, a larger water nanocluster is seen to have a larger average number of hydrogen bonds and a longer hydrogen-bond structural relaxation time.     

界面氢键对受限水结构和动态特性的影响

应用反应力场分子动力学方法, 模拟了水限制在全羟基化二氧化硅晶体表面间的弛豫过程, 研究了基底表面与水形成的界面氢键, 及其对受限水结构和动态特性行为的影响. 当基底表面硅醇固定时, 靠近基底表面水分子中的氧原子与基底表面的氢原子形成强氢键, 这使得靠近表面水分子中的氧原子比对应的氢原子更靠近基底表面, 从而水分子的偶极矩远离表面. 当基底表面硅醇可动时, 靠近基底表面水分子与基底表面原子形成两种强氢键, 一种是水分子中的氧原子与表面的氢原子形成的强氢键, 数量较少, 另一种是水分子中的氢原子与表面的氧原子形成的强氢键, 数量较多, 这使得靠近表面水分子中的氢原子比对应的氧原子更靠近表面, 从而水分子的偶极矩指向表面. 在相同几何间距下, 当基底表面硅醇可动时, 表面的活动性使得几何限制作用减弱, 导致了受限水分层现象没有固定表面限制下的明显. 此外, 固定表面比可动表面与水形成的界面氢键作用较强, 数量较多, 导致了可动表面限制下水的运动更为剧烈.     

Role of oxygen functional groups for structure and dynamics of interfacial water on low rank coal surface: a molecular dynamics simulation

Molecular dynamics simulations were employed to study the effects of oxygen functional groups for structure and dynamics properties of interfacial water molecules on the subbituminous coal surface. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of subbituminous coal according to XPS results, and the composing proportion for hydroxyl, carbonyl and carboxyl is 25:3:5. The hydration energy with ?386.28 kJ/mol means that the adsorption process between water and coal surface is spontaneous. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxyl and carbonyl is similar.     

Enhancement of water self-diffusion at super-hydrophilic surface with ordered water

It has been well acknowledged that molecular water structures at the interface play an important role in the surface properties, such as wetting behavior or surface frictions. Using molecular dynamics simulation, we show that the water self-diffusion on the top of the first ordered water layer can be enhanced near a super-hydrophilic solid surface. This is attributed to the fewer number of hydrogen bonds between the first ordered water layer and water molecules above this layer, where the ordered water structures induce much slower relaxation behavior of water dipole and longer lifetime of hydrogen bonds formed within the first layer.     

Molecular dynamics simulations of PVP kinetic inhibitor in liquid water and hydrate/liquid water systems

The possible effects of PVP (poly(N-vinylpyrrolidone)) on the properties of liquid and water in clathrate hydrate has been investigated using NVT molecular dynamics simulations. A model for a monomer of the PVP polymer is immersed in three systems, liquid water, a unit cell of a hydrate in liquid water with a hydrate former and a third system where some of the liquid water molecules of this last system are replaced by a PVP monomer. Both molecular dynamics simulation and integral equation theory predict hydrogen bonding between the double bonded oxygen in the PVP ring and hydrogens in water. For the composite system, the PVP monomer has a preference for hydrogen bonding to hydrogens from the water molecules at the surface of the hydrate lattice. The simulations indicate that the PVP monomer tends to orient perpendicular to the hydrate surface. For the model systems in this study PVP may form hydrogen bonds with liquid water through the double bonded oxygen in the ring. When a hydrate crystal is immersed in the liquid water phase this hydrogen bonding is shifted towards the hydrate due to a more favourable Coulomb interaction involving hydrogens from more than one water molecule at the hydrate surface. The PVP monomer has a preference for perpendicular orientation with respect to the hydrate surface. A scheme is suggested for the characterization of kinetic hydrate inhibitors based on molecular dynamics simulations and on three basic properties. In addition to the energy between the active groups of the inhibitor and hydrate water another point of focus is the free energy changes in the interactions between the inhibitor and water as the charges are changed from zero to the original model charges. In particular the difference between this integral for the (hydrate water)–(PVP monomer) interaction and the (liquid water)–(PVP inhibitor) interaction should reflect the driving forces in freezing the inhibitor out from the liquid water phase and onto the hydrate surface. The third property in the characterization scheme is the diffusivities of groups connecting to the hydrate crystal, relative to the diffusivities of the hydrate crystal. Results are presented from simulations where a small cavity with a methane model as a guest is immersed in liquid water with free methane molecules at a temperature of 150K. Changes in structure, diffusivities and energy indicate a tendency towards a more solid–like structurde around the cavity.     

Water in porous glasses. A computer simulation study

We report molecular dynamics simulations of water confined in a cylindrical silica pore. The pore geometry and size is similar to that of typical pores in porous Vycor glass. In the present study we focus on the dependence of microscopical structural and dynamical properties on the degree of hydration of the pore. We have performed five simulations of systems between 19 and 96 % hydration. In all cases, water adsorbs strongly on the pore surface, clearly demonstrating the hydrophilic nature of the Vycor surface. Two layers of water molecules are affected strongly by the interactions with the glass surface. With decreasing degree of hydration an increasing volume in the center of the pore is devoid of water molecules. At 96 % hydration the center is a continuous and homogeneous region that has, however, a lower density than bulk water at ambient conditions. A well-pronounced mobility profile exists, where molecules in the center of the pores have substantially higher self diffusion coefficients than molecules on the pore surface. The spectral densities of center of mass and hydrogen atom motion show the signature of confinement for the molecules close to the pore surface, while the spectral densities in the center of the pore are similar to those in bulk water. The molecular dynamics results are in good agreement with recent experiments. Our data indicate that the dependence of experimental data on the level of hydration of the Vycor sample is due to the different relative contribution of molecules adsorbed on the pore surface and bulk-like molecules in the interior of the pore to the experimental averages.     

表/界面水的扫描探针技术研究进展

表面和界面水在自然界、人们的日常生活以及现代科技中无处不在.它在物理、化学、环境学、材料学、生物学、地质学等诸多基础学科和应用领域起到至关重要的作用.因此,表面和界面水的功能与特性的研究,是水基础科学的一项核心任务.然而,由于水分子之间氢键相互作用的复杂性,及其与水-固界面相互作用的竞争,使得表(界)面水对于局域环境的影响非常敏感,往往需要深入到分子层次研究其微观结构和动力学过程.近年来,新型扫描探针技术的发展使得人们可以在单分子甚至亚分子尺度上对表(界)面水展开细致的实空间研究.本文着重介绍几种代表性的扫描探针技术及其在表(界)面水体系中的应用,包括:超高真空扫描隧道显微术、单分子振动谱技术、电化学扫描隧道显微术和非接触式原子力显微术.此外,本文还将对表(界)面水扫描探针技术研究面临的挑战和未来发展方向进行了展望.     

Oxide/water interfaces: how the surface chemistry modifies interfacial water properties

The organization of water at the interface with silica and alumina oxides is analysed using density functional theory-based molecular dynamics simulation (DFT-MD). The interfacial hydrogen bonding is investigated in detail and related to the chemistry of the oxide surfaces by computing the surface charge density and acidity. We find that water molecules hydrogen-bonded to the surface have different orientations depending on the strength of the hydrogen bonds and use this observation to explain the features in the surface vibrational spectra measured by sum frequency generation spectroscopy. In particular, 'ice-like' and 'liquid-like' features in these spectra are interpreted as the result of hydrogen bonds of different strengths between surface silanols/aluminols and water.     

Structure and dynamics of ordered water in a thick nanofilm on ionic surfaces

The structure and dynamics of water in a thick film on an ionic surface are studied by molecular dynamic simulations. We find that there is a dense monolayer of water molecules in the vicinity of the surface. Water molecules within this layer not only show an upright hydrogen-down orientation, but also an upright hydrogen-up orientation. Thus, water molecules in this layer can form hydrogen bonds with water molecules in the next layer. Therefore, the two-dimensional hydrogen bond network of the first layer is disrupted, mainly due to the O atoms in this layer, which are affected by the next layer and are unstable. Moreover, these water molecules exhibit delayed dynamic behavior with relatively long residence time compared with those bulk-like molecules in the other layers. Our study should be helpful to further understand the influence of water film thickness on the interfacial water at the solid-liquid interface.     

Bound water: Evidence from and implications for the dielectric properties of aqueous solutions

Broadband dielectric spectra of a variety of aqueous solutions are evaluated as to indications of water that may be considered bound. Static permittivity decrements due to depolarizing internal electric fields, from kinetic depolarization, as well as from dielectric saturation are discussed. The latter effect reflects the preferential orientation of water permanent dipole orientations within strong Coulombic field of small ions, especially multivalent cations. Such water may be considered bound even though rapid rotations around the orientation of the electric dipole moment are definitely possible and also a fast exchange of water molecules between the hydration region and the bulk may take place. Water exhibiting large dielectric relaxation times, as typical for regions with large local concentration of foreign matter, may also be named bound. However, no clear evidence for interaction energies exceeding the hydrogen bond energy of pure water has been found. Rather enhanced relaxation times at low water content reflect the small concentration of hydrogen bonding sites and thus low probability density for the formation of a new hydrogen bond. Potential interferences of the water relaxation with relaxations from other molecules or from ionic structures are mentioned briefly.     

Effects of co-solutes on the hydrogen bonding structure and dynamics in aqueous N-methylacetamide solution: a molecular dynamics simulations study

The effects of trimethylamine-N-oxide (TMAO), urea and tetramethyl urea (TMU) on the hydrogen bonding structure and dynamics of aqueous solution of N-methylacetamide (NMA) are investigated by classical molecular dynamics simulations. The modification of the water's hydrogen bonding structure and interactions is calculated in presence of these co-solutes. It is observed that the number of four-hydrogen-bonded water molecules in the solution decreases significantly in the presence of TMAO rather than urea and TMU. The lifetime and structural relaxation time of water–water and NMA–water hydrogen bonds show a strong increase with the addition of TMAO and TMU in the solution, whereas the change is nominal in case of urea solution. It is also found that the translational and rotational dynamics of water and NMA slowdown with increasing the concentration of these osmolytes. The slower dynamics of water and NMA is more pronounced in case of TMAO and TMU solution, as these co-solutes strengthen the average hydrogen bond energies between water–water and NMA–water, whereas urea has a little effect on the hydrogen bonding structure and dynamics of aqueous NMA solution. The calculated self-diffusion coefficient values for water and these co-solutes are in similar pattern with experimental observations.     

Comparison of reactivity on step and terrace sites of Pd (3 3 2) surface for the dissociative adsorption of hydrogen: A quantum chemical molecular dynamics study

The notion of “active sites” is fundamental to heterogeneous catalysis. However, the exact nature of the active sites, and hence the mechanism by which they act, are still largely a matter of speculation. In this study, we have presented a systematic quantum chemical molecular dynamics (QCMD) calculations for the interaction of hydrogen on different step and terrace sites of the Pd (3 3 2) surface. Finally the dissociative adsorption of hydrogen on step and terrace as well as the influence of surface hydrogen vacancy for the dissociative adsorption of hydrogen has been investigated through QCMD. This is a state-of-the-art method for calculating the interaction of atoms and molecules with metal surfaces. It is found that fully hydrogen covered (saturated) step sites can dissociate hydrogen moderately and that a monovacancy surface is suitable for significant dissociative adsorption of hydrogen. However in terrace site of the surface we have found that dissociation of hydrogen takes place only on Pd sites where the metal atom is not bound to any pre-adsorbed hydrogen atoms. Furthermore, from the molecular dynamics and electronic structure calculations, we identify a number of consequences for the interpretation and modeling of diffusion experiments demonstrating the coverage and directional dependence of atomic hydrogen diffusion on stepped palladium surface.     

Solvation properties of non-polar amino acids in water and methanol: a molecular dynamics study

Solvation properties of the amino acids alanine, valine, leucine and phenylalanine in water and in methanol at infinite dilution have been studied by using molecular dynamics simulations. The free energy transfer of the solutes from water to methanol has been calculated by means of a slow-growth technique and the results compared with diverse hydrophobicity scales available in the literature. Chothia's linear relation between the accessible surface area and the hydrophobicity is shown. The line obtained has a slope equivalent to 49cal mol^?1 Å^?2. The radial distribution of solvent molecules around the amino acids residues is considered, and a comparison is made of the number of hydrogen bonds formed between water molecules in the presence and in the absence, respectively, of the amino acids. The calculations show clearly that the solvation structure near the non-polar amino acids is richer for methanol than for water.     

A molecular dynamics simulation study on the cavitation inception of water with dissolved gases

The promotion/prevention mechanism of dissolved gases on cavitation inception is essential for many high-tech industries and research. In the present study, large-scale molecular dynamics simulations are performed to investigate the effects of water cavitation caused by different gas types by using nitrogen and oxygen gases with TIP4P/2005 water. The cavitation inception behaviour is analyzed via Mean First Passage Time method. Water with dissolved gases has a higher nucleation rate and is easier to cavitate than pure water. At the same gas concentration, the cavitation of water with nitrogen is promoted to a greater extent than that with oxygen. The number and energy of hydrogen bond (HB) are further calculated by the Acceptor-Hydrogen-Donor method to explain this promotion mechanism. The number and energy of HB in water with gases decrease compared with those in pure water. The introduction of gases weakens the HB network and promotes cavitation inception because of weaker interactions between gas and water molecules. A model is developed to describe the relationship between nucleation rate and HB energy. Gas molecules assemble on the surface of bubbles during water cavitation, which may decrease the free energy of bubble surface, maintain the existing bubble, and contribute to the growth process.     

Dynamical and structural properties of adsorbed water molecules at the TiO2 rutile-(110) surface: interfacial hydrogen bonding probed by ab-initio molecular dynamics

ABSTRACT

Ab-initio molecular dynamics (AIMD) simulations have been carried out to study a range of different and energetically-accessible adsorbed-water configurations and motifs for their vibrational and structural characteristics, in contact with rutile-(110) interfaces at 100?K. The radial pair distribution function between the titanium atoms at the interface and the hydrogen and oxygen atoms in the water monolayer show an orientation of the water molecules parallel to the surface of titania, and with hydrogen atoms pointed in the opposite direction to the surface. In some cases, a distinctive vibrational frequency region between 2500 and 3000?cm^?1 has also been observed, due to a strong dispersion interaction between water molecules. This behaviour is also seen in experimental studies of thin-film water coverage on TiO₂ surfaces.     

How polar hydroxyl groups affect surface hydrophobicity on model talc surfaces

Based on molecular dynamics simulations, we have studied the wetting behaviors of water on the talc-like surface with different surface polarity by modifying the charge distribution of surface hydroxyl (–OH) groups. With the change of the charge of the hydrogen atom (denoted as δ_q) in –OH group, the contact angle decreases from 91° to 50° and then remains constant. On the surfaces with the larger charge of hydrogen atoms (δ_q ≥ 0.2 e), a water droplet is formed above a water monolayer, which is exactly contacted on the surface. Each water molecule in the monolayer forms one hydrogen bond (H bond) with surface –OH groups, without participating in any H bond with the water molecules within the monolayer or with the water molecules above the monolayer. The polarity of the –OH group also has a great influence on the dynamic behaviors of the interface water, such as residence time, hydrogen bond lifetime and self-diffusion coefficient. The diffusion of water molecules in the water monolayer near the highly polar surface is greatly suppressed, and the residence time of water molecules in the water monolayer even exceeds 12 ns.     

Ice tessellation on a hydroxylated silica surface

The adsorption of water on a fully hydroxylated silica surface is studied by using density-functional total-energy and molecular dynamics calculations. The (100) surface of beta(alpha)-cristobalite covered by geminal hydroxyls has been taken as the substrate. A well-ordered and stable two-dimensional ice with quadrangular and octagonal patterns of hydrogen bond (H-bond) networks-an ice tessellation-is found on the surface for the first time. With the vibrational recognition, the four water molecules in the quadrangle are found to be bonded by strong H bonds while the quadrangles are connected to each other by weak H bonds. This configuration is the most stable, because all the water molecules are fully saturated with H bonds either to each other or to the surface hydroxyl groups.     

石墨烯表面的特征水分子排布及其湿润透明特性的分子动力学模拟

石墨烯因其独特的分子构型、卓越的物理化学性能而受到广泛关注.本文首先利用分子动力学模拟比较了单层石墨烯、铜、二氧化硅三者表面的浸润性,除了接触角的比较,还分析了基底表面的水分子排布,得到石墨烯表面的特征水分子排布为:表面有两层密集的水分子层,其中靠近基底的密集水分子层中O—H键与垂直基底方向夹角集中在90°附近,并且基底表面的氢键几乎都垂直于基底.另一方面,本文研究了石墨烯浸润透明特性,发现在铜和二氧化硅上添加一层石墨烯,对铜的浸润性影响较小,对二氧化硅的浸润性影响很大,不仅使其上接触角显著增大,还使得基底表面的水分子排布呈现出类似单层石墨烯上的规律.本文使用分子动力学模拟方法从微观尺度验证了文献的实验结果,从基底表面水分子排布角度分析了石墨烯独特的浸润透明特性,为进一步开发石墨烯在微结构设计上的应用提供了理论指导.     

Influence of the reversible adsorption of methyl alcohol on magnetization reversal in ferromagnets

A change in the dynamics of 180° domain walls on the surface of a soft amorphous ferromagnet in methyl alcohol atmosphere is established by means of a magnetooptic method. A reversible decrease in the relaxation frequency of the domain walls near the surface in the presence of methyl alcohol admolecules is observed. This effect is related to the magnetic defects resulting from methyl alcohol adsorption on the ferromagnet surface, which proceeds through hydrogen bonding. Based on earlier data for the influence of the reversible adsorption of water molecules on the domain wall dynamics in ferromagnets, it is concluded that reversible adsorption through the mechanism of hydrogen bond formation considerably affects the domain wall dynamics in soft ferromagnets.     

Direct observation of hydrogen-bond exchange within a single water dimer

The dynamics of water dimers was investigated at the single-molecule level by using a scanning tunneling microscope. The two molecules in a water dimer, bound on a Cu(110) surface at 6 K, were observed to exchange their roles as hydrogen-bond donor and acceptor via hydrogen-bond rearrangement. The interchange rate is approximately 60 times higher for (H2O)2 than for (D2O)2, suggesting that quantum tunneling is involved in the process. The interchange rate is enhanced upon excitation of the intermolecular mode that correlates with the reaction coordinate.