Improving on equilibrium isotropic nuclear shielding constants

Isotropic nuclear shielding constants at the equilibrium molecular structure σ_eq and zero‐point vibrational corrections (ZPVCs) to σ_eq are evaluated using the B3LYP/aug‐cc‐pVTZ level of theory, as well as the KT2/aug‐cc‐pVTZ level of theory. Various scaling factors and systematic corrections are obtained by linear regression to experimental shielding constants. Comparisons of the scaled and systematically corrected equilibrium and vibrationally averaged shielding constants reveal that, at the 99% confidence level, the ZPVCs via second‐order perturbation theory do not improve the agreement of B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ calculated shielding constants with experiment. This holds true when the same analysis is applied to CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq of Teale et al. [Journal of Chemical Physics 2013, 138, 024111]. In addition, at the 99% confidence level, B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ scaled and systematically corrected shielding constants are found to be statistically no different from CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq. The use of scaling factors and systematic corrections could thus provide a cheaper but yet reasonably accurate alternative for the study of nuclear shielding constants of larger systems.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Nonempirical calculations of the one‐bond 29Si–13C spin–spin coupling constants taking into account relativistic and solvent corrections

The computational study of the one‐bond ²⁹Si–¹³C spin–spin coupling constants has been performed at the second‐order polarization propagator approximation (SOPPA) level in the series of 60 diverse silanes with a special focus on the main factors affecting the accuracy of the calculation including the level of theory, the quality of the basis set, and the contribution of solvent and relativistic effects. Among three SOPPA‐based methods, SOPPA(MP2), SOPPA(CC2), and SOPPA(CCSD), the best result was achieved with SOPPA(CCSD) when used in combination with Sauer's basis set aug‐cc‐pVTZ‐J characterized by the mean absolute error of calculated coupling constants against the experiment of ca 2 Hz in the range of ca 200 Hz. The SOPPA(CCSD)/aug‐cc‐pVTZ‐J method is recommended as the most accurate and effective computational scheme for the calculation of ¹J(Si,C). The slightly less accurate but essentially more economical SOPPA(MP2)/aug‐cc‐pVTZ‐J and/or SOPPA(CC2)/aug‐cc‐pVTZ‐J methods are recommended for larger molecular systems. It was shown that solvent and relativistic corrections do not play a major role in the computation of the total values of ¹J(Si,C); however, taking them into account noticeably improves agreement with the experiment. The rovibrational corrections are estimated to be of about 1 Hz or 1–1.5% of the total value of ¹J(Si,C). Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of the nitrile and methyl torsional vibrations of n‐propyl cyanide in its S0 electronic state

This work presents a structural and vibrational theoretical study of n‐propyl cyanide as a function of the nitrile and methyl torsional modes. A potential energy hypersurface is built at the MP4(SDQ)/aug‐cc‐pVTZ//MP2/aug‐cc‐pVTZ theory level. The equilibrium structure is found in a gauche conformation. Another minimum is found for the trans form. The maximum appears in a cis conformation. For the first time, the interconversion barriers between the different forms are calculated. A two‐dimensional anharmonic vibrational Hamiltonian is built for the nitrile and methyl torsional modes. We find the vibrational energy levels to organize in two stacks associated to the gauche and trans forms. Fundamental frequencies of 113.12 and 220.54 cm^?1 are predicted for the nitrile and methyl torsions in the equilibrium, gauche, conformer. In addition, we find symmetry allowed transitions between the gauche and trans energy levels stacks. The lowest transition is predicted to appear at 24.49 cm^?1. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational thermochemistry of glycolaldehyde

We use a variant of the focal point analysis to refine estimates of the relative energies of the four low‐energy torsional conformers of glycolaldehyde. The most stable form is the cis‐cis structure which enjoys a degree of H‐bonding from hydroxyl H to carbonyl O; here dihedral angles τ₁ (O?C? C? O) and τ₂ (C? C? O? H) both are zero. We optimized structures in both CCSD(T)/aug‐cc‐pVDZ and aug‐cc‐pVTZ; the structures agree within 0.01 Å for bond lengths and 1.0 degrees for valence angles, but the larger basis brings the rotational constants closer to experimental values. According to our extrapolation of CCSD(T) energies evaluated in basis sets ranging to aug‐cc‐pVQZ the trans‐trans form (180°, 180°) has a relative energy of 12.6 kJ/mol. The trans‐gauche conformer (160°, ±75°) is situated at 13.9 kJ/mol and the cis‐trans form (0°, 180°) at 18.9 kJ/mol. Values are corrected for zero point vibrational energy by MP2/aug‐cc‐pVTZ frequencies. Modeling the vibrational spectra is best accomplished by MP2/aug‐cc‐pVTZ with anharmonic corrections. We compute the Watsonian parameters that define the theoretical vibrational‐rotational spectra for the four stable conformers, to assist the search for these species in the interstellar medium. Six transition states are located by G4 and CBS‐QB3 methods as well as extrapolation using energies for structures optimized in CCSD(T)/aug‐cc‐pVDZ structures. We use two isodesmic reactions with two well‐established thermochemical computational schemes G4 and CBS‐QB3 to estimate energy enthalpy and Gibbs energy of formation as well as the entropy of the gas phase system. Our extrapolated electronic energies of species appearing in the isodesmic reactions produce independent values of thermodynamic quantities consistent with G4 and CBS‐QB3. © 2013 Wiley Periodicals, Inc.     

F2 dimer: Improved intermolecular potential energy surface using ab initio calculations

A computational study on the intermolecular potential energy of 44 different orientations of F₂ dimers is presented. Basis set superposition error (BSSE) corrected potential energy surface is calculated using the supermolecular approach at CCSD(T) and QCISD(T) levels of theory. The interaction energies obtained using the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets are extrapolated to the complete basis set limit using the latest extrapolation scheme. The basis set effect is checked and it is found that the extrapolated intermolecular energies provide the best compromise between the accuracy and computational cost. Among 1320 energy points of F₂–F₂ system covering more relative orientations, the most stable structure of the dimers was obtained with a well depth of ?146.62 cm^?1 that related to cross configuration, and the most unstable structure is related to linear orientation with a well depth of ?52.63 cm^?1. The calculated second virial coefficients are in good agreement with experimental data. The latest extrapolation scheme of the complete basis set limit at the CCSD(T) level of theory is used to determine the intermolecular potential energy surface of the F₂ dimer. Comparing the results obtained by the latest scheme with those by older schemes show that the new approach provides the best compromise between accuracy and computational cost.     

Theoretical investigation of the halogen bonded complexes between carbonyl bases and molecular chlorine

The halogen bonded complexes between six carbonyl bases and molecular chlorine are investigated theoretically. The interaction energies calculated at the CCSD(T)/aug‐cc‐pVTZ level range between ?1.61 and ?3.50 kcal mol^?1. These energies are related to the ionization potential, proton affinity, and also to the most negative values (V_s,min) on the electrostatic potential surface of the carbonyl bases. A symmetry adapted perturbation theory decomposition of the energies has been performed. The interaction results in an elongation of the Cl? Cl bond and a contraction of the CF and CH bonds accompanied by a blue shift of the ν(CH) vibrations. The properties of the Cl₂ molecules are discussed as a function of the σ*(Cl? Cl) occupation, the hybridization, and the occupation of the Rydberg orbitals of the two chlorine atoms. Our calculations predict a large enhancement of the infrared and Raman intensities of the ν(Cl? Cl) vibration on going from isolated to complexed Cl₂. © 2015 Wiley Periodicals, Inc.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

Absorption,resonance, and near‐resonance Raman studies of the tetracyanoquinodimethane neutral and its monoanion in terms of density functional theory and complete active space self‐consistent field methods

The electronic structure of the 1¹B_1u and 1²B_3u excited electronic states of the tetracyanoquinodimethane (TCNQ) neutral and its charged derivative are studied within the framework of complete active space self‐consistent field (CASSCF) and Becke's three‐parameter hybrid method with Lee–Yang–Parr correlation functional (B3LYP) methods applied to the level aug‐cc‐p‐VDZ basis set. Both CASSCF/aug‐cc‐p‐VDZ and B3LYP/aug‐cc‐p‐VDZ treatments provide the ground‐state and the excited state geometries; these are then used to assess the Franck–Condon (FC) parameters in the 1¹B_1u state of the neutral TCNQ and in the 1²B_3u state of the TCNQ monoanion. The quality of numerical results is then tested on the base of available experimental near‐resonance and resonance Raman data. The studies are performed in terms of the vibronic model, which takes both FC and mode‐mixing (Dushinsky) effects into account. This somewhat simplified vibronic model leads to very good agreement between the theory and the Raman experiments concerning both neutral TCNQ and its monoanion. In particular, the calculated excitation profiles of the ν₂ = 2215 cm^?1, ν₄ = 1389 cm^?1, ν₅ = 1195 cm^?1, and ν₉ = 336 cm^?1 fundamentals are shown to be in excellent agreement with those for the TCNQ monoanion. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Computational NMR coupling constants: Shifting and scaling factors for evaluating 1JCH

Optimized shifting and/or scaling factors for calculating one‐bond carbon–hydrogen spin–spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97‐2 and M06‐L) and basis sets (TZVP, HIII‐su3, EPR‐III, aug‐cc‐pVTZ‐J, ccJ‐pVDZ, ccJ‐pVTZ, ccJ‐pVQZ, pcJ‐2 and pcJ‐3) using 68 organic molecular systems with 88 ¹J_CH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of ¹J_CH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root‐mean‐square deviations lie between 4.7 and 16.4 Hz and are reduced to 4–6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97‐2 functionals in combination with HIII‐su3, aug‐cc‐pVTZ‐J and pcJ‐2 basis sets. Copyright © 2013 John Wiley & Sons, Ltd.     

Theoretical study on the gas phase reaction of propargyl alcohol with hydroxyl radical

The reaction of propargyl alcohol with hydroxyl radical has been studied extensively at CCSD(T)/aug‐cc‐pVTZ//MP2/cc‐pVTZ level. This is the first time to gain a conclusive insight into the reaction mechanism and kinetics for this important reaction in detail. Two reaction mechanisms were revealed, namely addition/elimination and hydrogen abstraction mechanism. The reaction mechanism confirms that OH addition to C?C triple bond forms the chemically activated adducts, IM1 (·CHCOHCH₂OH) and IM2 (CHOH·CCH₂OH), and the hydrogen abstraction pathways (? CH₂OH bonded to the carbon atom and alcohol hydrogen) may occur via low barriers. Harmonic model of Rice–Ramsperger–Kassel–Marcus theory and variational transition state theory are used to calculate the overall and individual rate constants over a wide range of temperatures and pressures. The calculated rate constants are in good agreement with the experimental data. At atmospheric pressure with Ar as bath gas, IM1 (·CHCOHCH₂OH) and IM2 (CHOH·CCH₂OH) formed by collisional stabilization are dominant in the low temperature range. The production of CHCCHOH + H₂O via hydrogen abstraction becomes dominate at higher temperature. The fraction of IM3 (CH₂COHCH₂·O) is very significant over the moderate temperature range. © 2014 Wiley Periodicals, Inc.     

ChemInform Abstract: Smallest Deltahedra Silicon Dicarbide: C2Si32‐.

A total of 20 singlet and 36 triplet C₂Si₃^2‐ isomers are obtained by quantum chemical calculations at the B3LYP/aug‐cc‐pVTZ level.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Theoretical investigation of the structures and spectroscopic properties of (H2O4)n (n = 1–4) clusters

Density functional theory and ab initio calculations were performed to elucidate the hydrogen interactions in (H₂O₄)_n (n = 1–4) clusters. The optimized geometries, binding energies, and harmonic vibrational frequencies were predicted at various levels of theory. The trans conformer of the H₂O₄ monomer was predicted to be the most stable structure at the CCSD(T)/aug‐cc‐pVTZ level of theory. The binding energies per H₂O₄ monomer increased in absolute value by 9.0, 10.1, and 11.8 kcal/mol from n = 2 to n = 4 at the MP2/cc‐pVTZ level of theory (after the zero‐point vibrational energy and basis set superposition error corrections). This result implies that the intermolecular hydrogen bonds were stronger in the long‐chain clusters, that is, the formation of the longer chain in the (H₂O₄)_n clusters was more energetically favorable.     

Divergent Pathways and Competitive Mechanisms of Metathesis Reactions between 3‐Arylprop‐2‐ynyl Esters and Aldehydes: An Experimental and Theoretical Study

Mechanistic studies of the reaction between 3‐arylprop‐2‐ynyl esters and aldehydes catalyzed by BF₃ ? Et₂O were performed by isotopic labeling experiments and quantum chemical calculations. The reactions are shown to proceed by either a classical alkyne–carbonyl metathesis route or an unprecedented addition–rearrangement cascade. Depending on the structure of the starting materials and the reaction conditions, the products of these reactions can be Morita–Baylis–Hillman (MBH) adducts that are unavailable by traditional MBH reactions or E‐ and Z‐α,β‐unsaturated ketones. ¹⁸O‐Labeling studies suggested the existence of two different reaction pathways to the products. These pathways were further examined by quantum chemical calculations that employed the DFT(wB97XD)/6‐311+G(2d,p) method, together with the conductor‐like screening model for realistic solvation (COSMO‐RS). By using the wB97XD functional, the accuracy of the computed data is estimated to be 1–2 kcal mol^?1, shown by the careful benchmarking of various DFT functionals against coupled cluster calculations at the CCSD(T)/aug‐cc‐pVTZ level of theory. Indeed, most of the experimental data were reproduced and explained by theory and it was convincingly shown that the branching point between the two distinct mechanisms is the formation of the first intermediate on the reaction pathway: either the four‐membered oxete or the six‐membered zwitterion. The deep mechanistic understanding of these reactions opens new synthetic avenues to chemically and biologically important α,β‐unsaturated ketones.     

Prediction of water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) using correlation‐consistent and polarization‐consistent basis sets in the Kohn–Sham basis set limit

Density functional theory (DFT) was used to estimate water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) in the Kohn–Sham (KS) complete basis set (CBS) limit. Correlation‐consistent cc‐pVxZ and cc‐pCVxZ (x = D, T, Q, 5, and 6), and their modified versions (ccJ‐pVxZ, unc‐ccJ‐pVxZ, and aug‐cc‐pVTZ‐J) and polarization‐consistent pc‐n and pcJ‐n (n = 0, 1, 2, 3, and 4) basis sets were used, and the results fitted with a simple mathematical formula. The performance of over 20 studied density functionals was assessed from comparison with the experiment. The agreement between the CBS DFT‐predicted isotropic shieldings, spin–spin values, and the experimental values was good and similar for the modified correlation‐consistent and polarization‐consistent basis sets. The BHandH method predicted the most accurate ¹H, ¹⁷O isotropic shieldings and ¹J(OH) coupling constant (deviations from experiment of about ? 0.2 and ? 1 ppm and 0.6 Hz, respectively). The performance of BHandH for predicting water isotropic shieldings and ¹J(OH) is similar to the more advanced methods, second‐order polarization propagator approximation (SOPPA) and SOPPA(CCSD), in the basis set limit. Copyright © 2008 John Wiley & Sons, Ltd.     

Nature of the interaction between ammonia derivatives and carbon disulfide. A theoretical investigation

The interactions between NH₃, its methylated and chlorinated derivatives and CS₂ are investigated by ab initio CCSD(T) and density functional BLYP‐D3 methods. The CCSD(T)/aug‐cc‐pVTZ calculated interaction energies of complexes characterized by the S···N chalcogen bonds range between ?1.71 and ?2.78 kcal mol^?1. The S···N bonds are studied by atoms in molecules, natural bond orbital, and noncovalent interaction methods. The lack of correlation between the interaction energies of methylated amines complexes and the electrostatic potential results from the lone pair effect in aliphatic amines. Different structures of CS₂ complexed with ammonia derivatives, stabilized by other than the S···N chalcogen bonds, are also predicted. These structures are characterized by interaction energies ranging between 1.15 and 3.46 kcal mol^?1. The results show that the complexing ability of CS₂ is not very high but this molecule is able to attack the electrophilic or nucleophilic sites of a guest molecule.     

Theoretical investigation on the IR spectra of four‐membered silicon oxide ring

To test the mass effect on the ring vibrational frequencies (SiO₂) shift of the four‐membered silicon oxide ring, the deuterium and the tritium substituted cyclodisiloxanes on the hydrogen positions are examined at the CCSD(T)/cc‐pVTZ level of theory. The SiO₂ ring vibrations for a silicon oxide surface model compound, substituted cyclodisiloxane (Si₃O₅? O₂? Si₃O₅) with two six‐membered ring, are also calculated at the B3LYP/cc‐pVTZ level of theory. Our results of 909 and 920 cm^?1 are in good agreement with the experimental result. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

硫的取代对基态硫代乙酸单分子解离反应影响的理论研究

应用B3LYP方法,结合6-31G**、cc-pVDZ、aug-cc-pVDZ和cc-pVTZ基组对硫代乙酸的两种异构体CH3C(O)SH和CH3C(S)OH在基态势能面上的9个单分子反应进行了研究。本文计算预测硫代乙酸主要以CH3C(O)SH的形式存在,两种异构体均以顺式构象为优势构象。通过对比CH3C(O)SH、CH3C(S)OH和 CH3C(O)OH的反应性差异,我们可以得出结论：CH3C(O)OH中-OH基团的O被S取代后,只有当-SH作为一个整体参加反应时才对分子解离过程有较大影响;而C=O或C=S对反应性影响较小。