Δ ² machine learning for reaction property prediction

Reactive potentials serve as essential tools for investigating chemical reactions with moderate computational costs. However, traditional reactive potentials often depend on fixed, semiempirical parameters, which limits their accuracy and transferability. Overcoming these limitations can significantly expand the applicability of reactive potentials, enabling the simulation of a broader range of reactions under diverse conditions and the prediction of reaction properties, such as barrier heights. This work introduces ANI-1xBB, a novel ANI-based reactive ML potential trained on off-equilibrium molecular conformers generated through an automated bond-breaking workflow. ANI-1xBB significantly enhances the prediction of reaction energetics, barrier heights, and bond dissociation energies, surpassing those of conventional ANI models. Our results show that ANI-1xBB improves transition state modeling and reaction pathway prediction while generalizing effectively to pericyclic reactions and radical-driven processes. Furthermore, the automated data generation strategy supports the efficient construction of large-scale, high-quality reactive data sets, reducing reliance on expensive QM calculations. This work highlights ANI-1xBB as a practical model for accelerating the development of reactive machine learning potentials, offering new opportunities for modeling reaction phenomena.

We formulate a reaction prediction problem in terms of node-classification in a disconnected graph of source molecules and generalize a graph convolution neural network for disconnected graphs.

Machine learning interatomic potentials (MLIPs) have emerged as powerful tools for investigating atomistic systems with high accuracy and a relatively low computational cost. However, a common and unaddressed challenge with many current neural network (NN) MLIP models is their limited ability to accurately predict the relative energies of systems containing isolated or nearly isolated atoms, which appear in various reactive processes. To address this limitation, we present a mathematical technique for modifying any existing atom-centered NN architecture to account for the energies of isolated atoms. The result produces a consistent prediction of the atomization energy (AE) of a system using minimal constraints on the model. Using this technique, we build a model architecture that we call hierarchically interacting particle neural network (HIP-NN)-AE, an AE-constrained version of the HIP-NN, as well as ANI-AE, the AE-constrained version of the accurate NN engine for molecular energies (ANI). Our results demonstrate AE consistency of AE-constrained models, which drastically improves the AE predictions for the models. We compare the AE-constrained approach to unconstrained models as well as models from the literature in other scenarios, such as bond dissociation energies, bond dissociation pathways, and extensibility tests. These results show that the constraints improve the model performance in some of these tasks and do not negatively affect the performance on any tasks. The AE constraint approach thus offers a robust solution to the challenges posed by isolated atoms in energy prediction tasks.

Reaction mechanism search tools have demonstrated the ability to provide insights into likely products and rate-limiting steps of reacting systems. However, reactions involving several concerted bond changes─as can be found in many key steps of natural product syntheses─can complicate the search process. To mitigate these complications, we present a mechanism search strategy particularly suited to help expedite exploration of an exemplary family of such complex reactions, cyclizations. We provide a cost-effective strategy for identifying relevant elementary reaction steps by combining graph-based enumeration schemes and machine learning techniques for intermediate filtering. Key to this approach is our use of a neural network potential (NNP), AIMNet2-rxn, for computational evaluation of each candidate reaction pathway. In this article, we evaluate the NNP's ability to estimate activation energies, demonstrate the correct anticipation of stereoselectivity, and recapitulate complex enabling steps in natural product synthesis.