Comparison of semiempirical MO methods applied to large molecules

The heats of formation for 19 molecules have been calculated with PM3 and AM1 semiempirical methods. The values obtained have been compared with experimental heats of formation. With PM3 and AM1 the average differences between calculated and experimental heats of formation are 8.45 and 12.34 kcal mol^?1 respectively. There are significant differences when large molecules are considered: this suggests that the parameterization should be done including larger molecules.     

Quantum chemical study of the formation of acetylenic alcohols and 7-methylidene-6,8-dioxabicyclo[3.2.1]octanes from acetylene and carbonyl compounds

The geometric structures of the conformers of 2,3,4-trimethyl-7-methylidene-1,5-di-(thiophen-2-yl)-6,8-dioxabicyclo[3.2.1]octane were studied by the MP2/6-311++G**// B3LYP/6-31G* method. The average deviations of the calculated bond lengths from those determined by X-ray diffraction are no larger than 0.02 Å, and the deviations of the bond angles and dihedral angles are 0.7 and 3.5°, respectively. In the gas phase, the chair conformation is thermodynamically ～7.3 kcal mol^?1 more favorable than the boat conformation, and, taking into account the solvent effect of DMSO, the chair conformation is ～9.4 kcal mol^?1 more favorable. For formaldehyde and a series of ketones, including alkyl hetaryl ketones, the enthalpies of competitive formation of acetylenic alcohols and 7-methylidene-6,8-dioxabicyclo[3.2.1]octanes were calculated. The formation of the latter compounds is characterized by a considerable decrease in the enthalpy (to ?91 kcal mol^?1).     

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.     

Die Bildungsenthalpie von Zinndijodid,SnJ2 (c)

The heat of reaction for SnJ₂ (c)+J₂ (c)+4045 CS₂ (l)=[SnJ₄; 4045 CS₂] (sol) has been determined to be (?41.12±0.55) kJ mol^?1, [(?9.83±0.13) kcal mol^?1] by isoperibol solution calorimetry. Combining this result with the heat of formation of SnJ₄ in CS₂ determined in a previous investigation¹¹ the value (?153.9±1.40) kJ mol^?1, [(?36.9±0.33) kcal mol^?1] has been derived for the heat of formation, ΔH _f ^ι (SnJ₂;c; 298.15 K), of tin diiodide.     

A photoionization study of [C4H7]+ ions in the gas phase

The ionization and [C₄H₇]⁺ appearance energies for a series of C₄H₇CI and C₄H₇Br isomers have been measured by dissociative photoionization mass spectrometry. Cationic heats of formation, based on the stationary electron convention, are derived. No threshold ion is observed with a heat of formation corresponding to the trans-1-methylallyl cation, although there is evidence for formation of the less stable cis isomer. A 298 K heat of formation of 871 kJ mol^?1 is obtained for the cyclopropylcarbinyl cation, with the cyclobutyl cation having a higher value of 886 kJ mol^?1. At the HF/6-31G^** level, ab initio molecular orbital calculations show the 2-butenyl, isobutenyl and homoallyl cations to be stable forms of [C₄H₇]⁺, being less stable than the trans-1-methylallyl cation by 101 kJ mol^?1, 159 kJ mol^?1 and 164 kJ mol^?1, respectively. However, threshold formation is not observed for any of these ions, the fragmentation of appropriate precursor molecules producing [C₄H₇]⁺ ions with lower energy structures.     

Kinetics and thermochemistry of the methyl radical: Study of the CH3 + HCl reaction

The kinetics of the reaction between CH₃ and HCl was studied in a tubular reactor coupled to a photoionization mass spectrometer. Rate constants were measured as a function of temperature (296–495 K) and were fitted to an Arrhenius expression: k₁ = 5.0(±0.7) × 10^?13 exp{?1.4(±0.3) kcal mol^?1/RT} cm³ molecule^?1 s^?1. This information was combined with known kinetic parameters of the reverse reaction to obtain Second Law determinations of the methyl radical heat of formation {34.7(±0.6) kcal mol^?1} and entropy {46(±2) cal mol^?1 K^?1} at 298 K. Using the known entropy of CH₃, a more accurate Third Law determination of the CH₃ heat of formation at this temperature was also obtained {34.8(±0.3) kcal mol^?1}. The values of k₁ obtained in this study are between those reported in prior investigations. The results were also used to test the accuracy of the thermochemical information which can be obtained from kinetic studies of R + HX (X = Cl, Br, I) reactions of the type described here.     

Thermal properties of some cyclic disulfides: naphthalene disulfide and diphenylene disulfide

The specific heat, the melting heat and entropy, the vaporization heat of naphtalene disulfide (C₁₀H₆S₂) and of diphenylene disulfide (C₁₂H₈S₂) have been determined by differential scanning calorimetry (DSC).Over the temperature range examined the specific heat may be represented as follows:

where T is the temperature in degrees Kelvin, while melting heat, vaporization heat, melting entropy are for naphtalene disulfide: 3.10 kcal mol^?1, 6.42 kcal mol^?1, 7.87 cal deg^? mol^?1 and for diphenylene disulfide: 4.62 kcal mol^?1, 6.90 kcal mol^?1 and 11.87 cal deg^?1 mol^?1.     

Exploring new species on the [H,S, Se,Cl] potential energy surface

This investigation is concerned with the characterization of seleno‐sulfide‐halogen model systems, the isomerization processes, and the dissociation into diatomic fragment channels on the [H, S, Se, Cl] potential energy surface. Structural, energetic, and vibrational data were obtained at the CCSD(T) and MP2 levels of theory with the series of correlation consistent basis sets and extrapolated to the complete basis set (CBS) limit. For the frequencies, additional computations were performed to include the contribution of anharmonic effects, and for the determination of the heats of formation, important corrections incorporating core‐valence correlation effects and relativistic effects (scalar and spin‐orbit) were taken into account. CCSD(T)/CBS relative stability (kcal mol^?1) follows the order: HSSeCl (0.0), HSeSCl (8.80), SSeHCl (23.52), and SeSHCl (25.87). The cis‐rotational barrier for the two lowest isomers is practically identical (10.14 and 10.09 kcal mol^?1), whereas for the trans barrier, we obtained 9.25 (HSSeCl) and 8.45 (HSeSCl) kcal mol^?1. Dissociation of HSSeCl (HSeSCl) into HS (HSe) + SeCl (SCl) requires 59.70 (56.30) kcal mol^?1. For the most stable isomer, we predict a value of the heat of formation at 298.15 K of 2.53 kcal mol^?1. One of the outcomes of this research is that the MP2 results are consistent with those of CCSD(T). The MP2 method turns out to be a reliable alternative for a first exploration of larger catenated species, although it accounts for a lesser fraction of correlation effects. © 2012 Wiley Periodicals, Inc.     

Different catalysis role of in‐loop and out‐of‐loop waters in assisting HNS/HSN proton transfer isomerizations: Bridging vs. surrounding effect

In this work, a density function theory (DFT) study is presented for the HNS/HSN isomerization assisted by 1–4 water molecules on the singlet state potential energy surface (PES). Two modes are considered to model the catalytic effect of these water molecules: (i) water molecule(s) participate directly in forming a proton transfer loop with HNS/HSN species, and (ii) water molecules are out of loop (referred to as out‐of‐loop waters) to assist the proton transfer. In the first mode, for the monohydration mechanism, the heat of reaction is 21.55 kcal · mol^?1 at the B3LYP/6‐311++G** level. The corresponding forward/backward barrier lowerings are obtained as 24.41/24.32 kcal · mol^?1 compared with the no‐water‐assisting isomerization barrier T (65.52/43.87 kcal · mol^?1). But when adding one water molecule on the HNS, there is another special proton‐transfer isomerization pathway with a transition state 10T′ in which the water is out of the proton transfer loop. The corresponding forward/backward barriers are 65.89/65.89 kcal · mol^?1. Clearly, this process is more difficult to follow than the R–T–P process. For the two‐water‐assisting mechanism, the heat of reaction is 19.61 kcal · mol^?1, and the forward/backward barriers are 32.27/12.66 kcal · mol^?1, decreased by 33.25/31.21 kcal · mol^?1 compared with T. For trihydration and tetrahydration, the forward/backward barriers decrease as 32.00/12.60 (30T) and 37.38/17.26 (40T) kcal · mol^?1, and the heat of reaction decreases by 19.39 and 19.23 kcal · mol^?1, compared with T, respectively. But, when four water molecules are involved in the reactant loop, the corresponding energy aspects increase compared with those of the trihydration. The forward/backward barriers are increased by 5.38 and 4.66 kcal · mol^?1 than the trihydration situation. In the second mode, the outer‐sphere water effect from the other water molecules directly H‐bonded to the loop is considered. When one to three water molecules attach to the looped water in one‐water in‐loop‐assisting proton transfer isomerization, their effects on the three energies are small, and the deviations are not more than 3 kcal · mol^?1 compared with the original monohydration‐assisting case. When adding one or two water molecules on the dihydration‐assisting mechanism, and increasing one water molecule on the trihydration, the corresponding energies also are not obviously changed. The results indicate that the forward/backward barriers for the three in‐loop water‐assisting case are the lowest, and the surrounding water molecules (out‐of‐loop) yield only a small effect. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Remarks on the analysis of torsional energy surfaces of cycloheptane and cyclooctane by molecular mechanics

A modified version (MM 2′) of the Allinger's 1977 force field is checked against cycloheptane and cyclooctane. Cycloheptane is characterized by two pseudorotating itineraries, chair/twist-chair and boat/twist-boat, separated by a barrier of 8.5 kcal mol^?1. The activation energy in the C/TC pseudorotation is estimated to be 0.96 kcal mol^?1, while B and TB transform into each other freely at an energy level 3.8 kcal mol^?1 above the global energy minimum (TC). With cyclooctane the lowest energy is calculated for the boat-chair form which participates in a pseudorotational process with TBC through a saddle point lying 3.5 kcal mol^?1 above BC. The chair/chair and boat/boat families contain only one local minimum, crown and BB, respectively, on the MM 2′ surface. The results are presented as an illustration for quick coverage of torsional energy surface by two-bond driver calculation with the block-diagonal Newton–Raphson minimization, followed by the force search of stationary points by full-matrix Newton–Raphson optimization.     

Molecular mechanics (MM4) calculations on carbonyl compounds part I: aldehydes

Aliphatic aldehydes have been studied with the aid of the MM4 force field. The structures, moments of inertia, vibrational spectra, conformational energies, barriers to internal rotation, and dipole moments have been examined for six compounds (nine conformations). MM4 parameters have been developed to fit the indicated quantities to the wide variety of experimental data. Ab initio (MP2) and density functional theory (B3LYP) calculations have been used to augment and/or replace experimental data, as appropriate. Because more, and to some extent, better, data have become available since MM3 was developed, it was anticipated that the overall accuracy of the information calculated with MM4 would be better than with MM3. The best single measure of the overall accuracy of a force field is the accuracy to which the moments of inertia of a set of compounds (from microwave spectroscopy) can be reproduced. For all of the 20 moments (seven conformations) experimentally known for the aldehyde compounds, the MM4 rms error is 0.30%, while with MM3, the most accurate force field presently available, the rms error over the same set is 1.01%. The calculation of the vibrational spectra was also improved overall. For the four aldehydes that were fully analyzed (over a total of 78 frequencies), the rms errors with MM4 and MM3 are 18 and 38 cm^?1, respectively. These improvements came from several sources, but the major ones were separate parameters involving the carbonyl carbon for formaldehyde, the alkyl aldehydes and the ketones, and new crossterms featured in the MM4 force field that are not present in the MM3 version. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1396–1425, 2001     

Molecular mechanics calculations on carbonyl compounds. II. Open‐chain ketones

Open‐chain aliphatic ketones were studied with the molecular mechanics (MM4) force field. A total of seven compounds were examined. Structures were well fit, including moments of inertia. Rotational barriers, vibrational spectra, and dipole moments were also well fit. The overall root mean square errors for MM3 and MM4 were 0.27 and 0.18%, respectively, for the six moments of inertia (known experimentally for two compounds) and 31 and 20 cm^?1, respectively, for the vibrational frequencies (over 99 weighted modes). © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1426–1450, 2001     

The kinetics of olation of some hydroxoaquotetrammine-cobalt(III) complexes in the solid state

The solid-state kinetics for the olation reactions of [Co(NH₃₄(OH)(H₂O)]X₂ (where X  Cl^?, Br^?, or $\text{1}\text{2}$SO^2?₄ were determined by several different methods using dynamic and isothermal thermogravimetric data. For the reduced-time plot method, E values were 20, 43, and 25 kcal mol^?1 for the chloride-bromide, and sulfate complexes, respectively. For the Jacobs and Kureishy method, E values of 21, 37, and 17 kcal mol^?1, were obtained for the above three complexes, respectively. A possible reaction pathway is suggested for the olation reaction.     

Sels de Pyrylium. XVIII. Etude par RMN du Carbone-13 de Sels d'Aminopyrylium et des Barrières de Rotation dans quelques Cations N,N-Diméthylaminopyrylium

The C-2—N bond of 2-N,N-dimethylaminopyrylium cations has a partial π character due to the conjugation of the nitrogen lone-pair with the ring. The values of ΔG^≠, ΔH^≠, ΔS^≠ parameters related to the corresponding hindered rotation have been determined by ¹³C NMR total bandshape analysis. This conjugation decreases the electrophilic character of carbon C-4 so that the displacement of the alkoxy group is no longer possible. Such a hindered rotation also exists in 4-N,N-dimethylaminopyrylium cations and the corresponding ΔG^≠ parameters have been evaluated. Comparison of these two cationic species shows that hindered rotation around the C—N bond is larger in position 4 than in position 2. Furthermore, the barrier to internal rotation around the C-2? N bond decreases with increasing electron donating power of the substituent at position 4. ΔG^≠ values decreases from 19.1 kcal mol^?1 (79.9 kJ mol^?1) to 12.6 kcal mol^?1 (52.7 kJ mol^?1) according to the following sequence for the R-4 substituents: -C₆H₅, -CH₃, -OCH₃, -N(CH₃)₂.     

The enthalpies of formation of CH3Mn(CO)5 and of CH3Re(CO)5, and the strengths of the CH3Mn and CH3Re bonds

From measurements of the heats of iodination of CH₃Mn(CO)₅ and CH₃Re(CO)₅ at elevated temperatures using the ‘drop’ microcalorimeter method, values were determined for the standard enthalpies of formation at 25° of the crystalline compounds: ΔH^o_f[CH₃Mn(CO)₅, c] = ?189.0 ± 2 kcal mol^?1 (?790.8 ± 8 kJ mol^?1), ΔH^o_f[Ch₃Re(CO)₅,c] = ?198.0 ± kcal mol^?1 (?828.4 ± 8 kJ mo^?1). In conjunction with available enthalpies of sublimation, and with literature values for the dissociation energies of MnMn and ReRe bonds in Mn₂(CO)₁₀ and Re₂(CO)₁₀, values are derived for the dissociation energies: D(CH₃Mn(CO)₅) = 27.9 ± 2.3 or 30.9 ± 2.3 kcal mol^?1 and D(CH₃Re(CO)₅) = 53.2 ± 2.5 kcal mol^?1. In general, irrespective of the value accepted for D(MM) in M₂(CO)₁₀, the present results require that, D(CH₃Mn) = $\text{1}\text{2}$D(MnMn) + 18.5 kcal mol^?1 and D(CH₃Re) = $\text{1}\text{2}$D(ReRe) + 30.8 kcal mol^?1.     

Addition of ethyl radicals to carbon monoxide. Kinetic and thermochemical properties of the propionyl radical

Azoethane was irradiated in the presence of carbon monoxide in the temperature range of 238 to 378 K. Kinetic parameters for the addition of ethyl radicals to carbon monoxide and for the decomposition of propionyl radicals were determined. The rate constants were found to be log k(cm³ mol^?1 sec^?1) = 11.19 - 4.8/θ and log k(sec^?1) = 12.77 - 14.4/θ, respectively. Estimated thermochemical properties of the propionyl radical are ΔH_f⁰ = -10.6 ± 1.0 kcal mol^?1, S⁰ = 77.3 ± 1.0 cal K^?1 mol^?1, and D(C₂H₅CO? H) = 87.4 kcal mol^?1.     

A Highly Oxidizing and Isolable Oxoruthenium(V) Complex [RuV(N4O)(O)]2+: Electronic Structure,Redox Properties,and Oxidation Reactions Investigated by DFT Calculations

The electronic structure and redox properties of the highly oxidizing, isolable Ru^V?O complex [Ru^V(N₄O)(O)]²⁺, its oxidation reactions with saturated alkanes (cyclohexane and methane) and inorganic substrates (hydrochloric acid and water), and its intermolecular coupling reaction have been examined by DFT calculations. The oxidation reactions with cyclohexane and methane proceed through hydrogen atom transfer in a transition state with a calculated free energy barrier of 10.8 and 23.8 kcal mol^?1, respectively. The overall free energy activation barrier (ΔG^≠=25.5 kcal mol^?1) of oxidation of hydrochloric acid can be decomposed into two parts: the formation of [Ru^III(N₄O)(HOCl)]²⁺ (ΔG=15.0 kcal mol^?1) and the substitution of HOCl by a water molecule (ΔG^≠=10.5 kcal mol^?1). For water oxidation, nucleophilic attack on Ru^V?O by water, leading to O? O bond formation, has a free energy barrier of 24.0 kcal mol^?1, the major component of which comes from the cleavage of the H? OH bond of water. Intermolecular self‐coupling of two molecules of [Ru^V(N₄O)(O)]²⁺ leads to the [(N₄O)Ru^IV? O₂? Ru^III(N₄O)]⁴⁺ complex with a calculated free energy barrier of 12.0 kcal mol^?1.     

An ab initio study of the structures and energetics of CH3OCl and CH3ClO

The equilibrium structures, vibrational spectra, and heats of formation for CH₃OCl and CH₃ClO have been estimated using high levels of ab initio molecular orbital theory. The lowest energy isomer is found to be CH₃OCl, and its heat of formation is estimated to be −13.5±2 kcal mol⁻¹, in good agreement with bond additivity estimates. Results for the CH₃ClO isomer are presented for the first time, and it is found to be 40.5 kcal mol⁻¹ higher in energy relative to CH₃OCl. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 29–35, 1999     

Substituenteneffekte auf die C–C-Bindungsstärke, 8. Standardbildungsenthalpien und Spannungsenthalpien von 1,1,2,2-Tetraethylethylenglykol-dimethylether und D,L-1,2-Dimethyl-1,2-diphenylethylenglykol-dimethylether

Effects of Substituents on the Strength of C - C Bonds, 8¹. - Heats of Formation and Strain of 1,1,2,2-Tetraethylethylene Glycol Dimethyl Ether and D,L .-1,2-Dimethyl-l,2-diphenylethylene Glycol Dimethyl Ether The heats of combustion of the title compounds 1 and 2 were measured calorimetrically with the result (kcal mol ^-1, s. d. in parentheses) ΔH°_c = − 1880.1 (± 0.6) and − 2373.3 (± 1.4). The heat of vaporisation of 1 ΔH^v = 14.3 (± 0.3) and the heat of sublimation of 2 ΔH^sub = 27.2 (± 0.5) were derived from their temperature dependance of the vapor pressure. The latter were determined between 30 and 80°C using a flow method. The resulting standard heats of formation ΔH°_t(g) = −122.4 (± 0.7) and −43.8 (±1.5) for 1 and 2 correspond to a strain enthalpy (_s) of 15.9 and 8.0 kcal mol^-1, respectively. The steric strain of the dimethoxyethanes 1 and 2 is about one fourth lower than the strain of the corresponding dimethylethanes 3 and 4 bearing the same substituents. Thus, a methoxy group causes less steric stress than a methyl group.