HF,Ne等电子体系对内对间电子相关能的比较研究

采用MELD程序ROHF－OPT1方法在MP2/6-311++G(d)水平上,计算了HF分子基态1^∑和Ne原子基态1^S的对内对间电子相关能,并对两等电子体系的对相关能进行了分析和比较,深入研究两等电子体系对内对间相关能所具有的共同规律和存在的差异性,通过比较说明分子内的化学键是影响电子相关能的重要因素之一。HF分子和Ne原子两体系三重激发和四重激发对体系电子相关能的贡献的计算结果表明高激发项对电子相关能的贡献在精确量子化学计算中是不可以忽略的。     

HF,Ne等电子体系对内对间电子相关能的比较研究

采用MELD程序ROHF－OPT1方法在MP2/6-311++G(d)水平上,计算了HF分子基态1^∑和Ne原子基态1^S的对内对间电子相关能,并对两等电子体系的对相关能进行了分析和比较,深入研究两等电子体系对内对间相关能所具有的共同规律和存在的差异性,通过比较说明分子内的化学键是影响电子相关能的重要因素之一。HF分子和Ne原子两体系三重激发和四重激发对体系电子相关能的贡献的计算结果表明高激发项对电子相关能的贡献在精确量子化学计算中是不可以忽略的。     

多原子强离子键KX(X = OH, NC)及(KX)2体系电子相关能的简单计算

通过MELD程序MP2-OPT2方法在6-311++G(d)及6-311++G(3df, 3pd)基组下对KX(X = OH, NC)体系对内对间电子相关能的计算结果和分析, 并根据KX(X = OH, NC)分子中K和X组分的内层电子对相关能和内层电子相关效应的传递性, 定义并且计算了K^δ +和K^δ -对KX分子体系的电子相关能贡献. 基于K⁺, X-和KX(X = OH, NC)的对相关能的比较结果, 提出了利用其组成金属离子和负离子基团的电子相关能之和来计算多原子强离子键分子体系电子相关能的简单计算方法, 此方法对Gaussian98程序MP2(full)/6-311++G(d)水平下二聚体系(KOH))₂和(KNC))₂的电子相关能的近似计算结果误差很小, 应用此“化整为零”的方案计算大分子离子键体系电子相关能既可以基本达到化学计算要求的精度又能节约大量的计算工作量.     

应用Monte Carlo方法计算He原子包含电子相关作用的基态能量

应用Monte Carlo方法计算He原子包含电子相关波函数的基态能量,获得了与精确值非常接近的结果．实践表明,应用Monte Carlo方法有可能在多电子体系中直接采用包含任意2个电子间距离ry的函数作为变分函数来考虑电子相关作用．     

线性氟代烷烃电子相关能的基团加和性

应用Meld 程序在MP2-OPT2/6-311++G(d)水平计算了线性氟代烷烃分子CH3(CH2)mF(m=0-5)中—CH3、—CH2—和—F基团电子相关能贡献值. 计算结果表明, 在CH3(CH2)mF(m=0-5)体系中, 端基—F、—CH3基团电子相关能贡献值Ecorr(—F)和Ecorr(—CH3)随着m的增大而逐渐减小. 分子中α位置—CH2—基团电子相关能贡献值Ecorr(—CH2—)大于其他位置的贡献值. 通过计算结果可以推断, 在CH3(CH2)mF体系中随着m的逐渐增大, 远离端基—F的—CH2—基团电子相关能贡献值逐渐减小并将趋于不变, 此—CH2—基团可看作一个标准的亚甲基且其Ecorr(—CH2—)的数值在CH3(CH2)mF体系中具有传递性. 应用Meld程序在MP2-OPT2/6-311++G(d)水平对CH3(CH2)mF(m=2-5)体系的计算结果和应用Gaussian 98程序在MP2/6-311++G(d)//HF/6-311++G(d)水平对CH3(CH2)mF(m=2-10)体系的计算结果均表明, 体系总电子相关能与体系中(m-1)数值呈线性关系.     

CdH分子的量子化学研究

在相对论赝势从头算全有效空间多组态自洽场方法（CASSCF）基础上，对CdH分子进行了二级电子相关（SOCI）计算，得出此体系几个电子态的电子结构、势能曲线及光谱常数，认为，静态相关效应及动态相关效应对该体系的光谱性质等都有较大影响．     

基于块相关框架的多参考微扰理论和多参考耦合簇理论

传统单参考电子相关方法已经发展成熟,但很多时候无法正确描述共价键解离、双(多)自由基和激发态等电子之间相关性非常强的体系。近年来发展的多参考微扰理论和多参考耦合簇理论以多个行列式的线性组合为参考波函数,采用不同的方式有效考虑电子之间的动态相关,对强关联体系的描述取得了显著的改进。但根据理论出发点和精度要求的不同发展出了许多多参考理论,仍无一个公认的、令人满意的方案。本文将结合与常见电子相关方法在理论框架和计算精度上的比较,详细阐述块相关理论的基本原理,并介绍基于块相关的"另类"多参考电子相关方法。最后本文还简单展望了多参考电子相关方法今后的发展趋势。     

高能体系分子间相互作用研究: 含NNO~2和NH~2混合物

以abinitioHF/6-31G^*计算求得NH~3+NH~2NO~2的两种优化构型,经MP4电子相关能校正和Boys-Bernardi方案校正基组叠加误差求得精确的分子间相互作用能。还用PM3方法计算研究TATB(均三氨基三硝基苯)分别与HMX(奥克托金)和RDX(黑索金)的混合体系,经色散能校正电子相关近似地求得分子间相互作用能。结果表明,NH~3与NH~2NO~2之间的最大结合能为-38.32kJ/mol;分子间相互作用增强了N-NO~2键强度;TATB与HMX,RDX的结合能远大于石墨与HMX或RDX的结合能,表明TATB对HMX和RDX的润湿和钝感作用较石墨更强。     

考虑相关能效应的过渡态计算方法

提出了一种在X_a方法基础上,同时包含电子相关效应和考虑电子自相互作用的过渡态计算方法.由此计算的一系列原子电离势的结果表明,相关能对电离势的计算结果有很大影响,其数值为—0.41eV～0.94eV.引入相关能效应使计算结果明显趋于合理.此计算模型兼有简便和合理的特点,且能适用于分子体系的计算.     

DNA碱基的电子相关效应

分别在Hartree-Fock和二级Moller-Plesset微扰理论MP2的电子相关校正水平, 用6-31G^*^*基函数对四种DNA碱基胞嘧啶、胸腺嘧啶、鸟嘌呤和腺嘌呤的能量、偶极矩、电荷分布和分子静电势等性质进行了系统的从头计算研究。其中, 采用Z矢量方法在波函数中加入MP2级别的电子相关校正; 分别用Mulliken布居分析、静电势导出电荷CHELPG以及自然布居分析NPA方法计算分子中原子电荷。在上述计算结果的基础上, 系统地分析了DNA碱基的电子相关效应。     

Examination of numerical accuracy on fragment-DFT calculations with integral values of total electron density functions

Fragment molecular orbital (FMO) method gives a powerful tool to investigate electronic structures for biological substances, and ABINIT-MP program has been developed to implement ab initio FMO calculations effectively. We introduced DFT code into ABINIT-MP and applied fragment-DFT (F-DFT) calculations to model polypeptides. The total accuracy of numerical integrations employed in those calculations was examined by the total numbers of electrons in the molecules. It is shown that the numerical integral of the total density function under the fragment approximation works as an indicator for the numerically total accuracy on the F-DFT implementation.     

Development of a new density functional program for all-electron calculation of proteins

In this article, we propose a new molecular orbital program for all-electron calculation of proteins which is based on density functional theory. To carry it out in a fully analytical way, we adopted the (pure-) analytical Xα method and modified it for saving a lot of memories for large-scale calculations. The recent software technology sophisticated in information science is inevitably applied to achieve calculations of large molecular systems. The program is coded by the object-oriented language C + +, its output is shown graphically, and the most of the procedures in this program are controlled through an efficient graphical user interface developed by ourselves. Such technology supports the safe construction of the huge software, the tidy representation of enormous data, and the ready control of complex calculations. Test calculations with various sizes of glycine polypeptides indicate that the computation time is proportional to the 1.7 powers of the number of residues. This result suggests that the all-electron calculations of proteins consisting of over 1000 atoms could be performed with distributed and/or massively parallel computers. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 245–256, 1997     

A computational method for approximate evaluation of ionization constants from potentiometric data

A general equation is given for the evaluation of proton dissociation constants from potentiometric data. The method is applicable to monobasic acids, and also polybasic acids, where the successive dissociation steps may or may not overlap. A computer program was employed for the calculations. The method was tested by evaluation of the constants for 22 acids, from monobasic to octobasic, with satisfactory agreement with literature data.     

分子积分和分子对称性——分子局部对称性在abinitio程序中实现

对前人利用分子对称性简化分子积分计算之原理给予了严格的数学证明,并将利用分子局部对称性简化分子积分计算的方法在MQAB 80abinitio程序中实现.实算算例表明,最多可节省计算的CPU时间50%以上。该方法具有简便、有效之特点.     

Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide

Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree–Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by “breaking” bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. Results from tests of the method to find relative conformational energies and geometries of alanine dipeptides and alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1468–1494, 1999     

核磁AA'BB'氢谱体系的图谱解析程序

核磁AA′BB′氢谱体系图谱全部分辨时共有28条谱线,相互关系复杂,计算繁难。借助群论可大大简化计算,且便于识别谱线归属和获取谱参数,特别适于计算机解析。基本原理AA′BB′氢谱体系可组成16个初始基函数,其构型特征具有二重对称性。     

CASPER: a computer program used for structural analysis of carbohydrates.

The computer program CASPER and its algorithms are described. The program is aimed at facilitating the determination of structures of oligosaccharides and regular polysaccharides, requiring as input either the one-dimensional 1H or 13C NMR spectrum or the 2D C,H-correlation NMR spectrum together with information on components and linkages. The databases, the method of simulating spectra, options of the program, and techniques for faster calculations are described as well as an example of a structural determination.     

Graphical representation of a new algorithm for nonorthogonalab initio valence bond calculations

A pictorial representation of the algorithm using successive expansion method for the nonorthogonal VB calculations is given. With the help of this representation and the graph analysis, the efficiency of this algorithm is improved and theN! problem is reduced by a factor of about (N!)^1/2. Anab initio VB program for valence bond self-consistent-field (VBSCF) calculations has been implemented based on this algorithm. Some VBSCF calculations have been performed for systems of up to 14 electrons. The statistics of the CPU time of the calculations indicate that this new group-theoretical approach is quite practical.     

Simple evaluation of Franck–Condon factors and non‐Condon effects in the Morse potential

The calculation of Franck–Condon factors between different 1‐D Morse potential eigenstates using a formula derived from the Wigner function is discussed. Our numerical calculations using a simple program written in Mathematica are compared with other calculations. We show that our results have a similar accuracy as those calculations performed with more sophisticated methods. We discuss the extension of our method to include non‐Condon effects in the calculation. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem 88: 280–295, 2002     

The MMP2 calculational method

The molecular mechanics (MMP2) program and procedures for the treatment of conjugated hydrocarbons, and some of the results which they can achieve are described. The program is an updated version of the similar MMP1 program, but contains some differences. It is based on an SCF π system calculation, rather than on the VESCF method used earlier. All parameters are compatible with those in the MM2 program. Hence it is possible to calculate heats of formation, resonance energies, and structures for conjugated hydrocarbons in a way that is consistent with the calculations on non-conjugated molecules. The overall results as far as structure and energy are somewhat better than they were with the MMP1 program.