价键理论新进展

概要介绍了现代价键理论的几个主要方法,并讨论了它们各自的特点及其发展现状,并重点介绍了键表方法的基本理论、计算程序及一些应用。     

A programmable spin-free method for configuration interaction

In this note a method is presented for quick implementation of configuration interaction (CI) calculations in molecules. A spin-free Hamiltonian for anN electron system in a spin stateS, expressed in terms of the generators for the unitary group algebra, is diagonalized over orbital configurations forming a basis for the irreducible representation [2^1/2N-S 1^2S ] of the permutation group S _N . It has been found that the basic algebraic expressions necessary for the CI calculation involve a limited category of permutations. These have been displayed explicitly. On leave from the Indian Institute of Technology, Bombay, India.     

Simultaneous determination of nuclear and electronic wave functions without Born–Oppenheimer approximation: Ab initio NO+MO/HF theory

We develop a simultaneous determination method of nuclear and electronic wave functions without the Born–Oppenheimer approximation. We examine two expanding methods, namely, molecular orbital (MO)-type and valence bond (VB)-type expansions for a nuclear orbital, which is a one-particle wave function of a nucleus. The VB-type expansion is shown to be more accurate than the MO-type one because of the local nature of the nuclei. We also investigate the basis function expansion of the nuclear orbital and propose a scheme to determine the orbital exponent for the nuclear basis function. Numerical calculations confirm the accuracy and feasibility of the present method. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

MOLCAS 7: The Next Generation

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two‐electron integrals and in the generation of so‐called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self‐consistent field, density functional theory, 2nd order perturbation theory, complete‐active space self‐consistent field multiconfigurational reference 2nd order perturbation theory, and coupled‐cluster methods. The report further elaborates on the implementation of a restricted‐active space self‐consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas‐Kroll‐Hess transformation for one‐component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so‐called picture‐change‐free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Application of cholesky-like matrix decomposition methods to the evaluation of atomic orbital integrals and integral derivatives

When viewed as a square two-indexed matrix, the array of atomic orbital-based, two-electron integrals (ij|kl) is a positive semidefinite array. Beebe and Linderberg showed, in 1977, that actual or near linear dependencies often exist within the types of atomic orbital basis sets employed in conventional quantum chemical calculations. In fact, large (i.e., higher quality) bases were shown to be substantially more redundant than smaller or more spatially separated bases. In situations where there exists significant basis near redundancy, the rank (r) of the (ij|kl) ≡ V_l,J matrix of integrals will be significantly smaller than the matrix dimension M. When this occurs, it proves computationally tractable to decompose the M-dimensional matrix V into components L ( V = LL ^T) which contain all of the information needed to form the full V matrix. The Cholesky algorithm allow such a decomposition to be carried out and forms the basis of the work described here. The method is found to be highly successful in reducing the number of integrals and integral derivatives that must actually be calculated. In particular, results on the C₂ molecule indicate that the algorithm can be superior to traditional methods of integral derivative generation if the orbital basis is large enough to contain appreciable near redundancy. In contrast, results on benzene with a more spatially delocalized basis show that conventional methods are preferred whenever substantial basis (near) redundancy is not present.     

Configuration interaction matrix elements for atoms using permutation group algebra

A procedure is described for the efficient evaluation of the energy matrix elements necessary for atomic configuration-interaction calculations. With the orbital configurations of an N electron system in spin state S written as the irreducible representations [2^1/2N?S, 1^2S] of the permutation group S( N ), it is possible to evaluate readily the energy matrix elements of a spin-free Hamiltonian expressed in terms of the generators of the unitary group. We show how the use of angular momentum ladder operators permits the effective generation of a basis of eigenstates of ??², ??_z as well as ??² and ??_z, for which the energy matrix elements may be evaluated with ease.     

A general,recursive, and open‐ended response code

We present a new implementation of a recent open‐ended response theory formulation for time‐ and perturbation‐dependent basis sets (Thorvaldsen et al., J. Chem. Phys. 2008, 129, 214108) at the Hartree–Fock and density functional levels of theory. A novel feature of the new implementation is the use of recursive programming techniques, making it possible to write highly compact code for the analytic calculation of any response property at any valid choice of rule for the order of perturbation at which to include perturbed density matrices. The formalism is expressed in terms of the density matrix in the atomic orbital basis, allowing the recursive scheme presented here to be used in linear‐scaling formulations of response theory as well as with two‐ and four‐component relativistic wave functions. To demonstrate the new code, we present calculations of the third geometrical derivatives of the frequency‐dependent second hyperpolarizability for HSOH at the Hartree–Fock level of theory, a seventh‐order energy derivative involving basis sets that are both time and perturbation dependent. © 2014 Wiley Periodicals, Inc.     

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.     

Implementation of molecular symmetry in valence bond calculation

A novel strategy for the construction of many-electron symmetry-adapted wave function is proposed for ab initio valence bond (VB) calculations and is implemented for valence bond self-consistent filed (VBSCF) and breathing orbital valence bond (BOVB) methods with various orbital optimization algorithms. Symmetry-adapted VB functions are constructed by the projection operator of symmetry group. The many-electron symmetry-adapted wave function is expressed in terms of symmetry-adapted VB functions, and thus the VB calculations can be performed with the molecular symmetry restriction. Test results show that molecular symmetry reduces the computational cost of both the iteration numbers and CPU time. Furthermore, excited states with specific symmetry can be conveniently obtained in VB calculations by using symmetry-adapted VB functions.     

Ab initio calculation of maximum bond order hybrid orbitals

Summary The maximum bond order hybrid orbital (MBOHO) procedure is tested onab initio level by use of the density matrix in Löwdin orthogonalized atomic orbital basis. The direct MBOHO calculation based on the whole density matrix includes also the hybridization of the inner atomic orbitals, and the MBOHO calculation based on the valence orbital part of the density matrix considers only the hybridization of the valence atomic orbitals. The concrete MBOHO calculations based on theab initio calculation with STO-3G basis show that the components of the s atomic orbitals in MBOHOs and the maximum bond orders obtained from the two kinds of MBOHO calculations are very close to each other, and that the two kinds of MBOHOs all have the excellent correlativity with the nuclear spin-spin coupling constants.The project supported by National Natural Science Foundation of China and the Excellent Young University Teacher's Foundation of State Education Commission of China.     

XMVB: a program for ab initio nonorthogonal valence bond computations

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.     

Implementation of a general multireference configuration interaction procedure with analytic gradients in a semiempirical context using the graphical unitary group approach

The graphical unitary group approach has been applied in an efficient implementation of a general multireference configuration interaction (MRCI) method for use with small active molecular orbital spaces in a semiempirical framework. Gradients can be computed analytically for molecular orbitals from a closed-shell or a half-electron open-shell Hartree-Fock calculation. CPU times for single point energy and gradient calculations are reported. The code allows MRCI geometry optimizations of large molecules, as illustrated for the singlet ground state and the four lowest triplet states of fullerene C(76).     

A charge-conserving approximation method for ab initio calculations

A new method is presented for approximate ab initio calculations in quantum chemistry. It is called CCAM (charge conserving approximation method). The calculation method does not include the use of empirical parameters. We use Slater type orbitals as basis set, replacing STO's by STO-2G functions to evaluate three- and four-center integrals and making the STO-2G two-orbital charge distributions have the same total charge as STO. The results are presented for test calculations on five molecules. In view of these results, CCAM is better than ab initio calculations over STO-6G in the results on total energies, kinetic energies and occupied orbital energies. In atomic populations, dipole moments and unoccupied orbital energies, CCAM is also satisfactory. We estimate that CCAM would be as fast as ab initio calculations over STO-2G in evaluating molecular integrals.     

Correlation corrected band structures of homopolypeptides v. B3LYP band structures of 19 homopolypeptides

Ab initio B3LYP crystal orbital (CO) calculations have been performed on the 19 homopolypeptides (PolyGly, PolyAla, PolySer, PolyThre, PolyLeu, PolyiLeu, PolyVal, PolyAspAc, PolyAsp, PolyGlutAc, PolyGlut, PolyHist, PolyProl, PolyCyst, PolyMeth, PolyTyr, PolyPhenAla, PolyArg, and PolyLys) in their β pleated sheet conformation. Keeping the main chain conformation fixed as in PolyGly, the side chain geometries were optimized. For the calculation 2n+1 different k points were used with n = 8 for the case of simpler and n = 10 for more complicated side chains. The basis set applied was the double ζ one of Clementi. According to the results obtained, the conduction bands are shifted upward and the valance band downward, compared with the results of previous BLYP 1 and LDA 7 CO calculations. The bandwidths are similar to the previous cases. The band edges are in many cases not at the endpoints of the first Brillouin zones, causing nonmonotonous dispersion of both the conduction bands (CB) and the valance bands (VB), respectively. The fundamental gap values due to the upward shifts of the CB and downward shifts of the VB are substantially larger than in the case of our previous DFT CO calculations (values 6.0–7.0 eV). They are very close to the gap values, which can be estimated on the basis of experimental ultraviolet (UV) spectra of some homopolypeptides and on the basis of intermediate exciton theoretical calculations (6.5–7.5 eV). These surprisingly good results for the gaps are due to the compensation of errors (LDA or BLYP gives too small and simple HF provides too large gap values) in the B3LYP method. The admixture of the exact HF exchange with a weight of 0.19 obviously compensates the self interaction error occurring in the LDA or BLYP methods. This article discusses whether/how this result could be established by other B3LYP CO calculations on simple polymer chains and on stacked systems (e.g., nucleotide base stacks). Furthermore, a comparative analysis of the ground state DFT methods, the HF method and of the optimized effective potential method could throw more light on our successful theoretical results for the gaps of the homopolypeptides. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

An efficient algorithm for complete active space valence bond self‐consistent field calculation

This article describes a novel algorithm for the optimization of valence bond self‐consistent field (VBSCF) wave function for a complete active space (CAS), so‐called VBSCF(CAS). This was achieved by applying the strategies adopted in the optimization of CASSCF wave functions to VBSCF(CAS) wave functions, using an auxiliary orthogonal orbital set that generates the same configuration space as the original nonorthogonal orbital set. Theoretical analyses and test calculations show that the VBSCF(CAS) method shares the same computational scaling as CASSCF. The test calculations show the current capability of VBSCF method, which involves millions of VB structures. © 2012 Wiley Periodicals, Inc.     

Variational treatment on the energy of the He‐sequence ground state with weakest bound electron potential model theory

The concept of spin–orbital of the weakest bound electron is described used to construct the antisymmetric wave function of atomic or ionic systems within weakest bound electron potential model theory (WBEPM theory). The total energies of He‐sequence (Z = 2–9) in the ground states is calculated with a variational method. The effect of fixed orbital approximation is discussed quantitatively. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

A systematic preparation of new contracted Gaussian-type orbital sets. III. Second-row atoms from Li through ne

Four minimal Gaussian basis sets are generated for the second-row atoms Li through Ne. The first one, MINI-1, consists of a 3-term contraction of primitive Gaussian-type orbitals for 1s, 2s, and 2p atomic orbitals. The convenient shorthand notation would be (3,3) for Li? Be and (3,3/3) for B? Ne. The second one, MINI-2, can be represented by (3,3/4) for B? Ne. In the same way, MINI-3 is described as (4,3) for Li? Be, and MINI-3 and MINI-4 are represented by (4,3/3) and (4,3/4) for B? Ne, respectively. Although the four basis sets are the minimal type, they give the valence shell orbital energies which are close to those of DZ. These four and other sets derived from them are tested for the hetero- and homodiatomic molecules and some organic molecules. They are found to give the orbital energies that agree well with those given by extended calculations. Atomization energies and other spectroscopic constants are also calculated and compared with those of extended calculations. The results clearly indicate that the present basis sets can be used very effectively in the molecular calculations.     

Comparison of atomic charges derived via different procedures

Atomic charges were obtained from ab initio molecular orbital calculations using a variety of procedures to compare them and assess their utility. Two procedures based on the molecular orbitals were examined, the Mulliken population analysis and the Weinhold–Reed Natural Population Analysis. Two procedures using the charge density distribution were included; the Hirshfeld procedure and Bader's Atoms in Molecules method. Charges also were derived by fitting the electrostatic potential (CHELPG) and making use of the atomic polar tensors (GAPT). The procedures were first examined for basis set independence, and then applied to a group of hydrocarbons. The dipole moments for these molecules were computed from the various atomic charges and compared to the total SCF dipole moments. This was followed by an examination of a series of substituted methanes, simple hydrides, and a group of typical organic compounds such as carbonyl derivatives, nitriles, and nitro compounds. In some cases, the ability of the charges to reproduce electrostatic potentials was examined. © John Wiley & Sons, Inc.     

Unitary group tensor operator algebras for many-electron systems. III. Matrix elements in U(n 1 +n 2) ⊃ U(n 1) × U(n 2) partitioned basis

Exploiting our earlier results [J. Math. Chem. 4 (1990) 295–353 and 13 (1993) 273–316] on the unitary group U(n) Racah-Wigner algebra, specifically designed for quantum chemical calculations of molecular electronic structure, and the related tensor operator formalism that enabled us to introduce spin-free orbital equivalents of the second quantization-like creation and annihilation operators as well as higher rank symmetric, antisymmetric and adjoint tensors, we consider the problem of U(n) basis partitioning that is required for group-function type approaches to the many-electron problem. Using the U(n) U(n ₁) × U(n ₂),n =n ₁ +n ₂ adapted basis, we evaluate all required matrix elements of U(n) generators and their products that arise in one- and two-body components of non-relativistic electronic Hamiltonians. The formalism employed naturally leads to a segmented form of these matrix elements, with many of the segments being identical to those of the standard unitary group approach. Relationship with similar approaches described earlier is briefly pointed out.     

VB2000: Pushing valence bond theory to new limits

Full valence bond (VB) calculations for a system of N electrons have always been hindered by the rapidly growing value of N!, which effectively imposes a limit N < 20. Often, however, not all electrons in a molecule are of interest; if we focus on a “group” G of N_G electrons (e.g., in an “active” region), then it is N_G! that sets the limit. In this work, group function (GF) theory is used to represent a molecule as a collection of interacting electron groups, each with a many‐electron wave function of any chosen form (e.g., VB, MO‐SCF, MCSCF), and each GF is optimized individually in a step‐by‐step process. An efficient VB algorithm allows for up to 14 electrons in any VB group and this combination of GF and VB methods greatly extends the range of feasibility of molecular calculations with VB‐type wave functions: Thus, (1) a large system can be divided into any number of smaller subsystems (groups); (2) each group may contain any chosen number of electrons; (3) the form of any group function (including its level of accuracy) may be chosen at will by the program user. A number of sample calculations are briefly presented. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002