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.     

Thermochemistry of the reactions between CN+ and H2O in the gas phase

The [H₂, C, N, 0]⁺ potential energy surface (PES) has been explored by means of high-level ab initio calculations, carried out in the framework of the G2 theory. From this survey we concluded that the predominant products of the CN⁺ + H₂O reaction are the result of the dissociation of HNCOH⁺ species and to a much lesser extent of the CNHOH⁺ cation to yield CNH⁺ + OH. According to our results HCN⁺ should not be a product of this reaction because all pathways leading to its formation are unfavourable with regards to other competitive processes. Other reactive channels lead to the formation of the H₂ONC⁺ structure which dissociates into CN + H₂O⁺. The loss of NH(3σ) and O(³P) seems to take place following spin-forbidden reaction paths through an intersystem crossing between the singlet and the triplet PESs. The global minimum of the PES, H₂NCO⁺ is easily accessible and should lead to the loss of carbon monoxide which has not been experimentally observed in CN⁺ + H₂O reactions. We cannot offer a clear explanation for this disagreement between theory and experiment.     

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.     

Structure and bonding of the (C6H6)+ 2 radical

To elucidate the relative stability of various structures of the benzene dimer cation radical, (C₆H₆)⁺ ₂ in its ground and low-lying excited states, ab initio complete active space self-consistent field (CASSCF), multi-reference singly and doubly excited configuration interaction (MRSDCI), and multi-reference coupled pair approximation (MRCPA) calculations were performed. Full optimization was performed at the CASSCF level for various structures of the dimer cation, followed by MRSDCI and MRCPA calculations. It was found that the global minimum of the cation is at a slipped C_2h sandwich structure but there are some other sandwich structures with almost the same stability, being within about kcal mol^?1. T-shape structures are less stable than the sandwich structures, by more than 5 kcal mol^?1 by MRCPA calculations. Low lying electronic excited states in various structures are also discussed.     

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.     

Probing spin-orbit quenching in Cl (2P) + H2 via crossed molecular beam scattering

In our previous work we investigated electronically non-adiabatic effects in using crossed molecular beam scattering coupled with velocity mapped ion imaging. The prior experiments placed limits on the cross-section for electronically non-adiabatic spin-orbit excitation and electronically non-adiabatic spin-orbit quenching . In the present work, we investigate electronically non-adiabatic spin-orbit quenching for which is the required first step for the reaction of Cl^* to produce ground state HCl+H products. In these experiments we collide Cl (²P) with H₂ at a series of fixed collision energies using a crossed molecular beam machine with velocity mapped ion imaging detection. Through an analysis of our ion images, we determine the fraction of electronically adiabatic scattering in Cl^* +H₂, which allows us to place limits on the cross-section for electronically non-adiabatic scattering or quenching. We determine the following quenching cross-sections σ _quench(2.1 kcal/mol) = 26 ± 21 ?², σ _quench(4.0 kcal/mol) = 21 ± 49 ?², and σ _quench(5.6 kcal/mol) = 14 ± 41 ?².     

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.     

Theoretical rovibrational spectroscopy beyond fc-CCSD(T): the cation CNC+

An accurate near-equilibrium potential energy surface (PES) for CNC⁺ is constructed based on a high-level composite ab initio method. By combining explicitly correlated all-electron CCSD(T)-F12b with scalar relativistic effects and higher order correlation up to coupled cluster theory with singles, doubles, triples and quadruples (CCSDTQ) we achieve convergence in the wavenumbers of the fundamentals to ca. 1 cm^?1. Rovibrational energies are calculated in a variational approach and vibrational term energies and rotational constants are in excellent agreement with available experimental data. Accurate values for centrifugal distortion constants of CNC⁺ in different vibrational states are predicted. Especially the centrifugal distortion constants in the vibrational ground state of D₀ = 0.563 · 10^?6 cm^?1 and H₀ = 0.188 · 10^?10 cm^?1 should be superior to experimentally derived values. Reassignments of some experimentally observed transitions are suggested based on a comparison of experimental and calculated term differences. The bending part of the PES appears to be almost quartic and the band origin of the bending vibration is predicted at 94.2 cm^?1. Absolute line intensities are calculated for various transitions in CNC⁺. For the bending vibration, an intensity is predicted that is three orders of magnitude smaller than for the antisymmetric stretching vibration.     

Comparison of the Kinetics of the Exchange Mechanism of the Thallium Ion with 18C6 and with Pentaglyme in Acetonitrile Solutions

Abstract

In acetonitrile solutions, the exchange reaction is bimolecular in the Tl⁺ + 18C6 system, while in the Tl⁺ + pentaglyme system the associative-dissociative and the bimolecular mechanisms coexist at room temperature and the bimolecular exchange reaction dominates at 263° K. For the bimolecular mechanism in the case of Tl⁺ + 18C6 and the associative-dissociative mechanism in the case of Tl⁺ + pentaglyme, the activation energies of the exchange reactions change with temperature. At 298° K, in the Tl⁺ + 18C6 system the activation energy for the bimolecular exchange reaction is ≈ 2 kcal.mol^?1 and exchange rate constant (k₁) is (4.1 ± 0.1) × 10⁷ s^?1mol^?1; in the Tl⁺ + pentaglyme system, the activation energy for the associative-dissociative exchange reaction is ≈ 5 kcal mol^?1 and the decomplexation rate constant (k_?2) is (2.2 ± 0.4) X 10⁵ s^?1. The activation energy for the bimolecular exchange in the Tl⁺ + pentaglyme system was determined to be 3.00 ± 0.05 kcal.mol_?1 and the exchange rate constant (3.0 ± 0.1) X 10⁸ s^?1 mol^?1.     

Synthesis and Characterization of 2,2-Biimidazole Complexes of Oxocations of Molybdenum (VI,V) and Uranium(VI)

Abstract

2,2′-Biimidazole complexes of MoO₂ ⁺², MoO₂ ⁺ and UO₂ ⁺² have been prepared and characterized by elemental analysis, conductance; and ¹H NMR, IR and electronic spectra. Two types of complexes have been identified. Those obtained from slightly acidic solutions have the formulae MoO₂ (H₂bim)Cl₂.2H₂O 1, UO₂(H₂bim) (Ac)₂ 2 and UO₂(H₂bim)Cl₂.2H₂O 3; whereas those from alkaline solutions have the formulae Mo₂O₄(Hbim)₂.2H₂O 4, and MO₂(Hbim)₂ (M = Mo(VI) 5, U(VI) 6). The infrared spectra of these complexes show characteristic biimidazole frequencies in the 3200–2500, 1550–1000 and 750 cm^?1 regions as well as metal oxygen double bonds in the 900 cm^?1 region. The stoichiometries of the acetate complex has been confirmed from ¹H NMR signal ratios of bimidazole to acetate protons at 7.3 and 2.3 ppm, respectively. The electronic spectrum of molybdenum(V) complex showed d-d transition band at ?13,500 cm^?1 in accord with that reported for copper (d⁹) imidazole complexes; as well as peaks due to charge transfer bands at 30,000–26,000 cm^?1 Peaks assignable to BIM → U(VI) were located at ?26,600 cm^?1. The most probable structures of these complexes have been suggested.     

A global potential energy surface for H2S+(X 4A′′) and quasi-classical trajectory study of the S+(4S) + H2(X1Σ+g) reaction

A novel global potential energy surface for H₂S⁺(X?⁴A″) based on accurate ab initio calculations is presented. Energies are calculated at the multi-reference configuration interaction level with Davidson correction using aug-cc-pVQZ basis set plus core-polarisation high-exponent d functions. A grid of 4552 points is used for the least-square fitting procedure in the frame of a many-body expansion. The topographical features of the new potential energy surface are here discussed in detail. Such a surface is then employed for dynamic studies of the S(⁴S) + H₂(X?¹Σ⁺_g) →SH⁺(X?³Σ^?) + H(²S) reaction using the quasi-classical trajectory method. State specific trajectories are calculated, for both ground and ro-vibrationally excited initial states of H₂(X?¹Σ⁺_g). Corrections to the zero point energy leakage of the classical calculations are also presented. Calculated reaction cross sections and rate constants are here reported and compared with available literature.     

Investigation of the kinetics of OH∗ and CH∗ chemiluminescence in hydrocarbon oxidation behind reflected shock waves

The temporal variation of chemiluminescence emission from OH^?(A² Σ ⁺) and CH^?(A² Δ) in reacting Ar-diluted H₂/O₂/CH₄, C₂H₂/O₂ and C₂H₂/N₂O mixtures was studied in a shock tube for a wide temperature range at atmospheric pressures and various equivalence ratios. Time-resolved emission measurements were used to evaluate the relative importance of different reaction pathways. The main formation channel for OH^? in hydrocarbon combustion was studied with CH₄ as benchmark fuel. Three reaction pathways leading to CH^? were studied with C₂H₂ as fuel. Based on well-validated ground-state chemistry models from literature, sub-mechanisms for OH^? and CH^? were developed. For the main OH^?-forming reaction CH+O₂=OH^?+CO, a rate coefficient of k ₂=(8.0±2.6)×10¹⁰ cm³?mol^?1?s^?1 was determined. For CH^? formation, best agreement was achieved when incorporating reactions C₂+OH=CH^?+CO (k ₅=2.0×10¹⁴ cm³?mol^?1?s^?1) and C₂H+O=CH^?+CO (k ₆=3.6×10¹²exp(?10.9 kJ?mol^?1/RT) cm³?mol^?1?s^?1) and neglecting the C₂H+O₂=CH^?+CO₂ reaction.     

Theoretical study of mechanisms for NCO with CH3 reaction

A detailed computational study has been performed at the QCISD(T)/6-311++G(d,p)//B3LYP/6-311++G(d,p) level for the NCO with CH₃ reaction by constructing singlet and triplet potential energy surfaces (PES). The results show that the title reaction is more favorable for the singlet PES than the triplet PES. On the singlet PES, the dominant channel is the barrierless addition of the O or N atom to the C atom of the methyl group to form CH₃NCO (IM1) and CH₃OCN (IM2). On the triplet PES, the favorable channel is the barrierless addition of the N atom to the C atom of the methyl group to form an intermediate CH₃NCO (³IM2), which then undergoes a N–C bond scission process to give out CH₃N + CO.     

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.     

Separation of overlapping singlet and triplet band systems in Cl2CS by LFES and MRS techniques

The visible absorption band systems of Cl₂CS have been reexamined using laser fluorescence excitation and magnetic rotation spectroscopy. Since the former is sensitive only to S₁ ← S₀ transitions and the latter primarily to T₁ ← S₀ transition, the spin multiplicities of the features in the overlapping singlet and triplet systems could be unambiguously determined. Analyses of the spectra gave new values for ν′₆ = 189 cm^?1 in the $\text{A}\text{?}$ state and ν′₁ = 923 cm^?1 in the $\text{a}\text{?}$ state. Bands with types A and C polarization were found to occur only very weakly in the $\text{A}\text{?}\text{←}\text{X}\text{?}$ spectrum in marked contrast to the corresponding system in the prototype thiocarbonyl compound, H₂CS.     

Addition complex and insertion isomers on the potential energy surface of the reaction of indium dimer with water studied with relativistic ECP

Stationary points on the lowest singlet and triplet In₂ + H₂O potential energy surfaces (PESs) have been explored using the coupled cluster method, including single and double excitations with perturbative triples (CCSD(T)), and the density functional theory (DFT), employing the effective core potential (ECP) for indium (In), which accounts for scalar relativistic effects, with the triple-zeta quality basis set. The CCSD(T) calculated binding energy and anharmonic ν₂-bending mode frequency for the triplet ground-state addition complex, In₂^… OH₂(³B₁), are consistent with the complex detected in the matrix isolation infrared (IR) spectroscopic study under the thermal conditions. The two minimum energy crossing points between the triplet and the singlet PESs that have been located between the structures of In₂^…OH₂ and the transition state for the O–H bond breakage are not likely to be thermally accessible under the low-temperature matrix conditions. With the CCSD(T)-calculated In₂ + H₂O reaction profile and anharmonic vibrational frequencies for several In₂(H)(OH) insertion product isomers, we support the IR matrix isolation detection (by two experimental groups) of the lowest energy singlet double-bridged In(μ-H)(μ-OH)In isomer. For the proposed two-step mechanism of H₂ elimination from the In₂(H)(OH) species, the estimated energy barriers are also compatible with experiment.     

Computational studies of the halooxyhalocarbenes XOCX (X = F,Cl)

Ab initio calculations at the B3LYP, MP2, MP4 and CCSD(T) levels of theory were performed to predict the stability of the halooxyhalocarbenes, XOCX (X = F, Cl). The calculations indicate that the nonlinear FOCF molecule is stable with an energy 16 kJ mol^?1 below the energy of possible reacting fragments F₂ and CO. However, a nonlinear equilibrium structure for ClOCCl was located, but it was found to be about 192 kJ mol^?1 higher in energy than the energy of Cl₂ and CO. The charge distribution in these molecules was analysed using the atoms in molecules (AIM) method.     

P‐Heterocyclic silylenes: a survey of stability with density functional theory

The effects of phosphorous atom on the stability, multiplicity, and reactivity of six‐member cyclic silylenes are investigated at B3LYP/AUG‐cc‐pVTZ//B3LYP/6‐31+G* and MP2/6‐311++G**//B3LYP/6‐31+G* coupled with appropriate isodesmic reactions. From a thermodynamic point of view, 1H‐2‐silaphosphinine‐2‐ylidene ( 1_a ) and 1H‐4‐silaphosphinine‐4‐ylidene ( 2_a ) are relatively the most stable with singlet–triplet energy gaps (ΔE_S–T) of 37.0 and 28.1 kcal/mol, respectively. The calculated energy barrier for the 1,2‐H shift of 1_a to the corresponding 2‐silapyridine ( 1 ) is 26.5 kcal/mol, which is lower than the 28.8 kcal/mol required for the 1,4‐H shift of 2_a to the corresponding 4‐silapyridine ( 2 ). In contrast to the previous reports, isodesmic reactions indicate that π‐donor/σ‐donor phosphorous destabilizes the singlet while stabilizes the triplet state. Both 1_a and 2_a silylenes appear invulnerable to the head‐to‐head as well as the head‐to‐tail dimerization, inviting experimental explorations. Copyright © 2011 John Wiley & Sons, Ltd.