Abstract:Finite-temperature lattice free energy differences between polymorphs of molecular crystals are fundamental to understanding and predicting the relative stability relationships underpinning polymorphism, yet are computationally expensive to obtain. Here, we implement and critically assess machine-learning-enabled targeted free energy calculations derived from flow-based generative models to compute the free energy difference between two ice crystal polymorphs (Ice XI and Ic), modelled with a fully flexible empirical classical force field. We demonstrate that even when remapping through an analytical reference distribution, such methods enable a cost-effective and accurate calculation of free energy differences between disconnected metastable ensembles when trained on locally ergodic data sampled exclusively from the ensembles of interest. Unlike classical free energy perturbation methods, such as the Einstein crystal method, the targeted approach analysed in this work requires no additional sampling of intermediate perturbed Hamiltonians, offering significant computational savings in the system sizes compared in this work. To systematically assess the accuracy of the method, we monitored the convergence of free energy estimates during training by implementing an overfitting-aware weighted averaging strategy. By comparing our results with ground-truth free energy differences computed with the Einstein crystal method, we assess the accuracy and efficiency of two different model architectures, employing two different representations of the supercells degrees of freedom (Cartesian vs. quaternion-based). We conduct our assessment by comparing free energy differences between crystal supercells of different sizes and temperatures and assessing the accuracy in extrapolating lattice free energies to the thermodynamic limit. While at low temperatures and in small system sizes, the models perform with similar accuracy, we note that for larger systems and high temperatures, the choice of representation is key to obtaining generalisable results of quality comparable to that obtained from the Einstein crystal method. We hope the current work to be a useful stepping stone towards efficient free energy calculations in larger, more conformationally flexible systems.

From a Liquid to a Crystal Without Going Through a First-Order Phase Transition: Determining the Free Energy of Melting with Glassy Intermediates

Estimating the Energy State of Liquids from Crystal Melting

An Efficient Method for Computing Excess Free Energy of Liquid

Computation of the chemical potential and solubility of amorphous solids

Tailoring Tetrahedral and Pair-Correlation Entropies of Glass-Forming Liquids for Energy Storage Applications at Ultralow Temperatures

Free Energy Of Crystal-Liquid Interface

Ice Phases under Ambient and High Pressure: Insights from Density Functional Theory

Estimation of the ice melting point in molecular dynamics simulations based on the finite-size effects

Theory of amorphous ices

Crystal melting influenced by particle cooperativity of the liquid

Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models

Accurate Lattice Free Energies of Packing Polymorphs from Probabilistic Generative Models

Time-dependent Density Functional Theory for the Freezing/melting Transition in Interfacial Systems

Phase Transition Between Crystalline Variants of Ordinary Ice.

Potential Energy Landscape of the Apparent First-Order Phase Transition between Low-Density and High-Density Amorphous Ice

Predicting Glass-to-Glass and Liquid-to-Liquid Phase Transitions in Water using Classical Nucleation Temperature

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

Energy Gap and the Ideal Glass as a Defective Crystal: A Lattice Model of Monatomic Systems

Thermodynamic and Dynamical Heterogeneities During Glass Transition of Water

Dielectric Anomalies in Crystalline Ice: Indirect Evidence of the Existence of a Liquid-Liquid Critical Point in H2O

Homogeneous nucleation of phase transformations in supercooled water