No Arabic abstract
Four scenarios have been proposed for the low--temperature phase behavior of liquid water, each predicting different thermodynamics. The physical mechanism which leads to each is debated. Moreover, it is still unclear which of the scenarios best describes water, as there is no definitive experimental test. Here we address both open issues within the framework of a microscopic cell model by performing a study combining mean field calculations and Monte Carlo simulations. We show that a common physical mechanism underlies each of the four scenarios, and that two key physical quantities determine which of the four scenarios describes water: (i) the strength of the directional component of the hydrogen bond and (ii) the strength of the cooperative component of the hydrogen bond. The four scenarios may be mapped in the space of these two quantities. We argue that our conclusions are model-independent. Using estimates from experimental data for H bond properties the model predicts that the low-temperature phase diagram of water exhibits a liquid--liquid critical point at positive pressure.
Over the years, plenty of classical interaction potentials for water have been developed and tested against structural, dynamical and thermodynamic properties. On the other hands, it has been recently observed (F. Martelli et. al, textit{ACS Nano}, textbf{14}, 8616--8623, 2020) that the topology of the hydrogen bond network (HBN) is a very sensitive measure that should be considered when developing new interaction potentials. Here we report a thorough comparison of 11 popular non polarizable classical water models against their HBN, which is at the root of water properties. We probe the topology of the HBN using the ring statistics and we evaluate the quality of the network inspecting the percentage of broken and intact HBs. For each water model, we assess the tendency to develop hexagonal rings (that promote crystallization at low temperatures) and pentagonal rings (known to frustrate against crystallization at low temperatures). We then introduce the emph{network complexity index}, a general descriptor to quantify how much the topology of a given network deviates from that of the ground state, namely of hexagonal or cubic ice. Remarkably, we find that the network complexity index allows us to relate, for the first time, the dynamical properties of different water models with their underlying topology of the HBN. Our study provides a benchmark against which the performances of new models should be tested against, and introduces a general way to quantify the complexity of a network which can be transferred to other materials and that links the topology of the HBN with dynamical properties. Finally, our study introduces a new perspective that can help in rationalizing the transformations among the different phases of water and of other materials.
During compression of a water dimer calculated with high-precision first-principles methods, the trends of H-bond and O-H bond lengths show quantum effect of the electronic structure. We found that the H-bond length keeps decreasing, while the O-H bond length increases up to the stable point and decreases beyond it when the water dimer is further compressed. The remarkable properties are different from those observed in most previous researches which can be understood and extrapolated through classical simulation. The observations can be explained by the decrease in orbital overlap and change in the exchange repulsion interaction between water monomers. The dominant interaction between water monomers changes from electrostatic interaction to exchange repulsion at the turning point of the O-H bond length when the O...O distance is decreased. These findings highlight the quantum effect on the hydrogen bond in extreme conditions and play an important role in the recognition of the hydrogen bond structure and mechanism.
Numerical simulations for a number of water models have supported the possibility of a metastable liquid-liquid critical point (LLCP) in the deep super-cooled region. Here we consider a theoretical model for a supercooled liquid water monolayer and its mathematical mapping onto a percolation problem. The mapping allows us to identify the finite-size clusters at any state-point, and the infinite cluster at the critical point, with the regions of correlated hydrogen bonds (HBs). We show that the percolation line coincides with the first-order liquid-liquid phase transition ending at the LLCP. At pressures below the LLCP, the percolation line corresponds to the strong maxima of the thermodynamic response functions and to the locus of maximum correlation length (Widom line). At higher pressures, we find a percolation transition with a positive slope and we discuss its possible relation with the thermodynamics.
We have carried out a density functional theory study on the structures of DMSO clusters and analysed the structure and their stability using molecular electrostatic potential and quantum theory of atoms-in-molecules (QTAIM). The ground state geometry of the DMSO clusters, prefer to exist in ouroboros shape. Pair wise interaction energy calculation show the interaction between methyl groups of adjacent DMSO molecules and a destabilization is is created by the methyl groups which are away from each other. Molecular electrostatic potential analysis shows the existence of hole on the odd numbered clusters, which helps in their highly directional growth. QTAIM analysis show the existence of two intermolecular hydrogen bonds, of type SOC hydrogen bonds and methyl CHC dihydrogen bonds. The computed and Laplacian values were all positive for the intermolecular bonds, supporting the existence of noncovalent interactions. The computed ellipticity for the dihydrogen bonds have values > 2, which confirms the delocalization of electron, are mainly due to the hydrogen-hydrogen interactions of methyl groups. A plot of total hydrogen bonding energy vs the observed total local electron density shows linearity with correlation coefficient of near unity, which indicates the cooperative effects of intermolecular dihydrogen HH bonds.
The glass transition in hydrogen-bonded glass formers differs from the glass transition in other glass formers. The Eshelby rearrangements of the highly viscous flow are superimposed by strongly asymmetric hydrogen bond rupture processes, responsible for the excess wing. Their influence on the shear relaxation spectrum is strong in glycerol and close to zero in PPE, reflecting the strength of the hydrogen bond contribution to the high frequency shear modulus. An appropriate modification of a recent theory of the highly viscous flow enables a quantitative common description of the relaxation spectra in shear, linear and non-linear dielectrics, and heat capacity.