Theoretical study of long oligothiophene dications: bipolaron vs polaron pair vs triplet state

A series of oligothiophene dications (from the sexithiophene dication to the 50-mer oligothiophene dication, nT2+, n = 6-50) were studied. Density functional theory (DFT) at the B3LYP/6-31G(d) level and, in some cases, also at BLYP/6-31Gd, was applied to study the singlet and triplet states of the whole series. We found that the singlet state is the ground state for all oligothiophene dications up to the 20-mer, and that the singlet and triplet states are degenerate for longer oligomers. Thus, the triplet state is never a pure ground state for these dications. We found that, for short oligothiophenes dication (e.g., 6T2+), the bipolaron state is the more important state, with only a small contribution made by the polaron pair state. For medium size oligothiophene dications (e.g., 14T2+), both the bipolaron state and the polaron-pair state contribute to the electronic structure. Finally, in long oligothiophene dications, such as 30T2+ and 50T2+, the contribution from the polaron pair state becomes dominant, and these molecules can be considered as consisting of two independent cation radicals or a polaron pair. Results from isodesmic reactions show that the stability of oligothiophene cation radicals over dications is inversely proportional to chain length. Small oligothiophene dications (n = 6-12) were studied at the CASSCF(m,m)/6-31G(d) (m = 4, 6, and 10) level. The major conclusions of this paper regarding the relative energy of the singlet state versus the triplet state and regarding the relative stability of the bipolaron versus the polaron pair were also supported by CASSCF calculations.     

Role of singlet diradicals in reactions of 2-carbenabicyclo[3.2.1]octa-3,6-diene

The generation of 2-carbenabicyclo[3.2.1]octa-3,6-diene (1) results in the formation of C(8)H(8) hydrocarbons endo-6-ethynylbicyclo[3.1.0]hex-2-ene (4), semibullvalene (5), and 5-ethynyl-1,3-cyclohexadiene (6), and C(8)H(10) hydrocarbons bicyclo[3.2.1]octa-2,6-diene (7), tricyclo[3.2.1.0(4,6)]oct-2-ene (8), and tetracyclo[3.3.0.0(2,8)0(4,6)]octane (9). Focus is placed on three mechanistic pathways for the formation of the C(8)H(10) hydrocarbon fraction: (a) abstraction of hydrogen by triplet carbene 1T to produce an equilibrating set of monoradicals, (b) interconversion of triplet carbene 1T into tricyclic triplet diradical 19T and tetracyclic triplet diradical 20T, and (c) interconversion of singlet 1S with analogous singlet diradical 19S and 20S. Ab initio calculations at the (U)B3LYP/6-311+G(3df,2p)//(U)B3LYP/6-31G(d,p) and broken spin symmetry UBS B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) levels rule out choices (a) and (b) and are consistent with the singlet diradical process.     

Thermochemistry of 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide. Evaluation of the mean dissociation enthalpy of the (N-O) bond

The standard enthalpy of formation of the 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide compound in the gas-phase was derived from the enthalpies of combustion of the crystalline solid measured by static bomb combustion calorimetry and its enthalpy of sublimation determined by Knudsen mass-loss effusion at T= 298.15 K. This value is (383.8 +/- 5.4) kJ mol(-1) and was subsequently combined with the experimental gas-phase enthalpy of formation of atomic oxygen and with the computed gas-phase enthalpy of formation of 2-amino-3-quinoxalinecarbonitrile, (382.0 +/- 6.3) kJ mol(-1), in order to estimate the mean (N-O) bond dissociation enthalpy in the gas-phase of 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide. The result obtained is (248.3 +/- 8.3) kJ mol(-1), which is in excellent agreement with the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d) computed value.     

Homoaromaticity in carbene intermediates

Analyses of isodesmic reactions for singlet 7-carbenanorbornene (4S), 8-carbena-endo-tricyclo[3.2.1.0(2,4)]octane (5S), 3-carbenabicyclo[3.1.0]hexane (3S), 2-carbenanorbornene 6S, and 2-carbenabicyclooctadiene 7S at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level provide stabilization energies of 13.83, 13. 50, 3.00, -2.22, and -3.01 kcal/mol, respectively. The C7 carbene in 4S and the C8 carbene center in 5S are strongly bent toward the double bond and cyclopropane ring, respectively, in contrast to their related triplets, 4T and 5T and parent hydrocarbons. The geometric change for 3S compared to 3T or to parent bicyclo[3.1. 0]hexane is minimal. Comparison of the stability of 6S with 2-carbenanorbornane and the geometry of 6S with that of 6T and also with the singlet and triplet 2-carbenanorbornane suggests very modest bridging. The stabilization energy and geometry of 2-carbenabicyclooctadiene 7 resemble antihomoaromatic bicyclooctadienyl cation 9 rather than the related homoaromatic bicyclooctadienyl anion 8. The diamagnetic susceptibility exhaltations (Deltachi) of 3, 4, 5, and 7, calculated at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G level, are -0.7, 22.7, 26.0, and -10.3 cgs-ppm, respectively. The singlet-triplet energy differences, DeltaE(TS), for carbenacyclohexane, carbenacyclopentane, 3, 4, 5, 6, 7, and 2-carbenabicyclo[3.2.1]oct-3-ene at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level are 3.4, 10.3, 8.8, 27.1, 25.9, 12.7, -4.0, and -0.9 kcal/mol, respectively.     

Thermal rearrangements of 2-vinylcyclopropylidene to cyclopentadiene and vinylallene: a theoretical investigation

In an attempt to clarify the favored rearrangement reaction of vinylcyclopropylidenes, the prototype thermal rearrangements of singlet 2-vinylcyclopropylidene (1) leading to 1,3cyclopentadiene (2) and 1,2,4-pentatriene (vinylallene) (3) were investigated by means of ab initio quantum-mechanical electronic-structure calculations. The B3LYP functional with the 6-31G(d) basis set was employed for geometry optimization of the equilibrium and transition-state structures relevant to the two reaction pathways and for computing their harmonic vibrational frequencies. Final energies were evaluated by single-point calculations at the CCSD(T) level of theory with the 6-311 + G(3df,2p) basis set. The rearrangement of s-cis 1 to 2 is found to occur by a three-step pathway. The first step involves the formation of a nonclassical carbene (5), which is an internal pi complex between the pi molecular orbital of the double bond and the empty p atomic orbital of the carbene carbon. In the second step, the nonplanar five-membered ring geometry of 5 flattens to reach the planar structure of 3-cyclopentenylidene (4). The last step is the 1,2-migration of a alpha-hydrogen atom to the carbene center in 4. The rate-determining step for the rearrangement of s-cis 1 to 2 is the formation of 5, with a predicted global deltaG++(220 K) of only 0.6 kcalmol(-1). The rearrangement of s-trans 1 to 2 requires an initial conversion of s-trans 1 to the s-cis conformer, with a predicted deltaG++(220 K) of 1.8 kcalmol(-1). The transition structure for the ring-opening of s-trans 1 into s-trans 3 (deltaG++(220 K)=4.7 kcalmol(-1)) is more energetic than that for the ring-opening of s-cis 1 into s-cis 3 (deltaG++(220 K)=2.5 kcalmol(-2)) due to larger repulsive nonbonded H...H interactions in the former transition structure. On the basis of these results, it is suggested that if the reaction of 1,1-dibromo-2-vinylcyclopropane with methyllithium at -78 degrees C leads to the initial formation of carbene 1, then the reaction should yield 2 as the main product together with small amounts of 3. This theoretical prediction nicely agrees with experimental findings.     

烯丙基自由基(·C3H5)与一氧化氮(NO)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE水平上对反应·CHCHCH3+NO进行了计算, 并建立了其单重态的反应势能面. 在该反应中, 分别找到生成P1(CH3CHO+HCN), P2(CH3CHO+HNC), P3(CH3CN+HCHO), P4(CH3CCH+HNO)的4条产物通道, 其中·CHCHCH3和NO中的氮原子直接连接形成m1(trans-CH3CHCHNO), m1经过顺反异构形成m2(cis-CH3CHCHNO), m2再经过CCNO四元环合, 然后发生环解离, 最后生成产物P1(CH3CHO+HCN)是最可行的产物通道, 其余三条通道为次要产物通道. 该体系中生成P1的反应路径与同类体系·C2H3+NO的主要反应路径相类似, 两者的差别是前者为动力学可行的反应, 而后者为动力学不可行反应, 这使得·CHCHCH3+NO反应比·C2H3+NO反应更具有实际意义.     

The origin of the two-electron/four-centers C--C bond in pi-TCNE(2)2- dimers: electrostatic or dispersion?

The structure and stability of the pi-TCNE(2)2- dimers in K2TCNE2 aggregates is revisited trying to find if the origin of their two-electron/four-centers C--C bond are the electrostatic K+-TCNE- interactions or the dispersion interactions between the anions. The study is done at the HF, B3LYP, CASSCF (2,2), and MCQDPT/CASSCF (2,2) levels using the 6-31+G(d) basis set. Our results show that the only minima of this aggregate that preserves the pi-TCNE(2)2- structure has the two K+ atoms placed in equatorial positions in between the two TCNE- planes. When the K+ atoms are placed along the D2h axis of the anions the structure is not a minimum. The main energetic component responsible for the stability of these aggregates comes from the cation-anion interactions. However, a proper accounting of the dispersion component (as done in the MCQDPT/CASSCF (2,2) calculations) is needed to make the closed-shell singlet more stable than the open-shell singlet. Thus, the bond results from the combination of the electrostatic and dispersion components, being the first the dominant one. The optimum geometry of the closed-shell singlet is very similar to the experimental one found in crystals.     

Theoretical study of the reactions of CF3OCHF2 with the hydroxyl radical and the chlorine atom.

Theoretical investigations are carried out on the reaction mechanism of the reactions of CF3OCHF2 (HFOC-125) with the OH radials and Cl atoms, as well as the heats of formation of CF3OCHF2 and CF3OCF2. The electronic structure information on the potential energy surface for each reaction is obtained at the B3LYP/6-311G(d,p) level, and energetic information is further refined by G3(MP2) theory. The direct dynamics calculation of the hydrogen abstraction reactions are also performed at the G3(MP2)//B3LYP/ 6-311G(d,p) level. The classical energy profile is corrected by interpolated single-point energies (ISPE) approach, incorporating the small-curvature tunnelling effect (SCT) calculated by the variational transition-state theory (VTST). The rate constants are in good agreement with the experimental data and are found to be k1 = 4.95 x 10(-30)T(5.40)exp(-347/T) and k2 = 1.86 x 10(-23)T(3.43)exp(-1579/T) cm3 molecule(-1)s(-1) over the temperature range 200-2000 K. The rate constants at 298 K for these two reactions are 3.38 x 10(-16) and 2.80 x 10(-17) cm3 molecule(-1)s(-1), respectively. Using group-balanced isodesmic reactions as working chemical reactions, the standard enthalpies of formation for CF3OCHF2 and CF3OCF2 are -312.3 +/- 1.0 and -257.3 +/- 1.0 kcalmol(-1), respectively, evaluated by G3(MP2) theory based on the B3LYP/6-311G(d,p) geometries. The theoretical studies provide rate constants of the title reactions and the enthalpies of the species, which are important parameters in determining the atmospheric lifetime and the feasible pathways for the loss of HFOC-125.     

Theoretical study of the energetics of the reactions of triplet dioxygen with hydroquinone, semiquinone, and their protonated forms: relation to the mechanism of superoxide generation in the respiratory chain

One-electron reduction of the dioxygen molecule by the reduced form of mitochondrial ubiquinones (Q) of the NADH dehydrogenase (complex I) and mitochondrial cytochrome bc1 (complex III) is believed to be the main source of the superoxide anion radical O2*- and the hydroperoxide radical OOH*. In this work, we modeled the energetics of four possible reactions of the triplet ((3)Sigma(g)) dioxygen-molecule reduction by fully reduced and protonated ubiquinone (QH2; reaction 1), its deprotonated form (QH-; reaction 2), the semiquinone radical (QH*; reaction 3), and the semiquinone anion radical (Q*-; reaction 4), by means of ab initio calculations with the 6-31G(d) and 6-31+G(d) basis set in the restricted open-shell Hartree-Fock (ROHF), unrestricted Hartree-Fock (UHF), and complete active space self-consistent field (CASSCF) with dynamic correlation [at the second-order M?ller-Plesset (MP2) or multiple reference M?ller-Plesset (MRMP), respectively] schemes and the basis set superposition error (BSSE) correction included, as well as semiempirical AM1 and PM3 calculations in the UHF and ROHF schemes. 2-Butene-1,4-dione and p-benzoquinone were selected as model compounds. For the reduced forms of both compounds, reaction 1 turned out to be energetically unfavorable at all levels of theory, this agreeing with the experimentally observed diminished reductive properties of hydroquinone derivatives at low pH. For 2-butene-1,4-dione treated at the most advanced MRMP/CASSCF/6-31+G(d) level, the energies of reactions 1-4 are 4.7, -34.3, -15.0, and -4.1 kcal/mol, respectively. This finding suggests that reactions 2 and 3 are the most likely mechanisms of electron transfer to molecular oxygen in aprotic environments and that proton transfer is involved in this process. Nearly the same energies of reactions 2 and 3 were calculated at the MRMP/CASSCF/6-31+G(d) level for reduced forms of p-benzoquinone. Inclusion of diffuse functions in the basis set and dynamic correlation at the CASSCF level appears essential. Because deprotonated ubiquinol is unlikely to exist in physiological environments, reaction 3 appears to be the most likely mechanism of one-electron reduction of oxygen; however, if oxygen can penetrate cytochrome bc1 as far as the Q(o) center where ubiquinol can be deprotonated, reaction 2 can also come into play. The energies of reactions 2 and 3 calculated at the MRMP/CASSCF/6-31+G(d) level are most closely reproduced in the ab initio and semiempirical UHF PM3 calculations. Additional semiempirical calculations on more realistic models of ubiquinone, 2,3-dimethoxy-6-methyl-p-benzoquinone and 2,3-dimethoxy-5-isoprenyl-6-methyl-p-benzoquinone, gave qualitatively the same relations between the energies of reactions 2 and 3 as those carried out for p-benzoquinone species, thereby suggesting that this method could be used in studying electron-transfer reactions from reduced quinone derivatives to molecular oxygen in more complex systems, such as a model of the Q(o) site of cytochrome bc1, where applying ab initio methods is unfeasible.     

Theoretical study on the mechanism of the 1CHCl + NO reaction

The complex doublet potential energy surface of the CHClNO system, including 31 minimum isomers and 84 transition states, is investigated at the QCISD(T)/6-311G(d, p)//B3LYP/6-31G(d, p) level in order to explore the possible reaction mechanism of the singlet CHCl with NO. Various possible isomerization and dissociation channels are probed. The initial association between 1CHCl and NO at the terminal N-site can almost barrierlessly lead to the chainlike adducts HClCNO a (a1, a2) followed by the direct Cl-extrusion to product P9 Cl + HCNO, which is the most feasible channel. Much less competitively, a (a1, a2) undergoes a ring-closure leading to the cyclic isomer c-C(HCl)NO d followed by a concerted Cl-shift and N-O cleavage of d to form the branched isomers ClNC(H)O f (f1, f2). Eventually, f (f1, f2) may take a direct H-extrusion to produce P7 H + ClNCO or a concerted 1,2-H-shift and Cl-extrusion to form P1 Cl + HNCO. The low-lying products P2 HCl + NCO, P3 Cl + HOCN, P14 HCO + 3NCl, P6 ClO + HCN, and P13 ClNC + OH may have the lowest yields observed. Our calculations show that the product distributions of the title reaction are quite different from those of the analogous 1CHF + NO reaction, yet are similar to those of another analogous 3CH2 + NO reaction. The similarities and discrepancies among the three reactions are discussed in terms of the substitution effect. The present article may assist in future experimental identification of the product distributions for the title reaction and may be helpful for understanding the halogenated carbene chemistry.     

Theoretical Study on Photoisomerization of 2-Methylfuran to 3-Methylfuran

1 INTRODUCTION 2-Methylfuran belongs to the basic heteroaromatic compounds relevant to many fields of modern che- mistry, ranging from the study of natural products and biologically active substances to the develop- ment of building blocks for organic synthesis and conducting polymers[1]. Since the photochemistry ofR-furan was gradually recognized in 1960s[2～7], lots of interest has been aroused. Herein we only study one branch of photoche- mistry of R-furan: the isomerization of 2-methy…     

A theoretical characterization of the photoisomerization channels of 1,2-cyclononadienes on both singlet and triplet potential-energy surfaces

The ground-, (1)(pipi*)-, and (3)(pipi*)-state potential-energy surfaces of 1,2-cyclononadiene and isomeric C(9)H(14) species, as well as 1-methyl-1,2-cyclononadiene and isomeric C(10)H(16) species were all mapped using CASSCF and the 6-31G(d) basis set. Theoretical results were found to be in good agreement with the available experimental observations for both 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene isomerization reactions under singlet and triplet direct or sensitized irradiation. Extremely efficient decay occurs from the first singlet excited state to the ground state through at least three different conical intersections (surface crossings). The first of these crossing points is accessed by a one-bond ring closure. From this conical intersection point (CI-A or CI-C), some possible subsequent ground-state reaction paths have been identified: 1) intramolecular C--H bond insertion to form the bicyclic photoproduct and 2) intramolecular C--H bond insertion to form tricyclic photoproducts. An excited state [1,3]-sigmatropic shift leads to the second conical intersection (CI-B or CI-E), which can give a three-bond cyclononyne species. Besides these, in the singlet photochemical reactions of 1-methyl-1,2-cyclononadiene, excited-state, one allenic C--H bond insertion leads to a third conical intersection (CI-D). Possible ground-state reaction pathways from this structure lead to the formation of a diene photoproduct or to transannular insertion photoproducts. Moreover, in the case of triplet 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene photoisomerization reactions, both chemical reactions will adopt a 1,3-biradical (T(1)/S(0)-1, T(1)/S(0)-2, and T(1)/S(0)-3), which may undergo intersystem crossings leading to the formation of tricyclic or bicyclic photoproducts. The results obtained allow a number of predictions to be made.     

Ab initio study on the oxidation of NCN by O (3P): prediction of the total rate constant and product branching ratios

The reaction of NCN with O is relevant to the formation of prompt NO according to the new mechanism, CH+N2-->cyclic-C(H)NN- -->HNCN-->H+NCN. The reaction has been investigated by ab initio molecular orbital and transition state theory calculations. The mechanisms for formation of possible product channels involved in the singlet and triplet potential energy surfaces have been predicted at the highest level of the modified GAUSSIAN-2 (G2M) method, G2M (CC1). The barrierless association/dissociation processes on the singlet surface were also examined with the third-order Rayleigh-Schr?dinger perturbation (CASPT3) and the multireference configuration interaction methods including Davidson's correction for higher excitations (MRCI+Q) at the CASPT3(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) and MRCI+Q(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) levels. The rate constants for the low-energy channels producing CO+N2, CN+NO, and N(4S)+NCO have been calculated in the temperature range of 200-3000 K. The results show that the formation of CN+NO is dominant and its branching ratio is over 99% in the whole temperature range; no pressure dependence was noted at pressures below 100 atm. The total rate constant can be expressed by: kt=4.23x10(-11) T0.15 exp(17/T) cm3 molecule(-1) s(-1).     

Reaction Mechanism of 3,4‐Dinitrofuroxan Formation from Glyoxime: Dehydrogenation and Cyclization of Oxime

The reaction pathway of the formation of 3,4‐dinitrofuroxan from glyoxime is theoretically investigated under experimental conditions with 25 % nitric acid and dinitrogentetroxide reagents to clarify the mechanism of formation of a furoxan ring by glyoxime. The geometric configurations of minima and transition‐state species are optimized at the (U)B3LYP/6‐311++G** level. The CCSD(T) and CASSCF(10e,8o)/CASSCF(9e,8o) single‐point energy corrections at the same level are performed on top of the optimized geometries. A subsequent dynamic correlation by using NEVPT2/6‐311++G**‐level single‐point energy calculations based on the CASSCF results is also performed to obtain accurate energy values. The formation reaction is analyzed from two processes: glyoxime nitration and 3,4‐dinitroglyoxime (nitration product) oxidative cyclization. Calculation results indicate that the electrophilic substitution of nitronium ions from the protonated HNO₃ and the abstraction of hydrogen ions by HNO₃ molecules are requisites of glyoxime nitration. The formation of a furoxan ring from 3,4‐dinitroglyoxime involves two possible mechanisms: 1) oxydehydrogenation by NO₂ molecules and the subsequent torsion of NO radical groups to form a ring and 2) the alternation of dehydrogenation and cyclization. The intermediates and transition states in both routes exhibit monoradical and diradical characteristics. Singlet and triplet reactions are considered for the diradical species. Results show that the singlet reaction mechanism is more favorable for cyclization than the triplet reaction. The formation of a furoxan ring from oxime is in accordance with the stepwise intermolecular dehydrogenation and intramolecular torsion to the ring.     

Theoretical study of low-lying triplet states of aniline

Multireference complete active space self-consistent-field CASSCF(10,12)/ANO and second-order perturbation theory MS-CASPT2 calculations were performed to determine the vertical low-lying singlet and triplet states of aniline. The sequence of the seven lower lying triplet states is T1(1(3)A'), T2(1(3)A' '), T3(2(3)A'), T4(3(3)A'), T5(2(3)A' '), T6(4(3)A'), and T7(3(3)A' '). The 3(3)A', 4(3)A', and 3(3)A' ' states are assigned as 3s, 3py, and 3pz Rydberg states, respectively, while other states correspond to pi <-- pi excitations. Both the T1 and T2 states are found to be below at the lowest-lying singlet S1 (1(1)A' ') state. Geometry, vibrational modes, and electron distribution of the lowest lying T1 state were determined using UB3LYP calculations. The vertical and adiabatic singlet-triplet energy gaps DeltaE(S0-T1) amount to 3.7 and 3.5 +/- 0.2 eV, respectively. In clear contrast with the S0 state, the triplet aniline is no longer aromatic, and its protonation occurs preferentially at the ring meta-carbon site, with a proton affinity PA = 243 +/- 3 kcal/mol.     

Isomers of cyclo-heptasulfur and their coordination to Li(+): an ab initio molecular orbital study

The potential energy hypersurfaces (PESs) of heptasulfur (S7) and of [LiS7]+ have been investigated by ab initio molecular orbital calculations at the G3X(MP2) level of theory. Besides the chair-like seven-membered ring (1a) as the global minimum structure, eight local minimum structures and one transition state have been located on the PES of S7. The barrier for pseudorotation of 1a is only 5.6 kJ mol(-1). The boat-like S7 ring (1b) is 12.1 kJ mol(-1) less stable than 1a, followed by three isomers of connectivity S6=S and four open-chain isomers. On the basis of multireference calculations at the MRCI(4,4)+Q/6-311G(d) level, the most stable open-chain form of S7 is a triplet of relative energy 133.1 kJ mol(-1). Thus, the reaction energy (deltaE0) for the ring opening of 1a is 133.1 kJ mol(-1), halfway between those of the highly symmetrical rings S6 and S8. Because of their strong multireference characters, the stabilities of the biradicalic singlet chains are significantly overestimated by the single-reference-based G3X(MP2) method. The calculated vibrational spectrum of 1a is in good agreement with experimental data. The various isomers of S7 form stable complexes with Li+ with coordination numbers of 1-4 for the metal atom and binding energies in the range of -93.8 to -165.7 kJ mol(-1). A total of 15 isomeric complexes have been located, with 13 of them containing cyclic ligands. The global minimum structure (2a) is composed of 1a, with the Li+ cation linked to the four negatively charged sulfur atoms (symmetry C(s)). The corresponding complex 2c containing the ligand 1b is by 23.4 kJ mol(-1) less stable than 2a, and a bicyclic crown-shaped LiS7 cation (2e) is by 34.9 kJ mol(-1) less stable than 2a. Even less stable are four complexes with the branched S6=S ligand. SS bond activation by polarization of the valence electrons takes place on coordination of Li+ to cyclo-S7 (1a).     

The aromaticity of pyracylene: an experimental and computational study of the energetics of the hydrogenation of acenaphthylene and pyracylene

In this work, the aromaticity of pyracylene (2) was investigated from an energetic point of view. The standard enthalpy of hydrogenation of acenaphthylene (1) to acenaphthene (3) at 298.15 K was determined to be minus sign(114.5 +/- 4.2) kJ x mol(-1) in toluene solution and minus sign(107.9 +/- 4.2) kJ x mol(-1) in the gas phase, by combining results of combustion and reaction-solution calorimetry. A direct calorimetric measurement of the standard enthalpy of hydrogenation of pyracylene (2) to pyracene (4) in toluene at 298.15 K gave -(249.9 plus minus 4.6) kJ x mol(-1). The corresponding enthalpy of hydrogenation in the gas phase, computed from the Delta(f)H(o)m(cr) and DeltaH(o)m(sub) values obtained in this work for 2 and 4, was -(236.0 +/- 7.0) kJ x mol(-1). Molecular mechanics calculations (MM3) led to Delta(hyd)H(o)m(1,g) = -110.9 kJ x mol(-1) and Delta(hyd)H(o)m(2,g) = -249.3 kJ x mol(-1) at 298.15 K. Density functional theory calculations [B3LYP/6-311+G(3d,2p)//B3LYP/6-31G(d)] provided Delta(hyd)H(o)m(2,g) = -(244.6 +/- 8.9) kJ x mol(-1) at 298.15 K. The results are put in perspective with discussions concerning the "aromaticity" of pyracylene. It is concluded that, on energetic grounds, pyracylene is a borderline case in terms of aromaticity/antiaromaticity character.     

A systematic computational study on the reactions of HO2 with RO2: The HO2 + CH3O2(CD3O2) and HO2 + CH2FO2 reactions

A systematic theoretical study of the reactions of HO2 with RO2 has been carried out. The major concern of the present work is to gain insight into the reaction mechanism and then to explain experimental observations and to predict new product channels for this class of reactions of importance in the atmosphere. In this paper, the reaction mechanisms for two reactions, namely, HO2 + CH3O2 and HO2 + CH2FO2, are reported. Both singlet and triplet potential energy surfaces are investigated. The complexity of the present system makes it impossible to use a single ab initio method to map out all the reaction paths. Various ab initio methods including MP2, CISD, QCISD(T), CCSD(T), CASSCF, and density function theory (B3LYP) have been employed with the basis sets ranging from 6-31G(d) to an extrapolated complete basis set (CBS) limit. It has been established that the CCSD(T)/cc-pVDZ//B3LYP/6-311G(d,p) scheme represents the most feasible method for our systematic study. For the HO2 + CH3O2 reaction, the production of CH3OOH is determined to be the dominant channel. For the HO2 + CH2FO2 reaction, both CH2FOOH and CHFO are major products, whereas the formation of CHFO is dominant in the overall reaction. The computational findings give a fair explanation for the experimental observation of the products.     

Mechanisms and kinetics for the thermal decomposition of 2-azido-N,N-dimethylethanamine (DMAZ)

To gain insight into the mechanisms and kinetics of 2-azido-N,N-dimethylethanamine's (DMAZ's) thermal decomposition postulated reaction paths were simulated with ab initio and density functional theory quantum chemistry models. Four reaction types were modeled: (i) spin-allowed and spin-forbidden paths involving N-N(2) bond fission and nitrene formation, (ii) HN(3) elimination with the formation of (dimethylamino)ethylene, (iii) N-N(2) bond fission with the formation of molecules with three- or four-membered heterocyclic rings, and (iv) simple scission of C-H, C-N, and C-C bonds. The geometries of stationary points of the reactions were obtained with a MPWB1K/6-31+G(d,p) model. To locate and model the geometries of minimum energy intersystem crossing points for triplet nitrene formation and isomerization, unrestricted broken spin symmetry calculations were performed. Employed to model an analogous path for methyl azide's decomposition, this approach was found to yield results similar to those obtained with a CASSCF(10,8)/aug-cc-pVDZ model. Of the four reaction types studied, N-N(2) bond fissions with singlet or triplet nitrene formation were found to have the lowest barriers. Barriers for paths to cyclic products were found to be 2-4 kcal/mol higher. Kinetic rate expressions for individual paths were derived from the quantum chemistry results, and spin-allowed nitrene formation was found to be dominant at all temperatures and pressures examined. The expression 2.69 × 10(9) (s(-1))T(1.405) exp(-39.0 (kcal/mol)/RT), which was derived from QCISD(T)/6-31++G(3df,2p)//MPWB1K/6-31+G(d,p) results, was found to be representative of this reaction's gas-phase rate. Adjusted on the basis of results from self-consistent reaction field models to account for solvation by n-dodecane, the expression became 1.11 × 10(9) (s(-1))T(1.480) exp(-37.6 (kcal/mol)/RT). Utilizing this result and others derived in the study, a model of the decomposition of n-dodecane-solvated DMAZ was constructed, and it generated simulations that well-reproduce previously published measured data for the process.     

Au+(1S, 3D)与N2O(1Σ+)反应机理的理论研究

采用密度泛函B3LYP方法, O和N用6-311+G*基组, Au+用赝势基组(8s7p6d)/[6s5p3d], 研究了Au+(1S, 3D)离子和N2O(1Σ+)分子的反应机理. 报道了在基态单重态和激发三重态势能面上各反应物、中间体和过渡态的构型特征及能量. 结果表明, 两个主反应通道Au+(1S)+ N2O(1Σ+)→1NA-Complex-1→1NA-TS1→1NA-Complex-2→1NA-Crossing→[3OAuNN]+和Au+(1S)+ N2O(1Σ+)→1NB-Complex→1NB-Crossing→[AuNN(1Σ+)]++O(3P)都需经过反应交叉势能面, 出现“系间窜越”. 用内禀坐标单点计算垂直激发态的方法确定了势能面交叉点, 并用含时密度泛函TD-B3LYP方法进一步探讨了自旋翻转机理.