Accurate ab initio predictions of ionization energies and heats of formation for the 2-propyl, phenyl, and benzyl radicals

The ionization energies (IEs) for the 2-propyl (2-C(3)H(7)), phenyl (C(6)H(5)), and benzyl (C(6)H(5)CH(2)) radicals have been calculated by the wave-function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled cluster level with single and double excitations plus quasiperturbative triple excitation [CCSD(T)]. The zero-point vibrational energy correction, the core-valence electronic correction, and the scalar relativistic effect correction have been also made in these calculations. Although a precise IE value for the 2-C(3)H(7) radical has not been directly determined before due to the poor Franck-Condon factor for the photoionization transition at the ionization threshold, the experimental value deduced indirectly using other known energetic data is found to be in good accord with the present CCSD(T)/CBS prediction. The comparison between the predicted value through the focal-point analysis and the highly precise experimental value for the IE(C(6)H(5)CH(2)) determined in the previous pulsed field ionization photoelectron (PFI-PE) study shows that the CCSD(T)/CBS method is capable of providing an accurate IE prediction for C(6)H(5)CH(2), achieving an error limit of 35 meV. The benchmarking of the CCSD(T)/CBS IE(C(6)H(5)CH(2)) prediction suggests that the CCSD(T)/CBS IE(C(6)H(5)) prediction obtained here has a similar accuracy of 35 meV. Taking into account this error limit for the CCSD(T)/CBS prediction and the experimental uncertainty, the CCSD(T)/CBS IE(C(6)H(5)) value is also consistent with the IE(C(6)H(5)) reported in the previous HeI photoelectron measurement. Furthermore, the present study provides support for the conclusion that the CCSD(T)/CBS approach with high-level energy corrections can be used to provide reliable IE predictions for C(3)-C(7) hydrocarbon radicals with an uncertainty of +/-35 meV. Employing the atomization scheme, we have also computed the 0 K (298 K) heats of formation in kJ/mol at the CCSD(T)/CBS level for 2-C(3)H(7)/2-C(3)H(7) (+) ,C(6)H(5)/C(6)H(5) (+), and C(6)H(5)CH(2)/C(6)H(5)CH(2) (+) to be 105.2/822.7 (90.0/806.4), 351.4/1148.5 (340.4/1138.8), and 226.2/929.0 (210.3/912.7), respectively. Comparing these values with the available experimental values, we find that the discrepancies for the 0 and 298 K heats of formation values are < or =2.6 kJ/mol for 2-C(3)H(7)/2-C(3)H(7) (+),< or =4.1 kJ/mol for C(6)H(5)/C(6)H(5) (+), and < or =3.2 kJ//mol for C(6)H(5)CH(2)C(6)H(5)CH(2) (+).     

Accurate ab initio predictions of ionization energies of hydrocarbon radicals: CH2, CH3, C2H, C2H3, C2H5, C3H3, and C3H5

The ionization energies for methylene (CH2), methyl (CH3), ethynyl (C2H), vinyl (C2H3), ethyl (C2H5), propargyl (C3H3), and allyl (C3H5) radicals have been calculated by the wave-function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled-cluster level with single and double excitations plus a quasiperturbative triple excitation [CCSD(T)]. When it is appropriate, the zero-point vibrational energy correction, the core-valence electronic correction, the scalar relativistic effect correction, the diagonal Born-Oppenheimer correction, and the high-order correlation correction have also been made in these calculations. The comparison between the computed ionization energy (IE) values and the highly precise experimental IE values determined in previous pulsed field ionization-photoelectron (PFI-PE) studies indicates that the CCSD(T)/CBS method is capable of providing accurate IE predictions for these hydrocarbon radicals achieving error limits well within +/-10 meV. The benchmarking of the CCSD(T)/CBS IE predictions by the PFI-PE experimental results also lends strong support for the conclusion that the CCSD(T)/CBS approach with high-level energy corrections can serve as a valuable alternative for reliable IE determination of radicals, particularly for those radicals with very unfavorable Franck-Condon factors for photoionization transitions near their ionization thresholds.     

High-level ab initio predictions for the ionization energies and heats of formation of five-membered-ring molecules: thiophene, furan, pyrrole, 1,3-cyclopentadiene, and borole, C4H4X/C4H4X+ (X = S, O, NH, CH2, and BH)

The ionization energies (IEs) and heats of formation (ΔH°(f0)/ΔH°(f298)) for thiophene (C(4)H(4)S), furan (C(4)H(4)O), pyrrole (C(4)H(4)NH), 1,3-cyclopentadiene (C(4)H(4)CH(2)), and borole (C(4)H(4)BH) have been calculated by the wave function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled-cluster level with single and double excitations plus a quasi-perturbative triple excitation [CCSD(T)]. Where appropriate, the zero-point vibrational energy correction (ZPVE), the core-valence electronic correction (CV), and the scalar relativistic effect (SR) are included in these calculations. The respective CCSD(T)/CBS predictions for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2), being 8.888, 8.897, 8.222, and 8.582 eV, are in excellent agreement with the experimental values obtained from previous photoelectron and photoion measurements. The ΔH°(f0)/ΔH°(f298) values for the aforementioned molecules and their corresponding cations have also been predicted by the CCSD(T)/CBS method, and the results are compared with the available experimental data. The comparisons between the CCSD(T)/CBS predictions and the experimental values for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2) suggest that the CCSD(T)/CBS procedure is capable of predicting reliable IE values for five-membered-ring molecules with an uncertainty of ±13 meV. In view of the excellent agreements between the CCSD(T)/CBS predictions and the experimental values for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2), the similar CCSD(T)/CBS IE and ΔH°(f0)/ΔH°(f298) predictions for C(4)H(4)BH, whose thermochemical data are not readily available due to its reactive nature, should constitute a reliable data set. The CCSD(T)/CBS IE(C(4)H(4)BH) value is 8.868 eV, and ΔH°(f0)/ΔH°(f298) values for C(4)H(4)BH and C(4)H(4)BH(+) are 269.5/258.6 and 1125.1/1114.6 kJ/mol, respectively. The highest occupied molecular orbitals (HOMO) of C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, C(4)H(4)CH(2), and C(4)H(4)BH have also been studied by the natural bond orbital (NBO) method, and the extent of π-electron delocalization in these five-membered rings are discussed in correlation with their molecular structures and orbitals.     

Photodissociation of 1-bromo-2-butene, 4-bromo-1-butene, and cyclopropylmethyl bromide at 234 nm studied using velocity map imaging

We present photofragment imaging experiments to characterize potential photolytic precursors of three C4H7 radical isomers: 1-methylallyl, cyclopropylmethyl, and 3-buten-1-yl radicals. The experiments use 2+1 resonance enhanced multiphoton ionization (REMPI) with velocity map imaging to state-selectively detect the Br(2P(3/2)) and Br(2P(1/2)) atoms as a function of their recoil velocity imparted upon photodissociation of 1-bromo-2-butene, cyclopropylmethyl bromide, and 4-bromo-1-butene at 234 nm as well as the angular distributions of the photofragments. Energy and momentum conservation allows the internal energy distribution of the nascent momentum-matched radicals to be derived. The radicals are detected with single photon photoionization at 157 nm. In the case of the 1-methylallyl radical the photoionization cross section is expected to be independent of internal energy in the range of 7-30 kcal/mol. Thus, comparison of the product recoil kinetic energy distribution derived from the measurement of the 1-methylallyl velocity distribution, detecting the radicals with 157 nm photoionization, with a linear combination of the Br atom recoil kinetic energy distributions allows us to derive reliable REMPI line strength ratios for the detection of Br atoms and to test the assumption that the photoionization cross section does not strongly depend on the internal energy of the radical. This line strength ratio is then used to determine the branching to the Br(2P(3/2)) and Br(2P(1/2)) product channels for the other two photolytic systems and to determine the internal energy distribution of their momentum-matched radicals. (We also revisit earlier work on the photodissociation of cyclobutyl bromide which detected the Br atoms and momentum-matched cyclobutyl radicals.) This allows us to test whether the 157 nm photoionization of these radicals is insensitive to internal energy for the distribution of total internal (vibrational+rotational) energy produced. We find that 157 nm photoionization of cyclopropylmethyl radicals is relatively insensitive to internal energy, while 3-buten-1-yl radicals show a photoionization cross section that is markedly dependent on internal energy with the lowest internal energy radicals not efficiently detected by photoionization at 157 nm. We present electronic structure calculations of the radicals and their cations to understand the experimental results.     

Photodissociation of cyclobutyl bromide at 234 nm studied using velocity map imaging

This study investigates the 234 nm photodissociation dynamics of cyclobutyl bromide using a two-dimensional photofragment velocity imaging technique. The spin-orbit ground- and excited-state Br(2P) atoms are state-selectively detected via [2+1] resonance enhanced multiphoton ionization (REMPI), whereas the cyclobutyl radicals are ionized using 157 nm laser light. The Br(2P(3/2)) and the Br(2P(1/2)) atoms and their c-C4H7 radical cofragments evidence a single-peaked, Gaussian-shaped translational energy distribution ranging from approximately 14 to approximately 39 kcal/mol and angular distributions with significant parallel character. The Br(2P(1/2))/ Br(2P(3/2)) spin-orbit branching ratio is determined to be 0.11 +/- 0.07 by momentum match between the Br(2P) photofragments and the recoiling c-C4H7 fragments, assuming a uniform photoionization probability of the c-C4H7 radicals with an internal energy range of 10-35 kcal/mol. The REMPI line strength ratio for the detection of Br(2P(3/2)) and Br(2P(1/2)) atoms at 233.681 and 234.021 nm, respectively, is therefore derived to be 0.10 +/- 0.07. The measured recoil kinetic energies of the c-C4H7 radicals, and the resulting distribution of internal energies, indicates some of the radicals are formed with total internal energies above the barrier to isomerization and subsequent dissociation, but our analysis indicates they may be stable due to the substantial fraction of the internal energy which is partitioned to rotational energy of the radicals.     

Photoionization of propargyl and bromopropargyl radicals: a threshold photoelectron spectroscopic study

In this Article, we present mass-selected threshold photoelectron spectra of propargyl as well as the 1- and 3-bromopropargyl radicals. The reactive intermediates were produced by flash pyrolysis of suitable precursors and ionized by VUV synchrotron radiation. The TPES of the propargyl radical was simulated using data from a recent high-level computational study. An ionization energy (IE) of 8.71 ± 0.02 eV was obtained, in excellent agreement with computations, but slightly above previous experimental IEs. The pyrolysis of 1,3-dibromopropyne delivers both 1- and 3-bromopropargyl radicals that can be distinguished by their different ionization energies (8.34 and 8.16 eV). To explain the vibrational structure, a Franck-Condon simulation was performed, based on DFT calculations, which can account for all major spectral features. Bromopropargyl photoionizes dissociatively beginning at around 10.1 eV. Cationic excited states of 1- and 3-bromopropargyl were tentatively identified. The dissociative photoionization of the precursor (1,3-dibromopropyne) was also examined, delivering an AE(0K) (C(3)H(2)Br(+)/C(3)H(2)Br(2)) of 10.6 eV.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Theoretical study on the isomerization behavior between alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals

[reaction: see text] Ab initio calculations using 6-311G**, cc-pVDZ, aug-cc-pVDZ, and a (valence) double-zeta pseudopotential (DZP) basis set, with (QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional methods (BHandHLYP, B3LYP) predict that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals exist as isomers. At the CCSD(T)/cc-pVDZ//BHandHLY/cc-pVDZ level of theory, energy barriers of 15.1 and 17.7-21.7 kJ mol(-)(1) are calculated for the isomerization of s-trans-propenoyl and s-trans-crotonoyl radical to ketenylmethyl and 1-ketenylethyl radical, respectively. Similar results are obtained for the reactions of s-trans isomers involving silyl, germyl, and stannyl groups with energy barriers (DeltaE++) of 12.2-12.4, 13.1-13.9, and 12.9-18.2 kJ mol(-)(1) at the CCSD(T)/DZP//BHandHLYP/DZP calculation, respectively. These results suggest that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals are not canonical forms but are isomeric species that can rapidly interconvert.     

Photodissociation dynamics of the 2-methylallyl radical

Time- and frequency-resolved photoionization of the hydrogen atom product from a jet-cooled electronically excited 2-methylallyl radical, C4H7, provides information on the dissociation dynamics. The measured dissociation rates and kinetic energy release of 2-methylallyl and its isotopologue CD3C3H4 combined with high level ab initio calculations suggests unimolecular dissociation with methylenecyclopropane and hydrogen as the major C-H bond fission channel with no evidence for nonstatistical behavior in dissociation. Other possible dissociation and isomerization pathways are discussed based on the results of the coupled-cluster ab initio calculations.     

Electron super-rich radicals in the gas phase. A neutralization-reionization mass spectrometric and ab initio/RRKM study of diaminohydroxymethyl and triaminomethyl radicals

Diaminohydroxymethyl (1) and triaminomethyl (2) radicals were generated by femtosecond collisional electron transfer to their corresponding cations (1+ and 2+, respectively) and characterized by neutralization-reionization mass spectrometry and ab initio/RRKM calculations at correlated levels of theory up to CCSD(T)/aug-cc-pVTZ. Ion 1+ was generated by gas-phase protonation of urea which was predicted to occur preferentially at the carbonyl oxygen with the 298 K proton affinity that was calculated as PA = 875 kJ mol-1. Upon formation, radical 1 gains vibrational excitation through Franck-Condon effects and rapidly dissociates by loss of a hydrogen atom, so that no survivor ions are observed after reionization. Two conformers of 1, syn-1 and anti-1, were found computationally as local energy minima that interconverted rapidly by inversion at one of the amine groups with a <7 kJ mol-1 barrier. The lowest energy dissociation of radical 1 was loss of the hydroxyl hydrogen atom from anti-1 with ETS = 65 kJ mol-1. The other dissociation pathways of 1 were a hydroxyl hydrogen migration to an amine group followed by dissociation to H2N-C=O* and NH3. Ion 2+ was generated by protonation of gas-phase guanidine with a PA = 985 kJ mol-1. Electron transfer to 2+ was accompanied by large Franck-Condon effects that caused complete dissociation of radical 2 by loss of an H atom on the experimental time scale of 4 mus. Radicals 1 and 2 were calculated to have extremely low ionization energies, 4.75 and 4.29 eV, respectively, which belong to the lowest among organic molecules and bracket the ionization energy of atomic potassium (4.34 eV). The stabilities of amino group containing methyl radicals, *CH2NH2, *CH(NH2)2, and 2, were calculated from isodesmic hydrogen atom exchange with methane. The pi-donating NH2 groups were found to increase the stability of the substituted methyl radicals, but the stabilities did not correlate with the radical ionization energies.     

High-resolution threshold photoelectron study of the propargyl radical by the vacuum ultraviolet laser velocity-map imaging method

By employing the vacuum ultraviolet (VUV) laser velocity-map imaging (VMI) photoelectron scheme to discriminate energetic photoelectrons, we have measured the VUV-VMI-threshold photoelectrons (VUV-VMI-TPE) spectra of propargyl radical [C(3)H(3)(X?(2)B(1))] near its ionization threshold at photoelectron energy bandwidths of 3 and 7 cm(-1) (full-width at half-maximum, FWHM). The simulation of the VUV-VMI-TPE spectra thus obtained, along with the Stark shift correction, has allowed the determination of a precise value 70 156 ± 4 cm(-1) (8.6982 ± 0.0005 eV) for the ionization energy (IE) of C(3)H(3). In the present VMI-TPE experiment, the Stark shift correction is determined by comparing the VUV-VMI-TPE and VUV laser pulsed field ionization-photoelectron (VUV-PFI-PE) spectra for the origin band of the photoelectron spectrum of the X?(+)-X? transition of chlorobenzene. The fact that the FWHMs for this origin band observed using the VUV-VMI-TPE and VUV-PFI-PE methods are nearly the same indicates that the energy resolutions achieved in the VUV-VMI-TPE and VUV-PFI-PE measurements are comparable. The IE(C(3)H(3)) value obtained based on the VUV-VMI-TPE measurement is consistent with the value determined by the VUV laser PIE spectrum of supersonically cooled C(3)H(3)(X?(2)B(1)) radicals, which is also reported in this article.     

On the thermal unimolecular decomposition of the cyclohexoxy radical--an experimental and theoretical study

The kinetics of the thermal unimolecular decomposition of the cyclohexoxy radical (c-C(6)H(11)O) was experimentally studied, and the results were analyzed in terms of statistical rate theory with molecular and transition state data from quantum chemical calculations. Laser flash photolysis of cyclohexylnitrite at 351 nm was used to produce c-C(6)H(11)O radicals, and their concentration was monitored by laser-induced fluorescence after excitation at 356.2 or 365.2 nm. The experiments were performed at temperatures ranging from 293 to 341 K and pressures between 5 and 55 bar with helium as the bath gas. Over the whole temperature range, biexponential profiles were observed, which is an indication of a consecutive reaction with a pre-equilibrium. From our quantum chemical calculations, it follows that this pre-equilibrium corresponds to the reversible ring-opening via beta-C-C bond fission to form the 6-oxo-1-hexyl radical (l-C(6)H(11)O), c-C(6)H(11)O <--> l-C(6)H(11)O (1,-1). The following temperature-dependent rate coefficients were deduced with an estimated uncertainty of +/-30%: k(1)(T) = 3.80 x 10(13) exp(-50.1 kJ mol(-1)/RT) s(-1) and k(-1)(T) = 3.02 x 10(8) exp(-23.8 kJ mol(-1)/RT) s(-1); a pressure dependence was not observed. In our theoretical analysis, the different conformers of c-C(6)H(11)O were explicitly taken into account, and the C-C torsional motions in l-C(6)H(11)O were treated as hindered internal rotators using a recently suggested approach. This explicit consideration of the hindered internal rotators significantly improved the agreement between the experimentally determined rate coefficients and the results from the quantum chemical computations.     

Cis-trans isomerization of chemically activated 1-methylallyl radical and fate of the resulting 2-buten-1-peroxy radical

The cis-trans isomerization of chemically activated 1-methylallyl is investigated using RRKM/Master Equation methods for a range of pressures and temperatures. This system is a prototype for a large range of allylic radicals formed from highly exothermic (～35 kcal/mol) OH + alkene reactions. Energies, vibrational frequencies, anharmonic constants, and the torsional potential of the methyl group are computed with density functional theory for both isomers and the transition state connecting them. Chemically activated radicals are found to undergo rapid cis-trans isomerization leading to stabilization of significant amounts of both isomers. In addition, the thermal rate constant for trans → cis isomerization of 1-methylallyl is computed to be high enough to dominate reaction with O(2) in 10 atm of air at 700 K, so models of the chemistry of the (more abundant and more commonly studied) trans-alkenes may need to be modified to include the cis isomers of the corresponding allylic radicals. Addition of molecular oxygen to 1-methylallyl radical can form 2-butene-1-peroxy radical (CH(3)CH═CHCH(2)OO(?)), and quantum chemistry is used to thoroughly explore the possible unimolecular reactions of the cis and trans isomers of this radical. The cis isomer of the 2-butene-1-peroxy radical has the lowest barrier (via 1,6 H-shift) to further reaction, but this barrier appears to be too high to compete with loss of O(2).     

烷烃与氢过氧自由基氢提取反应类反应能垒与速率常数的精确计算

氢过氧自由基从烷烃分子中提取氢的反应是碳氢燃料中低温燃烧化学中非常重要的一类反应。本文用等键反应方法计算了这一类反应的动力学参数。所有反应物、过渡态、产物的几何结构均在HF/6-31+G(d)水平下优化得到。以反应中的过渡态反应中心的几何结构守恒为判据,该反应类可用等键反应处理。本文选取了乙烷和氢过氧自由基的氢提取反应为参考反应,其它反应作为目标反应,用等键反应方法对目标反应在HF/6-31+G(d)水平的近似能垒和反应速率常数进行了校正。为了验证方法的可靠性,选取C5以下的烷烃分子体系,对等键反应方法校正结果和高精度CCSD(T)/CBS直接计算结果进行了比较,最大绝对误差为5.58k J?mol~(-1),因此,采用等键反应方法只需用低水平HF从头算方法就可以再现高精度CCSD(T)/CBS计算结果,从而解决了该反应类中大分子体系的能垒的精确计算。本文的研究为碳氢化合物中低温燃烧模拟中重要的烷烃与氢过氧自由基氢提取反应提供了准确的动力学参数。     

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.     

State of the art theoretical study and comparison to experiment for the phenol...argon complex

The phenol...argon complex was studied by means of various high level ab initio quantum mechanics methods and high resolution threshold ionization spectroscopy. The structure and stabilization energy of different conformers were determined. Stabilization energy of van der Waals bonded and H-bonded PhOH...Ar complex determined at CCSD(T) complete basis set (CBS) level for CP-RI-MP2/cc-pVTZ/Ar aug-cc-pVTZ geometries amount to 434 and 285 cm(-1). The CCSD(T)/CBS were constructed either as a sum of MP2/CBS interaction energy and CCSD(T) correction term [difference between CCSD(T) and MP2 correlation energies determined with medium basis set] or directly from CCSD(T)/aug-cc-pVDZ and aug-cc-pVTZ energies. Both schemes provide very similar values. Harmonic vibrational analysis revealed that the H-bonded structure does not represent energy minimum but first order transition structure. The respective imaginary vibrational mode (16 cm(-1)) connects two possible argon locations -- above and below the phenol aromatic ring. Including the DeltaZPVE, we obtained stabilization enthalpy at 0 K of 389 cm(-1). This value is marginally higher (25-35 cm(-1), 0.07-0.10 kcal/mol) than the experimental value. The determination of DeltaZPVE constitutes the most significant error and possible improvements should come from more accurate evaluation of the (nonharmonic) vibrational frequencies.     

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3

Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y > or = 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin-orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to approximately 6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up to 3-4 kcal/mol. DFT using the B3LYP functional and 6-311+G** basis set predict isodesmic reaction energies to within a few kcal/mol compared with the CCSD(T)/CBS method so isodesmic reactions involving BN compounds and the analogous hydrocarbons can be used to estimate heats of formation. Heats of formation of c-B3N3H12 and c-B3N3H6 are -95.5 and -115.5 kcal/mol at 298 K, respectively, using our best calculated CCSD(T)/CBS approach. The experimental value for c-B3N3H6 appears to be approximately 7 kcal/mol too negative. Enthalpies, entropies, and free energies are calculated for many dehydrocoupling and dehydrogenation reactions that convert BNH6 to alicyclic and cyclic oligomers and H2(g). Generally, the reactions are highly exothermic and exergonic as well because of the release of 1 or more equivalents of H2(g). For c-B3N3H12 and c-B3N3H6, available experimental data for sublimation and vaporization lead to estimates of their condensed phase 298 K heats of formation: DeltaHf degrees [c-B3N3H12(s)] = -124 kcal/mol and DeltaHf degrees [c-B3N3H6(l)] = -123 kcal/mol. The reaction thermochemistries for the dehydrocoupling of BNH6(s) to c-B3N3H12(s) and the dehydrogenation of c-B3N3H12(s) to c-B3N3H6(l) are much less exothermic compared with the gas-phase reactions due to intermolecular forces which decrease in the order BNH6 > cyclo-B3N3H12 > cyclo-B3N3H6. The condensed phase reaction free energies are less negative compared with the gas-phase reactions but are still too favorable for BNH6 to be regenerated from either c-B3N3H12 or c-B3N3H6 by just an overpressure of H2.     

Benchmark calculations on the adiabatic ionization potentials of M-NH(3) (M=Na,Al,Ga,In,Cu,Ag)

The ground states of the M-NH(3) (M=Na,Al,Ga,In,Cu,Ag) complexes and their cations have been studied with density functional theory and coupled cluster [CCSD(T)] methods. The adiabatic ionization potentials (AIPs) of these complexes are calculated, and these are compared to results from high-resolution zero-electron kinetic energy photoelectron spectroscopy. By extrapolating the CCSD(T) energies to the complete basis set (CBS) limit and including the core-valence, scalar relativistic, spin-orbit, and zero-point corrections, the CCSD(T) method is shown to be able to predict the AIPs of these complexes to better than 6 meV or 0.15 kcal/mol. 27 exchange-correlation functionals, including one in the local density approximation, 13 in the generalized gradient approximation (GGA), and 13 with hybrid GGAs, were benchmarked in the calculations of the AIPs. The B1B95, mPW1PW91, B98, B97-1, PBE1PBE, O3LYP, TPSSh, and HCTH93 functionals give an average error of 0.1 eV for all the complexes studied, with the B98 functional alone yielding a maximum error of 0.1 eV. In addition, the calculated metal-ammonia harmonic stretching frequencies with the CCSD(T) method are in excellent agreement with their experimental values, whereas the B3LYP method tends to underestimate these stretching frequencies. The metal-ammonia binding energies were also calculated at the CCSD(T)/CBS level, and are in excellent agreement with the available experimental values considering the error limits, except for Ag-NH(3) and Ag(+)-NH(3), where the calculations predict stronger bond energies than measured by about 4 kcal/mol, just outside the experimental error bars of +/-3 kcal/mol.     

Ab initio based potential energy surfaces and Franck-Condon analysis of ionization thresholds of cyclic-C3H and linear-C3H

We report a Franck-Condon analysis in reduced dimensionality of the ionization thresholds of linear(l)-C3H and cyclic(c)-C3H using MP2-based potential energy surfaces and CCSD(T)/aug-cc-pVTZ calculations of electronic energies at selected geometries. The potential energy surfaces are fits to tens of thousands of MP2/aug-cc-pVTZ energies for the neutral and cation systems. These fits properly describe the invariance of the potential with respect to all permutations of the three C atoms. The realism of the potential surfaces is assessed by comparing stationary-point structures, energies, and normal-mode frequencies with previous high-level ab initio calculations. Several key vibrational modes in this ionization process are located at saddle points and so a numerical approach to obtain the Franck-Condon factors for those modes is done. On the basis of this analysis combined with a simple harmonic treatment of the energies of the remaining modes and key electronic energy differences obtained with CCSD(T)/aug-cc-pVTZ calculations, we find the threshold ionization energy of l-C3H to be 9.06 eV and for c-C3H we estimate the threshold to be in the range 9.70-9.76 eV. We estimate these values are accurate to within +/-0.05 eV.     

A comparison of orbital interactions in the additions of phosphonyl and acyl radicals to double bonds

Calculation of the barriers for addition of the H2P(=O) and HC(=O) radicals to alkenes, at the CCSD(T)/aug-cc-pVDZ//BHandHLYP/6-311G** level, indicates that both radicals display ambiphilic behaviour. For the HC(=O) radical this behaviour occurs because a secondary orbital interaction of the type pi*(C=O)<--HOMO acts in conjunction with the primary SOMO<--HOMO interaction to balance the SOMO-->LUMO interaction. For the H2P(=O) radical, on the other hand, the much higher-lying LUMO (the sigma*P-O orbital) allows for only minimal secondary interaction, and this radical's ambiphilic behaviour is therefore reflective of a balance between SOMO-->LUMO and SOMO<--HOMO interactions.