Enthalpies of formation of gas-phase N3, N3-, N5+, and N5- from Ab initio molecular orbital theory, stability predictions for N5(+)N3(-) and N5(+)N5(-), and experimental evidence for the instability of N5(+)N3(-)

Ab initio molecular orbital theory has been used to calculate accurate enthalpies of formation and adiabatic electron affinities or ionization potentials for N3, N3-, N5+, and N5- from total atomization energies. The calculated heats of formation of the gas-phase molecules/ions at 0 K are DeltaHf(N3(2Pi)) = 109.2, DeltaHf(N3-(1sigma+)) = 47.4, DeltaHf(N5-(1A1')) = 62.3, and DeltaHf(N5+(1A1)) = 353.3 kcal/mol with an estimated error bar of +/-1 kcal/mol. For comparison purposes, the error in the calculated bond energy for N2 is 0.72 kcal/mol. Born-Haber cycle calculations, using estimated lattice energies and the adiabatic ionization potentials of the anions and electron affinities of the cations, enable reliable stability predictions for the hypothetical N5(+)N3(-) and N5(+)N5(-) salts. The calculations show that neither salt can be stabilized and that both should decompose spontaneously into N3 radicals and N2. This conclusion was experimentally confirmed for the N5(+)N3(-) salt by low-temperature metathetical reactions between N5SbF6 and alkali metal azides in different solvents, resulting in violent reactions with spontaneous nitrogen evolution. It is emphasized that one needs to use adiabatic ionization potentials and electron affinities instead of vertical potentials and affinities for salt stability predictions when the formed radicals are not vibrationally stable. This is the case for the N5 radicals where the energy difference between vertical and adiabatic potentials amounts to about 100 kcal/mol per N5.     

Computational s-block thermochemistry with the correlation consistent composite approach

The correlation consistent composite approach (ccCA) is a model chemistry that has been shown to accurately compute gas-phase enthalpies of formation for alkali and alkaline earth metal oxides and hydroxides (Ho, D. S.; DeYonker, N. J.; Wilson, A. K.; Cundari, T. R. J. Phys. Chem. A 2006, 110, 9767).The ccCA results contrast to more widely used model chemistries where calculated enthalpies of formation for such species can be in error by up to 90 kcal mol-1. In this study, we have applied ccCA to a more general set of 42 s-block molecules and compared the ccCA DeltaHf values to values obtained using the G3 and G3B model chemistries. Included in this training set are water complexes such as Na(H2O)n+ where n = 1 - 4, dimers and trimers of ionic compounds such as (LiCl)2 and (LiCl)3, and the largest ccCA computation to date: Be(acac)2, BeC10H14O4. Problems with the G3 model chemistries seem to be isolated to metal-oxygen bonded systems and Be-containing systems, as G3 and G3B still perform quite well with a 2.7 and 2.6 kcal mol-1 mean absolute deviation (MAD), respectively, for gas-phase enthalpies of formation. The MAD of the ccCA is only 2.2 kcal mol-1 for enthalpies of formation (DeltaHf) for all compounds studied herein. While this MAD is roughly double that found for a ccCA study of >350 main group (i.e., p-block) compounds, it is commensurate with typical experimental uncertainties for s-block complexes. Some molecules where G3/G3B and ccCA computed DeltaHf values deviate significantly from experiment, such as (LiCl)3, NaCN, and MgF, are inviting candidates for new experimental and high-level theoretical studies.     

The correlation-consistent composite approach: application to the G3/99 test set

The correlation-consistent composite approach (ccCA), an ab initio composite technique for computing atomic and molecular energies, recently has been shown to successfully reproduce experimental data for a number of systems. The ccCA is applied to the G3/99 test set, which includes 223 enthalpies of formation, 88 adiabatic ionization potentials, 58 adiabatic electron affinities, and 8 adiabatic proton affinities. Improvements on the original ccCA formalism include replacing the small basis set quadratic configuration interaction computation with a coupled cluster computation, employing a correction for scalar relativistic effects, utilizing the tight-d forms of the second-row correlation-consistent basis sets, and revisiting the basis set chosen for geometry optimization. With two types of complete basis set extrapolation of MP2 energies, ccCA results in an almost zero mean deviation for the G3/99 set (with a best value of -0.10 kcal mol(-1)), and a 0.96 kcal mol(-1) mean absolute deviation, which is equivalent to the accuracy of the G3X model chemistry. There are no optimized or empirical parameters included in the computation of ccCA energies. Except for a few systems to be discussed, ccCA performs as well as or better than Gn methods for most systems containing first-row atoms, while for systems containing second-row atoms, ccCA is an improvement over Gn model chemistries.     

DFT calculations of the ionization potentials and electron affinities of serinamide

Adiabatic and vertical ionization potentials (IPs) and valence electron affinities (EAs) of serinamide in the gas phase have been determined using density functional theory (DFT) B3LYP, B3P86, and B3PW91 methods with the 6‐311++G** and 6‐311G** basis sets, respectively. IPs and EAs of serinamide in solution have been calculated with the B3LYP method using the 6‐311++G** and 6‐311G** basis sets. Eight possible conformers of serinamide and its charged states in the gas phase have been optimized employing the DFT B3LYP method with 6‐311++G** and 6‐311G** basis sets, respectively. All the adiabatic and vertical ionization potentials (AIPs and VIPs) of eight serinamide conformers in our work are positive values, whether in the gas phase or in solutions; the IPs in solutions are smaller than the results in the gas phase and decrease with increased dielectric constants in solutions. This finding indicates that the cationic states in solutions are more stable than those in the gas phase. All EAs of eight serinamide conformers are negative values in the gas phase, indicating that the anionic states are unstable with respect to electron autodetachment, both adiabatically and vertically. In contrast, all other adiabatic electron affinities (AEAs) are negative values in solutions except for 6S in water; 7S in chloroform, acetone, and water; and 8S in acetone and water, and increase with increasing of dielectric constants in solutions. All vertical electron affinities (VEAs) are negative values in solutions; however, no good rule has been found for these values in solutions. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Accurate heats of formation of the "Arduengo-type" carbene and various adducts including H2 from ab initio molecular orbital theory

The heats of formation of saturated and unsaturated diaminocarbenes (imadazol(in)-2-ylidenes) have been calculated by using high levels of ab initio electronic structure theory. The calculations were done at the coupled cluster level through noniterative triple excitations with augmented correlation consistent basis sets up through quadruple. In addition, four other corrections were applied to the frozen core atomization energies: (1) a zero point vibrational correction; (2) a core/valence correlation correction; (3) a scalar relativistic correction; (4) a first-order atomic spin-orbit correction. The value of DeltaHf( 298) for the unsaturated carbene 1 is calculated to be 56.4 kcal/mol. The value of DeltaHf( 298) for the unsaturated triplet carbene (3)1 is calculated to be 142.8 kcal/mol, giving a singlet-triplet splitting of 86.4 kcal/mol. Addition of a proton to 1 forms 3 with DeltaHf( 298)(3) = 171.6 kcal/mol with a proton affinity for 1 of 250.5 kcal/mol at 298 K. Addition of a hydrogen atom to 1 forms 4 with DeltaHf( 298)(4) = 72.7 kcal/mol and a C-H bond energy of 35.8 kcal/mol at 298 K. Addition of H- to 1 gives 5 with DeltaHf( 298)(5) = 81.2 kcal/mol and 5 is not stable with respect to loss of an electron to form 4. Addition of H2 to the carbene center forms 6 with DeltaHf( 298)(6) = 41.5 kcal/mol and a heat of hydrogenation at 298 K of -14.9 kcal/mol. The value of DeltaHf( 298) for the saturated carbene 7 (obtained by adding H2 to the C=C bond of 1) is 47.4 kcal/mol. Hydrogenation of 7 to form the fully saturated imidazolidine, 8, gives DeltaHf( 298)(8) = 14.8 kcal/mol and a heat of hydrogenation at 298 K of -32.6 kcal/mol. The estimated error bars for the calculated heats of formation are +/-1.0 kcal/mol.     

Heats of formation of boron hydride anions and dianions and their ammonium salts [BnHmy-][NH4+]y with y=1-2

Thermochemical parameters of the closo boron hydride BnHn2- dianions, with n=5-12, the B3H8- and B11H14- anions, and the B5H9 and B10H14 neutral species were predicted by high-level ab initio electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T)/complete basis set (CBS) extrapolated energies, plus additional corrections were used to predict the heats of formation of the simplest BnHmy- species in the gas phase in kcal/mol at 298 K: DeltaHf(B3H8-)=-23.1+/-1.0; DeltaHf(B5H52-)=119.4+/-1.5; DeltaHf(B6H62-)=64.1+/-1.5; and DeltaHf(B5H9)=24.1+/-1.5. The heats of formation of the larger species were evaluated by the G3 method from hydrogenation reactions (values at 298 K, in kcal/mol with estimated error bars of+/-3 kcal/mol): DeltaHf(B7H72-)=51.8; DeltaHf(B8H82-)=46.1; DeltaHf(B9H92-)=24.4; DeltaHf(B10H102-)=-12.5; DeltaHf(B11H112-)=-11.8; DeltaHf(B12H122-)=-86.3; DeltaHf(B11H14-)=-57.3; and DeltaHf(B10H14)=18.7. A linear correlation between atomization energies of the dianions and energies of the BH units was found. The heats of formation of the ammonium salts of the anions and dianions were predicted using lattice energies (UL) calculated from an empirical expression based on ionic volumes. The UL values (0 K) of the BnHn2- dianions range from 319 to 372 kcal/mol. The values of UL for the B3H8- and B11H14- anions are 113 and 135 kcal/mol, respectively. The calculated lattice energies and gas-phase heats of formation of the constituent ions were used to predict the heats of formation of the ammonium crystal salts [BnHmy-][NH4+]y. These results were used to evaluate the thermodynamics of the H2 release reactions from the ammonium hydro-borate salts.     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Theoretical study of the enthalpies of formation for C40H56 carotenes

Theoretical study of the enthalpies of formation (DeltaHf) for polyenes up to nine ethylene units and for several C40H56 carotenes including beta-carotene, alpha-carotene, lycopene, and prolycopene is presented. For polyenes and small branched alkenes, we used G2, G3, and G3MP2B3 theories, and the DeltaHf values were evaluated with the atomization, isodesmic bond separation, and homodesmic schemes. The applicability of six DFT functionals were evaluated by comparing their predictions with those obtained using G3 theory within the atomization scheme. Additivity approaches, including atom equivalents and group equivalents using DFT and semiempirical theories, were explored. We found that group equivalents associated with isodesmic reactions are able to provide the most accurate predictions within the test set. The predictions from the six functionals are in good agreement with the G3 results. Among them, B3LYP performs the best, with an average absolute deviation of only 0.30 kcal/mol. The application of DFT in the prediction for the DeltaHf value of C40H56 carotenes is promising.     

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.     

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

Accurate and efficient method for predicting thermochemistry of polycyclic aromatic hydrocarbons - bond-centered group additivity

A self-consistent estimation method for the thermochemical properties of polycyclic aromatic hydrocarbons (PAH) is presented. This method is based on enthalpies of formation (DeltaHf(degrees), entropies (S(degrees)298, and heat capacities (C(degrees)p obtained from B3LYP/6-31G(d) calculations of the total energies and frequencies for 139 PAHs, including C(60) and C(70) fullerenes. The enthalpies of formation were calculated using an optimized set of homodesmic reactions given the available experimental DeltaHf(degrees) of PAHs. The theoretical entropies were compared with the existing experimental entropies, and some inconsistencies in the experimental data were identified. The estimation method presented here is a systematic extension of the widely employed atom-centered group additivity method originally proposed by Benson. This new method is based on bond-centered groups that define bonds linking two atom-centered groups and specify the size of the rings to which they belong. In addition, a term to describe the resonance energy is included. The thermochemical properties of PAHs up to C(70) fullerene are estimated with a mean average deviation of 2.8 kcal mol(-1) in DeltaHf(degrees), 0.7 cal K(-1) mol(-1) in S(degrees)298, and about 0.5 cal K(-1) mol(-1) in the C(degrees)p. This bond-centered group additivity method for the thermochemical properties of PAHs significantly expands both the range of systems that can be estimated and the accuracy of the estimations. The results of this work also allow us to assess the quality of available experimental data. For example, there are strong indications that the literature DeltaHf(degrees)of benzo[k]fluoranthene is about 10 kcal mol(-1) too low.     

DFT study of ionization potentials and electron affinities of formamide and its methylation derivatives

The ionization potentials (IPs) and electron affinities (EAs) of formamide in the gas phase have been calculated using density functional theory (DFT), ab initio HF and Møller-Plesset perturbational theory (MP) at 6-311++G** basis set. The results indicate that the IPs of formamide obtained with DFT and MP are in agreement with the results obtained from experiment. And B3LYP has been confirmed to be the most accurate method in calculating the AIPs and VIPs of formamide through our work. IPs and EAs of formamide in solution are not known experimentally, therefore IPs and EAs of formamide in chloroform, acetone, and dimethylsulfoxide have been calculated using polarized continuum model (PCM) with B3LYP/6-311++G** level and have been compared with the values in the gas phase. The AIPs and VIPs of formamide have been compared with those of its methylation derivatives. All EAs of methylation derivatives of formamide are bigger than those of formamide conformers in the gas phase with BLYP, B3LYP, and B3P86 methods at 6-311++G** basis set. All these indicate that all anions of methylation derivatives of formamide are more stable than anions of formamide with respect to electron detachment adiabatically and vertically in the gas phase.     

Enthalpy of the gas-phase CO2 + Mg reaction from ab initio total energies

Various highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3(MP2), G3//B3LYP, G3(MP2)//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction enthalpies of the ground-state reaction of CO2 with Mg. All model chemistries predict highly endothermic reactions, with DeltaH(298) = 63.6-69.7 kcal x mol(-1). The difference between the calculated reaction enthalpies and the experimental value, evaluated with recommended experimental standard enthalpies of formation for products and reactants, is more than 20 kcal x mol(-1) for all methods. This difference originates in the incorrect experimental enthalpy of formation of gaseous MgO given in thermochemical databases. When the theoretical formation enthalpy for MgO calculated by a particular method is used, the deviation is reduced to 1.3 kcal x mol(-1). The performance of the methodologies used to calculate the heat of this particular reaction and the enthalpy of formation of MgO are discussed.     

Thermochemistry of disputed soot formation intermediates C4H3 and C4H5

Accurate isomeric energy differences and standard enthalpies of formation for disputed intermediates in soot formation, C(4)H(3) and C(4)H(5), have been determined through systematic extrapolations of ab initio energies. Electron correlation has been included through second-order Z-averaged perturbation theory (ZAPT2), and spin-restricted, open-shell coupled-cluster methods through triple excitations [ROCCSD, ROCCSD(T), and ROCCSDT] utilizing the correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, and 6), followed by extrapolations to the complete basis set limit via the focal point method of Allen and co-workers. Reference geometries were fully optimized at the ROCCSD(T) level with a TZ(2d1f,2p1d) basis set. Our analysis finds that the resonance-stabilized i-C(4)H(3) and i-C(4)H(5) isomers lie 11.8 and 10.7 kcal mol(-1) below E-n-C(4)H(3) and E-n-C(4)H(5), respectively, several kcal mol(-1) (more, less) than reported in recent (diffusion Monte Carlo, B3LYP density-functional) studies. Moreover, in these systems Gaussian-3 (G3) theory suffers from large spin contamination in electronic wave functions, poor reference geometries, and anomalous vibrational frequencies, but fortuitous cancellation of these sizable errors leads to isomerization energies apparently accurate to 1 kcal mol(-1). Using focal-point extrapolations for isodesmic reactions, we determine the enthalpies of formation (delta(f)H(0) (composite function)) for i-C(4)H(3), Z-n-C(4)H(3), E-n-C(4)H(3), i-C(4)H(5), Z-n-C(4)H(5), and E-n-C(4)H(5) to be 119.0, 130.8, 130.8, 78.4, 89.7, and 89.1 kcal mol(-1), respectively. These definitive values remove any remaining uncertainty surrounding the thermochemistry of these isomers in combustion models, allowing for better assessment of whether even-carbon pathways contribute to soot formation.     

DFT Calculations of the Ionization Potential and Electron Affinity of Alaninamide

Adiabatic and vertical ionization potentials (IPs) and valence electron affinities (EAs) of alaninamide in gas phase have been determined using density functional theory (BLYP,B3LYP,B3P86) methods with 6-311++G(d,p) basis set,respectively. IPs and EAs of alaninamide in solutions have been calculated at the B3LYP/6-311++G(d,p) level. Five possible conformers of alaninamide and their charged states have been optimized employing density functional theory B3LYP method with 6-311++(d,p) basis set,respectively.     

Supramolecular binding thermodynamics by dispersion-corrected density functional theory

The equilibrium association free enthalpies ΔG(a) for typical supramolecular complexes in solution are calculated by ab initio quantum chemical methods. Ten neutral and three positively charged complexes with experimental ΔG(a) values in the range 0 to -21?kcal?mol(-1) (on average -6?kcal?mol(-1) ) are investigated. The theoretical approach employs a (nondynamic) single-structure model, but computes the various energy terms accurately without any special empirical adjustments. Dispersion corrected density functional theory (DFT-D3) with extended basis sets (triple-ζ and quadruple-ζ quality) is used to determine structures and gas-phase interaction energies (ΔE), the COSMO-RS continuum solvation model (based on DFT data) provides solvation free enthalpies and the remaining ro-vibrational enthalpic/entropic contributions are obtained from harmonic frequency calculations. Low-lying vibrational modes are treated by a free-rotor approximation. The accurate account of London dispersion interactions is mandatory with contributions in the range -5 to -60?kcal?mol(-1) (up to 200?% of ΔE). Inclusion of three-body dispersion effects improves the results considerably. A semilocal (TPSS) and a hybrid density functional (PW6B95) have been tested. Although the ΔG(a) values result as a sum of individually large terms with opposite sign (ΔE vs. solvation and entropy change), the approach provides unprecedented accuracy for ΔG(a) values with errors of only 2?kcal?mol(-1) on average. Relative affinities for different guests inside the same host are always obtained correctly. The procedure is suggested as a predictive tool in supramolecular chemistry and can be applied routinely to semirigid systems with 300-400 atoms. The various contributions to binding and enthalpy-entropy compensations are discussed.     

Thinking out of the black box: accurate barrier heights of 1,3-dipolar cycloadditions of ozone with acetylene and ethylene

Accurate barriers for the 1,3-dipolar cycloadditions of ozone with acetylene and ethylene have been determined via the systematic extrapolation of ab initio energies within the focal point approach of Allen and co-workers. Electron correlation has been accounted for primarily via coupled cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [CCSD, CCSD(T), CCSDT, and CCSDT(Q)]. For the concerted [4 + 2] cycloadditions, the final recommended barriers are DeltaH(0K) = 9.4 +/- 0.2 and 5.3 +/- 0.2 kcal mol(-1) for ozone adding to acetylene and ethylene, respectively. These agree with recent results of Cremer et al. and Anglada et al., respectively. The reaction energy for O3 + C2H2 exhibits a protracted convergence with respect to inclusion of electron correlation, with the CCSDT/cc-pVDZ and CCSDT(Q)/cc-pVDZ values differing by 2.3 kcal mol-1. Recommended enthalpies of formation (298 K) for cycloadducts 1,2,3-trioxole and 1,2,3-trioxolane are +32.8 and -1.6 kcal mol(-1), respectively. Popular composite ab initio approaches [CBS-QB3, CBS-APNO, G3, G3B3, G3(MP2)B3, G4, G4(MP3), and G4(MP2)] predict a range of barrier heights for these systems. The CBS-QB3 computed barrier for ozone and acetylene, DeltaH(0K) = 4.4 kcal mol(-1), deviates by 5 kcal mol(-1) from the focal point value. CBS-QB3 similarly underestimates the barrier for the reaction of ozone and ethylene, yielding a prediction of only 0.7 kcal mol(-1). The errors in the CBS-QB3 results are significantly larger than mean errors observed in application to the G2 test set. The problem is traced to the nontransferability of MP2 basis set effects in the case of these reaction barriers. The recently published G4 and G4(MP2) approaches perform substantially better for O3 + C2H2, predicting enthalpy barriers of 9.0 and 8.4 kcal mol(-1), respectively. For the prediction of these reaction barriers, the additive corrections applied in the majority of the composite approaches considered lead to worse agreement with the reference focal point values than would be obtained relying only on single point energies evaluated at the highest level of theory utilized within each composite method.     

Halogenated 6,13-bis(triisopropylsilylethynyl)-5,7,12,14-tetraazapentacene: applications for ambipolar air-stable organic field-effect transistors

Density functional theory calculations were performed to explore the influence of halogenation on the reorganization energies (λ), adiabatic ionization potentials (IPs), adiabatic electron affinities (EAs), and air stabilities of a series of pentacene (PENT) and tetraceno[2,3-b]thiophene (TbTH) derivatives. According to calculated IP and EA values, all well-known PENT and TbTH derivatives in this paper are air-stable p-channel but not air-stable n-channel organic field-effect transistors (OFETs) due to insufficient EAs, consistent with experimental observations. The calculated results show that attaching two or more halogen atoms onto air-unstable 6,13-bis(triisopropylsilylethynyl)-5,7,12,14-tetraazapentacene (TIPS-N4PENT) is sufficient for promoting ambipolar air-stable properties. The electronic coupling and band structure calculations indicate that halogenated TIPS-N4PENT derivatives have potential applications in high-performance ambipolar air-stable OFETs. They also provide rational guidelines for the design of ambipolar air-stable organic semiconductors (OSCs).     

Thermochemical properties of selenium fluorides, oxides, and oxofluorides

The bond dissociation energies (BDEs), fluoride and fluorocation affinities, and electron affinities of SeF(n) (n = 1-6), SeOF(n) (n = 0-4), and SeO(2)F(n) (n = 0-2) have been predicted with coupled cluster CCSD(T) theory extrapolated to the complete basis set limit. To achieve near chemical accuracy, additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies. These included corrections for core-valence effects, scalar relativistic effects, for first-order atomic spin-orbit effects, and vibrational zero point energies. The adiabatic BDEs contain contributions from product reorganization energies and, therefore, can be much smaller than the diabatic BDEs and can vary over a wide range. For thermochemical calculations, the adiabatic values must be used, whereas for bond strength and kinetic considerations, the diabatic values should be used when only small displacements of the atoms without change of the geometry of the molecule are involved. The adiabatic Se-F BDEs of SeF(n) (n = 1-6) are SeF(6) = 90, SeF(5) = 27, SeF(4) = 93, SeF(3) = 61, SeF(2) = 86, and SeF = 76 kcal/mol, and the corresponding diabatic values are SeF(6) = 90, SeF(5) = 88, SeF(4) = 93, SeF(3) = 74, SeF(2) = 86, and SeF = 76 kcal/mol. The adiabatic Se-O BDEs of SeO(n) (n = 1-3), SeOF(n) (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 23 to 107 kcal/mol, whereas the diabatic ones range from 62 to 154 kcal/mol. The adiabatic Se-F BDEs of SeOF(n) (n = 1-4) and SeO(2)F(n) (n = 1,2) range from 20 to 88 kcal/mol, whereas the diabatic ones range from 73 to 112 kcal/mol. The fluoride affinities of SeF(n), (n = 1-6), SeO(n), (n = 1-3), SeOF(n), (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 15 to 121 kcal/mol, demonstrating that the Lewis acidity of these species covers the spectrum from very weak (SeF(6)) to very strong (SeO(3)) acids. The electron affinities which are a measure of the oxidizing power of a species, span a wide range from 1.56 eV in SeF(4) to 5.16 eV in SeF(5) and for the free radicals are much higher than for the neutral molecules. Another interesting feature of these molecules and ions stems from the fact that many of them possess both a Se free valence electron pair and a free unpaired valence electron, raising the questions of their preferred location and their influence on the Se-F and Se═O bond strengths.     

Thermochemistry of key soot formation intermediates: C3H3 isomers

Accurate standard enthalpies of formation for allene, propyne, and four C3H3 isomers involved in soot formation mechanisms have been determined through systematic focal point extrapolations of ab initio energies. Auxiliary corrections have been applied for anharmonic zero-point vibrational energy, core electron correlation, the diagonal Born-Oppenheimer correction (DBOC), and scalar relativistic effects. Electron correlation has been accounted for via second-order Z-averaged perturbation theory (ZAPT2) and primarily through coupled-cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [ROCCSD, ROCCSD(T), ROCCSDT, and UCCSDT(Q)]. The correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, 6), was employed. The CCSDT(Q) corrections do not exceed 0.12 kcal mol(-)1 for the relative energies of the systems considered here, indicating a high degree of electron correlation convergence in the present results. Our recommended values for the enthalpies of formation are as follows: Delta(f)H(o)(0)(propargyl) = 84.76, Delta(f)H(o)(0) (1-propynyl) = 126.60, Delta(f)H(o)(0) (cycloprop-1-enyl) = 126.28, Delta(f)H(o)(0)(cycloprop-2-enyl) = 117.36, Delta(f)H(o)(0)(allene) = 47.41, and Delta(f)H(o)(0)(propyne) = 46.33 kcal mol(-1), with estimated errors no larger than 0.3 kcal mol(-1). The corresponding C3H3 isomerization energies are about 1 kcal mol(-1) larger than previous coupled-cluster results and several kcal mol(-1) below those previously obtained using density functional theory.