Controlling product selectivity during dioxygen reduction with Mn complexes using pendent proton donor relays and added base

The catalytic reduction of dioxygen (O2) is important in biological energy conversion and alternative energy applications. In comparison to Fe- and Co-based systems, examples of catalytic O2 reduction by homogeneous Mn-based systems is relatively sparse. Motivated by this lack of knowledge, two Mn-based catalysts for the oxygen reduction reaction (ORR) containing a bipyridine-based non-porphyrinic ligand framework have been developed to evaluate how pendent proton donor relays alter activity and selectivity for the ORR, where Mn(p-tbudhbpy)Cl (1) was used as a control complex and Mn(nPrdhbpy)Cl (2) contains a pendent -OMe group in the secondary coordination sphere. Using an ammonium-based proton source, N,N'-diisopropylethylammonium hexafluorophosphate, we analyzed catalytic activity for the ORR: 1 was found to be 64% selective for H2O2 and 2 is quantitative for H2O2, with O2 binding to the reduced Mn(ii) center being the rate-determining step. Upon addition of the conjugate base, N,N'-diisopropylethylamine, the observed catalytic selectivity of both 1 and 2 shifted to H2O as the primary product. Interestingly, while the shift in selectivity suggests a change in mechanism for both 1 and 2, the catalytic activity of 2 is substantially enhanced in the presence of base and the rate-determining step becomes the bimetallic cleavage of the O-O bond in a Mn-hydroperoxo species. These data suggest that the introduction of pendent relay moieties can improve selectivity for H2O2 at the expense of diminished reaction rates from strong hydrogen bonding interactions. Further, although catalytic rate enhancements are observed with a change in product selectivity when base is added to buffer proton activity, the pendent relays stabilize dimer intermediates, limiting the maximum rate.

Dioxygen Reactivity of Copper(I)/Manganese(II)-Porphyrin Assemblies: Mechanistic Studies and Cooperative Activation of O2

The oxidation of transition metals such as manganese and copper by dioxygen (O₂) is of great interest to chemists and biochemists for fundamental and practical reasons. In this report, the O₂ reactivities of 1:1 and 1:2 mixtures of [(TPP)Mn^II] (1; TPP: Tetraphenylporphyrin) and [(tmpa)Cu^I(MeCN)]⁺ (2; TMPA: Tris(2-pyridylmethyl)amine) in 2-methyltetrahydrofuran (MeTHF) are described. Variable-temperature (-110 °C to room temperature) absorption spectroscopic measurements support that, at low temperature, oxygenation of the (TPP)Mn/Cu mixtures leads to rapid formation of a cupric superoxo intermediate, [(tmpa)Cu^II(O₂^•-)]⁺ (3), independent of the presence of the manganese porphyrin complex (1). Complex 3 subsequently reacts with 1 to form a heterobinuclear μ-peroxo species, [(tmpa)Cu^II-(O₂^2-)-Mn^III(TPP)]⁺ (4; λ_max = 443 nm), which thermally converts to a μ-oxo complex, [(tmpa)Cu^II-O-Mn^III(TPP)]⁺ (5; λ_max = 434 and 466 nm), confirmed by electrospray ionization mass spectrometry and nuclear magnetic resonance spectroscopy. In the 1:2 (TPP)Mn/Cu mixture, 4 is subsequently attacked by a second equivalent of 3, giving a bis-μ-peroxo species, i.e., [(tmpa)Cu^II-(O₂^2-)-Mn^IV(TPP)-(O₂^2-)-Cu^II(tmpa)]²⁺ (7; λ_max = 420 nm and δ_pyrrolic = -44.90 ppm). The final decomposition product of the (TPP)Mn/Cu/O₂ chemistry in MeTHF is [(TPP)Mn^III(MeTHF)₂]⁺ (6), whose X-ray structure is also presented and compared to literature analogs.

A DMRG/CASPT2 Investigation of Metallocorroles: Quantifying Ligand Noninnocence in Archetypal 3d and 4d Element Derivatives

Hybrid density functional theory (B3LYP) and density matrix renormalization group (DMRG) theory have been used to quantitatively compare the degree of ligand noninnocence (corrole radical character) in seven archetypal metallocorroles. The seven complexes, in decreasing order of corrole noninnocent character, are Mn[Cor]Cl > Fe[Cor]Cl > Fe[Cor](NO) > Mo[Cor]Cl2 > Ru[Cor](NO) ≈ Mn[Cor]Ph ≈ Fe[Cor]Ph ≈ 0, where [Cor] refers to the unsubstituted corrolato ligand. DMRG-based second-order perturbation theory calculations have also yielded detailed excited-state energetics data on the compounds, shedding light on periodic trends involving middle transition elements. Thus, whereas the ground state of Fe[Cor](NO) (S = 0) is best described as a locally S = 1/2 {FeNO}7 unit antiferromagnetically coupled to a corrole A' radical, the calculations confirm that Ru[Cor](NO) may be described as simply {RuNO}6-Cor3-, that is, having an innocent corrole macrocycle. Furthermore, whereas the ferromagnetically coupled S = 1{FeNO}7-Cor•2- state of Fe[Cor](NO) is only ∼17.5 kcal/mol higher than the S = 0 ground state, the analogous triplet state of Ru[Cor](NO) is higher by a far larger margin (37.4 kcal/mol) relative to the ground state. In the same vein, Mo[Cor]Cl2 exhibits an adiabatic doublet-quartet gap of 36.1 kcal/mol. The large energy gaps associated with metal-ligand spin coupling in Ru[Cor](NO) and Mo[Cor]Cl2 reflect the much greater covalent character of 4d-π interactions relative to analogous interactions involving 3d orbitals. As far as excited-state energetics is concerned, DMRG-CASPT2 calculations provide moderate validation for hybrid density functional theory (B3LYP) for qualitative purposes, but underscore the possibility of large errors (>10 kcal/mol) in interstate energy differences.

Reactivity of the Bicyclic Amido-Substituted Silicon(I) Ring Compound Si4 {N(SiMe3 )Mes}4 with FLP-Type Character

The bicyclic amido‐substituted silicon(I) ring compound Si₄{N(SiMe₃)Mes}₄2 (Mes=Mesityl=2,4,6‐Me₃C₆H₂) features enhanced zwitterionic character and different reactivity from the analogous compound Si₄{N(SiMe₃)Dipp}₄1 (Dipp=2,6‐ⁱPr₂C₆H₃) due to the smaller mesityl substituents. In a reaction with the N‐heterocyclic carbene NHCMe4 (1,3,4,5‐tetramethyl‐imidazol‐2‐ylidene), we observe adduct formation to give Si₄{N(SiMe₃)Mes}₄ ⋅ NHCMe4 (3). This adduct reacts further with the Lewis acid BH₃ to yield the Lewis acid–base complex Si₄{N(SiMe₃)Mes}₄ ⋅ NHCMe4  ⋅ BH₃ (4). Coordination of AlBr₃ to 2 leads to the adduct 5. Calculated proton affinities and fluoride ion affinities reveal highly Lewis basic and very weak Lewis acidic character of the low‐valent silicon atoms in 1 and 2. This is confirmed by protonation of 1 and 2 with Brookharts acid yielding 6 and 7. Reaction with diphenylacetylene only occurs at 111 °C with 2 in toluene and is accompanied by fragmentation of 2 to afford the silacyclopropene 8 and the trisilanorbornadiene species 9.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

Cycloadditions with a Stable Charge-Separated Cyclobutadiene-Type Amido-Substituted Silicon Ring Compound

Reductive debromination of {N(SiMe₃)₂}SiBr₃ with Rieke magnesium results in the formation of the five‐vertex silicon cluster with one bromine substituent Si₅{N(SiMe₃)₂}₅Br, 1, and the cyclobutadiene analogue 2 in a 1:1 ratio. The latter features a planar four‐membered silicon ring with a charge‐separated electronic situation. Two silicon atoms in 2 are trigonal planar and the other two trigonal pyramidal. In cycloadditions with ethylene, diethylacetylene, 1,5‐cyclooctadiene, and 2,3‐dimethyl‐1,3‐butadiene cyclic unsaturated ring compounds (3–6) were formed at room temperature in quantitative reactions. Two of the products (3 and 6) show photochemical isomerization with LED light (λ=405 nm) to afford saturated ring compounds 4 e and 6′.

Nitric Oxide Dioxygenation by O2 Adducts of Manganese Porphyrins

We describe here nitric oxide dioxygenation (NOD) by the dioxygen manganese porphyrin adducts Mn(Por)(η²-O₂) (Por^2- = the meso-tetra-phenyl or meso-tetra-p-tolylporphyrinato dianions, TPP^2- and TTP^2-). The Mn(Por)(η²-O₂) was assembled by adding O₂ to sublimed layers of Mn^II(Por). When NO was introduced and the temperature was slowly raised from 80 to 120 K, new IR bands with correlated intensities grew concomitant with depletion of the υ(O₂) band. Isotope labeling experiments with ¹⁸O₂, ¹⁵NO, and N¹⁸O combined with DFT calculations provide the basis for identifying the initial intermediates as the six-coordinate peroxynitrito complexes (ON)Mn(Por)(η¹-OONO). Further warming to room temperature led to formation of the nitrato complexes Mn(Por)(η¹-ONO₂), thereby demonstrating the ability of these metal centers to promote NOD. However, comparable quantities of the nitrito complexes Mn(Por)(η¹-ONO) are also formed. In contrast, when the analogous reactions were initiated with the weak σ-donor ligand tetrahydrofuran or dimethyl sulfide present in the layers, formation of Mn(Por)(η¹-ONO₂) is strongly favored (∼90%). The latter are formed via a 6-coordinate intermediate (L)Mn(Por)(η¹-ONO₂) (L = THF or DMS) that loses L upon warming. These reaction patterns are compared to those observed previously with analogous iron and cobalt porphyrin complexes.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

Oxidation of an indole substrate by porphyrin iron(iii) superoxide: relevance to indoleamine and tryptophan 2,3-dioxygenases

Reaction of FeIII(O2˙-)(TPP) with 2,3-dimethylindole at -40 °C gives the ring-opened, dioxygenated N-(2-acetyl-phenyl)-acetamide product. The reaction was monitored in situ by low-temperature UV-vis and 1H NMR spectroscopies. This work demonstrates that a discrete iron(iii)(superoxo) porphyrin is competent to carry out indole oxidation, as proposed for the tryptophan and indoleamine 2,3-dioxygenases.

Transition Metal Spin-State Energetics by MC-PDFT with High Local Exchange

The energetics of the spin states of transition metal complexes have been explored with a variety of electronic structure methods, but the calculations require a compromise between accuracy and affordability. In this work, the spin splittings of several iron complexes are studied with multiconfiguration pair-density functional theory (MC-PDFT). The results are compared to previously published results obtained by complete active space second-order perturbation theory (CASPT2) and CASPT2 with coupled-cluster semicore correlation (CASPT2/CC). In contrast to CASPT2's systematic overstabilization of high-spin states with respect to the CASPT2/CC reference, MC-PDFT with the tPBE on-top functional understabilizes high-spin states. This systematic understabilization is largely corrected by revising the exchange and correlation contributions to the on-top functional using the high local-exchange approximation (tPBE-HLE). Moreover, tPBE-HLE correctly predicts the spin of the ground state in most cases, while CASPT2 incorrectly predicts high-spin ground states in all cases. This is encouraging for practical work because tPBE and tPBE-HLE are faster than CASPT2 by a factor of 50 even in a moderately sized example.

Electronic Structure of N-Bridged High-Valent Diiron-Oxo

Density functional theory (DFT) and an advanced ab initio technique based on density matrix renormalization group (DMRG-CASPT2) were employed to investigate a reactive N-bridged high-valent diiron-oxo species involved in H-abstraction reactions. We studied in detail two important doublet states, the ground state with two iron(IV) centers and a mixed valence Fe^V -Fe^IV excited state. We found that the latter state is low-lying. Furthermore, its electronic structure and spin density imply that it has significantly higher H-abstraction reactivity than the ground state. This low-lying excited state might be the reason behind the high oxidation reactivity of this diiron-oxo species towards methane.