Theoretical calculations for cobaltoporphyrin with nitric oxide as axial ligand

Crystal Structure and Electron Configuration of {MNO}8 Heme Complexes

Although research on nitrosyl (NO) heme complexes and their one-electron reduced form, nitroxyl (or nitroxyl anion, NO-) derivatives, has been going on for decades, there are still disagreements about the electrical configuration of nitroxyl complexes, and the majority of the work on this topic is based on theoretical calculations. Following the initial nitroxyl iron porphyrin crystal structure, we present two further polymorphic forms of [CoCp2][Fe(TFPPBr8)(NO)]. Using the same completely halogenated porphyrin ligand, we also present two polymorphic forms of nitrosyl cobalt(II) complexes, which are another sort of {MNO}8 structure. In addition to the EXANES and EPR studies of these {FeNO}7 and {CoNO}8 complexes, the {FeNO}8 [CoCp2][Fe(TFPPBr8)(NO)] complex is also investigated by temperature-dependent Mössbauer experiments for the first time with the {FeNO}7 precursor as a control sample. The analysis of the Mössbauer and crystal structural parameters between these two types of {MNO}8 (M = Fe or Co) species and previously reported analogous ones allow us to conclude that the electronic configuration of [Fe(TFPPBr8)(NO)]- is best described as an intermediate between low-spin Fe(II)-NO- and Fe(I)-NO•.

Not Limited to Iron: A Cobalt Heme-NO Model Facilitates N-N Coupling with External NO in the Presence of a Lewis Acid to Generate N2 O

Some bacterial heme proteins catalyze the coupling of two NO molecules to generate N2 O. We previously reported that a heme Fe-NO model engages in this N-N bond-forming reaction with NO. We now demonstrate that (OEP)CoII (NO) similarly reacts with 1 equiv of NO in the presence of the Lewis acids BX3 (X=F, C6 F5 ) to generate N2 O. DFT calculations support retention of the CoII oxidation state for the experimentally observed adduct (OEP)CoII (NO⋅BF3 ), the presumed hyponitrite intermediate (P.+ )CoII (ONNO⋅BF3 ), and the porphyrin π-radical cation by-product of this reaction, and that the π-radical cation formation likely occurs at the hyponitrite stage. In contrast, the Fe analogue undergoes a ferrous-to-ferric oxidation state conversion during this reaction. Our work shows that cobalt hemes are chemically competent to engage in the NO-to-N2 O conversion reaction.

Multireference Perturbation Theory with Cholesky Decomposition for the Density Matrix Renormalization Group

We present a second-order N-electron valence state perturbation theory (NEVPT2) based on a density matrix renormalization group (DMRG) reference wave function that exploits a Cholesky decomposition of the two-electron repulsion integrals (CD-DMRG-NEVPT2). With a parameter-free multireference perturbation theory approach at hand, the latter allows us to efficiently describe static and dynamic correlation in large molecular systems. We demonstrate the applicability of CD-DMRG-NEVPT2 for spin-state energetics of spin-crossover complexes involving calculations with more than 1000 atomic basis functions. We first assess, in a study of a heme model, the accuracy of the strongly and partially contracted variant of CD-DMRG-NEVPT2 before embarking on resolving a controversy about the spin ground state of a cobalt tropocoronand complex.

Role of Spin States in Nitric Oxide Binding to Cobalt(II) and Manganese(II) Porphyrins. Is Tighter Binding Always Stronger?

Binding of nitric oxide (NO) to metalloporphyrins and heme groups is important in biochemistry while challenging to describe accurately by density functional theory (DFT) calculations. Here, the structural and thermochemical aspect of NO binding to Co(II) and Mn(II) porphyrins is investigated by DFT and DFT-D (dispersion-corrected) calculations, supported by reliable coupled-cluster methodology (CCSD(T)), and critically correlated with the experimental data. It is argued that whereas the bonding of NO to Co(II) porphyrin is a simple radical recombination, the bonding of NO to Mn(II) porphyrin is accompanied by a crossing of spin states. For this reason, the spin-state conversion energy contributes to the Mn-NO bond energy, and the paradigmatic correlation between bond length and bond energy is violated for the considered nitrosyl complexes: the Mn-NO bond is (structurally) shorter by ∼0.2 Å, albeit (energetically) weaker by a few kcal/mol, compared with the Co-NO bond. Moreover, none of the many tested DFT methods can reproduce the Co-NO and Mn-NO bond energies simultaneously, except for calculations with B3LYP*-D3, TPSSh-D3, and M06-D3 methods supplemented with the proposed spin-state energy correction (to compensate for an error on the calculated spin-state conversion energy). The results of this study are important to appreciate the role of spin-state changes in ligand binding properties of heme-related models. They also highlight the need for accurate calculations for correct interpretation of experimental data, including the qualitative structure-energy relationship.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Electronic structure and ligand vibrations in FeNO, CoNO, and FeOO porphyrin adducts

The gaseous ligands, CO, NO, and O2 interact with the Fe ion in heme proteins largely via backbonding of Fe electrons to the π* orbitals of the XO (X = C, N, O) ligands. In these FeXO adducts, the Fe-X stretching frequency varies inversely with the X-O stretching frequency, since increased backbonding strengthens the Fe-X bond while weakening the X-O bond. Inverse frequency correlations have been observed for all three ligands, despite differing electronic and geometric structures, and despite variable composition of the "FeX" vibrational mode, in which Fe-X stretching and Fe-X-O coordinates are mixed for bent FeXO adducts. We report experimental data for 5-coordinate Co(II)(NO) porphyrin adducts (isoelectronic with Fe(II)(OO) adducts), and the results of density functional theory (DFT) modeling for 5-coordinate Fe(II)(NO), Co(II)(NO), and Fe(II)(OO) adducts. Inverse ν(MX)/ν(XO) correlations are obtained computationally, using model porphyrins with graded electron-donating and -withdrawing substituents to modulate the backbonding. Computed slopes agree satisfactorily with experiment, provided nonhybrid functionals are used, which avoid overemphasizing high-spin states. The BP86 functional gives correct ground states, a closed-shell singlet for Co(II)(NO) and an open-shell singlet for the isoelectronic Fe(II)(OO), as corroborated by structural data for Co(II)(NO), and the ν(MX)/ν(XO) slope agreement with experiment for both adducts. However, for Fe(II)(OO) adducts, the computed inverse ν(MX)/ν(XO) correlation applies only to porphyrins with electron-donating and withdrawing substituents of moderate strength. For substituents more donating than -CH3, a direct correlation is obtained, the Fe-O and O-O bonds weakening in concert. This effect is ascribed to the dominance of σ bonding via the in-plane dxz(+dz(2))-π* orbital, when electron-donating substituents raise the d orbital energies sufficiently to render backbonding (dyz-π*) unimportant.

Cobalt-porphyrin catalyzed electrochemical reduction of carbon dioxide in water. 1. A density functional study of intermediates

The reduction of carbon dioxide by cobalt porphyrins is thought to be a multistep reaction with several possible intermediates and reaction pathways. We here investigate a number of possible intermediates in this reaction using density functional theory, including both hybrid (B3LYP) and pure (PBE and BP86) functionals. Optimum structures are located, and harmonic vibrational frequencies and thermal corrections are computed for the low-lying electronic states for all intermediates. Free energies of solvation are predicted for all species, providing a reaction profile in the aqueous phase, which enables identification of likely pathways. Finally, the reaction energy for the binding of carbon dioxide to the cobalt porphine cation is determined in the gas phase and in solution.

Cobalt-porphyrin catalyzed electrochemical reduction of carbon dioxide in water. 2. Mechanism from first principles

We apply first principles computational techniques to analyze the two-electron, multistep, electrochemical reduction of CO(2) to CO in water using cobalt porphyrin as a catalyst. Density functional theory calculations with hybrid functionals and dielectric continuum solvation are used to determine the steps at which electrons are added. This information is corroborated with ab initio molecular dynamics simulations in an explicit aqueous environment which reveal the critical role of water in stabilizing a key intermediate formed by CO(2) bound to cobalt. By use of potential of mean force calculations, the intermediate is found to spontaneously accept a proton to form a carboxylate acid group at pH < 9.0, and the subsequent cleavage of a C-OH bond to form CO is exothermic and associated with a small free energy barrier. These predictions suggest that the proposed reaction mechanism is viable if electron transfer to the catalyst is sufficiently fast. The variation in cobalt ion charge and spin states during bond breaking, DFT+U treatment of cobalt 3d orbitals, and the need for computing electrochemical potentials are emphasized.