Topological building blocks of hydrogen bond network in water

Structural effects of water clusters on viscosity at high shear rates

In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.

Mechanism of Cations Suppressing Proton Diffusion Kinetics for Electrocatalysis

Proton transfer is crucial for electrocatalysis. Accumulating cations at electrochemical interfaces can alter the proton transfer rate and then tune electrocatalytic performance. However, the mechanism for regulating proton transfer remains ambiguous. Here, we quantify the cation effect on proton diffusion in solution by hydrogen evolution on microelectrodes, revealing the rate can be suppressed by more than 10 times. Different from the prevalent opinions that proton transport is slowed down by modified electric field, we found water structure imposes a more evident effect on kinetics. FTIR test and path integral molecular dynamics simulation indicate that proton prefers to wander within the hydration shell of cations rather than to hop rapidly along water wires. Low connectivity of water networks disrupted by cations corrupts the fast-moving path in bulk water. This study highlights the promising way for regulating proton kinetics via a modified water structure.

Temperature Dependence of Hydrogen Bond Networks of Liquid Water: Thermodynamic Properties and Structural Heterogeneity from Topological Descriptors

Topological analyses of hydrogen bond networks were performed based on the complex network and island statistics of liquid water at different temperatures. The influence of temperature on the liquid water structures and the topological properties of the hydrogen bond networks was investigated by Metropolis Monte Carlo simulations with the TIP4P/2005 potential model. The bilinear behavior of the second peak in the radial distribution function with the temperature was properly reproduced by these simulations. The average connectivity also displayed a bilinear behavior consistent with being a local descriptor. The semiglobal average path length (or geodesic distance) descriptor showed an unprecedented trimodal distribution, whose areas were dependent on the temperature. Considering equilibrium between these three sets of networks, standard enthalpy and entropy of equilibrium were determined for the first time, providing new insights into the structural heterogeneities of liquid water with interesting perspectives for modeling these quantitative properties of hydrogen bond networks.

Water clusters and density fluctuations in liquid water based on extended hierarchical clustering methods

The microscopic structures of liquid water at ambient temperatures remain a hot debate, which relates with structural and density fluctuations in the hydrogen bond network. Here, we use molecular dynamics simulations of liquid water to study the properties of three-dimensional cage-like water clusters, which we investigate using extended graph-based hierarchical clustering methods. The water clusters can cover over 95% of hydrogen bond network, among which some clusters maximally encompass thousands of molecules extending beyond 3.0 nm. The clusters imply fractal behaviors forming percolating networks and the morphologies of small and large clusters show different scaling rules. The local favored clusters and the preferred connections between adjacent clusters correspond to lower energy and conformational entropy depending on cluster topologies. Temperature can destroy large clusters into small ones. We show further that the interior of clusters favors high-density patches. The water molecules in the small clusters, inside which are the void regarded as hydrophobic objects, have a preference for being more tetrahedral. Our results highlight the properties and changes of water clusters as the fundamental building blocks of hydrogen bond networks. In addition, the water clusters can elucidate structural and density fluctuations on different length scales in liquid water.

Effect of Cage Occupancy on Stability and Decomposition of Methane Hydrate

Gas hydrates usually contain a certain number of empty cages that will both affect the hydrate stability and reduce the gas storage capacity. In this work, by MD simulations, we found that the hydrate stability is related to the cage occupancy, the empty cage types, and especially the distribution of empty cages. With the decrease of overall occupancy, the stability of hydrate becomes worse. Under the same overall occupancy, the more concentrated the empty cages are, the more unstable the hydrate is and hence the faster it decomposes. The methane molecules may migrate between distorted cages during the decomposition, resulting in a temporary increase in the stability of hydrate. Hydrates with different empty cage distributions show different decomposition mechanisms: when empty cages are concentrated, the melting rate is fast first due to the rapid decomposition of empty cages, but the remaining filled cages will reduce the melting rate; when empty cages are separated on the contrary, the early melting is slow because of the high local occupancy, and the following melting will be accelerated because of the high melting surface area. It indicates that the empty cage distribution plays a controlling role in hydrate decomposition kinetics at different stages.

A Machine Learning Model to Classify Dynamic Processes in Liquid Water*

The dynamics of water molecules plays a vital role in understanding water. We combined computer simulation and deep learning to study the dynamics of H-bonds between water molecules. Based on ab initio molecular dynamics simulations and a newly defined directed Hydrogen (H-) bond population operator, we studied a typical dynamic process in bulk water: interchange, in which the H-bond donor reverses roles with the acceptor. By designing a recurrent neural network-based model, we have successfully classified the interchange and breakage processes in water. We have found that the ratio between them is approximately 1 : 4, and it hardly depends on temperatures from 280 to 360 K. This work implies that deep learning has the great potential to help distinguish complex dynamic processes containing H-bonds in other systems.

Molecular Origins of Deformation in Amorphous Methane Hydrates

Water and methane can stay together under low temperature and high pressure in the forms of liquid solutions and crystalline solids. From liquid and gaseous states to crystalline solids or the contrary processes, amorphous methane hydrates can occur in these evolution scenarios. Herein, mechanical properties of amorphous methane hydrates are explored for the first time to bridge the gap between mechanical responses of monocrystalline and polycrystalline methane hydrates. Our results demonstrate that mechanical properties of amorphous methane hydrates are strongly governed by our original proposed order parameter, namely, normalized hydrogen-bond directional order parameter. Followed by this important achievement, a multistep deformation mechanism core is proposed to explain mechanical properties of amorphous methane hydrates. Through an extensive detailed analysis of amorphous methane hydrates, our simulation results not only greatly enlarge our fundamental understanding for mechanical responses of amorphous methane hydrates in geological systems but also offer a fresh perspective in structure-property topics of solid materials in future science and technology.

Efficient Structural Relaxation of Polycrystalline Graphene Models

Large samples of experimentally produced graphene are polycrystalline. For the study of this material, it helps to have realistic computer samples that are also polycrystalline. A common approach to produce such samples in computer simulations is based on the method of Wooten, Winer, and Weaire, originally introduced for the simulation of amorphous silicon. We introduce an early rejection variation of their method, applied to graphene, which exploits the local nature of the structural changes to achieve a significant speed-up in the relaxation of the material, without compromising the dynamics. We test it on a 3200 atoms sample, obtaining a speed-up between one and two orders of magnitude. We also introduce a further variation called early decision specifically for relaxing large samples even faster, and we test it on two samples of 10,024 and 20,000 atoms, obtaining a further speed-up of an order of magnitude. Furthermore, we provide a graphical manipulation tool to remove unwanted artifacts in a sample, such as bond crossings.

A hierarchical clustering method of hydrogen bond networks in liquid water undergoing shear flow

Many properties of water, such as turbulent flow, are closely related to water clusters, whereas how water clusters form and transform in bulk water remains unclear. A hierarchical clustering method is introduced to search out water clusters in hydrogen bonded network based on modified Louvain algorithm of graph community. Hydrogen bonds, rings and fragments are considered as 1st-, 2nd-, and 3rd-level structures, respectively. The distribution, dynamics and structural characteristics of 4th- and 5th-level clusters undergoing non-shear- and shear-driven flow are also analyzed at various temperatures. At low temperatures, nearly 50% of water molecules are included in clusters. Over 60% of clusters remain unchanged between neighboring configurations. Obvious collective translational motion of clusters is observed. The topological difference for clusters is elucidated between the inner layer, which favors 6-membered rings, and the external surface layer, which contains more 5-membered rings. Temperature and shearing can not only accelerate the transformation or destruction of clusters at all levels but also change cluster structures. The assembly of large clusters can be used to discretize continuous liquid water to elucidate the properties of liquid water.

Iterative Cup Overlapping: An Efficient Identification Algorithm for Cage Structures of Amorphous Phase Hydrates

Molecular dynamics studies have revealed that the nucleation pathway of clathrate hydrates involves the evolution from amorphous to crystalline hydrates. In this study, complete cages are further classified into the standard edge-saturated cages (SECs) and nonstandard edge-saturated cages (non-SECs). Centered on studying the structure and evolution of non-SECs and SECs, we propose a novel and efficient algorithm, iterative cup overlapping (ICO), to monitor hydrate nucleation and growth in molecular simulations by identifying SECs and discuss possible causes of the instability of non-SECs. Manipulation of topological information makes it possible for ICO to avoid the repeated searches for identified cages and deduce all SECs with low time costs, improving the efficiency of identification significantly. The accuracy and efficiency of ICO were verified by comparing the identification results with other well-proven algorithms. Furthermore, it was found that non-SECs have short lifetimes and eventually decompose or reorganize into more stable structures. Some evidence suggests that the instability of non-SECs is closely related to the hydrogen-bonding configuration of water-ring aggregations that they contain. The spontaneous evolution of the hydrogen-bonding network into the tetrahedral network may be the main factor that causes the conversion of QWRAs and the evolution of non-SECs.

Unraveling nucleation pathway in methane clathrate formation

Methane clathrates are widespread on the ocean floor of the Earth. A better understanding of methane clathrate formation has important implications for natural-gas exploitation, storage, and transportation. A key step toward understanding clathrate formation is hydrate nucleation, which has been suggested to involve multiple evolution pathways. Herein, a unique nucleation/growth pathway for methane clathrate formation has been identified by analyzing the trajectories of large-scale molecular dynamics (MD) simulations. In particular, ternary water-ring aggregations (TWRAs) have been identified as fundamental structures for characterizing the nucleation pathway. Based on this nucleation pathway, the critical nucleus size and nucleation timescale can be quantitatively determined. Specifically, a methane hydration layer compression/shedding process is observed to be the critical step in (and driving) the nucleation/growth pathway, which is manifested through overlapping/compression of the surrounding hydration layers of the methane molecules, followed by detachment (shedding) of the hydration layer. As such, an effective way to control methane hydrate nucleation is to alter the hydration layer compression/shedding process during the course of nucleation.

Growth and arrest of topological cycles in small physical networks

The chordless cycle sizes of spatially embedded networks are demonstrated to follow an exponential growth law similar to random graphs if the number of nodes [Formula: see text] is below a critical value [Formula: see text] For covalent polymer networks, increasing the network size, as measured by the number of cross-link nodes, beyond [Formula: see text] results in a crossover to a new regime in which the characteristic size of the chordless cycles [Formula: see text] no longer increases. From this result, the onset and intensity of finite-size effects can be predicted from measurement of [Formula: see text] in large networks. Although such information is largely inaccessible with experiments, the agreement of simulation results from molecular dynamics, Metropolis Monte Carlo, and kinetic Monte Carlo suggests the crossover is a fundamental physical feature which is insensitive to the details of the network generation. These results show random graphs as a promising model to capture structural differences in confined physical networks.

Inhibition Effect of Kinetic Hydrate Inhibitors on the Growth of Methane Hydrate in Gas–Liquid Phase Separation State

The effect of kinetic hydrate inhibitors (KHIs) on the growth of methane hydrate in the gas–liquid phase separation state is studied at the molecular level. The simulation results show that the kinetic inhibitors, named PVP and PVP-A, show good inhibitory effects on the growth of methane hydrate under the gas–liquid phase separation state, and the initial position of the kinetic hydrate inhibitors has a major effect on the growth of methane hydrates. In addition, inhibitors at different locations exhibit different inhibition performances. When the inhibitor molecules are located at the gas–liquid phase interface, increasing the contact area between the groups of the inhibitor molecules and methane is beneficial to enhance the inhibitory performance of the inhibitors. When inhibitor molecules are located at the solid–liquid phase interface, the inhibitor molecules adsorbed on the surface of the hydrate nucleus and decreased the direct contact of hydrate nucleus with the surrounding water and methane molecules, which would delay the growth of hydrate nucleus. Moreover, the increase of hydrate surface curvature and the Gibbs–Thomson effect caused by inhibitors can also reduce the growth rate of methane hydrate.

Topological Water Network Analysis Around Amino Acids

Water molecules play a key role in protein stability, folding, function and ligand binding. Protein hydration has been studied using free energy perturbation algorithms. However, the study of protein hydration without free energy calculation is also an active field of research. Accordingly, topological water network (TWN) analysis has been carried out instead of free energy calculation in the present work to investigate hydration of proteins. Water networks around 20 amino acids in the aqueous solution were explored through molecular dynamics (MD) simulations. These simulation results were compared with experimental observations. Water molecules from the protein data bank structures showed TWN patterns similar to MD simulations. This work revealed that TWNs are effected by the surrounding environment. TWNs could provide valuable clues about the environment around amino acid residues in the proteins. The findings from this study could be exploited for TWN-based drug discovery and development.

Functional Hydration Behavior: Interrelation between Hydration and Molecular Properties at Lipid Membrane Interfaces

Water is an abundant commodity and has various important functions. It stabilizes the structure of biological macromolecules, controls biochemical activities, and regulates interfacial/intermolecular interactions. Common aspects of interfacial water can be obtained by overviewing fundamental functions and properties at different temporal and spatial scales. It is important to understand the hydrogen bonding and structural properties of water and to evaluate the individual molecular species having different hydration properties. Water molecules form hydrogen bonds with biomolecules and contribute to the adjustment of their properties, such as surface charge, hydrophilicity, and structural flexibility. In this review, the fundamental properties of water molecules and the methods used for the analyses of water dynamics are summarized. In particular, the interrelation between the hydration properties, determined by molecules, and the properties of molecules, determined by their hydration properties, are discussed using the lipid membrane as an example. Accordingly, interesting water functions are introduced that provide beneficial information in the fields of biochemistry, medicine, and food chemistry.

Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core

The primary driving force for protein folding is the formation of a well-packed, anhydrous core. However, recently, the crystal structure of an antifreeze protein, maxi, has been resolved where the core of the protein is filled with water, which apparently contradicts the existing notion of protein folding. Here, we have performed standard molecular dynamics (MD) simulation, replica exchange MD (REMD) simulation, and umbrella sampling using TIP4P water at various temperatures (300, 260, and 240 K) to explore the origin of this unusual structural feature. It is evident from standard MD and REMD simulations that the protein is found to be stable at 240 K in its unusual state. The core of protein has two layers of semi-clathrate water separating the methyl groups of alanine residues from different helical strands. However, with increasing temperature (260 and 300 K), the stability decreases as the core becomes dehydrated, and methyl groups of alanine are tightly packed driven by hydrophobic interactions. Calculation of the potential of mean force by an umbrella sampling technique between a pair of model hydrophobes resembling maxi protein at 240 K shows the stabilization of second solvent-separated minima (SSM), which provides a thermodynamic rationale of the unusual structural feature in terms of weakening of the hydrophobic interaction. Because the stabilization of SSMs is implicated for cold denaturation, it suggests that the maxi protein is so designed by nature where the cold denatured-like state becomes the biologically active form as it works near or below the freezing point of water.

Methane Hydrate Nucleation within Elastic Confined Spaces: Suitable Spacing and Elasticity Can Accelerate the Nucleation

Elastic materials are candidates for process intensification of gas storage by forming gas hydrate. In this work, molecular dynamics simulations of hydrate nucleation in elastic silica double layers were performed to study the effect of elastic confined spaces on hydrate formation. It is found that in narrow confined spaces, hexagonal rings dominated the hydrogen bond network of water molecules established rapidly by a multisite nucleation mechanism. With molecules added, a bilayer water structure was formed finally because elastic space can adapt the volume expansion. In medium and wide confined spaces, hydrates were formed from a series of "pseudo cages" which are considered as precursors of complete hydrate cages. Moreover, the induction time for nucleation was a minimum when the elasticity of the silica layer changes: nucleation is fastest in the weak-elastic system. When the elasticity increases, it becomes hard to adapt the volume expansion during nucleation and also difficult to nucleate in very weak-elastic systems because of the fluctuation of the layers.

Formation of a nanobubble and its effect on the structural ordering of water in a CH4-N2-CO2-H2O mixture

The replacement of methane (CH4) from its hydrate by a mixture of nitrogen (N2) and carbon dioxide (CO2) involves the dissociation of methane hydrate leading to the formation of a CH4-N2-CO2-H2O mixture that can significantly influence the subsequent steps of the replacement process. In the present work, we study the evolution of dissolved gas molecules in this mixture by applying classical molecular dynamics simulations. Our study shows that a higher CO2 : N2 ratio in the mixture enhances the formation of nanobubbles composed of N2, CH4 and CO2 molecules. To understand how the CO2 : N2 ratio affects nanobubble nucleation, the distribution of molecules in the bubble formed is examined. It is observed that unlike N2 and CH4, the density of CO2 in the bubble reaches a maximum at the surface of the bubble. The accumulation of CO2 molecules at the surface makes the bubble more stable by decreasing the excess pressure inside the bubble as well as surface tension at its interface with water. It is found that a frequent exchange of gas molecules takes place between the bubble and the surrounding liquid and an increase in concentration of CO2 in the mixture leads to a decrease in the number of such exchanges. The effect of nanobubbles on the structural ordering of water molecules is examined by determining the number of water rings formed per unit volume in the mixture. The role of nanobubbles in water structuring is correlated to the dynamic nature of the bubble arising from the exchange of gas molecules between the bubble and the liquid.

Pressure increases the ice-like order of water at hydrophobic interfaces

Ice-like order of water at hydrophobic interfaces is quantified on different length scales based on molecular dynamics simulations.

Enhanced heterogeneous ice nucleation by special surface geometry

The freezing of water typically proceeds through impurity-mediated heterogeneous nucleation. Although non-planar geometry generically exists on the surfaces of ice nucleation centres, its role in nucleation remains poorly understood. Here we show that an atomically sharp, concave wedge can further promote ice nucleation with special wedge geometries. Our molecular analysis shows that significant enhancements of ice nucleation can emerge both when the geometry of a wedge matches the ice lattice and when such lattice match does not exist. In particular, a 45° wedge is found to greatly enhance ice nucleation by facilitating the formation of special topological defects that consequently catalyse the growth of regular ice. Our study not only highlights the active role of defects in nucleation but also suggests that the traditional concept of lattice match between a nucleation centre and crystalline lattice should be extended to include a broader match with metastable, non-crystalline structural motifs.

Understanding ice nucleation is important for the development of accurate cloud models. Here Bi et al. show that sharp wedges can enhance ice nucleation both when the wedge geometry matches the ice lattice and when such matching is absent, in which case nucleation is promoted by topological defects.

Giant Pockels effect of polar organic solvents and water in the electric double layer on a transparent electrode

The Pockels effect of polar organic solvents and water within the electric double layer on an ITO electrode is studied to find that water has the largest Pockels coefficient, followed in order by methanol, ethanol, and dimethyl sulfoxide.

Phase behavior of a binary fluid mixture of quadrupolar molecules

We propose a model molecule to investigate microscopic properties of a binary mixture with a closed-loop coexistence region. The molecule is comprised of a Lennard-Jones particle and a uniaxial quadrupole. Gibbs ensemble Monte Carlo simulations demonstrate that the high-density binary fluid of the molecules with the quadrupoles of the same magnitude but of the opposite signs can show closed-loop immiscibility. We find that an increase in the magnitude of the quadrupoles causes a shrinkage of the coexistence region. Molecular dynamics simulations also reveal that aggregates with two types of molecules arranged alternatively are formed in the stable one-phase region both above and below the coexistence region. String structures are dominant below the lower critical solution temperature, while branched aggregates are observed above the upper critical solution temperature. We conclude that the anisotropic interaction between the quadrupoles of the opposite signs plays a crucial role in controlling these properties of the phase behavior.

Chiral Ordering in Supercooled Liquid Water and Amorphous Ice

The emergence of homochiral domains in supercooled liquid water is presented using molecular dynamics simulations. An individual water molecule possesses neither a chiral center nor a twisted conformation that can cause spontaneous chiral resolution. However, an aggregation of water molecules will naturally give rise to a collective chirality. Such homochiral domains possess obvious topological and geometrical orders and are energetically more stable than the average. However, homochiral domains cannot grow into macroscopic homogeneous structures due to geometrical frustrations arising from their icosahedral local order. Homochiral domains are the major constituent of supercooled liquid water and the origin of heterogeneity in that substance, and are expected to be enhanced in low-density amorphous ice at lower temperatures.

Origin of Self-preservation Effect for Hydrate Decomposition: Coupling of Mass and Heat Transfer Resistances

Gas hydrates could show an unexpected high stability at conditions out of thermodynamic equilibrium, which is called the self-preservation effect. The mechanism of the effect for methane hydrates is here investigated via molecular dynamics simulations, in which an NVT/E method is introduced to represent different levels of heat transfer resistance. Our simulations suggest a coupling between the mass transfer resistance and heat transfer resistance as the driving mechanism for self-preservation effect. We found that the hydrate is initially melted from the interface, and then a solid-like water layer with temperature-dependent structures is formed next to the hydrate interface that exhibits fractal feature, followed by an increase of mass transfer resistance for the diffusion of methane from hydrate region. Furthermore, our results indicate that heat transfer resistance is a more fundamental factor, since it facilitates the formation of the solid-like layer and hence inhibits the further dissociation of the hydrates. The self-preservation effect is found to be enhanced with the increase of pressure and particularly the decrease of temperature. Kinetic equations based on heat balance calculations is also developed to describe the self-preservation effect, which reproduces our simulation results well and provides an association between microscopic and macroscopic properties.

How Properties of Solid Surfaces Modulate the Nucleation of Gas Hydrate

Molecular dynamics simulations were performed for CO2 dissolved in water near silica surfaces to investigate how the hydrophilicity and crystallinity of solid surfaces modulate the local structure of adjacent molecules and the nucleation of CO2 hydrates. Our simulations reveal that the hydrophilicity of solid surfaces can change the local structure of water molecules and gas distribution near liquid-solid interfaces, and thus alter the mechanism and dynamics of gas hydrate nucleation. Interestingly, we find that hydrate nucleation tends to occur more easily on relatively less hydrophilic surfaces. Different from surface hydrophilicity, surface crystallinity shows a weak effect on the local structure of adjacent water molecules and on gas hydrate nucleation. At the initial stage of gas hydrate growth, however, the structuring of molecules induced by crystalline surfaces are more ordered than that induced by amorphous solid surfaces.

Nucleation rate analysis of methane hydrate from molecular dynamics simulations

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

Order Parameters and Algorithmic Approaches for Detection and Demarcation of Interfaces in Hydrate-Fluid and Ice-Fluid Systems

Some aspects of the use of order parameter fields in molecular dynamics simulations to delimit solid phases containing water, namely ice and hydrate, in both hydrophilic and hydrophobic fluids are examined; this includes the influences of rectangular meshes and of filtering on the quality of these parameters. Three order parameters are studied: the mass density, ρ; an angular tetrahedrality measure, Sg (Chau and Hardwick, Mol. Phys. 1998, 93, 511); and the water-dimer dihedral angle, F4 (Rodger et al. Fluid Phase Equilib. 1996, 116, 326). The parameters are studied to find their ability to distinguish between bulk phases, their consistency in different environments, their noise susceptibility, and their ability to demarcate the interface region. Spatial sampling and filtering are covered in detail, and some temporal features are illustrated by using autocorrelation maps. The parameters are employed to determine the position of interfaces as functions of time and, with the capillary wave fluctuation method (Hoyt et al. Phys. Rev. Lett. 2001, 86, 5530; Math. Comput. Simul. 2010, 80, 1382), to estimate solid-fluid interfacial stiffnesses, with partial success for the hydrophilic/hydrophobic-type interfaces.

Reaction coordinate of incipient methane clathrate hydrate nucleation

Nucleation from solution is a ubiquitous phenomenon with relevance to myriad scientific disciplines, including pharmaceuticals, biomineralization, and disease. One prominent example is the nucleation of clathrate hydrates, multicomponent crystalline inclusion compounds relevant to the energy industry where they block pipelines and also constitute a potential vast energy resource. Despite their importance, the molecular mechanism of incipient hydrate formation remains unknown. Herein, we employ advanced molecular simulation tools (pB histogram, equilibrium path sampling) to provide a statistical-mechanical basis for extracting physical insight into the molecular steps by which clathrates form. Through testing the Mutually Coordinated Guest (MCG) order parameter, we demonstrate that both guest (methane) and host (water) structuring are crucial to accurately describe the nucleation of hydrates and determine a critical nucleus size of MCG-1 = 16 at 255 K and 500 bar. Equipped with a validated (and novel) reaction coordinate, subsequent equilibrium path sampling simulations yield the free energy barrier and nucleation rate. The resulting quantitative nucleation process is described by the MCG clustering mechanism. This constitutes a significant advance in the field of hydrates research, as the fitness of a molecular descriptor has never been statistically verified. More broadly, this work has significance to a wide range of multicomponent nucleation contexts wherein the formation mechanism depends on contributions from both solute and solvent.

Two-component order parameter for quantifying clathrate hydrate nucleation and growth

Methane clathrate hydrate nucleation and growth is investigated via analysis of molecular dynamics simulations using a new order parameter. This order parameter (OP), named the Mutually Coordinated Guest (MCG) OP, quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. It is the first two-component OP used for quantifying hydrate nucleation and growth. The algorithm for calculating the MCG OP is described in detail. Its physical motivation and advantages compared to existing methods are discussed.

Proton transfer through the water gossamer

The diffusion of protons through water is understood within the framework of the Grotthuss mechanism, which requires that they undergo structural diffusion in a stepwise manner throughout the water network. Despite long study, this picture oversimplifies and neglects the complexity of the supramolecular structure of water. We use first-principles simulations and demonstrate that the currently accepted picture of proton diffusion is in need of revision. We show that proton and hydroxide diffusion occurs through periods of intense activity involving concerted proton hopping followed by periods of rest. The picture that emerges is that proton transfer is a multiscale and multidynamical process involving a broader distribution of pathways and timescales than currently assumed. To rationalize these phenomena, we look at the 3D water network as a distribution of closed directed rings, which reveals the presence of medium-range directional correlations in the liquid. One of the natural consequences of this feature is that both the hydronium and hydroxide ion are decorated with proton wires. These wires serve as conduits for long proton jumps over several hydrogen bonds.

Insights into the ultraviolet spectrum of liquid water from model calculations: the different roles of donor and acceptor hydrogen bonds in water pentamers

With a view toward a better understanding of changes in the peak position and shape of the first absorption band of water with condensation or temperature, results from electronic structure calculations using high level wavefunction based and time-dependent density functional methods are reported for water pentamers. Excitation energies, oscillator strengths, and redistributions of electron density are determined for the quasitetrahedral water pentamer in its C(2v) equilibrium geometry and for many pentamer configurations sampled from molecular simulation of liquid water. Excitations associated with surface molecules are removed in order to focus on those states associated with the central molecule, which are the most representative of the liquid environment. The effect of hydrogen bonding on the lowest excited state associated with the central molecule is studied by adding acceptor or donor hydrogen bonds to tetramer and trimer substructures of the C(2v) pentamer, and by sampling liquid-like configurations having increasing number of acceptor or donor hydrogen bonds of the central molecule. Our results provide clear evidence that the blueshift of excitation energies upon condensation is essentially determined by acceptor hydrogen bonds, and the magnitudes of these shifts are determined by the number of such, whereas donor hydrogen bonds do not induce significant shifts in excitation energies. This qualitatively different role of donor and acceptor hydrogen bonds is understood in terms of the different roles of the 1b(1) monomer molecular orbitals, which establishes an intimate connection between the valence hole and excitation energy shifts. Since the valence hole of the lowest excitation associated with the central molecule is found to be well localized in all liquid-like hydrogen bonding environments, with an average radius of gyration of ~1.6 Å that is much lower than the nearest neighbor O-O distance, a clear and unambiguous connection between hydrogen bonding environments and excitation energy shifts can be established. Based on these results, it is concluded that peak position of the first absorption band is mainly determined by the relative distribution of single and double acceptor hydrogen bonding environments, whereas the shape of the first absorption band is mainly determined by the relative distribution of acceptor and broken acceptor hydrogen bonding environments. The temperature dependence of the peak position and shape of the first absorption band can be readily understood in terms of changes to these relative populations.

Nucleation of the CO2 hydrate from three-phase contact lines

Using molecular dynamics simulations on the microsecond time scale, we investigate the nucleation and growth mechanisms of CO(2) hydrates in a water/CO(2)/silica three-phase system. Our simulation results indicate that the CO(2) hydrate nucleates near the three-phase contact line rather than at the two-phase interfaces and then grows along the contact line to form an amorphous crystal. In the nucleation stage, the hydroxylated silica surface can be understand as a stabilizer to prolong the lifetime of adsorbed hydrate cages that interact with the silica surface by hydrogen bonding, and the adsorbed cages behave as the nucleation sites for the formation of an amorphous CO(2) hydrate. After nucleation, the nucleus grows along the three-phase contact line and prefers to develop toward the CO(2) phase as a result of the hydrophilic nature of the modified solid surface and the easy availability of CO(2) molecules. During the growth process, the population of sI cages in the formed amorphous crystal is found to increase much faster than that of sII cages, being in agreement with the fact that only the sI hydrate can be formed in nature for CO(2) molecules.