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.

[Fe4S4] cubane in sulfite reductases: new insights into bonding properties and reactivity.

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe4S4] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe4S4] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe4S4] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe4S4] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe4S4] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

Electronic Structure of RhO2+, Its Ammoniated Complexes (NH3)1-5RhO2+, and Mechanistic Exploration of CH4 Activation by Them.

High-level electronic structure calculations are initially performed to investigate the electronic structure of RhO2+. The construction of potential energy curves for the ground and low-lying excited states allowed the calculation of spectroscopic constants, including harmonic and anharmonic vibrational frequencies, bond lengths, spin-orbit constants, and excitation energies. The equilibrium electronic configurations were used for the interpretation of the chemical bonding. We further monitored how the Rh-O bonding scheme changes with the gradual addition of ammonia ligands. The nature of this bond remains unaffected up to four ammonia ligands but adopts a different electronic configuration in the pseudo-octahedral geometry of (NH3)5RhO2+. This has consequences in the activation mechanism of the C-H bond of methane by these complexes, especially (NH3)4RhO2+. We show that the [2 + 2] mechanism in the (NH3)4RhO2+ case has a very low energy barrier comparable to that of a radical mechanism. We also demonstrate that methane can coordinate to the metal in a similar fashion to ammonia and that knowledge of the electronic structure of the pure ammonia complexes provides qualitative insights into the optimal reaction mechanism.

Scandium in Neutral and Positively Charged Ammonia Complexes: Balancing between Sc2+ and Sc3.

High-level quantum chemical approaches are performed to study the stability and electronic structure of tri-, di-, monocationic, and neutral scandium ammonia complexes. The calculated binding energies of all Sc(NH3)1-83+,2+,+,0 complexes reveal the higher stability of hexa- and octacoordinated systems. The ground states of Sc(NH3)6,82+ and Sc(NH3)6,8 have a Sc2+(3d1) center, while there are two competitive electronic states for Sc(NH3)6,8+ with a Sc2+(3d1) or a Sc3+ center. The remaining electrons occupy an outer diffuse s-type orbital (1s). The lower lying states involve 3d-3d transitions for Sc(NH3)6,82+ but outer 1s-1p transitions for Sc(NH3)6,8+,0. The addition of one electron to Sc(NH3)6,83+,2+,+ reduces the binding energies but shortens the Sc-N bond lengths. The comparison with the vanadium and yttrium ammonia complexes (studied earlier) reveal the unique identity of scandium balancing between a d- and s-block element.

Methane to Methanol Conversion Facilitated by Transition-Metal Methyl and Methoxy Units: The Cases of FeCH3+ and FeOCH3.

The conversion of methane to methanol (MTM) catalyzed by FeOCH3+ and FeCH3+ is investigated by means of multireference configuration interaction (MRCI), single-reference coupled clusters (CC), and density functional theory (DFT) approaches. Our dual purpose is the assessment of the applied methodologies and the performance of the proposed catalytic cycle, which involves both of the titled units. The investigated cycle aims to bypass the limitations of metal-oxide catalysts and offers an alternative promising method for efficient MTM transformation. From the technical viewpoint, we found that generally accurate electron correlation treatment is more important than accurately calculated geometries. The combination of optimal DFT geometries with MRCI and CC energetics provides a good compromise between accuracy and efficiency, although there are cases where multireference calculations must be used to obtain correct structures.

Ligand field effects on the ground and excited states of reactive FeO2+ species.

High-valent Fe(iv)-oxo species have been found to be key oxidizing intermediates in the mechanisms of mononuclear iron heme and non-heme enzymes that can functionalize strong C-H bonds. Biomimetic Fe(iv)-oxo molecular complexes have been successfully synthesized and characterized, but their catalytic reactivity is typically lower than that of the enzymatic analogues. The C-H activation step proceeds through two competitive mechanisms, named σ- and π-channels. We have performed high-level wave function theory calculations on bare FeO2+ and a series of non-heme Fe(iv)-oxo model complexes in order to elucidate the electronic properties and the ligand field effects on those channels. Our results suggest that a coordination environment formed by a weak field gives access to both competitive channels, yielding more reactive Fe(iv)-oxo sites. In contrast, a strong ligand environment stabilizes only the σ-channel. Our concluding remarks will aid the derivation of new structure-reactivity descriptors that can contribute to the development of the next generation of functional catalysts.

The role of O(1D) in the oxidation mechanism of ethylene by iodosobenzene and other hypervalent molecules.

The role of the first excited state of oxygen (1D) is proven essential for the description of terminal iodine-oxygen chemical bonds. The description of the I-O bond as a dative one from iodine to O(1D) provides a simple and accurate picture which explains the oxidation properties of iodosobenzene and similar in nature molecules.

A fresh perspective on metal ammonia molecular complexes and expanded metals: opportunities in catalysis and quantum information.

Recent advances in our comprehension of the electronic structure of metal ammonia complexes have opened avenues for novel materials with diffuse electrons. These complexes in their ground state can host peripheral "Rydberg" electrons which populate a hydrogenic-type shell model imitating atoms. Aggregates of such complexes form the so-called expanded or liquid metals. Expanded metals composed of d- and f-block metal ammonia complexes offer properties, such as magnetic moments and larger numbers of diffuse electrons, not present for alkali and alkaline earth (s-block) metals. In addition, tethering metal ammonia complexes via hydrocarbon chains (replacement of ammonia ligands with diamines) yields materials that can be used for redox catalysis and quantum computing, sensing, and optics. This perspective summarizes the recent findings for gas-phase isolated metal ammonia complexes and projects the obtained knowledge to the condensed phase regime. Possible applications for the newly introduced expanded metals and linked solvated electrons precursors are discussed and future directions are proposed.

Impact of Graphene Quantum Dot Edge Morphologies on Their Optical Properties.

The optoelectronic properties of functionalized graphene quantum dots (GQDs) have been explored by simulating electronic structure of three different shapes of GQDs containing exclusively zigzag or armchair edges in both pristine and functionalized forms. Absorption spectra and transition densities for the low-lying excited states are evaluated by using time-dependent density functional theory and compared for different functionalization species. The functionalization position dictates the optical properties of square GQDs, where isomers with CH2 in the intermediate positions (excluding corner and center positions) have higher electronic transition energies and exciton delocalization than other isomers. Rhombic GQDs with all armchair edges exhibit high steric flexibility, and their complete passivation results in the largest structural deformation from planarity and strongest red-shifts. A steady red-shift in the absorption energy is observed following the order F, CH3, Cl, and Br substitutions. This suggests that the steric effects due to large van der Waals radii overcome electronegative effects.

Aufbau Rules for Solvated Electron Precursors: Be(NH3)40,± Complexes and Beyond.

Tetra-amino beryllium complexes and ions, Be(NH3)40,±, have a tetrahedral Be(NH3)42+ core with one, two, or three outer electrons orbiting its periphery. Our calculations reveal a new class of molecular entities, solvated electron precursors, with Aufbau rules (1s, 1p, 1d, 2s, 1f, 2p, 2d) that differ from their familiar hydrogenic counterparts and resemble those of jellium or nuclear-shell models. The core's radial electrostatic potential suffices to reproduce the chief features of the ab initio results. Wave function and electron-propagator methods combined with diffuse basis sets are employed to calculate accurate geometries, ionization energies, electron affinities, and excitation energies.