Benson group additivity values of phosphines and phosphine oxides: Fast and accurate computational thermochemistry of organophosphorus species

Composite quantum chemical methods W1X-1 and CBS-QB3 are used to calculate the gas phase standard enthalpy of formation, entropy, and heat capacity of 38 phosphines and phosphine oxides for which reliable experimental thermochemical information is limited or simply nonexistent. For alkyl phosphines and phosphine oxides, the W1X-1, and CBS-QB3 results are mutually consistent and in excellent agreement with available G3X values and empirical data. In the case of aryl-substituted species, different computational methods show more variation, with G3X enthalpies being furthest from experimental values. The calculated thermochemical data are subsequently used to determine Benson group additivity contributions for 24 Benson groups and group pairs involving phosphorus, thereby allowing fast and accurate estimations of thermochemical data of many organophosphorus compounds of any complexity. Such data are indispensable, for example, in chemical process design or estimating potential hazards of new chemical compounds. © 2018 Wiley Periodicals, Inc.     

A theoretical study of the thermodynamics and kinetics of small organosulfur compounds

The performance of a large set of ab initio procedures in predicting geometries, thermochemical and kinetic data of small sulfur compounds is assessed. Geometries and thermochemical data for H₂S, (CH₃)₂S, H₂S₂, (CH₃)₂S₂ and H₂C=S are studied using the HF method, density functional theory methods (B3LYP, BHandHLYP, MPW1PW91 and BMK), post-HF methods [MP2, MP3, MP4, CCSD, CCSD(T) and QCISD] and composite techniques (G3, G3B3, CBS-QB3 and W1U). A set of five reactions involving these small organosulfur compounds is studied and the influence of the level of theory on transition state geometries, reaction barriers and rate coefficients is assessed. Independent of the level of theory used, accurate geometries are obtained with the 6-311G(2d,d,p) and cc-pVTZ basis sets, both reproducing experimental bond lengths and bond angles within 2 pm and 0.5°. Besides composite methods, the BMK/cc-pVTZ method is the only studied method that succeeds to predict standard enthalpies of formation within 10 kJ mol^?1 of the experimental data. The best agreement with experimental rate coefficients is obtained with the BHandHLYP/cc-pVTZ method, closely followed by the composite methods and the BMK/cc-pVTZ method. All these methods succeed to reproduce the experimental rate coefficients within a factor 4. To obtain an accurate prediction of both thermochemical and kinetic data for organosulfur compounds, the commonly used composite methods G3B3 and CBS-QB3 and the BMK/cc-pVTZ method prove to be valuable tools.     

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters, but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.     

Computational methods in organic thermochemistry. 1. Hydrocarbon enthalpies and free energies of formation

Standard state enthalpies and free energies of formation can be computed with reasonable accuracy (usually within 4 and often 2 kJ/mol) using high level model chemistries. A comparison set of nearly 300 organic compounds ranging from 1 to 10 carbon atoms having a variety of functional groups for which enthalpy and free energy literature values are available has been examined using G2, G2MP2, G3, G3MP2, G3B3, G3MP2B3, CBS-QB3, and density functional (B3LYP/6-311+G(3df,2p)) model chemistries. G3 gives an average mean absolute deviation of 3.0 and 13.4 kJ/mol for the enthalpies and free energies, respectively, using the atomization method and 3.1 and 3.7 kJ/mol when bond separation reactions are employed. G3 and G3B3 are the most accurate overall; the related G3MP2 and G3MP2B3 are nearly as accurate and can compute larger molecules. CBS-QB3 was also found to be accurate but is more limited in the size of molecules that can be computed. The density functional energies were found to have large deviations from the literature values using either the atomization or the bond separation method. Regardless of the model employed, the free energies are increasingly underestimated by computation as the size of the molecule increases. A series of corrections applied to the aliphatic hydrocarbons is presented, which usually reduces the deviations to less than 4 kJ/mol regardless of the size of the molecule.     

Thermochemistry of methyl and ethyl nitro, RNO2, and nitrite, RONO, organic compounds

Computational quantum theory is employed to determine the thermochemical properties of n-alkyl nitro and nitrite compounds: methyl and ethyl nitrites, CH3ONO and C2H5ONO, plus nitromethane and nitroethane, CH3NO2 and C2H5NO2, at 298.15 K using multilevel G3, CBS-QB3, and CBS-APNO composite methods employing both atomization and isodesmic reaction analysis. Structures and enthalpies of the corresponding aci-tautomers are also determined. The enthalpies of formation for the most stable conformers of methyl and ethyl nitrites at 298 K are determined to be -15.64 +/- 0.10 kcal mol-1 (-65.44 +/- 0.42 kJ mol-1) and -23.58 +/- 0.12 kcal mol-1 (-98.32 +/- 0.58 kJ mol-1), respectively. DeltafHo(298 K) of nitroalkanes are correspondingly evaluated at -17.67 +/- 0.27 kcal mol-1 (-74.1 +/- 1.12 kJ mol-1) and -25.06 +/- 0.07 kcal mol-1 (-121.2 +/- 0.29 kJ mol-1) for CH3NO2 and C2H5NO2. Enthalpies of formation for the aci-tautomers are calculated as -3.45 +/- 0.44 kcal mol-1 (-14.43 +/- 0.11 kJ mol-1) for aci-nitromethane and -14.25 +/- 0.44 kcal mol-1 (-59.95 +/- 1.84 kJ mol-1) for the aci-nitroethane isomers, respectively. Data are evaluated against experimental and computational values in the literature with recommendations. A set of thermal correction parameters to atomic (H, C, N, O) enthalpies at 0 K is developed, to enable a direct calculation of species enthalpy of formation at 298.15 K, using atomization reaction and computation outputs.     

A novel predictive model for formation enthalpies of Si and Ge hydrides with propane- and butane-like structures

Butane- and propane-like silicon-germanium hydrides and chlorinated derivatives represent a new class of precursors for the fabrication of novel metastable materials at low-temperature regimes compatible with selective growth and commensurate with the emerging demand for the reduced thermal budgets of complementary metal oxide semiconductor integration. However, predictive simulation studies of the growth process and reaction mechanisms of these new compounds, needed to accelerate their deployment and fine-tune the unprecedented low-temperature and low-pressure synthesis protocols, require experimental thermodynamic data, which are currently unavailable. Furthermore, traditional quantum chemistry approaches lack the accuracy needed to treat large molecules containing third-row elements such as Ge. Accordingly, here we develop a method to accurately predict the formation enthalpy of these compounds using atom-wise corrections for Si, Ge, Cl, and H. For a test set of 15 well-known hydrides of Si and Ge and their chlorides, such as Si(3)H(8), Ge(2)H(6), SiGeH(6), SiHCl(3), and GeCl(4), our approach reduces the deviations between the experimental and predicted formation enthalpies obtained from complete basis set (CBS-QB3), G2, and B3LPY thermochemistry to levels of 1-3 kcal/mol, or a factor of ～5 over the corresponding uncorrected values. We show that our approach yields results comparable or better than those obtained using homodesmic reactions while circumventing the need for thermochemical data of the associated reaction species. Optimized atom-wise corrections are then used to generate accurate enthalpies of formation for 39 pure Si-Ge hydrides and a selected group of 20 chlorinated analogs, of which some have recently been synthesized for the first time. Our corrected enthalpies perfectly reproduce the experimental stability trends of heavy butane-like compounds containing Ge. This is in contrast to the direct application of the CBS-QB3 method, which yields erratic predictions. Our approach also provides quantitative bond-additivity rules for the chlorination of these heavier species. Finally, we discuss structure and bonding trends across the entire sequence of butane-, propane-, and ethane-like molecules with a special focus on the isomeric variations.     

Thermochemical and kinetic analyses on oxidation of isobutenyl radical and 2-hydroperoxymethyl-2-propenyl radical

In recognition of the importance of the isobutene oxidation reaction in the preignition chemistry associated with engine knock, the thermochemistry, chemical reaction pathways, and reaction kinetics of the isobutenyl radical oxidation at low to intermediate temperature range were computationally studied, focusing on both the first and the second O2 addition to the isobutenyl radical. The geometries of reactants, important intermediates, transition states, and products in the isobutenyl radical oxidation system were optimized at the B3LYP/6-311G(d,p) and MP2(full)/6-31G(d) levels, and the thermochemical properties were determined on the basis of ab initio, density functional theory, and statistical mechanics. Enthalpies of formation for several important intermediates were calculated using isodesmic reactions at the DFT and the CBS-QB3 levels. The kinetic analysis of the first O2 addition to the isobutenyl radical was performed using enthalpies at the CBS-QB3 and G3(MP2) levels. The reaction forms a chemically activated isobutenyl peroxy adduct which can be stabilized, dissociate back to reactants, cyclize to cyclic peroxide-alkyl radicals, and isomerize to the 2-hydroperoxymethyl-2-propenyl radical that further undergoes another O2 addition. The reaction channels for isomerization and cyclization and further dissociation on this second O2 addition were analyzed using enthalpies at the DFT level with energy corrections based on similar reaction channels for the first O2 addition. The high-pressure limit rate constants for each reaction channel were determined as functions of temperature by the canonical transition state theory for further kinetic model development.     

Thermochemistry and bond dissociation energies of ketones

Ketones are a major class of organic chemicals and solvents, which contribute to hydrocarbon sources in the atmosphere, and are important intermediates in the oxidation and combustion of hydrocarbons and biofuels. Their stability, thermochemical properties, and chemical kinetics are important to understanding their reaction paths and their role as intermediates in combustion processes and in atmospheric chemistry. In this study, enthalpies (ΔH°(f 298)), entropies (S°(T)), heat capacities (C(p)°(T)), and internal rotor potentials are reported for 2-butanone, 3-pentanone, 2-pentanone, 3-methyl-2-butanone, and 2-methyl-3-pentanone, and their radicals corresponding to loss of hydrogen atoms. A detailed evaluation of the carbon-hydrogen bond dissociation energies (C-H BDEs) is also performed for the parent ketones for the first time. Standard enthalpies of formation and bond energies are calculated at the B3LYP/6-31G(d,p), B3LYP/6-311G(2d,2p), CBS-QB3, and G3MP2B3 levels of theory using isodesmic reactions to minimize calculation errors. Structures, moments of inertia, vibrational frequencies, and internal rotor potentials are calculated at the B3LYP/6-31G(d,p) density functional level and are used to determine the entropies and heat capacities. The recommended ideal gas-phase ΔH°(f 298), from the average of the CBS-QB3 and G3MP2B3 levels of theory, as well as the calculated values for entropy and heat capacity are shown to compare well with the available experimental data for the parent ketones. Bond energies for primary, secondary, and tertiary radicals are determined; here, we find the C-H BDEs on carbons in the α position to the ketone group decrease significantly with increasing substitution on these α carbons. Group additivity and hydrogen-bond increment values for these ketone radicals are also determined.     

Thermochemical properties and bond dissociation energies of C3-C5 cycloalkyl hydroperoxides and peroxy radicals: cycloalkyl radical + (3)O2 reaction thermochemistry

Cyclic aliphatic hydrocarbons are major components in modern fuels; they can be present in the reactants, and they can be formed during the gas-phase oxidation processes. In combustion and thermal oxidation processes, these cyclics will form radicals that react with (3)O(2) to form peroxy radicals. In this study, density functional theory and higher level ab initio calculations are used to calculate thermochemical properties and bond dissociation energies of 3-5-membered cycloalkanes, corresponding hydroperoxides, hydroperoxycycloalkyl radicals, and cycloalkyl radicals that occur in these reaction systems. Geometries, vibration frequencies, and thermochemical properties, ΔH(f 298)°, are calculated with the B3LYP/6-31 g(d,p), B3LYP/6-31 g(2d,2p), composite CBS-QB3, and G3MP2B3 methods. Standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Group additivity contributions are developed, and application of group additivity with comparison to calculated values is illustrated. Entropy and heat capacities, S°(T) and C(p)°(T) (5 K ≤ T ≤ 5000), are determined using geometric parameters and frequencies from the B3LYP/6-31 g(d,p) calculations.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

Calculation of the enthalpies of formation and proton affinities of some isoquinoline derivatives

Ab initio molecular orbital theory has been used to calculate enthalpies of formation of isoquinoline, 1-hydroxyisoquinoline, 5-hydroxyisoquinoline, and 1,5-dihydroxyisoquinoline as well as some pyridine and quinoline derivatives. The proton affinities of the four isoquinoline derivatives were also obtained. The high-level composite methods G3(MP2), G3(MP2)//B3LYP, G3//B3LYP, and CBS-QB3 have been used for this study, and the results have been compared with available experimental values. For six of the eight studied compounds, the theoretical enthalpies of formation were very close to the experimental values (to within 4.3 kJ · mol⁻¹); where comparison was possible, the theoretical and experimental proton affinities were also in excellent agreement with one another. However, there is an extraordinary discrepancy between theory and experiment for the enthalpies of formation of 1-hydroxyisoquinoline and 1,5-dihydroxyisoquinoline, suggesting that the experimental values for these two compounds should perhaps be re-examined. We also show that popular low cost computational methods such as B3LYP and MP2 show very large deviations from the benchmark values.     

Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

Quantum-chemical and group-additive calculations of the enthalpies of formation of hydroxylamines and oximes

The quality of calculations of the enthalpies of formation by the G2, G3, G3B3, and G3MP2B3 methods was estimated for substituted hydroxylamines and oximes. External correction of composite methods was performed to improve the accuracy of the results. The method of homodesmic reactions was used in calculations of large compounds comprising more than seven or eight nonhydrogen atoms. The enthalpies of formation of oximes and hydroxylamines were analyzed for their correspondence to the additive rule. The enthalpy increments of the thermochemical groups characteristic of the compounds under consideration were calculated, and certain trends in their variations revealed. The conclusion was drawn that the enthalpies of formation of organic derivatives of oximes and hydroxylamines of various structures could be calculated with high accuracy.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Theoretical study of the thermodynamics and kinetics of hydrogen abstractions from hydrocarbons

Thermochemical and kinetic data were calculated at four cost-effective levels of theory for a set consisting of five hydrogen abstraction reactions between hydrocarbons for which experimental data are available. The selection of a reliable, yet cost-effective method to study this type of reactions for a broad range of applications was done on the basis of comparison with experimental data or with results obtained from computationally demanding high level of theory calculations. For this benchmark study two composite methods (CBS-QB3 and G3B3) and two density functional theory (DFT) methods, MPW1PW91/6-311G(2d,d,p) and BMK/6-311G(2d,d,p), were selected. All four methods succeeded well in describing the thermochemical properties of the five studied hydrogen abstraction reactions. High-level Weizmann-1 (W1) calculations indicated that CBS-QB3 succeeds in predicting the most accurate reaction barrier for the hydrogen abstraction of methane by methyl but tends to underestimate the reaction barriers for reactions where spin contamination is observed in the transition state. Experimental rate coefficients were most accurately predicted with CBS-QB3. Therefore, CBS-QB3 was selected to investigate the influence of both the 1D hindered internal rotor treatment about the forming bond (1D-HR) and tunneling on the rate coefficients for a set of 21 hydrogen abstraction reactions. Three zero curvature tunneling (ZCT) methods were evaluated (Wigner, Skodje & Truhlar, Eckart). As the computationally more demanding centrifugal dominant small curvature semiclassical (CD-SCS) tunneling method did not yield significantly better agreement with experiment compared to the ZCT methods, CD-SCS tunneling contributions were only assessed for the hydrogen abstractions by methyl from methane and ethane. The best agreement with experimental rate coefficients was found when Eckart tunneling and 1D-HR corrections were applied. A mean deviation of a factor 6 on the rate coefficients is found for the complete set of 21 reactions at temperatures ranging from 298 to 1000 K. Tunneling corrections play a critical role in obtaining accurate rate coefficients, especially at lower temperatures, whereas the hindered rotor treatment only improves the agreement with experiment in the high-temperature range.     

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.     

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.     

Enthalpies of formation and substituent effects of ortho-, meta-, and para-aminotoluenes from thermochemical measurements and from ab initio calculations

The molar enthalpies of vaporization of 2-amino-, 3-amino-, and 4-aminotoluenes were obtained from the temperature dependence of the vapor pressure measured by the transpiration method. The molar enthalpy of sublimation of 4-aminotoluene was measured in the same way. The standard (p(o) = 0.1 MPa) molar enthalpy of formation delta fH(o)m(cr) at the temperature T = 298.15 K of crystalline 4-aminotoluene was measured using combustion calorimetry. The thermochemical investigations of aminotoluenes available in the literature were collected and combined with our own experimental results to obtain our own reliable standard molar enthalpies of formation at T = 298.15 K in the gaseous state. Ab initio calculations of aminotoluenes have been performed using the MP2Full/6-31G(d) and G3(MP2) basis sets, and the results from the bond separation method are in excellent agreement with the experiment. These new results help to resolve the uncertainty in the available thermochemical data on aminotoluenes. Weak mutual interactions of substituents in aminotoluenes have been realized using an isodesmic reaction procedure.     

Gas-phase thermochemical properties of some tri-substituted phenols: A density functional theory study

The study of the energetics of phenolic compounds has a considerable practical interest since this family of compounds includes numerous synthetic and naturally occurring antioxidants. In this work, density functional theory (DFT) has been used to investigate gas-phase thermochemical properties of the following tri-substituted phenols: 2,4,6-trimethylphenol, 2,6-dimethyl-4-tert-butylphenol, 2,6-dimethyl-4-methoxyphenol, 2,4,6-tri-tert-butylphenol, 2,6-di-tert-butyl-4-methylphenol, 2,6-di-tert-butyl-4-methoxyphenol, 2,4,6-trimethoxyphenol, 2,6-dimethoxy-4-methylphenol and 2,6-dimethoxy-4-tert-butylphenol. Molecular structures were computed with the B3LYP and the ωB97X-D functionals and the 6-31G(d) basis set. More accurate energies were obtained from single-point energy calculations with both functionals and the 6-311++G(2df,2pd) basis set. Standard enthalpies of formation of the phenolic molecules and phenoxyl radicals were derived using an appropriate homodesmotic reaction. The OH homolytic bond dissociation enthalpies, gas-phase acidities and adiabatic ionization enthalpies were also calculated. The general good agreement found between the calculated and the few existent experimental gas-phase thermochemical parameters gives confidence to the estimates concerning the phenolic compounds which were not yet experimentally studied.     

Computational study of the thermochemistry of organophosphorus(III) compounds

The enthalpies of formation of organophosphorus(III) compounds have been calculated at the G3X, G3X(MP2), and B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) levels of theory using the atomization energy procedure and the method of isodesmic reactions. The Delta f H298 degree values for 50 relatively large molecules with up to 10 non-hydrogen atoms, such as P(CH3)3, P(C2H5)3, P(OCH3)3, n-C4H9OPCl2, [(CH3)2N]2PCl, (C2H5)2NPCl2, and [(CH3)2N]2PCN, have been calculated directly from the G3X atomization energies. A good agreement between the known experimental values and G3X results for 14 compounds provides support to our predictions for remaining species whose experimental enthalpies of formation are unknown or known with relatively large uncertainties. On the basis of our calculations and sometimes conflicting experimental data a set of internally consistent enthalpies of formation has been recommended for organophosphorus(III) compounds. Our computational results call into question the experimental enthalpies of formation of P(C2H5)3 and P(n-C4H9)3. From comparison with most reliable experimental data, the accuracy of the theoretical enthalpies of formation is estimated as ranging from 5 to 10 kJ/mol. The recommended Delta f H298 degree values were used to derive the group additivity values (GAVs) for 45 groups involving the phosphorus(III) atom. These GAVs significantly extend the applicability of Benson's group additivity method and may be used to estimate the enthalpies of formation of larger organophosphorus(III) compounds, where high level quantum chemical calculations are impracticable.