Direct CI method in restricted configuration spaces

A configuration interaction (CI) method in restricted CI (RCI) space obtained by imposing inequality-type restrictions on the occupancies of groups of molecular orbitals (MOs) was studied. The direct CI approach in such space was analyzed, and some recommendations concerning practical implementation of the RCI method are given. The corresponding program has been written in FORTRAN 77 for an IBM 486 DX personal computer and has been used for electronic structure calculations on transition metal complexes using a valence MO basis with the INDO approximation. © 1996 by John Wiley & Sons, Inc.     

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients B(K) of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, Delta E(K) = (E-H(KK))B(K)2/(1-B(K)2), with B(K) determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Delta E(dis), is approximated by the sum of Delta E(K)s of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound E(S) is computed by CI in a selected space S, and E(M) = E(S) + Delta E(dis) + delta E, where delta E is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm(-1)) is achieved in a model space M of 1.4 x 10(9) CSFs (1.1 x 10(12) determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 x 10(12) CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since Delta E(dis) can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of E(S) is taken up in a companion paper.     

Fast calculation of excited-state potentials for rare-gas diatomic molecules: Ne2 and Ar2

A fast method for obtaining excited-state potentials of rare-gas diatomic molecules is described. Two types of excited orbitals are used: molecular orbitals calculated in the field of a singly charged molecular ion, and atomic orbitals (properly symmetrized) obtained in a similar atomic system. The RPA equations are solved within the manifold of excitations from the highest occupied orbital in each symmetry to the lowest excited orbital of either type in each symmetry. A simple model for estimating the dynamic correlation correction to excitation (and ionization) energies is given. Applications to excited states of Ne₂(^1,3Σ⁺_{g, u}, ^1,3Π_{g, u}) and Ar₂(^1,3Σ⁺_{g, u}) are described. Two-electron integral transformations involve only three orbitals of each symmetry, and the RPA matrices are four-dimensional. The computational effort required for all excited-state potentials adds less than one-tenths (in terms of computer time) to the effort involved in the preliminary ground state Hartree—Fock calculations. The resulting potentials compare favorably with more elaborate CI calculations and give good agreement with spectroscopic and scattering data. Potential curves for the molecular ions are also given.     

给体-桥-受体型系列C~6~0吡咯/二茂铁的电子结构及 二阶非线性光学性质的理 论研究

用INDO/2和INDO/CI方法,计算了系列给体-桥-受体(D-B-A)型C~6~0吡咯/二茂铁(C~6~0PY/Fc)的结构和电子光谱,计算结果和实验结果一致,在正确的光谱基础上,用INDO/CI-SOS公式计算了该系列分子的二阶非线性光学系数β~i~j~k,β~μ。考察了分子电子结构对β~μ影响的微观本质。     

An ab initio study of the ground and excited states of p-benzoquinone

The results of anab initio SCF calculation for the ground state and CI calculations for the excited states of p-benzoquinone are presented and discussed. A minimum basis set of Slater type orbitals was employed and the CI calculations were performed by considering single excitations from valence to virtual SCF molecular orbitals. The convergence of the calculated excitation energies is studied as a function of the number of orbitals used in the CI calculations. These calculations explain quite well the experimental results.     

Role of doubly excited states in calculations of optical activity

It is shown that doubly excited states play an important role in calculations of the optical activity of molecules with well-conjugated electron systems, such as the DNA bases. In some significant cases it is necessary to include a large number of excited states in the configuration interaction (CI ) to obtain a reliable, converging result. A new version of the CNDO/OPTIC method, which includes doubly excited states in the CI , is proposed. As an application, the electric transition moments in different pyrimidines are considered. The calculated results agree with experimental data and results obtained from ab initio calculations and INDO calculations using doubly excited states in the CI .     

Determination of one-centre core integrals from the average energies of atomic configurations

One-center core integrals for valence orbitals are determined from the experimental average energies of neutral atomic configurations from Li through Zn. These values are compared with those estimated from CNDO /1, “INDO /1”, CNDO /2, “INDO /2” and with theoretical values calculated from a pseudo-potential method. The agreement is good between values obtained from neutral atoms and from the psuedo-potential calculation except for the 3d orbitals of the transition elements where the theoretically calculated integrals over single ξ functions are not realistic. These two methods reproduce both term and average configuration energies for the first two rows of atoms; the semiempirical method reliably reproduces them for the third row. The CNDO /1 and INDO /1 methods underestimate atomic energies, while the CNDO /2 and INDO /2 procedures fail rather poorly. The propriety of using core integrals estimated semiempirically in molecular orbital calculations is discussed.     

Complete active-space configuration interaction with optimized orbitals: Application to Li2

We performed CAS –CI calculations on Li₂ using a set of molecular orbitals (MO ) optimized with a procedure that, in the case of highly symmetric molecules, permits extraction of a small set of MO out of a large set of atomic orbitals (AO ). The dimension of the CAS –CI space was of about 12 million symmetry-adapted determinants. We determined some spectroscopic constants of Li₂ with three different atomic basis sets of increasing quality. The values obtained with the largest atomic basis set are very close to the experimental results. © 1995 John Wiley & Sons, Inc.     

价键理论中的组态相互作用

近 2 0年来 ,从头计算水平的价键 (VB)方法得到了人们的重视 ,并广泛应用于化学反应等问题的研究[1～ 5] ,然而目前价键理论的计算方法仍然很不完善 .用 VBSCF方法进行计算虽然比较简单 ,能正确地描述化学反应的形成机理 ,但数值结果不理想 ;而用 BOVB方法[4 ] 进行计算虽然可以得到较好的计算结果 ,但存在收敛困难等问题 .分子轨道理论中的组态相互作用是一种简单直接的电子相关能计算方法 ,显然这一方法可以应用于价键方法中 .然而与分子轨道理论方法不同 ,在价键方法中 ,无法直接得到空轨道 ,此外如何选取激发价键函数使得计算结果…     

An analysis of the through-bond interaction using the localized molecular orbitals—II : Long-range proton hyperfine coupling in bridgehead alkyl radicals: bicyclo[1.1.1]pent-1-yl,bicyclo[2.1.1]hex-1-yl and bicyclo[2.2.1]hept-1-yl radicals

The INDO calculations were performed on three bridgehead alkyl radicals; bicyclo[1.1.1]pent-1-yl, bicyclo[2.1.1]hex-l-yl and bicyclo[2.2.1]hept-1-yl radicals. We have transformed the canonical molecular orbitals obtained by the INDO method into the localized molecular orbitals. With the use of the obtained localized molecular orbitals, the variation in the hyperfine coupling constant at the bridgehead proton in these radicals was pursued in terms of the through-bond (and/or the through-space) interaction according to the method by which we selectively can pick up a particular interaction between the specified localized molecular orbitals in a radical. As a result of this analysis, it was found that the hyperfine coupling constants in these radicals can be expressed by the summation of several terms; through-virtuals, through-space, through-bond, and some other coupling terms.     

Ab-initio calculations on large molecules and solids by desirable computational procedures

Desirable computational procedures developed here recently for ab-initio calculations on large molecules are outlined. Effective core model potentials (MODPOT) permit calculations of valence electrons only explicitly, yet accurately; a charge-conserving integral prescreening evaluation to decide whether a block of integrals will be larger than a preset threshold and thus be calculated explicitly is effective for spatially extended systems; an efficient MERGE technique to save and reuse common invariant skeletal integrals is useful for geometry variations and for adding basis fcuntions, substituent groups and molecules; and an effective configuration interaction (CI) Hamiltonian into which are folded the effects of the occupied molecular orbitals from which no excitations are allowed is useful for molecular decompositions and intermolecular reactions. These techniques have been extended for CI calculations on breaking a chemical bond in a molecule in a crystal or solid; atom-class/atomic-class potential functions and dispersion calculations have been added. In a new program, POLY-CRYST, all the integral strategies for large molecules are meshed.     

Basis set selections for relativistic self-consistent field calculations

Practical methods of generating reliable and economic basis sets for relativistic self-consistent fields (RSCF) calculations are developed. Large component basis sets are generated from constrained optimizations of exponents in the nonrelativistic atomic calculations for light atoms. For heavy atoms, large component basis sets for inner core orbitals are generated by fitting numerical atomic spinors of Dirac-Hartree-Fock calculations with appropriate number of Slater-type functions. Small component basis sets are obtained by using the kinetic balance condition and other computational criteria. With judicious selections of the basis sets, virtual orbitals in RSCF calculations become very similar to those in nonrelativistic calculations, implying that relativistic virtual orbitals can be used in electron correlation calculations in the same manner as the conventional nonrelativistic virtual orbitals. It is also evident that the Koopmans' theorem is also valid in RSCF results.     

An efficient atomic orbital based second-order Møller-Plesset gradient program

Based on the orbital-invariant atomic orbital formulation of the MP2 (M?ller-Plesset second-order perturbation theory) energy and gradient [P. Pulay and S. Saeb?, Theor. Chim. Acta 69, 357 (1986)], we have derived and programmed detailed working equations for closed-shell MP2 gradients. The orbital-invariant form avoids the difficulties of other formulations with frozen orbitals, and allows the use of arbitrary occupied orbitals, an important consideration for local correlation theories, although the present program uses canonical molecular orbitals. The atomic orbital formulation offers savings both in storage and computer time. Test calculations on systems containing up to approximately 100 atoms and approximately 1000 basis functions, performed on a single personal computer, are reported. Parallelization of the code is underway.     

Study of localized molecular orbitals using group theory methods and its approach to the many-electron correlation problem. IV. The symmetry-adaptation of many-center integrals and hamiltonian matrix elements in MCSCF calculations

Group theoretic methods are presented for the transformations of integrals and the evaluation of matrix elements encountered in multiconfigurational self-consistent field (MCSCF) and configuration interaction (CI) calculations. The method has the advantages of needing only to deal with a symmetry unique set of atomic orbitals (AO) integrals and transformation from unique atomic integrals to unique molecular integrals rather than with all of them. Hamiltonian matrix element is expressed by a linear combination of product terms of many-center unique integrals and geometric factors. The group symmetry localized orbitals as atomic and molecular orbitals are a key feature of this algorithm. The method provides an alternative to traditional method that requires a table of coupling coefficients for products of the irreducible representations of the molecular point group. Geometric factors effectively eliminate these coupling coefficients. The saving of time and space in integral computations and transformations is analyzed. © 1994 by John Wiley & Sons, Inc.     

A new fragment-based approach for calculating electronic excitation energies of large systems

We present a new fragment-based scheme to calculate the excited states of large systems without necessity of a Hartree-Fock (HF) solution of the whole system. This method is based on the implementation of the renormalized excitonic method [M. A. Hajj et al., Phys. Rev. B 72, 224412 (2005)] at ab initio level, which assumes that the excitation of the whole system can be expressed by a linear combination of various local excitations. We decomposed the whole system into several blocks and then constructed the effective Hamiltonians for the intra- and inter-block interactions with block canonical molecular orbitals instead of widely used localized molecular orbitals. Accordingly, we avoided the prerequisite HF solution and the localization procedure of the molecular orbitals in the popular local correlation methods. Test calculations were implemented for hydrogen molecule chains at the full configuration interaction, symmetry adapted cluster/symmetry adapted cluster configuration interaction, HF/configuration interaction singles (CIS) levels and more realistic polyene systems at the HF/CIS level. The calculated vertical excitation energies for lowest excited states are in reasonable accordance with those determined by the calculations of the whole systems with traditional methods, showing that our new fragment-based method can give good estimates for low-lying energy spectra of both weak and moderate interaction systems with economic computational costs.     

无机硼酸盐体系的二阶非线性光学晶体材料

总结了无机硼酸盐二阶非线性光学晶体材料的结构特征，即无中心对称性和阴离子基团的π共轭类芳香性结构；讨论了阴离子基团理论处理实际体系时遇到的困难及其假设中的不尽合理之处。以BBO（低温相偏硼酸钡β-BaB2O4）和LBO（三硼酸或LiB3O5）晶体为例，从含组态作用的INDO/S记结合态求和方法计算得到的结果分析后表明，BBO和LBO的价带具有相同的原子轨道组成，导带的底部，在BBO中主要是Ba原子轨道的贡献；在LBO中主要是（B3O5）-基团的贡献。前者对二阶极化率的贡献主要来自环外O原子的2p轨道到Ba原子的电荷转移，而后者的贡献则主要来自[B－O]阴离子基团内的电荷转移．计算得到的BBO和LBO的价带和导带间的电子跃迁能分别是6.28eV和7.26eV，相当好地接近于紫外吸收边的测量值6．43eV和7.78eV，最大的二阶非线性极化系数Xyyy（7．18X10-9esu）和Xzyy（2.38X10－9esu）相当接近实验测量值d22（5．30X10－9esu）和d32（2．79X10－9esu）。     

Ab initio CI calculations on benzene with an extended basis set

An extended basis set of triple zeta plus polarization quality is employed to carry out configuration interaction (CI ) calculations of the three lowest singlet and triplet excited states of benzene. The CI calculation is carried out by taking into account single and double excitations of π and σ electrons. In the CI , composite natural orbitals (CNO s), which are constructed from the natural orbitals of the ground state of ethylene, are used as virtual orbitals. The aim of using CNOs is to reduce the number of virtual orbitals to be used in constructing configuration-state functions, thus cutting down CI dimensions without losing reasonable accuracy. The excitation energies resulting from the CI are in fairly good agreement with experiment. The root mean square of the deviation is 0.22 eV for the six calculated energies and the largest disagreement is 0.37 eV for the third singlet excited state. To obtain better excitation energies by an ab initio calculation, it seems likely that we need to take into account more electron correlation than in the present calculation. © 1994 John Wiley & Sons, Inc.     

Modified all-valence INDO/spd method for ground and excited state properties of isolated molecules and molecular complexes

A new semiempirical all-valence method, GRINDOL (Ghost and Rydberg INDO ), based on the INDO approximation, is described. Linderberg–Seamans relation (extended to the d and Rydberg orbitals) for the resonance integrals and a new semitheoretical expression for the core-core repulsion term and energy correction including basis-set superposition error (intermolecular as well as intramolecular) has been applied. The proposed method enables calculation of ground and excited state properties. The ground state results (including intermolecular interactions) as well as the spectral properties are in reasonable agreement with relevant experimental (or ab initio) studies for isolated molecules, molecular complexes, and transition metal compounds. The method contains only one adjustable parameter, all two-center integrals and terms are only basis-set dependent. The one-center integrals are evaluated from the respective atomic terms.     

The modified all-valence INDO method with the inclusion of spin–orbit coupling

The all-valence INDO method has been modified for the inclusion of spin–orbit coupling effects. In the method presented, the Hamiltonian includes spin–orbit coupling and the basis set constitutes the singlet and triplet determinental wave functions constructed from molecular orbitals resulting from nonrelativistic calculations. Eigenvectors obtained are later used for the evaluation of transition probabilities among different states. The results presented include lifetimes of different states of organic molecules and transition energies for halogen molecules and they are in a good agreement with experimental results. © 1992 John Wiley & Sons, Inc.     

Deorthogonalization of atomic originals in the CNDO approach

The procedure for deorthogonalization (D) of atomic orbitals in the semiempirical CNDO approach is reviewed. For comparative studies, CNDO/2, CNDO/2D, and STO -3G calculations of molecular dipole moments and Mulliken populations are carried out on 35 prototype molecules containing H, C, N, O, and F atoms. The calculated values are assessed on the basis of how well they agree with experimental trends, chemical bonding theories, and ab initio molecular orbital (MO) values. Results of analyses indicate that the CNDO/2D values for dipole moments are in reasonable agreement with experimental values, and those for net atomic charges and electron populations bear greater resemblance to the ab initio (STO -3G and 6-31G**) values than the original CNDO/2 values. These findings, together with those of previous investigators, demonstrate unequivocally the advantages of incorporating deorthogonalization into routine CNDO/2 or INDO calculations as a means to obtain reasonable estimates of charge distributions.