CH3F分子间相互作用的ab initio研究

在MP2/6-311++G(3d,3p)电子相关校正水平上, 对CH3F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证, 发现3种势能面上有极小点的构型. 进一步在高级电子相关校正的MP4SDTQ、 CCSD(T)/6-311++G(3df,3pd)方法水平上, 对其中总能量最小的构型进行精确计算, 得到二聚体的结合能为-9.707 kJ/mol. 研究结果支持了由光谱实验结果推测的构型, 解释了CH3F二聚体的谐振频率的多样性.      

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

硝酸乙酯分子间相互作用的ab initio研究

在abinitio-HF/6-31G水平上求得硝酸乙酯二聚体势能面上的四种优化构型和电子结构。经MP2电子相关校正和基组叠加误差(BSSE)以及零点能(ZPE)校正，求得二聚体的最大结合能为11.46kJ.mol^-^1，还进行HF/6-311G和HF/6-311++G水平的总能量比较计算，发现6-31G基组对计算结合能比较适合，二子体系间的电荷转移很少，对优化构型进行振动分析，并基于统计热力学求得从单体形成二聚体的热力学性质变化。     

硫代甲酰胺双聚体的量子化学计算

在MP2/6 31G(d)和MP2(FC)/6 311 G(d,p)水平上,对硫代甲酰胺(HC-SNH2)及其3种构型双聚体进行几何全优化计算,经振动频率分析,确认为势能超曲面上的稳定驻点.然后在MP2/6 311 G(2df,2p)水平上进行单点能计算和基组重叠误差(BSSE)校正以获得相互作用能.并利用自然键轨道(NBO)理论和分子中的原子(AIM)理论探讨HCSNH2之间相互作用的本质.     

2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

甲醛与乙腈相互作用的理论研究

在HF/6-311G(d,p)、 MP2/6-311G(d,p)和B3LYP/6-311G(d,p)水平上,对H2CO和CH3CN以及设计的4种结构H2CO…CH3CN复合物等进行几何全优化和振动频率计算,排除振动频率为负值的非局域极小点结构,并对稳定的环状构型复合物结合能进行基组重叠误差校正和零点振动能校正.分子间相互作用的能量分解分析显示,静电能在H2CO...CH3CN相互作用能量中占主导地位,电荷转移能居第二位.     

GeH4与HX(X＝F,Cl,Br)间二氢键复合物的理论研究

对GeH4与HX形成的二氢键复合物的结构特征及本质进行了探讨.在MP2/6-311 ++G(3 d,3p)水平优化、频率验证得到复合物的分子结构,通过分子的几何参数及电子密度拓扑分析,确认GeH4与卤化氢已形成了二氢键复合物.MP2/6-311 ++ G(3d,3p)水平下进行BSSE校正后的结合能为3.281到4.5...     

硝酸甲酯分子间相互作用的DFT和ab initio比较

用密度泛函理论(DFT)和从头算(ab initio)方法,分别在B3LYP/6 31G和HF/6 31G水平上求得硝酸甲酯三种二聚体的全优化几何构型和电子结构,并用6 311G和6 311＋＋G基组进行总能量计算.对HF/6 31G计算结果进行MP4SDTQ电子相关校正.在各基组下均进行基组叠加误差(BSSE)和零点能(ZPE)校正求得结合能.对6 31G优化构型作振动分析并基于统计热力学求得200～600 K温度下单体和二聚体的热力学性质.详细比较两种方法的相应计算结果,发现DFT求得的分子间距离较短,分子内键长较长,所得结合能均小于相应ab initio计算值.     

亚甲基氰自由基(HCCN)与一氧化氮(NO)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE水平上对反应HCCN+NO的二重态反应势能面进行了计算,得到了4种产物:P1(HCN+NCO),P2(OH+NCCN),P₃[HCN+(CNO)]和P4(HCN+CNO).其中产物P1为主要产物,P2为次要产物,P₃和P4很难得到.在G2(B3LYP/MP2/CC)水平,对产物P1和P2的反应通道的单点能量进行了校正.     

乙烯、乙炔和苯与氨基锂间的π型锂键相互作用

分别在DFT-B3LYP和MP2/6-311++G**水平上求得乙烯、乙炔和苯与氨基锂锂键复合物势能面上的3个几何构型. 频率分析表明,3个构型均为稳定构型. 计算结果表明,形成锂键复合物后,质子供体N(2)- Li(4)的键长明显增大,且其伸缩振动的频率发生了不同程度的红移. 分别在乙烯…氨基锂、乙炔…氨基锂和苯…氨基锂三种复合物中,经MP2/6-311++G**水平计算的同时含基组重叠误差(BSSE)和零点振动能校正的单体间锂键相互作用能分别为-26.04、-24.86 和 -30.02 kJ·mol-1. 自然键轨道理论（NBO）分析表明单体间的弱相互作用属于π-s型锂键.     

Analysis of vibrational spectra of 3-halo-1-propanols CH(2)XCH(2)CH(2)OH (X is Cl and Br)

The conformational stability and the three rotor internal rotations in 3-chloro- and 3-bromo-1-propanols were investigated by DFT-B3LYP/6-311+G and ab initio MP2/6-311+G, MP3/6-311+G and MP4(SDTQ)//MP3/6-311+G levels of theory. On the calculated potential energy surface twelve distinct minima were located all of which were not predicted to have imaginary frequencies at the B3LYP level of theory. The calculated lowest energy minimum in the potential curves of both molecules was predicted to correspond to the Gauche-gauche-trans (Ggt) conformer in excellent agreement with earlier microwave and electron diffraction results. The equilibrium constants for the conformational interconversion of the two 3-halo-1-propanols were calculated at the B3LYP/6-311+G level of calculation and found to correspond to an equilibrium mixture of about 32% Ggt, 18% Ggg1, 13% Tgt, 8% Tgg and 8% Gtt conformations for 3-chloro-1-propanol and 34% Ggt, 15% Tgt, 13% Ggg1, 9% Tgg and 7% Gtt conformations for 3-bromo-1-propanol at 298.15K. The nature of the high energy conformations was verified by carrying out solvent experiments using formamide ( epsilon=109.5) and MP3 and MP4//MP3 calculations. The vibrational frequencies of each molecule in its three most stable forms were computed at the B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecules.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

Ab initio prediction of MnOH2/MnOH2 electron transfer reactivity at electron correlation level

An ab initio calculation of electron transfer reactivity of Mn²⁺OH₂/Mn³⁺OH₂ redox system, in gas and in solution, has been made in this paper. The detailed geometry optimization and the scanning of the potential energy surfaces have been carried out at UMP2 (full)/6-311+G^* level. The relevant energy quantities (such as the activation energy and the binding energy) are calculated at different levels of theory (HF, MP2, MP3, MP4 and QCISD and corresponding spin-projection PUHF, PMP2 and PMP3) with the same basis set (6-311+G^*). Both all-electrons and valence electrons have been correlated in energy calculations. The electronic transmission coefficient is calculated using the ab initio potential energy surface slopes and the coupling matrix element determined from the two-state model and the Slater-type d-electron wave functions. The pair distribution function is obtained using both a classical sphere approximation scheme and a novel scheme. The relevant kinetic parameters are obtained at different ab initio levels in terms of the new model. The contact-distance dependence of these parameters and the applicability of the presented models are also discussed.     

Electronic structure and vibrational spectra of X20H20 (X=C, Si, Ge, Sn): A theoretical study

Optimized geometries, HOMO–LUMO gaps, vertical ionization potentials and electron affinities are obtained using HF, and B3LYP methods with 6-311G** basis set for C₂₀H₂₀, Si₂₀H₂₀ and Ge₂₀H₂₀. For germanium and tin analogues, B3LYP calculations are performed with LANL2DZ effective core potential. Electron correlation is included by doing MP2 calculation. The harmonic frequencies of all the compounds are obtained using B3LYP with 6-311G** and/or LANL2DZ basis sets. The force field and vibrational spectra are analyzed and 74 symmetry unique non-redundant local force constants are evaluated. Probable assignments are proposed for all the fundamentals based on the potential energy distribution.     

Three rotor potential energy scans, conformational equilibrium constants and vibrational analysis of 3-fluoro-1-propanol CH(2)FCH(2)CH(2)OH

The conformational stability and the three rotor internal rotations in 3-fluoro-1-propanol were investigated by the DFT-B3LYP/6-311+G** and the ab initio MP2/6-311+G** levels of theory. The calculated potential energy curves of the molecule at both levels of theory were consistent with complex conformational equilibria of about 12 minima, all of which were predicted to have real frequencies at both the B3LYP and the MP2 levels. The lowest energy minimum in the potential curves of 3-fluoro-1-propanol was predicted to correspond to the Gauche-gauche-trans (Ggt) conformer in excellent agreement with microwave and electron diffraction results. The equilibrium constants for the conformational interconversion of the molecule were calculated and found to correspond to an equilibrium mixture of about 33% Ggt, 14% Ggg1 and 13% Gg1g and about 43% Ggt, 12% Ggg1 and 10% Gg1g distribution by the B3LYP/6-311+G** and the MP2/6-311+G** calculations, respectively, at 298.15K. The vibrational frequencies of each molecule in its three stable forms were computed at B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecule.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

Distance dependence analysis of the electron transfer reactivity for the metastable FeOH2/FeOH2 system: an ab initio investigation

A novel theoretical scheme and ab initio application in discussing the electron transfer (ET) reactivity are presented in this paper, and are also calibrated in terms of the mono-hydrated iron ion system, Fe²⁺–OH₂/Fe³⁺–OH₂. The detailed geometry optimizations have been made at UMP2(full)/6-311+G^* level, and the activation geometrical configuration and the energy have been obtained at this level of theory using the activation model and the ab initio potential energy surface fitted from MP2(full)/6-311+G^* single point energies. The corresponding energy quantities (such as the activation energy, and dissociation energy) have also been obtained at different levels of theory (HF, MP2, MP3, MP4, QCISD and PUHF, PMP2 and PMP3 with the spin-projection) and a same basis set (6-311+G^*). The electron correlation calculations include the all electron correlation and the valence electron correlation. The electronic transmission coefficient is calculated using the ab initio potential energy surface slopes and the coupling matrix element determined from the two-state model and the Slater-type d-electron wave functions. The pair distribution function is calculated using two different schemes. Taking the pair distribution function and the local ET rate into account, a statistically averaged overall observed ET rate scheme and a spherically averaged local ET rate scheme are proposed. The relevant kinetic parameters are obtained in terms of these new schemes at different ab initio calculational levels. The contact distance dependence of these parameters and the applicability of the presented models and ab initio calculational method are also discussed.     

Basis set effects on the intermolecular interaction of the H2-H2 system obtained using ab initio molecular orbital calculations with the Møller-Plesset perturbation correction

The intermolecular interaction potential of the H₂-H₂ system was calculated by an ab initio molecular orbital method using several basis sets (up to 6-31 lG(3pd)) with inclusion of the electron correlation correction of the Møller-Plesset perturbation method and the basis set superposition error (BSSE) correction of the counterpoise method in order to evaluate the basis set effect. The calculated interaction energies depend strongly on the basis set used. Whereas the interaction energies of the repulsive and coulombic energy components calculated at the Hartree-Fock level are not affected by a change of basis set, the dispersion energy component depends strongly on the basis set used. Parameters of an exp-6-1 type non-bonding interaction potential were optimized on the basis of the MP4(SDTQ)/6-311G(3p) level intermolecular interaction energies of the H₂-H₂ system.     

Basis set dependence of the vibrational wavenumbers of the out-of-plane modes of conjugated π-electron ring systems

Basis set dependence of the vibrational wavenumbers of out-of-plane modes calculated at the MP2 level of ab initio molecular orbital theory is examined for benzene, p-benzoquinodimethane, p-benzoquinone, furan, and thiophene. Various polarization functions up to (3df,p) are used in combination with the 6-31G and 6-311G basis sets. It is shown that, especially in the case of normal modes with alternate out-of-plane motions of the carbon atoms (such as the ν₄ (b_2g) mode of benzene), the calculated wavenumbers depend strongly on the exponents (_d) of the d functions on the carbon atoms. It is therefore necessary to include d functions with an optimum exponent (_d0.4) on the carbon atoms to obtain reasonable out-of-plane vibrational force fields. In a few cases (such as the ν₁₆ (a₂) mode of furan), inclusion of a set of f functions on the carbon atoms has some effects on the calculated wavenumbers of out-of-plane modes. However, unless the basis set contains an optimum set of d functions, inclusion of a set of f functions does not improve the agreement between the observed and calculated vibrational wavenumbers. As a case with an exaggerated effect of basis set, it is shown that the wavenumber of the ν₄₁ (b_2g) mode of the planar optimized structure of p-benzoquinodimethane is calculated to be imaginary by using the 6-311G(d,p) or 6-311G(df,p) basis set at the MP2 level. For all the molecules treated in the present study, reasonable out-of-plane vibrational force fields are obtained by using the 6-31G(2df,p) and 6-311G(2df,p) basis sets.