Thermal isomerization and decomposition of 1,2-butadiene in shock waves

1,2-Butadiene diluted with Ar was heated behind reflected shock waves over the temperature and the total density range of 1100–1600 K and 1.36 × 10^?5 ? 1.75 × 10^?5 mol/cm³. The major products were 1,3-butadiene, 1-butyne, 2-butyne, vinylacetylene, diacetylene, allene, propyne, C₂H₆, C₂H₄, CH₄, and benzene, which were analyzed by gas chromatography. The UV kinetic absorption spectroscopy at 230 nm showed that 1,2-butadiene rapidly isomerizes to 1,3-butadiene from the initial stage of the reaction above 1200 K. In order to interpret the formation of 1,3-butadiene, 1-butyne, and 2-butyne, it was necessary to include the parallel isomerizations of 1,2-butadiene to these isomers. The present data were successfuly modeled with a 82 reaction mechanism. From the modeling, rate constant expressions were derived for the isomerization 1,2-butadiene = 1,3-butadiene to be k₃ = 2.5 × 10¹³ exp(?63 kcal/RT) s^?1 and for the decomposition 1,2-butadiene = C₃H₃ + CH₃ to be k₆ = 2.0 × 10¹⁵ exp(?75 kcal/RT) s^?1, where the activation energies, 63 kcal/mol and 75 kcal/mol, were assumed. These rate constants are only applicable under the present experimental conditions, 1100–1600 K and 1.23–2.30 atm. © 1995 John Wiley & Sons, Inc.     

1-三氯锡烷基-2,3-丁二烯和2-三氯锡烷基-1,3-丁二烯与甲醛反应的B3LYP研究

用密度泛函方法在B3LYP/6-31G^**水平上研究了1-三氯锡烷基-2,3-丁二烯和2-三氯锡烷基-1,3-丁二烯与甲醛的反应.优化得到各驻点的几何构型,通过振动分析和内禀反应坐标对过渡态进行了确认,解析了反应路径.并用SCRF(PCM)方法在同一水平上对在CH₂Cl₂溶液中的两反应进行了研究.计算了两反应在气相和CH₂Cl₂溶液中的活化能垒、自由能和平衡常数.结果表明,反应具有很强的选择性,主要得到1-三氯锡烷基-2,3-丁二烯与甲醛反应的产物.该结果与实验事实一致.     

Rotation about the C?N bond in 2-aza-1,3-butadiene and the N?N bond in 2,3-diaza-1,3 butadiene: A molecular orbital study

The geometry and energy of 2-aza-1,3-butadiene and 2,3-diaza-1,3-butadiene have been calculated using the 6-31G* basis set as a function of the CNCC and CNNC dihedral angles, respectively. With the 2-aza derivative potential minima are located at 0° (trans) and at about 130° for a gauche structure approximately 9.5 kJ mol^?1 less stable than the trans. Potential maxima are at about 75° giving a gauche barrier height of approximately 19 kJ mol^?1 relative to the trans structure, and at 180° (cis) giving a barrier height of approximately 14.5 kJ mol^?1 relative to the 130° gauche structure. With the 2,3-diaza derivative the gauche barrier has disappeared and there are a series of gauche structures in the region 70°–100° of almost equal energy 12.5-15 kJ mol^?1 less stable than the trans. In addition the cis barrier is much greater, nearly 70 kJ mol^?1 relative to the trans structure. Inclusion of electron correlation, accounting for about 50% of the correlation energy, produces no significant changes in the shape of the potential energy curves. There are systematic and progressive changes in almost all the geometrical parameters as the ?CH? groups in butadiene are replaced by ?N? . The outward tilt and compression within the methylene groups show adverse steric interactions to be operative in the cis structures. The values of V_nn indicate that gauche structures of both the 2-aza and the 2,3-diaza derivatives near the cis structure are more compact (as with butadiene), and gauche structures of the 2-aza derivative near the trans structure are less compact (as with butadiene). Originating in the changes in bond lengths and bond angles, rotation-independent nuclear–nuclear interactions again play an important role.     

Ab initio study on the internal rotation of five π-conjugated hydrocarbons at MP2 level

The potential-energy curves of internal rotation were calculated for 1,3-butadiene at the MP2/6-311G** level, for isoprene and 1,3-pentadiene at the MP2/6-311G* level, and for 2,3-dimethyl-1,3-butadiene and styrene at the MP2/6-31G* level. The geometries of the energy minima (stable conformers) and maxima (transition states) on the curves are completely optimized. For butadiene and its methyl derivatives, two stable rotamers, s-trans and gauche conformers, are obtained. s-trans forms have the lowest energies and gauche conformers twisted by 39.9°–48.3° around the central bond of the butadiene skeleton are, on average, 9.8 kJ/mol above the trans forms. s-cis forms are rotational transition states. The computed gauche–cis barriers range from 4.30 to 11.70 kJ/mol. The regular effects of methyl substitutions at the end and central carbons are found. For styrene, the planar form is calculated to be a saddle point which is only about 1 kJ/mol higher in total energy than a twisted minimum, in which the torsional angle between the phenyl and vinyl planes is 27.4°. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 659–667, 1998     

Comparison of two-body and three-body decomposition of ethanedial, propanal, propenal, n-butane, 1-butene, and 1,3-butadiene

We investigated two-body (binary) and three-body (triple) dissociations of ethanedial, propanal, propenal, n-butane, 1-butene, and 1,3-butadiene on the ground potential-energy surfaces using quantum-chemical and Rice-Ramsperger-Kassel-Marcus calculations; most attention is paid on the triple dissociation mechanisms. The triple dissociation includes elimination of a hydrogen molecule from a combination of two separate terminal hydrogen atoms; meanwhile, the rest part simultaneously decomposes to two stable fragments, e.g., C(2)H(4), C(2)H(2), or CO. Transition structures corresponding to the concerted triple dissociation were identified using the B3LYP/6-311G(d,p) level of theory and total energies were computed using the method CCSD(T)/6-311+G(3df, 2p). The forward barrier height of triple dissociation has a trend of ethanedial < propanal < propenal < n-butane < 1-butene < 1,3-butadiene, pertaining to the reaction enthalpy. Ratios of translational energies of three separate fragments could be estimated from the transition structure of triple dissociation. The synchronous concerted dissociation of propanal, propenal, and 1-butene leading to three different types of molecular fragments by breaking nonequivalent chemical bonds is rare. The triple dissociation of propanal, n-butane, 1-butene, and 1,3-butadiene were investigated for the first time. To outline a whole picture of dissociation mechanisms, some significant two-body dissociation channels were investigated for the calculations of product branching ratios. The triple dissociation plays an important role in the three carbonyl compounds, but plays a minor or negligible role in the three hydrocarbons.     

Refined ab initio 6-31G split-valence basis set optimization of the molecular structures of biphenyl in twisted,planar, and perpendicular conformations

Improved full ab initio optimizations of the molecular structure of biphenyl in twisted minimum energy, coplanar, and perpendicular conformations by use of Poles's GAUSSIAN 82 program have been performed in the 6-31G basis set. These lead to geometries and energies of much higher reliability than our earlier STO-3G results. The torsional angle Φ_min obtained now is 45.41° in close agreement with the recent experimental value of 44.4° ± 1.2°. Calculated CC distances may be converted to experimental ED r_g-values by means of independently determined linear regression correlations with very high statistical confidence, although they agree better with experimental x ray data for coplanar biphenyl without this correction. Calculated intramolecular angles are very similar for both STO-3G and 6-31G basis sets. The calculated torsional energy barrier towards Φ = 90° (ΔE⁹⁰) is 6.76 kJ/mol in close agreement with the experimental-31G value of 6.5 ± 2.0 kJ/mol. For coplanar biphenyl with D_2h-symmetry the calculated torsional energy barrier ΔE⁰ is 13.26 kJ/mol which is surprisingly much higher than the experimental value of 6.0 ± 2.1 kJ/mol. This discrepancy could not be resolved by optimizations assumed for two kinds of distortions of planarity of orthohydrogens from the molecular plane of the coplanar carbon atoms. But for the twisted minimum energy conformation asymmetric bending of ortho-H atoms lead to a torsional angle Φ_min = 44.74° together with a dihedral angle towards ortho-H of 1.22°, and consequently even to an increase of torsional energy barriers to ΔE⁰ = 13.51 and ΔE⁹⁰ = 6.91 kJ/mol.     

Reaction‐path dynamics and theoretical rate constants for the reaction CH4 + O3 → HOOO + CH3

We present a direct ab initio dynamics study of thermal rate constants of the hydrogen abstraction reaction of CH₄ + O₃ → HOOO +CH₃. The geometries of all the stationary points are optimized at MPW1K/6‐31+G(d,p), MPWB1K/6‐31+G(d,p), and BHandHLYP/6‐31+G(d,p) levels of theory. The energies are refined at a multi‐high‐level method. The extended Arrhenius expression fitted from the CVT/SCT and μVT/Eckart rate constants of ozonolysis of methane in the temperature range 200–2500 K are k^CVT/SCT(T) = 5.96 × 10^?29T^4.49e^(?17321.3/T) and k^μVT/Eckart(T) = 7.92 × 10^?29T^4.46e^(?17301.7/T), respectively. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

<Emphasis Type="Italic">Ab initio</Emphasis> study of the reaction path of the reaction HNCO+CN

We report ab initio UMP2 calculations of the reaction of CN with HNCO using 6-311G(d,p) basis sets. The obtained results show that the reaction has two product channels: HNCO+CN→HCN+NCO (1) and HNCO+CN→HNCN+CO (2). Channel (1) is a hydrogen abstraction reaction, which is a concerted process. The calculated potential energy barrier is 20.80 kJ/mol at UMP2(full)/6-311G(d,p) level. In the range of reaction temperature (1000-2100 K), the conventional transition theory rate constant for channel (1) ranges from 0.32×10⁻¹¹ to 6.9×10⁻¹¹cm³· mol⁻¹· s⁻¹, which is close to the experimental value. Channel (2) is a stepwise reaction involving an intermediate during the process of reaction. The UMP2(full)/6-311G(d,p) potential energy barrier is 83.42 kJ/mol for the rate-controlling step, which is much higher than that of channel (1).     

A computational study of the reaction kinetics of methyl radicals with trifluorohalomethanes

Ab initio calculations have been used to characterize the transition states for halogen abstraction by CH₃ in reactions with CF₄, CF₃Cl, CF₃Br, and CF₃I (1–4). Geometries and frequencies were obtained at the HF/6-31G(d) and MP2=full/6-31G(d) levels of theory. Energy barriers were computed via the Gaussian-2 methodology, and the results were employed in transition state theory analyses to obtain the rate constants over 298–2500 K. There is good accord with literature measurements in the approximate temperature range 360–500 K for reactions (2–4), and the computed activation energies are accurate to within ±6 kJ mol⁻¹. Recommended rate constant expressions for use in combustion modeling are k;₁=1.6×10⁻¹⁹ (T/K)^2.41 exp(−13150 K/T), k₂=8.4×10⁻²⁰(T/K)^2.34 exp(−5000 K/T), k₃=4.6×10⁻¹⁹ (T/K)^2.05 exp(−3990 K/T), and k₄=8.3×10⁻¹⁹ (T/K)^2.18 exp(−1870 K/T) cm³ molecule⁻¹ s⁻¹. The results are discussed in the context of flame suppression chemistry. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet 30: 179–184, 1998.     

Ab initio potential energy surface and excited vibrational states for the electronic ground state of Li2H

A 285-point multi-reference configuration-interaction involving single and double excitations (MRS-DCI) potential energy surface for the electronic ground state of Li₂H is determined by using 6-311G (2df, 2pd) basis set. A Simons-Parr-Finlan polynomial expansion is used to fit the discrete surface with a X² of 4.64 × 10^-6. The equilibrium geometry occurs at R_e =0.172 nm and <LiHLi =94.10. The dissociation energy for reaction Li₂H(²A)⇑ Li₂(¹⌆_g)+H(²S) is 243.910 kJ/mol. and that for reaction Li₂H(²A)⇑HLi(¹be)+Li(²S) is 106.445 kJ/mol. The inversion barrier height is 50.388 kJ/mol. The vibrational energy levels are calculated using the discrete variable representation (DVR) method. Project supported by the National Natural Science Foundation of China (grant No. 29673029) and by the Special Doctoral Research Foundation of the State Education Commission of China.     

Ab initio study on the mechanism of reaction HNCO+NH2

Ab initio UMP2 and UQCISD(T) calculations, with 6-311G** basis sets, were performed for the titled reactions. The results show that the reactions have two product channels: NH₂+ HNCO?NH₃+NCO (1) and NH₂+HNCO?N₂H₃+CO (2), where reaction (1) is a hydrogen abstraction reaction via an H-bonded complex (HBC), lowering the energy by 32.48 kJ/mol relative to reactants. The calculated QCISD(T)//MP2(full) energy barrier is 29.04 kJ/mol, which is in excellent accordance with the experimental value of 29.09 kJ/mol. In the range of reaction temperature 2300–2700 K, transition theory rate constant for reaction (1) is 1.68×10¹¹–3.29×10¹¹ mL·mol^-1·s^-1, which is close to the experimental one of 5.0×10¹¹mL·mol^-1·s^-1or less. However, reaction (2) is a stepwise reaction proceeding via two orientation modes,cis andtrans, and the energy barriers for the rate-control step at our best calculations are 92.79 kJ/mol (forcis-mode) and 147.43 kJ/mol (fortrans-mode), respectively, which is much higher than reaction (1). So reaction (1) is the main channel for the titled reaction.     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Thermal dimerization of 1,3-butadiene: Kinetics of the formation of cis,cis-cycloocta-1,5-diene

The thermal reaction of 1,3-butadiene (BD) has been studied between 464 and 557°K at pressures between 49 and 450 torr. The products are 4-vinylcyclohexene (VCH) and cis, cis-cycloocta-1,5-diene (COD), and their formations are second order. The rate constant (in 1/mol · sec) for VCH is given by and that for COD by The thermal reaction of COD has also been studied. The temperature was varied from 505 to 586°K and the pressure from 15 to 51 torr. The rate constants (in sec^?1) for the formations of VCH and BD are given by A biradical mechanism seems to fit these results. The heat of formation and the entropy of COD are estimated.     

Eclipsed and staggered conformations of (SIH3)2F+: An ab initio study

Ab initio calculations at 6–31G**, 6–31++G**, and MP2/6–31G** levels were performed on disilyl–fluoronium, (SiH₃)₂F⁺, with the SiH₃ group eclipsed or staggered. Optimized geometries, total energies, dipole moments, atomic charges, electronic density, and vibrational frequencies were computed. The results were compared with calculated structural parameters and vibrational frequencies of H₃SiF, H₂SiF⁺, H₂SiF^?, and H₄SiF⁺ ions. The basis-set effects were studied. Several thermochemistry parameters—ZPE, thermal energy, rotational constants, and entropies—were also calculated. © 1994 John Wiley & Sons, Inc.     

Mechanism of thermal reactions of dimerization and trimerization of hexafluoro-1,3-butadiene: A quantum-chemical study

The mechanism of thermal oligomerization of hexafluoro-l,3-butadiene was examined using RHF, ROHF, and GVB/DH, and by B3LYP/6-31G* and /6-311G* quantum-chemical methods. The energies of highly reactive excited states of the monomer and of intermediate biradical species were estimated. Isomeric biradicals (excited monomers and dimers) with the free valences localized on two carbon atoms were shown to coexist in the reaction mixture. Recombination of such biradicals and their reactions with neutral stable molecules give rise to diverse products, depending on the reaction temperature.     

PAH formation under single collision conditions: reaction of phenyl radical and 1,3-butadiene to form 1,4-dihydronaphthalene

The crossed beam reactions of the phenyl radical (C(6)H(5), X(2)A(1)) with 1,3-butadiene (C(4)H(6), X(1)A(g)) and D6-1,3-butadiene (C(4)D(6), X(1)A(g)) as well as of the D5-phenyl radical (C(6)D(5), X(2)A(1)) with 2,3-D2-1,3-butadiene and 1,1,4,4-D4-1,3-butadiene were carried out under single collision conditions at collision energies of about 55 kJ mol(-1). Experimentally, the bicyclic 1,4-dihydronaphthalene molecule was identified as a major product of this reaction (58 ± 15%) with the 1-phenyl-1,3-butadiene contributing 34 ± 10%. The reaction is initiated by a barrierless addition of the phenyl radical to the terminal carbon atom of the 1,3-butadiene (C1/C4) to form a bound intermediate; the latter underwent hydrogen elimination from the terminal CH(2) group of the 1,3-butadiene molecule leading to 1-phenyl-trans-1,3-butadiene through a submerged barrier. The dominant product, 1,4-dihydronaphthalene, is formed via an isomerization of the adduct by ring closure and emission of the hydrogen atom from the phenyl moiety at the bridging carbon atom through a tight exit transition state located about 31 kJ mol(-1) above the separated products. The hydrogen atom was found to leave the decomposing complex almost parallel to the total angular momentum vector and perpendicularly to the rotation plane of the decomposing intermediate. The defacto barrierless formation of the 1,4-dihydronaphthalene molecule involving a single collision between a phenyl radical and 1,3-butadiene represents an important step in the formation of polycyclic aromatic hydrocarbons (PAHs) and their partially hydrogenated counterparts in combustion and interstellar chemistry.     

Argon-matrix-isolation Raman spectra and density functional study of 1,3-butadiene conformers

s-trans, s-cis and gauche conformers of 1,3-butadiene have been studied using density functional theory and the coupled-cluster method using double substitutions (CCD). Matrix isolation Raman and IR data for the minor conformer were obtained and are used in combination with the theoretical results to resolve earlier ambiguities in vibrational assignments. Based on high-quality Hessians, new harmonic stretching force constants are reported for the carbon backbone of s-trans-1,3-butadiene. For the minor conformer the best unscaled root mean square error of the calculated frequencies for the s-cis and gauche geometries are 17.5 cm⁻¹ and 7.4 cm⁻¹, respectively, primarily due to a better agreement of the gauche results for the vibrations at 983 cm⁻¹, 596 cm⁻¹ and 470 cm⁻¹ which depend strongly on the torsional angle. Although this points towards the gauche form rather than the s-cis form, the calculated transition dipole moment directions at the CCD/6-311G(d,p) level confirm the earlier conclusion that the minor conformer has C _{2 v} symmetry in the matrix. It is concluded that either the better agreement between the frequencies calculated for the gauche form and the observed values is coincidental, or that the molecule is indeed nonplanar in the matrix and tunnels very rapidly between the two mirror-image forms (or its lowest vibrational level lies above the barrier). Received: 1 July 1998 / Accepted: 26 October 1998 / Published online: 15 February 1999     

Ab initio study of possible intermediates in the no(4Π) catalyzed geometric isomerizations of olefins

The geometries of acyclic and three-membered ring (nitroxide) H₄C₂NO radicals in their ground ²Π electronic states have been optimized completely at ab initio UHF and ROHF theoretical levels with the STO-3G and the 6-31G** basis sets. The optimizations favour the cyclic nitroxide structure energetically. However ΔE(acyclic - cyclic) at the UHF and ROHF/6-31G** levels are only 3.2 and 1.9 kcal mol^-1, respectively. Incomplete MP2/6-311G** optimizations support these results. The zero-point energy computed at the ROHF/6-31G** level for the nitroxide radical is 2.5 kcal mol^-1 higher than that for the acyclic structure, thus reversing the relative energies by 0.6 kcal mol^-1. The energies of the two radical structures, relative to the sum of those for ethylene and NO, are very close to literature values of the activation energies for the thermal, NO catalyzed geometrical isomerizations of olefins. Thus cyclic nitroxide intermediates may play a role not only in the Hg 6(³P₁) photosensitized, but also in the thermal, NO-catalyzed geometric isomerizations of olefins. Paper dedicated to Professor Otto P. Strausz; presented in part at the 75th Canadian Chemical Congress and Exhibition, Edmonton, May 31 – June 4, 1992.     

A molecular orbital study of the rotation about the C-C bond in 1,3-butadiene

The geometry and energy of 1,3-butadiene have been calculated using the 6-311G** basis set as a function of the CCCC dihedral angle-0 ° (trans), 30 °, 60 °, 75 °, 90 °, 120 °, 135 °, 150 °, 165 ° and 180 ° (cis)-assuming that the vinyl groups remain planar. Potential minima are located at 0 ° and 141.4 °, with the trans structure more stable than the gauche by 13.2 kJ mol^–1. Potential maxima are located at 76.7 °, giving a barrier height of 25.4 kJ mol^–1 relative to the trans structure, and at 180 ° giving a barrier height of 3.0 kJ mol^–1 relative to the 141.4 °-gauche structure. Using the 6-31G* basis set the inclusion of electron correlation, accounting for about 52% of the correlation energy, was found to produce no significant change in the shape of the potential energy curve. The magnitude of the expectation energy differences is such that both barriers with respect to the 14l.4 °-gauche maximum structure can be categorized unequivocally as attractive-dominant, whereas the values for the energy barrier with respect to the trans structure, although characteristic of a repulsive-dominant barrier at the 6–311G** level, are sufficiently small that higher level calculations might give the opposite result. Analysis of V _nn for the conversion reactions cis 150 °-gauche, trans 60 °-gauche, and trans 90 °-gauche in terms of the individual contributions from the various internuclear interactions shows that nonbonded interactions are important, not only in initiating the destabilization of the crowded cis structure, but also through-out the entire range of CCCC dihedral angles, 0 ° to 180 °.