Intrinsic dinitrogen activation at bare metal atoms

There is ongoing interest in metal complexes which bind dinitrogen and facilitate either its reduction or oxidation under mild conditions. In nature, the enzyme nitrogenase catalyzes this process, and dinitrogen fixation occurs under mild and ambient conditions at a metal-sulfur cluster in the center of the MoFe protein, but the mechanism of this process remains largely unknown. In the last few years, new important discoveries have been made in this field. In this review are discussed recent findings on the interaction of N(2) with metal atoms and metal-atom dimers from all groups of the periodic table as provided by gas-phase as well as matrix-isolation experiments. Intrinsic dinitrogen activation at such bare metal atoms is then related to corresponding processes at complexes, clusters, and surfaces.     

Energetics and mechanism of ammonia synthesis through the Chatt Cycle: conditions for a catalytic mode and comparison with the Schrock Cycle

Through a series of DFT calculations the energy profile of the Chatt cycle is evaluated. This is the counterpiece of our earlier investigations of the Schrock cycle (Angew. Chem. 2005, 117, 5783; Angew. Chem. Int. Ed. 2005, 44, 5639), applying the same quantumchemical methodology and approximations. As for the Schrock cycle, decamethylchromocene acts as reductant. The protonation reactions are considered to be mediated by HBF4/diethyl ether or lutidinium. For all protonation and reduction steps the corresponding free reaction enthalpy changes are calculated. The derived energy profile and corresponding reaction mechanism bear strong similarities to the Schrock cycle. In particular, the most endergonic reaction is the first protonation of the N2 complex and the most exergonic reaction is the cleavage of the N--N bond. If lutidinium is employed as acid and Cp2*Cr as reductant, the reaction course involves steps that are not thermally allowed. For HBF4/diethyl ether as the acid and Cp2*Cr as reducant, however, a catalytic cycle consisting of thermally allowed reactions is principally feasible. This cycle involves a Mo I-fluoro complex as dinitrogen intermediate. It is shown that regeneration to the Mo 0-bis(dinitrogen) complex is thermally not accessible in this system. Moreover, the Mo I fluoro-dinitrogen complex is labile towards disproportionation. The implications of these results with respect to the realization of a catalytic system on the basis of Mo and W phosphine complexes are discussed.     

Elucidating the Coordination Chemistry and Mechanism of Biological Nitrogen Fixation

How does the enzyme nitrogenase reduce the inert molecule N₂ to NH₃ under ambient conditions that are so different from the energy‐expensive conditions of the best industrial practices? This review focuses on recent theoretical investigations of the catalytic site, the iron–molybdenum cofactor FeMo‐co, and the way in which it is hydrogenated by protons and electrons and then binds N₂. Density functional calculations provide reaction profiles and activation energies for possible mechanistic steps. This establishes a conceptual framework and the principles for the coordination chemistry of FeMo‐co that are essential to the chemical mechanism of catalysis. The model advanced herein explains relevant experimental data.     

Enzymatic and catalytic reduction of dinitrogen to ammonia: Density functional theory characterization of alternative molybdenum active sites

We used density functional calculations to model dinitrogen reduction by a FeMo cofactor containing a central nitrogen atom and by a Mo‐based catalyst. Plausible intermediates, reaction pathways, and relative energetics in the enzymatic and catalytic reduction of N₂ to ammonia at a single Mo center are explored. Calculations indicate that the binding of N₂ to the Mo atom and the subsequent multiple proton–electron transfer to dinitrogen and its protonated species involved in the conversion of N₂ are feasible energetically. In the reduction of N₂ the Mo atom experiences a cycled oxidation state from Mo(IV) to Mo(VI) by nitrogenase and from Mo(III) to Mo(VI) by the molybdenum catalyst, respectively, tuning the gradual reduction of N₂. Such a wide range of oxidation states exhibited by the Mo center is crucial for the gradual reduction process via successive proton–electron transfer. Present results suggest that the Mo atom in the N‐centered FeMo cofactor is a likely alternative active site for dinitrogen binding and reduction under mild conditions once there is an empty site available at the Mo site. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Car-Parrinello molecular dynamics study of the initial dinitrogen reduction step in Sellmann-type nitrogenase model complexes

We have studied reduction reactions for nitrogen fixation at Sellmann-type model complexes with Car-Parrinello simulation techniques. These dinuclear complexes are especially designed to emulate the so-called open-side FeMoco model. The main result of this work shows that in order to obtain the reduced species several side reactions have to be suppressed. These involve partial dissociation of the chelate ligands and hydrogen atom transfer to the metal center. Working at low temperature turns out to be one necessary pre-requisite in carrying out successful events. The successful events cannot be described by simple reaction coordinates. Complicated processes are involved during the initiation of the reaction. Our theoretical study emphasizes two experimental strategies which are likely to inhibit the side reactions. Clamping of the two metal fragments by a chelating phosphane ligand should prevent dissociation of the complex. Furthermore, introduction of tert-butyl substituents could improve the solubility and should thus allow usage of a wider range of (mild) acids, reductants, and reaction conditions.     

The use of quantum molecular calculations to guide a genetic algorithm: a way to search for new chemistry

The process of gene-based molecular evolution has been simulated in silico by using massively parallel density functional theory quantum calculations, coupled with a genetic algorithm, to test for fitness with respect to a target chemical reaction in populations of genetically encoded molecules. The goal of this study was the identification of transition-metal complexes capable of mediating a known reaction, namely the cleavage of N(2) to give the metal nitride. Each complex within the search space was uniquely specified by a nanogene consisting of an eight-digit number. Propagation of an individual nanogene into successive generations was determined by the fitness of its phenotypic molecule to perform the target reaction and new generations were created by recombination and mutation of surviving nanogenes. In its simplest implementation, the quantum-directed genetic algorithm (QDGA) quickly located a local minimum on the evolutionary fitness hypersurface, but proved incapable of progressing towards the global minimum. A strategy for progressing beyond local minima consistent with the Darwinian paradigm by the use of environmental variations coupled with mass extinctions was therefore developed. This allowed for the identification of nitriding complexes that are very closely related to known examples from the chemical literature. Examples of mutations that appear to be beneficial at the genetic level but prove to be harmful at the phenotypic level are described. As well as revealing fundamental aspects of molecular evolution, QDGA appears to be a powerful tool for the identification of lead compounds capable of carrying out a target chemical reaction.     

Reaction mechanism of apocarotenoid oxygenase (ACO): a DFT study

The mechanism of the oxidative cleavage catalyzed by apocarotenoid oxygenase (ACO) was studied by using a quantum chemical (DFT: B3 LYP) method. Based on the available crystal structure, relatively large models of the unusual active-site region, in which a ferrous ion is coordinated by four histidines and no negatively charged ligand, were selected and used in the computational investigation of the reaction mechanism. The results suggest that binding of dioxygen to the ferrous ion in the active site promotes one-electron oxidation of carotenoid leading to a substrate radical cation and a Fe-bound superoxide radical. Recombination of the two radicals, which can be realized in at least two different ways, yields a reactive peroxo species that subsequently evolves into either a dioxetane or an epoxide intermediate. The former easily decays into the final aldehyde products, whereas the oxidation of the epoxide to the proper products of the reaction requires involvement of a water molecule. The calculated activation barriers favor the dioxetane mechanism, yet the mechanism involving the epoxide intermediate cannot be ruled out.     

An electronic perspective on the reduction of an n=n double bond at a conserved dimolybdenum core

Density functional theory has been used to assess the role of the bimetallic core in supporting reductive cleavage of the N=N double bond in [Cp2Mo2(mu-SMe)3(mu-eta1:eta1-HN=NPh)]+. The HOMO of the complex, the Mo-Mo delta orbital, plays a key role as a source of high-energy electrons, available for transfer into the vacant orbitals of the N=N unit. As a result, the metal centres cycle between the Mo(III) and Mo(IV) oxidation states. The symmetry of the Mo-Mo delta "buffer" orbital has a profound influence on the reaction pathway, because significant overlap with the redox-active orbital on the N=N unit (pi* or sigma*) is required for efficient electron transfer. The orthogonality of the Mo-Mo delta and N-N sigma* orbitals in the eta1:eta1 coordination mode ensures that electron transfer into the N-N sigma bond is effectively blocked, and a rate-limiting eta1:eta1-->eta1 rearrangement is a necessary precursor to cleavage of the bond.     

Theoretical study of N2O reduction by CO in Fe-BEA Zeolite.

Quantum mechanical (QM) and QM/molecular mechanics (MM) studies of the full catalytic cycle of N(2)O reduction by CO in Fe-BEA zeolite, that is, oxidation of BEA-Fe by N(2)O and reduction of BEA-Fe-alphaO by CO, is presented. A large QM cluster, representing half of the channel of the BEA zeolite, is used. The contribution of the MM embedding to the calculated activation energies is found to be negligible. The minimum-energy paths for N(2)O decomposition and reduction with CO are calculated using the nudged elastic band (NEB) method. Calculated and experimental activation energies are in good agreement. The two possible orientations for the gaseous molecules adsorbing on the Fe site that are found lead to different activation energies.     

A photochemical activation scheme of inert dinitrogen by dinuclear Ru(II) and Fe(II) complexes

A general photochemical activation process of inert dinitrogen coordinated to two metal centers is presented on the basis of high-level DFT and ab initio calculations. The central feature of this activation process is the occupation of an antibonding pi* orbital upon electronic excitation from the singlet ground state S0 to the first excited singlet state S1. Populating the antibonding LUMO weakens the triple bond of dinitrogen. After a vertical excitation, the excited complex may structurally relax in the S1 state and approaches its minimum structure in the S1 state. This excited-state minimum structure features the dinitrogen bound in a diazenoid form, which exhibits a double bond and two lone pairs localized at the two nitrogen atoms, ready to be protonated. Reduction and de-excitation then yield the corresponding diazene complex; its generation represents the essential step in a nitrogen fixation and reduction protocol. The consecutive process of excitation, protonation, and reduction may be rearranged in any experimentally appropriate order. The protons needed for the reaction from dinitrogen to diazene can be provided by the ligand sphere of the complexes, which contains sulfur atoms acting as proton acceptors. These protonated thiolate functionalities bring protons close to the dinitrogen moiety. Because protonation does not change the pi*-antibonding character of the LUMO, the universal and well-directed character of the photochemical activation process makes it possible to protonate the dinitrogen complex before it is irradiated. The pi*-antibonding LUMO plays the central role in the activation process, since the diazenoid structure was obtained by excitation from various occupied orbitals as well as by a direct two-electron reduction (without photochemical activation) of the complex; that is, the important bending of N2 towards a diazenoid conformation can be achieved by populating the pi*-antibonding LUMO.     

Lewis adducts of the side-on end-on dinitrogen-bridged complex [{(NPN)Ta}2(mu-H)2(mu-eta 1:eta 2-N2)] with AlMe3, GaMe3, and B(C6F5)3: synthesis, structure, and spectroscopic properties

Reaction of the side-on end-on dinitrogen complex [{(NPN)Ta}(2)(mu-H)(2)(mu-eta(1):eta(2)-N(2))] (1; in which NPN=(PhNSiMe(2)CH(2))(2)PPh), with the Lewis acids XR(3) results in the adducts [{(NPN)Ta}(2)(mu-H)(2)(mu-eta(1):eta(2)-NNXR(3))], XR(3)=GaMe(3) (2), AlMe(3) (3), and B(C(6)F(5))(3) (4). The solid-state molecular structures of 2, 3, and 4 demonstrate that the N-N bond length increases relative to those found in 1 by 0.036, 0.043, and 0.073 A, respectively. In solution complexes 2-4 are fluxional as evidenced by variable-temperature (1)H NMR spectroscopy. The (15)N{(1)H} NMR spectra of 2-4 are reported; furthermore, their vibrational properties and electronic structures are evaluated. The vibrational structures are found to be closely related to that of the parent complex 1. Detailed spectroscopic analysis on 2-4 leads to the identification of the theoretically expected six normal modes of the Ta(2)N(2) core. On the basis of experimental frequencies and the QCB-NCA procedure, the force constants are determined. Importantly, the N-N force constant decreases from 2.430 mdyn A(-1) in 1 to 1.876 (2), 1.729 (3), and 1.515 mdyn A(-1) (4), in line with the sequence of N-N bond lengths determined crystallographically. DFT calculations on a generic model of the Lewis acid adducts 2-4 reveal that the major donor interaction between the terminal nitrogen atom and the Lewis acid is mediated by a sigma/pi hybrid molecular orbital of N(2), corresponding to a sigma bond. Charge analysis performed for the adducts indicates that the negative charge on the terminal nitrogen atom of the dinitrogen ligand increases with respect to 1. The lengthening of the N-N bond observed for the Lewis adducts is therefore explained by the fact that charge donation from the complex fragment into the pi* orbitals of dinitrogen is increased, while electron density from the N-N bonding orbitals p(sigma) and pi(h) is withdrawn due to the sigma interaction with the Lewis acid.     

Exploring Electronic and Steric Effects on the Insertion and Polymerization Reactivity of Phosphinesulfonato PdII Catalysts

Thirteen different symmetric and asymmetric phosphinesulfonato palladium complexes ([{( ^X 1 ‐Cl)‐μ‐M}_n], M=Na, Li, 1 =^X(P^O)PdMe) were prepared (see Figure 1). The solid‐state structures of the corresponding pyridine or lutidine complexes were determined for ^(MeO)2 1‐py , ^(iPrO)2 1‐lut , ^(MeO,Me2) 1‐lut , ^(MeO)3 1‐lut , ^CF3 1‐lut , and ^Ph 1‐lut . The reactivities of the catalysts ^X 1 , obtained after chloride abstraction with AgBF₄, toward methyl acrylate (MA) were quantified through determination of the rate constants for the first and the consecutive MA insertion and the analysis of β‐H and other decomposition products through NMR spectroscopy. Differences in the homo‐ and copolymerization of ethylene and MA regarding catalyst activity and stability over time, polymer molecular weight, and polar co‐monomer incorporation were investigated. DFT calculations were performed on the main insertion steps for both monomers to rationalize the effect of the ligand substitution patterns on the polymerization behaviors of the complexes. Full analysis of the data revealed that: 1) electron‐deficient catalysts polymerize with higher activity, but fast deactivation is also observed; 2) the double ortho‐substituted catalysts ^(MeO)2 1 and ^(MeO)3 1 allow very high degrees of MA incorporation at low MA concentrations in the copolymerization; and 3) steric shielding leads to a pronounced increase in polymer molecular weight in the copolymerization. The catalyst properties induced by a given P‐aryl (alkyl) moiety were combined effectively in catalysts with two different non‐chelating aryl moieties, such as ^{c HexO/(MeO)2} 1 , which led to copolymers with significantly increased molecular weights compared to the prototypical ^MeO 1 .     

Nitrogen fixation under mild ambient conditions: part I--the initial dissociation/association step at molybdenum triamidoamine complexes

In several recent studies Schrock and collaborators demonstrated for the first time how molecular dinitrogen can be catalytically transformed under mild and ambient conditions to ammonia by a molybdenum triamidoamine complex. In this work, we investigate the geometrical and electronic structures involved in this process of dinitrogen activation with quantum chemical methods. Density functional theory (DFT) has been employed to calculate the coordination energies of ammonia and dinitrogen relevant for the dissociation/association step in which ammonia is substituted by dinitrogen. In the DFT calculations the triamidoamine chelate ligand has been modeled by a systematic hierarchy of increasingly complex substituents at the amide nitrogen atoms. The most complex ligand considered is an experimentally known ligand with an HMT = 3,5-(2,4,6-Me3C6H2)2C6H3 substituent. Several assumptions by Schrock and collaborators on key reaction steps are confirmed by our calculations. Additional information is provided on many species not yet observed experimentally. Particular attention is paid to the role of the charge of the complexes. The investigation demonstrates that dinitrogen coordination is enhanced for the negatively charged metal fragment, that is, coordination is more favorable for the anionic metal fragment than for the neutral species. Coordination of N2 is least favorable for the cationic metal fragment. Furthermore, ammonia abstraction from the cationic complex is energetically unfavorable, while NH3 abstraction is less difficult from the neutral and easily feasible from the anionic low-spin complex.     

Insights into the Mechanism of Reaction of [(C5Me5)2SmII(thf)2] with CO2 and COS by DFT Studies

Reaction mechanisms for the oxidative reactions of CO₂ and COS with [(C₅Me₅)₂Sm] have been investigated by means of DFT methods. The experimental formation of oxalate and dithiocarbonate complexes is explained. Their formation involve the samarium(III) bimetallic complexes [(C₅Me₅)₂Sm‐CO₂‐Sm(C₅Me₅)₂] and [(C₅Me₅)₂Sm‐COS‐Sm(C₅Me₅)₂] as intermediates, respectively, ruling out radical coupling for the formation of the oxalate complex.     

Acid-catalyzed nucleophilic aromatic substitution: experimental and theoretical exploration of a multistep mechanism

The mechanism for the acid-mediated substitution of a phenolic hydroxyl group with a sulfur nucleophile has been investigated by a combination of experimental and theoretical methods. We conclude that the mechanism is distinctively different in nonpolar solvents (i.e., toluene) compared with polar solvents. The cationic mechanism, proposed for the reaction in polar solvents, is not feasible and the reaction instead proceeds through a multistep mechanism in which the acid (pTsOH) mediates the proton shuffling. From DFT calculations, we found a rate-determining transition state with protonation of the hydroxyl group to generate free water and a tight ion pair between a cationic protonated naphthalene species and a tosylate anion. Kinetic experiments support this mechanism and show that, at moderate concentrations, the reaction is first order with respect to 2-naphthol, n-propanethiol, and p-toluenesulfonic acid (pTsOH). Experimentally determined activation parameters are similar to the calculated values (Delta H exp not equal =105+/-9, Delta H calcd not equal =118 kJ mol(-1); Delta G exp not equal =112+/-18, Delta G calcd not equal =142 kJ mol(-1)).     

A theoretical study of surface reduction mechanisms of CeO(2)(111) and (110) by H(2).

Reaction mechanisms for the interactions between CeO(2)(111) and (110) surfaces are investigated using periodic density functional theory (DFT) calculations. Both standard DFT and DFT+U calculations to examine the effect of the localization of Ce 4f states on the redox chemistry of H(2)-CeO(2) interactions are described. For mechanistic studies, molecular and dissociative local minima are initially located by placing an H(2) molecule at various active sites of the CeO(2) surfaces. The binding energies of physisorbed species optimized using the DFT and DFT+U methods are very weak. The dissociative adsorption reactions producing hydroxylated surfaces are all exothermic; exothermicities at the DFT level range from 4.1 kcal mol(-1) for the (111) to 26.5 kcal mol(-1) for the (110) surface, while those at the DFT+U level are between 65.0 kcal mol(-1) for the (111) and 81.8 kcal mol(-1) for the (110) surface. Predicted vibrational frequencies of adsorbed OH and H(2)O species on the surfaces are in line with available experimental and theoretical results. Potential energy profiles are constructed by connecting molecularly adsorbed and dissociatively adsorbed intermediates on each CeO(2) surface with tight transition states using the nudged elastic band (NEB) method. It is found that the U correction method plays a significant role in energetics, especially for the intermediates of the exit channels and products that are partially reduced. The surface reduction reaction on CeO(2)(110) is energetically much more favorable. Accordingly, oxygen vacancies are more easily formed on the (110) surface than on the (111) surface.     

Oxygenolysis of flavonoid compounds: DFT description of the mechanism for quercetin.

Flavonoids are naturally occurring phenol derivatives present in substantial amounts in a large variety of plants, fruits and vegetables daily eaten by humans. Most of these compounds exhibit several interesting biological activities, such as antiradical and antioxidant actions. Indeed, by complexation with specific enzymes, flavonoids are notably liable to metabolize molecular dioxygen. On the basis of experimental results describing oxygenolysis of the flavonoid quercetin, activated by the enzyme quercetin 2,3-dioxygenase (2,3-QD),ur attention has focused on the role of metal center in the activation of the substrate quercetin. Thus, in the present study, by means of DFT calculations at the B3LYP/ 6-31(+)G* level on model molecular systems, we describe different mechanisms for dioxygen metabolization by quercetin. Stationary points are described, and energetic and structural analyses along the reaction paths are reported. Our calculations show that the copper cation must act as an oxidant towards the substrate and that the reaction proceeds through a 1,3-cycloaddition.