孪函数多电子基组在LiH分子方面的应用

选择从头算ＳＴＯ－６Ｇ基组，以孪函数Ｎ电子基组计算ＬｉＨ分子的基态和激发态。结果表明，无需全部孪函数Ｎ电子基组，只要适当选择少数基组即可达到所需精度。     

求解弱相互作用分子A—BC的振转能级和波函数的SCF—CI方法

给出了一种计算Ａ－ＢＣ型弱相互作用分子体系的分子间振转态的自洽场－组态相互作用（ＳＣＦ－ＣＩ）方法，首先使用ＳＣＦ方法优化径向伸缩振动和弯曲振动的基函数，再用ＣＩ方法确定精确的振转激发发态能级，在求解一维的伸缩振动了ＳＣＦ方程时，采用Ｎｕｍｅｒｏｖ－Ｊｏｈｎｓｏｎ算法在给定区间上求解获得振动波函数的数值解，具体计算了Ａｒ－ＨＣｌ和Ａｒ－Ｎ２体系的振动激发态的能级以验证该方法，结果表明，本方法可用较     

非线性光学极化率密度泛函理论计算的基组效应

由于分子的非线性光学性质与分子外层电子行为及激发性质密切相关，扩散函数对分子的非线性光学极化率计算非常重要．在ADF程序中的DZP基组基础上缀加扩散函数，构造出我们称之为DZP df的新基组．通过对5个模型分子的含频二阶非线性光学极化率的密度泛函理论计算，表明新基组可以得出较DZP基组更为准确可靠的结果，同时比较ADF程序内置带有扩散函数的大基组，计算量大为减少．     

层间化合物hNB—Cu稳定性的理论研究

以分子簇为模型，用从头算方法，在６－３１Ｇ＾＊基组水平上对六方氮化硼（ｈＢＮ）进行几何全优化，再用优化构型对Ｃｕ的六方氮化硼层间化合物（ｈＮＢ－Ｃｕ）做单点计算，根据计算结果，从层间距、Ｍｕｌｌｉｋｅｎ布居、轨道相互作用、原子净电荷及前线轨道等角度分析了ｈＢＮ－Ｃｕ的电子结构及其稳定性，通过劝分析计算确认了这一层间化合物只是一个介稳结构，阐明了Ｃｕ使六方氮化硼陶瓷致密化的作用机理以及ｈＮＢ－Ｃｕ的弱抗氧化性。     

OClO里德堡态激发能的准确预测及其阴离子低能激发态的从头算研究

采用全活化空间自洽场方法(CASSCF)研究了OClO阴离子7个低能电子态及其自由基的基态. 为了进一步考虑动态电子相关效应, 采用二级多组态微扰理论(CASPT2)获得更加可靠的能量值. 此外, 在ANO-L基组的基础上, 在OClO自由基的电荷中心增加了为研究里德堡态所建立的1s1p1d的波函数, 并应用多组态二级微扰理论(MS-CASPT2)方法获得了里德堡态的准确电子激发能.     

OClO里德堡态激发能的准确预测及其阴离子低能激发态的从头算研究

采用全活化空间自洽场方法(CASSCF)研究了OClO阴离子7个低能电子态及其自由基的基态.为了进一步考虑动态电子相关效应,采用二级多组态微扰理论(CASPT2)获得更加可靠的能量值.此外,在ANO-L基组的基础上,在OClO自由基的电荷中心增加了为研究里德堡态所建立的1s1p1d的波函数,并应用多组态二级微扰理论(MS-CASPT2)方法获得了里德堡态的准确电子激发能.     

HCB和HCB-的激发态光谱与态-态转化的理论研究

HCB是最简单的有机硼化合物，理论上已有少量研究，报导了部分激发态及其谐振动频率．关于HCB-及其异构体HBC-的研究未见报导．本文采用精确的量子化学计算方法，对它们的激发态性质、光谱和态一态转化进行了系统的研究，为实验上对这些物质的研究提供了一些有用的信息．计算采用极化的三Zeta基组（TZP）进行基态和激发态的构型优化与振动性质调查各电子态相关能估算应用基于多组态自洽场分子轨道的一级组态相互作用（FOCI／MCSCF）方法．计算过程中，考虑到对称性破裂（symmetrybrealling）的影响同，对存在对称性破裂的体系采用…     

建立在群对称轨道(SMO)及特征组态(CD)基础上的线性分子SMO-CDCI方法

定义和讨论了线性分子的群对称轨道（SMO）及特征组态函数（CD），它能够对线性分子的分子轨道与组态函数进行简单而完全的对称分类．SMO－CDCI方法是一种高效的计算方法，大大节省了线性分子CI计算的机时与内存．     

簇合物Mo2X4(X＝S,O)电子结构和光谱性质的Ab Initio研究

在ａｂ ｉｎｉｔｉｏ水平上对含ｃｉｓ／ｔｒａｎｓ－Ｍｏ２Ｓ４核和含Ｍｏ２Ｏ４核簇合物的电子结构进行了研究。对钼原子选取和构造了一组适合于含Ｍｏ－Ｍｏ金属键的双核钼簇合物从头算的基组，利用该基组并结合自然键轨道方法，对含上述簇胳的３个簇合物的电子结构进行了研究，并与相应的钨族合物进行了比较。     

完全活性空间组态相互作用能量的拟合和外推

完全活性空间组态相互作用计算与完全活性空间中的活性电子数和活性轨道数有关,但完全活性空间组态相互作用的能量不是活性电子数和活性轨道数的单调递减函数,因此活性轨道数和活性电子数不能用来外推完全活性空间组态相互作用的能量。为此,我们定义了一个新的变量:活性空间中的最大未占满轨道数。我们对一系列单重态、双重态和三重态分子进行了完全活性空间组态相互作用的计算,并利用活性空间中的活性电子数和最大未占满轨道数这两个变量,对这些基态能量进行了拟合和外推,拟合的均方根误差都在10~(-6)数量级。外推能量的精度优于MP4,对小分子体系,其精度高于CCSD。外推的完全的组态相互作用(FCI)能量值和实际计算的FCI值也很接近。另外,我们还利用外推能量来优化双原子分子的平衡键长,并计算谐振频率,其精度优于CASSCF。     

孪函数N电子基向Slater多电子反对称基的展开

发展了一套孪函数多电子基向Ｓｌａｔｅｒ行列式的展开方法，并得到展开系数的解析表达式，使得孪函数多电子基直接与Ｈａｒｔｅｒｒ－Ｆｏｃｋ从头计算相关联。     

Gaussian geminal basis set optimization with crude SCF reference state

Gaussian geminal basis functions for second-order correlation energy calculations according to the Sinanogˇlu method are optimized with reference to rather crude SCF functions. The optimized geminal basis set is then used in a one-step calculation of the correlation energy with respect to the near-Hartree-Fock reference State. The numerical results for the beryllium atom indicate the usefulness of the proposed technique.     

Electron correlation described by extended geminal models: The EXGEM2 and EXGEM3 models

The general extended geminal model is reviewed and two new approximate models EXGEM 2 and EXGEM 3 are introduced. In a test calculation on water using a double-zeta basis set, the approximate models recover 93.4% of the full CI correlation energy defined within the same basis set. A test calculation on the neon atom demonstrates that the performance of the models in comparison with the full CI , will improve as the basis set is increased. It is suggested that for basis sets of moderate and large size, and which include polarization functions, the extended geminal models are likely to recover 95%–97% of the correlation energy obtainable by the full CI .     

Variational formulation of perturbative explicitly-correlated coupled-cluster methods

We present a variational formulation of the recently-proposed CCSD(2)(R12) method [Valeev, Phys. Chem. Chem. Phys., 2008, 10, 106]. The centerpiece of this approach is the CCSD(2)(R12) Lagrangian obtained via L?wdin partitioning of the coupled-cluster singles and doubles (CCSD) Hamiltonian. Extremization of the Lagrangian yields the second-order basis set incompleteness correction for the CCSD energy. We also developed a simpler Hylleraas-type functional that only depends on one set of geminal amplitudes by applying screening approximations. This functional is used to develop a diagonal orbital-invariant version of the method in which the geminal amplitudes are fixed at the values determined by the first-order cusp conditions. Extension of the variational method to include perturbatively the effect of connected triples produces the method that approximates the complete basis-set limit of the standard CCSD plus perturbative triples [CCSD(T)] method. For a set of 20 small closed-shell molecules, the method recovered at least 94.5/97.3% of the CBS CCSD(T) correlation energy with the aug-cc-pVDZ/aug-cc-pVTZ orbital basis set. For 12 isogyric reactions involving these molecules, combining the aug-cc-pVTZ correlation energies with the aug-cc-pVQZ Hartree-Fock energies produces the electronic reaction energies with a mean absolute deviation of 1.4 kJ mol(-1) from the experimental values. The method has the same number of optimized parameters as the corresponding CCSD(T) model, does not require any modification of the coupled-cluster computer program, and only needs a small triple-zeta basis to match the precision of the considerably more expensive standard quintuple-zeta CCSD(T) computation.     

Quintuple-zeta quality coupled-cluster correlation energies with triple-zeta basis sets

The explicitly-correlated coupled-cluster method CCSD(T)(R12) is extended to include F12 geminal basis functions that decay exponentially with the interelectronic distance and reproduce the form of the average Coulomb hole more accurately than linear-r(12). Equations derived using the Ansatz 2 strong orthogonality projector are presented. The convergence of the correlation energy with orbital basis set for the new CCSD(T)(F12) method is studied and found to be rapid, 98% of the basis set limit correlation energy is typically recovered using triple-zeta orbital basis sets. The performance for reaction enthalpies is assessed via a test set of 15 reactions involving 23 molecules. The title statement is found to hold equally true for total and relative correlation energies.     

Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model

To approach the complete basis set limit of the "gold-standard" coupled-cluster singles and doubles plus perturbative triples [CCSD(T)] method, we extend the recently proposed perturbative explicitly correlated coupled-cluster singles and doubles method, CCSD(2)(R12) [E. F. Valeev, Phys. Chem. Chem. Phys. 8, 106 (2008)], to account for the effect of connected three-electron correlations. The natural choice of the zeroth-order Hamiltonian produces a perturbation expansion with rigorously separable second-order energy corrections due to the explicitly correlated geminals and conventional triple and higher excitations. The resulting CCSD(T)(R12) energy is defined as a sum of the standard CCSD(T) energy and an amplitude-dependent geminal correction. The method is technically very simple: Its implementation requires no modification of the standard CCSD(T) program and the formal cost of the geminal correction is small. We investigate the performance of the open-shell version of the CCSD(T)(R12) method as a possible replacement of the standard complete-basis-set CCSD(T) energies in the high accuracy extrapolated ab initio thermochemistry model of Stanton et al. [J. Chem. Phys. 121, 11599 (2004)]. Correlation contributions to the heat of formation computed with the new method in an aug-cc-pCVXZ basis set have mean absolute basis set errors of 2.8 and 1.0 kJmol when X is T and Q, respectively. The corresponding errors of the standard CCSD(T) method are 9.1, 4.0, and 2.1 kJmol when X=T, Q, and 5. Simple two-point basis set extrapolations of standard CCSD(T) energies perform better than the explicitly correlated method for absolute correlation energies and atomization energies, but no such advantage found when computing heats of formation. A simple Schwenke-type two-point extrapolation of the CCSD(T)(R12)aug-cc-pCVXZ energies with X=T,Q yields the most accurate heats of formation found in this work, in error on average by 0.5 kJmol and at most by 1.7 kJmol.     

Inclusion of hydrogen p orbitals in the semiempirical calculation of NMR parameters. I. Influence on the FPT INDO spin–spin coupling constants

Hydrogen 2p orbitals have been introduced into the basis set to calculate the Fermi contact term of spin–spin coupling constants using the FPT INDO method. Different coupling constants show different sensitivity to these hydrogen polarization functions. Some improvements are found for molecules containing N or F. Calculations of proton-proton geminal coupling constants give more negative results than those of FPT INDO , yielding a better agreement with experimental values. The π-transmission mechanism is notably exaggerated.     

Communication: Second-order multireference perturbation theory with explicit correlation: CASPT2-F12

An explicitly correlated complete active space second-order perturbation (CASPT2-F12) method is presented which strongly accelerates the convergence of CASPT2 energies and properties with respect to the basis set size. A Slater-type geminal function is employed as a correlation factor to represent the electron-electron cusp of the wave function. The explicitly correlated terms in the wave function are internally contracted. The required density matrix elements and coupling coefficients are the same as in conventional CASPT2, and the additional computational effort for the F12 correction is small. The CASPT2-F12 method is applied to the singlet-triplet splitting of methylene, the dissociation energy of ozone, and low-lying excited states of pyrrole.