The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.     

A coupled reference interaction site model/molecular dynamics study of the potential of mean force curve of the SN2 Cl- + CH3Cl reaction in water

An application of the coupled reference interaction site model (RISM)/simulation methodology to the calculation of the potential of mean force (PMF) curve in aqueous solution for the identity nucleophilic substitution reaction Cl(-) + CH(3)Cl is performed. The free energy of activation is calculated to be 27.1 kcal/mol which compares very well with the experimentally determined barrier height of 26.6 kcal/mol. Furthermore, the calculated PMF is almost superimposed with that previously calculated using the computationally rigorous Monte Carlo with importance sampling method (Chandrasekhar, J.; Smith, S. F.; Jorgensen, W. L. J. Am. Chem. Soc. 1985, 107, 154). Using the calculated PMF, a crude estimate of the solvated kinetic transmission coefficient also compares well with that of previous more accurate simulations. These results indicate that the coupled RISM/simulation method provides a cost-effective methodology for studying reactions in solution.     

The effect of the torsional and stretching vibrations of C2H6 on the H + C2H6 --> H2 + C2H5 reaction

We present a three-dimensional quantum scattering model to treat reactions of the type H + C2H6 --> H2 + C2H5. The model allows the torsional and the stretching degrees of freedom to be treated explicitly. Zero-point energies of the remaining modes are taken into account in electronic structure calculations. An analytical potential-energy surface was developed from a minimal number of ab initio geometry evaluations using the CCSD(T,full)/cc-pVTZ//MP2(full)/cc-pVTZ level of theory. The reaction is endothermic by 1.5 kcal mol(-1) and exhibits a vibrationally adiabatic barrier of 12.0 kcal mol(-1). The results show that the torsional mode influences reactivity when coupled with the vibrational C-H stretching mode. We also found that ethyl radical products are formed internally excited in the torsional mode.     

Arene-cation interactions of positive quadrupole moment aromatics and arene-anion interactions of negative quadrupole moment aromatics

Intermolecular interactions involving aromatic pi-electron density are widely believed to be governed by the aromatic molecular quadrupole moment, Theta(zz). Arene-cation binding is believed to occur primarily with negative Theta(zz) aromatics, and arene-anion binding is believed to occur largely with positive Theta(zz) aromatics. We have performed quantum mechanical computations that show the cation binding of positive Theta(zz) aromatics and the anion binding of negative Theta(zz) aromatics is quite common in the gas phase. The pi-electron density of hexafluorobenzene, the prototypical positive Theta(zz ) aromatic (experimental Theta(zz) = 9.5 +/- 0.5 DA), has a Li+ binding enthalpy of -4.37 kcal/mol at the MP2(full)/6-311G**level of theory. The RHF/6-311G** calculated Theta(zz) value of 1,4-dicyanobenzene is +11.81 DA, yet it has an MP2(full)/6-311G** Li+ binding enthalpy of -12.65 kcal/mol and a Na+ binding enthalpy of -3.72 kcal/mol. The pi-electron density of benzene, the prototypical negative Theta(zz) aromatic (experimental Theta(zz) = -8.7 +/- 0.5 DA), has a F- binding enthalpy of -5.51 kcal/mol. The RHF/6-311G** calculated Theta(zz) of C6H2I4 is -10.45 DA, yet it has an MP2(full)/6-311++G** calculated F- binding enthalpy of -20.13 kcal/mol. Our results show that as the aromatic Theta(zz) value increases the cation binding enthalpy decreases; a plot of cation binding enthalpies versus aromatic Theta(zz) gives a line of best of fit with R2 = 0.778. No such correlation exists between the aromatic Theta(zz) value and the anion binding enthalpy; the line of best fit has R2 = 0.297. These results are discussed in terms of electrostatic and polarizability contributions to the overall binding enthalpies.     

Ab initio and direct dynamics studies of the reaction of singlet methylene with acetylene and the lifetime of the cyclopropene complex

The energetics of the (1)CH(2) + C(2)H(2) --> H + C(3)H(3) reaction are accurately calculated using an extrapolated coupled-cluster/complete basis set (CBS) method based on the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. The reaction enthalpy (0 K) is predicted to be -20.33 kcal/mol. This reaction has no classical barrier in either the entrance or exit channel. However, there are several stable intermediates-cyclopropene (c-C(3)H(4)), allene (CH(2)CCH(2)), and propyne (CH(3)CCH)-along the minimum energy path. These intermediates with zero-point energy corrections lie below the reactants by 87.11 (c-C(3)H(4)), 109.69 (CH(2)CCH(2)), and 110.78 kcal/mol (CH(3)CCH). The vibrationally adiabatic ground-state (VAG) barrier height for c-C(3)H(4) isomerization to allene is obtained as 45.2 kcal/mol, and to propyne as 37.2 kcal/mol. In addition, the (1)CH(2) + C(2)H(2) reaction is investigated utilizing the dual-level "scaling all correlation" (SAC) ab initio method of Truhlar et al., i.e., the UCCSD(SAC)/cc-pVDZ theory. Results show that the reaction occurs via long-lived complexes. The lifetime of the cyclopropene intermediate is obtained as 3.2 +/- 0.4 ps. It is found that the intermediate propyne can be formed directly from reactants through the insertion of (1)CH(2) into a C-H bond of C(2)H(2). However, compared to the major mechanism in which the propyne is produced through a ring-opening of the cyclopropene complex, this reaction pathway is much less favorable. Finally, the theoretical thermal rate constant exhibits a negative temperature dependence, which is in excellent agreement with the previous results. The temperature dependence is consistent with the earlier RRKM results but weaker than the experimental observations at high temperatures.     

The reaction of tricarbon with acetylene: an ab initio/RRKM study of the potential energy surface and product branching ratios

Ab initio calculations of the potential energy surface for the C3(1Sigmag+)+C2H2(1Sigmag+) reaction have been performed at the RCCSD(T)/cc-pVQZ//B3LYP/6-311G(d,p) + ZPE[B3LYP/6-311G(d,p)] level with extrapolation to the complete basis set limit for key intermediates and products. These calculations have been followed by statistical calculations of reaction rate constants and product branching ratios. The results show the reaction to begin with the formation of the 3-(didehydrovinylidene)cyclopropene intermediate i1 or five-member ring isomer i7 with the entrance barriers of 7.6 and 13.8 kcal/mol, respectively. i1 rearranges to the other C5H2 isomers, including ethynylpropadienylidene i2, singlet pentadiynylidene i3, pentatetraenylidene i4, ethynylcyclopropenylidene i5, and four- and five-member ring structures i6, i7, and i8 by ring-closure and ring-opening processes and hydrogen migrations. i2, i3, and i4 lose a hydrogen atom to produce the most stable linear isomer of C5H with the overall reaction endothermicity of approximately 24 kcal/mol. H elimination from i5 leads to the formation of the cyclic C5H isomer, HC2C3, +H, 27 kcal/ mol above C3+C2H2. 1,1-H2 loss from i4 results in the linear pentacarbon C5+H2 products endothermic by 4 kcal/mol. The H elimination pathways occur without exit barriers, whereas the H2 loss from i4 proceeds via a tight transition state 26.4 kcal/mol above the reactants. The characteristic energy threshold for the reaction under single collision conditions is predicted be in the range of approximately 24 kcal/mol. Product branching ratios obtained by solving kinetic equations with individual rate constants calculated using RRKM and VTST theories for collision energies between 25 and 35 kcal/mol show that l-C5H+H are the dominant reaction products, whereas HC2C3+H and l-C5+H2 are minor products with branching ratios not exceeding 2.5% and 0.7%, respectively. The ethynylcyclopropenylidene isomer i5 is calculated to be the most stable C5H2 species, more favorable than triplet pentadiynylidene i3t by approximately 2 kcal/mol.     

Gaseous reaction mechanism of C2F radical with water

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.     

Quantum Monte Carlo study of the reaction: Cl+CH3OH-->CH2OH+HCl

A theoretical study is reported of the Cl+CH3OH-->CH2OH+HCl reaction based on the diffusion Monte Carlo (DMC) variant of the quantum Monte Carlo method. Using a DMC trial function constructed as a product of Hartree-Fock and correlation functions, we have computed the barrier height, heat of reaction, atomization energies, and heats of formation of reagents and products. The DMC heat of reaction, atomization energies, and heats of formation are found to agree with experiment to within the error bounds of computation and experiment. M?ller-Plesset second order perturbation theory (MP2) and density functional theory, the latter in the B3LYP generalized gradient approximation, are found to overestimate the experimental heat of reaction. Intrinsic reaction coordinate calculations at the MP2 level of theory demonstrate that the reaction is predominantly direct, i.e., proceeds without formation of intermediates, which is consistent with a recent molecular beam experiment. The reaction barrier as determined from MP2 calculations is found to be 2.24 kcal/mol and by DMC it is computed to be 2.39(49) kcal/mol.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Quantum chemical study of the mechanism of reaction between NH (X 3sigma-) and H2, H2O, and CO2 under combustion conditions

Reactions of ground-state NH (3sigma-) radicals with H2, H2O, and CO2 have been investigated quantum chemically, whereby the stationary points of the appropriate reaction potential energy surfaces, that is, reactants, products, intermediates, and transition states, have been identified at the G3//B3LYP level of theory. Reaction between NH and H2 takes place via a simple abstraction transition state, and the rate coefficient for this reaction as derived from the quantum chemical calculations, k(NH + H2) = (1.1 x 10(14)) exp(-20.9 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K, is found to be in good agreement with experiment. For reaction between triplet NH and H2O, no stable intermediates were located on the triplet reaction surface although several stable species were found on the singlet surface. No intersystem crossing seam between triplet NH + H2O and singlet HNO + H2 (the products of lowest energy) was found; hence there is no evidence to support the existence of a low-energy pathway to these products. A rate coefficient of k(NH + H2O) = (6.1 x 10(13)) exp(-32.8 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K for the reaction NH (3sigma-) + H2O --> NH2 (2B) + OH (2pi) was derived from the quantum chemical results. The reverse rate coefficient, calculated via the equilibrium constant, is in agreement with values used in modeling the thermal de-NO(x) process. For the reaction between triplet NH and CO2, several stable intermediates on both triplet and singlet reaction surfaces were located. Although a pathway from triplet NH + CO2 to singlet HNO + CO involving intersystem crossing in an HN-CO2 adduct was discovered, no pathway of sufficiently low activation energy was discovered to compare with that found in an earlier experiment [Rohrig, M.; Wagner, H. G. Proc. Combust. Inst. 1994, 25, 993.].     

HNNC radical and its role in the CH+N2 reaction

A previously unreported channel in the spin-allowed reaction path for the CH+N2 reaction that involves the HNNC radical is presented. The structures and energetics of the HNNC radical and its isomers HCNN and HNCN and the relevant intermediates and transition states that are involved in the proposed mechanism are obtained at the coupled cluster singles and doubles level of theory with noniterative triples correction (CCSD(T)) using a converging series of basis sets aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ. The aug-cc-pVQZ basis is used for all the final single point energy calculations using the CCSD(T)/aug-cc-pVTZ optimized geometries. We find the HNNC radical to have a heat of formation of DeltafH0 (HNNC)=116.5 kcal mol(-1). An assessment of the quality of computed data of the radical species HNCN and HCNN is presented by comparison with the available experimental data. We find that HNNC can convert to HNCN, the highest barrier in this path being 14.5 kcal mol(-1) above the energy of the CH+N2 reactants. Thus, HNNC can play a role in the high-temperature spin-allowed mechanism for the reaction of CH+N2 proposed by Moskaleva, Xia, and Lin (Chem. Phys. Lett. 2000, 331, 269).     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

Carboranes and baskets from reaction of B4H10 with allene

The reaction of the boron hydride B4H10 with allene was studied at the CCSD(T)/6-311+G(d)//MP2/6-31G(d) level. The mechanism is surprisingly complex with 44 transition states and several branching points located. The four carboranes and one basket that have been observed experimentally are all connected by pathways that have very similar free energies of activation. In addition, two new structures, a basket (2,4-(CH2CH2CH2)B4H8, 5a) and a "classical" structure (1,4-(Me2C)bisdiborane, 7), which might be obtained from the B4H10 + C3H4 reaction under the right conditions (hot/cold, quenched, etc.) have been identified. The first branch point in the reaction is the competition between H2 elimination from B4H10 (DeltaG(298 K) = 32.2 kcal/mol) and the hydroboration of allene by B4H10 (DeltaG(298 K) = 31.3 kcal/mol). The next branch point in the hydroboration mechanism controls the formation of 2,4-(MeCHCH2)B4H8 (1) (DeltaG(298 K) = 31.5 kcal/mol) and arachno-1,2/arachno-1,3-Me2-1-CB4H7 (8 and 8a) (DeltaG(298 K) = 34.3 kcal/mol). Another branch point in the H2-elimination mechanism controls the formation of 1-Me-2,5-micro-CH2-1-CB4H7 (29) (DeltaG(298 K) = 0.1 kcal/mol) and 2,5-micro-CHMe-1-CB4H7 (25/26) (DeltaG(298 K) = 7.3 kcal/mol). Formation of 2-Me-2,3-C2B4H7, a carborane observed in the reaction of methylacetylene with B4H10, is calculated to be blocked by a high barrier for H2 elimination. All free energies are relative to B4H10 + allene. An interesting reaction step discovered is the "reverse hydroboration step" in which a hydrogen atom is transferred from carbon back to boron, which allows a CH hydrogen to shuttle between the terminal and central carbon of allene.     

Reaction mechanism and tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution: a combined Monte Carlo and quantum mechanics study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order M?ller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.     

Studies of the kinetics and thermochemistry of the forward and reverse reaction Cl + C6H6 = HCl + C6H5

The laser flash photolysis resonance fluorescence technique was used to monitor atomic Cl kinetics. Loss of Cl following photolysis of CCl4 and NaCl was used to determine k(Cl + C6H6) = 6.4 x 10(-12) exp(-18.1 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) over 578-922 K and k(Cl + C6D6) = 6.2 x 10(-12) exp(-22.8 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) over 635-922 K. Inclusion of literature data at room temperature leads to a recommendation of k(Cl + C6H6) = 6.1 x 10(-11) exp(-31.6 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) for 296-922 K. Monitoring growth of Cl during the reaction of phenyl with HCl led to k(C6H5 + HCl) = 1.14 x 10(-12) exp(+5.2 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) over 294-748 K, k(C6H5 + DCl) = 7.7 x 10(-13) exp(+4.9 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) over 292-546 K, an approximate k(C6H5 + C6H5I) = 2 x 10(-11) cm(3) molecule(-1) s(-1) over 300-750 K, and an upper limit k(Cl + C6H5I) < or = 5.3 x 10(-12) exp(+2.8 kJ mol(-1)/RT) cm(3) molecule(-1) s(-1) over 300-750 K. Confidence limits are discussed in the text. Third-law analysis of the equilibrium constant yields the bond dissociation enthalpy D(298)(C6H5-H) = 472.1 +/- 2.5 kJ mol(-1) and thus the enthalpy of formation Delta(f)H(298)(C6H5) = 337.0 +/- 2.5 kJ mol(-1).     

A multilayered representation,quantum mechanical and molecular mechanics study of the CH3F + OH− reaction in water

The bimolecular nucleophilic substitution (S_N2) reaction of CH₃F + OH^? in aqueous solution was investigated using a combined quantum mechanical and molecular mechanics approach. Reactant complex, transition state, and product complex along the reaction pathway were analyzed in water. The potentials of mean force were calculated using a multilayered representation with the DFT and CCSD(T) level of theory for the reactive region. The obtained free energy activation barrier for this reaction at the CCSD(T)/MM representation is 18.3 kcal/mol which agrees well with the experimental value at ～21.6 kcal/mol. Both the solvation effect and solute polarization effect play key roles on raising the activation barrier height in aqueous solution, with the former raising the barrier height by 3.1 kcal/mol, the latter 1.5 kcal/mol. © 2013 Wiley Periodicals, Inc.     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Theoretical study of the C6H3 potential energy surface and rate constants and product branching ratios of the C2H(2Sigma+) + C4H2(1Sigma(g)+) and C4H(2Sigma+) + C2H2(1Sigma(g)+) reactions

Ab initio CCSD(T)cc-pVTZ//B3LYP6-311G(**) and CCSD(T)/complete basis set (CBS) calculations of stationary points on the C(6)H(3) potential energy surface have been performed to investigate the reaction mechanism of C(2)H with diacetylene and C(4)H with acetylene. Totally, 25 different C(6)H(3) isomers and 40 transition states are located and all possible bimolecular decomposition products are also characterized. 1,2,3- and 1,2,4-tridehydrobenzene and H(2)CCCCCCH isomers are found to be the most stable thermodynamically residing 77.2, 75.1, and 75.7 kcal/mol lower in energy than C(2)H + C(4)H(2), respectively, at the CCSD(T)/CBS level of theory. The results show that the most favorable C(2)H + C(4)H(2) entrance channel is C(2)H addition to a terminal carbon of C(4)H(2) producing HCCCHCCCH, 70.2 kcal/mol below the reactants. This adduct loses a hydrogen atom from the nonterminal position to give the HCCCCCCH (triacetylene) product exothermic by 29.7 kcal/mol via an exit barrier of 5.3 kcal/mol. Based on Rice-Ramsperger-Kassel-Marcus calculations under single-collision conditions, triacetylene+H are concluded to be the only reaction products, with more than 98% of them formed directly from HCCCHCCCH. The C(2)H + C(4)H(2) reaction rate constants calculated by employing canonical variational transition state theory are found to be similar to those for the related C(2)H + C(2)H(2) reaction in the order of magnitude of 10(-10) cm(3) molecule(-1) s(-1) for T = 298-63 K, and to show a negative temperature dependence at low T. A general mechanism for the growth of polyyne chains involving C(2)H + H(C[triple bond]C)(n)H --> H(C[triple bond]C)(n+1)H + H reactions has been suggested based on a comparison of the reactions of ethynyl radical with acetylene and diacetylene. The C(4)H + C(2)H(2) reaction is also predicted to readily produce triacetylene + H via barrierless C(4)H addition to acetylene, followed by H elimination.