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 electron affinity of Al13H cluster: high level ab initio study

Al₁₃H clusters have been considered candidates for cluster assembled materials. Here we have carried out benchmark calculations for the Al₁₃H cluster, both neutral and anionic, with the aim of verifying the nature of stationary points on the potential energy surface, studying dynamics of H atom and determining an adiabatic electron affinity. A range of correlated methods applied include second-order perturbation theory (MP2), spin-component-scaled MP2, coupled electron pair (CEPA) and coupled cluster singles and doubles with perturbative triple corrections (CCSD(T)). These methods are used in combination with the correlation consistent basis sets through aug-cc-pVTZ including extrapolation to the complete basis set (CBS) limit. Performance of several different flavours of density functional theory (DFT) such as generalised gradient approximation (GGA), hybrid GGA, meta-GGA and hybrid-meta-GGA is assessed with respect to the ab initio correlated reference data. The harmonic force constant analysis is systematically performed with the MP2 and DFT methods. The MP2 results show that for neutral Al₁₃H only the hollow structure is a potential energy minimum, with the bridged structure being a transition state for the H shift from the hollow site to the adjacent hollow site. The CCSD(T)/aug-cc-pVTZ (CCSD(T)/CBS) estimate of the energy barrier to this H shift is 2.6 (2.9) kcal/mol, implying that the H atom movement over the Al₁₃H cluster surface is facile. By contrast, the DFT force constant analysis results suggest additional terminal and bridged minima structures. For the anion Al₁₃H^?, exhibiting ‘stiffer’ potential energy surface compared to the neutral, the existence of the hollow and terminal isomers is consistent with the earlier photoelectron spectroscopy assignment. The adiabatic electron affinity of Al₁₃H is determined to be 2.00 and 1.95 eV (the latter including the ΔZPE correction) based on the CCSD(T) energies extrapolated to the CBS limit, whereas the respective CCSD(T)/CBS thermodynamic EA values are 2.79 and 2.80 eV.     

A theoretical study on the kinetics of multichannel Multiwell reaction of H2S(1A1) with HO2(2A′′)

The mechanism and kinetics of the reaction of hydrogen sulphide (H₂S(¹A₁)) with hydroperoxyl radical (HO₂(²A″)) on the lowest doublet potential energy surface have been theoretically studied. The potential energy surface for possible pathways has been investigated by employing Complete Basis Set (CBS), DFT, and CCSD(T) methods. Three possible pathways are suggested for the title reaction. The most probable entrance channel consists of formation of a hydrogen-bonded pre-reaction complex (vdw1) and two energised intermediates. Multichannel RRKM-Steady State Approximation and CVT calculations have been carried out to compute the rate constants over a broad range of temperature from 200?K to 3000?K to cover the atmospheric and combustion conditions and pressure from 0.1 to 2000?Torr. No sign of pressure dependence was observed for the title reaction over the stated range of pressure. We have shown that the major products of the title reaction are H₂O₂ and SH while at higher temperatures, formation of the other products such as H₂O, HOS, HSOH and OH are feasible, too. Our calculated overall rate constant is in agreement with the reported experimental data in the literature.     

Physisorption and dissociative chemisorption of H2 on sub-nanosized Al$_{13}^{-}$ anion cluster: ab initio study

We have studied the interaction of Al₁₃^-_{13}^{-} anion cluster with H₂. Both the long range interaction and dissociative adsorption have been examined using the established correlated ab initio methods, MP2 and CCSD(T), in conjunction with the augmented correlation consistent basis sets up to aug-cc-pVTZ. The formation of the weakly bound (physisorbed) end-on anion complex Al₁₃^-_{13}^{-}...H₂ is predicted for the interacting Al...H distances of 3.95 ? with the H-H axis pointing towards the ‘hollow’ site of Al₁₃^-_{13}^{-} and binding energy (D_e)D_{e}) of 0.7 kcal/mol at the estimated complete basis set (CBS) limit of CCSD(T). The barrier height for H₂ dissociation on Al₁₃^-_{13}^{-} of 41.6 (42.9) kcal/mol calculated at the ZPVE-corrected CCSD(T)/aug-cc-pVTZ (estimated CCSD(T)/CBS) level is at least twice as large as that evaluated by us for a dissociative adsorption of H₂ on an open-shell Al₁₃ neutral cluster. To our knowledge, this report presents the first “benchmark” quality study of the physisorption and dissociative chemisorption of molecular hydrogen on Al₁₃^-_{13}^{-} anion cluster.     

Weakly interacting molecular clusters of CO with H2O,SO2, and NO+

Intermolecular interactions in three dimers, CO···H₂O, CO···SO₂, and CO···NO⁺, were studied at the CCSD(T) level of theory, using a series of the augmented correlation consistent polarised basis sets. Interaction energy and its components as well as vibrational spectra for local minima were computed using both harmonic and anharmonic approximations. While CO···H₂O and CO···SO₂ are weakly bound with the binding energies ?7.4 and ?6.4 kJ/mol, CO···NO⁺ is much more stable with the binding energy of ?32.8 kJ/mol corresponding to ΔG = ?4.7 kJ/mol at 254 K.     

Isomers of OCS and their reaction with H2O on singlet potential energy surface

The isomers of the carbonyl sulfide (OCS) molecule are investigated in detail at CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level of theory. One cyclic isomer was identified along with three different linear minima of the OCS molecule. Three interconversion transition states were also located between cyclic and linear forms of OCS. Among these four isomers, the singlet potential energy surface (PES) for the molecule–molecule reaction between the three most energetically favoured isomers of OCS and H₂O has been explored theoretically at the CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level. This singlet PES comprises of three paths. Path 1 is the reaction of linear OCS molecule with water producing the major product P1 (CO₂?+?H₂S), minor product P2 (S?+?HCOOH) and two isomers via 14 minima and 15 transition states. The Path 2 is an isomerization process in which cyclic isomer of OCS reacts with water molecule via another initial barrierless aduct producing five isomers of the OCS–H₂O system through five interconversion transition states. The reaction of linear COS isomer with water is shown in Path 3. This path produces the radicals SH and COOH from another COS–H₂O complex via a transition state. Among these three products, the product P1 is energetically most favoured. The overall exothermicity of the product channels for the formation of major product P1 on PES is calculated to be about 10.60?kcal/mol possessing initial high entrance barriers of 45.48 and 55.47?kcal/mol in two possible pathways. As the process is favoured thermodynamically but not kinetically, the reaction is expected to be very slow.     

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.     

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.     

New trends along hydrogen polyoxides: unusually long oxygen–oxygen bonds in H2O6 and H2O7

By means of coupled cluster, G4 and B3LYP calculations, we characterised polyoxides H₂O₆ and H₂O₇. These two molecules behave very differently from lower polyoxides, given that some isomers present unusual bonding. In the case of H₂O₆, the central OO bond is predicted, to be 1.909 Å, at the CCSD(T)/cc-pVTZ level. Two conformational isomers of H₂O₆ display nearly the same enthalpy of formation, but only one of them has extremely long OO bonds. For H₂O₇, the most stable isomer also has two unusually long OO bonds. At the CCSD(T) level, and after extrapolation to the complete basis set limit, this isomer is predicted to be 5.37 kcal mol^–1 more stable than the one with short OO bonds, and the longest OO bond distance is expected to be close to 1.96 Å. Analysis of correlation energies indicated that the new isomers found for H₂O₆ and H₂O₇ are among the most strongly correlated molecules that can be formed with first-row atoms.     

A coupled cluster study of the magnetisability,rotational g-tensor and quadrupole moment of NF3, PF3 and AsF3

ABSTRACT

A coupled-cluster investigation of magnetic and electric properties of NF₃, PF₃ and AsF₃ provides for a comparison with recent experimental data. For PF₃, achieving reliable values for the magnetisability and rotational g-tensor of PF₃ has been particularly challenging. We report the most accurate calculations to date for PF₃; for the vibrationally corrected anisotropic magnetisability, our extrapolated CCSD(T)/CBS value of ?0.290 a.u is within the uncertainty limits of the most recent experimental value of ?0.286 ± 0.042 a.u. For the rotational g-tensor of PF₃, agreement between theory and experiment for the g_⊥ component is excellent (deviation of less than 0.0006 a.u.). However, the g_|| component remains problematic even though our vibrationally corrected CCSD(T)/CBS value of ?0.0387 a.u is in closer agreement with the recently revised experimental value of ?0.0470 ± 0.0020 a.u. than the original value of ?0.0815 ± 0.0020 a.u. The origin of the remaining discrepancy remains unclear. Dipole and quadrupole moments have also been investigated.     

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.     

Ab initio anharmonic force field and equilibrium structure of the sulfonium ion

The semi-experimental equilibrium structure of the sulfonium ion, , has been obtained from the experimental ground-state rotational constants available for five isotopologues and the corresponding vibrational corrections computed at the CCSD(T)/cc-pwCVQZ level of theory. This geometry has been found in very good agreement with the pure ab initio equilibrium structure calculated at the CCSD(T) level of theory using a basis set of sextuple-zeta quality and including core correlation corrections. The anharmonic force field has been used for deriving spectroscopic properties: in particular, in addition to the vibrational corrections, the rotational parameters of the SH₂D⁺ isotopic species, not yet experimentally observed, have been predicted to a guessed good accuracy.     

高次激发相关能在构建HCN-He亚光谱精度势能面中的角色:通过实验高分辨率光谱进行严格检测

Interpreting high-resolution rovibrational spectra of weakly bound complexes commonly requires spectroscopic accuracy (<1 cm^-1) potential energy surfaces (PES). Constructing high-accuracy ab initio PES relies on the high-level electronic structure approaches and the accurate physical models to represent the potentials. The coupled cluster approaches including single and double excitations with a perturbational estimate of triple excitations (CCSD(T)) have been termed the "gold standard" of electronic structure theory, and widely used in generating intermolecular interaction energies for most van der Waals complexes. However, for HCN-He complex, the observed millimeter-wave spectroscopy with high-excited resonance states has not been assigned and interpreted even on the ab initio PES computed at CCSD(T) level of theory with the complete basis set (CBS) limit. In this work, an effective three-dimensional ab initio PES for HCN-He, which explicitly incorporates dependence on the Q₁ (C-H) normal-mode coordinate of the HCN monomer has been calculated at the CCSD(T)/CBS level. The post-CCSD(T) interaction energy has been examined and included in our PES. Analytic two-dimensional PESs are obtained by least-squares fitting vibrationally averaged interaction energies for v₁(C-H)=0, and 1 to the Morse/Long-Range potential function form with root-mean-square deviations (RMSD) smaller than 0.011 cm^-1. The role and significance of the post-CCSD(T) interaction energy contribution are clearly illustrated by comparison with the predicted rovibrational energy levels. With or without post-CCSD(T) corrections, the value of dissociation limit (D₀) is 8.919 or 9.403 cm^-1, respectively. The predicted millimeter-wave transitions and intensities from the PES with post-CCSD(T) excitation corrections are in good agreement with the available experimental data with RMS discrepancy of 0.072 cm^-1. Moreover, the infrared spectrum for HCN-He complex is predicted for the first time. These results will serve as a good starting point and provide reliable guidance for future infrared studies of HCN doped in (He)_n clusters.     

Coupled cluster study of NMR shielding constants and spin-rotation constants in SiH4, PH3 and H2S molecules

Ab initio methods are applied to analyse the NMR shielding constants and spin-rotation constants in SiH₄, PH₃ and H₂S molecules. The electron correlation effects are studied applying the MP2 and coupled cluster perturbation approaches. The basis set convergence is examined at the same time, and the final results for the equilibrium geometries are obtained at the CCSD(T)/cc-pCVQZ level. Zero-point vibrational and temperature contributions are computed at the SCF, MP2 and CCSD level of approximation. In addition, for the shielding constants we also estimate the relativistic effects, to determine total values of the shielding of the third-row nuclei in the studied molecules. Our final results for the shielding constants at 300?K are σ (²⁹Si in SiH₄)?=?482.35?ppm, σ (³¹P in PH₃)?=?611.64?ppm and σ (³³S in H₂S)?=?736.13?ppm. These values, together with estimated corrections and error bars, can be used to determine absolute NMR shielding scales for the heavy nuclei.     

The catalytic effects of H2CO3, CH3COOH,HCOOH and H2O on the addition reaction of CH2OO + H2O → CH2(OH)OOH

The addition reaction of CH₂OO + H₂O → CH₂(OH)OOH without and with X (X = H₂CO₃, CH₃COOH and HCOOH) and H₂O was studied at CCSD(T)/6-311+ G(3df,2dp)//B3LYP/6-311+G(2d,2p) level of theory. Our results show that X can catalyse CH₂OO + H₂O → CH₂(OH)OOH reaction both by increasing the number of rings, and by adding the size of the ring in which ring enlargement by COOH moiety of X inserting into CH₂OO···H₂O is favourable one. Water-assisted CH₂OO + H₂O → CH₂(OH)OOH can occur by H₂O moiety of (H₂O)₂ or the whole (H₂O)₂ forming cyclic structure with CH₂OO, where the latter form is more favourable. Because the concentration of H₂CO₃ is unknown, the influence of CH₃COOH, HCOOH and H₂O were calculated within 0–30 km altitude of the Earth's atmosphere. The results calculated within 0–5 km altitude show that H₂O and HCOOH have obvious effect on enhancing the rate with the enhancement factors are, respectively, 62.47%–77.26% and 0.04%–1.76%. Within 5–30 km altitude, HCOOH has obvious effect on enhancing the title rate with the enhancement factor of 2.69%–98.28%. However, compared with the reaction of CH₂OO + HCOOH, the rate of CH₂OO···H₂O + HCOOH is much slower.     

Anharmonic force fields of cis- and trans-S1 C2H2

We calculate second-order vibrational perturbation theory (VPT2) anharmonic force fields for the cis and trans conformers of S₁ C₂H₂, and compare the results to experiment. The vibrational assignments of recently observed levels belonging to the cis well are of particular interest. A refined estimate of the cis origin position (44,870?±?10?cm^?1) is proposed, and preliminary low-energy fits to the global J?=?K?=?0 trans level structure are also described. The performance of perturbation theory in this isomerizing system is examined, and both surprising successes and failures are encountered. We examine these and their causes, and offer practical suggestions for avoiding the pitfalls of applying perturbation theory to systems with large amplitude motions.     

Semiexperimental equilibrium structure of 1-methylisatin from gas-phase electron diffraction data and structural changes in isatin due to 1-methyl and 5-fluoro substituents as predicted by coupled cluster computations

ABSTRACT

The molecular structure of 1-methylisatin (1) has been studied by gas-phase electron diffraction (GED) and quantum chemical computation up to the coupled cluster (CCSD(T)) level of theory. The semiexperimental equilibrium structure (C_s point group symmetry) has been determined from the GED data taking into account anharmonic vibrational corrections calculated in curvilinear coordinates with the B2PLYP force field. To observe fine structural effects due to the presence of different kinds of substituents, the accurate molecular structures of isatin (2) and 5-fluoroisatin (3) were also computed at the CCSD(T) level. The (O=)C–C(=O) carbon–carbon bonds of the pyrroline moiety in 1–3 are found to be remarkably longer than the typical single C–C bond. The electron donating methyl group causes a decrease of the C?N?C angle and an elongation of the N?C bond lengths in the pyrroline ring by 0.7° and up to 0.008?Å, respectively, whereas the electron withdrawing fluorine atom increases the ipso CCC angle by 2.5° in comparison to that in unsubstituted isatin.     

Additivity and non-additivity of dissociation energies in intermolecular interactions. Theoretical studies on (H2)n,n = 2-8, (CO2)n,n = 2-6 and (HF)n,n = 2-8

CCSD(T) and MP2 results using the aug-cc-pV5Z basis set are reported for chain, cyclic and other structures of the clusters (H₂)_n, n?=?2-8, (CO₂)_n, n?=?2-6 and (HF)_n, n?=?2-8. In chain-like structures of (H₂)_n and (CO₂)_n, with the bonding type of the dimer maintained, the dissociation energy D_e of the dimer doubles for the trimer, triples for the tetramer, and so on. Due to these systems being dominated by short-range forces, interactions are essentially restricted to neighbouring monomers. For other types of (H₂)_n and (CO₂)_n structures, the multipliers relative to the dimerisation energy can be much higher. Dissociation energies for the hexamers in S₆ symmetry of both H₂ (379?cm^?1) and CO₂ (4925?cm^?1) are over ten times the respective dimerisation energies. For the chain-like trimer of HF, however, D_e is in excess of double the dimer value. Mainly due to longer-range dipolar forces, the interactions reach beyond the neighbouring monomers. The interaction energy between HF monomers in chains follows an approximate R^?2 (R being the F–F distance) relationship, The calculated dissociation energies of the HF octamer are 15,985?cm^?1 (factor of 10.4) for the chain, and 21,003?cm^?1 (factor of 13.7) for the C_6h cyclic structure.     

用从头算方法研究GeCl2分子的非谐振力场和光谱常数

用从头算方法的MP2和CCSD（T）方法结合cc-pVTZ基组计算了二氯化锗同位素（⁷⁰GeCl₂、⁷²GeCl₂和⁷⁶GeCl₂）分子的平衡结构、光谱常数和非谐振力场.二氯化锗的几何结构、转动常数、振转相互作用常数、谐频、非谐振常数、四次和六次离心畸变常数、三次和四次力常数的计算结果与实验结果符合较好,二氯化锗分子的同位素效应较小,可能的原因是Ge同位素的质量变化相对较小.两种方法计算的结果均与实验结果符合,但CCSD（T）方法比MP2计算结果的偏差稍大一些,可能的原因是CCSD（T）方法在描述过共价Cl原子的电子相关时不够充分.