We use molecular dynamics simulations to compare and contrast the liquid-state anomalies in the Stillinger-Weber models of monatomic water (mW), silicon (Si), and germanium (Ge) over a fairly wide range of temperatures and densities. The relationships between structure, entropy, and mobility, as well as the extent of the regions of anomalous behavior, are discussed as a function of the degree of tetrahedrality. We map out the cascade of density, structural, pair entropy, excess entropy, viscosity, and diffusivity anomalies for these three liquids. Among the three liquids studied here, only mW displays anomalies in the thermal conductivity, and this anomaly is evident only at very low temperatures. Diffusivity and viscosity, on the other hand, show pronounced anomalous regions for the three liquids. The temperature of maximum density of the three liquids shows re-entrant behavior consistent with either singularity-free or liquid-liquid critical point scenarios proposed to explain thermodynamic anomalies. The order-map, which shows the evolution of translational versus tetrahedral order in liquids, is different for Ge than for Si and mW. We find that although the monatomic water reproduces several thermodynamic and dynamic properties of rigid-body water models (e.g., SPC/E, TIP4P/2005), its sequence of anomalies follows, the same as Si and Ge, the silica-like hierarchy: the region of dynamic (diffusivity and viscosity) anomalies encloses the region of structural anomalies, which in turn encloses the region of density anomaly. The hierarchy of the anomalies based on excess entropy and Rosenfeld scaling, on the other hand, reverses the order of the structural and dynamic anomalies, i.e., predicts that the three Stillinger-Weber liquids follow a water-like hierarchy of anomalies. We investigate the scaling of diffusivity, viscosity, and thermal conductivity with the excess entropy of the liquid and find that for dynamical properties that present anomalies there is no universal scaling of the reduced property with excess entropy for the whole range of temperatures and densities. Instead, Rosenfeld's scaling holds for all the three liquids at high densities and high temperatures, although deviations from simple exponential dependence are observed for diffusivity and viscosity at lower temperatures and intermediate densities. The slope of the scaling of transport properties obtained for Ge is comparable to that obtained for simple liquids, suggesting that this low tetrahedrality liquid, although it stabilizes a diamond crystal, is already close to simple liquid behavior for certain properties.

The total, triplet, and pair contributions to the entropy with increasing tetrahedrality are mapped out for the Stillinger-Weber liquids to demonstrate the qualitative and quantitative differences between triplet-dominated, tetrahedral liquids and pair-dominated, simple liquids with regard to supercooling and crystallization. The heat capacity anomaly of tetrahedral liquids originates in local ordering due to both pair and triplet correlations. The results suggest that structural correlations can be directly related to thermodynamic anomalies, phase changes, and self-assembly in other atomic and colloidal fluids.

Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice− liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice−liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice−water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.

In the present molecular dynamics study, we investigate the effects of increasing pressure on the structural morphology of trihexyl(tetradecyl)phosphonium bromide (P/Br) and trihexyl(tetradecyl)phosphonium dicyanamide (P/DCA) ionic liquids (ILs). Special attention was paid to how charge and polarity orderings, which are present in the microscopic structure of these ILs at ambient conditions, respond to very high external pressure. The simulated X-ray scattering structure functions, S(q)s, of the two systems reveal that both the characteristic orderings show appreciable responsiveness towards the applied pressure change. At a given pressure, a slight difference between the polarity ordering (PO), charge ordering (CO), and adjacency correlations (AC) for both the systems points towards different microscopic structure of the two ILs due to change in anion. Beyond a certain pressure, we observe emergence of a new low-q peak in the S(q)s of both the systems. The new peak is associated with formation of crystalline order in these systems at higher pressures and the real space length-scale corresponding to the crystalline order lies in between those of polarity- and charge-ordering. Beyond the transition pressure, the crystallinity of both the systems increases with increasing pressure and the corresponding length-scale shifts towards smaller values upon increasing pressure. We also observe that the extent of the usual polarity ordering decreases upon increasing pressure for both the P/Br and P/DCA systems. We demonstrate that the disappearance of the usual polarity peak is due to decreased polar-polar and apolar-apolar correlations and enhanced correlations between the charged and uncharged groups of the ions. This scenario is completely reversed for the components corresponding to the crystalline order, the polar-polar and apolar-apolar correlations are enhanced and polar-apolar correlations are diminished at higher pressure. In addition, the charge ordering peak, which is not so obvious from the total S(q) but from ionic and sub-ionic partial components of it, shifts towards lower q values for P/Br. Instead, for the P/DCA, at the highest pressure studied the CO peak occurs at a q-value higher than that at the ambient pressure.

Porous crystalsincluding zeolites, metal–organic frameworks, and clathratesare widely used as catalysts and molecular sieves and in energy applications. These materials are synthesized from solution through an intermediate amorphous phase. However, the structural gap between the amorphous and porous crystal phases can lead to large nucleation barriers and difficulties in the control of crystal polymorphs. Previous reports indicate that porous mesophases can facilitate the amorphous to porous crystal transformation, directing the crystallization toward specific crystal structures. To date, it is not known how mesophase stability and synthesis temperature impact the facilitation. Here, we use molecular simulations and nucleation theory to investigate the crystallization of a zeolite in a family of models with tunable mesophase stability. We find that the nucleation mechanism evolves from non-classical to classical as the mesophase stability approaches that of the crystal. The simulations reveal that even an unstable mesophase can facilitate porous crystal nucleation through the formation of a transient fluctuation within which the crystal is originated. We conclude that tapping into the medium-range order of mesophases that have ordered pores without crystalline tiling is promising to increase the crystallization rates of porous crystals while directing the synthesis toward specific polymorphs.