Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

Hydrogen bond dynamics and microscopic structure of confined water inside carbon nanotubes

We have investigated the density and temperature dependences of microscopic structure and hydrogen bond dynamics of water inside carbon nanotubes (CNTs) using molecular dynamics simulation. The CNTs are treated as rigid, and smoothly truncated extended simple point charge water model is adopted. The results show that as the overall density increases, the atomic density profiles of water inside CNTs become sharper, the peaks shift closer to the wall, and a new peak of hydrogen atomic density appears between the first (outermost) and second layer. The intermittent hydrogen bond correlation function C(HB)(t) of water inside CNTs decays slower than that of bulk water, and the rate of decay decreases as the tube diameter decreases. C(HB)(t) clearly decays more slowly for the first layer of water than for other regions inside CNTs. The C(HB)(t) of the interlayer hydrogen bonds decays faster than those of the other regions and even faster than that of the bulk water. On the other hand, the hydrogen bond lifetimes of the first layer are shorter than those of the inner layer(s). Interlayer hydrogen bond lifetimes are clearly shorter than those of the constituent layers. As a whole, the hydrogen bond lifetimes of water inside CNTs are shorter than those of bulk water, while the relaxation of C(HB)(t) is slower for the confined water than for bulk water. In other words, hydrogen bonds of water inside CNTs break more easily than those of bulk water, but the water molecules remain in each other's vicinity and can easily reform the bonds.     

Radiowave dielectric investigation of water confined in channels of carbon nanotubes

Structure and dynamics of water confined in channels of diameter of few nanometer in size strongly differ from the ones of water in the bulk phase. Here, we present radiowave dielectric relaxation measurements on water-filled single-walled carbon nanotubes, with the aim of highlighting some aspects on the molecular electric dipole organization of water responding to high spatial confinement in a hydrophobic environment. The observed dielectric spectra, resulting into two contiguous relaxation processes, allow us to separate the confined water in the interior of the nanotubes from external water, providing support for the existence in the confinement region of water domains held together by hydrogen bonds. Our results, based on the deconvolution of the dielectric spectra due to the presence of a bulk and a confined water phase, furnish a significantly higher Kirkwood correlation factor, larger than the one of water in bulk phase, indicating a strong correlation between water molecules inside nanotubes, not seen in bulk water.     

Adsorption of water molecules inside a Au nanotube: a molecular dynamics study

A molecular dynamics simulation of water molecules through a Au nanotube with a diameter of 20 A at bulk densities 0.8, 1, and 1.2 gcm(3) has been carried out. The water molecules inside a nanoscale tube, unlike those inside a bulk tube, have a confined effect. The interaction energy of the Au nanotube wall has a direct influence on the distribution of water molecules inside the Au tube in that the adsorption of the water molecules creates shell-like formations of water. Moreover, the high number of adsorbed molecules has already achieved saturation at the wall of the Au nanotube at three bulk densities. This work compares the distribution percentage profiles of hydrogen bonds for different regions inside the tube. The structural characteristics of water molecules inside the tube have also been studied. The results reveal that the numbers of hydrogen bonds per water molecule influence the orientational order parameter q. In addition, the phenomenon of a group of molecules bonded inside the tube can be observed as the number of hydrogen bonds increase.     

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

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

Description of collective effects in computer models of water

The dynamics of a system containing 3456 water molecules in a cubic cell with periodic boundary conditions at 297 K was simulated. The time dependence of distances between oxygen atoms was examined for many pairs of molecules. These distances often oscillate around a certain average value over long periods of time (10 ps and longer). These average values can be about 2.8 Å (hydrogen bond) or much larger, up to 12–13 Å and more. This suggests that big groups of molecules are involved in a concerted motion. Lists of hydrogen bonds in 50 configurations divided by an interval of about 1 ps are compared. The average lifetime of a hydrogen bond is about 7 ps. The network of hydrogen bonds is colored according to their lifetimes for one of the configurations. The bonds that live longer than 7 ps form an infinite cluster. The bonds that live longer than 8 ps join to form a great number of finite clusters including several hundreds of nodes (molecules). These clusters contain few closed cycles. Even the bonds that live longer than 20 ps are united into clusters each containing two or three nodes (molecules). The self-diffusion coefficient for molecules involved in long-lived bonds is likely to be slightly smaller than that for molecules which do not participate in these bonds.     

Dynamics of water confined in the interdomain region of a multidomain protein

Molecular dynamics simulations are performed to study the dynamics of interfacial water confined in the interdomain region of a two-domain protein, BphC enzyme. The results show that near the protein surface the water diffusion constant is much smaller and the water-water hydrogen bond lifetime is much longer than that in bulk. The diffusion constant and hydrogen bond lifetime can vary by a factor of as much as 2 in going from the region near the hydrophobic domain surface to the bulk. Water molecules in the first solvation shell persist for a much longer time near local concave sites than near convex sites. Also, the water layer survival correlation time shows that on average water molecules near the extended hydrophilic surfaces have longer residence times than those near hydrophobic surfaces. These results indicate that local surface curvature and hydrophobicity have a significant influence on water dynamics.     

Investigation of water diffusion process in ethyl cellulose-based films by attenuated total reflectance Fourier transform infrared spectroscopy and two-dimensional correlation analysis

Two-dimensional correlation spectroscopy (2Dcos) analysis on the time-resolved attenuated total reflectance Fourier transform infrared spectroscopy has been employed to investigate the diffusion behavior of water in ethyl cellulose/triethyl citrate (EC/TEC) films with varied TEC content. The diffusion coefficients of water in the EC-based films are calculated from the diffusion curves according to the Fickian Diffusion Model. And they are observed to increase with TEC content, possibly caused by the increased free volume in the film matrix. In the 2Dcos analysis, the broad O–H stretching vibration region is split into at least four bands, which can be assigned to four different states of water molecules, that is, bulk water, cluster water, relatively free water and free water. Moreover, it is found that as water molecules disperse into the EC-based films, cluster water with moderate hydrogen bonds diffuses faster than bulk water with strong hydrogen bonds. The relatively free water and free water are formed during the diffusion process due to their interactions with the films matrix, which makes water molecules confined and restricted in limited space.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation 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 with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

A first principles molecular dynamics study of lithium atom solvation in binary liquid mixture of water and ammonia: structural, electronic, and dynamical properties

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.     

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.     

Development of non‐bonded interaction parameters between graphene and water using particle swarm optimization

New Lennard‐Jones parameters have been developed to describe the interactions between atomistic model of graphene, represented by REBO potential, and five commonly used all‐atom water models, namely SPC, SPC/E, SPC/Fw, SPC/Fd, and TIP3P/Fs by employing particle swarm optimization (PSO) method. These new parameters were optimized to reproduce the macroscopic contact angle of water on a graphene sheet. The calculated line tension was in the order of 10⁻¹¹ J/m for the droplets of all water models. Our molecular dynamics simulations indicate the preferential orientation of water molecules near graphene–water interface with one O H bond pointing toward the graphene surface. Detailed analysis of simulation trajectories reveals the presence of water molecules with ≤∼1, ∼2, and ∼4 hydrogen bonds at the surface of air–water interface, graphene–water interface, and bulk region of the water droplet, respectively. Presence of water molecules with ≤∼1 and ∼2 hydrogen bonds suggest the existence of water clusters of different sizes at these interfaces. The trends observed in the libration, bending, and stretching bands of the vibrational spectra are closely associated with these structural features of water. The inhomogeneity in hydrogen bond network of water at the air–water and graphene–water interface is manifested by broadening of the peaks in the libration band for water present at these interfaces. The stretching band for the molecules in water droplet shows a blue shift as compared to the pure bulk water, which conjecture the presence of weaker hydrogen bond network in a droplet. © 2017 Wiley Periodicals, Inc.     

Absorption and Structural Property of Ethanol/Water Mixture with Carbon Nanotubes

Molecular dynamics (MD) simulations are performed to study the structure and adsorption of ethanol/water mixture within carbon nanotubes (CNTs). Inside the (6,6) and (10,10) CNTs, there are always almost full of ethanol molecules and hardly water molecules. Inside wider CNTs, there are some water molecules, while the ethanol mass fractions inside the CNTs are still much higher than the corresponding bulk values. A series of structural analysis for the molecules inside and outside the CNTs are performed, including the distributions of radial, axial, angular density, orientation, and the number of hydrogen bonds. The angular density distribution of the molecules in the first solvation shell outside the CNTs indicates that the methyl groups of ethanol molecules have the strongest interaction with the carbon wall, and are pinned to the centers of the hexagons of the CNTs. Based on the understanding of the microscopic mechanism of these phenomena, we propose that the CNTs prefer to contain ethanol rather than methanol.     

On the Problem of Criteria for Phase Transitions in Water Clusters (A Hexamer and Octamer Example)

The study concerns the dynamics of a network of hydrogen bonds in small water clusters (on the example of hexamer and octamer clusters). The numerical experiment is carried out with the molecular dynamics method using two analytical rigid potentials TIP4P and TIP5P. According to the obtained results, the criteria of phase transitions provided by the cluster theory are not reliable enough in the case of small water clusters. This is testified by the distributions of potential energies of individual molecules in the clusters. It is shown that the appearance of several peaks in these distributions reflects the number of hydrogen bonds between the molecule and other molecules of the cluster.     

Molecular dynamics simulation study of association in trifluoroethanol/water mixtures

Mixtures of Trifluoroethanol (TFE) and water with different proportions are studied using molecular dynamics simulations. The radial and spatial distribution functions, as well as the size distribution of TFE clusters are obtained from the trajectories. The variation of radial and spatial distribution functions with composition show that the addition of TFE enhances the water structure, but the hydrogen bonds between TFE molecules are broken as TFE is diluted with water. The TFE‐rich solutions have stronger TFE–water hydrogen bonds. The clustering of TFE molecules in low concentration region is attributed to the hydrophobic interactions between CF₃ groups. The distribution of cluster sizes in solution supports these conclusions. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

水中受拉伸负载下单壁纳米碳管机械性质的分子动力学研究

利用分子动力学方法研究了(5,5)扶手椅型和(10,10)锯齿型纳米碳管在水中受拉伸负载下的机械性质.通过计算纳米碳管中氧和氢原子的局部密度分布研究了限制效应.结果表明,碳管在水中的杨式系数与在真空下相同,而碳管在水中的拉伸应力小于在真空中的.     

受限于纳米碳管中的乙醇分子的结构和扩散的分子模拟研究

利用分子模拟研究了常温常压下受限于(8,8) (管径1.081 nm)和(15,15) (管径 2.035 nm)单壁纳米碳管中的乙醇分子. 对受限分子的径向密度分布和氢键等静态性质以及扩散性质进行了分析. 结果显示在管内乙醇分子的平均氢键数目和主体相一致. 乙醇分子在(8,8)碳管内具有高度有序的结构, 而在(15,15)碳管内由于空间的增大导致结构有序度的降低, 其中分子取向已呈随机分布. 进一步对扩散系数的分析发现, 在管内乙醇分子的轴向扩散系数低于主体相, 特别在(8,8)碳管内乙醇分子几乎丧失了轴向扩散能力.     

Molecular dynamics simulations of atomistic hydration structures of poly(vinyl methyl ether)

Molecular dynamics simulations have been performed on the aqueous solutions of poly(vinyl methyl ether) (PVME) at various concentrations. Both radial and spatial distribution functions are used to investigate the detailed hydration structures. The structures of water are found to get increasingly concentrated when polymers are introduced and the water motions are severely hindered by the polymer matrix. At low concentrations, larger populations of tt conformers in meso dyads than those at higher concentrationsare found and this phenomenon is believed to be due to the increasing in bonding of water molecule to two ether oxygens in meso dyad. At higher concentrations, the size and conformations of polymers are quite similar to those in bulk. A transition of hydrogen bond fractions between PVME and water at around the concentration of 0.3 is observed and this value is perfectly in agreement with the results of conformational analysis and Raman spectra. Second neighbor hydrogen bond statistics revealed the domination of complicated hydrogen bond networks at low concentrations, but single hydrogen bonds as well as isolated clusters composed of 2-4 water molecules are usual around each polymer repeat unit.     

In situ study of diffusion and interaction of water and mono- or divalent anions in a positively charged membrane using two-dimensional correlation FT-IR/attenuated total reflection spectroscopy

Two-dimensional (2D) correlation analysis based on time-resolved FT-IR/attenuated total reflection (ATR) spectroscopy has been used to study the diffusion behavior of water and mono- or divalent anions in the positively charged membranes of different charge density. In 2D FT-IR/ATR spectra, the splitting of the water delta(OH) bending band in the spectral range 1700-1500 cm-1 indicates that there are three different states of water in the positively charged membrane, that is, the water molecules forming strong or weak hydrogen bonds with hydrophilic groups of the membrane and water molecularly dispersed with weaker hydrogen bonds. The wavenumber difference of the delta(OH) band in the low- and high-charge-density membrane indicates that water molecules form much stronger hydrogen bonds with hydrophilic groups in the high-charge-density membrane. The sequential order of the three water bands intensity changes shows that, in the process of water diffusion into the high-charge-density membrane, the hydrogen-bonding interaction between hydrophilic groups of the membrane and water molecules takes place gradually due to the highly cross-linked network structure of the membrane; in the process of water diffusion into the low-charge-density membrane, the strong hydrogen-bonding interaction between hydrophilic groups of the membrane and water molecules takes place instantaneously and this type of water easily diffuses due to the weak interactions between the water molecules and the membrane polymer. Furthermore, the diffusion processes of the electrolyte solution such as NaAc and Na2SO4 aqueous solutions in the positively charged membrane have also been examined.     

Long-range influence of carbohydrates on the solvation dynamics of water--answers from terahertz absorption measurements and molecular modeling simulations

We present new terahertz (THz) spectroscopic measurements of solvated sugars and compare the effect of two disaccharides (trehalose and lactose) and one monosaccharide (glucose) with respect to the solute-induced changes in the sub-picosecond network dynamics of the hydration water. We found that the solute affects the fast collective network motions of the solvent, even beyond the first solvation layer. For all three carbohydrates, we find an increase of 2-4% in the THz absorption coefficient of the hydration water in comparison to bulk water. Concentration-dependent changes in the THz absorption between 2.1 and 2.8 THz of the solute-water mixture were measured with a precision better than 1% and were used to deduce a dynamical hydration shell, which extends from the surface up to 5.7 +/- 0.4 and 6.5 +/- 0.9 A for the disaccharides lactose and trehalose, respectively, and 3.7 +/- 0.9 A for the glucose. This exceeds the values for the static hydration shell as determined, for example, by scattering, where the long-range structure was found to be not significantly affected by the solute beyond the first hydration shell. When comparing all three carbohydrates, we found that the solute-induced change in the THz absorption depends on the product of molar concentration of the solute and the number of hydrogen bonds between the carbohydrate and water molecules. We can conclude that the long-range influence on the sub-picosecond collective water network motions of the hydration water is directly correlated with the average number of hydrogen bonds between the molecule and adjacent water molecules for carbohydrates. This implies that monosaccharides have a smaller influence on the surrounding water molecules than disaccharides. This could explain the bioprotection mechanism of sugar-water mixtures, which has been found to be more effective for disaccharides than for monosaccharides.