Thermochemistry,bond energies and internal rotor barriers of methyl sulfinic acid,methyl sulfinic acid ester and their radicals

Sulfur–Oxygen containing hydrocarbons are formed in oxidation of sulfides and thiols in the atmosphere, on aerosols and in combustion processes. Understanding their thermochemical properties is important to evaluate their formation and transformation paths. Structures, thermochemical properties, bond energies, and internal rotor potentials of methyl sulfinic acid CH₃S(?O)OH, its methyl ester CH₃S(?O)OCH₃ and radicals corresponding to loss of a hydrogen atom have been studied. Gas phase standard enthalpies of formation and bond energies were calculated using B3LYP/6‐311G (2d, p) individual and CBS‐QB3 composite methods employing work reactions to further improve accuracy of the ${\Delta} _{{\bf f}} H_{{\bf 298}}^{{\bf o}} $ . Molecular structures, vibration frequencies, and internal rotor potentials were calculated. Enthalpies of the parent molecules CH₃S(?O)OH and CH₃S(?O)OCH₃ are evaluated as ?77.4 and ?72.7 kcal mol^?1 at the CBS? QB3 level; Enthalpies of radicals C^?H₂? S(?O)? OH, CH₃? S^?(?O)₂, C^?H₂? S(?O)? OCH₃ and CH₃? S(?O)? OC^?H₂ (CBS‐QB3) are ?25.7, ?52.3, ?22.8, and ?26.8 kcal mol^?1, respectively. The CH₃C(?O)O—H bond dissociation energy is of 77.1 kcal mol^?1. Two of the intermediate radicals are unstable and rapidly dissociate. The CH₃S(?O)? O. radical obtained from the parent CH₃? S(?O)? OH dissociates into methyl radical (${\bf CH}_{{\bf 3}}^{{\bf .}} $ ) plus SO₂ with endothermicity (ΔH_rxn) of only 16.2 kcal mol^?1. The CH₃? S(?O)? OC^?H₂ radical dissociates into CH₃? S^?=O and CH₂=O with little or no barrier and an exothermicity of ?19.9 kcal mol^?1. DFT and the Complete Basis Set‐QB3 enthalpy values are in close agreement; this accord is attributed to use of isodesmic work reactions for the analysis and suggests this combination of B3LYP/work reaction approach is acceptable for larger molecules. Copyright © 2010 John Wiley & Sons, Ltd.     

Ab initio molecular dynamics study of aqueous formaldehyde and methanediol

The solvation shell of aqueous formaldehyde has been studied by ab initio molecular dynamics. Two different DFT approaches using BLYP and PBE functionals were explored. The results show only a slightly different mobility in the solvation shells and allow characterization of the hydrogen bonded structure with a H₂C?=?O··HOH hydrogen bond lifetime of ca. 3 ps. Formaldehyde hydrolysis was studied by following the reverse process, methanediol decomposition, by Blue Moon constrained MD showing that four water molecules are directly involved in the reaction and assisted by the whole hydration shell. The total energy of the aqueous methanediol to formaldehyde inter-conversion process is calculated with a barrier height of ca. 95?kJ?mol^?1 while the corresponding free energy barrier is only ΔG^??=?46?kJ?mol^?1 at 300?K.     

Computational study on structures,thermochemical properties,and bond energies of disulfide oxygen (S–S–O)‐bridged CH3SSOH and CH3SS(=O)H and radicals

Cleavage of disulfide bonds is a common method used in linking peptides to proteins in biochemical reactions. The structures, internal rotor potentials, bond energies, and thermochemical properties (Δ_fH°, S°, and Cp(T)) of the S–S bridge molecules CH₃SSOH and CH₃SS(=O)H and the radicals CH₃SS?=O and C?H₂SSOH that correspond to H‐atom loss are determined by computational chemistry. Structure and thermochemical parameters (S° and Cp(T)) are determined using density functional Becke, three‐parameter, Lee–Yang–Parr (B3LYP)/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p). The enthalpies of formation for stable species are calculated using the total energies at B3LYP/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p), and the higher level composite CBS–QB3 levels with work reactions that are close to isodesmic in most cases. The enthalpies of formation for CH₃SSOH, CH₃SS(=O)H are ?38.3 and ?16.6 kcal mol^?1, respectively, where the difference is in enthalpy RSO–H versus RS(=O)–H bonding. The C–H bond energy of CH₃SSOH is 99.2 kcal mol^?1, and the O–H bond energy is weaker at 76.9 kcal mol^?1. Cleavage of the weak O–H bond in CH₃SSOH results in an electron rearrangement upon loss of the CH₃SSO–H hydrogen atom; the radical rearranges to form the more stable CH₃SS· = O radical structure. Cleavage of the C–H bond in CH₃SS(=O)H results in an unstable [CH₂SS(=O)H]* intermediate, which decomposes exothermically to lower energy CH₂ = S + HSO. The CH₃SS(=O)–H bond energy is quite weak at 54.8 kcal mol^?1 with the H–C bond estimated at between 91 and 98 kcal mol^?1. Disulfide bond energies for CH₃S–SOH and CH3S–S(=O)H are low: 67.1 and 39.2 kcal mol^?1. Copyright © 2011 John Wiley & Sons, Ltd.     

Ab Initio Investigation of the Microspecies and Energy in Hydrated Strontium Ion Clusters

Quantum chemistry calculations were used to study the structure and energy of strontium (Sr) ion hydrated clusters [Sr(H₂O)_1?25]²⁺. The saturated hydration number of the first hydration layer of Sr²⁺ was 8, and the hydration distance was 2.58 Å. The second hydration layer had 1–9 hydration numbers, and the hydration distance was in the range of 4.4–4.6 Å. This work also developed the relationship between the thermodynamic data (average water binding energy E_n and successive water binding energy ΔE_n,n?1, etc.) of the aforementioned low-energy structure and the hydration structures. The first hydration layer was formed by the strong electrostatic interaction between Sr²⁺ and water molecules, and the decrease in ΔE_n,n?1 was relatively large. Hydrogen bonds were formed between water molecules of the second hydration layer and water molecules of the inner layer, and the decrease in ΔE_n,n?1 was relatively small. When one water molecule was added beyond the second hydration layer, ΔE_n,n?1 was close to the hydrogen bond energy 8.88 kcal/mol (37.1 kJ/mol) of dimer water molecule, indicating that there was very weak interaction between Sr²⁺ and the water molecules beyond the second hydration layer.     

Structures,vibrational spectra,and relative energetics of FC(O)ONO and FC(O)NO2 isomers at DFT and ab initio levels

The structural and vibrational parameters of FC(O)ONO and FC(O)NO₂ isomers were examined theoretically using the B3LYP/6-311+G(3df) and CCSD(T)/6-311G(d) methods. Four conformers of FC(O)ONO isomer and one FC(O)NO₂ isomer are found here. Among them, the trans–cis and cis–cis FC(O)ONO configuration are new conformers. The energetics were refined with G3//B3LYP and CBS-QB3 calculations. The trans–trans conformer of the FC(O)ONO isomer is found to be the lowest energy structure, with an estimated heat of formation of ?104.9 kcal mol^?1 at 0 K as determined from CBS-QB3 theory. The next lowest structure is the cis–trans FC(O)ONO lying 1.7 kcal mol^?1 above the trans–trans structural form. The highest energy structure is the FC(O)NO₂ isomer with a predicted heat of formation of ?84.8 kcal mol^?1. A comparison of the relative stability of the FCNO₃ isomers with the isomers of ClCNO₃ shows that the Cl analogues follow the same pattern of stability, as do the F isomers. However, the chlorine isomers are unstable relative to their fluorine analogues.     

The heats of formation of the haloacetylenes XCCY [X,Y = H,F, Cl]: basis set limit ab initio results and thermochemical analysis

The heats of formation of haloacetylenes are evaluated using the recent W1 and W2 ab initio computational thermochemistry methods. These calculations involve CCSD and CCSD(T) coupled cluster methods, basis sets of up to spdfgh quality, extrapolations to the one-particle basis set limit, and contributions of inner-shell correlation, scalar relativistic effects. and (where relevant) first-order spin-orbit coupling. The heats of formation determined using W2 theory are: δH₁ ²⁹⁸(HCCH) = 54.48 kcal mol^?1, δH_f ²⁹⁸(HCCH) = 25.15 kcal mol, δH_f ²⁹⁸(FCCF) = 1.38 kcal mol^?1, δH_f ²⁹⁸(HCCC1) = 54.83 kcal mol^?1, δH_f ²⁹⁸(CICCC1) = 56.21 kcal mol^?1, and δH_f ²⁹⁸(FCCC1) = 28.47 kcal mo1^?1. Enthalpies of hydrogenation and destabilization energies relative to acetylene were obtained at the WI level of theory. So doing we find the following destabilization order for acetylenes: FCCF > ClCCF > HCCF > ClCCCl > HCCCI > HCCH. By a combination of WI theory and isodesmic reactions. we show that the generally accepted heat of formation of 1,2-dichloroethane should be revised to ?31.8 ± 0.6 kcal mol^?1, in excellent agreement with a very recent critically evaluated review. The performance of compound thermochemistry schemes, such as G2, G3, G3X and CBS-QB3 theories, has been analysed.     

Accurate ab initio anharmonic force field and heat of formation for silane

From large basis set coupled cluster calculations and a minor empirical adjustment, an anharmonic force field for silane has been derived that is consistently of spectroscopic quality (±1 cm^?1 on vibrational fundamentals) for all isotopomers of silane studied. Inner-shell polarization functions have an appreciable effect on computed properties and even on anharmonic corrections. From large basis set coupled cluster calculations and extrapolations to the infinite-basis set limit, we obtain TAE₀ = 303.80 ± 0.18 kcal mol^?1, which includes an anharmonic zero-point energy (19.59 kcal mol^?1), inner-shell correlation (—0.36 kcal mol^?1), scalar relativistic corrections (— 0.70 kcal mol^?1) and atomic spin-orbit corrections (—0.43 kcal mol^?1). In combination with the recently revised ΔH ^o _{f, o}[Si(g)], we obtain ΔH ^o _f.o[SiH₄(g)] = 9.9 ± 0.4 kcal mol^?1 in between the two established experimental values.     

The keto-enol equilibrium in substituted acetaldehydes: focal-point analysis and ab initio limit

High-level ab initio electronic structure calculations up to the CCSD(T) theory level, including extrapolations to the complete basis set (CBS) limit, resulted in high precision energetics of the tautomeric equilibrium in 2-substituted acetaldehydes (XH₂C-CHO). The CCSD(T)/CBS relative energies of the tautomers were estimated using CCSD(T)/aug-cc-pVTZ, MP3/aug-cc-pVQZ, and MP2/aug-cc-pV5Z calculations with MP2/aug-cc-pVTZ geometries. The relative enol (XHC?=?CHOH) stabilities (ΔE _{e,CCSD(T)/CBS}) were found to be 5.98?±?0.17, ?1.67?±?0.82, 7.64?±?0.21, 8.39?±?0.31, 2.82?±?0.52, 10.27?±?0.39, 9.12?±?0.18, 5.47?±?0.53, 7.50?±?0.43, 10.12?±?0.51, 8.49?±?0.33, and 6.19?±?0.18?kcal?mol^?1 for X?=?BeH, BH₂, CH₃, Cl, CN, F, H, NC, NH₂, OCH₃, OH, and SH, respectively. Inconsistencies between the results of complex/composite energy computations methods Gn/CBS (G2, G3, CBS-4M, and CBS-QB3) and high-level ab initio methods (CCSD(T)/CBS and MP2/CBS) were found. DFT/aug-cc-pVTZ results with B3LYP, PBE0 (PBE1PBE), TPSS, and BMK density functionals were close to the CCSD(T)/CBS levels (MAD?=?1.04?kcal?mol^?1).     

The reactions of Cl+ (3P) and Cl+ (1D) with hydrogen sulphide. A G2 molecular orbital study

G2 ab initio molecular orbital calculations have been performed to study the potential energy surfaces (PESs) associated with the reactions of Cl⁺ in its ³P ground state and in its ¹D first excited state with hydrogen sulphide. [H₂, Cl, S]⁺ singlet and triplet state cations present very different bonding characteristics. The latter are systematically ion-dipole or hydrogen-bonded weakly bound species, while the former are covalent molecular ions. As a consequence, although the Cl⁺(³P) is 34.5 kcal mol^?1 more stable than Cl⁺(¹D), the global minimum of the singlet PES lies 37.3 kcal mol^?1 below the global minimum of the triplet PES. Both singlet and triplet potential energy surfaces show significant differences with respect to those associated with Cl⁺ + H₂O reactions as well as with SH₂ reactions with F⁺. In both cases, the major product should be SH⁺ ₂; SH⁺ and HCl⁺ being the minor products, in agreement with the experimental evidence. The estimated heat of formation for the most stable H₂SCl⁺ singlet state species is 198 ± 1 kcal mol^?1.     

硫酸催化乙二醛气体相水化反应机理和速率常数理论计算研究

采用M06-2X和CCSD(T)高阶量化计算和传统过渡态理论研究硫酸催化乙二醛气体相水化反应.对HCOCHO+H₂O, HCOCHO+H₂O+H₂O, HCOCHO+H₂O+H₂O, HCOCHO+H₂O...H₂SO₄和HCOCHO+H₂O+H₂SO₄五个路径的反应机理和速率常数进行了研究.计算结果表明硫酸具有较强的催化能力,能显著减小乙二醛水化反应的能垒,在CCSD(T)/6-311++G(3df,3pd)//M06-2X/6-311++G(3df,3pd)理论水平,当硫酸分子参与乙二醛水化反应时,反应能垒从37.15 kcal/mol减少至7.08 kcal/mol.在室温条件下,硫酸催化乙二醛水化反应的反应速率1.34×10^-11 cm³/(molecule.s),是等量水分子参与乙二醛水化反应的速率的10¹²倍,大于乙二醛与OH自由基反应的反应速率1.10×10^-11 cm³/(molecule.s).这表明大气条件下,硫酸催化乙二醛水化反应可以发生,同乙二醛与OH自由基反应相竞争.     

氧化石墨烯负载的Pt单原子催化硼胺烷水解机理的理论研究

本文研究了氧化石墨烯负载Pt单原子(Pt₁/Gr-O)催化硼胺烷(NH₃BH₃)全水解反应机理,即一分子的NH₃BH₃生成三分子的氢气(H₂)的过程. 在水解路径中,首先吸附的硼胺烷连续断裂两个B-H键生成第一分子的H₂. 接着,一个H₂O分子与^*BHNH₃基团(^*表示吸附态)反应生成^*BH(H₂O)NH₃,其中伸长的O-H键断裂后形成^*BH(OH)NH₃. 然后,第二个H₂O与^*BH(OH)NH₃反应生成^*BH(OH)(H₂O)NH₃,在指向Pt₁/Gr-O表面的O-H断裂后,生成BH(OH)₂NH₃并脱附到水溶液中. 两个水分子脱氢产生的两个H原子脱附生成第二个H₂分子,且Pt₁/Gr-O催化剂恢复. 脱附后的BH(OH)₂NH₃在水溶液中水解生成第三个H₂分子. 纵观整个水解反应,H₂O分子和^*BHNH₃基团的结合是反应速控步,其反应能垒是16.1 kcal/mol. 因此,Pt₁/Gr-O有希望成为室温催化NH₃BH₃全水解催化剂.     

A high level computational study of the CH4/CF4 dimer: how does it compare with the CH4/CH4 and CF4/CF4 dimers?

The interaction within the methane–methane (CH₄/CH₄), perfluoromethane–perfluoromethane (CF₄/CF₄) methane–perfluoromethane dimers (CH₄/CF₄) was calculated using the Hartree–Fock (HF) method, multiple orders of Møller–Plesset perturbation theory [MP2, MP3, MP4(DQ), MP4(SDQ), MP4(SDTQ)], and coupled cluster theory [CCSD, CCSD(T)], as well as the PW91, B97D, and M06-2X density functional theory (DFT) functionals. The basis sets of Dunning and coworkers (aug-cc-pVxZ, x?=?D, T, Q), Krishnan and coworkers [6-311++G(d,p), 6-311++G(2d,2p)], and Tsuzuki and coworkers [aug(df, pd)-6-311G(d,p)] were used. Basis set superposition error (BSSE) was corrected via the counterpoise method in all cases. Interaction energies obtained with the MP2 method do not fit with the experimental finding that the methane–perfluoromethane system phase separates at 94.5?K. It was not until the CCSD(T) method was considered that the interaction energy of the methane–perfluoromethane dimer (?0.69?kcal?mol^?1) was found to be intermediate between the methane (?0.51?kcal?mol^?1) and perfluoromethane (?0.78?kcal?mol^?1) dimers. This suggests that a perfluoromethane molecule interacts preferentially with another perfluoromethane (by about 0.09?kcal?mol^?1) than with a methane molecule. At temperatures much lower than the CH₄/CF₄ critical solution temperature of 94.5?K, this energy difference becomes significant and leads perfluoromethane molecules to associate with themselves, forming a phase separation. The DFT functionals yielded erratic results for the three dimers. Further development of DFT is needed in order to model dispersion interactions in hydrocarbon/perfluorocarbon systems.     

Homo‐Diels–Alder reaction of a very inactive diene,bicyclo[2,2,1]hepta‐2,5‐diene,with the most active dienophile, 4‐phenyl‐1,2,4‐triazolin‐3,5‐dione. Solvent,temperature, and high pressure influence on the reaction rate

Solvent, temperature, and high pressure influence on the rate constant of homo‐Diels–Alder cycloaddition reactions of the very active hetero‐dienophile, 4‐phenyl‐1,2,4‐triazolin‐3,5‐dione (1), with the very inactive unconjugated diene, bicyclo[2,2,1]hepta‐2,5‐diene (2), and of 1 with some substituted anthracenes have been studied. The rate constants change amounts to about seven orders of magnitude: from 3.95^.10^?3 for reaction (1+2) to 12200 L mol^?1 s^?1 for reaction of 1 with 9,10‐dimethylanthracene (4e) in toluene solution at 298 K. A comparison of the reactivity (ln k₂) and the heat of reactions (?_r‐nH) of maleic anhydride, tetracyanoethylene and of 1 with several dienes has been performed. The heat of reaction (1+2) is ?218 ± 2 kJ mol^?1, of 1 with 9,10‐dimethylanthracene ?117.8 ± 0.7 kJ mol^?1, and of 1 with 9,10‐dimethoxyanthracene ?91.6 ±0.2 kJ mol^?1. From these data, it follows that the exothermicity of reaction (1+2) is higher than that with 1,3‐butadiene. However, the heat of reaction of 9,10‐dimethylanthracene with 1 (?117.8 kJ mol^?1) is nearly the same as that found for the reaction with the structural C=C counterpart, N‐phenylmaleimide (?117.0 kJ mol^?1). Since the energy of the N=N bond is considerably lower (418 kJ/bond) than that of the C=C bond (611 kJ/bond), it was proposed that this difference in the bond energy can generate a lower barrier of activation in the Diels–Alder cycloaddition reaction with 1. Linear correlation (R = 0.94) of the solvent effect on the rate constants of reaction (1+2) and on the heat of solution of 1 has been observed. The ratio of the volume of activation (?V^≠) and the volume of reaction (?V_r‐n) of the homo‐Diels–Alder reaction (1+2) is considered as “normal”: ?V^≠/?V_r‐n = ?25.1/?30.95 = 0.81. Copyright © 2012 John Wiley & Sons, Ltd.     

Infrared spectra and isomeric structures of hydroxide ion-water clusters OH- (H2O)1-5: a comparison with H3O (H2O)1-5

Size-selected hydroxide ion water tetramers and pentamers [OH-(H₂O)_4,5], produced by a supersonic expansion, have been investigated using vibrational predissociation spectroscopy in conjunction with ab initio calculations based on density functional theory (DFT). The observed spectra in the frequency range 2650–3850 cm^?1 show some broad absorption bands attributed to the free and hydrogen bonded OH stretches of OH-(H₂O)_4,5 at an estimated cluster temperature of 170 K. DFT calculations performed at the B3LYP/6-31 + G* level reveal five and eight possible low lying isomeric forms for OH-(H₂O)₄ and OH-(H₂O)₅, respectively. The global minimum isomer of the tetramer is tri-solvated cyclic, which is energetically more stable than the tetra-solvated wheel-shaped form with an OH- ion at its centre. Compact cage-like lowest energy structures are found for the pentamer, in which the water molecules can act either as a single-donor-single-acceptor, as a double-proton-donor, or as a double-donor-single-acceptor in both the firs1183t a1192nd the second solvation shell of the OH- ion core. Interconversion among the isomers appears to be rapid as manifested in the observed spectra dominated by broad and congested absorptions. To understand the nature of spectral broadening and congestion, systemic comparisons of the results are made against those of the corresponding protonated cations, H+(H₂O)_n|1 and the corresponding halide anions, X-(H²O)ⁿ X = F, Cl, Br, and I. It is suggested that the spectral complexities observed for OH-(H₂O)_4,5 are predominantly a result of sampling configurations with a large distribution of O_solvent-O_ion-O_solvent angles and O_solvent … H-O_solvent distances between water molecules in the firs1175t a1181nd/or second hydration shells, together with the existence of more than one isomer in the supersonic expanson and rapid isomeric interconversion among them.     

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C₂H₂–Kr and C₂H₂–Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n–1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C₂H₂ moiety corresponds to distances of 4.064 and 4.229?Å, and angles of 65.22° and 68.67° for C₂H₂–Kr and C₂H₂–Xe, respectively. The interaction energy of both complexes is estimated to be ?151.88 (1.817?kJ?mol^?1) and ?182.76?cm^?1 (2.186?kJ?mol^?1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed.     

Phase transition in hexacelsian at about 580 K

Temperature-induced structure and microstructure changes in hexacelsians (BaAl₂Si₂O₈) that have been synthesised from the Ba-exchanged LTA and FAU zeolites (hexacelsian_LTA and hexacelsian_FAU) show that the phase transition near 580?K exists only in hexacelsian_LTA. The X-ray powder diffraction method has been used to follow the evolution of the structure during the phase transition, as described here. The excess thermodynamic quantities Gibbs free energy (G), entropy (S) and enthalpy (H) are obtained through the Landau theory of phase transition. The constants of proportionality between the G and ordering parameter (Q) are: h?=??170345?J?mol^?1, a?=??66.6?J?mol^?1?K^?1 and b?=??410534?J?mol^?1. The abrupt change in the trigonal distortion of the single six-member tetrahedral [SiO₄]^4? and [AlO₄]^5? ring near 580?K is responsible for the phase transition. The phase transition is non-convergent, ferroelastic, pure and proper.     

Effect of protonation on individual hydrogen bonds in the 8-oxoguanine-cytosine base pair: NMR,NBO and AIM analyses

Protonation increases the total binding energy of the 8-oxoguanine-cytosine (8OG:C) base pair by 60–70% at the B3LYP/6-311++G(d,?p) level of theory. It changes the individual H-bond energies, estimated from electron charge densities at bond critical points, by 1.16 to ?16.41?kcal?mol^?1. The individual H-bond energies and the two bond X–Y spin–spin coupling constants (^2hJ_X–Y) increase with protonation where 8OG behaves as an H-bond donor; the reverse is true for the H-bonds in which the 8OG unit acts as an H-bond acceptor. Similar to ^2hJ_X–Y, the value of ^1hJ_O–H (a one-bond H?···?Y spin–spin coupling constant) is distance dependent and in linear correlation with the O?···?H distance, but the ^1hJ_N–H values are independent of the N–H distance and the PSO term is the predominant portion in it. The ¹J_X–H spin–spin coupling constant is dominated by the negative FC term for all hydrogen bonds, although the PSO term is the best to investigate the behaviour of ¹J_X–H across the X–H?·?Y H-bond.     

Theoretical studies of ethylene addition to transition metal compounds with carbene and oxo groups LnM(CH2)(O)

Quantum chemical calculations using density functional theory at the B3LYP level in combination with relativistic effective core potentials for the metals and TZ2P valence basis sets have been carried out for elucidating the reaction pathways of ethylene addition to MeReO₂(CH₂) ( C1 ). The results are compared with our previous studies of ethylene addition to OsO₂(CH₂)₂ ( A1 ) and OsO₃(CH₂) ( B1 ). Significant differences have been found between the ethylene additions to the osmium compounds A1 and B1 and the rhenium compound C1 . Seven pathways for the reaction C1 +C₂H₄ were studied, but only the [2+2]_Re,C addition yielding rhenacyclobutane C5 is an exothermic process with a high activation barrier of 48.9 kcal mol^?1. The lowest activation energy (27.7 kcal mol^?1) is calculated for the [2+2]_Re,C addition, which leads to the isomeric form C5 ′. Two further concerted reactions [3+2]_O,C, [3+2]_O,O, and [2+2]_Re,O and the addition/hydrogen migration of ethylene to one oxo ligand are endothermic processes which have rather high activation barriers (>35 kcal mol^?1). Four isomerization processes of C1 have very large activation energies of >65 kcal mol^?1. The ethylene addition to the osmium compounds A1 and B1 are much more exothermic and have lower activation barriers than the C₂H₄ addition to C1 . Copyright © 2007 John Wiley & Sons, Ltd.     

Calculation of anharmonic effects in the unimolecular dissociation of M2+(H2O)2 (M = Be,Mg, and Ca)

The anharmonic and harmonic rate constants of the unimolecular dissociation of M²⁺(H₂O)₂ (M = Be, Mg, and Ca) were calculated using the Rice–Ramsperger–Kassel–Marcus theory. The anharmonic effects of the reactions were investigated. The results show that the energy barrier of the dissociation of Be²⁺(H₂O)₂ is 68.47 kcal/mol, and the anharmonic (T_4000K = 4.28×10⁸ s^?1) and harmonic (T_4000K = 4.22×10⁸ s^?1) rate constants were close in value in both the canonical and microcanonical systems. The energy barriers of the two steps for the dissociation, Mg²⁺(H₂O)₂ → MgOH⁺+H₃O⁺, were 37.41 and 11.39 kcal/mol, and those for the dissociation, Ca²⁺(H₂O)₂ → CaOH⁺+H₃O⁺, were 21.15 and 26.42 kcal/mol. The anharmonic effect of the two reactions is significant and cannot be neglected in both the canonical and microcanonical systems. The comparison also shows that the rate constants of the dissociation of Ca²⁺(H₂O)₂ have the maximum values, while those of Be²⁺(H₂O)₂ have the minimum values in the three reactions; however, the anharmonic effect also shows the similar trend in the comparison.