-
Anomalous behaviour of transport properties in a supercooled water-glycerol mixture
Authors:
Alberto Zaragoza,
Rajat Kumar,
Jose Martín Roca,
Prabir Khatua,
Valeria Molinero,
Frédéric Caupin,
Chantal Valeriani
Abstract:
Glycerol acts as a natural cryoprotectant by depressing the temperature of ice nucleation and slowing down the dynamics of water mixtures. In this work we characterize dynamics -- diffusion, viscosity, and hydrogen-bond dynamics -- as well as density anomaly and structure of water mixtures with 1\% to 50\% w/w glycerol at low temperatures via molecular dynamics simulations using all-atom and coars…
▽ More
Glycerol acts as a natural cryoprotectant by depressing the temperature of ice nucleation and slowing down the dynamics of water mixtures. In this work we characterize dynamics -- diffusion, viscosity, and hydrogen-bond dynamics -- as well as density anomaly and structure of water mixtures with 1\% to 50\% w/w glycerol at low temperatures via molecular dynamics simulations using all-atom and coarse-grained models. Simulations reveal distinct violations of the Stokes-Einsten relation in the low temperature regime for water and glycerol. Deviations are positive for water at all concentrations, and positive for glycerol in very dilute solutions but turning negative in concentrated ones. The all-atom and coarse-grained models reveal an unexpected crossover in the dynamics of the 1% and 10 % w/w glycerol at the lowest simulated temperatures. This crossover manifests in the diffusion coefficients of water and glycerol, as well as in the viscosity and lifetime of hydrogen-bonds in water. We interpret that the crossover originates on the opposing dependence with glycerol concentration of the two factors controlling the solution's slow-down: the increase in tetrahedrally coordinated water and the dynamics and clustering of the glycerol molecules. We anticipate that this dynamic crossover will also occur for solution of water with other polyols.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
CEGANN: Crystal Edge Graph Attention Neural Network for multiscale classification of materials environment
Authors:
Suvo Banik,
Debdas Dhabal,
Henry Chan,
Sukriti Manna,
Mathew Cherukara,
Valeria Molinero,
Subramanian KRS Sankaranarayanan
Abstract:
Machine learning models and applications in materials design and discovery typically involve the use of feature representations or "descriptors" followed by a learning algorithm that maps them to a user-desired property of interest. Most popular mathematical formulation-based descriptors are not unique across atomic environments or suffer from transferability issues across different application do…
▽ More
Machine learning models and applications in materials design and discovery typically involve the use of feature representations or "descriptors" followed by a learning algorithm that maps them to a user-desired property of interest. Most popular mathematical formulation-based descriptors are not unique across atomic environments or suffer from transferability issues across different application domains and/or material classes. In this work, we introduce the Crystal Edge Graph Attention Neural Network (CEGANN) workflow that uses graph attention-based architecture to learn unique feature representations and perform classification of materials across multiple scales (from atomic to mesoscale) and diverse classes ranging from metals, oxides, non-metals and even hierarchical materials such as zeolites and semi ordered materials such as mesophases. We first demonstrate a case study where the classification is based on a global, structure-level representation such as space group and structural dimensionality (e.g., bulk, 2D, clusters etc.). Using representative materials such as polycrystals and zeolites, we next demonstrate the transferability of our network in successfully performing local atom-level classification tasks, such as grain boundary identification and other heterointerfaces. We also demonstrate classification in (thermal) noisy dynamical environments using a representative example of crystal nucleation and growth of a zeolite polymorph from an amorphous synthesis mixture. Finally, we characterize the formation of a binary mesophase and its phase transitions and the growth of ice, demonstrating the performance of CEGANN in systems with thermal noise and compositional diversity. Overall, our approach is agnostic to the material type and allows for multiscale classification of features ranging from atomic-scale crystal structures to heterointerfaces to microscale grain boundaries.
△ Less
Submitted 15 November, 2022; v1 submitted 20 July, 2022;
originally announced July 2022.
-
Nature of the anomalies in the supercooled liquid state of the mW model of water
Authors:
Vincent Holten,
David T. Limmer,
Valeria Molinero,
Mikhail A. Anisimov
Abstract:
The thermodynamic properties of the supercooled liquid state of the mW model of water show anomalous behavior. Like in real water, the heat capacity and compressibility sharply increase upon supercooling. One of the possible explanations of these anomalies, the existence of a second (liquid-liquid) critical point, is not supported by simulations for this model. In this work, we reproduce the anoma…
▽ More
The thermodynamic properties of the supercooled liquid state of the mW model of water show anomalous behavior. Like in real water, the heat capacity and compressibility sharply increase upon supercooling. One of the possible explanations of these anomalies, the existence of a second (liquid-liquid) critical point, is not supported by simulations for this model. In this work, we reproduce the anomalies of the mW model with two thermodynamic scenarios: one based on a non-ideal "mixture" with two different types of local order of the water molecules, and one based on weak crystallization theory. We show that both descriptions accurately reproduce the model's basic thermodynamic properties. However, the coupling constant required for the power laws implied by weak crystallization theory is too large relative to the regular backgrounds, contradicting assumptions of weak crystallization theory. Fluctuation corrections outside the scope of this work would be necessary to fit the forms predicted by weak crystallization theory. For the two-state approach, the direct computation of the low-density fraction of molecules in the mW model is in agreement with the prediction of the phenomenological equation of state. The non-ideality of the "mixture" of the two states never becomes strong enough to cause liquid-liquid phase separation, also in agreement with simulation results.
△ Less
Submitted 1 April, 2013; v1 submitted 22 February, 2013;
originally announced February 2013.
-
The Rise and Fall of Anomalies in Tetrahedral Liquids
Authors:
Waldemar Hujo,
B. Shadrack Jabes,
Varun K. Rana,
Charusita Chakravarty,
Valeria Molinero
Abstract:
The thermodynamic liquid-state anomalies and associated structural changes of the Stillinger-Weber family of liquids are mapped out as a function of the degree of tetrahedrality of the interaction potential, focusing in particular on tetrahedrality values suitable for modeling C, H2O, Si, Ge and Sn. We show that the density anomaly, associated with a rise in molar volume on isobaric cooling, emerg…
▽ More
The thermodynamic liquid-state anomalies and associated structural changes of the Stillinger-Weber family of liquids are mapped out as a function of the degree of tetrahedrality of the interaction potential, focusing in particular on tetrahedrality values suitable for modeling C, H2O, Si, Ge and Sn. We show that the density anomaly, associated with a rise in molar volume on isobaric cooling, emerges at intermediate tetrahedralities (e.g. Ge, Si and H2O) but is absent in the low (e.g. Sn) and high (e.g. C) tetrahedrality liquids. The rise in entropy on isothermal compression associated with the density anomaly is related to the structural changes in the liquid using the pair correlation entropy. An anomalous increase in the heat capacity on isobaric cooling exists at high tetrahedralities but is absent at low tetrahedralities (e.g. Sn). Structurally, this heat capacity anomaly originates in a sharp rise in the fraction of four-coordinated particles and local tetrahedral order in the liquid as its structure approaches that of the tetrahedral crystal.
△ Less
Submitted 28 July, 2011;
originally announced July 2011.
-
Structural transformation in supercooled water controls the crystallization rate of ice
Authors:
Emily B. Moore,
Valeria Molinero
Abstract:
One of water's unsolved puzzles is the question of what determines the lowest temperature to which it can be cooled before freezing to ice. The supercooled liquid has been probed experimentally to near the homogeneous nucleation temperature TH{\approx}232 K, yet the mechanism of ice crystallization - including the size and structure of critical nuclei - has not yet been resolved. The heat capacity…
▽ More
One of water's unsolved puzzles is the question of what determines the lowest temperature to which it can be cooled before freezing to ice. The supercooled liquid has been probed experimentally to near the homogeneous nucleation temperature TH{\approx}232 K, yet the mechanism of ice crystallization - including the size and structure of critical nuclei - has not yet been resolved. The heat capacity and compressibility of liquid water anomalously increase upon moving into the supercooled region according to a power law that would diverge at Ts{\approx}225 K,(1,2) so there may be a link between water's thermodynamic anomalies and the crystallization rate of ice. But probing this link is challenging because fast crystallization prevents experimental studies of the liquid below TH. And while atomistic studies have captured water crystallization(3), the computational costs involved have so far prevented an assessment of the rates and mechanism involved. Here we report coarse-grained molecular simulations with the mW water model(4) in the supercooled regime around TH, which reveal that a sharp increase in the fraction of four-coordinated molecules in supercooled liquid water explains its anomalous thermodynamics and also controls the rate and mechanism of ice formation. The simulations reveal that the crystallization rate of water reaches a maximum around 225 K, below which ice nuclei form faster than liquid water can equilibrate. This implies a lower limit of metastability of liquid water just below TH and well above its glass transition temperature Tg{\approx}136 K. By providing a relationship between the structural transformation in liquid water, its anomalous thermodynamics and its crystallization rate, this work provides a microscopic foundation to the experimental finding that the thermodynamics of water determines the rates of homogeneous nucleation of ice.(5)
△ Less
Submitted 27 September, 2011; v1 submitted 6 July, 2011;
originally announced July 2011.
-
Water modeled as an intermediate element between carbon and silicon
Authors:
Valeria Molinero,
Emily B. Moore
Abstract:
Water and silicon are chemically dissimilar substances with common physical properties. Their liquids display a temperature of maximum density, increased diffusivity on compression, they form tetrahedral crystals and tetrahedral amorphous phases. The common feature to water, silicon and carbon is the formation of tetrahedrally coordinated units. We exploit these similarities to develop a coarse-…
▽ More
Water and silicon are chemically dissimilar substances with common physical properties. Their liquids display a temperature of maximum density, increased diffusivity on compression, they form tetrahedral crystals and tetrahedral amorphous phases. The common feature to water, silicon and carbon is the formation of tetrahedrally coordinated units. We exploit these similarities to develop a coarse-grained model of water (mW) that is essentially an atom with tetrahedrality intermediate between carbon and silicon. mW mimics the hydrogen-bonded structure of water through the introduction of a nonbond angular dependent term that encourages tetrahedral configurations. The model departs from the prevailing paradigm in water modeling: the use of long-ranged forces (electrostatics) to produce short-ranged (hydrogen-bonded) structure. mW has only short-range interactions yet it reproduces the energetics, density and structure of liquid water, its anomalies and phase transitions with comparable or better accuracy than the most popular atomistic models of water, at less than 1% of the computational cost. We conclude that it is not the nature of the interactions but the connectivity of the molecules that determines the structural and thermodynamic behavior of water. The speedup in computing time provided by mW makes it particularly useful for the study of slow processes in deeply supercooled water, the mechanism of ice nucleation, wetting-drying transitions, and as a realistic water model for coarse-grained simulations of biomolecules and complex materials.
△ Less
Submitted 16 September, 2008;
originally announced September 2008.
-
Tuning of tetrahedrality in a silicon potential yields a series of monatomic (metal-like) glassformers of very high fragility
Authors:
Valeria Molinero,
Srikanth Sastry,
C. Austen Angell
Abstract:
We obtain monatomic glass formers in simulations by modifying the tetrahedral character in a silicon potential to explore a triple point zone between potentials favoring diamond (dc) and bcc crystals. dc crystallization is always preceded by a polyamorphic transformation of the liquid, and is frustrated when the Kauzmann temperature of the high temperature liquid intersects the liquid-liquid coe…
▽ More
We obtain monatomic glass formers in simulations by modifying the tetrahedral character in a silicon potential to explore a triple point zone between potentials favoring diamond (dc) and bcc crystals. dc crystallization is always preceded by a polyamorphic transformation of the liquid, and is frustrated when the Kauzmann temperature of the high temperature liquid intersects the liquid-liquid coexistence line. The glass forming liquids are extraordinarily fragile. Our results suggest that Si and Ge liquids may be vitrified at a pressure close to the diamond-beta-tin-liquid triple point.
△ Less
Submitted 15 August, 2006; v1 submitted 11 October, 2005;
originally announced October 2005.