Theoretical study of two states reactivity of methane activation on iron atom and iron dimer

Density functional theory (DFT) calculations have been carried out to explore the catalytic activation of C–H bonds in methane by the iron atom, Fe, and the iron dimer, Fe2. For methane activation on an Fe atom, the calculations suggest that the activation of the first C–H bond is mediated via the triplet excited-state potential energy surface (PES), with initial excitation of Fe to the triplet state being necessary for the reaction to be energetically feasible. Compared with the breaking of the first C–H bond, the cleavage of the second C–H bond is predicted to involve a significantly higher barrier, which could explain experimental observations of the HFeCH3 complex rather than CH2FeH2 in the activation of methane by an Fe atom. For methane activation on an iron dimer, the cleavage of the first C–H bond is quite facile with a barrier only 11.2, 15.8 and 8.4 kcal/mol on the septet state energy surface at the B3LYP/6-311+G(2df,2dp), BPW91/6-311+G(2df,2dp) and M06/B3LYP level, respectively. Cleavage of the second C–H bond from HFe2CH3 involves a barrier calculated respectively as 18.0, 10.7 and 12.4 kcal/mol at the three levels. The results suggest that the elimination of hydrogen from the dihydrogen complex is a rate-determining step. Overall, our results indicate that the iron dimer Fe2 has a stronger catalytic effect on the activation of methane than the iron atom.

Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.

Mechanistic Aspects of Nucleophilic Substitution at Half-Sandwich Metal Complexes

The hydrolysis reactions of organometallic ruthenium(II) piano-stool complexes of the type Ru-II(eta(6)-cymene)(L)Cl](0/+) (1-5, where L = kappa(1)- or kappa(2)-1,1-bis(diphenylphosphino)methane,1,1bis-(diphenylphosphino)methane oxide, kappa(1)-mercaptobenzothiazole) have been studied using density functional theory at the B3LYP level. In addition to considering a syn attack in an associative fashion, where the nucleophile approaches from the same side as the leaving group, we have explored alternative paths such as an anti attack in an associative manner, where the nucleophile attacks from the opposite side of the leaving group. During the anti attack, an intermediate is formed and there is a coordination mode change of the arene ring from eta(6) to eta(2) along with its rotation. When the intermediate goes to the product, the arene ring slips back from eta(2) to eta(6) coordination. This coordinated movement of the arene ring makes the associative anti attack an accessible pathway for the substitution process. Our calculations predict very similar activation barriers for both syn and anti attacks. In the dissociative path, the rate-determining step is the generation of a coordinatively unsaturated 16-electron ruthenium species. This turns out to be viable once solvent effects are included. The large size of the ancillary ligands on Ru makes the dissociative process as favorable as the associative process. Activation energy calculations reveal that although the dissociative path is favorable for kappa(1) complexes, both dissociative and associative processes can have significant contribution to the hydrolysis reaction in kappa(2) complexes. Once activated by hydrolysis, these complexes react with guanine and adenine bases of DNA. The thermodynamic stabilities of complexes formed with the nucleobases are also presented.

Isomeric identification of methylated naphthalenes using gas phase infrared spectroscopy

We report gas phase mid-infrared spectra of 1- and 2- methyl naphthalenes at 0.2 cm(-1) resolution. Assignment of observed bands have been made using scaled quantum mechanical (SQM) calculations where the force fields rather the frequencies are scaled to find a close fit between observed and calculated bands. The structure of the molecules has been optimized using B3LYP level of theory in conjunction with standard 6-311G** basis set to obtain the harmonic frequencies. Using the force constants in Cartesian coordinates from the Gaussian output, scaled force field calculations are carried out using a modified version of the UMAT program in the QCPE package. Potential energy distributions of the normal modes obtained from such calculations helped us assign the observed bands and identify the unique features of the spectra of 1- and 2-MNs which are important for their isomeric identification.

Computational Study of the Migration of Rhenium from One Enantioface of an Olefin to the Other Facilitated by (C-H)center dot center dot center dot Re Interactions

The migration of a metal atom in a metal olefin complex from one pi face of the olefin to the opposite pi face has been rarely documented. Gladysz and co-workers showed that such a movement is indeed possible in monosubstituted chiral Re olefin complexes, resulting in diastereomerization. Interestingly, this isomerization occurred without dissociation, and on the basis of kinetic isotope effects, the involvement of a trans C-H bond was indicated. Either oxidative addition or an agostic interaction of the vinylic C-H(D) bond with the metal could account for the experimentally observed kinetic isotope effect. In this study we compute the free energy of activation for the migration of Re from one enantioface of the olefin to the other through various pathways. On the basis of DFT calculations at the B3LYP level we show that a trans (C-H)center dot center dot center dot Re interaction and trans C-H oxidative addition provide a nondissociative path for the diastereomerization. The trans (C-H)center dot center dot center dot Re interaction path is computed to be more favorable by 2.3 kcal mol(-1) than the oxidative addition path. While direct experimental evidence was not able to discount the migration of the metal through the formation of a eta(2)-arene complex (conducted tour mechanism), computational results at the B3LYP level show that it is energetically more expensive. Surprisingly, a similar analysis carried out at the M06 level computes a lower energy path for the conducted tour mechanism and is not consistent with the experimental isotope effects observed. Metal-(C-H) interactions and oxidative additions of the metal into C-H bonds are closely separated in energy and might contribute to unusual fluxional processes such as this diastereomerization.

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYP

Theoretical investigation of LaCn+ (n=2-8) clusters

LaCn+ (n = 2-8) have been studied using Hartree-Fock (HF) and B3LYP density functional method. The results indicated that at both levels, isomers with C-2v, C-s symmetry for n = 2, and edge insertion isomer for n = 4, 6, 8, as well as edge binding isomer for n = 3, 5, 7 were found as ground states. This is in good agreement with experimental results. The exceptional case is for n = 6 at B3LYP level, in which edge insertion and edge binding isomers were computed to be near isoenergetic. (C) 1997 Elsevier Science B.V.

Interactions of 1-Methylimidazole with UO2(CH3CO2)2 and UO2(NO3)2: Structural, Spectroscopic and Theoretical Evidence for Imidazole Binding to the Uranyl Ion

The first definitive high-resolution single-crystal X-ray structure for the coordination of the 1-methylimidazole (Meimid) ligand to UO2(Ac)2 (Ac = CH3CO2) is reported. The crystal structure evidence is confirmed by IR, Raman, and UV-vis spectroscopic data. Direct participation of the nitrogen atom of the Meimid ligand in binding to the uranium center is confirmed. Structural analysis at the DFT (B3LYP) level of theory showed a conformational difference of the Meimid ligand in the free gas-phase complex versus the solid state due to small energetic differences and crystal packing effects. Energetic analysis at the MP2 level in the gas phase supported stronger Meimid binding over H2O binding to both UO2(Ac)2 and UO2(NO3)2. In addition, self-consistent reaction field COSMO calculations were used to assess the aqueous phase energetics of combination and displacement reactions involving H2O and Meimid ligands to UO2R2 (R = Ac, NO3). For both UO2(NO3)2 and UO2(Ac)2, the displacement of H2O by Meimid was predicted to be energetically favorable, consistent with experimental results that suggest Meimid may bind uranyl at physiological pH. Also, log(Knitrate/KAc) calculations supported experimental evidence that the binding stoichiometry of the Meimid ligand is dependent upon the nature of the reactant uranyl complex. These results clearly demonstrate that imidazole binds to uranyl and suggest that binding of histidine residues to uranyl could occur under normal biological conditions.

The reaction path of CO and Fe2O3 in a chemical-looping combustion system

The reaction mechanism of CO and Fe₂O₃ in a chemical-looping combustion (CLC) was studied based on density functional theory (DFT) at B3LYP level in this paper. The structures of all reactants, intermediate, transition structures and products of this reaction had been optimized and characterized. The reaction path was validated by means of the intrinsic reaction coordinate (IRC) approach. The result showed that the reaction was divided into two steps, the adsorbed CO molecule on Fe ₂O₃ surface formed a medium state with one broken Fe-O bond in step1, and in step2, O atom broken here oxidized a subsequent CO molecule in the fuel reactor. Thus, Fe₂O₃ molecule transport O from air to oxide CO continually in the CLC process. The activation energy and rate coefficients of the two steps were also obtained.

A comparison of the reactivity of germylene and dimethylgermylene with some methylgermanes. Direct kinetic and quantum chemical studies

Time-resolved studies of germylene, GeH2, and dimethygermylene, GeMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to try to obtain rate coefficients for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas-phase. GeH2 + Me2GeH2 was studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 296-553 K. Only slight pressure dependences were found (at 386, 447 and 553 K). RRKM modelling was carried out to fit these pressure dependences. The high pressure rate coefficients gave the Arrhenius parameters: log(A/cm(3) molecule(-1)s(-1)) = -10.99 +/- 0.07 and E-a = -(7.35 +/- 0.48) kJ mol(-1). No reaction could be found between GeMe2 + Me2GeH2 at any temperature up to 549 K, and upper limits of ca. 10(-14) cm(3) molecule(-1)s(-1) were set for the rate coefficients. A rate coefficient of (1.33 +/- 0.04) x 10(-11)cm(3) molecule(-1)s(-1) was also obtained for GeH2 + MeGeH3 at 296 K. No reaction was found between GeMe2 and MeGeH3. Rate coefficient comparisons showed, inter alia, that in the substrate germane Me-for-H substitution increased the magnitudes of rate coefficients significantly, while in the germylene Me-for-H substitution decreased the magnitudes of rate coefficients by at least four orders of magnitude. Quantum chemical calculations (G2(MP2,SVP)// B3LYP level) supported these findings and showed that the lack of reactivity of GeMe2 is caused by a positive energy barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper.

Surprisingly slow reaction of dimethylsilylene with dimethylgermane: time-resolved kinetic studies and related quantum chemical calculations

Time-resolved studies of silylene, SiH2, and dimethylsilylene, SiMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to obtain rate constants for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas phase. SiMe2 + Me2GeH2 was studied at five temperatures in the range 299-555 K. Problems of substrate UV absorption at 193 nm at temperatures above 400 K meant that only three temperatures could be used reliably for rate constant measurement. These rate constants gave the Arrhenius parameters log(A/cm(3) molecule(-1) s(-1)) = -13.25 +/- 0.16 and E-a = -(5.01 +/- 1.01) kJ mol(-1). Only room temperature studies of SiH2 were carried out. These gave values of (4.05 +/- 0.06) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + Me2GeH2 at 295 K) and also (4.41 +/- 0.07) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + MeGeH3 at 296 K). Rate constant comparisons show the surprising result that SiMe2 reacts 12.5 times slower with Me2GeH2 than with Me2SiH2. Quantum chemical calculations (G2(MP2,SVP)//B3LYP level) of the model Si-H and Ge-H insertion processes of SiMe2 with SiH4/MeSiH3 and GeH4/MeGeH3 support these findings and show that the lower reactivity of SiMe2 with Ge-H bonds is caused by a higher secondary barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper. Other, related, comparisons of silylene reactivity are also presented.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

Joint Experimental and Theoretical Analysis of Order Disorder Effects in Cubic BaZrO(3) Assembled Nanoparticles under Decaoctahedral Shape

Periodic first-principles calculations based on density functional theory at the B3LYP level has been carried out to investigate the photoluminescence (PL) emission of BaZrO(3) assembled nanoparticles at room temperature. The defect created in the nanocrystals and their resultant electronic features lead to a diversification of electronic recombination within the BaZrO(3) band gap. Its optical phenomena are discussed in the light of photoluminescence emission at the green-yellow region around 570 nm. The theoretical model for displaced atoms and/or angular changes leads to the breaking of the local symmetry, which is based on the refined structure provided by Rietveld methodology. For each situation a band structure, charge mapping, and density of states were built and analyzed. X-ray diffraction (XRD) patterns, UV-vis measurements, and field emission scanning electron microscopy (FE-SEM) images are essential for a full evaluation of the crystal structure and morphology.

Effect of Different Solvent Ratios (Water/Ethylene Glycol) on the Growth Process of CaMoO(4) Crystals and Their Optical Properties

In this paper, calcium molybdate (CaMoO(4)) crystals (meso- and nanoscale) were synthesized by the coprecipitation method using different solvent volume ratios (water/ethylene glycol). Subsequently, the obtained suspensions were processed in microwave-assisted hydrothermal/solvothermal systems at 140 degrees C for 1 h. These meso- and nanocrystals processed were characterized by X-ray diffraction (X R I)), Fourier transform Raman (FT-Raman), Fourier transform infrared (FT-IR). ultraviolet visible (UV-vis) absorption spectroscopies, held-emission gun scanning electron microscopy (FEG-SEM). transmission electron microscopy (TEM). and photoluminescence (PL) measurements. X RI) patterns and FT-Raman spectra showed that these meso- and nanocrystals have a scheelite-type tetragonal structure without the presence of deleterious phases. FT-IR spectra exhibited a large absorption band situated at around 827 cm(-1), which is associated with the Mo-O anti-symmetric stretching vibrations into the [MoO(4)] clusters. FEG-SEM micrographs indicated that the ethylene glycol concentration in the aqueous solution plays an important role in the morphological evolution of CaMoO(4) crystals. High-resolution TEM micrographs demonstrated that the mesocrystals consist of several aggregated nanoparticles with electron diffraction patterns of monocrystal. In addition, the differences observed in the selected area electron diffraction patterns of CaMoO(4) crystals proved the coexistence of both nano- and mesostructures, First-principles quantum mechanical calculations based on the density functional theory at the B3LYP level were employed in order to understand the band structure find density of states For the CaMoO(4). UV-vis absorption measurements evidenced a variation in optical band gap values (from 3.42 to 3.72 cV) for the distinct morphologies. The blue and green PI. emissions observed in these crystals were ascribed to the intermediary energy levels arising from the distortions on the [MoO(4)] clusters clue to intrinsic defects in the lattice of anisotropic/isotropic crystals.