No Arabic abstract
Understanding the phase behaviors of nanoconfined water has driven notable research interests recently. In this work, we examine the structures and thermodynamics of water encapsulated under a graphene cover. We find layered water structures up to ~1000 molecules, which is stabilized by the spatial confinement and pressure induced by the adhesion between graphene cover and substrate. For monolayer encapsulations, we identify both crystalline lattices and defects. Free energy analysis shows that these low- entropy orders are compensated by high formation energies. There exists an order- disorder transition for this condensed phase at ~480-490 K, with a sharp reduction in the number of hydrogen bonds and increase in the entropy. These findings offer fundamental understandings of the encapsulated water, and provide guidance for practical applications with its presence, for example, in the design of nanoelectronic devices.
We evaluate the grand potential of a cluster of two molecular species, equivalent to its free energy of formation from a binary vapour phase, using a nonequilibrium molecular dynamics technique where guide particles, each tethered to a molecule by a harmonic force, move apart to disassemble a cluster into its components. The mechanical work performed in an ensemble of trajectories is analysed using the Jarzynski equality to obtain a free energy of disassembly, a contribution to the cluster grand potential. We study clusters of sulphuric acid and water at 300 K, using a classical interaction scheme, and contrast two modes of guided disassembly. In one, the cluster is broken apart through simple pulling by the guide particles, but we find the trajectories tend to be mechanically irreversible. In the second approach, the guide motion and strength of tethering are modified in a way that prises the cluster apart, a procedure that seems more reversible. We construct a surface representing the cluster grand potential, and identify a critical cluster for droplet nucleation under given vapour conditions. We compare the equilibrium populations of clusters with calculations reported by Henschel et al. [J. Phys. Chem. A 118, 2599 (2014)] based on optimised quantum chemical structures.
Deeply supercooled water exhibits complex dynamics with large density fluctuations, ice coarsening and characteristic time scales extending from picoseconds to milliseconds. Here, we discuss implications of these time scales as they pertain to two-phase coexistence and to molecular simulations of supercooled water. Specifically, we argue that it is possible to discount liquid-liquid criticality because the time scales imply that correlation lengths for such behavior would be bounded by no more than a few nanometers. Similarly, it is possible to discount two-liquid coexistence because the time scales imply a bounded interfacial free energy that cannot grow in proportion to a macroscopic surface area. From time scales alone, therefore, we see that coexisting domains of differing density in supercooled water can be no more than nano-scale transient fluctuations.
Using concepts from perturbation and local molecular field theories of liquids we divide the potential of the SPC/E water model into short and long ranged parts. The short ranged parts define a minimal reference network model that captures very well the structure of the local hydrogen bond network in bulk water while ignoring effects of the remaining long ranged interactions. This deconstruction can provide insight into the different roles that the local hydrogen bond network, dispersion forces, and long ranged dipolar interactions play in determining a variety of properties of SPC/E and related classical models of water. Here we focus on the anomalous behavior of the internal pressure and the temperature dependence of the density of bulk water. We further utilize these short ranged models along with local molecular field theory to quantify the influence of these interactions on the structure of hydrophobic interfaces and the crossover from small to large scale hydration behavior. The implications of our findings for theories of hydrophobicity and possible refinements of classical water models are also discussed.
Water modeling is a challenging problem. Its anomalies are difficult to reproduce, promoting the proliferation of a large number of computational models, among which researchers select the most appropriate for the property they study. In this chapter, we introduce a coarse-grained model introduced by Franzese and Stanley (FS) that accounts for the many-body interactions of water. We review mean-field calculations and Monte Carlo simulations on water monolayers for a wide range of pressures and temperatures, including extreme conditions. The results show the presence of two dynamic crossovers and explain the origin of diffusion anomalies. Moreover, the model shows that all the different scenarios, proposed in the last decades as alternative explanations of the experimental anomalies of water, can be related by the fine-tuning of the many-body (cooperative) interaction. Once this parameter is set from the experiments, the FS model predicts a phase transition between two liquids with different densities and energies in the supercooled water region, ending in a liquid-liquid critical point. From this critical point stems a liquid-liquid Widom line, i.e., the locus of maxima of the water correlation length, that in the FS model can be directly calculated. The results are consistent with the extrapolations from experiments. Furthermore, they agree with those from atomistic models but make predictions over a much wider thermodynamic region, allowing for a better interpretation of the available experimental data. All these findings provide a coherent picture of the properties of water and confirm the validity of the FS model that has proved to be useful for large-scale simulations of biological systems.
The free energetics of water density fluctuations near a surface, and the rare low-density fluctuations in particular, serve as reliable indicators of surface hydrophobicity; the easier it is to displace the interfacial waters, the more hydrophobic the underlying surface. However, characterizing the free energetics of such rare fluctuations requires computationally expensive, non-Boltzmann sampling methods like umbrella sampling. This inherent computational expense associated with umbrella sampling makes it challenging to investigate the role of polarizability or electronic structure effects in influencing interfacial fluctuations. Importantly, it also limits the size of the volume, which can be used to probe interfacial fluctuations. The latter can be particularly important in characterizing the hydrophobicity of large surfaces with molecular-level heterogeneities, such as those presented by proteins. To overcome these challenges, here we present a method for the sparse sampling of water density fluctuations, which is roughly two orders of magnitude more efficient than umbrella sampling. We employ thermodynamic integration to estimate the free energy differences between biased ensembles, thereby circumventing the umbrella sampling requirement of overlap between adjacent biased distributions. Further, a judicious choice of the biasing potential allows such free energy differences to be estimated using short simulations, so that the free energetics of water density fluctuations are obtained using only a few, short simulations. Leveraging the efficiency of the method, we characterize water density fluctuations in the entire hydration shell of the protein, ubiquitin; a large volume containing an average of more than six hundred waters.