The inner-sphere reorganization energies for AH2 + AH (A = Al,Si, P,S) electron self-exchange reactions in electron-transfer processes from ab initio calculations

This article presents an application of the accurate calculation scheme proposed recently for the inner-sphere reorganization energies of molecules of the type AH₂ (A = Al, Si, P, and S). A reasonable extension has been made. The inner-sphere reorganization energies for the title thermal electron self-exchange reactions are calculated in terms of ab initio MO self-consistent field method (HFSCF ) at different basis-set levels (6-31G **, 6-31 + G **, DZ , and DZP ) and the involved parameters are also determined. These calculated results have been calibrated by comparing optimized molecular geometrical parameters and corresponding energy properties with the experimental findings or other theoretical values. An approximation, in which the contribution from the bond length–bond angle to the potential energy surface is neglected, is adopted in constructing the calculation formulas via the function model. Its adequacy is discussed. Agreement among different calculation schemes is analyzed. © 1995 John Wiley & Sons, Inc.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.     

Nonempirical ab initio studies on inner-sphere reorganization energies of M2+(H2O)6/M3+(H2O)6 redox couples at valence basis level

On the basis of the basic feature of the electron transfer reactions, a new theoretical scheme and application of a nonempirical ab initio method in computing the inner-sphere reorganization energies (RE) of hydrated ions in electron transfer processes in solution are presented at valence STO basis (VSTO) level. The potential energy surfaces and the various molecular structural parameters for transition metal complexes are obtained using nonempirical molecular orbital (MO) calculations, and the results agree very well with experimentally observed ones from vibrational spectroscopic data. The results of inner-sphere REs obtained from these calculations via this new scheme give a good agreement with photoemission experimental findings and those from the improved self-exchange model proposed early for M²⁺(H₂O)₆/M³⁺(H₂O)₆(M = V, Cr, Mn, Fe, and Co) redox couple systems and are better than those from semiempirical INDO/II MO method and other classical methods. Further, the observed agreement of the optimized structural data and the results of inner-sphere REs of complexes with experimental findings confirms the following: (1) the validity of nonempirical MO calculation method to get accurate structural parameters and inner-sphere RE for the redox systems for which reliable vibrational spectroscopic data are not available, (2) the validity of the improved self-exchange model proposed early for inner-sphere RE, and (3) the reasonableness of some approximations adopted in this study. © 1997 John Wiley & Sons, Inc.     

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.     

离子簇合物HeN_2~+的电子光谱的从头算研究

本文报道了离子簇合物HeN2＋的电子光谱的从头算结果。在MCSCF／6－31G（d，p）水平上，对其基态进行了几何优化，用该构型．在SDC1／6－31（d，p）水平上计算了基态的总能量。用SDCI方法计算得到HeN2＋从基态到选择激发态的垂直跃迁能、振子强度、跃迁频率、辐射寿命以及Einstain系数。该结果可以较好的验证maier的实验。     

Structures and energies of D-galactose and galabiose conformers as calculated by ab initio and semiempirical methods

Optimized geometries and total energies of some conformers of alpha- and beta-D-galactose have been calculated using the RHF/6-31G* ab initio method. Vibrational frequencies were computed at the 6-31G* level for the conformers that favor internal hydrogen bonding, in order to evaluate their enthalpies, entropies, Gibbs free energies, and then their structural stabilities. The semiempirical AM1, PM3, MNDO methods have also been performed on the conformers GG, GT, and TG of alpha- and beta-D-galactose. In order to test the reliability of each semiempirical method, the obtained structures and energies from the AM1, PM3, and MNDO methods have been compared to those achieved using the RHF/6-31G* ab initio method. The MNDO method has not been investigated further, because of the large deviation in the structural parameters compared with those obtained by the ab initio method for the galactose. The semiempirical method that has yielded the best results is AM1, and it has been chosen to perform structural and energy calculations on the galabiose molecule (the disaccharides constituted by two galactose units alpha 1,4 linked). The goal of such calculations is to draw the energy surface maps for this disaccharide. To realize each map, 144 different possible conformations resulting from the rotations of the two torsional angles psi and phi of the glycosidic linkage are considered. In each calculation, at each increment of psi and phi, using a step of 30 degrees from 0 to 330 degrees, the energy optimization is employed. In this article, we report also calculations concerning the galabiose molecule using different ab initio levels such as RHF/6-31G*, RHF/6-31G**, and B3Lyp/6-31G*.     

Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths

The minimum energy path (MEP) of the reaction, CF₃CHFCF₃ + H → transition state (TS) → CF₃CFCF₃ + H₂, has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6‐31++G**, BH&HLYP/cc‐pVDZ, BMK/6‐31++G**, M05/6‐31+G**, M05‐2X/6‐31+G**, UMP2/6‐31++G**, PUMP2/6‐31++G**//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVDZ//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVTZ(spd,sp)//UMP2//6‐31++G**, RCCSD(T)/CBS//M05/6‐31+G**, and RCCSD(T)/CBS//UMP2/6‐31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero‐curvature, and small‐curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000–1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6‐31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (～1500 K), SCT corrections are significant at low temperatures (～300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections. © 2012 Wiley Periodicals, Inc.     

CH and Chalogen bond dissociation energies for fluorinated and chlorinated methane evaluated with hybrid B3LYP density functional theory methods and their comparison with experimental data and the CBS-Q ab initio computational approach

The CH and Chalogen bond dissociation energies (BDEs) were computed with the hybrid B3LYP/6-311 + G(2d,2p) theory model for chlorinated and fluorinated methane. All computed values were substantially lower (5–10 kcal mol⁻¹) than the experimental values. To obtain better agreement, a correlation factor was introduced. When this factor was applied, excellent agreement between the B3LYP/6-311 + G(2d,2p) computed energies and the experimental BDEs was observed. On the other hand, the CBS-Q ab initio computational approach generated BDEs which are in good agreement with experimental values without a correction factor.     

HeLi^n^+(n=0, 1)簇合物的稳定性和垂直激发态光谱

本文用ab initio研究了簇合物HeLi^n^+(n=0, 1)的几何构型和成键性质。在MP2(FULL)/6-31G**, 水平优化所得LeLi^+的平衡键长为0.2062nm, 与实验值0.205nm十分吻合。比较了HeLi^+(X^1∑^+和a^3∑^+), HeLi(X^2∑^+和a^4II)以及HLi(X^1∑^+)的稳定性, 计算了HeLi^+基态的相关能, 势能曲线和垂直激发态光谱。计算采用了6-31G**, 6-311G**,6-311G(2df, 2pd), 6-311G(3df, 2pd)和6-311+G(3df, 2pd)基组; 采用的方法包括MP2(FULL), MP4, MCSCF, MRSDCI, CCD和ST4CCD。计算表明, 同价HeLi^n^+中激发态的离解能均远比基态的大, 其中HeLi^+(a^3∑^+)的离解能最高(60.49kj/mol),说明激发态是稳定束缚态。HeLi^+基态比等电子体HLi分子基态的稳定性小得多。HeLi^+由A^1∑^+到B^1II的垂直跃迁(3σ→1π)振子强度较大而垂直跃迁能较小。     

An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

Theoretical study of the energetics of the reactions of triplet dioxygen with hydroquinone, semiquinone, and their protonated forms: relation to the mechanism of superoxide generation in the respiratory chain

One-electron reduction of the dioxygen molecule by the reduced form of mitochondrial ubiquinones (Q) of the NADH dehydrogenase (complex I) and mitochondrial cytochrome bc1 (complex III) is believed to be the main source of the superoxide anion radical O2*- and the hydroperoxide radical OOH*. In this work, we modeled the energetics of four possible reactions of the triplet ((3)Sigma(g)) dioxygen-molecule reduction by fully reduced and protonated ubiquinone (QH2; reaction 1), its deprotonated form (QH-; reaction 2), the semiquinone radical (QH*; reaction 3), and the semiquinone anion radical (Q*-; reaction 4), by means of ab initio calculations with the 6-31G(d) and 6-31+G(d) basis set in the restricted open-shell Hartree-Fock (ROHF), unrestricted Hartree-Fock (UHF), and complete active space self-consistent field (CASSCF) with dynamic correlation [at the second-order M?ller-Plesset (MP2) or multiple reference M?ller-Plesset (MRMP), respectively] schemes and the basis set superposition error (BSSE) correction included, as well as semiempirical AM1 and PM3 calculations in the UHF and ROHF schemes. 2-Butene-1,4-dione and p-benzoquinone were selected as model compounds. For the reduced forms of both compounds, reaction 1 turned out to be energetically unfavorable at all levels of theory, this agreeing with the experimentally observed diminished reductive properties of hydroquinone derivatives at low pH. For 2-butene-1,4-dione treated at the most advanced MRMP/CASSCF/6-31+G(d) level, the energies of reactions 1-4 are 4.7, -34.3, -15.0, and -4.1 kcal/mol, respectively. This finding suggests that reactions 2 and 3 are the most likely mechanisms of electron transfer to molecular oxygen in aprotic environments and that proton transfer is involved in this process. Nearly the same energies of reactions 2 and 3 were calculated at the MRMP/CASSCF/6-31+G(d) level for reduced forms of p-benzoquinone. Inclusion of diffuse functions in the basis set and dynamic correlation at the CASSCF level appears essential. Because deprotonated ubiquinol is unlikely to exist in physiological environments, reaction 3 appears to be the most likely mechanism of one-electron reduction of oxygen; however, if oxygen can penetrate cytochrome bc1 as far as the Q(o) center where ubiquinol can be deprotonated, reaction 2 can also come into play. The energies of reactions 2 and 3 calculated at the MRMP/CASSCF/6-31+G(d) level are most closely reproduced in the ab initio and semiempirical UHF PM3 calculations. Additional semiempirical calculations on more realistic models of ubiquinone, 2,3-dimethoxy-6-methyl-p-benzoquinone and 2,3-dimethoxy-5-isoprenyl-6-methyl-p-benzoquinone, gave qualitatively the same relations between the energies of reactions 2 and 3 as those carried out for p-benzoquinone species, thereby suggesting that this method could be used in studying electron-transfer reactions from reduced quinone derivatives to molecular oxygen in more complex systems, such as a model of the Q(o) site of cytochrome bc1, where applying ab initio methods is unfeasible.     

2-丁基-四氢噻吩亚砜13C-NMR的理论研究

在RHF/6-31G和B3LYP/6-31G水平上对顺式(Cis-)与反式(Trans-)2-丁基-四氢噻吩亚砜(BTHTO)进行几何优化,应用规范不变原子轨道法(GIAO)在6-31G、6-31+G、6-31++G和6-31+G(2d,p)水平上计算了Cis-和Trans-BTHTO的13C-NMR,对13C-NMR谱进行了归属。结果表明,BTHTO噻吩五元环的稳定构象呈半椅式,Cis-和Trans-BTHTO中与硫原子直接碳原子13C-NMR的显著差异主要是由于空间构型不同引起分子的静电势场对相应碳原子的屏蔽作用不同所致。     

Proton affinities of molecules containing nitrogen and oxygen: comparing ab initio and semiempirical results to experiments

We report a comparison of theoretical and experimental proton affinities at nitrogen and oxygen sites within a series of small molecules. The calculated proton affinities are determined using the semiempirical methods AM 1, MNDO , and PM 3; the ab initio Hartree–Fock method at the following basis levels: 3-21G //3-21G , 3-21+G //3-21G , 6-31G *//6-31G *, and 6-31+G (d, p)//6-31G *; and Møller–Plesset perturbation calculations: MP 2/6-31G *//6-31G *, MP 3/6-31G *//6-31G *, MP 2/6-31G +(d, p)//6-31G *, MP 3/6-31G +(d, p)//6-31G *, and MP 4(SDTQ )/6-31G +G (d, p)//6-31G *. The semiempirical methods have more nonsystematic scatter from the experimental values, compared to even the minimal 3-21G level ab initio calculations. The thermodynamically corrected 6-31G *//6-31G * proton affinities provide acceptable results compared to experiment, and we see no significant improvement over 6-31G *//6-31G * in the proton affinities with any of the higher-level calculations. © 1992 John Wiley & Sons, Inc.     

Theoretical studies on the S-N interactions in sulfoximine

The potential energy surface of sulfoximines has been searched using ab initio MO and Density Functional Calculations. The electronic structures of the isomers of sulfoximine have been studied using HF/6-31+G*, MP2(full)/6-31+G* and B3LYP/6-31+G* levels. Final energies of these molecules have been calculated at the high accuracy G2 and CBS-Q levels. Though a formal SN double bond is generally considered between sulfur and nitrogen in these systems, theoretical studies do not show any π interaction between them. S-N rotational barriers, bond dissociation energies, atomic charge analysis, and NBO analysis all indicate only a single bond across S-N with a very strong ionic interaction.     

Density functional efficiency in the calculations of vibrational frequencies and molecular structures of β-diketones

Density functional theory (DFT) levels are employed to calculate the vibrational frequencies and geometrical data of β-diketones. We evaluate the relative performance of the different levels by comparing theoretical results to experimental values. The applied DFT levels in this work are B3LYP, BLYP, B3P86, B3PW91, BPW91, G96LYP, BP86, and G96PW91 with the standard 6-31G, 6-31G*, 6-31G**, 6-31+G**, 6-31++G**, 6-311G**, 6-311++G** basis sets. The best results are obtained at the B3LYP, B3PW91, and B3P86 levels.     

Force-field parameters for the simulation of tetrahedral intermediates of serine hydrolases

CHARMM force-field parameters are reported for the tetrahedral intermediate of serine hydrolases. The fitting follows the standard protocol proposed for CHARMM22. The reference data include ab initio (RHF/6-31G*) interaction energies of complexes between water and the model compound 1,1-dimethoxyethoxide, torsional profiles of related model compounds from correlated ab initio (MP2/6-311+G*//B3LYP/6-31+G*) calculations, as well as molecular geometries and vibrational frequencies from density functional theory (B3LYP/6-31+G*). The optimized parameters reproduce the target data well. Their utility is demonstrated by a QM/MM study of the tetrahedral intermediate in Bacillus subtilis lipase A, and by classical molecular modeling of enantioselectivity in Pseudomonas aeruginosa lipase and its mutants.     

HCO＋HCN→HCHO＋CN反应机理研究

用ａｂｉｎｉｔｉｏ方法在ＵＨＦ／６－３１Ｇ基组上研究了ＨＣＯ　＋ＨＣＮ→ＨＣＨＯ＋ＣＮ　反应的反应机理。经寻找反应的过渡态和计算内禀反应坐标（ＩＲＣ）表明，该反应经一过渡态而形成产物，具有较宽的位垒，并属于后位垒反应。用ＵＭＰ２∥６－３１Ｇ方法计算的活化位垒为１４９．６０２ｋＪ／ｍｏｌ，与实验值（１４３．００１ｋＪ／ｍｏｌ）相当一致。     

Quantum-Chemical Study of Electronic and Steric Structure of Vinyltetrazoles: I. 2-Substituted 5-Vinyltetrazoles

The total Mulliken charges on the carbon atoms of the vinyl group, populations of S-trans-(N¹)conformers, and internal rotation energies were calculated ab initio (HF/6-31G**, MP2/6-31G**, and MP2/6-31G**//AM1) for a series of 2R-5-vinyltetrazoles (R = CH₃, C₂H₅, i-C₃H₇, t-C₄H₉, C₆H₅). The calculation results were compared to the available experimental data.     

Scaling factors for vibrational frequencies and zero-point vibrational energies of some recently developed exchange-correlation functionals

The scaling factors for the vibrational frequencies and zero-point vibrational energies evaluated at various combinations of recently developed exchange-correlation functionals and various basis sets are reported. The exchange-correlation functionals considered are B972, B98, HCTH, OLYP, O3LYP, G96LYP, PBE0 and VSXC functionals; the basis sets employed are 3-21G, 6-31G*, 6-31G**, 6-31+G, 6-311G*, 6-311G**, 6-311G(df,p), 6-311+G(df,p), cc-pVDZ and aug-cc-pVDZ. The experimental harmonic frequencies of 122 small molecules and the zero-point vibrational energies of 39 small molecules are used to determine the scaling factors through the least-square fitting procedure. It was found that the scaling factors do not depend significantly on the basis sets considered. The vibrational frequency scaling factors evaluated by using the B98 and PBE0 functionals are recommended on the basis of smallest root mean square error. The zero-point vibrational energy scaling factors evaluated from the B972 functional with Pople's double-zeta basis set and the HCTH functional with Pople's triple-zeta basis set are recommended on the basis of smallest root mean square error.     

Theoretical study of the reactions of the 1Σ+ ground state of MS+ (M = Sc,Y, and La) with oxygen‐transfer reagent MS+ + CO → ScO+ + CS in the gas phase

The reaction mechanisms of the ¹Σ⁺ ground state of MS⁺ (M = Sc, Y, and La) with oxygen‐transfer reagent MS⁺ + CO → MO⁺ + CS in the gas phase has been proposed and investigated by ab initio methods with the 6‐31G* basis set for nonmetal atoms and the effective core potentials of Lanl2dz for the metal atoms. A carbon migration from oxygen atom to sulfur atom via a four‐center transition state is involved on the reaction potential surface. The activation energies of the reactions are 34.0, 24.1, and 36.7 kcal/mol relative to their corresponding reactants and the reaction heats are 15.7, 18.6, and 18.0 kcal/mol (respectively, for M = Sc, Y, and La) at the MP4 (SDTQ)/6‐31G*//MP2/6‐31G* level plus zero‐point energy, which indicates that the cationic yttrium sulfide is more favorable for this type of reaction. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003