Deconvolution and Analysis of the 1H NMR Spectra of Crude Reaction Mixtures.

Nuclear magnetic resonance (NMR) spectroscopy is an important analytical technique in synthetic organic chemistry, but its integration into high-throughput experimentation workflows has been limited by the necessity of manually analyzing the NMR spectra of new chemical entities. Current efforts to automate the analysis of NMR spectra rely on comparisons to databases of reported spectra for known compounds and, therefore, are incompatible with the exploration of new chemical space. By reframing the NMR spectrum of a reaction mixture as a joint probability distribution, we have used Hamiltonian Monte Carlo Markov Chain and density functional theory to fit the predicted NMR spectra to those of crude reaction mixtures. This approach enables the deconvolution and analysis of the spectra of mixtures of compounds without relying on reported spectra. The utility of our approach to analyze crude reaction mixtures is demonstrated with the experimental spectra of reactions that generate a mixture of isomers, such as Wittig olefination and C-H functionalization reactions. The correct identification of compounds in a reaction mixture and their relative concentrations is achieved with a mean absolute error as low as 1%.

CoeffNet: predicting activation barriers through a chemically-interpretable, equivariant and physically constrained graph neural network.

Activation barriers of elementary reactions are essential to predict molecular reaction mechanisms and kinetics. However, computing these energy barriers by identifying transition states with electronic structure methods (e.g., density functional theory) can be time-consuming and computationally expensive. In this work, we introduce CoeffNet, an equivariant graph neural network that predicts activation barriers using coefficients of any frontier molecular orbital (such as the highest occupied molecular orbital) of reactant and product complexes as graph node features. We show that using coefficients as features offer several advantages, such as chemical interpretability and physical constraints on the network's behaviour and numerical range. Model outputs are either activation barriers or coefficients of the chosen molecular orbital of the transition state; the latter quantity allows us to interpret the results of the neural network through chemical intuition. We test CoeffNet on a dataset of SN2 reactions as a proof-of-concept and show that the activation barriers are predicted with a mean absolute error of less than 0.025 eV. The highest occupied molecular orbital of the transition state is visualized and the distribution of the orbital densities of the transition states is described for a few prototype SN2 reactions.

Assessing the Accuracy of Density Functional Approximations for Predicting Hydrolysis Reaction Kinetics.

Hydrolysis reactions are ubiquitous in biological, environmental, and industrial chemistry. Density functional theory (DFT) is commonly employed to study the kinetics and reaction mechanisms of hydrolysis processes. Here, we present a new data set, Barrier Heights for HydrOlysis - 36 (BH2O-36), to enable the design of density functional approximations (DFAs) and the rational selection of DFAs for applications in aqueous chemistry. BH2O-36 consists of 36 diverse organic and inorganic forward and reverse hydrolysis reactions with reference energy barriers ΔE‡ calculated at the CCSD(T)/CBS level. Using BH2O-36, we evaluate 63 DFAs. In terms of mean absolute error (MAE) and mean relative absolute error (MRAE), ωB97M-V is the best-performing DFA tested, while MN12-L-D3(BJ) is the best-performing pure (nonhybrid) DFA. Broadly, we find that range-separated hybrid DFAs are necessary to approach chemical accuracy (0.043 eV). Although the best-performing DFAs include a dispersion correction to account for long-range interactions, we find that dispersion corrections do not generally improve MAE or MRAE for this data set.

Machine Learning Full NMR Chemical Shift Tensors of Silicon Oxides with Equivariant Graph Neural Networks.

The nuclear magnetic resonance (NMR) chemical shift tensor is a highly sensitive probe of the electronic structure of an atom and furthermore its local structure. Recently, machine learning has been applied to NMR in the prediction of isotropic chemical shifts from a structure. Current machine learning models, however, often ignore the full chemical shift tensor for the easier-to-predict isotropic chemical shift, effectively ignoring a multitude of structural information available in the NMR chemical shift tensor. Here we use an equivariant graph neural network (GNN) to predict full 29Si chemical shift tensors in silicate materials. The equivariant GNN model predicts full tensors to a mean absolute error of 1.05 ppm and is able to accurately determine the magnitude, anisotropy, and tensor orientation in a diverse set of silicon oxide local structures. When compared with other models, the equivariant GNN model outperforms the state-of-the-art machine learning models by 53%. The equivariant GNN model also outperforms historic analytical models by 57% for isotropic chemical shift and 91% for anisotropy. The software is available as a simple-to-use open-source repository, allowing similar models to be created and trained with ease.