Thermodynamic implications of confinement for a waterlike fluid

Fast mass transport through carbon nanotube membranes

The May 19, 2006 issue of Science included a paper by Holt et al. on "Fast Mass Transport Through Sub-2-Nanometer Carbon Nanotubes". The paper was also featured on the cover, showing methane molecules translating inside a carbon nanotube (CNT). The authors explained how they prepared 2-6-mum thin membranes consisting of double-walled carbon nanotubes (DWNTs) all aligned perpendicular to the apparent membrane surface. These tubes are open at both ends and the space between the tubes is filled with dense Si(3)N(4). Pure gas and water fluxes were measured at room temperature with the application of a small pressure difference. Interpretation of the results led to the conclusion that the membranes showed much higher fluxes than what was estimated from Knudsen gas diffusion and Poiseuille viscous flow models. The membranes have a straight-channel morphology with a narrow pore-size distribution and exceptionally smooth pore walls. The unusual geometry and surface properties make it difficult to compare the membrane's properties with common membranes but there is no question that the mass transport in the aligned DWNTs is fast indeed. To appreciate how fast, we will consider their transport properties starting from the perspective of "conventional" porous membrane technology. Recent molecular dynamics simulations suggest that none of the classic models for gas (Knudsen) and water (Poiseuille) permeation work in a meaningful way for these nanotube membranes, and new models are needed.

Structural properties and liquid spinodal of water confined in a hydrophobic environment

We present the results of a computer simulation study of thermodynamical properties of TIP4P water confined in a hydrophobic disordered matrix of soft spheres upon supercooling. The hydrogen-bond network of water appears preserved in this hydrophobic confinement. Nonetheless a reduction in the average number of hydrogen bonds due to the geometrical constraints is observed. The liquid branch of the spinodal line is calculated from 350 K down to 210 K. The same thermodynamic scenario of the bulk is found: the spinodal curve is monotonically decreasing. The line of maximum density bends avoiding a crossing of the spinodal. There is, however, a shift both of the line of maximum density and of the spinodal toward higher pressures and lower temperatures with respect to bulk.

The effect of local defects on water adsorption in silicalite-1 zeolite: a joint experimental and molecular simulation study

We report a joint experimental and molecular simulation study of water condensation in silicalite-1 zeolite. A sample was synthesized using the fluoride route and was found to contain essentially no defects. A second sample synthesized using the hydroxide route was found to contain a small amount of silanol groups. The thermodynamics of water condensation was studied in these two samples, as well as in a commercial sample, in order to understand the effect of local defects on water adsorption. The molecular simulation study enabled us to qualitatively reproduce the experimentally observed condensation thermodynamics features. A shift and a rounding of the condensation transition was observed with an increasing hydrophilicity of the local defect, but the condensation transition was still observed above the water saturation vapor pressure P0. Both experiments and simulations agree on the fact that a small water uptake can be observed at very low pressure, but that the bulk liquid does not form from the gas phase below P0. The picture that emerges from the observed water condensation mechanism is the existence of a heterogeneous internal surface that is overall hydrophobic, despite the existence of hydrophilic "patches". This heterogeneous surface configuration is thermodynamically stable in a wide range of reduced pressures (from P/P0 = 0.2 to a few thousands), until the condensation transition takes place.

Large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction and nanobubble formation

We performed large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction manifested by the formation of nanobubble between nanometer-sized hydrophobic clusters at constrained equilibrium. Particular attention is placed on the tendency of formation and stability of nanobubbles in between model nanoassemblies which are composed of hydrophobic clusters (or patches) embedded in a hydrophilic substrate. On the basis of physical behavior of nanobubble formation, we observed a change from short-range molecular hydrophobic interaction to midrange nanoscopic interaction when the length scale of hydrophobe approaches to about 1 nm. We investigated the behavior of nanobubble formation with several different patterns of nonpolar-site distribution on the nanoassemblies but always keeping a constant ratio of nonpolar to polar monomer sites. Dynamical properties of confined water molecules in between nanoassemblies are also calculated.

Effect of Surface Polarity on Water Contact Angle and Interfacial Hydration Structure

We perform molecular dynamics simulations of water in the presence of hydrophobic/hydrophilic walls at T = 300 K and P = 0 GPa. For the hydrophilic walls, we use a hydroxylated silica model introduced in previous simulations [Lee, S. H.; Rossky, P. J. J. Chem. Phys. 1994, 100, 3334. Giovambattista, N.; Rossky, P. J.; Debenedetti, P. G.; Phys. Rev. E 2006, 73, 041604.]. By rescaling the physical partial atomic charges by a parameter 0 <or= k <or= 1, we can continuously transform the hydrophilic walls (hydroxylated silica, k = 1) into hydrophobic apolar surfaces (k = 0). From a physical point of view, k is the normalized magnitude of a surface dipole moment, and thus it quantifies the polarity of the surface. We calculate the contact angle of water for 0 <or= k <or= 1. We find that, at least for the present homogeneous, atomically flat, and defect-free surface model, the magnitude of the surface dipole correlates with the contact angle in a one-to-one correspondence. In particular, we find that polar surfaces with 0 < k <or= kc = 0.4 are macroscopically hydrophobic; that is, the contact angle is larger than 90 degrees . For the cutoff value k = kc, the magnitude of the dipole moment of the polar silica surface unit is 41% that of the water molecule dipole moment. We also study the water orientation distributions next to the walls (a microscopic property). We find that these distributions also correlate with the contact angle in a one-to-one correspondence. Thus, the structure of confined water, the surface polarity, and the contact angle are in a direct correspondence to each other, and therefore, each quantifies the hydrophobicity/hydrophilicity of the surface.

Structure of coexisting liquid phases of supercooled water: Analogy with ice polymorphs

The structural changes occurring in supercooled liquid water upon moving from one coexisting liquid phase to the other have been investigated by computer simulation using a polarizable interaction potential model. The obtained results favorably compare with recent neutron scattering data of high and low density water. In order to assess the physical origin of the observed structural changes, computer simulation of several ice polymorphs has also been carried out. Our results show that there is a strict analogy between the structure of various disordered (supercooled) and ordered (ice) phases of water, suggesting that the occurrence of several different phases of supercooled water is rooted in the same physical origin that is responsible for ice polymorphism.

Effect of confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.

The dewetting transition and the hydrophobic effect

A molecular-level description of the behavior of water in hydrophobic spaces is presented in terms of the coupled effects of solute size and atomic solute-solvent interactions. For model solutes with surface areas near those of protein contacts, we identify three different regions of solute-water interaction to be associated with three distinctly different structural characteristics of water in the intersolute region: dry, oscillating, and wet. A first orderlike phase transition is confirmed from the wet to dry state bridged by a narrow region with liquid-vapor oscillations in the intersolute region as the strength of the solute-water attractive dispersion interaction decreases. We demonstrate that the recent idea that cavitation in the intersolute region of nanoscopic solutes is preceded by the formation of a vapor layer around an individual solute is not the general case. The appearance of density waves pulled up around and outside of a nanoscopic plate occurs at lower interaction strengths than are required to obtain a wet state between such plates. We further show that chemically reasonable estimates of the interaction strength lead to a microscopically wet state and a hydrophobic interaction characterized by traps and barriers to association and not by vacuum induced collapse.

Analytical model for studying how environmental factors influence protein conformational stability in solution

We introduce an analytical modeling strategy for probing the conformational stability of globular proteins in aqueous solution. In this approach, the intrinsic (i.e., infinite dilution) thermodynamic stability and coarse structural properties of the proteins, as well as the effective protein-protein interactions, derive from a heteropolymer collapse theory that incorporates predicted temperature- and pressure-dependent hydrophobic interactions. Protein concentration effects are estimated by integrating this information into a molecular thermodynamic model, which is an ad hoc generalization of the exact equilibrium theory of a one-dimensional binary mixture of square-well particles that interconvert through an isomerization (i.e., folding) reaction. The end result is an analytical multiscale modeling approach which, although still schematic, can predict that folded proteins exhibit a closed-loop region of stability in the pressure-temperature plane and that protein concentration has a nonmonotonic effect on protein stability, results consistent with qualitative trends observed in both experiments of protein solutions and simulations of coarse-grained protein models.

Effect of field direction on electrowetting in a nanopore

We manifest a significant influence of field direction and polarity on surface wetting, when the latter is tuned by application of an external electric field. Thermodynamics of field-induced filling of hydrocarbon-like nanopores with water is studied by open ensemble molecular simulation. Increased field strength consistently results in water-filling and electrostriction in hydrophobic nanopores. A threshold field commensurate with surface charge density of about one elementary charge per 10 nm2 suffices to render prototypical paraffin surfaces hydrophilic. When a field is applied in the direction perpendicular to the confining walls, the competition between orientational polarization and angle preferences of interfacial water molecules relative to the walls results in an asymmetric wettability of opposing surfaces (Janus interface). Reduction of surface free energy observed upon alignment of confinement walls with field direction suggests a novel mechanism whereby the applied electric field can operate selectively on water-filled nanotubes while empty ones remain unaffected.

Effect of water-wall interaction potential on the properties of nanoconfined water

Much of the understanding of bulk liquids has progressed through study of the limiting case in which molecules interact via purely repulsive forces, such as a hard-core or "repulsive ramp" potential. In the same spirit, we report progress on the understanding of confined water by examining the behavior of waterlike molecules interacting with planar walls via purely repulsive forces and compare our results with those obtained for Lennard-Jones (LJ) interactions between the molecules and the walls. Specifically, we perform molecular dynamics simulations of 512 waterlike molecules interacting via the TIP5P potential and confined between two smooth planar walls that are separated by 1.1nm . At this separation, there are either two or three molecular layers of water, depending on density. We study two different forms of repulsive confinement, when the water-wall interaction potential is either (i) 1r;{9} or (ii) a WCA-like repulsive potential. We find that the thermodynamic, dynamic, and structural properties of the liquid in purely repulsive confinements qualitatively match those for a system with a pure LJ attraction to the wall. In previous studies that include attractions, freezing into monolayer or trilayer ice was seen for this wall separation. Using the same separation as these previous studies, we find that the crystal state is not stable with 1r;{9} repulsive walls but is stable with WCA-like repulsive confinement. However, by carefully adjusting the separation of the plates with 1r;{9} repulsive interactions so that the effective space available to the molecules is the same as that for LJ confinement, we find that the same crystal phases are stable. This result emphasizes the importance of comparing systems only using the same effective confinement, which may differ from the geometric separation of the confining surfaces.

Fluid in a closed narrow slit

The behavior of a fluid inside a closed narrow slit between solid walls is examined on the basis of the density functional theory. It is shown that the constraint of constant number of molecules leads to interesting effects which are absent when the slit is open and in contact with a reservoir. If the slit walls are identical, the density profiles at low temperatures or at high average densities rhoav of the fluid molecules in the slit have a sharp maximum in the middle of the slit, the value of the density at maximum being comparable to that of a liquid. The density of fluid at the walls is in this case comparable to the density of a vapor phase. At high temperatures or at low rhoav the fluid density in the middle of the slit is of the same order of magnitude as at the walls. For nonidentical walls the density maximum is shifted towards the wall with a stronger wall-fluid interaction. The transition between the two types (with and without the sharp maximum) of density profiles with the change of temperature in the slit occurs in a narrow range of temperatures, this range being larger for narrower slits. The pressures which the fluid exerts on the walls as well as the forces per unit area arising due to stresses in the sidewalls of the system can decrease with increasing rhoav. Such a behavior is not possible for homogeneous systems and can be explained by analyzing the fluid density at the walls when rhoav increases. The normal and transversal components of the pressure tensor were calculated as functions of the distance from the wall using the equation of hydrostatic equilibrium and direct calculation of the forces between molecules, respectively. The normal component decreases with increasing distance near the wall in contrast to the normal component near the liquid-vapor interface reported previously in the literature. The behavior of the transverse component does not depend on the fluid-solid interaction and is comparable to that for a liquid-vapor interface.

Heteropolymer collapse theory for protein folding in the pressure-temperature plane

We revisit a heteropolymer collapse theory originally introduced to explore how the balance between hydrophobic interactions and configurational entropy determines the thermal stability of globular proteins at ambient pressure. We generalize the theory by introducing a basic statistical mechanical treatment for how pressure impacts the solvent-mediated interactions between hydrophobic amino-acid residues. In particular, we estimate the strength of the hydrophobic interactions using a molecular thermodynamic model for the interfacial free energy between liquid water and a curved hydrophobic solute. The model, which also reproduces many of the distinctive thermodynamic properties of aqueous solutions in bulk and interfacial environments, predicts that the water-solute interfacial free energy is significantly reduced by the application of high hydrostatic pressures. This allows water to penetrate into folded heteropolymers at high pressure and break apart their hydrophobic cores, a scenario suggested earlier by information theory calculations. As a result, folded heteropolymers are predicted to display the kind of closed region of stability in the pressure-temperature plane exhibited by native proteins. We compare predictions of the collapse theory with experimental data for several proteins.

Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates

We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T = 300K. We investigate the effect of pressure (-0.15 GPa< or = P< or =0.2GPa) and plate separation (0.4 nm < or =d < or =1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between the plates at low enough P. The threshold value of d at which this transition occurs decreases with P (e.g., 1.6 nm at P approximately equal to -0.05 GPa, 0.5 nm at P approximately equal to 0.1 GPa), until, at high P, no capillary evaporation occurs. For d approximately equal to 0.6 nm and P > or =0.1 GPa, the system crystallizes into a bilayer ice. A P-d phase diagram showing the vapor, liquid, and bilayer ice phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains in the liquid phase at all P and d studied. Interestingly, we observe for this case that even at the P at which bulk water cavitates, the confined water remains in the liquid state. We also study systematically the state of hydration at different P for both kinds of plates. For the range of conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to enhance water structure and to push water molecules toward the plates. The average orientation of water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that upon pressurization, hydrophobic plates behave as "soft" surfaces (in the sense of accommodating pressure-dependent changes in water structure) while hydrophilic walls behave as "hard" surfaces.

Molecular simulation of liquid water confined inside graphite channels: Thermodynamics and structural properties

We carried out molecular dynamics simulations to describe the properties of water inside a narrow graphite channel. Two stable phases were found: a low-density one made of water clusters adsorbed on the graphite sheets and a liquid one that fills the entire channel, forming several layers around a bulk-like region. We analyzed the interfacial structure, orientational order, water residence times in several regions, and hydrogen bonding of this last water phase, calculating also a quantity of electrochemical interest, the probability of electron tunneling through interfacial water. The results are in good qualitative agreement with the available experimental data.

Gas Solubility in Hydrophobic Confinement

Measured forces between apolar surfaces in water have often been found to be sensitive to exposure to atmospheric gases despite low gas solubilities in bulk water. This raises questions as to how significant gas adsorption is in hydrophobic confinement, whether it is conducive to water depletion at such surfaces, and ultimately if it can facilitate the liquid-to-gas phase transition in the confinement. Open Ensemble molecular simulations have been used here to determine saturated concentrations of atmospheric gases in water-filled apolar confinements as a function of pore width at varied gas fugacities. For paraffin-like confinements of widths barely exceeding the mechanical instability threshold (spinodal) of the liquid-to-vapor transition of confined water (aqueous film thickness between three and four molecular diameters), mean gas concentrations in the pore were found to exceed the bulk values by a factor of approximately 30 or approximately 15 in cases of N2 and CO2, respectively. At ambient conditions, this does not result in visible changes in the water density profile next to the surfaces. Whereas the barrier to capillary evaporation has been found to decrease in the presence of dissolved gas (Leung, K.; Luzar, A.; and Bratko, D. Phys. Rev. Lett. 2003, 90, 065502), gas concentrations much higher than those observed at normal atmospheric conditions would be needed to produce noticeable changes in the kinetics of capillary evaporation. In simulations, dissolved gas concentrations corresponding to fugacities above approximately 40 bar for N2, or approximately 2 bar for CO2, were required to trigger expulsion of water from a hydrocarbon slit as narrow as 1.4 nm. For nanosized pore widths corresponding to the mechanical instability threshold or above, no significant coupling between adsorption layers at opposing confinement walls was observed. This finding explains the approximately linear increase in gas solubility with inverse confinement width and the apparent validity of Henry's law in the pores over a broad fugacity range.

Thermodynamics, structure, and dynamics of water confined between hydrophobic plates

We perform molecular dynamics simulations of 512 waterlike molecules that interact via the TIP5P potential and are confined between two smooth hydrophobic plates that are separated by 1.10 nm. We find that the anomalous thermodynamic properties of water are shifted to lower temperatures relative to the bulk by approximately 40 K. The dynamics and structure of the confined water resemble bulk water at higher temperatures, consistent with the shift of thermodynamic anomalies to lower temperature. Because of this T shift, our confined water simulations (down to T=220 K) do not reach sufficiently low temperature to observe a liquid-liquid phase transition found for bulk water at T approximately 215 K using the TIP5P potential, but we see inflections in isotherms at lower temperatures presumably due to the presence of a liquid-liquid critical point. We find that the different crystalline structures that can form for two different separations of the plates, 0.7 and 1.10 nm, have no counterparts in the bulk system, and we discuss the relevance to experiments on confined water.

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.

Hydrophobic hydration from small to large lengthscales: Understanding and manipulating the crossover

Small and large hydrophobic solutes exhibit remarkably different hydration thermodynamics. Small solutes are accommodated in water with minor perturbations to water structure, and their hydration is captured accurately by theories that describe density fluctuations in pure water. In contrast, hydration of large solutes is accompanied by dewetting of their surfaces and requires a macroscopic thermodynamic description. A unified theoretical description of these lengthscale dependencies was presented by Lum, Chandler, and Weeks [(1999) J. Phys. Chem. B 103, 4570-4577]. Here, we use molecular simulations to study lengthscale-dependent hydrophobic hydration under various thermodynamic conditions. We show that the hydration of small and large solutes displays disparate dependencies on thermodynamic variables, including pressure, temperature, and additive concentration. Understanding these dependencies allows manipulation of the small-to-large crossover lengthscale, which is nanoscopic under ambient conditions. Specifically, applying hydrostatic tension or adding ethanol decreases the crossover length to molecular sizes, making it accessible to atomistic simulations. With detailed temperature-dependent studies, we further demonstrate that hydration thermodynamics changes gradually from entropic to enthalpic near the crossover. The nanoscopic lengthscale of the crossover and its sensitivity to thermodynamic variables imply that quantitative modeling of biomolecular self-assembly in aqueous solutions requires elements of both molecular and macroscopic hydration physics. We also show that the small-to-large crossover is directly related to the Egelstaff-Widom lengthscale, the product of surface tension and isothermal compressibility, which is another fundamental lengthscale in liquids.

Modeling Water, the Hydrophobic Effect, and Ion Solvation

Water plays a central role in the structures and properties of biomolecules--proteins, nucleic acids, and membranes--and in their interactions with ligands and drugs. Over the past half century, our understanding of water has been advanced significantly owing to theoretical and computational modeling. However, like the blind men and the elephant, different models describe different aspects of water's behavior. The trend in water modeling has been toward finer-scale properties and increasing structural detail, at increasing computational expense. Recently, our labs and others have moved in the opposite direction, toward simpler physical models, focusing on more global properties-water's thermodynamics, phase diagram, and solvation properties, for example-and toward less computational expense. Simplified models can guide a better understanding of water in ways that complement what we learn from more complex models. One ultimate goal is more tractable models for computer simulations of biomolecules. This review gives a perspective from simple models on how the physical properties of water-as a pure liquid and as a solvent-derive from the geometric and hydrogen bonding properties of water.

On the Mechanism of Hydrophobic Association of Nanoscopic Solutes

The hydration behavior of two planar nanoscopic hydrophobic solutes in liquid water at normal temperature and pressure is investigated by calculating the potential of mean force between them at constant pressure as a function of the solute-solvent interaction potential. The importance of the effect of weak attractive interactions between the solute atoms and the solvent on the hydration behavior is clearly demonstrated. We focus on the underlying mechanism behind the contrasting results obtained in various recent experimental and computational studies on water near hydrophobic solutes. The length scale where crossover from a solvent separated state to the contact pair state occurs is shown to depend on the solute sizes as well as on details of the solute-solvent interaction. We find the mechanism for attractive mean forces between the plates is very different depending on the nature of the solute-solvent interaction which has implications for the mechanism of the hydrophobic effect for biomolecules.

Drying and Hydrophobic Collapse of Paraffin Plates

We perform molecular dynamics simulations of the hydrophobic collapse of two paraffin plates to examine how the collapse is mediated by realistic paraffin-water attractive van der Waals forces. We explore several aspects of the drying transition between the plates, including the critical separation for drying and the critical size of the vapor bubble required for the nucleation of the drying event. We also investigate the kinetics of hydrophobic collapse and find that the hydrophobic collapse occurs in about 100 ps. We compare these results with the simulations with the plate-water van der Waals attractions turned off and with recent results on the hydrophobic collapse of multidomain proteins. Last, we discuss the relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on solute surface using a simple macroscopic theory.

Water between Plates in the Presence of an Electric Field in an Open System

Molecular dynamics simulations of water at 298 K and 1 atm of pressure are used to investigate the electric-field dependence of the density and polarization density of water between two graphite-like plates of different sizes (9.8 x 9.2 and 17.7 x 17.2 A) in an open system for plate separations of 8.0, 9.5, and 16.4 A. The interactions with water were tuned to "hard-wall-like" and "normal" C-O hydrophobic potentials. Water between the larger plates at 16.4 A separation is layered but is metastable with respect to capillary evaporation at zero field (Bratko, D.; Curtis, R. A.; Blanch, H. W.; Prausnitz, J. M. J. Chem. Phys. 2001, 115, 3873). Applying a field decreases the density of the water between the plates, in apparent contradiction to thermodynamic and integral equation theories of bulk fluid electrostriction that ignore surface effects, rendering them inapplicable to finite-sized films of water between hydrophobic plates. This suggests that the free energy barrier for evaporation is lowered by the applied field. Water, between "hard-wall-like" plates at narrower separations of 9.5 A and less, shows a spontaneous but incomplete evaporation at zero field within the time scale of our simulation. Evaporation is further enhanced by an electric field. No such evaporation occurs, on these time scales, for the smaller plates with the "hard-wall-like" potential at a separation of 8.0 A at zero field, signaling a crossover in behavior as the plate dimension decreases, but the water density still diminishes with increasing field strength. These observations could have implications for the behavior of thin films of water between surfaces in real physical and biological systems.

Nanoscale hydrophobic interaction and nanobubble nucleation

We report large-scale atomistic simulation of midrange nanoscale hydrophobic interaction, manifested by the nucleation of nanobubble between nanometer-sized hydrophobes at constrained equilibrium. When the length scale of the hydrophobes is greater than 2 nm, the nanobubble formation shows hysteresis behavior resembling the first-order transition. Calculation of the potential of mean force versus interhydrophobe distance provides a quantitative measure of the strength of the nanoscale hydrophobic interaction.

A computer simulation study of water drying at the interface of protein chains

This study investigated the water drying (cavitation) in the interfacial region of two chains of a dimeric protein by nanosecond molecular dynamics simulations using explicit water representation. Separation-induced cavity of water was directly observed in the region. We evaluated the separation length scale of two chains on which the drying transition occurs, and the average number of water molecules that are expelled from the interfacial region during the transition. The obtained values can be rationalized by Kelvin equation for finite lateral size of confinement [K. Lum and A. Luzar, Phys. Rev. E 56, R6283 (1997)]. Also, we found that the drying transition is accompanied by an exponential reduction in the average hydrogen-bond number per interfacial water molecule. The results of this study may deepen the understanding of how hydrophobic interaction drives the assembly of protein chains.

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: A detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.

Dynamical behavior of one-dimensional water molecule chains in zeolites: Nanosecond time-scale molecular dynamics simulations of bikitaite

Nanosecond scale molecular dynamics simulations of the behavior of the one-dimensional water molecule chains adsorbed in the parallel nanochannels of bikitaite, a rare lithium containing zeolite, were performed at different temperatures and for the fully and partially hydrated material. New empirical potential functions have been developed for representing lithium-water interactions. The structure and the vibrational spectrum of bikitaite were in agreement both with experimental data and Car-Parrinello molecular dynamics results. Classical molecular dynamics simulations were extended to the nanosecond time scale in order to study the flip motion of water molecules around the hydrogen bonds connecting adjacent molecules in the chains, which has been observed by NMR experiments, and the dehydration mechanism at high temperature. Computed relaxation times of the flip motion follow the Arrhenius behavior found experimentally, but the activation energy of the simulated system is slightly underestimated. Based on the results of the simulations, it may be suggested that the dehydration proceeds by a defect-driven stepwise diffusion. The diffusive mechanism appears as a single-file motion: the molecules never pass one another, even at temperatures as high as about 1000 K, nor can they switch between different channels. However, the mean square displacement (MSD) of the molecules, computed with respect to the center of mass of the simulated system, shows an irregular trend from which the single-file diffusion cannot be clearly evidenced. If the MSDs are evaluated with respect to the center of mass of the molecules hosted in each channel, the expected dependence on the square root of time finally appears.

A simple analytical model of water

Water is an unusual liquid. It expands upon freezing, has minima in its volume, heat capacity, and isothermal compressibility with temperature, and shows signs of a first-order phase transition when supercooled. These anomalies disappear at high pressures. We review a recent analytical theory that predicts water's thermal properties and the main features of its phase diagram, including multiple crystalline phases and a fluid-fluid transition in the supercooled liquid. It also predicts a fragile-to-strong crossover in supercooled water's temperature-dependent relaxation processes. The theory is based on a simplified model for how triplets of waters interact via hydrogen bonds, steric repulsions, and dispersion attractions. It is designed to give simple insights into the microscopic origins of water's properties.

Dynamics of capillary drying in water

We use atomistic simulations to address the question when capillary evaporation of water confined in a hydrocarbonlike slit is kinetically viable. Activation barriers and absolute rates of evaporation are estimated using open ensemble Monte Carlo-umbrella sampling and molecular dynamics simulations. At ambient conditions, the evaporation rate in a water film four molecular diameters thick is found to be of the order 10(5)(nm(2) s)(-1), meaning that water readily evaporates. Films more than a few nanometers thick will persist in a metastable liquid state. Dissolved atmospheric gas molecules do not significantly decrease the activation barrier.

Prediction of beta-sheet topology and disulfide bridges in polypeptides

An ab initio method has been developed to predict beta architectures in polypeptides. The approach predicts the topology of beta-sheets and disulfide bridges through a novel superstructure-based mathematical framework originally established for chemical process synthesis problems. Two types of superstructure are introduced, both of which emanate from the principle that hydrophobic interactions drive the formation of a beta-structure. The mathematical formulation of the problem results in a set of integer linear programming (ILP) problems that can be solved to global optimality to identify the optimal beta-configuration. These (ILP) models can also predict a ranked ordered list of the best, second-best, third-best, etc., topologies of beta-sheets and disulfide bridges. The approach is shown to perform very well for several benchmark polypeptide systems, as well as polypeptides exhibiting challenging nonsequential beta-sheet topologies folds (56 to 187 amino acids).

Connecting local structure to interface formation: a molecular scale van der Waals theory of nonuniform liquids

This article reviews a new and general theory of nonuniform fluids that naturally incorporates molecular scale information into the classical van der Waals theory of slowly varying interfaces. The method optimally combines two standard approximations, molecular (mean) field theory to describe interface formation and linear response (or Gaussian fluctuation) theory to describe local structure. Accurate results have been found in many different applications in nonuniform simple fluids and these ideas may have important implications for the theory of hydrophobic interactions in water.