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.

A computationally efficient quasi-harmonic study of ice polymorphs using the FFLUX force field

Application of the FFLUX Force Field to Molecular Crystals: A Study of Formamide

DMC-ICE13: ambient and high pressure polymorphs of ice from Diffusion Monte Carlo and Density Functional Theory

Accurate Lattice Free Energies of Packing Polymorphs from Probabilistic Generative Models

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

Force field refinement for reproducing experimental infrared spectra of ionic liquids

Toward Gaussian Process Regression Modeling of a Urea Force Field

On-the-fly machine learning force field generation: Application to melting points

Quantum melting of charge ice and non-Fermi-liquid behavior: An exact solution for the extended Falicov-Kimball model in the ice-rule limit

High-accuracy thermodynamic properties to the melting point from ab initio calculations aided by machine-learning potentials

Data-efficient fine-tuning of foundational models for first-principles quality sublimation enthalpies

Stability-Ranking of Crystalline Ice Polymorphs Using Density-Functional Theory

Hydrogen mean force and anharmonicity in polycrystalline and amorphous ice

A Scalable Molecular Force Field Parameterization Method Based on Density Functional Theory and Quantum-Level Machine Learning

Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

Molecular modelling of ionic liquids: Perfluorinated anionic species with enlarged halogen substitutions

Machine learned Force-Fields for an ab-initio Quality Description of Metal-Organic Frameworks

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional