Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

叠氮化氰的光解离机理

采用多参考态方法,在MRCI+Q//CAS(10,9)/6-311+G(2df)水平上对叠氮化氰(N3CN)的光解离机理进行理论研究.优化得到基态(S0)和低激发态(S1、S2、T1)势能面上的极小点、过渡态、内转换交叉点(IC-S1/S0)和隙间窜跃交叉点(ISC-S1/T1)的结构和能量,构建反应势能面.在MRCI+Q//CAS(10,9)水平上计算N3CN的垂直激发能,并和实验值进行对比.结果表明,在S0、S1、S2和T1态势能面上,N—N键断裂生成N2+NCN是主要解离途径,而C—N键断裂通道是次要通道.实验观测到220 nm处的吸收峰对应分子由S0态到S1态的激发,对应主要光解离产物为NCN[a1△g];而在275 nm处的吸收峰则对应分子被激发到T1态,然后直接生成基态产物NCN[X3Σg-].我们的理论结果与实验测量符合得很好.     

Theoretical mechanistic study on the reaction of CN radical with HNCN

The mechanism for the reaction of the cyanogen radical (CN) with the cyanomidyl radical (HNCN) has been investigated theoretically. The electronic structure information of the singlet and triplet potential energy surfaces (PESs) is obtained at the B3LYP/6-311+G(3df,2p) level, and the single-point energies are refined at the CCSD(T)/6-311+G(3df,2p) level as well as by multilevel MCG3-MPWB method. The calculations show that the C atom of CN additions to middle- and end-N atoms of HNCN are two barrierless association processes leading to the energy-rich intermediates IM1 HN(CN)CN and IM2 HNCNCN, respectively, on the singlet PES. The higher barriers of the subsequent isomerization and dissociation channels from IM1 and IM2 indicate that these two intermediates, which have considerably thermodynamic and kinetic stability, are the dominant product at high pressure. While at low pressure, the most favorable product is P(2) H + NCNCN, which will be formed from both IM1 and IM2 via direct dissociation processes by the H-N bond rupture, and the secondary feasible product is P(4) HCN + (1) NCN, while P(5) HCCN + N(2) and P(6) HCNC + N(2) are the least competitive products. On the triplet PES, P(14) NCNC + HN may be a comparable competitive product at high temperature. In addition, the comparison between the mechanisms of the CN + HNCN and OH + HNCN reactions is made. The present results will enrich our understanding of the chemistry of the HNCN radical in combustion processes and interstellar space.     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

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.     

Theoretical study of the equilibrium structure, vibrational spectrum, and thermochemistry of the peroxynitrate CF2BrCFBrOONO2

The results of a theoretical study of the molecular structure and conformational mobilities of the peroxynitrate CF(2)BrCFBrOONO(2) and its radical decomposition product CF(2)BrCFBrOO are reported in this paper. The most stable structures were calculated from ab initio G3(MP2)B3 and G4(MP2) methods and from density functional theory at the B3LYP/6-311+G(d) and B3LYP/6-311+G(3df) levels of theory. The equilibrium conformation of CF(2)BrCFBrOONO(2) indicates that the bromine atoms lie in position anti to each other and possess a COON dihedral angle of 114°. A quantum statistical analysis shows that about 40% of the internal rotors can freely rotate at room temperature. Our best values for the standard enthalpies of formation of CF(2)BrCFBrOONO(2) and CF(2)BrCFBrOO at 298 K obtained from isodesmic reactions at the G3(MP2)//B3LYP/6-311+G(3df) level of theory are -144.7 and -127.0 kcal mol(-1). From these values and the enthalpy of formation of the NO(2) radical, a CF(2)BrCFBrOO-NO(2) bond dissociation enthalpy of 26.0 ± 2 kcal mol(-1) was estimated.     

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.     

B4O分子异构化稳定性的理论研究

采用CCSD(T)/6-311+G(2df)//B3LYP/6-311+G(d)方法, 系统研究了B₄O体系各个异构体的结构和能量, 以及重要异构体的解离和异构化稳定性. 结果表明, 单态平面三角形含-BO单元的异构体cB₃-BO ¹1能量最低, 其次是带状的异构体B₄O ¹2(10.9 kJ/mol), 并且¹1和¹2结构都具有良好的动力学稳定性. 因此对于B₄O体系, ¹1和¹2都是有可能存在的. 而文献报道的三态直线型结构BBBBO(146.0 kJ/mol)的能量比异构体¹1和¹2高得多.     

Consideration of spin states in determining the structure and decomposition of the transition metal pentazoles FeClN5, Fe(N5)2, Fe(H2O)4ClN5, and Fe(NH3)4ClN5

Optimised geometries and energies have been calculated for the unidentate, bidentate, and ferrocene-like structures of Fe(II) pentazoles of the form MN(5) and M'(N(5))(2), (M = FeCl, M'= Fe) using the UB3LYP DFT method with basis sets up to 6-311+G(3df) in order to study the effect of an unfilled d shell on the stability of the singlet, triplet, and quintet states; the ground state is a quintet and the quintet has a lower barrier to decomposition than the singlet.     

Reaction mechanism of the CCN radical with nitric oxide

To investigate the possibility of the carbyne radical CCN in removal of nitric oxide, a detailed computational study is performed at the Gaussian-3//B3LYP/6-31G(d) level on the CCN + NO reaction by constructing the singlet and triplet electronic state [C(2)N(2)O] potential energy surfaces (PESs). The barrierless formation of the chain-like isomers NCCNO (singlet at -106.5, triplet cis at -48.2 and triplet trans at -47.6 kcal/mol) is the most favorable entrance attack on both singlet and triplet PESs. Subsequently, the singlet NCCNO takes an O-transfer to form the branched intermediate singlet NCC(O)N (-85.6), which can lead to the fragments CN + NCO (-51.2) via the intermediate singlet NCOCN (-120.3). The simpler evolution of the triplet NCCNO is the direct N-O rupture to form the weakly bound complex triplet NCCN...O (-56.2) before the final fragmentation to NCCN + (3)O (-53.5). However, the lower lying products (3)NCN + CO (-105.6) and (3)CNN + CO (-74.6) are kinetically much less competitive. All the involved transition states for generation of CN + NCO and NCCN + (3)O lie much lower than the reactants. Thus, the novel reaction CCN + NO can proceed effectively even at low temperatures and is expected to play a role in both combustion and interstellar processes. Significant differences are found on the singlet PES between the CCN + NO and CH + NO reaction mechanisms.     

Transition states for deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of acrylonitrile

The geometries and energetics of transition states (TS) for radical deactivation reactions, including competitive combination and disproportionation reactions, have been studied for the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO)‐mediated free‐radical polymerization of acrylonitrile with quantum mechanical calculations at the DFT/UB3‐LYP/6‐311+G(3df,2p)//(U)AM1 level of theory (where DFT is density functional theory, AM1 is Austin model 1, and UAM1 is unrestricted Austin model 1). A method providing reasonable starting geometries for an effective search for TS between the TEMPO radical and 1‐cyanopropyl radical mimicking the growing polyacrylonitrile macroradical is shown. For the hydrogen atom abstraction reaction by the TEMPO radical from the 1‐cyanopropyl radical, practically one TS has been found, whereas for the combination reaction of the radicals, several TS have been found, mainly differing in out‐of‐plane angle α of the N? O bond in the TEMPO structure. α in the TS is correlated with the activation energy, ΔE, determined from the single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 level for the combination reaction of CH₃AN· with the TEMPO radical. The theoretical activation energy for the coupling reaction from DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 calculations has been estimated to be 11.6 kcal mol^?1, that is, only about 4.5 times smaller than ΔE for the disproportionation reaction obtained with the DFT UB3‐LYP/6‐311+G(3df, 2p)//(U)AM1 approach. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 914–927, 2006     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

Theoretical study on the mechanism of the NCO + HCNO reaction

The complex doublet potential surface of the NCO + HCNO reaction has been investigated at the QCISD(T)/6-311g(d,p)//UB3LYP/6-31G(d,p) level. We have found 29 isomers on the potential surface, which are connected by 38 transition states. The single-point energy calculations are performed at the high-level QCISD(T)/6-311G(d,p) for more accurate energy values. In various possible initial association ways, the end-N attack leading to HC2N2O2 a1 and a2 is the most favorable association way through a barrierless process. Through the thermodynamic and kinetic analyses, the product NO + CO + HCN should be the major product in both the low- and high-temperature conditions for its low-energy determination transition state. Our calculation is consistent with the available data in low-temperature condition and expected to be confirmed in the high-temperature condition.     

The nature and role of 8-deltacyclyl and 7-isodeltacyclyl cations in solvolytic reactions of parent brosylates

Solvolysis of exo-8-deltacyclyl brosylate proceeds directly through a C(2) delocalized cation to exo-8-deltacyclyl acetate. The solvolysis of the endo epimer presents a more complex picture, reacting via a classical deltacyclyl cation, the nonclassical C(2) delocalized cation, and the isodeltacyclyl cation. The solvolysis of exo-7-isodeltacyclyl brosylate generates the 7-isodeltacyclyl ion and subsequently the C(2) delocalized deltacyclyl cation forming 8-deltacyclyl acetate and exo-7-isodeltacyclyl acetate. The cationic intermediates in these three related reaction manifolds are characterized using density functional methods at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) and BPW91/6-311+G(3df, 2p)//BPW91/6-31G(d) levels.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

The Dowd-Beckwith ring expansion: a theoretical study

UB3LYP/6-311++G(d,p) and ROMP2/6-311++G(d,p)//UB3LYP/6-311++G(d,p) calculations including the effect of benzene solvent through the PCM-UAHF method render a concerted mechanism without fragmentation as the most favourable one for the Dowd-Beckwith radical ring expansion of the bromomethyl adduct of methyl cyclopentanone-2-carboxylate to yield methyl cyclohexanone-3-carboxylate. The corresponding concerted TS is a bicyclic alcoxy radical.     

Remarks on GIAO-DFT predictions of 13C chemical shifts

Predicting (13)C chemical shifts by GIAO-DFT calculations appears to be more accurate than frequently expected provided that: (a) the comparison between experimental and theoretical data is performed using the linear regression method, (b) a sufficiently high level of theory [e.g. B3LYP/6-311 + + G(2d,p)//B3LYP/6-311 + + G(2d,p) or PBE1PBE/6-311 + G(2df,p)//B3LYP/6-311 + + G(2d,p)] is used, (c) the experimental data originate from the measurements performed in one solvent whose influence is taken into account at the molecular geometry optimization step and, first of all, during the shielding calculation, (d) the experimental data are free of heavy atom effects or such effects are appropriately treated in calculations, and finally (e) the conformational compositions of the investigated objects are known.     

烯丙基自由基(·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反应更具有实际意义.