Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Treatment of dilute clusters of methanol and water by ab initio quantum mechanical calculations

Large molecular clusters can be considered as intermediate states between gas and condensed phases, and information about them can help us understand condensed phases. In this paper, ab initio quantum mechanical methods have been used to examine clusters formed of methanol and water molecules. The main goal was to obtain information about the intermolecular interactions and the structure of methanol/water clusters at the molecular level. The large clusters (CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10)) containing one molecule of one component (methanol or water) and many (12, 10) molecules of the other component were considered. M?ller-Plesset perturbation theory (MP2) was used in the calculations. Several representative cluster geometries were optimized, and nearest-neighbor interaction energies were calculated for the geometries obtained in the first step. The results of the calculations were compared to the available experimental information regarding the liquid methanol/water mixtures and to the molecular dynamics and Monte Carlo simulations, and good agreement was found. For the CH(4)O...(H(2)O)(12) cluster, it was shown that the molecules of water can be subdivided into two classes: (i) H bonded to the central methanol molecule and (ii) not H bonded to the central methanol molecule. As expected, these two classes exhibited striking energy differences. Although they are located almost the same distance from the carbon atom of the central methanol molecule, they possess very different intermolecular interaction energies with the central molecule. The H bonding constitutes a dominant factor in the hydration of methanol in dilute aqueous solutions. For the H(2)O...(CH(4)O)(10) cluster, it was shown that the central molecule of water has almost three H bonds with the methanol molecules; this result differs from those in the literature that concluded that the average number of H bonds between a central water molecule and methanol molecules in dilute solutions of water in methanol is about two, with the water molecules being incorporated into the chains of methanol. In contrast, the present predictions revealed that the central water molecule is not incorporated into a chain of methanol molecules, but it can be the center of several (2-3) chains of methanol molecules. The molecules of methanol, which are not H bonded to the central water molecule, have characteristics similar to those of the methane molecules around a central water molecule in the H(2)O...(CH(4))(10) cluster. The ab initio quantum mechanical methods employed in this paper have provided detailed information about the H bonds in the clusters investigated. In particular, they provided full information about two types of H bonds between water and methanol molecules (in which the water or the methanol molecule is the proton donor), including information about their energies and lengths. The average numbers of the two types of H bonds in the CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10) clusters have been calculated. Such information could hardly be obtained with the simulation methods.     

Dielectric Properties of Solutions Isomorphous with Water

The dielectric properties of aqueous solutions of ammonia, hydrogen fluoride, hydrogen peroxide, and ammonium halides are determined by the number of hydrogen bonds formed by the solute species with the molecules forming the water framework. For such solutions, the product of the concentration coefficient of the dielectric permittivity by the ratio of the number of hydrogen bonds per solute molecule and per water molecule is constant. In aqueous solutions with a water-like structure, the mechanism of the dielectric polarization is the same as in pure water (proton-activation mechanism) and is determined by proton fluctuations on the line of hydrogen bonds.__________Translated from Zhurnal Obshchei Khimii, Vol. 75, No. 1, 2005, pp. 39–45.Original Russian Text Copyright © 2005 by Potapov, Parkhomenko.     

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.     

Dependence of the number of hydrogen bonds per water molecule on its distance to a hydrophobic surface and a thereupon-based model for hydrophobic attraction

A water molecule in the vicinity of a hydrophobic surface forms fewer hydrogen bonds than a bulk molecule because the surface restricts the space available for other water molecules necessary for its hydrogen-bonding. In this vicinity, the number of hydrogen bonds per water molecule depends on its distance to the surface. Considering the number of hydrogen bonds per bulk water molecule (available experimentally) as the only reference quantity, we propose an improved probabilistic approach to water hydrogen-bonding that allows one to obtain an analytic expression for this dependence. (The original version of this approach [Y. S. Djikaev and E. Ruckenstein, J. Chem. Phys. 130, 124713 (2009)] provides the number of hydrogen bonds per water molecule in the vicinity of a hydrophobic surface as an average over all possible locations and orientations of the molecule.) This function (the number of hydrogen bonds per water molecule versus its distance to a hydrophobic surface) can be used to develop analytic models for the effect of hydrogen-bonding on the hydration of hydrophobic particles and their solvent-mediated interaction. Presenting a model for the latter, we also examine the temperature effect on the solvent-mediated interaction of two parallel hydrophobic plates.     

Study of molecular behavior in a water nanocluster: size and temperature effect

Temperature and size effects on the behavior of nanoscale water molecule clusters are investigated by molecular dynamics simulations. The flexible three-centered (F3C) water potential is used to model the inter- and intramolecular interactions of the water molecule. The differences between the structural properties for the surface region and those for the interior region of the cluster are also investigated. It is found that as the temperature rises, the average number of hydrogen bonds per water molecule decreases, but the ratio of surface water molecules increases. After comparing the water densities in interior regions and the average number of hydrogen bonds in those regions, we find there is no apparent size effect on water molecules in the interior region, whereas the size of the water cluster has a significant influence on the behavior of water molecules at the surface region.     

Very strong hydrogen bonds in a bent chain structure of fluorohydrogenate anions in liquid Cs(FH)2.3F

The structure of liquid Cs(FH)(2.3)F was revealed using a combination of high-energy x-ray and neutron diffraction measurements. We found that the strongest intermolecular H-F hydrogen bonds at an average distance of 1.36 A are accompanied by the formation of a high degree of bending of the oligomer chain in the melt, with [angle]FHF=150 degrees . A reverse Monte Carlo simulation showed that the average number of atoms per chain is 4.4. A detailed chain analysis of the atomic configuration revealed that (FH)(2)F(-) oligomer chains are the major entities in the liquid, and asymmetrical FHF(-) are formed owing to the strong H-F hydrogen bonds. The results suggest that an average of one or two HF molecules bond to each of the 11 fluorine atoms surrounding a cesium ion.     

Structure of methanol confined in MCM-41 investigated by large-angle X-ray scattering technique

Large-angle X-ray scattering (LAXS) measurements over a temperature range from 223 to 298 K have been made on methanol confined in mesoporous silica MCM-41 with two different pore diameters, 28 A (C14) and 21 A (C10), under both monolayer and capillary-condensed adsorption conditions. To compare the structure of methanol in the MCM-41 pores with that of bulk methanol, X-ray scattering intensities for bulk methanol in the same temperature range have also been measured. The radial distribution functions (RDFs) for the monolayer methanol samples showed that methanol molecules are strongly hydrogen bonded to the silanol groups on the MCM-41 surface, resulting in no significant change in the structure of adsorbed methanol with respect to the pore size and temperature. On the other hand, the RDFs for the capillary-condensed methanol samples showed that hydrogen-bonded chains of methanol molecules are formed in both pores. However, the distance and number of hydrogen bonds estimated from the RDFs suggested that hydrogen bonds between methanol molecules in the pores are significantly distorted or partly disrupted. It has been found that the hydrogen bonds are more distorted in the smaller pores of MCM-41. With decreasing temperature, however, the hydrogen-bonded chains of methanol in the pores were gradually ordered. A comparison of the present results on methanol in MCM-41 pores with those on water in the same pores revealed that the structural change with temperature is less significant for confined methanol than for confined water.     

The structure of liquid alcohols and the temperature dependence of vibrational bandwidth

To obtain the information about the structure peculiarities of liquid alcohols the temperature dependence of Raman profiles widths for these objects was studied. The results of these investigations have shown that in the region 150–340 K the widths of Raman bands in liquid methanol and ethanol are constant. From the data on the constancy of the number of hydrogen bonds per one molecule (in the same temperature range) the conclusion on the dynamic stability of the cluster structure of alcohols was made. The broadening of vibrational bands is determined by inhomogeneous broadening and by the dephasing of intramolecular vibrations due to the hydrogen bond dissociation.     

Hydrogen-bonded assembly of methanol on Cu(111)

Investigation of methanol's surface chemistry on metals is a crucial step towards understanding the reactivity of this important chemical feedstock. Cu is a relevant metal for methanol synthesis and reforming, but due to the weak interaction of methanol with Cu, an atomic scale view of methanol's coverage-dependent ordering and self-assembly on Cu(111), the most abundant facet of most nanoparticles, has not yet been possible. Low and variable temperature scanning tunneling microscopy coupled with density functional theory reveal a coverage-dependent range of highly ordered structures stabilized by two hydrogen bonds per molecule. While extended chains that resemble the hydrogen-bonded zigzag structures reported for solid methanol are an efficient way to pack methanol at higher coverages, lower surface coverages yield isolated hexamer units. These hexamers form the same number of hydrogen bonds as the chains but appear to repel one another on the surface. Annealing treatments lead to the desorption of methanol with almost no decomposition. This data serves as a useful guide to both the preferred adsorption geometries and energies of a variety of methanol structures on Cu(111) surfaces as a function of surface coverage.     

High-density hydration layer of lysozymes: molecular dynamics decomposition of solution scattering data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.     

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

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

Calculating the radial distribution functions of supercritical methanol by means of Car-Parrinello and classical molecular dynamics

Radial distribution functions and the average number of hydrogen bonds per methanol molecule under standard, subcritical, and supercritical conditions are obtained via classical molecular dynamics and Car-Parrinello nonempirical molecular dynamics. It is shown that independent methods of modeling yield close results. It is noted that the calculated radial distribution functions agree well with the experimental data only at T = 298 K and P = 0.1 MPa, while at high temperatures and pressures, considerable divergence from the experimental functions known from the literature is observed. It is concluded that both modeling methods reproduce the degree of hydrogen bonding in methanol and its variations depending on the state parameters and correspond closely to the experimental results.     

水、甲醇和乙醇液体微结构性质的Car-Parrinello 分子动力学模拟

采用Car-Parrinello 分子动力学(CPMD)方法分别研究了水、甲醇和乙醇的液体微结构性质.研究结果显示:在水、甲醇和乙醇三个体系中O…O径向分布函数曲线的第一个峰位置分别为0.278、0.276 和0.275nm; O…H径向分布函数曲线的第一个峰位置分别为0.178、0.176和0.177 nm.表明基团(氢原子、甲基、乙基)的差异对O…O第一个峰的位置影响很小.但基团的差异对径向分布函数峰高的影响却很显著,由水到乙醇第一个峰的高度逐渐变高.空间分布函数表明氧原子和氢原子在溶剂分子周围有取向地分布,这与径向分布函数所表现出尖锐的第一个峰相一致.氢键分布分析显示,水、甲醇和乙醇的平均氢键数分别为3.62、1.99 和1.87,表明水形成了网状氢键结构,而甲醇、乙醇形成链状氢键结构.     

The structure of glycerol in the liquid state: a neutron diffraction study

Neutron diffraction coupled with hydrogen/deuterium isotopic substitution has been used to investigate the structure of the pure cryoprotectant glycerol in the liquid state at 298 K and 1 atm. The neutron diffraction data were used to constrain a 3 dimensional computational model that is experimentally relevant using the empirical potential structure refinement (EPSR) technique. These simulations lead to a model structure of the glycerol molecule that is consistent with the experimental data. Interestingly, from interrogation of this structure, it is found that the number of hydrogen bonds per molecule is larger than had previously been suggested. Furthermore, converse to previous work, no evidence for intra-molecular hydrogen bonds is found. These results highlight the importance and relevance of using experimental data to inform computational modelling of even simple liquid systems.     

Influence of electric field on the hydrogen bond network of methanol

The understanding of the structure of hydrogen (H) bonding liquids in electric (E) fields is important in the context of several areas of research, such as electrochemistry, surface science, and thermodynamics of electrolyte solutions. We had earlier presented a general thermodynamic framework for this purpose, and had shown that the application of E field enhances H-bond interactions among water molecules. The present investigation with methanol suggests a different result-the H-bond structure, as indicated by the average number of H bonds per molecule, goes through a maxima with increasing field strength. This result is explained based on the symmetry in the location of the H-bonding sites in the two types of molecules.     

Whole-molecule calculation of log p based on molar volume, hydrogen bonds, and simulated 13C NMR spectra

The prediction of Log P is usually accomplished using either substructure or whole-molecule approaches. However, these methods are complicated, and previous whole-molecule approaches have not been successful for the prediction of Log P in very complex molecules. The observed chemical shifts in nuclear magnetic resonance (NMR) spectroscopy are related to the electrostatics at the nucleus, which are influenced by solute-solvent interactions. The different solvation effects on a molecule by either water or methanol have a strong effect on the NMR chemical shift value. Therefore, the chemical shift values observed in an aqueous and organic solvent should correlate to Log P. This paper develops a rapid, objective model of Log P based on molar volume, hydrogen bonds, and differences in calculated 13C NMR chemical shifts for a diverse set of compounds. A partial least squares (PLS) model of Log P built on the sum of carbon chemical shift differences in water and methanol, molar volume, number of hydrogen bond donors and acceptors in 162 diverse compounds gave an r2 value of 0.88. The average r2 for 10 training models of Log P made from 90% of the data was 0.87+/-0.01. The average q2 for 10 leave-10%-out cross-validation test sets was 0.87+/-0.05.     

甘油水溶液氢键特性的分子动力学模拟

为了研究低温保护剂溶液的结构和物理化学特性, 以甘油为保护剂, 采用分子动力学方法, 对不同浓度的甘油和水的二元体系进行了模拟. 得到了不同浓度的甘油水溶液在2 ns内的分子动力学运动轨迹, 通过对后1 ns内运动轨迹的分析, 得到了各个原子对的径向分布函数和甘油分子的构型分布. 根据氢键的图形定义, 分析了氢键的结构和动力学特性. 计算了不同浓度下体系中平均每个原子(O和H)和分子(甘油和水)参与氢键个数的百分比分布及其平均值. 同时还计算了所有氢键、水分子之间的氢键以及甘油与水分子之间的氢键的生存周期.     

Water-induced ethanol dewetting transition

The dewetting transitions of two hydrophobic plates immersed in pure water, aqueous ethanol solutions with concentrations from 25% to 90%, and pure ethanol were investigated by molecular dynamics simulations, where the dewetting transition was analogous to a first-order phase transition from liquid to vapor. It was found that the dewetting transitions occurred except that in the pure ethanol system. Although the ethanol molecules prefer to locate in the vicinity of the two plates, the inter-plate region is unfavorable for water molecules, due to losing more than one hydrogen bond. Moreover, each inter-plate water molecule forms hydrogen bonds on average with about two ethanol molecules. These intermolecular hydrogen bonds cause water and ethanol to cooperatively fill or exit the inter-plate region. Thus, water molecules play a more important role in the inter-plate filling/empty process, and induce the ethanol dewetting transition. Our results provide insight into the effect of water on the ethanol dewetting phenomena.     

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.