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.     

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

Theoretical prediction of the ionization energies of the C4H7 radicals: 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl radicals

The ionization energies (IEs) for the 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl 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 are included in these calculations. The present CCSD(T)/CBS results are then compared with the IEs determined in the photoelectron experiment by Schultz et al. [J. Am. Chem. Soc. 106, 7336 (1984)] The predicted IE value (7.881 eV) of 2-methylallyl radical is found to compare very favorably with the experimental value of 7.90+/-0.02 eV. Two ionization transitions for cis-1-methylallyl and trans-1-methylallyl radicals have been considered here. The comparison between the predicted IE values and the previous measurements shows that the photoelectron peak observed by Schultz et al. likely corresponds to the adiabatic ionization transition for the trans-1-methylallyl radical to form trans-1-methylallyl cation. Although a precise IE value for the cyclopropylmethyl radical has not been directly determined, the experimental value deduced indirectly using other known energetic data is found to be in good accord with the present CCSD(T)/CBS prediction. We expect that the Franck-Condon factor for ionization transition of c-C4H7-->bicyclobutonium is much less favorable than that for ionization transition of c-C4H7-->planar-C4H7+, and the observed IE in the previous photoelectron experiment is likely due to the ionization transition for c-C4H7-->planar-C4H7+. Based on our CCSD(T)/CBS prediction, the ionization transition of c-C4H7-->bicyclobutonium with an IE value around 6.92 eV should be taken as the adiabatic ionization transition for the cyclobutyl radical. The present study provides support for the conclusion that the CCSD(T)/CBS approach with high-level energetic corrections can be used to provide reliable IE predictions for C4 hydrocarbon radicals with an uncertainty of +/-22 meV. The CCSD(T)/CBS predictions to the heats of formation for the aforementioned radicals and cations are also presented.     

Structures, internal rotor potentials, and thermochemical properties for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals

Structures, enthalpy (Δ(f)H°(298)), entropy (S°(T)), and heat capacity (C(p)(T)) are determined for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals using the density functional B3LYP and composite CBS-QB3 calculations. Enthalpies of formation (Δ(f)H°(298)) are determined at the B3LYP/6-31G(d,p), B3LYP/6-31+G(2d,2p), and composite CBS-QB3 levels using several work reactions for each species. Entropy (S) and heat capacity (C(p)(T)) values from vibration, translational, and external rotational contributions are calculated using the rigid-rotor-harmonic-oscillator approximation based on the vibration frequencies and structures obtained from the density functional studies. Contribution to Δ(f)H(T), S, and C(p)(T) from the analysis on the internal rotors is included. Recommended values for enthalpies of formation of the most stable conformers of nitroacetone cc(═o)cno2, acetonitrite cc(═o)ono, nitroacetate cc(═o)no2, and acetyl nitrite cc(═o)ono are -51.6 kcal mol(-1), -51.3 kcal mol(-1), -45.4 kcal mol(-1), and -58.2 kcal mol(-1), respectively. The calculated Δ(f)H°(298) for nitroethylene c═cno2 is 7.6 kcal mol(-1) and for vinyl nitrite c═cono is 7.2 kcal mol(-1). We also found an unusual phenomena: an intramolecular transfer reaction (isomerization) with a low barrier (3.6 kcal mol(-1)) in the acetyl nitrite. The NO of the nitrite (R-ONO) in CH(3)C(═O')ONO moves to the C═O' oxygen in a motion of a stretching frequency and then a shift to the carbonyl oxygen (marked as O' for illustration purposes).     

Amorphous/crystalline (A/C) thermodynamic "rules of thumb": estimating standard thermodynamic data for amorphous materials using standard data for their crystalline counterparts

Standard thermochemical data (in the form of Δ(f)H° and Δ(f)G°) are available for crystalline (c) materials but rarely for their corresponding amorphous (a) counterparts. This paper establishes correlations between the sets of data for the two material forms (where known), which can then be used as a guideline for estimation of missing data. Accordingly, Δ(f)H°(a)/kJ mol(-1) ≈ 0.993Δ(f)H°(c)/kJ mol(-1) + 12.52 (R(2) = 0.9999; n = 50) and Δ(f)G°/kJ mol(-1) ≈ 0.988Δ(f)H°(c)/kJ mol(-1) + 0.70 (R(2) = 0.9999; n = 10). Much more tentatively, we propose that S°(298)(c)/J K(-1) mol(-1) ≈ 1.084S°(298)(c)/J K(-1) mol(-1) + 6.54 (R(2) = 0.9873; n = 11). An amorphous hydrate enthalpic version of the Difference Rule is also proposed (and tested) in the form [Δ(f)H°(M(p)X(q)·nH(2)O,a) - Δ(f)H°(M(p)X(q),a)]/kJ mol(-1) ≈ Θ(Hf)n ≈ -302.0n, where M(p)X(q)·nH(2)O represents an amorphous hydrate and M(p)X(q) the corresponding amorphous anhydrous parent salt.     

Thermochemical data and additivity group values for ten species of o-xylene low-temperature oxidation mechanism

o-Xylene could be a good candidate to represent the family of aromatic hydrocarbons in a surrogate fuel. This study uses computational chemistry to calculate standard enthalpies of formation at 298 K, Δ(f)H°(298 K), standard entropies at 298 K, S°(298 K), and standard heat capacities C(p)°(T) over the temperature range 300 K to 1500 K for ten target species present in the low-temperature oxidation mechanism of o-xylene: o-xylene (1), 2-methylbenzyl radical (2), 2-methylbenzylperoxy radical (3), 2-methylbenzyl hydroperoxide (4), 2-(hydroperoxymethyl)benzyl radical (5), 2-(hydroperoxymethyl)benzaldehyde (6), 1-ethyl-2-methylbenzene (7), 2,3-dimethylphenol (8), 2-hydroxybenzaldehyde (9), and 3-hydroxybenzaldehyde (10). Δ(f)H°(298 K) values are weighted averages across the values calculated using five isodesmic reactions and five composite calculation methods: CBS-QB3, G3B3, G3MP2, G3, and G4. The uncertainty in Δ(f)H°(298 K) is also evaluated. S°(298 K) and C(p)°(T) values are calculated at B3LYP/6-311G(d,p) level of theory from molecular properties and statistical thermodynamics through evaluation of translational, rotational, vibrational, and electronic partition functions. S°(298 K) and C(p)°(300 K) values are evaluated using the rigid-rotor-harmonic-oscillator model. C(p)°(T) values at T ≥ 400 K are calculated by treating separately internal rotation contributions and translational, external rotational, vibrational, and electronic contributions. The thermochemical properties of six target species are used to develop six new additivity groups taking into account the interaction between two substituents in ortho (ORT/CH2OOH/ME, ORT/ET/ME, ORT/CHO/OH, ORT/CHO/CH2OOH) or meta (MET/CHO/OH) positions, and the interaction between three substituents (ME/ME/OH123) located one beside the other (positions numbered 1, 2, 3) for two- or three-substituted benzenic species. Two other additivity groups are also developed using the thermochemical properties of benzenic species taken from the literature: the C/CB/H2/OO and the CB/CO groups. These groups extend the capacities of the group additivity method to deal with substituted benzenic species.     

On the nature of stabilization in weak, medium, and strong charge-transfer complexes: CCSD(T)/CBS and SAPT calculations

Weak, medium, and strong charge-transfer (CT) complexes containing various electron donors (C(2)H(4), C(2)H(2), NH(3), NMe(3), HCN, H(2)O) and acceptors (F(2), Cl(2), BH(3), SO(2)) were investigated at the CCSD(T)/complete basis set (CBS) limit. The nature of the stabilization for these CT complexes was evaluated on the basis of perturbative NBO calculations and DFT-SAPT/CBS calculations. The structure of all of the complexes was determined by the counterpoise-corrected gradient optimization performed at the MP2/cc-pVTZ level, and most of complexes possess a linear-like contact structure. The total stabilization energies lie between 1 and 55 kcal/mol and the strongest complexes contain BH(3) as an electron acceptor. When ordering the electron donors and electron acceptors on the basis of these energies, we obtain the same order as that based on the perturbative E2 charge-transfer energies, which provides evidence that the charge-transfer term is the dominant energy contribution. The CCSD(T) correction term, defined as the difference between the CCSD(T) and MP2 interaction energies, is mostly small, which allows the investigation of the CT complexes of this type at the "cheap" MP2/CBS level. In the case of weak and medium CT complexes (with stabilization energy smaller than about 15 kcal/mol), the dominant stabilization originates in the electrostatic term; the dispersion as well as induction and δ(HF) terms covering the CT energy contribution are, however, important as well. For strong CT complexes, induction energy is the second (after electrostatic) most important energy term. The role of the induction and δ(HF) terms is unique and characteristic for CT complexes. For all CT complexes, the CCSD(T)/CBS and DFT-SAPT/CBS stabilization energies are comparable, and surprisingly, it is true even for very strong CT complexes with stabilization energy close to 50 kcal/mol characteristic by substantial charge transfer (more than 0.3 e). It is thus possible to conclude that perturbative DFT-SAPT analysis is robust enough to be applied even for dative-like complexes with substantial charge transfer.     

Ab initio chemical kinetic study on the reactions of ClO with C2H2 and C2H4

The mechanisms for the reactions of ClO with C(2)H(2) and C(2)H(4) have been investigated at the CCSD(T)/CBS level of theory. The results show that in both systems, the interaction between the Cl atom of the ClO radical and the triple and double bonds of C(2)H(2) and C(2)H(4) forms prereaction van der Waals complexes with the O-Cl bond pointing perpendicularly toward the π-bonds, both with 2.1 kcal/mol binding energies. The mechanism is similar to those of the HO-C(2)H(2)/C(2)H(4) systems. The rate constants for the low energy channels have been predicted by statistical theory. For the reaction of ClO and C(2)H(2), the main channels are the production of CH(2)CO + Cl (k(1a)) and CHCO + HCl (k(1b)), with k(1a) = 1.19 × 10(-15)T(1.18) exp(-5814/T) and k(1b) = 6.94 × 10(-21) × T(2.60) exp(-6587/T) cm(3) molecule(-1) s(-1). For the ClO + C(2)H(4) reaction, the main pathway leads to C(2)H(4)O + Cl (k(2a)) with the predicted rate constant k(2a) = 2.13 × 10(-17)T(1.52) exp(-3849/T) in the temperature range of 300-3000 K. These rate constants are pressure-independent below 100 atm.     

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.     

Combined vacuum ultraviolet laser and synchrotron pulsed field ionization study of CH2BrCl

The pulsed field ionization-photoelectron (PFI-PE) spectrum of bromochloromethane (CH2BrCl) in the region of 85,320-88,200 cm-1 has been measured using vacuum ultraviolet laser. The vibrational structure resolved in the PFI-PE spectrum was assigned based on ab initio quantum chemical calculations and Franck-Condon factor predictions. At energies 0-1400 cm-1 above the adiabatic ionization energy (IE) of CH2BrCl, the Br-C-Cl bending vibration progression (nu1+=0-8) of CH2BrCl+ is well resolved and constitutes the major structure in the PFI-PE spectrum, whereas the spectrum at energies 1400-2600 cm-1 above the IE(CH2BrCl) is found to exhibit complex vibrational features, suggesting perturbation by the low lying excited CH2BrCl+(A 2A") state. The assignment of the PFI-PE vibrational bands gives the IE(CH2BrCl)=85,612.4+/-2.0 cm-1 (10.6146+/-0.0003 eV) and the bending frequencies nu1+(a1')=209.7+/-2.0 cm-1 for CH2BrCl+(X2A'). We have also examined the dissociative photoionization process, CH2BrCl+hnu-->CH2Cl++Br+e-, in the energy range of 11.36-11.57 eV using the synchrotron based PFI-PE-photoion coincidence method, yielding the 0 K threshold or appearance energy AE(CH2Cl+)=11.509+/-0.002 eV. Combining the 0 K AE(CH2Cl+) and IE(CH2BrCl) values obtained in this study, together with the known IE(CH2Cl), we have determined the 0 K bond dissociation energies (D0) for CH2Cl+-Br (0.894+/-0.002 eV) and CH2Cl-Br (2.76+/-0.01 eV). We have also performed CCSD(T, full)/complete basis set (CBS) calculations with high-level corrections for the predictions of the IE(CH2BrCl), AE(CH2Cl+), IE(CH2Cl), D0(CH2Cl+-Br), and D0(CH2Cl-Br). The comparison between the theoretical predictions and experimental determinations indicates that the CCSD(T, full)/CBS calculations with high-level corrections are highly reliable with estimated error limits of <17 meV.     

Thermodynamic properties of molecular borane phosphines, alane amines, and phosphine alanes and the [BH(4)(-)][PH(4)(+)], [AlH(4)(-)][NH(4)(+)], and [AlH(4)(-)][PH(4)(+)] salts for chemical hydrogen storage systems from ab initio electronic structure theory

The heats of formation for the molecules BH(3)PH(3), BH(2)PH(2), HBPH, AlH(3)NH(3), AlH(2)NH(2), HAlNH, AlH(3)PH(3), AlH(2)PH(2), HAlPH, AlH(4)(-), PH(3), PH(4), and PH(4)(+), as well as the diatomics BP, AlN, and AlP, have been calculated by using ab initio molecular orbital theory. The coupled cluster with single and double excitations and perturbative triples method (CCSD(T)) was employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-zeta, to extrapolate to the complete basis set limit. Additional d core functions were used for Al and P. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation of the salts [BH(4)(-)][PH(4)(+)](s), [AlH(4)(-)][NH(4)(+)](s), and [AlH(4)(-)][PH(4)(+)](s) have been estimated by using an empirical expression for the lattice energy and the calculated heats of formation of the two component ions. The calculations show that both AlH(3)NH(3)(g) and [AlH(4)(-)][NH(4)(+)](s) can serve as good hydrogen storage systems that release H(2) in a slightly exothermic process. In addition, AlH(3)PH(3) and the salts [AlH(4)(-)][PH(4)(+)] and [BH(4)(-)][PH(4)(+)] have the potential to serve as H(2) storage systems. The hydride affinity of AlH(3) is calculated to be -70.4 kcal/mol at 298 K. The proton affinity of PH(3) is calculated to be 187.8 kcal/mol at 298 K in excellent agreement with the experimental value of 188 kcal/mol. PH(4) is calculated to be barely stable with respect to loss of a hydrogen to form PH(3).     

Infrared spectrum of NH4+(H2O): evidence for mode specific fragmentation

The gas phase infrared spectrum (3250-3810 cm-1) of the singly hydrated ammonium ion, NH4+(H2O), has been recorded by action spectroscopy of mass selected and isolated ions. The four bands obtained are assigned to N-H stretching modes and to O-H stretching modes. The N-H stretching modes observed are blueshifted with respect to the corresponding modes of the free NH4+ ion, whereas a redshift is observed with respect to the modes of the free NH3 molecule. The O-H stretching modes observed are redshifted when compared to the free H2O molecule. The asymmetric stretching modes give rise to rotationally resolved perpendicular transitions. The K-type equidistant rotational spacings of 11.1(2) cm-1 (NH4+) and 29(3) cm-1 (H2O) deviate systematically from the corresponding values of the free molecules, a fact which is rationalized in terms of a symmetric top analysis. The relative band intensities recorded compare favorably with predictions of high level ab initio calculations, except on the nu3(H2O) band for which the observed value is about 20 times weaker than the calculated one. The nu3(H2O)/nu1(H2O) intensity ratios from other published action spectra in other cationic complexes vary such that the nu3(H2O) intensities become smaller the stronger the complexes are bound. The recorded ratios vary, in particular, among the data collected from action spectra that were recorded with and without rare gas tagging. The calculated anharmonic coupling constants in NH4+(H2O) further suggest that the coupling of the nu3(H2O) and nu1(H2O) modes to other cluster modes indeed varies by orders of magnitude. These findings together render a picture of a mode specific fragmentation dynamic that modulates band intensities in action spectra with respect to absorption spectra. Additional high level electronic structure calculations at the coupled-cluster singles and doubles with a perturbative treatment of triple excitations [CCSD(T)] level of theory with large basis sets allow for the determination of an accurate binding energy and enthalpy of the NH4+(H2O) cluster. The authors' extrapolated values at the CCSD(T) complete basis set limit are De [NH4+-(H2O)]=-85.40(+/-0.24) kJ/mol and DeltaH(298 K) [NH4+-(H2O)]=-78.3(+/-0.3) kJ/mol (CC2), in which double standard deviations are indicated in parentheses.     

Oxidative addition of the Cα-Cβ bond in β-O-4 linkage of lignin to transition metals using a relativistic pseudopotential-based ccCA-ONIOM method

A multi-level multi-layer QM/QM method, the relativistic pseudopotential correlation-consistent composite approach within an ONIOM framework (rp-ccCA-ONIOM), was applied to study the oxidative addition of the C(α)-C(β) bond in an archetypal arylglycerol β-aryl ether (β-O-4 linkage) substructure of lignin to Ni, Cu, Pd and Pt transition metal atoms. The chemically active high-level layer is treated using the relativistic pseudopotential correlation-consistent composite approach (rp-ccCA), an efficient methodology designed to reproduce an accuracy that would be obtained using the more computationally demanding CCSD(T)/aug-cc-pCV∞Z-PP, albeit at a significantly reduced computational cost, while the low-level layer is computed using B3LYP/cc-pVTZ. The thermodynamic and kinetic feasibilities of the model reactions are reported in terms of enthalpies of reactions at 298 K (ΔH°(298)) and activation energies (ΔH-act). The results obtained from the rp-ccCA:B3LYP hybrid method are compared to the corresponding values using CCSD(T) and several density functionals including B3LYP, M06, M06 L, B2PLYP, mPWPLYP and B2GP-PLYP. The energetics of the oxidative addition of C?C bond in ethane to Ni, Cu, Pd and Pt atoms are also reported to demonstrate that the rp-ccCA method effectively reproduces the accuracy of the CCSD(T)/aug-cc-pCV∞Z method. Our results show that in the catalytic activation of the C(α)-C(β) bond of β-O-4, the use of platinum metal catalysts will lead to the most thermodynamically favored reaction with the lowest activation barrier.     

Thermodynamic properties of molecular borane amines and the [BH4-][NH4+] salt for chemical hydrogen storage systems from ab initio electronic structure theory

The heats of formation for the borane amines BH3NH3, BH2NH2, and HBNH, tetrahedral BH4-, and the BN molecule have been calculated by using ab initio molecular orbital theory. Coupled cluster calculations with single and double excitations and perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-zeta, to extrapolate to the complete basis set limit. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation (in kcal/mol) at 0 K in the gas phase are Delta Hf(BH3NH3) = -9.1, Delta Hf(BH2NH2) = -15.9, Delta Hf(BHNH) = 13.6, Delta Hf(BN) = 146.4, and Delta Hf(BH4-) = -11.6. The reported experimental value for Delta Hf(BN) is clearly in error. The heat of formation of the salt [BH4-][NH4+](s) has been estimated by using an empirical expression for the lattice energy and the calculated heats of formation of the two component ions. The calculations show that both BH3NH3(g) and [BH4-][NH4+](s) can serve as good hydrogen storage systems which release H2 in a slightly exothermic process. The hydride affinity of BH3 is calculated to be 72.2 kcal/mol, in excellent agreement with the experimental value at 298 K of 74.2 +/- 2.8 kcal/mol.     

Structures, electron affinities, and harmonic vibrational frequencies of C6H5X/C6H5X- (X = N, S, NH, PH, CH2, and SiH2)

The molecular structures and electron affinities of the C6H5X/C6H5X- (X = N, S, NH, PH, CH2, and SiH2) species have been determined using seven different density functional or hybrid Hartree-Fock density functional methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. These methods have been carefully calibrated (Chem. Rev. 2002, 102, 231). The geometries are fully optimized with each density functional theory (DFT) method, and discussed. Harmonic vibrational frequencies were found to be within 3.2% of available experimental values for most functionals. Three different types of the neutral-anion energy separations reported in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The most reliable adiabatic electron affinities, obtained at the DZP++ BPW91 level of theory, are 1.45 (C6H5N), 2.29 (C6H5S), 1.57 (C6H5NH), 1.51 (C6H5PH), 0.91 (C6H5CH2), and 1.48 eV (C6H5SiH2), respectively. Compared with the experimental values, the average absolute error of the BPW91 method is 0.04 eV. The B3LYP and B3PW91 functionals also gave excellent predictions, with average absolute errors of 0.06 and 0.07 eV, respectively.     

Heats of formation of diphosphene, phosphinophosphinidene, diphosphine, and their methyl derivatives, and mechanism of the borane-assisted hydrogen release

The heats of formation of diphosphene (cis- and trans-P2H2), phopshinophosphinidene (singlet and triplet H2PP) and diphosphine (P2H4), as well as those of the P2H and P2H3 radicals resulting from PH bond cleavages, have been calculated by using high-level ab initio electronic structure theory. Energies were calculated using coupled-cluster theory with a perturbative treatment for triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets with additional tight d-functions on P (aug-cc-pV(n+d)Z) up to quadruple- or quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries and vibrational frequencies were determined with the CCSD(T) method. Core-valence and scalar relativistic corrections were included, as well as scaled zero-point energies. We find the following heats of formation (kcal/mol) at 298 [0] K: DeltaH(degree)(f)(P2H) = 53.4 [54.4]; DeltaH(degree)(f)(cis-P2H2) = 32.0 [33.9]; DeltaH(degree)(f)(trans-P2H2) = 28.7 [30.6]; DeltaH(degree)(f)(H2PP) = 53.7 [55.6]; DeltaH(degree)(f)(3H2PP) = 56.5 [58.3]; DeltaH(degree)(f)(P2H3) = 32.3 [34.8]; DeltaH(degree)(f)(P2H4) = 5.7 [9.1] (expt, 5.0 +/- 1.0 at 298 K); and DeltaH(degree)(f)(CH3PH2) = -5.0 [-1.4]. We estimate these values to have an accuracy of +/-1.0 kcal/mol. In contrast to earlier results, we found a singlet ground state for phosphinophosphinidene (H2PP) with a singlet-triplet energy gap of 2.8 kcal/mol. We calculated the heats of formation of the methylated derivatives CH3PPH, CH3HPPH2, CH3PPCH3, CH3HPP, (CH3)2PP, (CH3)2PPH2, and CH3HPPHCH3 by using isodesmic reactions at the MP2/CBS level. The calculated results for the hydrogenation reactions RPPR + H2 --> RHPPHR and R2PP + H2 --> R2PPH2 show that substitution of an organic substituent for H improves the energetics, suggesting that secondary diphosphines and diphosphenes are potential candidates for use in a chemical hydrogen storage system. A comparison with the nitrogen analogues is given. The mechanism for H2-generation from diphosphine without and with BH3 as a catalyst was examined. Including tunneling corrections, the rate constant for the catalyzed reaction is 4.5 x 1015 times faster than the uncatalyzed result starting from separated catalyst and PH2PH2.     

Photoion photoelectron coincidence spectroscopy of primary amines RCH2NH2 (R = H, CH3, C2H5, C3H7, i-C3H7): alkylamine and alkyl radical heats of formation by isodesmic reaction networks

Alkylamines (RCH(2)NH(2), R = H, CH(3), C(2)H(5), C(3)H(7), i-C(3)H(7)) have been investigated by dissociative photoionization by threshold photoelectron photoion coincidence spectroscopy (TPEPICO). The 0 K dissociation limits (9.754 +/- 0.008, 9.721 +/- 0.008, 9.702 +/- 0.012, and 9.668 +/- 0.012 eV for R = CH(3), C(2)H(5), C(3)H(7), i-C(3)H(7), respectively) have been determined by preparing energy-selected ions and collecting the fractional abundances of parent and daughter ions. All alkylamine cations produce the methylenimmonium ion, CH(2)NH(2)+, and the corresponding alkyl free radical. Two isodesmic reaction networks have also been constructed. The first one consists of the alkylamine parent molecules, and the other of the alkyl radical photofragments. Reaction heats within the isodesmic networks have been calculated at the CBS-APNO and W1U levels of theory. The two networks are connected by the TPEPICO dissociation energies. The heats of formation of the amines and the alkyl free radicals are then obtained by a modified least-squares fit to minimize the discrepancy between the TPEPICO and the ab initio values. The analysis of the fit reveals that the previous experimental heats of formation are largely accurate, but certain revisions are suggested. Thus, Delta(f)Ho(298K)[CH(3)NH(2)(g)] = -21.8 +/- 1.5 kJ mol-1, Delta(f)Ho(298K)[C(2)H(5)NH(2)(g)] = -50.1 +/- 1.5 kJ mol(-1), Delta(f)Ho(298K)[C(3)H(7)NH(2)(g)] = -70.8 +/- 1.5 kJ mol(-1), Delta(f)Ho(298K)[C(3)H(7)*] = 101.3 +/- 1 kJ mol(-1), and Delta(f)Ho(298K)[i-C(3)H(7)*] = 88.5 +/- 1 kJ mol(-1). The TPEPICO and the ab initio results for butylamine do not agree within 1 kJ mol-1; therefore, no new heat of formation is proposed for butylamine. It is nevertheless indicated that the previous experimental heats of formation of methylamine, propylamine, butylamine, and isobutylamine may have been systematically underestimated. On the other hand, the error in the ethyl radical heat of formation is found to be overestimated and can be decreased to +/- 1 kJ mol(-1); thus, Delta(f)Ho(298K)[C(2)H(5).] = 120.7 +/- 1 kJ mol(-1). On the basis of the data analysis, the heat of formation of the methylenimmonium ion is confirmed to be Delta(f)Ho(298K)[CH(2)NH(2)+] = 750.3 +/- 1 kJ mol(-1).     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

The heats of formation of diazene, hydrazine, N2H3+, N2H5+, N2H, and N2H3 and the Methyl Derivatives CH3NNH, CH3NNCH3, and CH3HNNHCH3

The heats of formation of N(2)H, diazene (cis- and trans-N(2)H(2)), N(2)H(3), and hydrazine (N(2)H(4)), as well as their protonated species (diazenium, N(2)H(3)(+), and hydrazinium, N(2)H(5)(+)), have been calculated by using high level electronic structure theory. Energies were calculated by using coupled cluster theory with a perturbative treatment of the triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets (aug-cc-pVnZ) up to quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries were optimized at the CCSD(T) level with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. Core-valence and scalar relativistic corrections were included, as well as scaled zero point energies. We find the following heats of formation (kcal/mol) at 0 (298) K: DeltaH(f)(N(2)H) = 60.8 (60.1); DeltaH(f)(cis-N(2)H(2)) = 54.9 (53.2); DeltaH(f)(trans-N(2)H(2)) = 49.9 (48.1) versus >/=48.8 +/- 0.5 (exptl, 0 K); DeltaH(f)(N(2)H(4)) = 26.6 (23.1) versus 22.8 +/- 0.2 (exptl, 298 K); DeltaH(f)(N(2)H(3)) = 56.2 (53.6); DeltaH(f)(N(2)H(3)(+)) = 231.6 (228.9); and DeltaH(f)(N(2)H(5)(+)) = 187.1 (182.7). In addition, we calculated the heats of formation of CH(3)NH(2), CH(3)NNH, and CH(3)HNNHCH(3) by using isodesmic reactions and at the G3(MP2) level. The calculated results for the hydrogenation reaction RNNR + H(2) --> RHNNHR show that substitution of an organic substituent for H improved the energetics, suggesting that these types of compounds may be possible to use in a chemical hydrogen storage system.     

Electronic structure and thermochemical properties of silicon-doped lithium clusters Li(n)Si0/+, n = 1-8: new insights on their stability

A theoretical investigation on small silicon-doped lithium clusters Li(n)Si with n = 1-8, in both neutral and cationic states is performed using the high accuracy CCSD(T)/complete basis set (CBS) method. Location of the global minima is carried out using a stochastic search method and the growth pattern of the clusters emerges as follows: (i) the species Li(n)Si with n ≤ 6 are formed by directly binding one Li to a Si of the smaller cluster Li(n-1)Si, (ii) the structures tend to have an as high as possible symmetry and to maximize the coordination number of silicon. The first three-dimensional global minimum is found for Li(4)Si, and (iii) for Li(7)Si and Li(8)Si, the global minima are formed by capping Li atoms on triangular faces of Li(6)Si (O(h)). A maximum coordination number of silicon is found to be 6 for the global minima, and structures with higher coordination of silicon exist but are less stable. Heats of formation at 0 K (Δ(f)H(0)) and 298 K (Δ(f)H(298)), average binding energies (E(b)), adiabatic (AIE) and vertical (VIE) ionization energies, dissociation energies (D(e)), and second-order difference in total energy (Δ(2)E) of the clusters in both neutral and cationic states are calculated from the CCSD(T)/CBS energies and used to evaluate the relative stability of clusters. The species Li(4)Si, Li(6)Si, and Li(5)Si(+) are the more stable systems with large HOMO-LUMO gaps, E(b), and Δ(2)E. Their enhanced stability can be rationalized using a modified phenomenological shell model, which includes the effects of additional factors such as geometrical symmetry and coordination number of the dopant. The new model is subsequently applied with consistency to other impure clusters Li(n)X with X = B, Al, C, Si, Ge, and Sn.