Hydration force between mica surfaces in aqueous KCl electrolyte solution

Liquid-vapor molecular dynamics simulations are performed to study the interaction forces between two mica surfaces in an aqueous KCl electrolyte solution. Strong repulsive hydration force is obtained within a distance of ~2 nm between the two mica surfaces, which cannot be explained by the continuum theory of double-layer repulsion. We find that this short-range repulsive hydration force is much stronger than the double-layer force between mica surfaces. Whereas the simulation system is much smaller than the surface force measurement system, fundamental mechanisms of repulsive hydration force are revealed. In particular, important features of the step-like force oscillatory behavior during normal compression and force hysteresis during retraction are observed. Detailed analysis of the ionic density distributions shows that the "forced adsorption" of diffusive K(+) ions onto mica surfaces during compression and the subsequent "slow desorption" of the absorbed K(+) ions from mica surfaces upon retraction are responsible for the hysteresis phenomenon. From a mechanics point of view, we attribute the load bearing capacity of the dense electrolyte to the very hard hydration shells of K(+) metal ions under confinement. We find that the hydrated K(+) ions and Cl(-) co-ions remain very diffusive in the aqueous film. Water molecules in the hydration layer are also very fluidic, in the sense that the diffusion constant of water molecules is less than its bulk value by at most 3 orders of magnitude under the extreme confinement.     

Structure and dynamics of water in nonionic reverse micelles: A combined time-resolved infrared and small angle x-ray scattering study

We study the structure and reorientation dynamics of nanometer-sized water droplets inside nonionic reverse micelles (water/Igepal-CO-520/cyclohexane) with time-resolved mid-infrared pump-probe spectroscopy and small angle x-ray scattering. In the time-resolved experiments, we probe the vibrational and orientational dynamics of the O-D bonds of dilute HDO:H(2)O mixtures in Igepal reverse micelles as a function of temperature and micelle size. We find that even small micelles contain a large fraction of water that reorients at the same rate as water in the bulk, which indicates that the polyethylene oxide chains of the surfactant do not penetrate into the water volume. We also observe that the confinement affects the reorientation dynamics of only the first hydration layer. From the temperature dependent surface-water dynamics, we estimate an activation enthalpy for reorientation of 45 ± 9 kJ mol(-1) (11?±?2 kcal mol(-1)), which is close to the activation energy of the reorientation of water molecules in ice.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

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.     

Vibrational density of states of hydration water at biomolecular sites: hydrophobicity promotes low density amorphous ice behavior

Inelastic neutron scattering experiments and molecular dynamics simulations have been used to investigate the low frequency modes, in the region between 0 and 100 meV, of hydration water in selected hydrophilic and hydrophobic biomolecules. The results show changes in the plasticity of the hydrogen-bond network of hydration water molecules depending on the biomolecular site. At 200 K, the measured low frequency density of states of hydration water molecules of hydrophilic peptides is remarkably similar to that of high density amorphous ice, whereas, for hydrophobic biomolecules, it is comparable to that of low density amorphous ice behavior. In both hydrophilic and hydrophobic biomolecules, the high frequency modes show a blue shift of the libration mode as compared to the room temperature data. These results can be related to the density of water molecules around the biological interface, suggesting that the apparent local density of water is larger in a hydrophilic environment.     

Formation of ice-like water structure on the surface of an antifreeze protein

Antifreeze proteins (AFPs) are found in different species from polar, alpine, and subarctic regions where they serve to inhibit ice crystal growth by adsorption to ice surfaces. Computational methods have the power to investigate the antifreeze mechanism in atomic detail. Molecular dynamics simulations of water under different conditions have been carried out to test our water model for simulations of biological macromolecules in extreme conditions: very low temperatures (200 K) and at the ice/liquid water interface. We show that the flexible F3C water model reproduces properties of water in the solid phase (ice I(h)), the supercooled liquid phase, and at the ice/liquid water interface. Additionally, the hydration of the type III AFP from ocean pout was studied as a function of temperature. Hydration waters on the ice-binding surface of the AFP were less distorted and more tetrahedral than elsewhere on the surface. More ice-like hydrating water structures formed on the ice-binding surface of the protein such that it created an ice-like structure in water within its first hydration layer but not beyond, suggesting that this portion of the protein has high affinity for ice surfaces.     

Novel ice structures in carbon nanopores: pressure enhancement effect of confinement

We report experimental results on the structure and melting behavior of ice confined in multi-walled carbon nanotubes and ordered mesoporous carbon CMK-3, which is the carbon replica of a SBA-15 silica template. The silica template has cylindrical mesopores with micropores connecting the walls of neighboring mesopores. The structure of the carbon replica material CMK-3 consists of carbon rods connected by smaller side-branches, with quasi-cylindrical mesopores of average pore size 4.9 nm and micropores of 0.6 nm. Neutron diffraction and differential scanning calorimetry have been used to determine the structure of the confined ice and the solid-liquid transition temperature. The results are compared with the behavior of water in multi-walled carbon nanotubes of inner diameters of 2.4 nm and 4 nm studied by the same methods. For D(2)O in CMK-3 we find evidence of the existence of nanocrystals of cubic ice and ice IX; the diffraction results also suggest the presence of ice VIII, although this is less conclusive. We find evidence of cubic ice in the case of the carbon nanotubes. For bulk water these crystal forms only occur at temperatures below 170 K in the case of cubic ice, and at pressures of hundreds or thousands of MPa in the case of ice VIII and IX. These phases appear to be stabilized by the confinement.     

Melting behavior of water in cylindrical pores: carbon nanotubes and silica glasses

We report a study of the effects of confinement in multi-walled carbon nanotubes and mesoporous silica glasses (SBA-15) on the solid structure and melting of both H(2)O and D(2)O ice, using differential scanning calorimetry, dielectric relaxation spectroscopy, and neutron diffraction. Multi-walled nanotubes of 2.4, 3.9 and 10 nm are studied, and the SBA-15 studied has pores of mean diameter 3.9 nm; temperatures ranging from approximately 110 to 290 K were studied. We find that the melting point is depressed relative to the bulk water for all systems studied, with the depression being greater in the case of the silica mesopores. These results are shown to be consistent with molecular simulation studies of freezing in silica and carbon materials. The neutron diffraction data show that the cubic phase of ice is stabilized by the confinement in carbon nanotubes, as well as in silica mesopores, and persists up to temperatures of about 240 K, above which there is a transition to the hexagonal ice structure.     

Effect of solute size and solute-water attractive interactions on hydration water structure around hydrophobic solutes.

Using Monte Carlo simulations, we investigated the influence of solute size and solute-water attractive interactions on hydration water structure around spherical clusters of 1, 13, 57, 135, and 305 hexagonally close-packed methanes and the single hard-sphere (HS) solute analogues of these clusters. We obtain quantitative results on the density of water molecules in contact with the HS solutes as a function of solute size for HS radii between 3.25 and 16.45 A. Analysis of these results based on scaled-particle theory yields a hydration free energy/surface area coefficient equal to 139 cal/(mol A2), independent of solute size, when this coefficient is defined with respect to the van der Waals surface of the solute. The same coefficient defined with respect to the solvent-accessible surface decreases with decreasing solute size for HS radii less than approximately 10 A. We also find that solute-water attractive interactions play an important role in the hydration of the methane clusters. Water densities in the first hydration shell of the three largest clusters are greater than bulk water density and are insensitive to the cluster size. In contrast, contact water densities for the HS analogues of these clusters decrease with solute size, falling below the bulk density of water for the two largest solutes. Thus, the large HS solutes dewet, while methane clusters of the same size do not.     

Role of electromechanical and mechanoelectric effects in protein hydration under hydrostatic pressure

Recent measurements of lysozyme hydration water density under non-denaturing pressure show that it is higher than that of bulk water in the same conditions. High protein hydration layer density has earlier been observed at ambient conditions and ascribed to electrostriction. We calculate the pressure-induced protein mean surface charge density increment Δσ. Within the hydration layer, the higher fields due to Δσ lead to an additional water compression via electrostriction. The increment Δσ is considered as due to a mechanoelectric effect in protein molecules. The mean value of the effective mechanoelectric coefficient d is calculated and compared with piezoelectric coefficients of amino acids and their compounds.     

Ion permeation dynamics in carbon nanotubes

Molecular dynamics simulations are carried out to investigate the permeation of ions and water in a membrane consisting of single wall carbon nanotubes possessing no surface charges connecting two reservoirs. Our simulations reveal that there are changes in the first hydration shell of the ions upon confinement in tubes of 0.82 or 0.90 nm effective internal diameter. Although the first minimum in the g(r) is barely changed in the nanotube compared to in the bulk solution, the hydration number of Na(+) ion is reduced by 1.0 (from 4.5 in bulk to 3.5 in the 0.90 nm tube) and the hydration number is reduced further in the 0.82 nm tube. The changes in the hydration shell of Cl(-) ion are negligible, within statistical errors. The water molecules of the first hydration shell of both ions exchange less frequently inside the tube than in the bulk solution. We compare ion trajectories for ions in the same tube under identical reservoir conditions but with different numbers of ions in the tubes. This permits investigation of changes in structure and dynamics which arise from multiple ion occupancy in a carbon nanotube possessing no surface charges. We also investigated the effects of tube flexibility. Ions enter the tubes so as to form a train of ion pairs. We find that the radial distribution profiles of Na(+) ions broaden significantly systematically with increasing number of ion pairs in the tube. The radial distribution profiles of Cl(-) ions change only slightly with increasing number of ions in the tube. Trajectories reveal that Na(+) ions do not pass each other in 0.90 nm tubes, while Cl(-) ions pass each other, as do ions of opposite charge. An ion entering the tube causes the like-charged ions preceding it in the tube to be displaced along the tube axis and positive or negative ions will exit the tube only when one or two other ions of the same charge are present in the tube. Thus, the permeation mechanism involves multiple ions and Coulomb repulsion among the ions plays an essential role.     

Ab Initio H2O in Realistic Hydrophilic Confinement

A protocol for the ab initio construction of a realistic cylindrical pore in amorphous silica, serving as a geometric nanoscale confinement for liquids and solutions, is presented. Upon filling the pore with liquid water at different densities, the structure and dynamics of the liquid inside the confinement can be characterized. At high density, the pore introduces long‐range oscillations into the water density profile, which makes the water structure unlike that of the bulk across the entire pore. The tetrahedral structure of water is also affected up to the second solvation shell of the pore wall. Furthermore, the effects of the confinement on hydrogen bonding and diffusion, resulting in a weakening and distortion of the water structure at the pore walls and a slowdown in diffusion, are characterized.     

Negative thermal expansion of water in hydrophobic nanospaces

The density and intermolecular structure of water in carbon micropores (w = 1.36 nm) are investigated by small-angle X-ray scattering (SAXS) and X-ray diffraction (XRD) measurements between 20 K and 298 K. The SAXS results suggest that the density of the water in the micropores increased with increasing temperature over a wide temperature range (20-277 K). The density changed by 10%, which is comparable to the density change of 7% between bulk ice (I(c)) at 20 K and water at 277 K. The results of XRD at low temperatures (less than 200 K) show that the water forms the cubic ice (I(c)) structure, although its peak shape and radial distribution functions changed continuously to those of a liquid-like structure with increasing temperature. The SAXS and XRD results both showed that the water in the hydrophobic nanospaces had no phase transition point. The continuous structural change from ice I(c) to liquid with increasing temperature suggests that water shows negative thermal expansion over a wide temperature range in hydrophobic nanospaces. The combination of XRD and SAXS measurements makes it possible to describe confined systems in nanospaces with intermolecular structure and density of adsorbed molecular assemblies.     

Dynamics of surface water in ZrO2 studied by quasielastic neutron scattering

A quasielastic neutron scattering experiment has revealed the dynamics of surface water in a high surface area zirconium oxide in the temperature range of 300-360 K. The characteristic times of the rotational (picoseconds) and translational (tens of picoseconds) components of diffusion motion are well separated. The rotational correlation time shows an Arrhenius-type behavior with an activation energy of 4.48 kJ/mol, which is lower compared to bulk water. The rotational diffusion at room temperature is slower by about a factor of 2 compared to bulk water, whereas the translational diffusion slows down by a factor of 40. In contrast to bulk water, the translational correlation time exhibits an Arrhenius-type temperature dependence with an activation energy of 11.38 kJ/mol. Comparison of different models for jump diffusion processes suggests that water molecules perform two-dimensional jumps at a well-defined, almost temperature-independent distance of 4.21-4.32 A. Such a large jump distance indicates a low molecular density of the layer of diffusing molecules. We argue that undissociated water molecules on an average form two hydrations layers on top of the surface layer of hydroxyl groups, and all the layers have similar molecular density. Quasielastic neutron scattering experiment assesses the dynamics of the outermost hydration layer, whereas slower motion of the water molecules in the inner hydration layer contributes to the elastic signal.     

Hydration of NaCl on glassy, supercooled-liquid, and crystalline water

Interactions of sodium chloride with amorphous and crystalline water films, leading to the possible formation of a dilute NaCl solution, were investigated using time-of-flight secondary ion mass spectrometry as a function of temperature. A monolayer of NaCl tends to remain on the surface or in subsurface sites of thick amorphous solid water films (200 monolayers); the Na+ ion is hydrated preferentially, whereas the Cl- ion is segregated at the surface. The hydration structure of NaCl is fundamentally unchanged for viscous liquid water that appears at temperatures higher than 136 K. The solubility of NaCl increases abruptly at 160 K because of the evolution of supercooled liquid water, which can hydrate the Cl- ion efficiently. However, the diffusion of the ions toward the bulk of supercooled liquid water is interrupted by crystallization; therefore, the dilute NaCl solution that is characterized by completely separated Na+-Cl- pairs may not be formed. When NaCl is deposited on the crystalline ice film, hydration of NaCl is enhanced above 160 K as well, indicating that a liquidlike phase coexists with crystals.     

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

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

Molecular simulation studies of nanoscale friction between phosphorylcholine self-assembled monolayer surfaces: correlation between surface hydration and friction

We performed all-atom molecular dynamics simulations to study the friction between surfaces covered with two phosphorylcholine self-assembled monolayers (PC-SAM) under shear. PC-SAM surfaces with a sqrt7 X sqrt7R19 degrees lattice structure and a parallel arrangement of the head groups were used as model zwitterionic surfaces. They provide a full representation of the zwitterionic nature of phospholipid surfaces, which are believed to play an important role in the lubrication of biological joints such as knees and hips. The surfaces were immersed in aqueous solutions and kept in contact with two regions of bulk water. Sodium chloride and potassium chloride solutions at various concentrations were employed to study the effects of the presence of ions on friction. The results show a strong relationship between surface hydration and friction. Higher ionic concentrations or ions with shorter Debye lengths cause a larger disruption to the hydration around the zwitterionic surfaces, leading to larger friction forces. In addition, the results show that under nanoscale confinement, the friction coefficients of PC-SAM surfaces in pure water are directly proportional to both shear velocity and surface separation distance. These results are comparable to previously published experimental studies.     

Hydration water and bulk water in proteins have distinct properties in radial distributions calculated from 105 atomic resolution crystal structures

Water plays a critical role in the structure and function of proteins, although the experimental properties of water around protein structures are not well understood. The water can be classified by the separation from the protein surface into bulk water and hydration water. Hydration water interacts closely with the protein and contributes to protein folding, stability, and dynamics, as well as interacting with the bulk water. Water potential functions are often parametrized to fit bulk water properties because of the limited experimental data for hydration water. Therefore, the structural and energetic properties of the hydration water were assessed for 105 atomic resolution (相似文献   

Hydration layer of a cationic micelle, C(10)TAB: structure, rigidity, slow reorientation, hydrogen bond lifetime, and solvation dynamics

We report a theoretical study of the structure and dynamics of the water layer (the hydration layer) present at the surface of the cationic micelle decyltrimethylammonium bromide (DeTAB) by using atomistic molecular dynamics simulations. The simulated micelle consisted of 47 surfactant molecules (and an equal number of bromide ions), in good agreement with the pioneering light scattering experiments by Debye which found an aggregation number of 50. In this micelle, three partially positively charged methyl groups of each surfactant headgroup face the surrounding water. The nature of the cationic micellar surface is found to play an important role in determining the arrangement of water which is quite different from that in the bulk or on the surface of an anionic micelle, like cesium perfluorooctanoate. Water molecules present in the hydration layer are found to be preferentially distributed in the region between the three partially charged methyl headgroups. It is found that both the translational and rotational motions of water exhibit appreciably slower dynamics in the layer than those in the bulk. The solvation time correlation function (TCF) of bromide ions exhibits a long time component which is found to originate primarily from the interaction of the probe with the micellar headgroups. Thus, the decay of the solvation TCF is controlled largely by the residence time of the probe in the surface. The residence time distribution of the water molecules also exhibits a slow time component. We also calculate the collective number density fluctuation in the layer and find a prominent slow component compared to the similar quantity in the bulk. This slow component demonstrates that water structure in the hydration layer is more rigid than that in the bulk. These results demonstrate that the slow dynamics of hydration layer water is generic to macromolecular surfaces of either polarity.     

Structure, phase behavior, and inhomogeneous fluid properties of binary dendrimer mixtures

The effective pair potentials between different kinds of dendrimers in solution can be well approximated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and strength of the effective interactions depend strongly upon the specific dendrimer architecture. We consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, to determine the bulk fluid structure and phase behavior. Using a simple mean field density functional theory (DFT) we find good agreement between theory and simulation results for the bulk fluid structure. Depending on the mixture, we find bulk fluid-fluid phase separation (macrophase separation) or microphase separation, i.e., a transition to a state characterized by undamped periodic concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in spherical cavities. Again, we find good agreement between the DFT and simulation results. For the dendrimer mixture exhibiting microphase separation, we observe a rather striking pattern formation under confinement.