Nature of Lone-Pair-Surface Bonds and Their Scaling Relations.

We investigate the (surface) bonding of a class of industrially and biologically important molecules in which the chemically active orbital is a 2 p electron lone pair located on an N or O atom bound via single bonds to H or alkyl groups. This class includes water, ammonia, alcohols, ethers, and amines. Using extensive density functional theory (DFT) calculations, we discover scaling relations (correlations) among molecular binding energies of different members of this class: the bonding energetics of a single member can be used as a descriptor for other members. We investigate the bonding mechanism for a representative (H2O) and find the most important physical surface properties that dictate the strength and nature of the bonding through a combination of covalent and noncovalent electrostatic effects. We describe the importance of surface intrinsic electrostatic, geometric, and mechanical properties in determining the extent of the lone-pair-surface interactions. We study systems including ionic materials in which the surface positive and negative centers create strong local surface electric fields, which polarize the dangling lone pair and lead to a strong "electrostatically driven bond". We emphasize the importance of noncovalent electrostatic effects and discuss why a fully covalent picture, common in the current first-principles literature on surface bonding of these molecules, is not adequate to correctly describe the bonding mechanism and energy trends. By pointing out a completely different mechanism (charge transfer) as the major factor for binding N- and O-containing unsaturated (radical) adsorbates, we explain why their binding energies can be tuned independently from those of the aforementioned species, having potential implications in scaling-driven catalyst discovery.

Mechanistic insights into heterogeneous methane activation.

While natural gas is an abundant chemical fuel, its low volumetric energy density has prompted a search for catalysts able to transform methane into more useful chemicals. This search has often been aided through the use of transition state (TS) scaling relationships, which estimate methane activation TS energies as a linear function of a more easily calculated descriptor, such as final state energy, thus avoiding tedious TS energy calculations. It has been shown that methane can be activated via a radical or surface-stabilized pathway, both of which possess a unique TS scaling relationship. Herein, we present a simple model to aid in the prediction of methane activation barriers on heterogeneous catalysts. Analogous to the universal radical TS scaling relationship introduced in a previous publication, we show that a universal TS scaling relationship that transcends catalysts classes also seems to exist for surface-stabilized methane activation if the relevant final state energy is used. We demonstrate that this scaling relationship holds for several reducible and irreducible oxides, promoted metals, and sulfides. By combining the universal scaling relationships for both radical and surface-stabilized methane activation pathways, we show that catalyst reactivity must be considered in addition to catalyst geometry to obtain an accurate estimation for the TS energy. This model can yield fast and accurate predictions of methane activation barriers on a wide range of catalysts, thus accelerating the discovery of more active catalysts for methane conversion.

Understanding trends in C-H bond activation in heterogeneous catalysis.

While the search for catalysts capable of directly converting methane to higher value commodity chemicals and liquid fuels has been active for over a century, a viable industrial process for selective methane activation has yet to be developed. Electronic structure calculations are playing an increasingly relevant role in this search, but large-scale materials screening efforts are hindered by computationally expensive transition state barrier calculations. The purpose of the present letter is twofold. First, we show that, for the wide range of catalysts that proceed via a radical intermediate, a unifying framework for predicting C-H activation barriers using a single universal descriptor can be established. Second, we combine this scaling approach with a thermodynamic analysis of active site formation to provide a map of methane activation rates. Our model successfully rationalizes the available empirical data and lays the foundation for future catalyst design strategies that transcend different catalyst classes.

Predicting Promoter-Induced Bond Activation on Solid Catalysts Using Elementary Bond Orders.

In this Letter, we examine bond activation induced by nonmetal surface promoters in the context of dehydrogenation reactions. We use C-H bond activation in methane dehydrogenation on transition metals as an example to understand the origin of the promoting or poisoning effect of nonmetals. The electronic structure of the surface and the bond order of the promoter are found to establish all trends in bond activation. On the basis of these results, we develop a predictive model that successfully describes the energetics of C-H, O-H, and N-H bond activation across a range of reactions. For a given reaction step, a single data point determines whether a nonmetal will promote bond activation or poison the surface and by how much. We show how our model leads to general insights that can be directly used to predict bond activation energetics on transition metal sulfides and oxides, which can be perceived as promoted surfaces. These results can then be directly used in studies on full catalytic pathways.