Reorientational Dynamics of Water Confined in Zeolites

The use of collective variables and enhanced sampling in the simulations of existing and emerging microporous materials

Microporous materials, including zeolites, metal-organic frameworks, and cage compounds, offer diverse functionalities due to their unique dynamics and guest confinement properties. These materials play a significant role in separation, catalysis, and sensing, but their complexity hinders exploration using traditional atomistic simulations. This review explores collective variables (CVs) paired with enhanced sampling as a powerful approach to enable efficient investigation of key features in microporous materials. We highlight successful applications of CVs in studying adsorption, diffusion, phase transitions, and mechanical properties, demonstrating their crucial role in guiding material design and optimisation. The future of CVs lies in integration with techniques like machine learning, allowing for enhanced efficiency and accuracy. By tailoring CVs to specific materials and developing multi-scale approaches we can further unlock the intricacies of these fascinating materials. Simulations are a cornerstone in unravelling the complexities of microporous materials and are crucial for our future understanding.

The influence of force fields on the structure and dynamics of water confined in ZIF-8 from atomistic simulations

The complexity of modeling flexible crystals, such as ZIF-8, mainly stems from the handling of intramolecular interactions. Numerous force fields have been proposed in the literature to describe the interactions between atoms in ZIF-8. We employ seven force fields to examine the structure and dynamic behavior of water molecules confined in ZIF-8, with the aim of investigating the impact of force fields on simulation results. Various structural characterization methods consistently indicate that the choice of different force fields has quantitative effects but no qualitative effects on the structural characteristics of confined water. Additionally, the force fields do not impact the qualitative description of the diffusion mechanism. Both mean-square displacement and van Hove autocorrelation function reveal two characteristic movements of water molecules diffusing in ZIF-8: a short-time intra-cavity hopping process and a long-time inter-cavity hopping process. However, the framework flexibility is found to play a crucial role in determining the order of spatial arrangement and local structure, self-diffusion coefficient and reorientational dynamics of confined water. Specifically, the DREIDING force field gives rise to an unrealistic stiff framework, enhancing the order of spatial arrangement and diminishing the local ordered structure of confined water. Meanwhile, it results in much slower translational and reorientational dynamics. Hence, the general DREIDING force field cannot be considered for providing a quantitative description of the water structure and dynamics.

Competing Nuclear Quantum Effects and Hydrogen-Bond Jumps in Hydrated Kaolinite

Recent work has shown that the dynamics of hydrogen bonds in pure clays are affected by nuclear quantum fluctuations, with different effects for the hydrogen bonds holding different layers of the clay together and for those within the same layer. At the clay-water interface there is an even wider range of types of hydrogen bond, suggesting that the quantum effects may be yet more varied. We apply classical and thermostated ring polymer molecular dynamics simulations to show that nuclear quantum effects accelerate hydrogen-bond dynamics to varying degrees. By interpreting the results in terms of the extended jump model of hydrogen-bond switching, we can understand the origins of these effects in terms of changes in the quantum kinetic energy of hydrogen atoms during an exchange. We also show that the extended jump mechanism is applicable not only to the hydrogen bonds involving water, but also those internal to the clay.

Quantifying the Effect of Spatial Confinement on the Diffusion and Nuclear Magnetic Resonance Relaxation of Water/Hydrocarbon Mixtures: A Molecular Dynamics Simulation Study

We implement molecular dynamics (MD) simulations to explore the relaxation mechanisms involved in the description of translational diffusion and rotational dynamics in water/hydrocarbon interfaces. The analysis of density profiles, self-diffusion coefficients, and nuclear magnetic resonance (NMR) relaxation properties as a function of the confinement layer width and type of hydrocarbons improves the understanding of confined water properties at water/oil interfaces. Density profile fluctuations reveal the presence of water-oil interactions close to the interface. MD results show that self-diffusion coefficients and NMR relaxation times of planar and cylindrical water/oil interfaces are strongly influenced by layer thickness and geometry. Shorter (between 20 and 60%) self-diffusion coefficients and ¹H NMR relaxation times were obtained for water/n-pentane, water/n-decane, and water/n-hexadecane systems than bulk diffusion coefficients. An increase in Larmor frequency from 2.3 MHz to 400 MHz shows that longitudinal relaxation time (T₁) of confined oil has slightly larger differences at higher frequencies than the transverse relaxation time (T₂). At 400 MHz, n-alkanes (n-pentane, n-decane, and n-hexadecane) exhibit longer relaxation times than at smaller frequency values (2.3 and 22 MHz). Analysis of spin-spin and spin-lattice times provides relevant information about inter- and intramolecular relaxation mechanisms of water and oil as a function of geometry and width of the interface layer. These MD results suggest that the strength of confinement and geometry play a vital role in the diffusion and NMR relaxation properties of water/oil interfaces.

Structure and dynamics of a water/methanol mixture confined in zeolitic imidazolate framework ZIF-8 from atomistic simulations

A classical atomistic simulation study is reported for the microscopic structure and dynamics of a water/methanol mixture confined in flexible nanoporous zeolitic imidazolate framework ZIF-8. Both the radial density distribution and vivid two-dimensional density profile demonstrate that methanol molecules can roughly be viewed as "embedded" between two layers of water molecules to form a "sandwich" structure. The reason for the formation of such a specific structure is explained based on the hydrogen-bonding state and the strength of various hydrogen bonds. The investigation of guest molecular diffusion shows that the self-diffusion coefficient of confined water is generally one to two orders of magnitude smaller than that of bulk water. In addition, the dependence of the self-diffusion coefficient on loading is non-monotonic: the self-diffusion coefficient firstly shows a significant increase and then decreases at higher loading. Moreover, both the structure and dynamics of the hydrogen bond (HB) network of confined water molecules are investigated in a spatially resolved manner. The results indicate that both the HB structure and dynamics of water molecules near the ZIF-8 surface deviate significantly from those of bulk water. However, while water molecules located at the pore center are relatively similar to bulk water molecules with respect to the HB structure, they exhibit strong slowdown in HB dynamics when compared with bulk water. This simulation study elucidates in detail the structural and dynamical properties of a water/methanol mixture in nanoscopic ZIF-8 confinement, which is expected to provide a deep insight into the role of porous fillers, such as ZIF-8, in improving the performance of the dehydration of alcohols via pervaporation and other related processes.

Adsorption of water in Na-LTA zeolites: an ab initio molecular dynamics investigation

The very wide range of applications of LTA zeolites, including the storage of tritiated water, implies that a detailed and accurate atomic-scale description of the adsorption processes taking place in their structure is crucial. To unravel with an unprecedented accuracy the mechanisms behind the water filling in NaA, we have conducted a systematic ab initio molecular dynamics investigation. Two LTA structural models, the conventional Z4A and the reduced one ZK4, have been used for static and dynamic ab initio calculations, respectively. After assessing this reduced model with comparative static DFT calculations, we start the filling of the α and β cages by water, molecule by molecule. This allowed us to thoroughly study the interaction of water molecules with the zeolite structure and between water molecules, progressively forming H-bond chains and ring patterns as the cage is being filled. The adsorption energies could then be calculated with an unprecedented accuracy, which showed that the interaction of the molecules with the zeolite weakens as their number increases. By these methods, we have been able to highlight the primary role of Na⁺ cations in the interaction of water with zeolite, and inversely, the role of water in the displacement of cations when it is sufficiently solvated, allowing the passage between the α and β cages. This phenomenon is possible thanks to the inhomogeneous distribution of water molecules on the cationic sites, as shown by our AIMD simulations, which allows the formation of water clusters. These results are important because they help in understanding how the coverage of cationic sites by water will affect the adsorption of other molecules inside the Na-LTA zeolite.

Diffusion of confined fluids in microporous zeolites and clay materials

Fluids exhibit remarkable variation in their structural and dynamic properties when they are confined at the nanoscopic scale. Various factors, including geometric restriction, the size and shape of the guest molecules, the topology of the host, and guest-host interactions, are responsible for the alterations in these properties. Due to their porous structures, aluminosilicates provide a suitable host system for studying the diffusion of sorbates in confinement. Zeolites and clays are two classes of the aluminosilicate family, comprising very ordered porous or layered structures. Zeolitic materials are important due to their high catalytic activity and molecular sieving properties. Guest molecules adsorbed by zeolites display many interesting features including unidimensional diffusion, non-isotropic rotation, preferred orientation and levitation effects, depending on the guest and host characteristics. These are useful for the separation of hydrocarbons which commonly exist as mixtures in nature. Similarly, clay materials have found application in catalysis, desalination, enhanced oil recovery, and isolation barriers used in radioactive waste disposal. It has been shown that the bonding interactions, level of hydration, interlayer spacing, and number of charge-balancing cations are the important factors that determine the nature of diffusion of water molecules in clays. Here, we present a review of the current status of the diffusion mechanisms of various adsorbed species in different microporous zeolites and clays, as investigated using quasielastic neutron scattering and classical molecular dynamics simulation techniques. It is impossible to write an exhaustive review of the subject matter, as it has been explored over several decades and involves many research topics. However, an effort is made to cover the relevant issues specific to the dynamics of different molecules in microporous zeolites and clay materials and to highlight a variety of interesting features that are important for both practical applications and fundamental aspects.

Structure and chemistry of graphene oxide in liquid water from first principles

Graphene oxide is a rising star among 2D materials, yet its interaction with liquid water remains a fundamentally open question: experimental characterization at the atomic scale is difficult, and modeling by classical approaches cannot properly describe chemical reactivity. Here, we bridge the gap between simple computational models and complex experimental systems, by realistic first-principles molecular simulations of graphene oxide (GO) in liquid water. We construct chemically accurate GO models and study their behavior in water, showing that oxygen-bearing functional groups (hydroxyl and epoxides) are preferentially clustered on the graphene oxide layer. We demonstrated the specific properties of GO in water, an unusual combination of both hydrophilicity and fast water dynamics. Finally, we evidence that GO is chemically active in water, acquiring an average negative charge of the order of 10 mC m-2. The ab initio modeling highlights the uniqueness of GO structures for applications as innovative membranes for desalination and water purification.

Shared hydrogen bonds: water in aluminated faujasite

Water confined in faujasite, a zeolite, with aluminium content, exhibits properties different from those of bulk water as well as water confined in siliceous faujasite. The RDF between oxygen of water (OW) and oxygen of aluminium (OAl) shows a prominent first peak near to 2.9 Å similar to any oxygen-oxygen RDF seen in bulk water and unlike water confined in siliceous faujasite. Further, HW-OAl shows a peak near 1.9 Å suggesting hydrogen bonding between hydrogen of water and OAl. The water satisfies the hydrogen bond criteria with both O1Al and O2Al indicating that it is participating in a shared hydrogen bond. The hydrogen bond exchange between such a water forming a shared hydrogen bond to OAl and another water molecule H2Ob is investigated through the changes in the distances and appropriate angles. The O-Al-O angle of the zeolite increases by about 7 degrees on the formation of the shared hydrogen bond. The jump dynamics of the shared hydrogen bond when the two bonds break simultaneously has been obtained and this is reported. This jump reorientation dynamics is different compared to normal hydrogen bonding reported by Laage and Hynes: it has a short lifetime, around 50-100 fs computed from SHB(t). The intermittent and continuous hydrogen bond correlation functions are also reported.

Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes

UiO‐66, composed by Zr‐oxide inorganic bricks [Zr₆(μ₃‐O)₄(μ₃‐OH)₄] and organic terephthalate linkers, is one of the most studied metal–organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties. In this work, it is shown, by means of first‐principle molecular‐dynamics simulations at operating conditions, how protic solvents may facilitate such changes in the metal coordination. Solvent can induce structural rearrangements in the material that can lead to undercoordinated but also overcoordinated metal sites. This is demonstrated by simulating activation processes along well‐chosen collective variables. Such enhanced MD simulations are able to track the intrinsic dynamics of the framework at realistic conditions.

Interfacial properties modulated by the water confinement in reverse micelles created by the ionic liquid-like surfactant bmim-AOT

The behavior of the interfacial water entrapped in reverse micelles (RMs) that were formed by the ionic liquid-like surfactant 1-butyl-3-methylimidazolium 1,4-bis-2-ethylhexylsulfosuccinate (bmim-AOT) was investigated with the use of UV-Vis absorption spectroscopy and nuclear magnetic resonance (NMR) relaxometry. The solvatochromism of two molecular probes, namely, 1-methyl-8-oxyquinolinium betaine (QB) and N,N,N',N'-tetramethylethylenediamine copper(ii)acetylacetonate tetraphenylborate ([Cu(acac)(tmen)][B(C6H5)4]), was investigated. As a comparison, the analog RMs formed by sodium 1,4-bis-2-ethylhexylsulfosuccinate (Na-AOT) were also explored. By varying the water content inside the RMs and consequently the different magnitude of the water-surfactant interactions at the interface, interesting properties were observed by comparing bmim-AOT and Na-AOT RMs. From the solvatochromic behavior of ([Cu(acac)(tmen)][B(C6H5)4]), we found that the interface in bmim-AOT RMs shows a smaller electron donating capacity than that in Na-AOT RMs. QB revealed that the interfacial region is a weaker hydrogen bond donor and less polar than the corresponding Na-AOT RMs. NMR experiments showed that the molecular motion of water in bmim-AOT RMs is less restricted than that of the water molecules confined in Na-AOT RMs. In summary, the results show how the nature of the bmim+ cation affects the interaction between the entrapped water and the RM interface, greatly modifying the interfacial water structure in comparison with the results known for Na-AOT.

Structure and Dynamics of Water Confined in Imogolite Nanotubes

We have studied the properties of water adsorbed inside nanotubes of hydrophilic imogolite, an aluminum silicate clay mineral, by means of molecular simulations. We used a classical force field to describe the water and the flexible imogolite nanotube and validated it against the data obtained from first-principles molecular dynamics. With it, we observe a strong structuration of the water confined in the nanotube, with specific adsorption sites and a distribution of hydrogen bond patterns. The combination of number of adsorption sites, their geometry, and the preferential tetrahedral hydrogen bonding pattern of water leads to frustration and disorder. We further characterize the dynamics of the water, as well as the hydrogen bonds formed between water molecules and the nanotube, which is found to be more than 1 order of magnitude longer than water-water hydrogen bonds.

Ab Initio Simulations of Water Dynamics in Aqueous TMAO Solutions: Temperature and Concentration Effects

We use ab initio molecular dynamics simulation to study the effect of hydrophobic groups on the dynamics of water molecules in aqueous solutions of trimethylamine N-oxide (TMAO). We show that hydrophobic groups induce a moderate (<2-fold) slowdown of water reorientation and hydrogen-bond dynamics in dilute solutions, but that this slowdown rapidly increases with solute concentration. In addition, the slowdown factor is found to vary very little with temperature, thus suggesting an entropic origin. All of these results are in quantitative agreement with prior classical molecular dynamics simulations and with the previously suggested excluded-volume model. The hydrophilic TMAO headgroup is found to affect water dynamics more strongly than the hydrophobic moiety, and the magnitude of this slowdown is very sensitive to the strength of the water-solute hydrogen-bond.

Photochemistry and Photophysics in Silica-Based Materials: Ultrafast and Single Molecule Spectroscopy Observation

Silica-based materials (SBMs) are widely used in catalysis, photonics, and drug delivery. Their pores and cavities act as hosts of diverse guests ranging from classical dyes to drugs and quantum dots, allowing changes in the photochemical behavior of the confined guests. The heterogeneity of the guest populations as well as the confinement provided by these hosts affect the behavior of the formed hybrid materials. As a consequence, the observed reaction dynamics becomes significantly different and complex. Studying their photobehavior requires advanced laser-based spectroscopy and microscopy techniques as well as computational methods. Thanks to the development of ultrafast (spectroscopy and imaging) tools, we are witnessing an increasing interest of the scientific community to explore the intimate photobehavior of these composites. Here, we review the recent theoretical and ultrafast experimental studies of their photodynamics and discuss the results in comparison to those in homogeneous media. The discussion of the confined dynamics includes solvation and intra- and intermolecular proton-, electron-, and energy transfer events of the guest within the SBMs. Several examples of applications in photocatalysis, (photo)sensors, photonics, photovoltaics, and drug delivery demonstrate the vast potential of the SBMs in modern science and technology.

Forced intrusion of water and aqueous solutions in microporous materials: from fundamental thermodynamics to energy storage devices

We review the high pressure forced intrusion studies of water in hydrophobic microporous materials such as zeolites and MOFs, a field of research that has emerged some 15 years ago and is now very active. Many of these studies are aimed at investigating the possibility of using these systems as energy storage devices. A series of all-silica zeolites (zeosil) frameworks were found suitable for reversible energy storage because of their stability with respect to hydrolysis after several water intrusion-extrusion cycles. Several microporous hydrophobic zeolite imidazolate frameworks (ZIFs) also happen to be quite stable and resistant towards hydrolysis and thus seem very promising for energy storage applications. Replacing pure water by electrolyte aqueous solutions enables to increase the stored energy by a factor close to 3, on account of the high pressure shift of the intrusion transition. In addition to the fact that aqueous solutions and microporous silica materials are environmental friendly, these systems are thus becoming increasingly interesting for the design of new energy storage devices. This review also addresses the theoretical approaches and molecular simulations performed in order to better understand the experimental behavior of nano-confined water. Molecular simulation studies showed that water condensation takes place through a genuine first-order phase transition, provided that the interconnected pores structure is 3-dimensional and sufficiently open. In an extreme confinement situations such as in ferrierite zeosil, condensation seem to take place through a continuous supercritical crossing from a diluted to a dense fluid, on account of the fact that the first-order transition line is shifted to higher pressure, and the confined water critical point is correlatively shifted to lower temperature. These molecular simulation studies suggest that the most important features of the intrusion/extrusion process can be understood in terms of equilibrium thermodynamics considerations.

Subtleties of catanionic surfactant reverse micelle assemblies revealed by a fluorescent molecular probe

In this work, the absorption and emission behavior of the cationic hemicyanine trans-4-[4-(dimethylamino)styryl]-N-methylpyridinium iodide (HC) in reverse micelles (RMs) formed by the catanionic surfactants benzyl-n-hexadecyldimethylammonium-1,4-bis-2-ethylhexylsulfosuccinate (AOT-BHD) and cetyltrimethylammonium-1,4-bis-2-ethylhexylsulfosuccinate (AOT-CTA) have been investigated. Our results show that the spectroscopic behavior of HC changes when the dye is dissolved in AOT-BHD or in AOT-CTA RMs. While HC undergoes an intramolecular charge-transfer process upon excitation in AOT-CTA RMs, in AOT-BHD RMs this process is inhibited due to a specific interaction between HC and the polar head group of the BHD⁺ cation. This implies that the chemical structure of CTA⁺ and BHD⁺ cations has a large impact on the excited stated from which HC emission occurs. Additionally, the structural difference between the two cations impacts on the water-RM interface interaction, which provides a way of controlling the solvation process in these RMs. Furthermore, differences in the interfacial fluidity between the two catanionic RMs is observed, a result that is particularly interesting with regard to these systems being used as nanoreactors.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

A review of the structure and dynamics of nanoconfined water and ionic liquids via molecular dynamics simulation

During the past decade, the research on fluids in nanoconfined geometries has received considerable attention as a consequence of their wide applications in different fields. Several nanoconfined systems such as water and ionic liquids, together with an equally impressive array of nanoconfining media such as carbon nanotube, graphene and graphene oxide have received increasingly growing interest in the past years. Water is the first system that has been reviewed in this article, due to its important role in transport phenomena in environmental sciences. Water is often considered as a highly nanoconfined system, due to its reduction to a few layers of water molecules between the extended surface of large macromolecules. The second system discussed here is ionic liquids, which have been widely studied in the modern green chemistry movement. Considering the great importance of ionic liquids in industry, and also their oil/water counterpart, nanoconfined ionic liquid system has become an important area of research with many fascinating applications. Furthermore, the method of molecular dynamics simulation is one of the major tools in the theoretical study of water and ionic liquids in nanoconfinement, which increasingly has been joined with experimental procedures. In this way, the choice of water and ionic liquids in nanoconfinement is justified by applying molecular dynamics simulation approaches in this review article.

How the cation 1-butyl-3-methylimidazolium impacts the interaction between the entrapped water and the reverse micelle interface created with an ionic liquid-like surfactant

The behavior of the interfacial water entrapped in reverse micelles (RMs) formed by the ionic liquid-like surfactant 1-butyl-3-methylimidazolium 1,4-bis-2-ethylhexylsulfosuccinate (bmim-AOT) dissolved in benzene (or chlorobenzene) was investigated using noninvasive techniques such as dynamic light scattering (DLS), static light scattering (SLS), FT-IR and (1)H NMR. The DLS and SLS results reveal the formation of discrete spherical and non-interacting water droplets stabilized by the bmim-AOT surfactant. Moreover, since the droplet size increases as the W0 (W0 = [water]/[surfactant]) value increases, water interacts with the RM interface. From FT-IR and (1)H NMR data, a weaker water-surfactant interaction in bmim-AOT RMs in comparison with the RMs created by sodium 1,4-bis-2-ethylhexylsulfosuccinate (Na-AOT) is detected. Consequently, there are less water molecules interacting with the interface in bmim-AOT RMs, and their hydrogen bond network is not completely disrupted as they are in Na-AOT RMs. The results show how the nature of the new cation impacts the interaction between the entrapped water and the RM interface, modifying the interfacial water structure in comparison with the results known for Na-AOT.

Water Dynamics in Metal-Organic Frameworks: Effects of Heterogeneous Confinement Predicted by Computational Spectroscopy

The behavior of water confined in MIL-53(Cr), a flexible metal-organic framework (MOF), is investigated through computational infrared spectroscopy. As the number of molecules adsorbed inside of the pores increases, the water OH stretch band of the linear infrared spectrum grows in intensity and approaches that of bulk water. To assess whether the water confined in MIL-53(Cr) becomes liquid-like, two-dimensional infrared spectra (2DIR) are also calculated. Confinement effects result in distinct chemical environments that appear as specific features in the 2DIR spectra. The evolution of the 2DIR line shape as a function of waiting time is well described in terms of the orientational dynamics of the water molecules, with chemical exchange cross peaks appearing at a time scale similar to the hydrogen bond rearrangement lifetime. The confining environment considerably slows the hydrogen bond dynamics relative to bulk as a result of the competition between water-framework and water-water interactions.

Effect of the cationic surfactant moiety on the structure of water entrapped in two catanionic reverse micelles created from ionic liquid-like surfactants

The behavior of water entrapped in reverse micelles (RMs) formed by two catanionic ionic liquid-like surfactants, benzyl-n-hexadecyldimethylammonium 1,4-bis-2-ethylhexylsulfosuccinate (AOT-BHD) and cetyltrimethylammonium 1,4-bis-2-ethylhexylsulfosuccinate (AOT-CTA), was investigated by using dynamic (DLS) and static (SLS) light scattering, FTIR, and (1)H NMR spectroscopy techniques. To the best of our knowledge, this is the first report in which AOT-CTA has been used to create RMs and encapsulate water. DLS and SLS results revealed the formation of RMs in benzene and the interaction of water with the RM interface. From FTIR and (1)H NMR spectroscopy data, a difference in the magnitude of the water-catanionic surfactant interaction at the interface is observed. For the AOT-BHD RMs, a strong water-surfactant interaction can be invoked whereas for AOT-CTA this interaction seems to be weaker. Consequently, more water molecules interact with the interface in AOT-BHD RMs with a completely disrupted hydrogen-bond network, than in AOT-CTA RMs in which the water structure is partially preserved. We suggest that the benzyl group present in the BHD(+) moiety in AOT-BHD is responsible for the behavior of the catanionic interface in comparison with the interface created in AOT-CTA. These results show that a simple change in the cationic component in the catanionic surfactant promotes remarkable changes in the RMs interface with interesting consequences, in particular when using the systems as nanoreactors.

Temperature dependence of hydrophobic hydration dynamics: from retardation to acceleration

The perturbation induced by a hydrophobic solute on water dynamics is essential in many biochemical processes, but its mechanism and magnitude are still debated. A stringent test of the different proposed pictures is provided by recent NMR measurements by Qvist and Halle (J. Am. Chem. Soc. 2008, 130, 10345-10353) which showed that, unexpectedly, the perturbation changes in a non-monotonic fashion when the solution is cooled below room temperature. Here we perform and analyze molecular dynamics simulations of a small paradigm amphiphilic solute, trimethylamine N-oxide (TMAO), in dilute aqueous solutions over the 218-350 K temperature range. We first show that our simulations properly reproduce the non-monotonic temperature dependence. We then develop a model which combines our previously suggested entropic excluded-volume effect with a perturbation factor arising from the difference between local structural fluctuations in the shell and the bulk. Our model provides a detailed molecular understanding of the hydrophobic perturbation over the full temperature range investigated. It shows that the excluded-volume factor brings a dominant temperature-independent contribution to the perturbation at all temperatures, and provides a very good approximation at room temperature. The non-monotonic temperature dependence of the perturbation is shown to arise from the structural factor and mostly from relative shifts between the shell and bulk distributions of local structures, whose amplitude remains very small compared to the widths of those distributions.