Mechanisms of difluoroethylene ozonolysis: a density functional theory study

We present a density functional theory (DFT) study on the mechanisms of gas-phase ozonolysis of three isomers of difluoroethylene, namely, cis-1,2-difluoroethylene, trans-1,2-difluoroethylene, and 1,1-difluoroethylene. MPW1K/cc-pVDZ and BHandHLYP/cc-pVDZ methods are employed to optimize the geometries of stationary points as well as the points on the minimum energy path (MEP). The energies of all the points were further refined at the QCISD(T)/cc-pVDZ and QCISD(T)/6-31+G(df,p) levels of theory with zero-point energy (ZPE) corrections. The ozone-cis-1,2-difluoroethylene reaction is predicted to be slower than the ozone-trans-1,2-difluoroethylene reaction. The enhanced reactivity of trans-1,2-difluoroethylene relative to the cis isomer is similar to the reactions of ozone with cis- and trans-dichloroethylene. The ozone-1,1-difluoroethylene reaction is predicted to be slower than the ozone-trans-1,2-difluoroethylene reaction. These results are in agreement with experimental studies. The calculated mechanisms indicate that in ozone-difluoroethylene reactions the yields of OH might be trivial, which is different from the reactions of ozone with unsaturated hydrocarbons.     

烯丙醇与臭氧反应机理的密度泛函理论研究

应用密度泛函理论的MPW1K,BHandHLYP和MPWB1K方法,结合6-31+G(d,p)基组优化了烯丙醇与臭氧反应势能面上各驻点的几何构型,通过同一水平的振动频率分析确认了中间体和过渡态.反应路径上的驻点都在HL理论水平下进行单点能量校正,并进行了MPW1K/6-31+G(d,p)水平下的零点振动能校正(ZPE).对反应机理的详尽分析表明臭氧抽取烯丙醇羟基基团中H的通道的反应势垒比臭氧加合烯丙醇双键基团通道的反应势垒高,臭氧与烯丙醇双键加合生成臭氧化物为最可几反应路径.在加合反应历程中,氢迁移通道需经过氢迁移和离解等复杂过程,最终要产生少量的OH自由基,与烃烯类臭氧化反应产生大量OH自由基的结果相反.     

Methane activation on Pt and Pt4: a density functional theory study

The activation mechanisms of a methane molecule on a Pt atom (CH4-Pt) and on a Pt tetramer (CH4-Pt4) were investigated using density functional theory (B3LYP and PW91) calculations. The results from these two functionals are different mostly in predicting the reaction barrier, in particular for the CH4-Pt system. A new lower energy pathway was identified for the CH4 dehydrogenation on a Pt atom. In the new pathway, the PtCH2 + H2 products were formed via a transition state, in which the Pt atom forms a complex with carbene and both dissociated hydrogen atoms. We report here the first theoretical study of methane activation on a Pt4 cluster. Among the five single steps toward dehydrogenation, our results show that the rate-limiting step is the third step, that is, breaking the second C-H bond, which requires overcoming an energy barrier of 28 kcal/mol. On the other hand, the cleavage of the first C-H bond, that is, the first reaction step, requires overcoming an energy barrier of 4 kcal/mol.     

Catalyzing aldehyde hydrosilylation with a molybdenum(VI) complex: a density functional theory study

[MoCl(2)O(2)] catalyzes the hydrosilylation reaction of aldehydes and ketones, as well as the reduction of other related groups, in apparent contrast to its known behavior as an oxidation catalyst. In this work, the mechanism of this reaction is studied by means of density functional theory calculations using the B3LYP functional complemented by experimental data. We found that the most favorable pathway to the first step, the Si--H activation, is a [2+2] addition to the Mo=O bond, in agreement with previous and related work. The stable intermediate that results is a distorted-square-pyramidal hydride complex. In the following step, the aldehyde approaches this species and coordinates weakly through the oxygen atom. Two alternative pathways can be envisaged: the classical reduction, in which a hydrogen atom migrates to the carbon atom to form an alkoxide, which then proceeds to generate the final silyl ether, or a concerted mechanism involving migration of a hydrogen atom to a carbon atom and of a silyl group to an oxygen atom to generate the silyl ether weakly bound to the molybdenum atom. In this Mo(VI) system, the gas-phase free energies of activation for both approaches are very similar, but if solvent effects are taken into account and HSiMe(3) is used as a source of silicon, the classical mechanism is favored. Several unexpected results led us to search for still another route, namely a radical path. The energy involved in this and the classical pathway are similar, which suggests that hydrosilylation of aldehydes and ketones catalyzed by [MoCl(2)O(2)] in acetonitrile may follow a radical pathway, in agreement with experimental results.     

C-H bond activation over metal oxides: a new insight into the dissociation kinetics from density functional theory

The C-H activation on metal oxides is a fundamental process in chemistry. In this paper, we report a density functional theory study on the process of the C-H activation of CH(4) on Pd(111), Pt(111), Ru(0001), Tc(0001), Cu(111), PdO(001), PdO(110), and PdO(100). A linear relationship between the C-H activation barrier and the chemisorption in the dissociation final state on the metal surfaces is obtained, which is consistent with the work in the literature. However, the relationship is poor on the metal oxide surfaces. Instead, a strong linear correlation between the barrier and the lattice O-H bond strength is found on the oxides. The new linear relationship is analyzed and the physical origin is identified.     

A systematic study of CO oxidation on metals and metal oxides: density functional theory calculations

CO oxidation on Ru(0001), Rh(111), Pd(111), Os(0001), Ir(111), Pt(111), and their corresponding metal oxides is studied using density functional theory. It is found that (i) the reactivity of metal oxide is generally higher than that of the corresponding metal, and (ii) on both metals and metal oxides, the higher the chemisorption energy is in the initial state, the larger the reaction barrier. The barriers are further analyzed by decomposing them into electronic and geometric effects, and the higher reactivity of metal oxides is attributed mainly to the surface geometric effect. Moreover, the electronic effect on both metals and metal oxides follows the same pattern: the shorter the OC-O bond distance in the TS, the higher the barrier.     

Crossed-beam study of Co+(3F4)+propane: experiment and density functional theory

A pulsed beam of Co+(3F4) crosses a pulsed beam of C3H8 or C3D8 gas under single collision conditions at collision energies of 0.01 eV and 0.21 eV. After a variable time delay t(ext) = 1-8 micros a fast high voltage pulse extracts product ions into a field-free flight tube for mass analysis. Consistent with earlier work, we observe prompt CoC3H6+ +H2 elimination products in 3:1 excess over CoC2H4+ +CH4 products at 0.21 eV on a 2-10 micros time scale. Long-lived CoC3H8+ complexes fragment predominantly back to Co+ +C3H8 reactants and to H2 elimination products on a 6-24 micros time scale. Density functional theory (B3LYP) calculations provide energies, geometries, and harmonic vibrational frequencies at key stationary points for use in a statistical rate model of the reaction. By adjusting two key multicenter transition state (MCTS) energies downward by 4-7 kcal mol(-1), we obtain good agreement with our decay time results and with the cross section versus collision energy of Armentrout and co-workers from 0.1-1.0 eV. B3LYP theory succeeds in finding relative energies of the MCTSs leading to CH4 and H2 in the proper order to explain the different product branching ratio for Co+ (which favors H2 over CH4) compared with its nearest neighbors Fe+ and Ni+ (which favor CH4 over H2).     

Dimerizations of nitrile oxides to furoxans are stepwise via dinitrosoalkene diradicals: a density functional theory study

Density functional theory calculations at the B3LYP/6-31G* level on the dimerization reactions of acetonitrile oxide and para-chlorobenzonitrile oxide to form furoxans indicate that these processes are stepwise involving dinitrosoalkene intermediates that have considerable diradical character. The rate-determining steps for these two reactions correspond to C-C bond formation. The retardation of dimerization in aromatic nitrile oxides arises from the interruption of conjugation between the nitrile oxide and aryl groups in the C-C bond formation step. The present study also suggests that the isomerization of single-ring furoxans occurs via a diradical intermediate mechanism.     

Substituent effects on benzdiyne: a density functional theory study

Density functional theory (DFT) studies were performed to investigate the effect of substituents on the properties of benzdiyne derivatives. Twelve substituted benzdiynes-C(6)X(2), where X = F, Cl, Br, Me, CF(3), CN, OH, NO(2), NH(2), OMe, NMe(2), and Ph-were considered along with the unsubstituted 1,4-benzdiyne. The structures, vibrational frequencies, and IR intensities of these benzdiynes were studied with a popular three-parameter hybrid density functional (B3LYP) combined with the split-valence 6-31G(d) basis set and Dunning's correlation-consistent polarized triple-zeta (cc-pVTZ) basis set. The relative stabilities of the substituted benzdiynes were studied with the help of reaction energies of isodesmic reactions, which showed that the electron-withdrawing groups destabilized the benzdiynes more than they did the corresponding benzenes, whereas the electron-donating groups stabilized the benzdiynes more than they did their benzene counterparts. Correlation analyses revealed that field/inductive effects played a more important role than did resonance effects. The changes in atomic charges and spin populations due to the substituents were also studied. The asymmetric nu(Ctbd1;C) stretching modes obtained were close to the 1500-cm(-)(1) mark. Reinvestigation of the experimental results supported these results; a weak IR band at 1486 cm(-)(1) was assigned to this asymmetric stretching mode in C(6)(CF(3))(2) F. Some other benzdiynes also had large IR intensity values for their asymmetric nu(Ctbd1;C) vibrational modes due to the coupling with other vibrational modes. Heats of formation for the substituted benzdiynes were obtained from the reaction energies calculated at the B3LYP/cc-pVTZ level of theory.     

Nitrogen activation via three-coordinate molybdenum complexes: comparison of density functional theory performance with wave function based methods

Obtaining an accurate theoretical model for the activation of dinitrogen by three-coordinate molybdenum amide complexes (e.g. Mo(NH2)3) is difficult due to the interaction of various high- and low-spin open-shell complexes along the reaction coordinate which must be treated with comparable levels of accuracy in order to obtain reasonable potential energy surfaces. Density functional theory with present-day functionals is a popular choice in this situation; however, the dinitrogen activation reaction energetics vary substantially with the choice of functional. An assessment of the reaction using specialized wave function based methods indicates that although current density functionals in general agree qualitatively on the mechanistic details of the reaction, a variety of high-level electron correlation methods (including CCSD(T), OD(T), CCSD(2), KS-CCSD(T), and spin-flip CCSD) provide a consistent but slightly different representation of the system.     

Mechanisms of Staudinger reactions within density functional theory

The Staudinger reactions of substituted phosphanes and azides have been investigated by using density functional theory. Four different initial reaction mechanisms have been found. All systems studied go through a cis-transition state rather than a trans-transition state or a one-step transition state. The one-step pathway of the phosphorus atom attacking the substituted nitrogen atom is always unfavorable energetically. Depending on the substituents on the azide and the phosphane, the reaction mechanism with the lowest initial reaction barrier can be classified into three categories: (1). like the parent reaction, PH(3) + N(3)H, the reaction goes through a cis-transition state, approaches a cis-intermediate, overcomes a PN-bond-shifting transition state, reaches a four-membered ring intermediate, dissociates N(2) by overcoming a small barrier, and results in the final products: N(2) and a phosphazene; (2). once reaching the cis-intermediate, the reaction goes through the N(2)-eliminating transition state and produces the final products; (3). the reaction has a concerted initial cis-transition state, in which the phosphorus atom attacks the first and the third nitrogen atoms of the azide simultaneously and reaches an intermediate, and then the reaction goes through similar steps of the first reaction mechanism. In contrast to the previous predictions on the relative stability of the unsubstituted cis-configured phosphazide intermediate and the unsubstituted trans-configured phosphazide intermediate, the total energy of the substituted trans-configured phosphazide intermediate is close to that of the substituted cis-configured phosphazide intermediate. The preference of the initial cis-transition state reaction pathway is thoroughly discussed. The relative stability of the cis- and the trans-intermediates is explored and analyzed with the aid of molecular orbitals. The effects of substituents and solvent effects on the reaction mechanisms of the Staudinger reactions are discussed in detail.     

Woodward-Hoffmann rules in density functional theory: initial hardness response

The Woodward-Hoffmann rules for pericyclic reactions, a fundamental set of reactivity rules in organic chemistry, are formulated in the language of conceptual density functional theory (DFT). DFT provides an elegant framework to introduce chemical concepts and principles in a quantitative manner, partly because it is formulated without explicit reference to a wave function, on whose symmetry properties the Woodward-Hoffmann [J. Am. Chem. Soc. 87, 395 (1965)] rules are based. We have studied the initial chemical hardness response using a model reaction profile for two prototypical pericyclic reactions, the Diels-Alder cycloaddition of 1,3-butadiene to ethylene and the addition of ethylene to ethylene, both in the singlet ground state and in the first triplet excited state. For the reaction that is thermally allowed but photochemically forbidden, the initial hardness response is positive along the singlet reaction profile. (By contrast, for the triplet reaction profile, a negative hardness response is observed.) For the photochemically allowed, thermally forbidden reaction, the behavior of the chemical hardness along the initial stages of the singlet and triplet reaction profiles is reversed. This constitutes a first step in showing that chemical concepts from DFT can be invoked to explain results that would otherwise require invoking the phase of the wave function.     

Adsorption of atoms on cu surfaces: a density functional theory study

The chemisorption of atoms (H, N, S, O, and C) on Cu surfaces has been systematically studied by the density functional theory generalized gradient approximation method with the slab model. Our calculated results indicate that the orders of the adsorption energy are H < N < S < O < C on Cu(111) and H < N < O < S < C on Cu(110) and Cu(100). Furthermore, the adsorption energies of the given atoms on Cu(100) are larger than those on Cu(111) and Cu(110). The preferred adsorption sites are a 3-fold hollow site on Cu(111) and a 4-fold hollow site on Cu(100), but the preferred adsorption sites on Cu(110) are different for different adatoms. The energy, as well as the geometry, is in good agreement with the experimental and other theoretical data. In addition, this study focuses on the electronic and geometric properties of the metal-atom (M-A) bond to explain the difference in adsorption energies among adatoms. A detailed investigation of the density of states curves explains the nature of the most stable site. Finally, we test the effect of the coverage and find that the surface coverage has no influence on the preferred adsorption sites of the given adatoms on Cu(110) with the exception of hydrogen and oxygen, but has much influence on the value of the adsorption energy.     

A density functional theory study on pelargonidin

A complete conformational analysis on the isolated and polarizable continuum model (PCM) modeled aqueous solution cation, quinonoidal, and anion forms of pelargonidin, comprising the diverse tautomers of the latter forms, was carried out at the B3LYP/6-31++G(d,p) level. The results indicate that the most stable conformer of cationic and quinonoidal forms of pelargonidin are completely planar in the gas phase, whereas that of the anionic form is not planar. In contrast, PCM calculations show that the plane of the B ring is slightly rotated with regard to the AC bicycle in the most stable conformer of the cation and quinonoidal form. The most stable conformers of the cation, both in gas phase and aqueous solution, display anti and syn orientations for, respectively, C2-C3-O-H and C6-C5-O-H dihedral angles, whereas syn and anti orientation of hydroxyls at 7 and 4' positions are nearly isoenergetic. The most stable tautomer of quinonoidal pelargonidin is obtained by deprotonating hydroxyl at C5 in gas phase but at C7 according to PCM. Also, the most stable tautomer of the anion is different in gas phase (hydrogens are abstracted from hydroxyls at C5 and C4') and PCM simulation (C3 and C5). Tautomeric equilibria affect substantially the geometries of the AC-B backbone providing bond length variations that basically agree with the predictions of the resonance model. Most of the conformers obtained display an intramolecular hydrogen bond between O3 and H6'. Nevertheless, this interaction is not present in the most stable anions. Ionization potentials and O-H bond dissociation energies computed for the most stable conformers of cation, quinonoidal, and anion forms are consistent with an important antioxidant activity.     

Matrix infrared spectra and density functional theory calculations of molybdenum hydrides

Laser-ablated Mo atoms react with H2 upon condensation in excess argon, neon, and hydrogen. The molybdenum hydrides MoH, MoH2, MoH4, and MoH6 are identified by isotopic substitution (H2, D2, HD, H2 + D2) and by comparison with vibrational frequencies calculated by density functional theory. The MoH2 molecule is bent, MoH4 is tetrahedral, and MoH6 appears to have the distorted trigonal prism structure.     

Density functional theory study of CO adsorption on molybdenum sulfide

CO adsorption on four MoSx (stoichiometric and nonstoichiometric) clusters has been investigated by using density functional method. It is found that CO prefers adsorption on the coordinatively unsaturated (1010) surface. The adsorption energy of high coverage shows the additivity as compared with that of one CO adsorption, and there is no significant repulsive interaction between the end-on adsorbed CO probes. The computed CO stretching frequencies (2000-2080 cm(-1)) agree perfectly with the experimental data (a broad band centered at 2070 cm(-1) with a tail extent to 2000 cm(-1)). No bridged CO adsorption is favored energetically under high CO concentration, and this might explain the catalytic ability of MoSx for C1 products instead of higher hydrocarbons and alcohols.     

C-N coupling on transition metal surfaces: a density functional theory study

We have investigated the formation of C-N bonds from individual atoms and single hydrogenated moieties on a series of transition metals. These reactions play a role in HCN formation at high oxygen coverage, also known as Andrussow oxidation, and they are fundamental to understand the ability of other materials to form part of alloys where Pt is the major component. Dehydrogenations take place quite easily under these high oxygen conditions and thus, the C+N, HC+N, and N+CH recombinations to form HCN or its isomer CNH might represent the rate-limiting steps for the reaction. For all the metals in the present study we have found that the activation energy for the reactions between H(x)C and NH(y) (x,y = 0,1) involved in C-N formation follow a linear relationship with the adsorption energy of the N atom. This is due to the common nature of all these transition states, where N-containing fragments get activated from three-fold hollow sites to bridge positions. The slopes of the linear dependence, though, depend on the valence of the N fragment, i.e., smaller slopes are found for NH moieties with respect to N ones.     

Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

Atmospheric oxidation mechanism of p-xylene: a density functional theory study

We report an investigation of the mechanistic features of OH-initiated oxidation reactions of p-xylene using density function theory (DFT). Reaction energies for the formation of the aromatic intermediate radicals have been obtained to determine their relative stability and reversibility, and their activation barriers have been analyzed to assess the energetically favorable pathways to propagate the p-xylene oxidation. OH addition is predicted to occur dominantly at the ortho position, with branching ratios of 0.8 and 0.2 for ortho and ipso additions, respectively, and the calculated overall rate constant is in agreement with available experimental studies. Under atmospheric conditions, the p-xylene peroxy radicals arising from initial OH and subsequent O(2) additions to the ring are shown to cyclize to form bicyclic radicals, rather than to react with NO to lead to ozone formation. With relatively low barriers, isomerization of the p-xylene bicyclic radicals to more stable epoxide radicals likely occurs, competing with O(2) addition to form bicyclic peroxy radicals. The study provides thermochemical and kinetic data for assessment of the photochemical production potential of ozone and formation of toxic products and secondary organic aerosol from p-xylene oxidation.     

Laser ionization of solid RDX: a density functional theory study

Structural Chemistry - Structural changes induced via ionization in an RDX lattice have been studied by using optimized [(RDX)2]0 conformers comprising eight combinations of four RDX isomers using...