Understanding the role of frustrated Lewis pairs as ligands in transition metal-catalyzed reactions

Heterogeneous Frustrated Lewis Pair Catalysts: Rational Structure Design and Mechanistic Elucidation Based on Intrinsic Properties of Supports

ConspectusThe discovery of reversible hydrogenation using metal-free phosphoborate species in 2006 marked the official advent of frustrated Lewis pair (FLP) chemistry. This breakthrough revolutionized homogeneous catalysis approaches and paved the way for innovative catalytic strategies. The unique reactivity of FLPs is attributed to the Lewis base (LB) and Lewis acid (LA) sites either in spatial separation or in equilibrium, which actively react with molecules. Since 2010, heterogeneous FLP catalysts have gained increasing attention for their ability to enhance catalytic performance through tailored surface designs and improved recyclability, making them promising for industrial applications. Over the past 5 years, our group has focused on investigating and strategically modifying various types of solid catalysts with FLPs that are unique from classic solid FLPs. We have explored systematic characterization techniques to unravel the underlying mechanisms between the active sites and reactants. Additionally, we have demonstrated the critical role of catalysts' intrinsic electronic and geometric properties in promoting FLP formation and stimulating synergistic effects. The characterization of FLP catalysts has been greatly enhanced by the use of advanced techniques such as synchrotron X-ray diffraction, neutron powder diffraction, X-ray photoelectron spectroscopy, extended X-ray absorption fine structure, elemental mapping in scanning transmission electron microscopy, electron paramagnetic resonance spectroscopy, diffuse-reflectance infrared Fourier transform spectroscopy, and solid-state nuclear magnetic resonance spectroscopy. These techniques have provided deeper insights into the structural and electronic properties of FLP systems for the future design of catalysts.Understanding electron distribution in the overlapping orbitals of LA and LB pairs is essential for inducing FLPs in operando in heterogeneous catalysts through target electron reallocation by external stimuli. For instance, in silicoaluminophosphate-type zeolites with weak orbital overlap, the adsorption of polar gas molecules leads to heterolytic cleavage of the Alδ+-Oδ- bond, creating unquenched LA-LB pairs. In a Ru-doped metal-organic framework, the Ru-N bond can be polarized through metal-ligand charge transfer under light, forming Ru+-N- pairs. This activation of FLP sites from the framework represents a groundbreaking innovation that expands the catalytic potential of existing materials. For catalysts already employing FLP chemistry to dynamically generate products from substrates, a complete mechanistic interpretation requires a thorough examination of the surface electronic properties and the surrounding environment. The hydrogen spillover ability on the Ru-doped FLP surfaces improves conversion efficiency by suppressing hydrogen poisoning at metal sites. In situ H2-H2O conditions enable the production of organic chemicals with excellent activity and selectivity by creating new bifunctional sites via FLP chemistry. By highlighting the novel FLP systems featuring FLP induction and synergistic effects and the selection of advanced characterization techniques to elucidate reaction mechanisms, we hope that this Account will offer innovative strategies for designing and characterizing FLP chemistry in heterogeneous catalysts to the research community.

Insights into the Reactivity of the Ring-Opening Reaction of Tetrahydrofuran by Intramolecular Group-13/P- and Al/Group-15-Based Frustrated Lewis Pairs

A theoretical study concerning key factors affecting activation energies for ring-opening reactions of tetrahydrofuran (THF) by G13/P-based (G13 = B, Al, Ga, In, and Tl) and Al/G15-based (G15 = N, P, As, Sb, and Bi) frustrated Lewis pairs (FLPs) featuring the dimethylxanthene scaffold was performed using density functional theory. Our theoretical findings indicate that only dimethylxanthene backbone Al/P-Rea (Rea = reactant) FLP-type molecules can be energetically favorable to undergo the ring-opening reaction with THF. Our theoretical evidence reveals that the shorter the separating distance between Lewis acidic (LA) and Lewis basic (LB) centers of the dimethylxanthene backbone FLP-type molecules, the greater the orbital overlaps between the FLP and THF and the lower the activation barrier for such a ring-opening reaction. Energy decomposition analysis (EDA) evidence suggests that the bonding interaction for such a ring-opening reaction is predominated by the donor-acceptor interaction (singlet-singlet interaction) compared to the electron-sharing interaction (triplet-triplet interaction). In addition, the natural orbitals for chemical valence (NOCV) evidence demonstrate that the bonding situations of such ring-opening reactions can be best described as FLP-to-THF forward bonding (the lone pair (G15) → the empty σ*(C-O)) and THF-to-FLP back bonding (the empty σ*(G13) ← filled p-π(O)). The EDA-NOCV observations show that the former plays a predominant role and the latter plays a minor role in such bonding conditions. The activation strain model reveals that the deformation energy of THF is the key factor in determining the activation energy of their ring-opening reactions. Comparing the geometrical structures of the transition states with their corresponding reactants, a linear relationship between them can be rationally explained by the Hammond postulate combined with the respective activation barriers calculated in this work.

Bonding Situation of σ-E-H Complexes in Transition Metal and Main Group Compounds

The ambiguous bonding situation of σ-E-H (E=Si, B) complexes in transition metal compounds has been rationalized by means of Density Functional Theory calculations. To this end, the combination of the Energy Decomposition Analysis (EDA) method and its Natural Orbital for Chemical Valance (NOCV) extension has been applied to representative complexes described in the literature where the possible η¹ versus η² coordination mode is not unambiguously defined. Our quantitative analyses, which complement previous data based on the application of the Quantum Theory of Atoms in Molecules (QTAIM) approach, indicate that there exists a continuum between genuine η¹ and η² modes depending mainly on the strength of the backdonation. Finally, we also applied this EDA-NOCV approach to related main-group species where the backdonation is minimal.

Rationalizing the Substituent Effects in Diels-Alder Reactions of Triazolinediones with Anthracene

In this work we tackle the problem of the substituent effects in the Diels-Alder cycloadditions between triazolinediones (TADs) and anthracene. Experiments showed that aryl TADs substituted with electron-withdrawing groups (EWG) are more reactive than those substituted with electron-donating (EDG) or alkyl groups. However, the molecular origin of this preference is not yet understood. By a combination of methods including the activation strain model (ASM), energy decomposition analysis (EDA), molecular orbital (MO) theory, and conceptual density functional theory (CDFT), we disclosed the substituent effects of TADs. First, ASM/EDA analysis revealed that the reactivity of alkyl and aryl-substituted TADs is controlled by interaction energies, ΔE_int, which are ultimately defined by orbital interactions between frontier molecular orbitals. Moreover, alkyl-TADs are also controlled by the extent of strain at the transition state. The MO analysis suggested that the rate acceleration for EWG-substituted TADs is due to a more favorable orbital interaction between the HOMO of anthracene and the LUMO of the TADs, which is corroborated by calculations of charge transfer at the transition states. From CDFT, the chemical potential of anthracene is higher than those of TADs, indicating a flow of electron density from anthracene to TADs, in agreement with the results from the electrophilicity index.

A neutral, acyclic, borataalkene-like ligand for group 11 metals: L- and Z-type ligands side by side

The overall neutral α-borylated phosphorus ylide Ph₃PC(Me)BEt₂ behaves like a polar borataalkene and can act as acyclic, ambiphilic π-type ligand with L- and Z-type functionalities side by side. In the complexes [MX{η²-Ph₃PC(Me)BEt₂}] (M = Cu, (Ag), Au; X = Cl, NTf₂), the bonding is dominated by the highly nucleophilic ylidic carbon atom (L-type ligand). The Lewis acidic boron atom furnishes nonetheless a small but significant bonding contribution (Z-type ligand).

Rationalizing the influence of α-cationic phospholes on π-catalysis

The physical factors behind the experimentally observed high activity of gold(I)-catalysts having an α-cationic phosphole as a ligand have been computationally explored. To this end, the gold(I)-catalysed hydroarylation reactions of phenylacetylene and mesitylene involving both neutral and cationic phosphole as well as phosphine ligands have been quantitatively analyzed in detail with the help of the activation strain model of reactivity in combination with the energy decomposition analysis method. It is found that the cationic phosphole ligands induce a dramatic change in both the geometry and the electronic structure of the initially formed π-complex which significantly enhances its electrophilicity. This results in an enhancement of the key π(mesitylene) → π*(LAu-acetylene complex) molecular orbital interaction which is the main factor responsible for the activating effect of these cationic ligands.

Significant Insight into the Origin of Reaction Barriers Determining Dihydrogen Activation by G13-P-P (G13 = Group 13 Element) and G15-P-Ga (G15 = Group 15 Element) Frustrated Lewis Pairs

The heterolytic cleavage of H₂ by multiply bonded phosphorus-bridged G13-P-P-Rea (G13 = B, Al, Ga, In, and Tl) and G15-P-Ga-Rea (G15 = N, P, As, Sb, and Bi) frustrated Lewis pairs (FLPs) has been theoretically investigated using density functional theory calculations. For the above nine FLP-type molecules, our theoretical findings suggest that only Al-P-P-Rea, Ga-P-P-Rea, and In-P-P-Rea can undergo the energetically feasible H₂ activation reaction from kinetic and thermodynamic viewpoints. Our study based on the activation strain model (ASM) reveals that gaining a better orbital overlap between G13-P-P-Rea and G15-P-Ga-Rea molecules and H₂ affected the reaction barriers through the atomic radius of G13 and G15. According to our energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) results, the bonding of these H₂ activation reactions involving G13-P-P-Rea and G15-P-Ga-Rea is dominated by the donor-acceptor interaction (singlet-singlet interaction) rather than the electron-sharing interaction (triplet-triplet interaction). Moreover, our EDA-NOCV evidence reveals that the best description for the above bonding situations is the lone pair(G15) → σ*(H₂) interaction rather than the empty p-π-orbital(G13) ← σ(H₂) interaction. In particular, the findings in this work based on theoretically calculated geometries and the corresponding relative free energies of the stationary points combined with the results from the above sophisticated methods nicely agree with the famous Hammond postulate.

Chemical reactivity from an activation strain perspective

Chemical reactions are ubiquitous in the universe, they are at the core of life, and they are essential for industrial processes. The drive for a deep understanding of how something occurs, in this case, the mechanism of a chemical reaction and the factors controlling its reactivity, is intrinsically valuable and an innate quality of humans. The level of insight and degree of understanding afforded by computational chemistry cannot be understated. The activation strain model is one of the most powerful tools in our arsenal to obtain unparalleled insight into reactivity. The relative energy of interacting reactants is evaluated along a reaction energy profile and related to the rigidity of the reactants' molecular structure and the strength of the stabilizing interactions between the deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). Owing to the connectedness between the activation strain model and Kohn-Sham molecular orbital theory, one is able to obtain a causal relationship between both the sterics and electronics of the reactants and their mutual reactivity. Only when this is accomplished one can eclipse the phenomenological explanations that are commonplace in the literature and textbooks and begin to rationally tune and optimize chemical transformations. We showcase how the activation strain model is the ideal tool to elucidate fundamental organic reactions, the activation of small molecules by metallylenes, and the cycloaddition reactivity of cyclic diene- and dipolarophiles.