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.

Accelerated structure-stability energy-free calculator

Accelerating Metadynamics-Based Free-Energy Calculations with Adaptive Machine Learning Potentials.

A generalized model for estimating adsorption energies of single atoms on doped carbon materialsâ€

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

Accelerated prediction of atomically precise cluster structures using on-the-fly machine learning

Prediction and Classification of Formation Energies of Binary Compounds by Machine Learning: An Approach without Crystal Structure Information

Accurate Free Energy Calculation via Multiscale Simulations Driven by Hybrid Machine Learning and Molecular Mechanics Potentials

Accelerating high-throughput searches for new alloys with active learning of interatomic potentials

Machine-Learned Fragment-Based Energies for Crystal Structure Prediction

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred Via Statistical Learning and Genetic Algorithm-Based Optimization

A critical examination of compound stability predictions from machine-learned formation energies

Atomic Energies from a Convolutional Neural Network.

Crystal Structure Prediction of Binary Alloys via Deep Potential

Accelerating the discovery of low-energy structure configurations: a computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Machine learning model for non-equilibrium structures and energies of simple molecules

Artificial intelligence real-time prediction and physical interpretation of atomic binding energies in nano-scale metal clusters

Accelerating crystal structure prediction by machine-learning interatomic potentials with active learning

FINETUNA: Fine-tuning Accelerated Molecular Simulations

Accurate Lattice Free Energies of Packing Polymorphs from Probabilistic Generative Models

Neural Network Potential with Multi-Resolution Approach Enables Accurate Prediction of Reaction Free Energies in Solution

Data-efficient modeling of catalytic reactions via enhanced sampling and on-the-fly learning of machine learning potentials