Molecular structures and electronic properties of helical thiophene carbon–sulfur oligomers, H2(C2S)nC2H2 (n?=?1–20)

Structure optimizations of the thiophene carbon–sulfur H₂(C₂S) _n C₂H₂ (n = 1–20) were carried out using density functional theory calculations at the B3LYP/6-31G(d) level. The B3LYP/6-31G(d) geometrical data for heptamer H₂(C₂S)₇C₂H₂ and undecamer H₂(C₂S)₁₁C₂H₂ are in good agreement with the X-ray crystallographic data for the helical (C₂S)_n β-heptathiophene and β-undecathiophene, respectively. Structural and electronic properties of helical oligothiophenes obtained at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level are reported. The strain energy formula of n oligothiophenes as a linear function of their molecular length was obtained.     

Molecular Structure and Conformational Composition of 1-[(Methylthio)methyl]-2-nitrobenzene (MTMNB): A Theoretical and Experimental Study

The conformational composition of gaseous MTMNB and the molecular structures of the rotational forms have been studied by electron diffraction at 130^∘C aided by results from ab initio and density functional theory calculations. The conformational potential energy surface has been investigated by using the B3LYP/6-31G(d,p) method. As a result, six minimum-energy conformers have been identified. Geometries of all conformers were optimized using MP2/6-31G(d,p), B3LYP/6-31G(d,p), and B3LYP/cc-pVTZ methods. These calculations resulted in accurate geometries, relative energies, and harmonic vibrational frequencies for all conformers. The B3LYP/cc-pVTZ energies were then used to calculate the Boltzmann distribution of conformers. The best fit of the electron diffraction data to calculated values was obtained for the six conformer model, in agreement with the theoretical predictions. Average parameter values (r_a in angstroms, angle α in degrees, and estimated total errors given in parentheses) weighted for the mixture of six conformers are r(C–C) = 1.507(5), r(C–C)_{ring, av} = 1.397(3), r(C–S)_av = 1.814(4), r(C–N) = 1.495(4), r(N–O)_av = 1.223(3), ∠(C–C–C)_ring = 116.0–122.5, ∠ C6–C4–C7 = 118.2(4), ∠ C–C–S = 113.6(6), ∠ C–S–C = 98.5(12), ∠ N–C–C4 = 121.9(3), ∠(O–N–C)_av = 116.8(3), ∠ O–N–O = 127.0(4). Torsional angles could not be refined. Theoretical B3LYP/cc-pVTZ torsional angles for the rotation about C–N bond, φ_C−N, were found to be 30.5–36.5^∘ for different conformers. As to internal rotation about C–C and C–S bonds, values of φ_C−C = 68–118^∘ and φ_C−S = 66–71^∘ were obtained for the three most stable conformers with gauche orientation with respect to these bonds. Some conclusions of this work were presented in a short communication in Russ. J. Phys. Chem. 2005, 79, 1701.     

Conformational properties of <Emphasis Type="Italic">ortho</Emphasis>-nitrobenzenesulfonamide in gas and crystalline phases. Intra- and intermolecular hydrogen bond

A combined gas-phase electron diffraction and quantum chemical (B3LYP/6-311+G**, B3LYP/cc-pvtz, MP2/cc-pvtz) study of molecular structure of 2-nitrobenzenesulfonamide (2-NBSA) was carried out. Quantum chemical calculations showed that 2-NBSA has four conformers, two of which are stabilized by intramolecular hydrogen bond. The latter (with the S–N bond in a close to orthogonal position around the phenyl ring and differing from each other by staggered or eclipsed positions of the N–H and S=O bonds in the SO₂NH₂ group) presented in a saturated vapor over 2-NBSA at T = 433 (3) K in commensurable amounts. Experimental internuclear distances (Ǻ) for the staggered conformer are (?): r _h1(C–H)_av. = 1.071(9), r _h1(C–C)_av. = 1.390(4), r _h1(C–S) = 1.789(8), r _h1(S=O)_av. = 1.427(6), r _h1(S–N) = 1.644(6), r _h1(N–O)_av. = 1.221(4), r _h1(C′–N) = 1.487(8), r _h1(N–H)_av. = 1.014. Calculations at B3LYP/cc-pvtz level were performed to determine the structure and the energies of the transition states between conformers. It was shown that the conformer structures of free molecule differ from those of a molecule stabilized by intermolecular hydrogen bonds in a crystal. Influence of a substituent X (X = –CH₃, –NO₂) on conformational features of the ortho-substituted benzenesulfonamide was established.     

Crystallographic and theoretical studies of 4,4′-dimethyl-7,7′-bi([1,2,5]thiadiazolo[3,4-b]pyridylidene)–chloranilic acid (1/1) with intermolecular O–H···N hydrogen bonds and S···O heteroatom interactions

Abstract  The molecular and crystal structure of a 1:1 co-crystal of 4,4′-dimethyl-7,7′-bi([1,2,5]thiadiazolo[3,4-b]pyridylidene)–chloranilic acid, (1), has been determined by X-ray diffraction at the monoclinic space group P2₁/c with cell parameters of a = 8.422(6), b = 7.343(4), c = 16.112(7) ?, β = 104.988(8)°, V = 962.5(10) ?³ and Z = 2. In the crystal structure, two components connect via the intermolecular O–H···N hydrogen bonds [2.804(4) ?] and S···O heteroatom interaction [2.945(3) ?] with R ₂ ²(7) couplings to form a unique and infinite one-dimensional supramolecular tape structure. The calculations of (1) at the HF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels can almost reproduce X-ray geometry. In addition, the distances of the intermolecular O–H···N and S···O interactions by MP2/6-31G(d) and B3LYP/6-31G(d) levels agree well with those in the crystal. The calculated binding energies corrected BSSE and ZPE are −4.487 (HF), −7.473 (MP2), and −5.640 (B3LYP) kcal/mol. The results suggest that the complex (1) is very stable and the dispersion interaction is significantly important for the attractive intermolecular interaction in (1). The NBO analysis has revealed that the n(N) → σ*(O–H) interaction gives the strongest stabilization to the system and the major interaction for the intermolecular S···O contact is n(O) → σ*(S–N). Index Abstract  In the crystal structure of the title compound, the molecules are linked by intermolecular O–H···N hydrogen bonds and short S···O heteroatom interactions with R ₂ ²(7) couplings to construct a unique and infinite one-dimensional supramolecular tape structure.      

Modeling Ar and Kr matrix effect on the ν s (Cl–H) and ν l (Cl–H) of Cl–H···NH3 by the IEF-PCM method

Ar and Kr matrix effect on the geometry and Cl–H stretching (ν _s (Cl–H)) and librational (ν _l (Cl–H)) frequencies of the hydrogen-bonded complex Cl–H···NH₃ are simulated within the framework of polarizable continuum model with integral equation formalism (IEF-PCM) at B3LYP and MP2 levels of theory with the basis set 6-311++G(2df,2pd). Within the framework of B3LYP and IEF-PCM, the simulated gas phase, Ar, and Kr matrix ν _s (Cl–H) of the complex are 2140, 1684, and 1550 cm⁻¹, respectively, which deviate from the experimental values (~2200, 1371, and 1218 cm⁻¹) by −60, 313, and 332 cm⁻¹. Within the framework of MP2 and IEF-PCM, the gas phase, Ar, and Kr matrix ν _s (Cl–H) are calculated as 2366, 2037, and 1957 cm⁻¹ by the harmonic approximation, and as 2177, 1876, and 1665 cm⁻¹ by the full-dimensional anharmonic correction. The matrix effect modeling is of greater importance than the anharmonic correction in accounting for the large experimental gas phase to Ar or Kr matrix shift of the ν _s (Cl–H) (−829 or −982 cm⁻¹). Our calculations do not support the assignment of the 733.8 and 736.9 cm⁻¹ bands to the Ar and Kr matrix ν _l (Cl–H).     

Reactions of cationic rhodium(I) carbonyl compounds with nucleophiles to form alkoxo,carboalkoxy and carboxylato complexes. Part II

Summary The rhodium(I) carbonyl compounds [Rh(CO)L₂2] [BF₄]. 1/2CH₂Cl_nn2 (L = PPh₂ or AsPh₃) react with the nucleophiles OMe⁻, RCOO⁻ (R = Me, Et) under nitrogen to form [Rh(OR)(CO)L₂] (1)–(2) and [Rh(OOCR)(CO)L₂] (7)–(₁₀), respectively. Addition of [Rh(CO)₂(PPh₃)₂]-[BF ₄] to OMe⁻ under nitrogen produces [Rh(COOMe)-(CO) (PPh₃)₂]-MeOH (3), whilst reactions of [Rh(CO)-(PPh₃)₂] [BF₄]·1/2CH₂Cl₂ and [Rh(CO)₂(PPh₃)₂] [BF₄] with OR_- (R = Me, Et or n-Pr) in the presence of CO produce [Rh(COOR)(CO)₂(PPh₃)₂] (4)–(6). The products have been characterised by i.r., ¹H, ³¹P, ¹³Cn.m.r. spectroscopy and elemental analysis.     

An ab initio and density functional study of microsolvation of carbon dioxide in water clusters and formation of carbonic acid

Geometry of the CO₂–H₂O complex and reaction barriers leading to the formation of H₂CO₃were studied at the RHF/6-311++G**, MP2/6-311++G**, B3LYP/AUG-cc-pVDZ, B3LYP/AUG-cc-pVTZ, MP2/AUG-cc-pVDZ and CCD/AUG-cc-pVDZ levels of theory. The rotational barrier of the CO₂–H₂O complex and the reaction barrier leading to the formation of H₂CO₃–H₂O from CO₂–(H₂O)₂ were studied using the first three of the above-mentioned methods. Microsolvation of CO₂ in water clusters having upto eight water molecules was studied using the B3LYP/AUG-cc-pVDZ method. Various methods except MP2/AUG-cc-pVDZ predict the equilibrium structure of the CO₂–H₂O complex to be symmetric while the MP2/AUG-cc-pVDZ method predicts it to be unsymmetric. Formation of H₂CO₃ from CO₂–H₂O is strongly catalyzed by the presence of a second water molecule. Atomic orbitals are strongly rehybridized in going from the equilibrium structures of the CO₂–H₂O and CO₂–(H₂O)₂ complexes to the transition states involved in the formation of H₂CO₃ and H₂CO₃–H₂O, respectively, as shown by hybridization displacement charges.     

The influence of steric hindrance on conformation properties and molecular structure of 2,4,6-trinitrobenzenesulfonic acid: gas electron diffraction and quantum chemical calculations

A combined gas-phase electron diffraction and quantum chemical (B3LYP/cc-pVTZ, MP2/cc-pVDZ) study of molecular structure of 2,4,6-trinitrobenzenesulfonic acid (2,4,6-tri-NBSA) was carried out. Quantum chemical calculations showed that 2,4,6-tri-NBSA possesses six conformers, which form three pairs of enantiomers with the relative energy of 0, 4.4/3.9, and 2.5/2.5 kcal/mol. It was experimentally established that at T = 444(5) K a saturated vapor over 2,4,6-tri-NBSA is, predominantly (up to 93 mol.%), represented by a low-energy enantiomers II and II′ characterized by intramolecular hydrogen bond between an H atom of the hydroxyl group and one of the O atoms of the NO₂ group. Experimental internuclear distances for the low-energy enantiomers are (?): r _h1(C–C)_av. = 1.387(4), r _h1(C–S) = 1.811(6), r _h1(S=O)_av. = 1.424(4), r _h1(S–O) = 1.579(4), r _h1(N–O)_av. = 1.214(3), r _h1(C–N)_av. = 1.491(5). Geometry of the conformer II points on existance of strong steric interactions between SO₂OH group and two ortho-nitro groups. Analysis of the orbital interactions between the substituents and benzene ring was carried out. Geometric parameters and energies of transition states between conformers were calculated (B3LYP).     

Conformational analysis and interpretation of ν(OH) band in the IR spectrum of o-vinylphenol: a DFT study

Two-dimensional cyclic potential energy surface for internal rotation of vinyl and hydroxyl substituents in o-vinylphenol molecule was constructed by the B3LYP/6-311G(d) method. It was shown that o-vinylphenol molecule exists in the gas phase as a mixture of seven rotamers denoted as A (A′), B (B′), C (C′) and D. The B3LYP/cc-pVTZ calculated percentage of the rotamers A and A′ in which OH…π intramolecular interaction occurs, is at most 24%. The height of barriers t interconversions between o-vinylphenol rotamers varies from 0.1 to 5.2 kcal mol⁻¹. According to B3LYP/cc-pVTZ calculations, the inclusion of solvent effect in the framework of the polarizable continuum model for a solution of o-vinylphenol in CCl₄ leads to a decrease in theoretical values of ν(OH) frequencies by about 4–9 cm⁻¹ and to an increase in the percentage of the rotamers without intramolecular hydrogen bond by about 4.3% compared to the corresponding gas-phase values. The simulated IR spectral contours of ν(OH) bands are in good agreement with experimental one Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 499–508, March, 2008.     

Theoretical study of reaction mechanism for Se + O<Subscript>3</Subscript>

The reaction mechanism of Se + O₃ on the singlet potential energy has been investigated at CCSD(T)/6-311++G(2df,2pd) level of theory based on the geometric parameters optimized at the B3LYP/6-311++G(3df,3pd) level of theory. The calculated results show that the reactants are firstly associated into the adduct Se–O₃ with any intrinsic barrier. Subsequently, through a variety of transformations of isomer Se–O₃, two kinds of products P₁(SeO₃(D_3h)) and P₂(SeO + ³O₂) are obtained. The breakage and formation of the chemical bonds in the reaction have been studied by the topological analysis of electronic density. The topological analysis results show that the ring transitional structure region does not only occur in cis-OSeOO → SeO₃(C_s) process but also occur in SeO₃(C_s) → SeO₃(D_3h).     

Theoretical study of structure and properties of hexuronic acid and <Emphasis Type="SmallCaps">d</Emphasis>-glucosamine structural units of glycosaminoglycans

The Becke3LYP functional of DFT theory was used to investigate molecular structure and sodium affinity of the systems CH₃CO₂Na (1), CH₃–O–SO₃Na (2), CH₃–NH–SO₃Na (3), saccharide_1Na₂ (4), saccharide_2Na (5), saccharide_3Na₃ (6), saccharide_4Na₂ (7), and saccharide_5Na₂ (8), respectively, which are models of N- and O-sulfate glycosaminoglycans. Interaction enthalpies, entropies and Gibbs energies of the sodium-coordinated systems in the gas phase were determined with the B3LYP/6-311+G(d,p) and B3LYP/6-31+G(d) methods. The computed Gibbs energies, ΔG ^o , of model systems 1–3 are negative and span a rather broad energy interval (from −500 to −1,500 kJ mol⁻¹). Gibbs interaction energies for sodium acetate, sodium sulfate and sodium sulfamate functions of the five saccharides, systems 4–8 are always lower than those values found for the model compounds 1–3. The ionization of sodium salts of saccharides studied in gas phase is in most cases connected with considerable conformational rearrangement of the ionic species. This rearrangement causes an additional energetic stabilization of anionic species and is connected with the substantial release of entropy.     

DFT studies of the phosphinidene derivative HPNaF and its insertion reaction with R–H(R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>)

The formations of the phosphinidene derivative HPNaF and its insertion reactions with R–H (R=F, OH, NH₂, CH₃) have been systematically investigated employing the density functional theory (DFT), such as the B3LYP and MPW1PW91 methods. A comparison with the results of MP2 calculations shows that MPW1PW91 underestimates the barrier heights for the four reactions considered. Similarly, the same is also true for the B3LYP method depending on the selected reactions, but by much less than MPW1PW91, where the barrier heights of the four reactions are 25.2, 85.7, 119.0, and 142.4 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. All the mechanisms of the four reactions are identical to each other, i.e., an intermediate has been located during the insertion reaction. Then, the intermediate could dissociate to substituted phosphinidane(H₂RP) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are −92.2, −68.1, −57.2, and −44.3 kJ/mol at the B3LYP/6-311+G* level of theory, respectively, where both the B3LYP and MPW1PW91 methods underestimate the reaction energies compared with the MP2 results. The linear correlations between the calculated barrier heights and the reaction energies have also been observed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Influence of diffuse and polarization functions on the second-order Møller–Plesset optimized dihedral angle of biphenyl

 Using 6-31G and 6-311G basis sets to which diffuse and polarization functions were added in a stepwise fashion (a total of 16 basis sets), Hartree–Fock (HF), MP2 and B3LYP geometry optimizations were performed on biphenyl. With the MP2 method, diffuse functions raise the dihedral angle φ, for example, from 46.3° for 6-31G to 54.1° for 6-311++G, while polarization functions lower it, for example, from 54.1° for 6-311++G to 42.1° for 6-311++G(2d,2p). For a single set of polarization functions, φ(MP2) lies close to or above φ(HF) (44–47°), but for a double set it is below φ(HF) and is close to B3LYP values (38–42°) which show little basis set dependence. The most reliable value for φ, 42.1° [MP2/6-311++G(2d,2p)], is expected to increase slightly by adding more diffuse functions. The corresponding best calculated energy barrier at 0° (coplanar conformation) is 2.83 kcal/mol, much higher than the experimental estimate (1.4 ± 0.5 kcal/mol). The barrier at 90° is 1.82 kcal/mol, in line with the experimental estimate (1.6 ± 0.5 kcal/mol) and with previous theoretical results. Received: 9 September 2002 / Accepted: 15 November 2002 / Published online: 1 April 2003 Correspondence to: Friedrich Grein e-mail: fritz@unb.ca Acknowledgement. The author would like to thank NSERC (Canada) for financial support.     

The molecular structure of N-methylsuccinimide studied by gas-phase electron diffraction (GED) and quantum-chemical methods

The molecular structure of N-methylsuccinimide was studied by the GED method at a nozzle temperature of 69–73°C. Anharmonic vibrational corrections to thermal-average r _a bond lengths, Δ(r _a − r _e), were calculated using the quadratic and cubic force constants from B3LYP/6-31G(df, p) calculations. The molecular skeleton was found to be planar within measurement errors. Some structural effects were likely caused by the conjugation of the N atom with two C=O bonds. The equilibrium geometric parameters derived from the experimental data and those from MP2/cc-pVTZ(seg-opt) calculations were in close agreement.     

Effects of fused benzene rings on tautomerizations and inversions of benzo, azabenzo, and oxabenzocycloheptatrienes at theoretical levels

Biologically important bicyclic species, including 6H-, 6H-6-aza-, and 6-oxabenzocycloheptatrienes (in which the benzene moiety is fused meta with respect to the tetrahedral constituents: –CH₂–, –NH–, and –O–, respectively), show strong shifts of tautomerizations in favor of the corresponding tricyclic benzonorcaradienes (with ΔH values of −11.49, −14.55, and −19.20 kcal mol⁻¹, respectively), at B3LYP/6-311++G**//B3LYP/6-31G*, and MP2/6-311++G**//MP2/6-31G* levels, and at 298 K. In contrast, such shifts are strongly disfavored by the isomeric bicyclic species in which the benzene moieties are fused ortho or para with respect to –CH₂–, –NH–, and –O–, respectively. Hence for species with ortho benzene rings including 5H-, 5H-5-aza- and 5-oxabenzocycloheptatrienes, tautomerization ΔH values are 30.76, 31.89, and 25.27 kcal mol⁻¹, respectively, while for species with para fused benzene moieties including 7H-, 7H-7-aza-, and 7-oxabenzocycloheptatrienes, tautomerization ΔH values are 24.12, 26.00, and 19.55 kcal mol⁻¹, respectively. NICS calculations are successfully used to rationalize these results. The calculated energy barriers for inversion of the seven-membered rings of bicyclic species predict a dynamic nature for all the structures except for the virtually planar 6H-6-aza- and 6-oxabenzocycloheptatrienes. Finally, our theoretical data are compared to the experimental results where available. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Experimental and theoretical studies on the reactions of 1-amino-5-benzoyl-4-phenyl-1<Emphasis Type="Italic">H</Emphasis>-pyrimidine-2-one/-thione compounds with ethyl acetoacetate

The reaction of 1-amino-5-(4-methylbenzoyl)-4-(4-methylphenyl)pyrimidin-2(1H)-one/-thione (1a,b) with ethyl acetoacetate (EA) afforded moderate to good yields (59–63%) of ethyl 2-methyl-4-(4-methylbenzoyl)-5-(4-methylphenyl)-7-oxo/-thioxo-3,3a-dihydropyrazolo[1,5-c]pyrimidine-3-carboxylate (2a,b). The newly synthesized compounds were characterized by elemental analyses, IR, ¹H and ¹³C NMR spectral data. All were compared with their previous analogues. The reaction mechanism of 1 with EA was studied by means of the B3LYP/6-31G(d,p) method. In addition, for reactants Fukui functions were performed using the data calculated with the Becke3–Lee–Yang–Parr (B3LYP) hybrid function.     

Molecular structure, barriers to internal rotation, and the enthalpy of formation of cumene hydroperoxide PhCMe2OOH: a DFT study

Conformational analysis of cumene hydroperoxide PhCMe₂OOH (1) has been carried out using the density functional methods B3LYP/6-31G(d,p) and B3LYP/6-311+G(3df,2p). Ignoring rotation of methyl groups, molecule 1 has seven conformers differing in orientation of the — CMe₂OOH fragment relative to the benzene ring and in mutual position of atoms in this fragment. The molecular structures, relative energies, and statistical distribution of the conformers were determined, and intramolecular rotational barriers were estimated. The enthalpies of formation of all conformers of molecule 1 were calculated using two approximations with inclusion of zero-point vibrational energy and temperature correction. Calculations using the isodesmic reaction (IDR) scheme made it possible to reduce the systematic error of the determination of the enthalpy of reactions. The total enthalpy of formation of compound 1 calculated with inclusion of statistical distribution of rotamers equals −19.7±3.6 kcal mol⁻¹. The combination of the B3LYP/6-31G(d,p) approximation and the IDR scheme gives fairly accurate results (relative error is ±0.4 kcal mol⁻¹) as compared to those obtained with the extended basis set 6-311+G(3df,2p). Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1157–1164, June, 2008.     

[C5Li5]Mgn[C5Li5] (n?=?2–8): novel sandwich complexes containing –Mg–Mg– chain

A density functional theory (DFT) investigation on novel sandwich-type D ₅ [C₅Li₅]Mg _n [C₅Li₅] (n = 2–8) complexes containing –Mg–Mg– chain has been performed in this work. The equilibrium geometries, electronic structures, vibrational frequencies, and stabilities of these complexes are researched by B3LYP and BP86 methods at 6-311+G(d) levels of theory. The Mg _n ²⁺ sandwich complexes with D ₅ symmetry are all true minima on the potential energy surface. NBO analyses for the series of complexes reveal that the Mg–Mg bond is a weak σ covalent bond. There are mainly electrostatic interactions between C₅Li₅ ⁻ ligands and Mg _n ²⁺(n = 2–8) nuclear in these complexes. The NICS and NICS_zz computed with GIAO-B3LYP/6-311+G(d) indicates that the C₅Li₅ ⁻ rings in the series of complexes are aromatic. These novel complexes turn out to be strongly thermodynamically favored in the gas phases and may be targeted in future experiments to expand the structural domain of sandwich-type complexes.     

Theoretical studies on the mechanisms of NCCO + O2 reaction

The mechanisms of the reaction of NCCO with molecular oxygen are investigated at the G3MP2//B3LYP/6-311G(d,p) levels for the first time. The calculation results show that two mechanisms are involved, namely, O attack on α atom mechanism and O attack on β atom mechanism, with six products yielded. The most feasible channel is the addition of O₂ to β atom in NCCO radical leading to the energy-rich intermediate IM1, NCC(O)OO, which can isomerize to a four-center-structure IM3, and then undergoes C–C and O–C bond fission to form P1(NCO + CO₂) finally. The barriers are 27.3 and 25.4 kcal/mol, respectively. For other channels involved in the two mechanisms, with less stable initial adducts and higher barrier, they are less conceivable dynamically and thermochemically.     

Theoretical study of normal and inverse electron demand cycloaddition reactions between cyclohexadiene and 2-aryl-4,6-dinitrobenzotriazole 1-oxides

In this work, we theoretically study Diels-Alder reactions between cyclohexadiene and a series of 2-aryl-4,6-dinitrobenzotriazole 1-oxides (5a-d). Based on various calculation methods such as SCF/6–31G* and DFT/B3LYP with the 6–31G* standard basis set, we discuss the thermodynamic possibility of these reactions. The reactions between a series of 6-nitrobenzotriazole 1-oxides(5aa ₁ -ab ₁ ) and cyclohexadiene are also studied. The text was submitted by the authors in English.