Enthalpy of formation and bond energies on unsaturated oxygenated hydrocarbons using G3MP2B3 calculation methods

Enthalpies of unsaturated oxygenated hydrocarbons and radicals corresponding to the loss of hydrogen atoms from the parent molecules are intermediates and decomposition products in the oxidation and combustion of aromatic and polyaromatic species. Enthalpies (Δ_fH⁰₂₉₈) are calculated for a set of 27 oxygenated and nonoxygenated, unsaturated hydrocarbons and 12 radicals at the G3MP2B3 level of theory and with the commonly used B3LYP/6‐311g(d,p) density functional theory (DFT) method. Standard enthalpies of formation (Δ_fH⁰₂₉₈) are determined from the calculated enthalpy of reaction (ΔH⁰_rxn,298) using isodesmic work reactions with reference species that have accurately known Δ_fH⁰₂₉₈ values. The deviation between G3MP2B3 and B3LYP methods is under ±0.5 kcal mol^?1 for 9 species, 18 other species differs by less than ±1 kcal mol^?1 , and 11 species differ by about 1.5 kcal mol^?1. Under them are 11 radicals derived from the above‐oxygenated hydrocarbons that show good agreement between G3MP2B3 and B3LYP methods. G3 calculations have been performed to further validate enthalpy values, where a discrepancy of more than 2.5 kcal mol^?1 exists between the G3MP3B3 and density functional results. Surprisingly the G3 calculations support the density functional calculations in these several nonagreement cases. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 633–648, 2005     

Calculation of Some Thermodynamic Properties and Detonation Parameters of 1‐Ethyl‐3‐Methyl‐H‐Imidazolium Perchlorate, [Emim][ClO4], on the Basis of CBS‐4M and CHEETAH Computations Supplemented by VBT Estimates

This paper estimates some thermochemical (in kcal mol^–1) and detonation parameters for the ionic liquid, [emim][ClO₄] and its associated solid in view of its investigation as an energetic material. The thermochemical values estimated, employing CBS‐4M computational methodology and volume‐based thermodynamics (VBT) include: lattice energy, U_POT([emim][ClO₄]) ≈? 123 ± 16 kcal · mol^–1; enthalpy of formation of the gaseous cation, Δ_fH°([emim]⁺, g) = 144.2 kcal · mol^–1 and anion, Δ_fH°([ClO₄]^–, g) = –66.1 kcal · mol^–1; the enthalpy of formation of the solid salt, Δ_fH°([emim][ClO₄],s) ≈? –55 ± 16 kcal · mol^–1 and for the associated ionic liquid, Δ_fH^o([emim][ClO₄],l) = –52 ± 16 kcal · mol^–1 as well as the corresponding Gibbs energy terms: Δ_fG°([emim][ClO₄],s) ≈? +29 ± 16 kcal · mol^–1 and Δ_fG^o([emim][ClO₄],l) = +24 ± 16 kcal · mol^–1 and the associated standard absolute entropies, of the solid [emim][ClO₄], S°₂₉₈([emim][ClO₄],s) = 83 ± 4 cal · K^–1 · mol^–1. The following combustion and detonation parameters are assigned to [emim][ClO₄] in its (ionic) liquid form: specific impulse (I_sp) = 228 s (monopropellant), detonation velocity (V_oD) = 5466 m · s^–1, detonation pressure (p_C–J) = 99 kbar, explosion temperature (T_ex) = 2842 K.     

High levelab initio stabilization energies of benzene

Summary G2 theory is shown to be reliable for calculating isodesmic and homodesmotic stabilization energies (ISE and HSE, respectively) of benzene. G2 calculations give HSE and ISE values of 92.5 and 269.1 kJ mol^–1 (298 K), respectively. These agree well with the experimental HSE and ISE values of 90.5±7.2 and 268.7±6.3 kJ mol^–1, respectively. We conclude that basis set superposition error corrections to the enthalpies of the homodesmotic or isodesmic reactions are not necessary in calculations of the stabilization energies of benzene using G2 theory. The calculated values of the enthalpies of formation of such molecules containing multiple bonds such as benzene ands-trans 1,3-butadiene, which are found from the enthalpies of isodesmic and homodesmotic reactions rather than of atomization reactions, demonstrate good performance of G2 theory. Estimates of theH _f ^o value for benzene from the G2 calculated enthalpies of homodesmotic reaction (2) and isodesmic reaction (3) are 80.9 and 82.5 kJ mol^–1 (298 K), respectively. These are very close to the experimentalH _f ^o value of 82.9±0.3 kJ mol^–1. TheH _f ^o value ofs-trans 1,3-butadiene calculated using the G2 enthalpy of isodesmic reaction (4) is 110.5 kJ mol^–1 and is in excellent agreement with the experimentalH _f ^o value of 110.0±1.1 kJ mol^–1.     

Heat of formation of Benzophenone Oxide

The heat of formation of benzophenone oxide, Ph₂CO₂, was measured using photoacoustic calorimetry. The enthalpy of the reaction Ph₂CN₂ + O₂ → Ph₂CO₂ + N₂ was found to be ?48.0 ±0.8 kcal mol^?1 and ΔH_f(Ph₂CN₂) was determined by measuring the reaction enthalpy for Ph₂CN₂ + EtOH → Ph₂CHOEt + N₂ (?53.6 ±1.0 kcal mol^?1). Taking ΔH_f(PhCHOEt) = ?10.6 kcal mol^?1 led to ΔH_f(Ph₂CN₂) = 99.2 ± 1.5 kcal mol^?1 and hence to ΔH_f(Ph₂CO₂) = 51.1 ± 2.0 kcal mol^?1. The results imply that the self-reaction of benzophenone oxide i.e., 2Ph₂CO₂ → 2Ph₂CO + O₂ is exothermic by ?76.0 ±4.0 kcal mol^?1.     

Bond dissociation energies and radical heats of formation in CH3Cl,CH2Cl2, CH3Br,CH2Br2, CH2FCl,and CHFCl2

An analysis of thermochemical and kinetic data on the bromination of the halomethanes CH_4–nX_n (X = F, Cl, Br; n = 1–3), the two chlorofluoromethanes, CH₂FCl and CHFCl₂, and CH₄, shows that the recently reported heats of formation of the radicals CH₂Cl, CHCl₂, CHBr₂, and CFCl₂, and the C? H bond dissociation energies in the matching halomethanes are not compatible with the activation energies for the corresponding reverse reactions. From the observed trends in CH₄ and the other halomethanes, the following revised ΔH°_f,298 (R) values have been derived: ΔH°_f(CH₂Cl) = 29.1 ± 1.0, ΔH°_f(CHCl₂) = 23.5 ± 1.2, ΔH_f(CH₂Br) = 40.4 ± 1.0, ΔH°_f(CHBr₂) = 45.0 ± 2.2, and ΔH°_f(CFCl₂) = ?21.3 ± 2.4 kcal mol^?1. The previously unavailable radical heat of formation, ΔH°_f(CHFCl) = ?14.5 ± 2.4 kcal mol^?1 has also been deduced. These values are used with the heats of formation of the parent compounds from the literature to evaluate C? H and C? X bond dissociation energies in CH₃Cl, CH₂Cl₂, CH₃Br, CH₂Br₂, CH₂FCl, and CHFCl₂.     

Evaluation of DFT Methods and Implicit Solvation Models for Anion‐Binding Host‐Guest Systems

Although supramolecular chemistry is traditionally an experimental discipline, computations have emerged as important tools for the understanding of supramolecules. We have explored how well commonly used density functional theory quantum mechanics and polarizable continuum solvation models can calculate binding affinities of host‐guest systems. We report the calculation of binding affinities for eight host–guest complexes and compare our results to experimentally measured binding free energies that span the range from ?2.3 to ?6.1 kcal mol^?1. These systems consist of four hosts (biotin[6]uril, triphenoxymethane, cryptand, and bis‐thiourea) with different halide ions (F^?, Cl^?, Br^?) in various media including organic and aqueous. The mean average deviation (MAD) of calculated from measured ΔG_a is 2.5 kcal mol^?1 when using B3LYP‐D3 with either CPCM or PCM. This MAD value lowers even more by eliminating two outliers: 1.1 kcal mol^?1 for CPCM and 1.2 kcal mol^?1 for PCM. The best DFT and implicit solvation model combination that we have studied is B3LYP?D3 with either CPCM or PCM.     

Data on the thermochemistry of azoalkanes

The rate constant of the primary decomposition step was determined for four symmetrical and four unsymmetrical azoalkanes. From the experimental activation energies and some literature enthalpy data, the following enthalpies of formation of radicals and group contributions were calculated: ΔH_? (CH₃N₂) = 51.5 ± 1.8 kcal mol^?1, ΔH_? (C₂H₅N₂) = 44.8 ± 2.5 kcal mol^?1, ΔH_? (2?C₃H₇N₂) = 37.9 ± 2.2 kcal mol^?1, [N_A-(C)] = 27.6 ± 3.7 kcal mol^?1, [N_A-(?_A) (C)] = 61.2 ± 3.1 kcal mol^?1.     

Thermochemical properties for n‐propyl,iso‐propyl,and tert‐butyl nitroalkanes,alkyl nitrites,and their carbon‐centered radicals

Density functional theory (DFT) based calculations are performed on a series of alkyl nitrites and nitroalkanes representing large‐scale primary, secondary, and tertiary nitro compounds and their radicals resulting from the loss of their skeletal hydrogen atoms. Geometries, vibration frequencies, and thermochemical properties [S°(T) and C°_p(T) (10 K ? T ? 5000 K)] are calculated at the B3LYP/6‐31G(d,p) DFT level. Δ_fH°₂₉₈ values are from B3LYP/6‐31G(d,p), B3LYP/6‐31+G(2d,2p), and the composite CBS‐QB3 levels. Potential energy barriers for the internal rotations have been computed at the B3LYP/6‐31G(d,p) level of theory, and the lower barrier contributions are incorporated into entropy and heat capacity data. The standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Recommended values derived from the most stable conformers of respective nitro‐ and nitrite isomers include ?30.57 and ?28.44 kcal mol^?1 for n‐propane‐, ?33.89 and ?32.32 kcal mol^?1 for iso‐propane‐, ?42.78 and ?41.36 kcal mol^?1 for tert‐butane‐nitro compounds and nitrites, respectively. Entropy and heat capacity values are also reported for the lower homologues: nitromethane, nitroethane, and corresponding nitrites. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 181–199, 2010     

The Heat of Formation of the Uranyl Dication: Theoretical Evaluation Based on Relativistic Density Functional Calculations

By using a set of model reactions, we estimated the heat of formation of gaseous UO₂²⁺ from quantum‐chemical reaction enthalpies and experimental heats of formation of reference species. For this purpose, we performed relativistic density functional calculations for the molecules UO₂²⁺, UO₂, UF₆, and UF₅. We used two gradient‐corrected exchange‐correlation functionals (revised Perdew–Burke–Ernzerhof (PBEN) and Becke–Perdew (BP)) and we accounted for spin‐orbit interaction in a self‐consistent fashion. Indeed, spin‐orbit interaction notably affects the energies of the model reactions, especially if compounds of U^IV are involved. Our resulting theoretical estimates for Δ_f (UO₂²⁺), 365±10 kcal mol^?1 (PBEN) and 370±12 kcal mol^?1 (BP), are in quantitative agreement with a recent experimental result, 364±15 kcal mol^?1. Agreement between the results of the two different exchange‐correlation functionals PBEN and BP supports the reliability of our approach. The procedure applied offers a general means to derive unknown enthalpies of formation of actinide species based on the available well‐established data for other compounds of the element in question.     

Enthalpy of formation of aqueous tungstate ion and of crystalline tungstic acid (H2WO4)

Calorimetric measurements of the enthalpy of reaction of WO₃(c) with excess OH^?(aq) have been made at 85°C. Similar measurements have been made with MoO₃(c) at both 85 and 25°C, to permit estimation of ΔH°=?13.4 kcal mol^?1 for the reaction WO₃(c)+2OH^?(aq)=WO^2?₄(aq)+H₂O(liq) at 25°C. Combination of this ΔH° with ΔH°_f for WO₃(c) leads to ΔH°_f=?256.5 kcal mol^?1 for WO^2?₄(aq). We also obtain ΔH°_f=?269.5 kcal mol^?1 for H₂WO₄(c). Both of these values are discussed in relation to several earlier investigations.     

Heterolysis of H2 Across a Classical Lewis Pair, 2,6‐Lutidine⋅BCl3: Synthesis,Characterization, and Mechanism

We report that 2,6‐lutidine?trichloroborane (Lut?BCl₃) reacts with H₂ in toluene, bromobenzene, dichloromethane, and Lut solvents producing the neutral hydride, Lut?BHCl₂. The mechanism was modeled with density functional theory, and energies of stationary states were calculated at the G3(MP2)B3 level of theory. Lut?BCl₃ was calculated to react with H₂ and form the ion pair, [LutH⁺][HBCl₃^?], with a barrier of ΔH^≠=24.7 kcal mol^?1 (ΔG^≠=29.8 kcal mol^?1). Metathesis with a second molecule of Lut?BCl₃ produced Lut?BHCl₂ and [LutH⁺][BCl₄^?]. The overall reaction is exothermic by 6.0 kcal mol^?1 (Δ_rG°=?1.1). Alternate pathways were explored involving the borenium cation (LutBCl₂⁺) and the four‐membered boracycle [(CH₂{NC₅H₃Me})BCl₂]. Barriers for addition of H₂ across the Lut/LutBCl₂⁺ pair and the boracycle B?C bond are substantially higher (ΔG^≠=42.1 and 49.4 kcal mol^?1, respectively), such that these pathways are excluded. The barrier for addition of H₂ to the boracycle B?N bond is comparable (ΔH^≠=28.5 and ΔG^≠=32 kcal mol^?1). Conversion of the intermediate 2‐(BHCl₂CH₂)‐6‐Me(C₅H₃NH) to Lut?BHCl₂ may occur by intermolecular steps involving proton/hydride transfers to Lut/BCl₃. Intramolecular protodeboronation, which could form Lut?BHCl₂ directly, is prohibited by a high barrier (ΔH^≠=52, ΔG^≠=51 kcal mol^?1).     

Kinetic and thermochemical parameters of chlorine atom transfer reactions from CF3Cl,CF3CF2Cl,CF2ClCF2Cl,and CF2ClCFCl2 in gas phase

The following reactions: (1) were studied over the temperature ranges 533–687 K, 563–663 K, and 503–613 K for the forward reactions respectively and over 683–763 K, for the back reaction. Arrhenius parameters for chlorine atom transfer were determined relative to the combination of the attacking radicals. The ΔH_r°(1) = ?3.95 ± 0.45 kcal mol^?1 was calculated and from this value the ΔH∮(C₂F₅Cl) = ?2.66.3 ± 2.5 kcal mol^?1 and D(C₂F₅-Cl) = 82.0 ± 1.2 kcal mol^?1 were obtained. Besides, the ΔH_r°(2) was estimated leading to D(CF₂ClCF₂Cl) = 79.2 ± 5 Kcal mol^?1. The bond dissociation energies and the heat of formation are compared with those of the literature. The effect of the halogen substitutents as well as the importance of the polar effects for halogen transfer processes are discussed.     

Chiral discrimination in hydrogen‐bonded complexes of 2‐methylol oxirane with hydrogen peroxide

A systematic quantum chemical study reveals the effects of chirality on the intermolecular interactions between two chiral molecules bound by hydrogen bonds. The methods used are second‐order Møller–Plesset perturbation theory (MP2) with the 6‐311++g(d,p) basis set. Complexes via the O? H···O hydrogen bond formed between the chiral 2‐methylol oxirane (S) and chiral HOOH (P and M) molecules have been investigated, which lead to four diastereomeric complexes. The nomenclature of the complexes used in this article is enantiomeric configuration sign corresponding to English letters. Such as: sm, sp. The relative positions of the methylol group and the hydrogen peroxide are designated as syn (same side) and anti (opposite side). The largest chirodiastaltic energy was ΔE_chir = ?1.329 kcal mol^?1 [9% of the counterpoise correct average binding energy D_e(corr)] between the sm‐syn and sp‐anti in favor of sm‐syn. The largest diastereofacial energy was ?1.428 kcal mol^?1 between sm‐syn and sm‐anti in favor of sm‐syn. To take into account solvents effect, the polarizable continuum model (PCM) method has been used to evaluate the chirodiastaltic energies, and diastereofacial energies of the 2‐methylol oxirane···HOOH complexes. The chiral 2,3‐dimethylol oxirane (S, S) is C₂ symmetry which offers two identical faces. Hence, the chirodiastaltic energy is identical to the diastereomeric energy, and is ΔE_chir = 0.563 kcal mol^?1 or 5.3% of the D_e(corr) in favor of s,s‐p. The optimized structures, interaction energies, and chirodiastaltic energies for various isomers were estimated. The harmonic frequencies, IR intensities, rotational constants, and dipole moments were also reported. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

London Dispersion Effects in a Distannene/Tristannane Equilibrium: Energies of their Interconversion and the Suppression of the Monomeric Stannylene Intermediate

Reaction of {LiC₆H₂−2,4,6-Cyp₃⋅Et₂O}₂ (Cyp=cyclopentyl) ( 1 ) of the new dispersion energy donor (DED) ligand, 2,4,6-triscyclopentylphenyl with SnCl₂ afforded a mixture of the distannene {Sn(C₆H₂−2,4,6-Cyp₃)₂}₂ ( 2 ), and the cyclotristannane {Sn(C₆H₂−2,4,6-Cyp₃)₂}₃ ( 3 ). 2 is favored in solution at higher temperature (345 K or above) whereas 3 is preferred near 298 K. Van't Hoff analysis revealed the 3 to 2 conversion has a ΔH=33.36 kcal mol⁻¹ and ΔS=0.102 kcal mol⁻¹ K⁻¹, which gives a ΔG^{300 K}=+2.86 kcal mol⁻¹, showing that the conversion of 3 to 2 is an endergonic process. Computational studies show that DED stabilization in 3 is −28.5 kcal mol⁻¹ per {Sn(C₆H₂−2,4,6-Cyp₃)₂ unit, which exceeds the DED energy in 2 of −16.3 kcal mol⁻¹ per unit. The data clearly show that dispersion interactions are the main arbiter of the 3 to 2 equilibrium. Both 2 and 3 possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).     

Outer‐Sphere 2 e−/2 H+ Transfer Reactions of Ruthenium(II)‐Amine and Ruthenium(IV)‐Amido Complexes

A diverse set of 2 e⁻/2 H⁺ reactions are described that interconvert [Ru^II(bpy)(en*)₂]²⁺ and [Ru^IV(bpy)(en‐H*)₂]²⁺ (bpy=2,2′‐bipyridine, en*=H₂NCMe₂CMe₂NH₂, en*‐H=H₂NCMe₂CMe₂NH⁻), forming or cleaving different O−H, N−H, S−H, and C−H bonds. The reactions involve quinones, hydrazines, thiols, and 1,3‐cyclohexadiene. These proton‐coupled electron transfer reactions occur without substrate binding to the ruthenium center, but instead with precursor complex formation by hydrogen bonding. The free energies of the reactions vary over more than 90 kcal mol⁻¹, but the rates are more dependent on the type of X−H bond involved than the associated ΔG °. There is a kinetic preference for substrates that have the transferring hydrogen atoms in close proximity, such as ortho ‐tetrachlorobenzoquinone over its para ‐isomer and 1,3‐cyclohexadiene over its 1,4‐isomer, perhaps hinting at the potential for concerted 2 e⁻/2 H⁺ transfers.     

Definitive heat of formation of methylenimine,CH2NH,and of methylenimmonium ion,CH2NH2+, by means of W2 theory

A long‐standing controversy concerning the heat of formation of methylenimine has been addressed by means of the W2 (Weizmann‐2) thermochemical approach. Our best calculated values, ΔH°_f,298(CH₂NH) = 21.1±0.5 kcal/mol and ΔH°_f,298(CH₂NH₂⁺) = 179.4±0.5 kcal/mol, are in good agreement with the most recent measurements but carry a much smaller uncertainty. As a byproduct, we obtain the first‐ever accurate anharmonic force field for methylenimine: upon consideration of the appropriate resonances, the experimental gas‐phase band origins are all reproduced to better than 10 cm^?1. Consideration of the difference between a fully anharmonic zero‐point vibrational energy and B3LYP/cc‐pVTZ harmonic frequencies scaled by 0.985 suggests that the calculation of anharmonic zero‐point vibrational energies can generally be dispensed with, even in benchmark work, for rigid molecules. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1297–1305, 2001     

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using complete basis set and Gaussian‐n models combined with CPCM continuum solvation methods

Complete basis set and Gaussian‐n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pK_a values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of ?264.61 kcal/mol for the free energy of solvation of H⁺, ΔG_s(H⁺), was combined with a value for G_gas(H⁺) of ?6.28 kcal/mol, to calculate pK_a values with cycle 1. The complete basis set gas‐phase methods used to calculate gas‐phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed‐phase free energies are slightly less accurate than the gas‐phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pK_a values accurate to less than half a pK_a unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Thermochemistry and Kinetics of the Thermal Decomposition of 1‐Chlorohexane

A theoretical kinetic study of the thermal decomposition of 1‐chlorohexane in gas phase between 600 and 1000 K was performed. Transition‐state theory and unimolecular reaction rate theory were combined with molecular information provided by quantum chemical calculations. Particularly, the B3LYP, BMK, M05–2X, and M06–2X formulations of the density functional theory (DFT) and the high‐level ab initio methods G3B3 and G4 were employed. The possible reaction channels for the thermal decomposition of 1‐chlorohexane were investigated, and the reaction takes place through the elimination of HCl with the formation of 1‐hexene. The derived high‐pressure limit rate coefficients are k _∞(600–1000 K) = (8 ± 5) × 10¹³ exp[‐((56.7 ± 0.4) kcal mol⁻¹/RT )] s⁻¹. The pressure effect over the reaction was analyzed from the calculation of the low‐pressure limit rate coefficients and the falloff curves. In addition, the standard enthalpies of formation at 298 K of −46.9 ± 1.5 kcal mol⁻¹ for 1‐chlorohexane and 5.8 ± 1.5 kcal mol⁻¹ for C₆H₁₃ radical were derived from isodesmic and isogiric reactions at high levels of theory.     

Photochemical Generation of Strained Cycloalkynes from Methylenecyclopropanes

The hydrocarbons 1‐cyclopentylidene‐1a,9b‐dihydro‐1H ‐cyclopropa[l ]phenanthrene and 1‐cyclobutylidene‐1a,9b‐dihydro‐1H ‐cyclopropa[l ]phenanthrene undergo photolysis in solution at ambient temperature to produce cyclohexyne and cyclopentyne, respectively. These strained cycloalkynes, formed via the putative cycloalkylidenecarbenes, were intercepted as Diels–Alder adducts. Calculations at the CCSD(T)/cc‐pVTZ//B3LYP/6‐31+G* level of theory show that singlet cyclopentylidenecarbene has to overcome a barrier of 9.1 kcal mol⁻¹ to rearrange into cyclohexyne (with ΔE for ring expansion=−15.1 kcal mol⁻¹). By contrast, cyclobutylidenecarbene only needs to surmount a barrier of 1.6 kcal mol⁻¹ to rearrange into cyclopentyne (with ΔE for ring expansion=−6.2 kcal mol⁻¹).