Automated reaction path searches for spin-forbidden reactions

Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as "spin-forbidden" reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin-forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed-spin eigenstate obtained from the diagonalization of the spin-coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the ^6,4 FeO⁺  + H₂ → ^6,4 Fe⁺  + H₂ O, ^6,4 FeO⁺  + CH₄ → ^6,4 Fe⁺  + CH₃ OH, and ⁷ Mn⁺  + OCS → ⁵ MnS⁺  + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.

Magnetic and structural correlations in [Fe(nsal2trien)] salts: the role of cation–anion interactions in the spin crossover phenomenon

Cation–anion and cation–solvent–anion interactions determine the SCO behaviour of six [Fe^III(nsal₂trien)] salts.

Electronic Effects on Room-Temperature, Gas-Phase C-H Bond Activations by Cluster Oxides and Metal Carbides: The Methane Challenge

This Perspective discusses a story of one molecule (methane), a few metal-oxide cationic clusters (MOCCs), dopants, metal-carbide cations, oriented-electric fields (OEFs), and a dizzying mechanistic landscape of methane activation! One mechanism is hydrogen atom transfer (HAT), which occurs whenever the MOCC possesses a localized oxyl radical (M-O^•). Whenever the radical is delocalized, e.g., in [MgO]_n^•+ the HAT barrier increases due to the penalty of radical localization. Adding a dopant (Ga₂O₃) to [MgO]₂^•+ localizes the radical and HAT transpires. Whenever the radical is located on the metal centers as in [Al₂O₂]^•+ the mechanism crosses over to proton-coupled electron transfer (PCET), wherein the positive Al center acts as a Lewis acid that coordinates the methane molecule, while one of the bridging oxygen atoms abstracts a proton, and the negatively charged CH₃ moiety relocates to the metal fragment. We provide a diagnostic plot of barriers vs reactants' distortion energies, which allows the chemist to distinguish HAT from PCET. Thus, doping of [MgO]₂^•+ by Al₂O₃ enables HAT and PCET to compete. Similarly, [ZnO]^•+ activates methane by PCET generating many products. Adding a CH₃CN ligand to form [(CH₃CN)ZnO]^•+ leads to a single HAT product. The CH₃CN dipole acts as an OEF that switches off PCET. [MC]⁺ cations (M = Au, Cu) act by different mechanisms, dictated by the M⁺-C bond covalence. For example, Cu⁺, which bonds the carbon atom mostly electrostatically, performs coupling of C to methane to yield ethylene, in a single almost barrier-free step, with an unprecedented atomic choreography catalyzed by the OEF of Cu⁺.

Implementation and performance of the artificial force induced reaction method in the GRRM17 program

This article reports implementation and performance of the artificial force induced reaction (AFIR) method in the upcoming 2017 version of GRRM program (GRRM17). The AFIR method, which is one of automated reaction path search methods, induces geometrical deformations in a system by pushing or pulling fragments defined in the system by an artificial force. In GRRM17, three different algorithms, that is, multicomponent algorithm (MC-AFIR), single-component algorithm (SC-AFIR), and double-sphere algorithm (DS-AFIR), are available, where the MC-AFIR was the only algorithm which has been available in the previous 2014 version. The MC-AFIR does automated sampling of reaction pathways between two or more reactant molecules. The SC-AFIR performs automated generation of global or semiglobal reaction path network. The DS-AFIR finds a single path between given two structures. Exploration of minimum energy structures within the hypersurface in which two different electronic states degenerate, and an interface with the quantum mechanics/molecular mechanics method, are also described. A code termed SAFIRE will also be available, as a visualization software for complicated reaction path networks. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

On the Electronic Origin of Remarkable Ligand Effects on the Reactivities of [NiL]+ Complexes (L=C6 H5 , C5 H4 N, CN) towards Methane

The gas-phase reactions of [NiL]⁺ (L=C₆ H₅ , C₅ H₄ N, CN) with methane have been explored by using electrospray-ionization mass spectrometry (ESI-MS) complemented by quantum chemical calculations. Though the phenyl Ni complex [Ni(C₆ H₅ )]⁺ exclusively abstracts one hydrogen atom from methane at ambient conditions, the cyano Ni complex [Ni(CN)]⁺ brings about both H-atom abstraction and ligand exchange to generate [Ni(CH₃ )]⁺ . In contrast, the complex 2-pyridinyl Ni [Ni(C₅ H₄ N)]⁺ is inert towards this substrate. The presence of the empty 4s(Ni) orbital dominates the proton-coupled electron transfer (PCET) processes for the investigated systems.

Quantum-Mechanical Study of the Reaction Mechanism for 2π-2π Cycloaddition of Fluorinated Methylene Groups

Perfluorocyclobutyl polymers are thermally and chemically stable, may be produced without a catalyst via thermal 2π-2π cycloaddition, and can form block structures, making them suitable for commercialization of specialty polymers. Thermal 2π-2π cycloaddition is a rare reaction that begins in the singlet state and proceeds through a triplet intermediate to form an energetically stable four-membered ring in the singlet state. This reaction involves two changes in spin state and, thus, two spin-crossover transitions. Presented here are density functional theory calculations that evaluate the energetics and reaction mechanisms for the dimerizations of two different polyfluorinated precursors, 1,1,2-trifluoro-2-(trifluoromethoxy)ethane and hexafluoropropylene. The spin-crossover transition states are thoroughly investigated, revealing important kinetics steps and an activation energy for the gas-phase cycloaddition of two hexafluoropropene molecules of 36.9 kcal/mol, which is in good agreement with the experimentally determined value of 34.3 kcal/mol. It is found that the first carbon-carbon bond formation is the rate-limiting step, followed by a rotation about the newly formed bond in the triplet state that results in the formation of the second carbon-carbon bond. Targeting the rotation of the C-C bond, a set of parameters were obtained that best produce high molecular weight polymers using this chemistry.

Facile and Reversible Formation of Iron(III)-Oxo-Cerium(IV) Adducts from Nonheme Oxoiron(IV) Complexes and Cerium(III)

Ceric ammonium nitrate (CAN) or Ce^IV (NH₄ )₂ (NO₃ )₆ is often used in artificial water oxidation and generally considered to be an outer-sphere oxidant. Herein we report the spectroscopic and crystallographic characterization of [(N4Py)Fe^III -O-Ce^IV (OH₂ )(NO₃ )₄ ]⁺ (3), a complex obtained from the reaction of [(N4Py)Fe^II (NCMe)]²⁺ with 2 equiv CAN or [(N4Py)Fe^IV =O]²⁺ (2) with Ce^III (NO₃ )₃ in MeCN. Surprisingly, the formation of 3 is reversible, the position of the equilibrium being dependent on the MeCN/water ratio of the solvent. These results suggest that the Fe^IV and Ce^IV centers have comparable reduction potentials. Moreover, the equilibrium entails a change in iron spin state, from S=1 Fe^IV in 2 to S=5/2 in 3, which is found to be facile despite the formal spin-forbidden nature of this process. This observation suggests that Fe^IV =O complexes may avail of reaction pathways involving multiple spin states having little or no barrier.

Nonclassical Single-State Reactivity of an Oxo-Iron(IV) Complex Confined to Triplet Pathways

C-H bond activation mediated by oxo-iron (IV) species represents the key step of many heme and nonheme O₂-activating enzymes. Of crucial interest is the effect of spin state of the Fe^IV(O) unit. Here we report the C-H activation kinetics and corresponding theoretical investigations of an exclusive tetracarbene ligated oxo-iron(IV) complex, [L^NHCFe^IV(O)(MeCN)]²⁺ (1). Kinetic traces using substrates with bond dissociation energies (BDEs) up to 80 kcal mol^-1 show pseudo-first-order behavior and large but temperature-dependent kinetic isotope effects (KIE 32 at -40 °C). When compared with a topologically related oxo-iron(IV) complex bearing an equatorial N-donor ligand, [L^TMCFe^IV(O) (MeCN)]²⁺ (A), the tetracarbene complex 1 is significantly more reactive with second order rate constants k'₂ that are 2-3 orders of magnitude higher. UV-vis experiments in tandem with cryospray mass spectrometry evidence that the reaction occurs via formation of a hydroxo-iron(III) complex (4) after the initial H atom transfer (HAT). An extensive computational study using a wave function based multireference approach, viz. complete active space self-consistent field (CASSCF) followed by N-electron valence perturbation theory up to second order (NEVPT2), provided insight into the HAT trajectories of 1 and A. Calculated free energy barriers for 1 reasonably agree with experimental values. Because the strongly donating equatorial tetracarbene pushes the Fe-d_x²-y² orbital above d_z², 1 features a dramatically large quintet-triplet gap of ∼18 kcal/mol compared to ∼2-3 kcal/mol computed for A. Consequently, the HAT process performed by 1 occurs on the triplet surface only, in contrast to complex A reported to feature two-state-reactivity with contributions from both triplet and quintet states. Despite this, the reactive Fe^IV(O) units in 1 and A undergo the same electronic-structure changes during HAT. Thus, the unique complex 1 represents a pure "triplet-only" ferryl model.

Sequential Gas-Phase Activation of Carbon Dioxide and Methane by [Re(CO)2]+: The Sequence of Events Matters!

The potential of carbonyl rhenium complexes in activating and coupling carbon dioxide and methane has been explored by using a combination of gas-phase experiments (FT-ICR mass spectrometry) and high-level quantum chemical calculations. While the complexes [Re(CO)_x]⁺ (x = 0, 1, 3) are thermally unreactive toward CO₂, [Re(CO)₂]⁺ abstracts one oxygen atom from this substrate spontaneously at ambient conditions. Based on ¹³C and ¹⁸O labeling experiments, the newly generated CO ligand is preferentially eliminated, and two mechanistic scenarios are considered to account for this unexpected finding. The oxo complex [ORe(CO)₂]⁺ reacts further with CH₄ to produce the dihydridomethylene complex [ORe(CO)(CH₂)(H)₂]⁺. However, coupling of the CO and CH₂ ligands to form CH₂═C═O does not take place. Further, the complexes [Re(CO)_x(CH₂)]⁺ (x = 1, 2), generated in the thermal reaction of [Re(CO)_x]⁺ (x = 1, 2) with CH₄, are inert toward CO₂. Mechanistic insight on the origin of this remarkable reactivity pattern has been derived from detailed quantum chemical calculations.

Spin-Forbidden Branching in the Mechanism of the Intrinsic Haber-Weiss Reaction

The mechanism of the O₂^⋅− and H₂O₂ reaction (Haber–Weiss) under solvent‐free conditions has been characterized at the DFT and CCSD(T) level of theory to account for the ease of this reaction in the gas phase and the formation of two different set of products (Blanksby et al., Angew. Chem. Int. Ed. 2007, 46, 4948). The reaction is shown to proceed through an electron‐transfer process from the superoxide anion to hydrogen peroxide, along two pathways. While the O₃^⋅− + H₂O products are formed from a spin‐allowed reaction (on the doublet surface), the preferred products, O^⋅−(H₂O)+³O₂, are formed through a spin‐forbidden reaction as a result of a favorable crossing point between the doublet and quartet surface. Plausible reasons for the preference toward the latter set are given in terms of the characteristics of the minimum energy crossing point (MECP) and the stability of an intermediate formed (after the MECP) in the quartet surface. These unique results show that these two pathways are associated with a bifurcation, yielding spin‐dependent products.

On the Origin of Reactivity Enhancement/Suppression upon Sequential Ligation: [Re(CO)x ]+ /CH4 (x=0-3) Couples

The thermal gas-phase reactions of rhenium carbonyl complexes [Re(CO)_x ]⁺ (x=0-3) with methane have been explored by using FT-ICR mass spectrometry complemented by high-level quantum chemical calculation. While it had been concluded in previous studies that addition of closed-shell ligands in general decreases the reactivity of metal ions, the current work provides an exception: the previously demonstrated inertness of atomic Re⁺ towards methane is completely changed upon ligation with CO. Both [Re(CO)]⁺ and [Re(CO)₂ ]⁺ bring about efficient dehydrogenation of the hydrocarbon under ambient conditions. However, addition of a third ligand to form [Re(CO)₃ ]⁺ completely quenches the reactivity.

Ménage-à-trois: single-atom catalysis, mass spectrometry, and computational chemistry

Genuine, single-atom catalysis can be realized in the gas phase and probed by mass spectrometry combined with computational chemistry.

The open-cubane oxo–oxyl coupling mechanism dominates photosynthetic oxygen evolution: a comprehensive DFT investigation on O–O bond formation in the S4state

How is O₂created in nature? Comprehensive DFT investigations determine the dominance of the open-cubane oxo–oxyl coupling mechanism over alternative possibilities.

The reaction of methyl peroxy and hydroxyl radicals as a major source of atmospheric methanol

Methyl peroxy, a key radical in tropospheric chemistry, was recently shown to react with the hydroxyl radical at an unexpectedly high rate. Here, the molecular reaction mechanisms are elucidated using high-level quantum chemical methodologies and statistical rate theory. Formation of activated methylhydrotrioxide, followed by dissociation into methoxy and hydroperoxy radicals, is found to be the main reaction pathway, whereas methylhydrotrioxide stabilization and methanol formation (from activated and stabilized methylhydrotrioxide) are viable minor channels. Criegee intermediate formation is found to be negligible. Given the theoretical uncertainties, useful constraints on the yields are provided by atmospheric methanol measurements. Using a global chemistry-transport model, we show that the only explanation for the high observed methanol abundances over remote oceans is the title reaction with an overall methanol yield of ∼30%, consistent with the theoretical estimates given their uncertainties. This makes the title reaction a major methanol source (115 Tg per year), comparable to global terrestrial emissions.

Penetrating the Elusive Mechanism of Copper-Mediated Fluoromethylation in the Presence of Oxygen through the Gas-Phase Reactivity of Well-Defined [LCuO](+) Complexes with Fluoromethanes (CH(4-n)Fn, n = 1-3)

Traveling wave ion mobility spectrometry (TWIMS) isomer separation was exploited to react the particularly well-defined ionic species [LCuO](+) (L = 1,10-phenanthroline) with the neutral fluoromethane substrates CH(4-n)Fn (n = 1-3) in the gas phase. Experimentally, the monofluoromethane substrate (n = 1) undergoes both hydrogen-atom transfer, forming the copper hydroxide complex [LCuOH](•+) and concomitantly a CH2F(•) radical, and oxygen-atom transfer, yielding the observable ionic product [LCu](+) plus the neutral oxidized substrate [C,H3,O,F]. DFT calculations reveal that the mechanism for both product channels relies on the initial C-H bond activation of the substrate. Compared to nonfluorinated methane, the addition of fluorine to the substrate assists the reactivity through a lowering of the C-H bond energy and reaction preorganization (through noncovalent interaction in the encounter complex). A two-state reactivity scenario is mandatory for the oxidation, which competitively results in the unusual fluoromethanol product, CH2FOH, or the decomposed products, CH2O and HF, with the latter channel being kinetically disfavored. Difluoromethane (n = 2) is predicted to undergo the analogous reactions at room temperature, although the reactions are less favored than those of monofluoromethane. The reaction of trifluoromethane (n = 3, fluoroform) through C-H activation is kinetically hindered under ambient conditions but might be expected to occur in the condensed phase upon heating or with further lowering of reaction barriers through templation with counterions, such as potassium. Overall, formation of CH(3-n)Fn(•) and CH(3-n)FnOH occurs under relatively gentle energetic conditions, which sheds light on their potential as reactive intermediates in fluoromethylation reactions mediated by copper in the presence of oxygen.

Photocycloaddition reaction of atropisomeric maleimides: mechanism and selectivity

We report a density functional study on the mechanism of the [2+2] photocyclization of atropisomeric maleimides.

Oxygen Reduction with a Bifunctional Iridium Dihydride Complex

The iridium dihydride [Ir(H)2 (HPNP)](+) (PNP=N(CH2 CH2 PtBu2 )2 ) reacts with O2 to give the unusual, square-planar iridium(III) hydroxide [Ir(OH)(PNP)](+) and water. Regeneration of the dihydride with H2 closes a quasi-catalytic synthetic oxygen-reduction reaction (ORR) cycle that can be run several times. Experimental and computational examinations are in agreement with an oxygenation mechanism via rate-limiting O2 coordination followed by H-transfer at a single metal site, facilitated by the cooperating pincer ligand. Hence, the four electrons required for the ORR are stored within the two covalent MH bonds of a mononuclear metal complex.

Dioxygen Activation by Non-Adiabatic Oxidative Addition to a Single Metal Center

A chromium(I) dinitrogen complex reacts rapidly with O2 to form the mononuclear dioxo complex [Tp(tBu,Me)Cr(V)(O)2] (Tp(tBu,Me) = hydrotris(3-tert-butyl-5-methylpyrazolyl)borate), whereas the analogous reaction with sulfur stops at the persulfido complex [Tp(tBu,Me)Cr(III)(S2)]. The transformation of the putative peroxo intermediate [Tp(tBu,Me)Cr(III)(O2)] (S = 3/2) into [Tp(tBu,Me)Cr(V)(O)2] (S = 1/2) is spin-forbidden. The minimum-energy crossing point for the two potential energy surfaces has been identified. Although the dinuclear complex [(Tp(tBu,Me)Cr)2(μ-O)2] exists, mechanistic experiments suggest that O2 activation occurs on a single metal center, by an oxidative addition on the quartet surface followed by crossover to the doublet surface.