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.

Predicting Crystal Morphology Using a Geometric Descriptor: A Comparative Study of Elemental Crystals with High-Throughput DFT Calculations

Predicting Intermetallic Surface Energies with High-Throughput DFT and Convolutional Neural Networks

Two descriptors avoiding explicit bond consideration on relative surface energies of various terminations in a crystal

New Generalized Crystallographic Descriptors for Structural Machine Learning

Interfacial Structure Analysis for the Prediction of Morphology of Crystals and Implications for the Design of Tailor-Made Additives

A Theoretical Method to Calculate the Surface Free Energies of Crystals

Predicting surface-energy anisotropy of metals with geometric properties of surfaces and atoms

Application of Electron Back Scattering Diffraction in Facet Crystalline Orientation Analysis

A computational strategy in consideration of the crystal surface structures for the crystal morphology prediction of energetic material

Quantitative Prediction of Surface Energy of High-Entropy-alloys Based on Intrinsic Descriptors

Computer simulation of surface and adatom properties of Lennard-Jones solids: A comparison between face-centered-cubic and hexagonal-close-packed structures

Accurate Lattice Free Energies of Packing Polymorphs from Probabilistic Generative Models

Morphology Prediction Methods and Their Applications in Energetic Crystals

Modelling of Surface Energies of Elemental Crystals

Anharmonic Terms of the Potential Energy Surface: A Group Theoretical Approach

Composition based crystal materials symmetry prediction using machine learning with enhanced descriptors

Noncovalent Interactions and Crystal Structure Prediction of Energetic Materials

Discovery of Stable Surfaces with Extreme Work Functions by High‐Throughput Density Functional Theory and Machine Learning

Physics guided deep learning for generative design of crystal materials with symmetry constraints

Numerical Method for Calculating Nanocrystals’ Edge Energies from Experimentally Observed Shape Evolution

Inverse Wulff Construction for Surface Energies of Coexisting and Missing Surfaces of Crystal Particles