Activation of Small Molecules by Transition Metal Complexes via Computational Methods
Authors:   , Najafian Ahmad 
Published: February 2025

The first study project is based on modeling Earth abundant 3d transition-metal methoxide complexes with potentially redox-noninnocent ligands for methane C–H bond activation to form methanol (LnM-OMe + CH4 → LnM–Me + CH3OH). Three types of complex consisting of tridentate pincer terpyridine-like ligands, and different first-row transition metals (M = Ti, V, Cr, Mn, Fe, Co, Ni, and Cu) were modeled to elucidate the reaction mechanism as well as the effect of the metal identity on the thermodynamics and kinetics of a methane activation reaction. The calculations showed that the d electron count of the metal is a more significant factor than the metal's formal charge in controlling the thermodynamics and kinetics of C–H activation. These researches suggest that late 3d-metal methoxide complexes that favor σ-bond metathesis pathways for methane activation will yield lower barriers for C–H activation, and are more profitable catalyst for future studies. Second, subsequently, on the basis of the first project, density functional theory is used to analyze methane C−H activation by neutral and cationic nickel-methoxide complexes. This study identifies strategies to further lower the barriers for methane C−H activation through evaluation of supporting ligand modifications, solvent polarity, overall charge of complex, metal identity and counterion effects. Overall, neutral low coordinate complexes (e.g. bipyridine) are calculated to have lower activation barriers than the cationic complexes. For both neutral and cationic complexes, the methane C−H activation proceed via a σ-bond metathesis rather than an oxidative addition/reductive elimination pathway. Neutralizing the cationic catalyst models by a counterion, BF4-, has a considerable impact on reducing the methane activation barrier free energy. Third, theoretical studies were performed to explore the effects of appended s-block metal ion crown ethers upon the redox properties of nitridomanganese(V) salen complexes, [(salen)MnV(N)(Mn+-crown ether)]n+, where, M = Na+, K+, Ba2+, Sr2+ for 1Na, 1K, 1Ba, 1Sr complexes respectively; A = complex without Mn+-crown ether and B = without Mn+). The results of the calculations reveal that ΔGrxn(e ̶ ) and thus reduction potentials are quite sensitive to the point charge (q) of the s-block metal ions. Methane activation by A, 1K and 1Ba complexes proceeds via a hydrogen atom abstraction (HAA) pathway with reasonable barriers for all complexes with ~ 4 kcal/mol difference in energy, more favorable free energy barrier for the complexes with higher point charge of metal ion. Changes in predicted properties as a function of continuum solvent dielectric constant suggest that the primary effect of the appended s-block ion is via "through space" interactions. Finally, a comprehensive DFT study of the electrocatalytic oxidation of ammonia to dinitrogen by a ruthenium polypyridyl complex, [(tpy)(bpy)RuII(NH3)]2+ (complex a), and its NMe2-substituted derivative (b), is presented. The thermodynamics and kinetics of electron (ET) and proton transfer (PT) steps and transition states are calculated. NMe2 substitution on bpy reduces the ET steps on average 8 kcal/mol for complex b as compared to a. The calculations indicate that N–N formation occurs by ammonia nucleophilic attack/H-transfer via a nitrene intermediate, rather than a nitride intermediate. Comparison of the free energy profiles of Ru-b with its first-row Fe congener reveals that the thermodynamics are less favorable for the Fe-b model, especially for ET steps. The N-H bond dissociation free energies (BDFEs) for NH3 to form N2 show the following trend: Ru-b <Ru-a <Fe-b, indicating the lowest and most favorable BDFEs for Ru-b complex.