New Strategies for Direct Methane-to-Methanol Conversion from Active Learning Exploration of 16 Million Catalysts

Directional multiobjective optimization of metal complexes at the billion-system scale

The discovery of transition metal complexes (TMCs) with optimal properties requires large ligand libraries and efficient multiobjective optimization algorithms. Here we provide the tmQMg-L library, containing 30k diverse and synthesizable ligands with robustly assigned charges and metal coordination modes. tmQMg-L enabled the generation of 1.37 million palladium TMCs, which were used to develop and benchmark the Pareto-Lighthouse multiobjective genetic algorithm (PL-MOGA). With fine control over aim and scope, this algorithm maximized both the polarizability and highest occupied molecular orbital-lowest unoccupied molecular orbital gap of the TMCs within selected regions of the Pareto front, without requiring prior knowledge on the objective limits. Instead of genetic operations on small ligand fragments, the PL-MOGA did whole-ligand mutation and crossover operations, which in chemical spaces containing billions of systems, yielded thousands of highly diverse TMCs in an interpretable manner.

Automated de Novo Design of Olefin Metathesis Catalysts: Computational and Experimental Analysis of a Simple Thermodynamic Design Criterion

Methods for computational de novo design of inorganic molecules have paved the way for automated design of homogeneous catalysts. Such studies have so far relied on correlation-based prediction models as fitness functions (figures of merit), but the soundness of these approaches has yet to be tested by experimental verification of de novo-designed catalysts. Here, a previously developed criterion for the optimization of dative ligands L in ruthenium-based olefin metathesis catalysts RuCl2(L)(L')(═CHAr), where Ar is an aryl group and L' is a phosphine ligand dissociating to activate the catalyst, was used in de novo design experiments. These experiments predicted catalysts bearing an N-heterocyclic carbene (L = 9) substituted by two N-bound mesityls and two tert-butyl groups at the imidazolidin-2-ylidene backbone to be promising. Whereas the phosphine-stabilized precursor assumed by the prediction model could not be made, a pyridine-stabilized ruthenium alkylidene complex (17) bearing carbene 9 was less active than a known leading pyridine-stabilized Grubbs-type catalyst (18, L = H2IMes). A density functional theory-based analysis showed that the unsubstituted metallacyclobutane (MCB) intermediate generated in the presence of ethylene is the likely resting state of both 17 and 18. Whereas the design criterion via its correlation between the stability of the MCB and the rate-determining barrier indeed seeks to stabilize the MCB, it relies on RuCl2(L)(L')(═CH2) adducts as resting states. The change in resting state explains the discrepancy between the prediction and the actual performance of catalyst 17. To avoid such discrepancies and better address the multifaceted challenges of predicting catalytic performance, future de novo catalyst design studies should explore and test design criteria incorporating information from more than a single relative energy or intermediate.

Economically viable co-production of methanol and sulfuric acid via direct methane oxidation

The direct oxidation of methane to methanol has been spotlighted research for decades, but has never been commercialized. This study introduces cost-effective process for co-producing methanol and sulfuric acid through a direct oxidation of methane. In the initial phase, methane oxidation forms methyl bisulfate (CH3OSO3H), then transformed into methyl trifluoroacetate (CF3CO2CH3) via esterification, and hydrolyzed into methanol. This approach eliminates the need for energy-intensive separation of methyl bisulfate from sulfuric acid by replacing the former with methyl trifluoroacetate. Through the superstructure optimization, our sequential process reduces the levelized cost of methanol to nearly two-fold reduction from the current market price. Importantly, this process demonstrates adaptability to smaller gas fields, assuring its economical operation across a broad range of gas fields. The broader application of this process could substantially mitigate global warming by utilizing methane, leading to a significantly more sustainable and economically beneficial methanol industry.

Toward De Novo Catalyst Discovery: Fast Identification of New Catalyst Candidates for Alcohol-Mediated Morita-Baylis-Hillman Reactions

Recently we have demonstrated how a genetic algorithm (GA) starting from random tertiary amines can be used to discover a new and efficient catalyst for the alcohol-mediated Morita-Baylis-Hillman (MBH) reaction. In particular, the discovered catalyst was shown experimentally to be eight times more active than DABCO, commonly used to catalyze the MBH reaction. This represents a breakthrough in using generative models for catalyst optimization. However, the GA procedure, and hence discovery, relied on two important pieces of information; 1) the knowledge that tertiary amines catalyze the reaction and 2) the mechanism and reaction profile for the catalyzed reaction, in particular the transition state structure of the rate-determining step. Thus, truly de novo catalyst discovery must include these steps. Here we present such a method for discovering catalyst candidates for a specific reaction while simultaneously proposing a mechanism for the catalyzed reaction. We show that tertiary amines and phosphines are potential catalysts for the MBH reaction by screening 11 molecular templates representing common functional groups. The method relies on an automated reaction discovery workflow using meta-dynamics calculations. Combining this method for catalyst candidate discovery with our GA-based catalyst optimization method results in an algorithm for truly de novo catalyst discovery.

Theoretical Perspective of Promoting Direct Methane-to-Methanol Conversion at Complex Metal Oxide-Metal Interfaces

Direct methane conversion to methanol has been considered as an effective and economic way to address greenhouse effects and the current high demand for methanol in industry. However, the process has long been challenging due to lack of viable catalysts to compromise the activation of methane that typically occurs at high temperatures and retaining of produced methanol that requires mild conditions. This Perspective demonstrates an effective strategy to promote direct methane to methanol conversion by engineering the active sites and chemical environments at complex metal oxide - copper oxide - copper interfaces. Such effort strongly depends on extensive theoretical studies by combining density functional theory (DFT) calculations and kinetic Monte Carlo (KMC) simulations to provide in-depth understanding of reaction mechanism and active sites, which build a strong basis to enable the identification of design principles and advance the catalyst optimization for selective CH₄-to-CH₃OH conversion.

Computational Discovery of Stable Metal-Organic Frameworks for Methane-to-Methanol Catalysis

The challenge of direct partial oxidation of methane to methanol has motivated the targeted search of metal-organic frameworks (MOFs) as a promising class of materials for this transformation because of their site-isolated metals with tunable ligand environments. Thousands of MOFs have been synthesized, yet relatively few have been screened for their promise in methane conversion. We developed a high-throughput virtual screening workflow that identifies MOFs from a diverse space of experimental MOFs that have not been studied for catalysis, yet are thermally stable, synthesizable, and have promising unsaturated metal sites for C-H activation via a terminal metal-oxo species. We carried out density functional theory calculations of the radical rebound mechanism for methane-to-methanol conversion on models of the secondary building units (SBUs) from 87 selected MOFs. While we showed that oxo formation favorability decreases with increasing 3d filling, consistent with prior work, previously observed scaling relations between oxo formation and hydrogen atom transfer (HAT) are disrupted by the greater diversity in our MOF set. Accordingly, we focused on Mn MOFs, which favor oxo intermediates without disfavoring HAT or leading to high methanol release energies─a key feature for methane hydroxylation activity. We identified three Mn MOFs comprising unsaturated Mn centers bound to weak-field carboxylate ligands in planar or bent geometries with promising methane-to-methanol kinetics and thermodynamics. The energetic spans of these MOFs are indicative of promising turnover frequencies for methane to methanol that warrant further experimental catalytic studies.

Cyclic iron tetra N-heterocyclic carbenes: synthesis, properties, reactivity, and catalysis

Cyclic iron tetracarbenes are an emerging class of macrocyclic iron N-heterocyclic carbene (NHC) complexes. They can be considered as an organometallic compound class inspired by their heme analogs, however, their electronic properties differ, e.g. due to the very strong σ-donation of the four combined NHCs in equatorial coordination. The ligand framework of iron tetracarbenes can be readily modified, allowing fine-tuning of the structural and electronic properties of the complexes. The properties of iron tetracarbene complexes are discussed quantitatively and correlations are established. The electronic nature of the tetracarbene ligand allows the isolation of uncommon iron(III) and iron(IV) species and reveals a unique reactivity. Iron tetracarbenes are successfully applied in C-H activation, CO₂ reduction, aziridination and epoxidation catalysis and mechanisms as well as decomposition pathways are described. This review will help researchers evaluate the structural and electronic properties of their complexes and target their catalyst properties through ligand design.

Active Learning Exploration of Transition-Metal Complexes to Discover Method-Insensitive and Synthetically Accessible Chromophores

Transition-metal chromophores with earth-abundant transition metals are an important design target for their applications in lighting and nontoxic bioimaging, but their design is challenged by the scarcity of complexes that simultaneously have well-defined ground states and optimal target absorption energies in the visible region. Machine learning (ML) accelerated discovery could overcome such challenges by enabling the screening of a larger space but is limited by the fidelity of the data used in ML model training, which is typically from a single approximate density functional. To address this limitation, we search for consensus in predictions among 23 density functional approximations across multiple rungs of "Jacob's ladder". To accelerate the discovery of complexes with absorption energies in the visible region while minimizing the effect of low-lying excited states, we use two-dimensional (2D)efficient global optimization to sample candidate low-spin chromophores from multimillion complex spaces. Despite the scarcity (i.e., ∼0.01%) of potential chromophores in this large chemical space, we identify candidates with high likelihood (i.e., >10%) of computational validation as the ML models improve during active learning, representing a 1000-fold acceleration in discovery. Absorption spectra of promising chromophores from time-dependent density functional theory verify that 2/3 of candidates have the desired excited-state properties. The observation that constituent ligands from our leads have demonstrated interesting optical properties in the literature exemplifies the effectiveness of our construction of a realistic design space and active learning approach.

Transition metal oxide complexes as molecular catalysts for selective methane to methanol transformation: any prospects or time to retire?

Transition metal oxides have been extensively used in the literature for the conversion of methane to methanol. Despite the progress made over the past decades, no method with satisfactory performance or economic viability has been detected. The main bottleneck is that the produced methanol oxidizes further due to its weaker C-H bond than that of methane. Every improvement in the efficiency of a catalyst to activate methane leads to reduction of the selectivity towards methanol. Is it therefore prudent to keep studying (both theoretically and experimentally) metal oxides as catalysts for the quantitative conversion of methane to methanol? This perspective focuses on molecular metal oxide complexes and suggests strategies to bypass the current bottlenecks with higher weight on the computational chemistry side. We first discuss the electronic structure of metal oxides, followed by assessing the role of the ligands in the reactivity of the catalysts. For better selectivity, we propose that metal oxide anionic complexes should be explored further, while hydrophylic cavities in the vicinity of the metal oxide can perturb the transition-state structure for methanol increasing appreciably the activation barrier for methanol. We also emphasize that computational studies should target the activation reaction of methanol (and not only methane), the study of complete catalytic cycles (including the recombination and oxidation steps), and the use of molecular oxygen as an oxidant. The titled chemical conversion is an excellent challenge for theory and we believe that computational studies should lead the field in the future. It is finally shown that bottom-up approaches offer a systematic way for exploration of the chemical space and should still be applied in parallel with the recently popular machine learning techniques. To answer the question of the title, we believe that metal oxides should still be considered provided that we change our focus and perform more systematic investigations on the activation of methanol.

Molecular Machine Learning for Chemical Catalysis: Prospects and Challenges

ConspectusIn the domain of reaction development, one aims to obtain higher efficacies as measured in terms of yield and/or selectivities. During the empirical cycles, an admixture of outcomes from low to high yields/selectivities is expected. While it is not easy to identify all of the factors that might impact the reaction efficiency, complex and nonlinear dependence on the nature of reactants, catalysts, solvents, etc. is quite likely. Developmental stages of newer reactions would typically offer a few hundreds of samples with variations in participating molecules and/or reaction conditions. These "observations" and their "output" can be harnessed as valuable labeled data for developing molecular machine learning (ML) models. Once a robust ML model is built for a specific reaction under development, it can predict the reaction outcome for any new choice of substrates/catalyst in a few seconds/minutes and thus can expedite the identification of promising candidates for experimental validation. Recent years have witnessed impressive applications of ML in the molecular world, most of them aimed at predicting important chemical or biological properties. We believe that an integration of effective ML workflows can be made richly beneficial to reaction discovery.As with any new technology, direct adaptation of ML as used in well-developed domains, such as natural language processing (NLP) and image recognition, is unlikely to succeed in reaction discovery. Some of the challenges stem from ineffective featurization of the molecular space, unavailability of quality data and its distribution, in making the right choice of ML model and its technically robust deployment. It shall be noted that there is no universal ML model suitable for an inherently high-dimensional problem such as chemical reactions. Given these backgrounds, rendering ML tools conducive for reactions is an exciting as well as challenging endeavor at the same time. With the increased availability of efficient ML algorithms, we focused on tapping their potential for small-data reaction discovery (a few hundreds to thousands of samples).In this Account, we describe both feature engineering and feature learning approaches for molecular ML as applied to diverse reactions of high contemporary interest. Among these, catalytic asymmetric hydrogenation of imines/alkenes, β-C(sp³)-H bond functionalization, and relay Heck reaction employed a feature engineering approach using the quantum-chemically derived physical organic descriptors as the molecular features─all designed to predict the enantioselectivity. The selection of molecular features to customize it for a reaction of interest is described, along with emphasizing the chemical insights that could be gathered through the use of such features. Feature learning methods for predicting the yield of Buchwald-Hartwig cross-coupling, deoxyfluorination of alcohols, and enantioselectivity of N,S-acetal formation are found to offer excellent predictions. We propose a transfer learning protocol, wherein an ML model such as a language model is trained on a large number of molecules (10⁵-10⁶) and fine-tuned on a focused library of target task reactions, as an effective alternative for small-data reaction discovery (10²-10³ reactions). The exploitation of deep neural network latent space as a method for generative tasks to identify useful substrates for a reaction is demonstrated as a promising strategy.

Machine Learning Models Predict Calculation Outcomes with the Transferability Necessary for Computational Catalysis

Virtual high-throughput screening (VHTS) and machine learning (ML) have greatly accelerated the design of single-site transition-metal catalysts. VHTS of catalysts, however, is often accompanied with a high calculation failure rate and wasted computational resources due to the difficulty of simultaneously converging all mechanistically relevant reactive intermediates to expected geometries and electronic states. We demonstrate a dynamic classifier approach, i.e., a convolutional neural network that monitors geometry optimizations on the fly, and exploit its good performance and transferability in identifying geometry optimization failures for catalyst design. We show that the dynamic classifier performs well on all reactive intermediates in the representative catalytic cycle of the radical rebound mechanism for the conversion of methane to methanol despite being trained on only one reactive intermediate. The dynamic classifier also generalizes to chemically distinct intermediates and metal centers absent from the training data without loss of accuracy or model confidence. We rationalize this superior model transferability as arising from the use of electronic structure and geometric information generated on-the-fly from density functional theory calculations and the convolutional layer in the dynamic classifier. When used in combination with uncertainty quantification, the dynamic classifier saves more than half of the computational resources that would have been wasted on unsuccessful calculations for all reactive intermediates being considered.