Ionization energy of KOH and the dissociation energies of KOH and KOH+

High level ab initio, up to RCCSD(T), and B3LYP calculations were employed to calculate thermochemical properties for KOH and KOH⁺. Basis sets were of both all-electron and effective core potential (ECP) types: in both cases large, flexible valence basis sets were used, and the largest basis sets were of quintuple-zeta quality. Both KOH and KOH⁺ were found to be linear; in the latter case, the Renner-Teller effect is discussed. The results are close to convergence with regard to both basis sets and levels of theory. The most reliable quantities are: first AIE(KOH) = 7.38±0.02eV; D₀(K- -OH) = 82 ± 1 kcal mol-¹; D₀(K+ …OH) = 11.4 ± 1 kcalmol-¹; δH_f(KOH,298K) = ?53 ± 1 kcalmol^?1; and δH_f(KOH⁺,298K) = 119 ± 1 kcalmol-¹.     

The enthalpy of formation of 2II CH

The standard enthalpy of formation, δ_fH^o, of² II CH has been determined at converged levels of ab initio electronic structure theory, including high order coupled cluster and full configuration interaction benchmarks. The atomic Gaussian basis sets employed include the (aug)-cc-p(C)VXZ family with X = 3, 4, 5 and 6. Extrapolations to the complete one-particle basis set and the full configuration interaction limits, where appropriate, have been performed to reduce remaining computational errors. Additional improvements in the enthalpy of formation of ²II CH were achieved by appending the valence-only treatment with core-valence correlation, relativistic effects including spin-orbit correlation, and the diagonal Born-Oppenheimer correction. The recommended values for δfH^o ₀ and δA_f H _o ²⁹⁸ of ²II CH are 592.48^+0.47 _?0.56 kJ mol_?1 and 595.93 ^+0.47 _?0.56 kJ mol_?1, respectively.     

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.     

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.     

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.     

Quantum chemical studies of the solvation of the hydroxide ion

Abstract

To understand and model the solvation of the hydroxide ion, OH(H₂O)^? _n clusters, n = 1?5, are studied using ab initio quantum chemical techniques, largely at the MP2 level of theory using a double zeta plus polarization functions basis extended by diffuse functions. Energies and vibrational frequencies, together with thermodynamic quantities such as enthalpies, entropies and Gibbs free energies, are computed. This permits comparison with experimental estimates of the successive thermodynamic changes associated with the reaction OH(H₂O)^? _n + H₂O → OH(H₂O)^? _n+1. The theoretical values are in good agreement with experiment. The free energy of hydration of OH^? is modelled by a composite discrete-continuum method where the effects of the first hydration shell (n = 3) are obtained from the gas phase cluster calculation, while the long-range effects are modelled using self consistent reaction field theory, namely by calculating the solvation energy of OH(H₂O)^? _n in a dielectric continuum. The best estimate of the solvation (free) energy at 298 K is ?84·5 kcal mol^?1, compared to the experimental value of ?102·8 kcal mol^?1.     

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).     

Explicitly correlated benchmark calculations on C8H8 isomer energy separations: how accurate are DFT,double-hybrid,and composite ab initio procedures?

Accurate isomerization energies are obtained for a set of 45 C₈H₈ isomers by means of the high-level, ab initio W1-F12 thermochemical protocol. The 45 isomers involve a range of hydrocarbon functional groups, including (linear and cyclic) polyacetylene, polyyne, and cumulene moieties, as well as aromatic, anti-aromatic, and highly-strained rings. Performance of a variety of DFT functionals for the isomerization energies is evaluated. This proves to be a challenging test: only six of the 56 tested functionals attain root mean square deviations (RMSDs) below 3?kcal?mol^?1 (the performance of MP2), namely: 2.9 (B972-D), 2.8 (PW6B95), 2.7 (B3PW91-D), 2.2 (PWPB95-D3), 2.1 (ωB97X-D), and 1.2 (DSD-PBEP86) kcal?mol^?1. Isomers involving highly-strained fused rings or long cumulenic chains provide a ‘torture test’ for most functionals. Finally, we evaluate the performance of composite procedures (e.g. G4, G4(MP2), CBS-QB3, and CBS-APNO), as well as that of standard ab initio procedures (e.g. MP2, SCS-MP2, MP4, CCSD, and SCS-CCSD). Both connected triples and post-MP4 singles and doubles are important for accurate results. SCS-MP2 actually outperforms MP4(SDQ) for this problem, while SCS-MP3 yields similar performance as CCSD and slightly bests MP4. All the tested empirical composite procedures show excellent performance with RMSDs below 1?kcal?mol^?1.     

Characterization of the [Xtilde] 1Σ+ à 3Π and à 1Π electronic states of BBO

The three lowest-lying electronic states, [Xtilde] ¹Σ⁺, à ³II and à ¹II, of the linear BBO molecule have been systematically investigated using ab initio electronic structure theory. The equilibrium structures and physical properties including dipole moments, vibrational frequencies and associated infrared intensities, Renner parameters and energetics for the three states of BBO have been determined employing SCF, CISD, CCSD and CCSD(T) levels of theory and a wide range of basis sets. The ground state of BBO presents a degenerate real bending frequency, while the à ³II and à ¹II states show two distinct real bending frequencies due to the Renner-Teller interaction. The bending motion of the à ¹II state was analysed using the equation-of-motion (EOM)-CCSD and EOM-CC3 techniques in order to avoid possible variational collapse to a lower-lying state. The [Xtilde] ¹Σ⁺-à ³II separation was predicted to be T ₀ = 16.6 kcal mol^?1 (5800 cm^?1, 0.719 eV) at the cc-pVQZ CCSD(T) level of theory. With the cc-pVQZ EOM-CC3 method the [Xtilde] ¹Σ⁺-à ¹II splitting was predicted to be T ₀ = 48.0 kcal mol^?1 (16 800 cm^?1, 2.08 eV), which is in good agreement with the experimental value of T ⁰ = 46.6 kcal mol^?1 (16 300 cm^?1, 2.02 eV). The Renner parameters and averaged harmonic frequencies of the bending mode were determined to be ? = 0.184 and ω₂ = 363 cm^?1 for the à ³II state, and ? = 0.246 and ω₂ = 383cm^?1 for the à ¹II state. The theoretical [Xtilde] ¹Σ⁺ state harmonic B-B stretching frequency ω₃ = 636 cm^?1 is somewhat higher than the experimental estimate of 582 cm^?1 and the predicted à ¹II state harmonic B-B stretching frequency ω₃ = 861 cm^?1 is significantly higher than the experimental estimate of 440 cm^?1     

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.     

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.     

Theoretical study of the CH3 + NS and related reactions: mechanism of HCN formation

More than thirty equilibrium and transition structures on the [CH₃NS] potential energy surface have been located using the B3LYP/6-311 + + G(d,p) method. Thioformaldoxime (3) turns out to be the most stable isomer followed by thionitrone (2), a three-membered ring, thionitrosyl methane (1) and thiazyl methane. These isomers are connected to each other by 1,2H and 1,3H shifts and ring—chain rearrangement, but the associated energy barriers are rather high, making most of them stable with respect to unimolecular transformations. Starting from CH₃ + NS, a possible initial atmospheric reaction, HCN formation appears to be the most favoured process through a cascade involvement of 1, 2 and 3. The standard heats of formation, ΔH ⁰ _f,298 are calculated to be: 3, 149kJmol^?1; and 1, 218kJmol^?1; using the CCSD(T)/6-311+ +G(3df,2p) method, with an error of ±10kJ mol^?1.     

Gas‐phase acidity of bis[(perfluoroalkyl)sulfonyl]imides. Effects of the perfluoroalkyl group on the acidity

The gas‐phase acidity (GA) values were determined for a number of perfluoroalkyl‐substituted sulfonylimides by measuring proton‐transfer equilibria using a Fourier transform ion cyclotron resonance (FT‐ICR) mass spectrometer. The GA scale below 286.5 kcal mol^?1 for (CF₃SO₂)₂NH was extended and partially revised. The GA value of (C₄F₉SO₂)₂NH which is currently the strongest acid was revised from 284.1 to 278.6 kcal mol^?1. The effect of fluorine atoms on the acidity of perfluoroalkyl‐substituted sulfonylimides was described with the following model where N_(α), N_(β), N_(γ), and N_(δ) are the numbers of fluorine atoms at α, β, γ, and δ position in R_fSO₂ (R_f = perfluoroalkyl group), respectively. This correlation indicates that the electron‐withdrawing ability of the R_fSO₂ group can be described in terms of the number of fluorine atoms in the perfluoroalkyl group corrected by taking into account their positions. Copyright © 2014 John Wiley & Sons, Ltd.     

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN+ and CNN−: a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN⁺ and CNN^? were investigated at a high level of theory. Very accurate structural parameters for the states X ³Σ^? and A ³Π of CNN, and X ²Π of both CNN⁺ and CNN^? were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Π states. At the UCCSD(T)/CBS_T-5 level of theory, Δ_f H(0 K) = 138.89 kcal/mol and Δ_f H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS_T-5 approach for the states X ³Σ^?, A ³Π, and a ¹Δ of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b ¹Σ⁺, c ¹Π, d ¹Σ^?, and B ³Σ^?, and excitation energies (T_e) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X ³Σ^? ← E ³Π.     

The concentration dependence of the proton chemical shift and the deuterium quadrupole coupling parameter for binary solutions of ethanol

Concentration dependent experimental measurements of the ethanol hydroxyl proton chemical shift σ_H for binary solutions were carried out. The solvents used were carbon tetrachloride (CCl₄), benzene, chloroform, acetonitrile, acetone and dimethylsulphoxide (DMSO). The chemical shift values range from 0.69 ppm (relative to TMS) for dilute ethanol (extrapolated to infinite dilution) in CCl₄ to 5.34 ppm for neat liquid ethanol. Ab initio calculations of the ethanol-solvent hydrogen bond energies show a correlation with the values for the chemical shift. The hydrogen bond energies for ethanol-solvent dimers range from 0.63 kcal mol^?1 for ethanol-CCl₄ to 9.34 kcal mol^?1 for ethanol-DMSO. Theoretical calculations show a linear correlation between the deuterium quadrupole coupling parameter XD ^ar d the isotropic proton chemical shift σ_H: X_D(kHz) = 291.48 ? 14.96 σ_H, where σ_H is the proton chemical shift in ppm relative to TMS (R ² = 0.99). Using the concentration dependent chemical shift data and this equation, X_D ia observed to range from 280 kHz for very dilute concentrations in CCl₄, where the primary species is ethanol monomer, to 210 kHz for the neat liquid that is comprised primarily of cyclic pentamers.     

The 1,2-hydrogen shift reaction for monohalogenophosphanes PH2X and HPX (X = F,Cl)

The aim of the present study was to perform a quantum chemical investigation in the 1,2-hydrogen shift reaction for the PH₂X and HPX molecules (X = F,Cl). Several phosphorus–halogen-bearing molecules were studied, including PH₂F, PH₂Cl, HPF, HPCl, HPFH, HPClH, PFH and PClH. The energies of stationary and saddle points on the ground electronic potential energy surface were investigated with post-Hartree–Fock methods [CCSD(T), MP2, QCISD] and different DFT functionals. The PH₂F 1,2-hydrogen shift energy barrier was 75 kcal mol^?1 at the CCSD(T) level and only a small increase in this value was observed for the HPF isomerisation. In contrast, the HPCl 1,2-hydrogen shift barrier is higher than the PH₂Cl one, which presented a barrier height of 69 kcal mol^?1 among CCSD(T) and composite methods. The rate constants of these unimolecular rearrangements varied from 10^?44 to 10^?38 s^?1, and these isomerisation channels exhibited large half-lives. In addition, the heat of formation of each monohalogenophosphane was also calculated. The Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbital (NBO) analysis were also employed to characterise the differences between the phosphorous–halogen bonds.     

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.     

High pressure 17O NMR study of solvent exchange on square planar complexes: Water exchange on [Pd(dien)H2O]2+

Abstract

The water exchange reaction of [Pd(dien)H₂O]²⁺ (dien = diethylenetriamine) was studied as function of temperature (268-308 K) and pressure 0.1-197 MPa) using ₁₇O NMR techniques. The rate and activation parameters are: k_cx = 5100 s^?1 at 298 K; ΔH^# =38 kJ mol^?1; ΔS^# = -47 JK^?1 mol^?1; ΔV^# = -2.8 cm³ mol^?1 at 296 K. The results are discussed in reference to solvent exchange data for other Pt(II) and Pd(II) complexes, and are interpreted in terms of an associatively activated substitution process.