Search

We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems. [ABSTRACT FROM AUTHOR]

The transport of small penetrants through disordered materials with glassy dynamics is encountered in applications ranging from drug delivery to chemical separations. Nonetheless, understanding the influence of the matrix structure and fluctuations on penetrant motions remains a persistent challenge. Here, we use event-driven molecular dynamics to investigate the transport of small, hard-sphere tracers embedded in matrices of square-well particles. Short-range attractions between matrix particles give rise to reentrant dynamics in the supercooled regime, in which the liquid's relaxation time increases dramatically upon heating or cooling. Heating results in a "repulsive" supercooled liquid where relaxations are frustrated by steric interactions between particles, whereas cooling produces an "attractive" liquid in which relaxations are hindered by long-lived interparticle bonds. Further cooling or heating, or compression, of the supercooled liquids results in the formation of distinct glasses. Our study reveals that tracer transport in these supercooled liquids and glasses is influenced by the matrix structure and dynamics. The relative importance of each factor varies between matrices and is examined in detail by analyzing particle mean-square displacements, caging behavior, and trajectories sampled from the isoconfigurational ensemble. We identify features of tracer dynamics that reveal the spatial and temporal heterogeneity of the matrices and show that matrix arrest is insufficient to localize tracers. [ABSTRACT FROM AUTHOR]

Two-phase simulations are commonly used to evaluate coexistence conditions, interfacial tensions, and other thermodynamic properties associated with first-order phase transitions. Calculation of these properties is often simplified when the interfaces between the two phases are flat or planar. Here, we derive a general thermodynamic criterion for selecting simulation cell dimensions to stabilize planar interfaces in phase-separated fluid-fluid systems with respect to homogeneous, single-phase states. The resulting expression is validated by analyzing the effects of simulation cell dimensions on the formation of planar liquid-vapor interfaces in the Lennard-Jones fluid and in the TIP4P/2005 model of water. We also perform large scale molecular dynamics simulations to study metastable liquid-liquid phase separation in the ST2 and TIP4P/2005 models of water under deeply supercooled conditions. Our simulations confirm the stability of a liquid-liquid interface in ST2, and they demonstrate that the corresponding interface for TIP4P/2005 can be stabilized by judiciously choosing the simulation cell aspect ratio in a manner consistent with the thermodynamic criterion. We posit that this sensitivity to the simulation cell aspect ratio may explain discrepancies between previous studies examining liquid-liquid separation in models of supercooled water. [ABSTRACT FROM AUTHOR]

Liquid water's isothermal compressibility (1) and isobaric heat capacity (2), and the magnitude of its thermal expansion coefficient (3), increase sharply on cooling below the equilibrium freezing point. Many experimental [...]

Recent equation of state calculations [E. Lascaris, Phys. Rev. Lett. 116, 125701 (2016)] for an ionic model of silica suggest that it undergoes a density-driven, liquid--liquid phase transition (LLPT) similar to the controversial transition hypothesized to exist in deeply supercooled water. Here, we perform extensive free energy calculations to scrutinize the model's low-temperature phase behavior and confirm the existence of a first-order phase transition between two liquids with identical compositions but different densities. The low-density liquid (LDL) exhibits tetrahedral order, which is partially disrupted in the high-density liquid (HDL) by the intrusion of additional particles into the primary neighbor shell. Histogram reweighting methods are applied to locate conditions of HDL-- LDL coexistence and the liquid spinodals that bound the two-phase region. Spontaneous liquid--liquid phase separation is also observed directly in large-scale molecular dynamics simulations performed inside the predicted two-phase region. Given its clear LLPT, we anticipate that this model may serve as a paradigm for understanding whether similar transitions occur in water and other tetrahedral liquids. [ABSTRACT FROM AUTHOR]

Recent equation of state calculations [E. Lascaris, Phys. Rev. Lett. 116, 125701 (2016)] for an ionic model of silica suggest that it undergoes a density-driven, liquid-liquid phase transition (LLPT) similar to the controversial transition hypothesized to exist in deeply supercooled water. Here, we perform extensive free energy calculations to scrutinize the model's low-temperature phase behavior and confirm the existence of a first-order phase transition between two liquids with identical compositions but different densities. The low-density liquid (LDL) exhibits tetrahedral order, which is partially disrupted in the high-density liquid (HDL) by the intrusion of additional particles into the primary neighbor shell. Histogram reweighting methods are applied to locate conditions of HDL- LDL coexistence and the liquid spinodals that bound the two-phase region. Spontaneous liquid-liquid phase separation is also observed directly in large-scale molecular dynamics simulations performed inside the predicted two-phase region. Given its clear LLPT, we anticipate that this model may serve as a paradigm for understanding whether similar transitions occur in water and other tetrahedral liquids. [ABSTRACT FROM AUTHOR]