Theoretical study on the unimolecular decomposition of thiophenol

The potential energy surface for the unimolecular decomposition of thiophenol (C(6)H(5)SH) is mapped out at two theoretical levels; BB1K/GTlarge and QCISD(T)/6-311+G(2d,p)//MP2/6-31G(d,p). Calculated reaction rate constants at the high pressure limit indicate that the major initial channel is the formation of C(6)H(6)S at all temperatures. Above 1000 K, the contribution from direct fission of the S-H bond becomes important. Other decomposition channels, including expulsion of H(2) and H(2)S are of negligible importance. The formation of C(6)H(6)S is predicted to be strong-pressure dependent above 900 K. Further decomposition of C(6)H(6)S produces CS and C(5)H(6). Overall, despite the significant difference in bond dissociation, i.e., 8-9 kcal/mol between the S-H bond in thiophenol and the O-H bond in phenol, H migration at the ortho position in the two molecules represents the most accessible initial channel.     

Energy barriers to gas-phase unimolecular decomposition of mono- and dinitrotoluenes

Alternative paths of gas-phase unimolecular decomposition of three nitrotoluenes and six dinitrotoluenes, in particular homolytic dissociation of the C–N bond, nitro–nitrite rearrangement, intramolecular hydrogen transfer from the methyl to nitro group with formation of isomeric aci-nitrotoluenes, and various paths involving formation of bicyclic intermediates, have been studied at the DFT B3LYP/6-31+G(2df,p) level of theory using GAUSSIAN 09 software package. The most energetically favorable path for o-nitrotoluene and 2,3-, 2,4-, 2,5-, and 2,6-dinitrotoluenes is the formation of aci-nitrotoluenes. The effect of the substrate structure on the competition between different mechanisms of these reactions has been analyzed.     

Theoretical studies of the NTO unimolecular decomposition

This work studies 39 decomposition paths among 18 intermediates and 14 transition states. Three types of intra-molecular proton migration and the direct scission of C–NO₂ were regarded as the initial steps of the unimolecular decomposition of NTO. The activation energies of the radicalization C–NO₂ homolysis step are 79.158, 79.781 and 80.652 kcal mol⁻¹. The activation energies of the ionization C–NO⁻¹₂ scission step are 262.488, 263.138 and 272.278 kcal mol⁻¹. The bottle neck activation energies of the C–NO₂H cleavage are 54.936, 63.257 and 71.247 kcal mol⁻¹. Two paths have the smallest bottle neck activation energy. Both of them have two proton migration steps and one internal rotation step prior to C–NO₂H cleavage. At lower temperatures, energy accumulated slowly. When the energy is high enough and reaction time is long enough for structure transformation, these two mechanisms should be the most probable decomposition paths. At high temperatures, the shortest (four steps) mechanism which goes through radicalization C–NO₂ scission should be the dominant path. There are five tautomers found in this study. Four of them are intra-molecular proton migration tautomers. The other one is an internal rotational tautomer. Their energy barriers for structure transfer are lower than any of the activation energies of the decomposition reactions. It may be regarded as one explanation of the insensitive property of NTO.     

Theoretical kinetic study of thermal unimolecular decomposition of cyclic alkyl radicals

Whereas many studies have been reported on the reactions of aliphatic hydrocarbons, the chemistry of cyclic hydrocarbons has not been explored extensively. In the present work, a theoretical study of the gas-phase unimolecular decomposition of cyclic alkyl radicals was performed by means of quantum chemical calculations at the CBS-QB3 level of theory. Energy barriers and high-pressure-limit rate constants were calculated systematically. Thermochemical data were obtained from isodesmic reactions, and the contribution of hindered rotors was taken into account. Classical transition state theory was used to calculate rate constants. The effect of tunneling was taken into account in the case of CH bond breaking. Three-parameter Arrhenius expressions were derived in the temperature range of 500-2000 K at atmospheric pressure, and the CC and CH bond breaking reactions were studied for cyclic alkyl radicals with a ring size ranging from three to seven carbon atoms, with and without a lateral alkyl chain. For the ring-opening reactions, the results clearly show an increase of the activation energy as the pi bond is being formed in the ring (endo ring opening) in contrast to the cases in which the pi bond is formed on the side chain (exo ring opening). These results are supported by analyses of the electronic charge density that were performed with Atoms in Molecules (AIM) theory. For all cycloalkyl radicals considered, CH bond breaking exhibits larger activation energies than CC bond breaking, except for cyclopentyl for which the ring-opening and H-loss reactions are competitive over the range of temperatures studied. The theoretical results compare rather well with the experimental data available in the literature. Evans-Polanyi correlations for CC and CH beta-scissions in alkyl and cycloalkyl free radicals were derived. The results highlight two different types of behavior depending on the strain energy in the reactant.     

Theoretical study of the unimolecular decomposition mechanisms of energetic TNAD and TNAZ explosives

Calculation methods based on hybrid Density Functional Theory (DFT) with the basis sets of the B3LYP/6‐31+G(d)//B3LYP/4‐31G(d) method and the differential overlap (INDO) program were used to derive reasonable decomposition mechanisms of 1,4,5,8‐tetranitro‐1,4,5,8‐tetraazadecalin (TNAD) and 1,3,3‐trinitroazetidine (TNAZ) explosives. All possible decomposition species and transition states, including intermediates and products, were identified and their corresponding enthalpy of formation and Gibbs free energy of formation were obtained using polyparametric modification equations. INDO bond energy calculation results reveal the weakest bonding site for reference and determine where cleavage can occur easily. This work is concerned mainly with eliminating HONO (cis or trans form). The activation energy for trans‐form HONO elimination is lower than that of cis‐form HONO elimination in the initial steps of both TNAD and TNAZ decomposition, being 18.5 kJ/mol and 33.3 kJ/mol, respectively. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical studies on the unimolecular decomposition of propanediols and glycerol

Polyols, a typical type of alcohol containing multiple hydroxyl groups, are being regarded as a new generation of a green energy platform. In this paper, the decomposition mechanisms for three polyol molecules, i.e., 1,2-propanediol, 1,3-propanediol, and glycerol, have been investigated by quantum chemistry calculations. The potential energy surfaces of propanediols and glycerol have been built by the QCISD(T) and CBS-QB3 methods, respectively. For the three molecules studied, the H(2)O-elimination and C-C bond dissociation reactions show great importance among all of the unimolecular decomposition channels. Rate constant calculations further demonstrate that the H(2)O-elimination reactions are predominant at low temperature and pressure, whereas the direct C-C bond dissociation reactions prevail at high temperature and pressure. The temperature and pressure dependence of calculated rate constants was demonstrated by the fitted Arrhenius equations. This work aims to better understand the thermal decomposition process of polyols and provide useful thermochemical and kinetic data for kinetic modeling of polyols-derived fuel combustion.     

Theoretical investigation of unimolecular decomposition channels of furan4

Density functional theory and high-level ab initio calculations were carried out to investigate three unimolecular decomposition channels of furan. All equilibrium and transition state structures along the proposed decomposition channels are fully optimized by B3LYP/6-31G** and characterized at the same level of theory by vibrational and intrinsic reaction coordinate analyses. Relative energies of the optimized structures were evaluated at theoretical levels up to QCISD(T)/6-311++G**. The theoretical results suggest that the unimolecular decomposition channel of isoxazole, proposed in an experimental study and implied to be the main decomposition channel of furan, is responsible only for the formation of HC(TRIPLE BOND)CH and H₂O(DOUBLE BOND)C(DOUBLE BOND)O, minor products of furan thermal decomposition. A new decomposition mechanism, proposed in the present study, is shown to be more likely responsible for the formation of CH₃C(TRIPLE BOND)CH and CO, major products of furan thermal decomposition. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 240–249, 1998     

H2ClCS单分子分解反应机理的理论研究

用量子化学B3LYP/6 - 311+G(d,p)方法优化了H2ClCS单分子分解反应驻点物种的几何构型,并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证.用QCISD(T)/6-311++G(d,p)方法计算各物种的单点能,并对总能量进行了零点能校正.利用经典过渡态理论(TST)与...     

Theoretical studies on the unimolecular decomposition of ethylene glycol

The unimolecular decomposition processes of ethylene glycol have been investigated with the QCISD(T) method with geometries optimized at the B3LYP/6-311++G(d,p) level. Among the decomposition channels identified, the H(2)O-elimination channels have the lowest barriers, and the C-C bond dissociation is the lowest-energy dissociation channel among the barrierless reactions (the direct bond cleavage reactions). The temperature and pressure dependent rate constant calculations show that the H(2)O-elimination reactions are predominant at low temperature, whereas at high temperature, the direct C-C bond dissociation reaction is dominant. At 1 atm, in the temperature range 500-2000 K, the calculated rate constant is expressed to be 7.63 × 10(47)T(-10.38) exp(-42262/T) for the channel CH(2)OHCH(2)OH → CH(2)CHOH + H(2)O, and 2.48 × 10(51)T(-11.58) exp(-43593/T) for the channel CH(2)OHCH(2)OH → CH(3)CHO + H(2)O, whereas for the direct bond dissociation reaction CH(2)OHCH(2)OH → CH(2)OH + CH(2)OH the rate constant expression is 1.04 × 10(71)T(-16.16) exp(-52414/T).     

Theoretical calculation of the unimolecular decomposition of the fluorine fluorosulfonate molecule

The RRKM theory is applied to the unimolecular decomposition of fluorine fluorosulfonate. Theoretical results are compared with the available experimental values for different temperatures and a good agreement is found.

---

PPKM . ; .

     

Theoretical determinations of the ambient conformational distribution and unimolecular decomposition of n-propylperoxy radical

The conformational distribution and unimolecular decomposition pathways for the n-propylperoxy radical have been generated at the CBS-QB3, B3LYP/6-31+G and mPW1K/6-31+G levels of theory. At each of the theoretical levels, the 298 K Boltzmann distributions and rotational profiles indicate that all five unique rotamers of the n-propylperoxy radical can be expected to be present in significant concentrations at thermal equilibrium. At the CBS-QB3 level, the 298 K distribution of rotamers is predicted to be 28.1, 26.4, 19.6, 14.0, and 11.9% for the gG, tG, gT, gG', and tT conformations, respectively. The CBS-QB3 C-OO bond dissociation energy (DeltaH298 K) for the n-propylperoxy radical has been calculated to be 36.1 kcal/mol. The detailed CBS-QB3 potential energy surface for the unimolecular decomposition of the n-propylperoxy radical indicates that important bimolecular products could be derived from two 1,4-H transfer mechanisms available at T < 500 K, primarily via an activated n-propylperoxy adduct.     

Nucleophilic reactions of charge delocalised carotenoid mono- and dications

In the present study insight was gained on the larger complexity of cationic mixtures of diaryl (phi,phi-carotene, isorenieratene) and aliphatic (psi,psi-carotene, lycopene) carotenes, prepared by reaction with BF3-etherate, compared with beta,beta-carotene. Chemical reactions of the mono- and dications prepared in situ from the allylic carotenols beta,beta-caroten-4-ol (isocryptoxanthin) and beta,beta-carotene-4,4'-diol (isozeaxanthin), and from isorenieratene and lycopene were investigated using selected O, N and S nucleophiles; water, methanol, azide and thioacetate. In total 22, including 18 new, neutral carotenoid products were isolated and identified by VIS, MS and NMR (in part) spectroscopy. Their structures were compatible with the structures of the cationic intermediates. The formal addition of hydride to the various dications, required to rationalise minor reaction products, is discussed in terms of more likely hydrogen radical or proton transfer in cationic reactions. Extensive E/Z isomerisation was observed for all quenching products. The potential use of carotenoid cations for the synthesis of 4,(4')-substituted beta,beta-carotenes and 7-oxabicyclo[2,2,1]heptane derivatives is discussed.     

Electron capture dissociation and infrared multiphoton dissociation of oligodeoxynucleotide dications

We report electron capture dissociation (ECD) and infrared multiphoton dissociation (IRMPD) of doubly protonated and protonated/alkali metal ionized oligodeoxynucleotides. Mass spectra following ECD of the homodeoxynucleotides polydC, polydG, and polydA contain w or d "sequence" ions. For polydC and polydA, the observed fragments are even-electron ions, whereas radical w/d ions are observed for polydG. Base loss is seen for polydG and polydA but is a minor fragmentation pathway in ECD of polydC. We also observe fragment ions corresponding to w/d plus water in the spectra of polydC and d(GCATGC). Although the structure of these ions is not clear, they are suggested to proceed through a pentavalent phosphorane intermediate. The major fragment in ECD of d(GCATGC) is a d ion. Radical a- or z-type fragment ions are observed in most cases. IRMPD primarily results in base loss, but backbone fragmentation is also observed. IRMPD provides more sequence information than ECD, but the spectra are more complex due to extensive base and water losses. It is proposed that the smaller degree of sequence coverage in ECD, with fragmentation mostly occurring close to the ends of the molecules, is a consequence of a mechanism in which the electron is captured at a P=O bond, resulting in a negatively charged phosphate group. Consequently, at least two protons (or alkali metal cations) must be present to observe a w or d fragment ion, a requirement that is less likely for small fragments.     

A theoretical study of the unimolecular decomposition of difluorochloromethane

Nonempirical quantum-chemical calculations have been performed in the 4-31G basis on the equilibrium geometry and vibrational frequencies of difluorochloromethane, as well as on the transition state in the decomposition CHF₂ClCF₂+HCl. The potentialenergy surface has been constructed and a simplified 1D dynamic model has been formulated for the reaction, which is based on dividing the degrees of freedom into active (dynamic) and adiabatic ones. Measurements on the reaction kinetics have been analyzed from calculations on the rate constant from the RRKM theory and within the framework of the 1D model.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 23, No. 3, pp. 274–281, May–June 1987.I am indebted to M. V. Bazilevskii, M. Ya. Gol'denberg, S. Ya. Umanskii, A. I. Voronin, and M. A. Teitel'boim for useful discussions.     

Competition between multiphoton fluorescence and multiphoton ionization in NO

A nitrogen laser pumped tunable dye laser has been used to observe the three-photon ionization of NO through a two-photon resonance with the C²II state. Fluorescence is also observed from this state. The wavelength dependence of both signals have been measured. A reaction mechanism is postulated, which includes the initial two-photon excitation of the C²II state as the rate-limiting step. This mechanism predicts the observed second-order intensity dependence of the ionization signal and shows that the simple rate equation treatment is valid in this system.     

Theoretical study of the thermal unimolecular rearrangement of fluoroethylidenes

The thermal unimolecular isomerization of fluoroethylidenes to the corresponding fluoroethylenes has been studied by the MNDO method. It has been shown that fluorine substitution on the carbene carbon increases the activation energy in comparison with the ethylidene rearrangement. To understand the reason for this increase in the activation energy, the charge-transfer effects have been analyzed. Fluorine substitution at other positions does not significantly affect the activation energies. The thermodynamic parameters for the reaction have been evaluated, using vibrational and rotational spectral data calculated in this work. RRKM calculations have been performed and high-pressure Arrhenius parameters calculated. Hydrogen–deuterium kinetic isotope effects indicate that the reaction rates are altered considerably on isotopic substitution, and the change in reaction rates depends upon the position of deuterium substitution, as well as on the number of hydrogens replaced by deuterium atoms. © 1992 John Wiley & Sons, Inc.     

Theoretical study of mono- and dithiocarbazic acid derivatives

The patterns of charge distribution in a series of derivatives of thiocarbazic acid are studied by means of several theoretical LCAO-MO methods: ab initio minimal basis-set, CNDO/2, INDO and EHT. Correlations are made between the electronic structures obtained from the calculations, the nature of the different substituents and some experimental IR, X-ray and XPS results.     

Theoretical study of unimolecular rearrangements of vinylidenes to acetylenes

The rearrangement of vinylidene to acetylene has been studied in detail by the density functional method, using Becke's three‐parameter exchange functional and the gradient‐corrected functional of Lee, Yang, and Parr. The rearrangement of the anion, as well as that of fluoro‐substituted systems, has also been investigated, in order to determine the effect of fluorine substitution on the activation barrier to the 1,2‐hydrogen shift, as well as the relative migratory aptitudes of hydrogen and fluorine. Natural bond orbital analysis is invoked to gain insight into the mechanisms of the rearrangements. Basis size effects are also discussed, particularly in relation to anionic systems. The need to include diffuse functions in geometry optimizations of anionic systems is reinforced by the present calculations. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Effective Hamiltonian models and unimolecular decomposition

Partitioning Hilbert space into two subspaces by using orthogonal projection operators yields compact forms for effective Hamiltonians for each of the subspaces. When one (the Q space) contains molecular bound states and the other (the P space) contains dissociative continua, a simple form for the non-Hermitian Q-space effective Hamiltonian, H(eff), can be obtained, subject to reasonable approximations. Namely, H(eff) = H0 - ivariant Planck's/2pi Gamma/2, where H0 is Hermitian, and the width operator variant Planck's/2pi Gamma accounts for couplings of the Q-space levels to the P-space continua. The P/Q partitioning procedure has been applied in many areas of atomic, molecular, and nuclear physics with widespread success. Inputting into this formalism ideas from random matrix theory in order to model independent open channels yields the random matrix H(eff) model. Despite numerous efforts, this model has failed to model satisfactorily the statistical transition-state theory of unimolecular decomposition (hereafter referred to as TST) in the regime of overlapping resonances, where nearly all such reactions occur. All statistical models of unimolecular decomposition are premised on rapid intramolecular vibrational redistribution (IVR) for a given set of good quantum numbers. The phase space thus accessed results in a threshold reaction rate of 1/h rho, and for K independent open channels, the rate is K/h rho. This reaction rate corresponds to a resonance width of K/2pi rho, and when K increases, the resonances (which are rho(-1) apart) overlap. In this regime, the random matrix H(eff) model fails because it does not introduce independent open channels. To illustrate the source of the problem, an analysis is carried out of a simple model that is obviously and manifestly inconsistent with TST. This model is solved exactly, and it is then put in the form of the random matrix H(eff) model, illustrating the one-to-one correspondence. This reveals the deficiencies of the latter. In manipulating this model into the form H0 - ivariant Planck's/2pi Gamma/2, it becomes clear that the independent open channels in the random matrix H(eff) model are inconsistent with TST. Rather, this model is one of gateway states (i.e., bound states that are coupled to their respective continua as well as to a manifold of zero-order bound states, none of which are coupled directly to the continua). Despite the fact that the effective Hamiltonian method is, by itself, beyond reproach, the random matrix H(eff) model is flawed as a model of unimolecular decomposition in several respects, most notably, bifurcations of the distributions of resonance widths in the regime of overlapping resonances.     

Theoretical study of the structure and unimolecular decomposition pathways of ethyloxonium, [CH3CH2OH2]+

Ab initio molecular orbital calculations with split-valence plus polarization basis sets and incorporating valence-electron correlation have been performed to determine the equilibrium structure of ethyloxonium ([CH₃CH₂OH₂]⁺) and examine its modes of unimolecular dissociation. An asymmetric structure (1) is predicted to be the most stable form of ethyloxonium, but a second conformational isomer of C_s symmetry lies only 1.4 kJ mol^?1 higher in energy than 1. Four unimolecular decomposition pathways for 1 have been examined involving loss of H₂, CH₄, H₂O or C₂H₄. The most stable fragmentation products, lying 65 kJ mol^?1 above 1, are associated with the H₂ elimination reaction. However, large barriers of 257 and 223 kJ mol^?1 have to be surmounted for H₂ and CH₄ loss, respectively. On the other hand, elimination of either C₂H₄ or H₂O from ethyloxonium can proceed without a barrier to the reverse associations and, with total endothermicities of 130 and 160 kJ mol^?1, respectively, these reactions are expected to dominate at lower energies. A second important equilibrium structure on the surface is a hydrogen-bridged complex, lying 53 kJ mol^?1 above 1. This complex is involved in the C₂H₄ elimination reaction, acts as an intermediate in the proton-transfer reaction connecting [C₂H₅]⁺ +H₂O and C₂H₄ + [H₃O]⁺ and plays an important role in the isotopic scrambling that has been observed experimentally in the elimination of either H₂O or C₂H₄ from ethyloxonium. The proton affinity of ethanol was calculated as 799 kJ mol^?1, in close agreement with the experimental value of 794 kJ mol^?1.