Methyl radical also reacts by the frontside mechanism: an ab initio study of some homolytic substitution reactions of methyl radical at silicon,germanium and tin

Ab initio calculations using 6-311G**, cc-pVDZ, aug-cc-pVDZ, and a (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional (B3LYO) calculations predict that homolytic substitution reactions of the methyl radical at the silicon atom in disilane can proceed via both backside and frontside attack mechanisms. At the highest level of theory (CCSD(T)/aug-cc-pVDZ//MP2/aug-cc-pVDZ), energy barriers (delta E) of 47.4 and 48.6 kJ mol-1 are calculated for the backside and frontside reactions respectively. Similar results are obtained for reactions involving germanium and tin with energy barriers (delta E) of between 46.5 and 67.3, and 41.0 and 73.3 kJ mol-1 for the backside and frontside mechanisms, respectively. These data suggest that homolytic substitution reactions of methyl radical at silicon, germanium, and tin can proceed via either homolytic substitution mechanism.     

An ab-initio study of some homolytic substitution reactions of acyl radicals at silicon, germanium and tint

Ab initio calculations using the 6-311G**, cc-pVDZ, and (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional (BHandHLYP, B3LYP) calculations predict that homolytic substitution reactions of acetyl radicals at the silicon atoms in disilane can proceed via both backside and frontside attack mechanisms. At the highest level of theory (CCSD(T)/cc-pVDZ//MP2/cc-pVDZ), energy barriers (deltaE double dagger) of 77.2 and 81.9 kJ mol(-1) are calculated for the backside and frontside reactions respectively. Similar results are obtained for reactions involving germanium and tin with energy barriers (deltaE double dagger) of 53.7-84.2, and 55.0-89.7 kJ mol(-1) for the backside and frontside mechanisms, respectively. These data suggest that both homolytic substitution mechanisms are feasible for homolytic substitution reactions of acetyl radicals at silicon, germanium, and tin. BHandHLYP calculations provide geometries and energy barriers for backside and frontside transition states in good agreement with those obtained by traditional ab initio techniques.     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.     

Is the CH2OH + O2 → CH2 = O + HO2 Reaction Barrierless? An Ab Initio Study on the Reaction Mechanism

The reaction mechanism of the CH₂OH + O₂ gas-phase reaction was investigated by means of ab initio calculations. MP2 and QCISD methodologies were used to obtain the stationary points on the potential energy surface. Single-point high-level QCISD(T) calculations were performed over the QCISD results in order to refine the energy of the transition states and the minima found. A new transition state concerning the initial O₂ addition to the CH₂OH radical was found, not reported so far for this reaction. Extra CCSD optimisation and single-point high-level CCSD(T) calculations upon the QCISD results confirm this TS. Additional RASSCF calculations show that its wave function has no significant multireferential character.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

An ab initio and DFT study of some halogen atom transfer reactions from alkyl groups to acyl radical

Ab initio calculations using the 6-311G**, cc-pVDZ, and (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (HF) the inclusion of electron correlation, and density functional (BHandHLYP, B3LYP) calculations predict that the transition states for the reaction of acetyl radical with several alkyl halides adopt an almost collinear arrangement of attacking and leaving radicals at the halogen atom. Energy barriers (DeltaE(double dagger)) for these halogen transfer reactions of between 89.2 (chlorine transfer from methyl group) and 25.3 kJ mol(-1) (iodine transfer from tert-butyl group) are calculated at the BHandHLYP/DZP level of theory. While the difference in forward and reverse energy barriers for iodine transfer to acetyl radical is predicted to be 15.1 kJ mol(-1) for primary alkyl iodide, these values are calculated to be 6.7 and -4.2 kJ mol(-1) for secondary and tertiary alkyl iodide respectively. These data are in good agreement with available experimental data in that atom transfer radical carbonylation reactions are sluggish with primary alkyl iodides, but proceed smoothly with secondary and tertiary alkyl iodides. These calculations also predict that bromine transfer reactions involving acyl radical are also feasible at moderately high temperature.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

An evaluation of harmonic vibrational frequency scale factors

Scale factors for obtaining fundamental vibrational frequencies, low-frequency vibrational frequencies, zero-point vibrational energies (ZPVEs), and thermal contributions to enthalpy and entropy have been derived through a least-squares approach from harmonic frequencies determined at more than 100 levels of theory. Wave function procedures (HF, MP2, QCISD, QCISD(T), CCSD, and CCSD(T)) and a large and representative range of density functional theory (DFT) approaches (B3-LYP, BMK, EDF2, M05-2X, MPWB1K, O3-LYP, PBE, TPSS, etc.) have been examined in conjunction with basis sets such as 6-31G(d), 6-31+G(d,p), 6-31G(2df,p), 6-311+G(d,p), and 6-311+G(2df,p). The vibrational frequency scale factors were determined by a comparison of theoretical harmonic frequencies with the corresponding experimental fundamentals utilizing a standard set of 1066 individual vibrations. ZPVE scale factors were generally obtained from a comparison of the computed ZPVEs with experimental ZPVEs for a smaller standard set of 39 molecules, though the effect of expansion to a 48 molecule data set was also examined. In addition to evaluating the scale factors for a wide range of levels of theory, we have also probed the effect on scale factors of varying the percentage of incorporated exact exchange in hybrid DFT calculations using a modified B3-LYP functional. This has revealed a near-linear relationship between the magnitude of the scale factor and the proportion of exact exchange. Finally, we have investigated the effect of basis set size on HF, MP2, B3-LYP, and BMK scale factors by deriving values with basis sets ranging from 6-31G(d) up to 6-311++G(3df,3pd) as well as with basis sets in the cc-pVnZ and aug-cc-pVnZ series and with the TZV2P basis.     

A gauge-independent zeroth-order regular approximation to the exact relativistic Hamiltonian--formulation and applications

A simple modification of the zeroth-order regular approximation (ZORA) in relativistic theory is suggested to suppress its erroneous gauge dependence to a high level of approximation. The method, coined gauge-independent ZORA (ZORA-GI), can be easily installed in any existing nonrelativistic quantum chemical package by programming simple one-electron matrix elements for the quasirelativistic Hamiltonian. Results of benchmark calculations obtained with ZORA-GI at the Hartree-Fock (HF) and second-order Moller-Plesset perturbation theory (MP2) level for dihalogens X(2) (X=F,Cl,Br,I,At) are in good agreement with the results of four-component relativistic calculations (HF level) and experimental data (MP2 level). ZORA-GI calculations based on MP2 or coupled-cluster theory with single and double perturbations and a perturbative inclusion of triple excitations [CCSD(T)] lead to accurate atomization energies and molecular geometries for the tetroxides of group VIII elements. With ZORA-GI/CCSD(T), an improved estimate for the atomization energy of hassium (Z=108) tetroxide is obtained.     

Symmetry breaking and electron correlation in O2−, O2, and O2+: A comparison of coupled cluster and quadratic configuration interaction approaches

Potential energy curves are calculated for O₂⁻, O₂, and O₂⁺ at the CCSD, QCISD, CCSD(T), and QCISD(T) levels of theory using aug-cc-pVDZ and aug-cc-pVTZ basis sets with electron correlation built onto inversion symmetry constrained and relaxed UHF wave functions. The spectroscopic constant r_e, w_e, w_e, x_e, D_j, and α_e, are determined from the potential curves using standard second-order perturbation theory expressions and are compared with experimental values to assess the relative accuracy of the theoretical approaches. Comparison of corresponding symmetry-constrained and symmetry-relaxed calculations indicates that the CCSD method is generally superior to CCSD(T), QCISD, and QCISD(T) in recovering from a symmetry-broken reference function. © 1996 John Wiley & Sons, Inc.     

A coupled-cluster study of the structure and vibrational spectra of pyrazole and imidazole

Harmonic force fields were calculated at the corresponding optimized geometries for pyrazole and imidazole at the HF, B3LYP, MP2, CCSD and CCSD(T) levels using the 6-31G* basis set and at the HF and B3LYP levels using the cc-pVTZ basis set. The agreement between the calculated and experimental geometries by the CCSD and CCSD(T) methods was generally similar to that obtained with the B3LYP and MP2 methods. The force fields were scaled using one-scale-factor (1SF), 3SF and 7SF scaling schemes. The scale factors were varied with respect to the experimental frequencies. Using 7SF scaling, the root-mean-square (RMS) deviation of the calculated frequencies from the experimental frequencies by the HF, B3LYP, MP2, CCSD and CCSD(T) methods and the 6-31G* basis set was 16, 7, 13, 11 and 11 cm(-1), respectively. This shows that the B3LYP method is preferred for force field calculations over the perturbative MP2, CCSD and CCSD(T) methods. Using 1SF scaling, the CCSD(T) scale factor was 0.931, the highest among the five methods used but close to that obtained with the B3LYP method and the cc-pVTZ basis set with lower RMS deviation.     

Interaction of dihydrogen with small and light molecules

Second-order M?ller-Plesset (MP2) calculations (using the approximate resolution of the identity, RI-MP2), explicitly correlated MP2 (MP2-R12) calculations, and coupled-cluster calculations including all single and double excitations with a perturbative estimate of triple excitations [CCSD(T)] are performed to study the interaction of molecular hydrogen with the small molecules HF, H2O, NH3, and LiOH. Different adsorption positions are studied. In the cases of H2O and NH3, the most favorable configuration places H2 in an end-on fashion on the O or N atom, respectively. In the cases of HF and LiOH, the H2 molecule takes a side-on position on the H atom of HF or the Li atom. With respect to MP2 calculations in a triple-zeta basis, both the enlargement of the basis set and the extension of the correlation treatment (CCSD(T) vs MP2) increase the interaction energy. The basis set limit CCSD(T) estimates of the interaction energy of H2 with the HF, H2O, NH3, and LiOH molecules amount to 4.40, 2.67, 3.02, and 10.74 kJ mol-1, respectively. The interaction energy for the simultaneous interaction of H2 with two LiOH molecules does not significantly exceed the value obtained for the interaction with a single LiOH molecule. Furthermore, the interaction energies (by MP2) of H2 with glycine, the glycine dimer, and imidazolium chloride amount to 2.78, 5.00, and 6.30 kJ mol-1, respectively.     

Benzooxirene. Ab initio calculations

Ab initio calculations have been performed on benzooxirene, the corresponding oxo carbene (“ketocarbene”), and the transition state linking the two. At the highest level used, QCISD(T)/6-31G^*//MP2(FULL)/6-1G^* with MP2(FULL)/ 6-31G^* zero point energy corrections, the relative energies of the oxirene, the transition state and the carbene are 0, 24.6, and −17.8 kJ mol⁻¹. Correlation energy effects are very important in this system: at the QCISD(T) level the oxirene lies above the carbene, as at the MP4 and HF levels, but at the MP2 level the ordering is reversed. Benzooxirene is probably slightly nonplanar: the HF/6-31G^* geometry is C_2v but the MP2(Fermi contact)/6-31G^* geometry is C_s with a 6-/3-ring coplanarity deviation of about 6.9 °, although in the MP2(FULL)/6-31G^* geometry this is reduced to about 3.1 °.     

Ab Initio Calculations on the Vinyl Fluoride Fragmentation Reactions

High level ab initio calculations for the fragmentation reactions of vinyl fluoride were performed. The relative energies calculated at the QCISD(T)/6-311G(2d,2p) level of theory, corrected with MP2/6-311G(2d,2p) zero-point energies (ZPEs), differ significantly from those obtained previously at a lower level of theory. The calculations suggest that both the threeand four-center HF elimination processes are likely to occur, with the three-center elimination favored over the four-center at high energies.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

Theoretical conformational analysis for neurotransmitters in the gas phase and in aqueous solution. Norepinephrine

The natural neurotransmitter (R)-norepinephrine takes the monocationic form in 93% abundance at the physiological tissue pH of 7.4. Ab initio and DFT/B3LYP calculations were performed for 12 protonated conformers of (R)-norepinephrine in the gas phase with geometry optimizations up to the MP2/6-311++G level, and with single-point calculations up to the QCISD(T) level at the HF/6-31G-optimized geometries. Four monohydrates were studied at the MP2/6-31G//HF/6-31G level. In the gas phase, the G1 conformer is the most stable with phenyl.NH(3)(+) gauche and HO(alc).NH(3)(+) gauche arrangements. A strained intramolecular hydrogen bond was found for conformers (G1 and T) with close NH(3)(+) and OH groups. Upon rotation of the NH(3)(+) group as a whole unit about the C(beta)-C(alpha) axis, a 3-fold potential was calculated with free energies for barriers of 3-12 kcal/mol at the HF/6-31G level. Only small deviations were found in MP2/6-311++G single-point calculations. A 2-fold potential was calculated for the phenyl rotation with free energies of 11-13 kcal/mol for the barriers at T = 310 K and p = 1 atm. A molecular mechanics docking study of (R)-norepinephrine in a model binding pocket of the beta-adrenergic receptor shows that the ligand takes a conformation close to the T(3) arrangement. The effect of aqueous solvation was considered by the free energy perturbation method implemented in Monte Carlo simulations. There are 4-5 strongly bound water molecules in hydrogen bonds to the conformers. Although hydration stabilizes mostly the G2 form with gauche phenyl.NH(3)(+) arrangement and a water-exposed NH(3)(+) group, the conformer population becomes T > G1 > G2, in agreement with the PMR spectroscopy measurements by Solmajer et al. (Z. Naturforsch. 1983, 38c, 758). Solvent effects reduce the free energies for barriers to 3-6 and 9-12 kcal/mol for rotations about the C(beta)-C(alpha) and the C(1)(ring)-C(beta) axes, respectively.     

Peculiar basis set dependence of the energetics of C2S2H2 isomers. In search of adequate and affordable basis set for routine calculations

Quantum mechanical calculations using the Hartree-Fock (HF), post-Hartree-Fock {MP2 and CCSD(T)}and gradient-corrected hybrid density functional variant, B3LYP methodology in conjugation with basis sets like 6-31G*, 6-311G** and cc-pVTZ, were employed to critically probe the right quality of basis set and the appropriate level of theory that can be applied in modeling the organo-sulfur compounds. The relative stability ordering of 1,2-dithiete and dithioglyoxal reveal the delicate basis set dependence especially on ‘S’ atom. Single point calculations at the B3LYP, MP2 and CCSD(T) levels using a series of basis set ranging from double ζ quality all the way up to the quintuple ζ quality were done to illustrate the effect of electron correlation and the basis set dependency for these compounds. Basis set requirements are much higher than affordable for medium sized molecules as very slow convergence is seen even when the calculations are carried out with basis set as high as cc-pV5Z at the B3LYP and MP2 level of theories.     

Intermolecular potential energy surface for CS2 dimer

A new four‐dimensional intermolecular potential energy surface for CS₂ dimer is obtained by ab initio calculation of the interaction energies for a range of configurations and center‐of‐mass separation distances for the first time. The calculations were performed using the supermolecular approach at the Møller–Plesset second‐order perturbation (MP2) level of theory with the augmented correlation consistent basis sets (aug‐cc‐pVxZ, x = D, T) and corrected for the basis‐set superposition error using the full counterpoise correction method. A two‐point extrapolation method was used to extrapolate the calculated energy points to the complete basis set limit. The effect of using the higher levels of theory, quadratic configuration interaction containing single, double, and perturbative triple excitations QCISD(T) and coupled cluster singles, doubles and perturbative triples excitations CCSD(T), on the shape of potential energy surface was investigated. It is shown that the MP2 level of theory apparently performs extremely poorly for describing the intermolecular potential energy surface, overestimating the total energy by a factor of nearly 1.73 in comparison with the QCISD(T) and CCSD(T) values. The value of isotropic dipole–dipole dispersion coefficient (C₆) of CS₂ fluid was obtained from the extrapolated MP2 potential energy surface. The MP2 extrapolated energy points were fitted to well‐known analytical potential functions using two different methods to represent the potential energy surface analytically. The most stable configuration of the dimer was determined at R = 6.23 au, α = 90°, β = 90°, and γ = 90°, with a well depth of 3.980 kcal mol^?1 at the MP2 level of theory. Finally, the calculated second virial coefficients were compared with experimental values to test the quality of the presented potential energy surface. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011.     

NH3…H2O体系氢键团簇结构的从头计算方案比较

对H2O，NH3单分子以及反式线性NH3…H2O双分子团簇应用HF，DFT，MR2，CCD ，QCISD等从头计算方法在不同的基组水平上进行几何结构优化和氨—水分子间相 互作用能计算。比较实验值和计算结果表明：对氨—水体系，MP2／6-31+G^**方法 能够满足氢键N-0键长误差为0．004nm，键角∠NOH误差为1°，NH3…H2O相互作用 能△E与CCD和QCISD偏差2．1kJ／mol计算精度的要求，并有较合理的计算代价。     

O原子与HNCO反应机理的量子化学及电子密度拓扑研究

采用MP2(Full),B3LYP,QCISD/MP2和CCSD(T)/MP2方法在6-311G(d,p)水平上对O原子与HNCO反应的微观机理进行了理论研究．采用MP2(Full)和B3LYP对反应位能面上的各驻点进行几何构型的全优化,振动分析和IRC计算证实了中间体和过渡态的真实性和相互连接关系．四种方法计算得到了四个反应通道的反应活化能．研究表明,O原子进攻HNCO中的H原子为反应的主要反应通道,该通道中形成了两个分子复合物,电子密度拓扑分析表明这两个分子复合物均为氢键复合物．