Interaction of petroleum-relevant organosulfur compounds with TiO2(110)

The interaction of two sets of structurally related molecules, thiophenol/thioanisole, and thiophene/tetrahydrothiophene, with vacuum-annealed and ion-bombarded TiO(2)(110) surfaces has been studied using a combination of temperature-programmed reaction spectroscopy (TPRS) and X-ray photoelectron spectroscopy (XPS). All thioethers studied were observed to adsorb and desorb from both surfaces without producing reaction products, while thiophenol, the only species studied containing a S-H bond, reacted with both surfaces. Approximately 25% of surface bound thiophenol decomposed over the vacuum-annealed surface. On the bombarded surface, thiophenol both decomposed into surface-bound C(x)H(y)/S fragments, and reacted to form benzene, which desorbed from the surface at 400 K. We propose that phenylthiolate formation on the bombarded surface leads to the observed production of benzene. These results highlight the importance of defects in the reactivity of titania, and lay the foundation for the study of larger, refractory sulfur compounds present in fuel.     

Thermochemical and kinetic analysis of the thermal decomposition of monomethylhydrazine: an elementary reaction mechanism

The reaction kinetics for the thermal decomposition of monomethylhydrazine (MMH) was studied with quantum Rice-Ramsperger-Kassel (QRRK) theory and a master equation analysis for pressure falloff. Thermochemical properties were determined by ab initio and density functional calculations. The entropies, S degrees (298.15 K), and heat capacities, Cp degrees (T) (0 < or = T/K < or = 1500), from vibrational, translational, and external rotational contributions were calculated using statistical mechanics based on the vibrational frequencies and structures obtained from the density functional study. Potential barriers for internal rotations were calculated at the B3LYP/6-311G(d,p) level, and hindered rotational contributions to S degrees (298.15 K) and Cp degrees (T) were calculated by solving the Schr?dinger equation with free rotor wave functions, and the partition coefficients were treated by direct integration over energy levels of the internal rotation potentials. Enthalpies of formation, DeltafH degrees (298.15 K), for the parent MMH (CH3NHNH2) and its corresponding radicals CH3N*NH2, CH3NHN*H, and C*H2NHNH2 were determined to be 21.6, 48.5, 51.1, and 62.8 kcal mol(-1) by use of isodesmic reaction analysis and various ab initio methods. The kinetic analysis of the thermal decomposition, abstraction, and substitution reactions of MMH was performed at the CBS-QB3 level, with those of N-N and C-N bond scissions determined by high level CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) calculations. Rate constants of thermally activated MMH to dissociation products were calculated as functions of pressure and temperature. An elementary reaction mechanism based on the calculated rate constants, thermochemical properties, and literature data was developed to model the experimental data on the overall MMH thermal decomposition rate. The reactions of N-N and C-N bond scission were found to be the major reaction paths for the modeling of MMH homogeneous decomposition at atmospheric conditions.     

A computational study on the kinetics and mechanism for the unimolecular decomposition of o-nitrotoluene

The kinetics and mechanism for the unimolecular decomposition of o-nitrotoluene (o-CH(3)C(6)H(4)NO(2)) have been studied computationally at the G2M(RCC, MP2)//B3LYP/6-311G(d, p) level of theory in conjunction with rate constant predictions with RRKM and TST calculations. The results of the calculations reveal 10 decomposition channels for o-nitrotoluene and its six isomeric intermediates, among them four channels give major products: CH(3)C(6)H(4) + NO(2), C(6)H(4)C(H)ON (anthranil) + H(2)O, CH(3)C(6)H(4)O (o-methyl phenoxy) + NO, and C(6)H(4)C(H(2))NO + OH. The predicted rate constants in the 500-2000 K temperature range indicate that anthranil production, taking place initially by intramolecular H-abstraction from the CH(3) group by NO(2) followed by five-membered ring formation and dehydration, dominates at temperatures below 1000 K, whereas NO(2) elimination becomes predominant above 1100 K and CH(3)C(6)H(4)O formation by the nitro-nitrite isomerization/decomposition process accounts for only 5-11% of the total product yield in the middle temperature range 800-1300 K. The branching ratio for CH(2)C(6)H(4)NO formation by the decomposition process of CH(2)C(6)H(4)N(O)OH is negligible. The predicted high-pressure-limit rate constants with the rate expression of 4.10 x 10(17) exp[-37000/T] s(-1) for the NO(2) elimination channel and 9.09 x 10(12) exp[-25800/T] s(-1) for the H(2)O elimination channel generally agree reasonably with available experimental data. The predicted high-pressure-limit rate constants for the NO and OH elimination channels are represented as 1.49 x 10(14) exp[-30000/T] and 1.31 x 10(15) exp[-38000/T] s(-1), respectively.     

Computational study on the mechanism and rate constant for the C6H5 + C6H5NO reaction

The reaction mechanism of C6H5 + C6H5NO involving four product channels on the doublet-state potential energy surface has been studied at the B3LYP/6-31+G(d, p) level of theory. The first reaction channel occurs by barrierless association forming (C6H5)2NO (biphenyl nitroxide), which can undergo isomerization and decomposition. The second channel takes place by substitution reaction producing C12H10 (biphenyl) and NO. The third and fourth channels involve direct hydrogen abstraction reactions producing C6H4NO + C6H6 and C6H5NOH + C6H4, respectively. Bimolecular rate constants of the above four product channels have been calculated in the temperature range 300-2000 K by the microcanonical Rice-Ramsperger-Kassel-Marcus theory and/or variational transition-state theory. The result shows the dominant reactions are channel 1 at lower temperatures (T < 800 K) and channel 3 at higher temperatures (T > 800 K). The total rate constant at 7 Torr He is predicted to be k(t) = 3.94 x 10(21) T(-3.09) exp(-699/T) for 300-500 K, 2.09 x 10(20) T(-3.56) exp(2315/T) for 500-1000 K, and 1.51 x 10(2) T(3.30) exp(-3043/T) for 1000-2000 K (in units of cm3 mol(-1) s(-1)), agreeing reasonably with the experimental data within their reported errors. The heats of formation of key products including biphenyl nitroxide, hydroxyl phenyl amino radical, and N-hydroxyl carbazole have been estimated.     

Near-ultraviolet photodissociation of thiophenol

H(D) Rydberg atom photofragment translational spectroscopy has been used to investigate the dynamics of H(D) atom loss C6H5SH(C6H5SD) following excitation at many wavelengths lambda phot in the range of 225-290 nm. The C6H5S cofragments are formed in both their ground (X(2)B1) and first excited ((2)B2) electronic states, in a distribution of vibrational levels that spreads and shifts to higher internal energies as lambda(phot) is reduced. Excitation at lambda(phot) > 275 nm populates levels of the first (1)pi pi* state, which decay by tunnelling to the dissociative (1)pi sigma* state potential energy surface (PES). S-H torsional motion is identified as a coupling mode facilitating population transfer at the conical intersection (CI) between the diabatic (1)pi pi* and (1)pi sigma* PESs. At shorter lambda(phot), the (1)pi sigma* state is deduced to be populated either directly or by efficient vibronic coupling from higher (1)pipi* states. Flux evolving on the (1)pi sigma* PES samples a second CI, at longer R(S-H), between the diabatic (1)pi sigma* and ground ((1)pi pi) PESs, where the electronic branching between ground and excited state C6H5S fragments is determined. The C6H5S(X(2)B1) and C6H5S((2)B2) products are deduced to be formed in levels with, respectively, a' and a' vibrational symmetry-behavior that reflects both Franck-Condon effects (both in the initial photoexcitation step and in the subsequent in-plane forces acting during dissociation) and the effects of the out-of-plane coupling mode(s), nu11 and nu16a, at the (1)pi sigma*/(1)pi pi CI. The vibrational state assignments enabled by the high-energy resolution of the present data allow new and improved estimations of the bond dissociation energies, D0(C6H5S-H) < or = 28,030 +/- 100 cm(-1) and D0(C6H5S-D) < or = 28,610 +/- 100 cm(-1), and of the energy separation between the X(2)B1 and (2)B2 states of the C6H5S radical, T(00) = 2800 +/- 40 cm(-1). Similarities, and differences, between the measured energy disposals accompanying UV photoinduced X-H (X = S, O) bond fission in thiophenol and phenol are discussed.     

Experimental and theoretical study of the photodissociation reaction of thiophenol at 243 nm: intramolecular orbital alignment of the phenylthiyl radical

The photoinduced hydrogen (or deuterium) detachment reaction of thiophenol (C(6)H(5)SH) or thiophenol-d(1) (C(6)H(5)SD) pumped at 243 nm has been investigated using the H (D) ion velocity map imaging technique. Photodissociation products, corresponding to the two distinct and anisotropic rings observed in the H (or D) ion images, are identified as the two lowest electronic states of phenylthiyl radical (C(6)H(5)S). Ab initio calculations show that the singly occupied molecular orbital of the phenylthiyl radical is localized on the sulfur atom and it is oriented either perpendicular or parallel to the molecular plane for the ground (B(1)) and the first excited state (B(2)) species, respectively. The experimental energy separation between these two states is 2600+/-200 cm(-1) in excellent agreement with the authors' theoretical prediction of 2674 cm(-1) at the CASPT2 level. The experimental anisotropy parameter (beta) of -1.0+/-0.05 at the large translational energy of D from the C(6)H(5)SD dissociation indicates that the transition dipole moment associated with this optical transition at 243 nm is perpendicular to the dissociating S-D bond, which in turn suggests an ultrafast D+C(6)H(5)S(B(1)) dissociation channel on a repulsive potential energy surface. The reduced anisotropy parameter of -0.76+/-0.04 observed at the smaller translational energy of D suggests that the D+C(6)H(5)S(B(2)) channel may proceed on adiabatic reaction paths resulting from the coupling of the initially excited state to other low-lying electronic states encountered along the reaction coordinate. Detailed high level ab initio calculations adopting multireference wave functions reveal that the C(6)H(5)S(B(1)) channel may be directly accessed via a (1)(n(pi),sigma(*)) photoexcitation at 243 nm while the key feature of the photodissociation dynamics of the C(6)H(5)S(B(2)) channel is the involvement of the (3)(n(pi),pi(*))-->(3)(n(sigma),sigma(*)) profile as well as the spin-orbit induced avoided crossing between the ground and the (3)(n(pi),sigma(*)) state. The S-D bond dissociation energy of thiophenol-d(1) is accurately estimated to be D(0)=79.6+/-0.3 kcalmol. The S-H bond dissociation energy is also estimated to give D(0)=76.8+/-0.3 kcalmol, which is smaller than previously reported ones by at least 2 kcalmol. The C-H bond of the benzene moiety is found to give rise to the H fragment. Ring opening reactions induced by the pi-pi(*)n(pi)-pi(*) transitions followed by internal conversion may be responsible for the isotropic broad translational energy distribution of fragments.     

Pyrolysis of tert-butyl tert-butanethiosulfinate, t-BuS(O)St-Bu: a computational perspective of the decomposition pathways

A systematic theoretical study has been performed on the low pressure thermal decomposition pathways of t-BuS(O)St-Bu using the CCSD(T)/cc-pV(D+d)Z//B3LYP/6-311++G(2d,2p), CCSD(T)/cc-pV(D+d)Z//PBEPBE/6-311++G(2d,2p), and G3B3 level of theories. Rate constants for the unimolecular decomposition pathways are calculated using Rice?Ramsperger?Kassel?Marcus (RRKM) theory. On the basis of the experimental observation and theoretical predictions, the pyrolysis channels are considered as primary and secondary pyrolysis reactions. The primary decomposition via a five-membered transition state leads to the formation of tert-butanethiosulfoxylic acid (t-BuSSOH) and 2-methylpropene (C4H8) almost exclusively having low-pressure limit rate constant k(1)(0) = 4.67 × 10(?6)T(?4.67) exp(?11.64 kcal mol(?1)/RT) cm3 mol(?1) s(?1) (T = 500?800 K). The primary decomposition via a six-membered transition state is also identified, and that leads to the tert-butanethiosulfinic acid t-BuS(OH)S, which is the branched chain isomer of t-BuSSOH. The formation of t-BuSSOH is thermodynamically as well as kinetically favorable over t-BuS(OH)S formation, and therefore the second product could not be found experimentally. Furthermore, calculation on secondary pyrolysis pathways involving the decomposition of t-BuSSOH leads to the formation of 1-oxatrisulfane (trans-HSSOH and cis-HSSOH) and their branched isomer S(SH)OH. These three secondary product formation rates are competitive, but thermodynamics do not favor the formation of the branched isomer. Among the secondary pyrolysis products, trans-HSSOH is the most stable one, and its formation rate constant at low pressure is calculated to be k(3)(0) = 5.49 × 10(28)T(?10.70) exp(?36.22 kcal mol(?1)/RT) cm3 mol(?1) s(?1) (T = 800?1500 K). Finally, the secondary pyrolysis pathway from less stable product t-BuS(OH)S is also predicted, and that leads to trans-HSSOH and cis-HSSOH products with almost equal rates. A bond-order analysis using Wiberg bond indexes obtained by natural bond orbital (NBO) calculation predicts that the primary and secondary pyrolysis of t-BuS(O)St-Bu occur via E1-like mechanism.     

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).     

Effect of ring substitution on the S-H bond dissociation enthalpies of thiophenols. An experimental and computational study

There are conflicting reports on the origin of the effect of Y substituents on the S-H bond dissociation enthalpies (BDEs) in 4-Y-substituted thiophenols, 4-YC(6)H(4)S-H. The differences in S-H BDEs, [4-YC(6)H(4)S-H] - [C(6)H(5)S-H], are known as the total (de)stabilization enthalpies, TSEs, where TSE = RSE - MSE, i.e., the radical (de)stabilization enthalpy minus the molecule (de)stabilization enthalpy. The effects of 4-Y substituents on the S-H BDEs in thiophenols and on the S-C BDEs in phenyl thioethers are expected to be almost identical. Some S-C TSEs were therefore derived from the rates of homolyses of a few 4-Y-substituted phenyl benzyl sulfides, 4-YC(6)H(4)S-CH(2)C(6)H(5), in the hydrogen donor solvent 9,10-dihydroanthracene. These TSEs were found to be -3.6 +/- 0.5 (Y = NH(2)), -1.8 +/- 0.5 (CH(3)O), 0 (H), and 0.7 +/- 0.5 (CN) kcal mol(-1). The MSEs of 4-YC(6)H(4)SCH(2)C(6)H(5) have also been derived from the results of combustion calorimetry, Calvet-drop calorimetry, and computational chemistry (B3LYP/6-311+G(d,p)). The MSEs of these thioethers were -0.6 +/- 1.1 (NH(2)), -0.4 +/- 1.1 (CH(3)O), 0 (H), -0.3 +/- 1.3 (CN), and -0.8 +/- 1.5 (COCH(3)) kcal mol(-1). Although all the enthalpic data are rather small, it is concluded that the TSEs in 4-YC(6)H(4)SH are largely governed by the RSEs, a somewhat surprising conclusion in view of the experimental fact that the unpaired electron in C(6)H(5)S(*) is mainly localized on the S. The TSEs, RSEs, and MSEs have also been computed for a much larger series of 4-YC(6)H(4)SH and 4-YC(6)H(4)SCH(3) compounds by using a B3P86 methology and have further confirmed that the S-H/S-CH(3) TSEs are dominated by the RSEs. Good linear correlations were obtained for TSE = rho(+)sigma(p)(+)(Y), with rho(+) (kcal mol(-1)) = 3.5 (S-H) and 3.9 (S-CH(3)). It is also concluded that the SH substituent is a rather strong electron donor with a sigma(p)(+)(SH) of -0.60, and that the literature value of -0.03 is in error. In addition, the SH rotational barriers in 4-YC(6)H(4)SH have been computed and it has been found that for strong electron donating (ED) Ys, such as NH(2), the lowest energy conformer has the S-H bond oriented perpendicular to the aromatic ring plane. In this orientation the SH becomes an electron withdrawing (EW) group. Thus, although the OH group in phenols is always in-plane and ED irrespective of the nature of the 4-Y substituent, in thiophenols the SH switches from being an ED group with EW and weak ED 4-Ys, to being an EW group for strong ED 4-Ys.     

Theoretical studies on the intramolecular hydrogen bond and tautomerism of 8-mercaptoquinoline in the gaseous phase and in solution using modern DFT methods

A theoretical quantum chemical study of the intramolecular hydrogen bonding interactions in 8-mercaptoquinoline has been carried out. Special attention has been paid to the rotation of S-H bond and intramolecular proton-transfer reactions. Therewith, the B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p), MPW1K/6-311++G(d,p), MPW1K/6-31+G(2d,2p), BH&HLYP/6-311++G(d,p), and G96LYP/6-311++G(d,p) methods have been used. By means of the Onsager and PCM reaction field methods, the effects of solvent on hydrogen-bond energies, conformational equilibria, rotational barriers, and tautomerism in aqueous solution have been studied. These simulations were done at the MPW1K/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. Natural-bond orbital analysis has been performed to study the intramolecular hydrogen bond (IHB) in the gaseous phase and in aqueous medium. The stability of forms under consideration in solution does not coincide with that in the gaseous phase, underlining a great importance of the electrostatic influence of solvent. Double-proton transfer in the prototropic tautomerization of 8-mercaptoquinoline, one water molecule complex in the gaseous phase and in solution, has been systematically studied. The double-proton transfer occurs concertedly and synchronously. The water-assisted tautomerization is kinetically less, but thermodynamically more favorable, compared to that of the single-proton transfer. As in the case with single-proton transfer, for water-assisted reaction, the tautomerization energies and barrier heights decrease with the increase in dielectric constant, which implies faster and more complete tautomerization of 8-mercaptoquinoline in a polar solvent.     

Analysis of the reactivity on the c7h6 potential energy surface

The reactivity and decomposition kinetics on the C(7)H(6) potential energy surface (PES) were investigated, determining structures of stationary points at the B3LYP/6-31+G(d,p) level and energies at the CCSD(T)/cc-pVTZ level with extension to the complete basis set limit. For the reactions characterized by a significant multireference character, the energies were calculated at the CASPT2/cc-pVTZ level. The portion of the PES investigated consisted of 27 wells connected by 39 saddle points. Of the 27 wells, 16 can be accessed through transition states having activation energies smaller than the dissociation threshold. In agreement with previous theoretical studies, it was found that the main interconversion channel takes place on the singlet PES and connects phenylcarbene, cycloheptatetrane, spiroheptatriene, fulvenallene, and three ethynylcyclopentadiene isomers. Two new mechanisms are proposed for the formation of 5-ethynylcyclopentadiene and for the conversion of spiroheptatriene to fulvenallene. The unimolecular decomposition kinetics was thoroughly investigated. It was found that the fastest high pressure decomposition channel, at the temperatures at which C(7)H(6) undergoes unimolecular decomposition (1500--2000 K), leads to the formation of cyclopentadienylidene and acetylene. The rate of crossing from the singlet to the triplet PES may affect considerably this reaction channel, as it is formally spin forbidden. The alternative pathway, which is the decomposition to fulvenallenyl, is however only a factor of 2--3 slower and significantly less activated (82 vs 96 kcal/mol).     

C2H+H2CO: a new route for formaldehyde removal

The title unknown reaction is theoretically studied at various levels to probe the interaction mechanism between the ethynyl radical (HC triple bond C) and formaldehyde (H(2)C double bond O). The most feasible pathway is a barrier-free direct H-abstraction process leading to acetylene and formyl radical (C(2)H(2)+HCO) via a weakly bound complex, and then the product can take secondary dissociation to the final product C(2)H(2)+CO+H. The C-addition channel leading to propynal plus H-atom (HCCCHO+H) has the barrier of only 3.6, 2.9, and 2.1 kcal/mol at the CCSD(T)/6-311+G(3df,2p)MP2//6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE, and G3//MP2 levels, respectively [CCSD(T)--coupled cluster with single, double, and triple excitations; ZPVE--zero-point vibrational energy; QCISD--quadratic configuration interaction with single and double excitations; G3//MP2-Gaussian-3 based on Moller-Plesset geometry]. The O addition also leading to propynal plus H atom needs to overcome a higher barrier of 5.3, 8.7, and 3.0 kcalmol at the three corresponding levels. The title no-barrier reaction presents a new efficient route to remove the pollutant H(2)CO, and should be included in the combustion models of hydrocarbons. It may also represent the fastest radical-H(2)CO reaction among the available theoretical data. Moreover, it could play an important role in the interstellar chemistry where the zero- or minute-barrier reactions are generally favored. Discussions are also made on the possible formation of the intriguing propynal in space via the title reaction on ice surface.     

氟原子与反式1,3-丁二烯分子反应动力学的理论研究

Theoretical studies of F atom reaction with trans-1,3-butadiene were carried out at the CCSD(T)/6-311G(d,p)/B3LYP/6-311G(d,p) levels. Energies and structures for all reactants, products and transition states were determined. Two reaction pathways involving the formation of the complexes CH2CHCHFCH2 and CH2CHCHCH2F were found in this reaction. Theoretical results suggest that the H atom channel observed in previous crossed beam experiment occurs likely via these two long-lived complex formation pathways. For the complex CH2CHCHFCH2 pathway, another reaction channel (C2H3+C2H3F) is also accessible. Relative importance of the C2H3+C2H3F channel versus the H formation channel via the same reaction pathway has also been estimated, suggesting that it would be difficult to observe the C2H3+C2H3F channel in a crossed molecular beam experiment. Theoretical analysis also shows that the HF formation proceeds via direct abstraction mechanisms, though it is likely a minor process in this reaction.     

An ignored but most favorable channel for NCO+C2H2 reaction

The NCO+C(2)H(2) reaction has been considered as a prototype for understanding the chemical reactivity of the isocyanate radical towards unsaturated hydrocarbons in fuel-rich combustion. It has also been proposed to provide an effective route for formation of oxazole-containing compounds in organic synthesis, and might have potential applications in interstellar processes. Unfortunately, this reaction has met mechanistic controversy both between experiments and between experiments and theoretical calculations. In this paper, detailed theoretical investigations at the Becke's three parameter Lee-Yang-Parr-B3LYP6-31G(d), B3LYP6-311++G(d,p), quadratic configuration interaction with single and double excitations QCISD6-31G(d), and Gaussian-3 levels are performed for the NCO+C(2)H(2) reaction, covering various entrance, isomerization, and decomposition channels. Also, the highly cost-expensive coupled-cluster theory including single and double excitations and perturbative inclusion of triple excitations CCSD(T)/aug-cc-pVTZ single-point energy calculation is performed for the geometries obtained at the Becke's three parameter Lee-Yang-Parr-B3LYP6-311++G(d,p) level. A previously ignored yet most favorable channel via a four-membered ring intermediate with allyl radical character is found. However, formation of P(3) H+HCCNCO and the five-membered ring channel predicted by previous experimental and theoretical studies is kinetically much less competitive. With the new channel, master equation rate constant calculations over a wide range of temperatures (298-1500 K) and pressures (10-560 Torr) show that the predicted total rate constants exhibit a positive-temperature dependence and no distinct pressure dependence effect. This is in qualitative agreement with available experimental results. Under the experimental conditions, the predicted values are about 50% lower than the latest experimental results. Also, the branching ratio variations of the fragments P(2) HCN+HCCO and P(5) OCCHCN+H as well as the intermediates L1 HCHCNCO, r4 cCHCHNC-O, and L5 NCHCHCO are discussed with respect to the temperature and pressure. Future experimental reinvestigations are strongly desired to test the newly predicted channel for the model NCO+C(2)H(2) reaction. Implications of the present results in various fields are discussed.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

Reaction of ketenyl radical with acetylene: a promising route for cyclopropenyl radical

The reaction of the ketenyl radical (HCCO) with acetylene (C(2)H(2)) is very relevant to the oxygen-acetylene flames and fuel-rich combustion process for nitrogen-containing compounds. Unfortunately, except for several rate constant measurements, the mechanism is completely unknown for this reaction. In this paper, detailed theoretical investigations are performed for the HCCO + C(2)H(2) reaction at the G3B3 level using the B3LYP/6-31G(d), B3LYP/6-311++G(d,p), and QCISD/6-31G(d) geometries. The exclusive fragmentation channel is the formation of the cyclopropenyl radical (c-C(3)H(3)) and carbon monoxide (CO) via the chainlike OCCHCHCH and three-membered ring OC-cCHCHCH intermediates. Thus, the mass spectroscopic peak of C(3)H(3)(+) in a previous experiment can be explained. The calculated overall reaction barrier is 4.4, 4.4, and 5.3 kcal/mol at the G3B3//B3LYP/6-31G(d), G3B3//B3LYP/6-311++G(d,p), and G3B3//QCISD/6-31G(d) levels, respectively. The title reaction may provide an effective route for generating the long-sought cyclopropenyl radical in the laboratory, which has been the long-standing subject of numerous theoretical studies as the simplest cyclic conjugate radical, and its bulky derivatives were already known. Future experimental investigations for the HCCO + C(2)H(2) reaction are greatly desired to test the predicted fragmentation channel. The implication of the present study in combustion and interstellar processes is discussed.     

自由基-分子反应:F+Propene(CH2CHCH3)的从头算研究

The radical-molecule reaction F+propene (CH2CHCH3) was studied in detail by using the Becke's three parameter Lee-Yang-Parr-B3LYP/6-311G(d,p) and coupled cluster with single, double, and triple excitationsCCSD(T)/6-311+G(2d,2p). It is shown that F+propene reaction mainly occurs through complex-formation mechanism: F attacks the double bond of propene leading to the formation of complex 1 and complex 2. As the two radical complexes are metastable, they can quickly dissociate to H+C3H5F, CH3+C2H3F and HF+C3H5. Based on the ab initio calculations, the CH3+C2H3F is the main channel, and the H elimination and HF forming channels also provide some contribution to products. The calculated values are in good agreement with the recently reported experimental results.     

Ab initio kinetics for the unimolecular reaction C6H5OH --> CO + C5H6

The unimolecular decomposition of C(6)H(5)OH on its singlet-state potential energy surface has been studied at the G2M//B3LYP/6-311G(d,p) level of theory. The result shows that the most favorable reaction channel involves the isomerization and decomposition of phenol via 2,4-cyclohexadienone and other low-lying isomers prior to the fragmentation process, producing cyclo-C(5)H(6) + CO as major products, supporting the earlier assumption of the important role of the 2,4-cyclohexadienone intermediate. The rate constant predicted by the microcanonical RRKM theory in the temperature range 800-2000 K at 1 Torr--100 atm of Ar pressure for CO production agrees very well with available experimental data in the temperature range studied. The rate constants for the production of CO and the H atom by O-H dissociation at atmospheric Ar pressure can be represented by k(CO) = 8.62 x 10(15) T(-0.61) exp(-37,300/T) s(-1) and k(H) = 1.01 x 10(71) T(-15.92) exp(-62,800/T) s(-1). The latter process is strongly P-dependent above 1000 K; its high- and low-pressure limits are given.     

乙醛二聚体内蓝移型氢键的理论研究

利用理论方法研究了乙醛二聚体内的氢键. 在MP2/6-31＋G(d), B3LYP/6-31＋G(d), B3LYP/6-311＋＋G(d,p)和B3LYP/6-311＋＋G(3df,2p)水平上, 利用常规方法和均衡校正方法对3种稳定的乙醛二聚体进行了几何优化和振动频率计算. 计算结果表明: 在二聚体A和C中乙醛中C—H键强烈收缩, 存在显著的C—H…O蓝移型氢键. 自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、分子内电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 分子内电子密度重排、轨道再杂化和电子受体内部结构重组属于键收缩效应. 在二聚体A和C中, 由于键收缩效应处于优势地位导致C—H…O蓝移氢键存在.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.