Ab initio and density functional theory study of the enthalpies of formation of F2SOx and FClSOx (x=1, 2)

Enthalpies of formation of F₂SO, F₂SO₂, FClSO and FClSO₂ molecules have been determined using ab initio molecular orbital theory and density functional theory (DFT) calculations. Different DFT approaches and levels of the Gaussian-3 and the complete basis set (CBS) ab initio model chemistries have been employed to calculate enthalpies of formation from both total atomization energies and isodesmic reaction schemes. The best values at 298 K for F₂SO, F₂SO₂, FClSO and FClSO₂ as derived from an average of G3, G3B3, CBS-Q and CBS-QB3 isodesmic energies are −140.6, −181.1, −92.6 and −132.3 kcal mol⁻¹, respectively. The results obtained suggest that the accumulated small component errors found in the DFT-based methods are significantly reduced at the ab initio levels employed. Structural properties, harmonic vibrational frequencies, mode assignations and infrared intensities derived from B3LYP and mPW1PW91 functional with the 6-311+G(3df) basis set are presented.     

Thermochemical database of Halomethanes, Halosilanes, Halophosphines, and Haloamines

An accurate thermochemical database for 28 halides of carbon, silicon, nitrogen, and phosphorus is presented. The database provides improved standard enthalpies of formation for several compounds of ecological importance (CH3F, CF2Cl2, CFCl3) together with enthalpies of other compounds which are not known due to experimental difficulties in measuring their enthalpies. We also present a comparison of the latest ab initio methods (CBS-QB3 and G3) which are used for thermochemical predictions. The comparison shows that the G3 method consistently underestimates delta H degree f by 1-2 kJ/mol (relative to CBS-QB3).     

Thermochemistry of 2- and 3-acetylthiophenes: calorimetric and computational study

The relative stabilities of 2- and 3-acetylthiophenes have been evaluated by experimental thermochemistry and the results compared to high-level ab initio calculations. The enthalpies of combustion, vaporization, and sublimation were measured by rotating-bomb combustion calorimetry, Calvet microcalorimetry, correlation gas chromatography, and Knudsen effusion techniques and the gas-phase enthalpies of formation, at T = 298.15 K, were determined. Standard ab initio molecular orbital calculations at the G2 and G3 levels were performed, and a theoretical study on the molecular and electronic structures of the compounds studied has been conducted. Calculated enthalpies of formation using atomization and isodesmic reactions are compared with the experimental data. Experimental and theoretical results show that 2-acetylthiophene is thermodynamically more stable than the 3-isomer. A comparison of the substituent effect of the acetyl group in benzene and thiophene rings has been carried out.     

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.     

Ab initio vs molecular mechanics thermochemistry: homocubanes

The standard enthalpies of formation and strain energies for a series of homocubanes have been investigated by high-level ab initio G3(MP2)/B3LYP method. The relative stabilities of isomers are discussed. The comparison is made between the results of ab initio and molecular mechanics methods with the aim of assessing their performances. The usefulness of high-level calculations for generating thermochemical databases of relatively large molecules (e.g. C(11)H(14)) was also demonstrated.     

Ab initio thermochemistry of some halogenated cyclopropanes

The standard enthalpies of formation for a series of chloro- and fluoro-substituted cyclopropanes have been calculated by using high-level ab initio G3/B3LYP methods. The relative stabilities of isomers and the influence of substituents on thermochemistry in several classes of substituted derivatives are discussed.     

Calorimetric and computational study of the thermochemistry of phenoxyphenols

Thermodynamic properties of 3- and 4-phenoxyphenol have been determined by using a combination of calorimetric and effusion techniques as well as by high-level ab initio molecular orbital calculations. The standard (p° = 0.1 MPa) molar enthalpies of formation in the condensed and gas states, Δ(f)H(m)°(cr or l) and Δ(f)H(m)°(g), at T = 298.15 K, of 3- and 4-phenoxyphenol were derived from their energies of combustion in oxygen, measured by a static bomb calorimeter, and from the enthalpies of vaporization or sublimation derived respectively by Calvet microcalorimetry for the 3-phenoxyphenol and by Knudsen effusion technique for the 4-phenoxyphenol. The theoretically estimated gas-phase enthalpies of formation were calculated from high-level ab initio molecular orbital calculations at the G3(MP2)//B3LYP level of theory. Furthermore, this composite approach was also used to obtain information about the gas-phase acidities, gas-phase basicities, proton and electron affinities, adiabatic ionization enthalpies, and, finally, O?H bond dissociation enthalpies. The good agreement between the G3MP2B3-derived values and the experimental gas-phase enthalpies of formation for the 3- and 4-phenoxyphenol gives confidence to the estimate concerning the 2-phenoxyphenol isomer, which was not experimentally studied, and to the estimates concerning the radical and the anion. Additionally, the experimental values of gas-phase enthalpies of formation were also compared with estimates based on the empirical scheme developed by Cox.     

Structures and energies of D-galactose and galabiose conformers as calculated by ab initio and semiempirical methods

Optimized geometries and total energies of some conformers of alpha- and beta-D-galactose have been calculated using the RHF/6-31G* ab initio method. Vibrational frequencies were computed at the 6-31G* level for the conformers that favor internal hydrogen bonding, in order to evaluate their enthalpies, entropies, Gibbs free energies, and then their structural stabilities. The semiempirical AM1, PM3, MNDO methods have also been performed on the conformers GG, GT, and TG of alpha- and beta-D-galactose. In order to test the reliability of each semiempirical method, the obtained structures and energies from the AM1, PM3, and MNDO methods have been compared to those achieved using the RHF/6-31G* ab initio method. The MNDO method has not been investigated further, because of the large deviation in the structural parameters compared with those obtained by the ab initio method for the galactose. The semiempirical method that has yielded the best results is AM1, and it has been chosen to perform structural and energy calculations on the galabiose molecule (the disaccharides constituted by two galactose units alpha 1,4 linked). The goal of such calculations is to draw the energy surface maps for this disaccharide. To realize each map, 144 different possible conformations resulting from the rotations of the two torsional angles psi and phi of the glycosidic linkage are considered. In each calculation, at each increment of psi and phi, using a step of 30 degrees from 0 to 330 degrees, the energy optimization is employed. In this article, we report also calculations concerning the galabiose molecule using different ab initio levels such as RHF/6-31G*, RHF/6-31G**, and B3Lyp/6-31G*.     

Pairwise-substitution effects and intramolecular hydrogen bonds in nitrophenols and methylnitrophenols. Thermochemical measurements and ab initio calculations

The standard molar enthalpies of formation in the gaseous state of a series of nitrophenols, 2-nitrophenol, 3-nitrophenol, 4-nitrophenol, 5-methyl-2-nitrophenol, 2-methyl-5-nitrophenol, and 3-methyl-4-nitrophenol, have been obtained from combustion calorimetry and results from the temperature dependence of the vapor pressure measured by the transpiration method. To verify the experimental data, ab initio calculations of all compounds have been performed using MP, DFT, and G3 methods. Enthalpies of formation derived from the G3 methods are in a good agreement with the experimental results. The quantitative analysis of ortho, meta, and para pairwise-substituent effects in nitrophenols has been performed, and the strength of intramolecular hydrogen bonding in o-nitrophenol has been derived from thermochemical results and compared with those obtained from spectroscopic experiments and ab initio calculations. The new results help to resolve uncertainties in the available thermochemical data on extended series of nitrophenols.     

Calorimetric and computational study of 1,3,5-trithiane

To understand the differences in conformational behavior and reactivity of oxygen- and sulfur-containing 1,3,5-heterocyclohexanes, the enthalpies of formation and sublimation of 1,3,5-trithiane, 1, have been measured. The numerical value of the enthalpy of formation for this compound in the solid state is -8.6 +/- 2.6 kJ mol(-1), while the corresponding value in the gaseous state is 84.6 +/- 2.6 kJ mol(-1). The value for the enthalpy of sublimation is 93.2 +/- 0.2 kJ mol(-1). Standard ab initio molecular orbital calculations at the G2(MP2), G2, and G3 levels were performed, and the calculated enthalpies of formation are compared with the experimental data. These experimental and theoretical studies support the relevance of through-space lone pair-lone pair electronic repulsion in the sulfur heterocycle.     

Ionic liquids. Combination of combustion calorimetry with high-level quantum chemical calculations for deriving vaporization enthalpies

In this work, the molar enthalpies of formation of the ionic liquids [C2MIM][NO3] and [C4MIM][NO3] were measured by means of combustion calorimetry. The molar enthalpy of fusion of [C2MIM][NO3] was measured using differential scanning calorimetry. Ab initio calculations of the enthalpy of formation in the gaseous phase have been performed for the ionic species using the G3MP2 theory. We have used a combination of traditional combustion calorimetry with modern high-level ab initio calculations in order to obtain the molar enthalpies of vaporization of a series of the ionic liquids under study.     

DFT and ab initio calculations on two reactions between hydrogen atoms and the fire suppressants 2-H heptafluoropropane and CF3Br

Reaction enthalpies and barrier heights of the reactions CF3Br+H-->CF3+HBr {reaction (1)} and CF3CHFCF3+H-->CF3CFCF3+H2 {reaction (2)} have been calculated at the near state-of-the-art ab initio level, and also by employing the B3LYP, BH&HLYP, BB1K, MPW1K, MPWB1K and TPSS1KCIS functionals. In addition, the integrated molecular orbital+molecular orbital (IMOMO) method has been used to study reaction (2). The ab initio benchmark values of the reaction enthalpy (298 K) and barrier height (0 K) of reaction (2) are reported for the first time {-(0.7+/-0.7) and 13.3+/-0.5 kcal/mole respectively}. When density functional theory (DFT) results are compared with ab initio benchmarks for both reactions (1) and (2), the MPWB1K functional is found to have the best performance of the six functionals used. The IMOMO method with the RCCSD/aug-cc-pVTZ and/or RCCSD(T)/aug-cc-pVTZ levels, as the high levels of calculation on the model system, gives reaction enthalpies and barrier heights of reaction (2), which agree with ab initio benchmark values to within 1 kcal/mole. Computed key geometrical parameters and imaginary vibrational frequencies of the transition state structures of reactions (1) and (2) obtained at different levels of calculation are compared. The magnitudes of the computed imaginary vibrational frequencies of the transition states of both reactions considered are found to be very sensitive to the levels of calculation used to obtain them. The heat of formation (298 K) of CF3CFCF3 calculated at the near state-of-the-art level has a value of -(318+/-3) kcal/mole.     

The fragment molecular orbital method for geometry optimizations of polypeptides and proteins

The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.     

Thermochemistry of ionic liquid-catalyzed reactions. experimental and theoretical study of chemical equilibria of isomerization and transalkylation of tert-butylbenzenes

The chemical equilibrium of mutual interconversions of tert-butylbenzenes was studied in the temperature range 286 to 423 K using chloroaluminate ionic liquids as a catalyst. Enthalpies of five reactions of isomerization and transalkylation of tert-butylbenzenes were obtained from temperature dependences of the corresponding equilibrium constants in the liquid phase. Molar enthalpies of vaporization of methyl-tert-butylbenzenes and 1,4-di tert-butylbenzene were obtained by the transpiration method and were used for a recalculation of enthalpies of reactions and equilibrium constants into the gaseous phase. Using these experimental results, ab initio methods (B3LYP and G3MP2) have been tested for prediction thermodynamic functions of the five reactions under study successfully. Thermochemical investigations of tert-butylbenzenes available in the literature combined with experimental results have helped to resolve contradictions in the available thermochemical data for tert-butylbenzene and to recommend consistent and reliable enthalpies of formation for this compound in the liquid and the gaseous state.     

One- and three-parameter density functional theory and high level ab initio theory study of the CH+CH disproportionation reaction

Complete basis set (CBS) ab initio computational studies were performed with the target being to explore the CH+CH potential energy surface. Several closed and open shell intermediates were located on the potential energy surface. Computed enthalpies for the branching reactions, as well as heats of formation are in excellent agreement. Although CBS computed energies are of high quality, this computational study is not capable of predicting the branching product ratio due to fact that neither the MP2 nor the 6-311G(2d,2p) basis set are sufficient to locate the reactant complexes and the transition state structures for the hydrogen and carbon transfer reactions in the reaction complexes. To properly explore the CH+CH potential energy surface a much higher ab initio theory level is required.     

反应类等键反应方法及类反应热动力学参数计算

本文提出了反应类等键反应方法,将通常用于热力学性质计算的等键反应方法推广用于类反应中反应势垒和反应焓变的计算. 对碳氢燃料低温燃烧反应机理中的一类重要反应类：?烷基自由基过氧化氢裂解生成烯烃和HO2自由基的反应势垒和反应焓变进行了计算. 通过对该类16个反应中的5个代表反应分别在不同计算水平HF、DFT、MP2、CCSD(T)的比较计算发现,采用等键反应方法可在较低从头算级别计算得到类反应较高精度的反应势垒,提高了计算的效率和精度. 本文采用反应类等键反应方法在B3LYP/6-311G(d,p)计算水平对该类16个反应进行了反应势垒和反应焓变的计算,并建立了反应势垒和反应能的线性自由能关系：delta V=71.02+0.41?delta E.     

Theoretical thermochemistry: Enthalpies of formation of a set of nitrogen‐containing compounds

Gas‐phase standard state formation enthalpies of 63 nitrogen‐containing compounds (NCCs) were studied by computational methods. Gaussian‐n and complete basis set composite methods were applied. After the calculation of the set of NCCs, the results were analyzed in several ways. All the seven selected methods depicted various precisions in this work. According to the calculations and further data processing, G4 was proven to be appropriate and capable of formation enthalpy calculations on NCCs with a mean absolute deviation of 0.63 kcal/mol. Thus, G4 calculation may help us on choosing experimental values and conducting predictions. Isodesmic reactions using G4 methods are also conducted and further correction was given. The high accuracy of G4 method makes it reliable for calculating gas‐phase enthalpies of formations and allows them to serve as a valuable check on the accuracy of reported experimental data. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

铀与氢初始反应产物形成热力学的理论研究

用量子化学从头计算的方法对铀与氢同位素分子或原子相互作用的某些初始反应产物的稳定几何构型进行全优化的基础上,研究铀与氢同位素相互作用的热力学,进而探讨了铀与氢初始反应产物的可能形成途径。     

Thermal features and thermochemistry of hexachlorohafnates of nitrogen aromatic bases

Hexachlorohafnates of pyridine and its three methyl-substituted derivatives were synthesized and examined by the thermoanalytical methods. The van't Hoff equation employed for the thermogravimetric αvs. T dependencies enabled evaluation of the heats of the thermal dissociation and subsequently enthalpies of formation and crystal lattice energies of the salts. Geometry and energy of formation of HfCl ₆ ^2? was determined at the ab initio Hartree-Fock SCF level, using all electron MINI basis set augmented with standard polarization functions (MINI^*). Electron correlation was considered at the MP2 level. Thermodynamic characteristics for the latter species were also obtained combining ab initio results with those of statistical thermodynamics. The usefulness of theoretical methods in examination of solid state energetics is briefly discussed.     

Experimental and computational thermochemical study of 2- and 3-thiopheneacetic acid methyl esters

Thiophene-based compounds have widespread use in modern drug design, biodiagnostics, electronic and optoelectronic devices, and conductive polymers. The present study reports an experimental and computational thermochemical study on the relative stabilities of 2- and 3-thiopheneacetic acid methyl esters. The enthalpies of combustion and vaporization were measured by a rotating-bomb combustion calorimeter, Calvet microcalorimetry, and correlation gas chromatography, and the gas-phase enthalpies of formation at T=298.15 K were determined. Standard ab initio molecular orbital calculations at the G3 level were performed, and a theoretical study of the molecular and electronic structure of the compounds studied was carried out. Calculated enthalpies of formation, using atomization and isodesmic reactions are in very good agreement with the experimental results.