The loge theory as a starting point for variational calculations. I. General formalism

It is shown that any expectation value of any observable associated with a molecule is the sum of loge contributions and of loge pair contributions. This result provides a rigorous theoretical basis for the study of additive properties of molecules. It is demonstrated that molecular wave functions (exact or approximate) can be expressed as a sum of functions corresponding to the various electronic events. Furthermore any of these event functions can be expressed in terms of correlated loge functions. This expression suggests many kinds of variational procedures of calculating wave functions (known methods and new ones). The case in which noncorrelated completely localized loge functions are used is discussed. If continuous functions are used the variational equation reduces to a sum of independent variational equations, each one corresponding to a particular electronic event. This is not so when discontinuous functions are used or when a delocalized function is added to replace the correlation interloge function. The noncorrelated completely localized loge model is analyzed in more detail. It is seen that local spin operators can be introduced and that each event density operator is the product of the loge density operators. Therefore that model is an independent loge model. The corresponding generalized self-consistent field equations are derived. This treatment helps us to understand how a localized state of a molecule can produce an ion containing a delocalized region, a phenomenon which is sometimes at the origin of some misunderstanding in photoelectron spectroscopy. Finally it is seen how virtual loge functions can be introduced to describe excited states.     

Variational principles for discontinuous wave functions and the independent particle model of electronic structure

Using the variational method advanced by McCavert and Rudge, we obtain the independent particle model variational principle for loge localized discontinuous wave functions. The transformation of the variational expressions into matrix form when the loge localized discontinuous orbitals are expanded in finite basis sets is discussed. The simplifications brought about by this new method in the evaluation of molecular integrals are indicated.     

NMR shielding tensors for density fitted local second-order Mo̸ller-Plesset perturbation theory using gauge including atomic orbitals

An efficient method for the calculation of nuclear magnetic resonance (NMR) shielding tensors is presented, which treats electron correlation at the level of second-order Mo?ller-Plesset perturbation theory. It uses spatially localized functions to span occupied and virtual molecular orbital spaces, respectively, which are expanded in a basis of gauge including atomic orbitals (GIAOs or London atomic orbitals). Doubly excited determinants are restricted to local subsets of the virtual space and pair energies with an interorbital distance beyond a certain threshold are omitted. Furthermore, density fitting is employed to factorize the electron repulsion integrals. Ordinary Gaussians are employed as fitting functions. It is shown that the errors in the resulting NMR shielding constant, introduced (i) by the local approximation and (ii) by density fitting, are very small or even negligible. The capabilities of the new program are demonstrated by calculations on some extended molecular systems, such as the cyclobutane pyrimidine dimer photolesion with adjacent nucleobases in the native intrahelical DNA double strand (ATTA sequence). Systems of that size were not accessible to correlated ab initio calculations of NMR spectra before. The presented method thus opens the door to new and interesting applications in this area.     

Formulation of the coupled cluster theory with localized orbitals in correlation calculations on polymers

It is shown that in the framework of coupled cluster theory both the correlation energy per unit cell and quasi-particle band structures of polymers can be computed directly from matrix elements of the excitation operator and the two-electron integrals calculated in localized orbital basis. Further, it is described how to thake advantage of the localized nature of the orbitals applied. Ab initio test calculations on a finite model system similar to the Pariser–Parr–Pople Hamiltonian are presented.     

On the use of spatial symmetry in ab initio calculations. Transformation of the two-electron integrals from atomic orbital to localized molecular orbital basis

Anm ⁵-dependent integral transformation procedure from atomic orbital basis to localized molecular orbitals is described for spatially extended systems with some Abelian symmetry groups. It is shown that exploiting spatial symmetry, the number of non-redundant integrals for normal saturated hydrocarbons can be reduced by a factor of 2.5-3.5, depending on the size of the system and on the basis. Starting from a list of integrals over basis functions in canonical order, the number of multiplications of the four-index transformation is reduced by a factor of 2.8-3.5 as compared to that of Diercksen's algorithm. It is pointed out that even larger reduction can be achieved if negligible integrals over localized molecular orbitals are omitted from the transformation in advance.     

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.     

Localizability of dynamic electron correlation

The convergence of the intrapair correlation energy for a localized internal orbital is investigated as the virtual subspace is enlarged. At variance with previous investigations of this kind, the virtual subspace is represented in atomic orbitals. This allows to define spatial relations between the orbitals involved. Typically, over 98% of the pair correlation energy is recovered by a small local basis set, consisting of the valence orbitals of the atoms with which the electron pair is associated. This opens the possibility of an efficient Cl procedure based on localized pairs.     

Solvation and ion-pairing properties of the aqueous sulfate anion: explicit versus effective electronic polarization

We assessed the relative merits of two approaches for including polarization effects in classical force fields for the sulfate anion. One of the approaches is the explicit shell model for atomic polarization and the other is an implicit dielectric continuum representation of the electronic polarization, wherein the polarizability density is spatially uniform. Both the solvation and ion association properties of sulfate were considered. We carried out an ab initio molecular dynamics simulation for a single sulfate anion in aqueous solution to obtain a benchmark for the solvation structure. For the ion-pairing properties, the models were compared to experimental thermodynamic data through Kirkwood-Buff theory, which relates the integrals of the pair correlation functions to measurable properties. While deficiencies were found for both of the approaches, the continuum polarization model was not systematically worse than the shell model. The shell model was found to give a more structured solution than the continuum polarization model, both with respect to solvation and ion pairing.     

Hartree–Fock MO–LCAO equations with charge-dependent atomic integrals

The Hartree–Fock equations are derived in the MO -LCAO approximation for the case when the integrals (except overlap integrals) over the atomic orbitals are charge-dependent. It is shown that inclusion of the overlap matrix in the iterative procedure gives equations which are too complicated for the simple model under consideration. The approach is applied to the VESCF method in the PPP scheme.     

Coupled-cluster theory in a projected atomic orbital basis

We present a biorthogonal formulation of coupled-cluster (CC) theory using a redundant projected atomic orbital (PAO) basis. The biorthogonal formulation provides simple equations, where the projectors involved in the definition of the PAO basis are absorbed in the integrals. Explicit expressions for the coupled-cluster singles and doubles equations are derived in the PAO basis. The PAO CC equations can be written in a form identical to the standard molecular orbital CC equations, only with integrals that are related to the atomic orbital integrals through different transformation matrices. The dependence of cluster amplitudes, integrals, and correlation energy contributions on the distance between the participating atomic centers and on the number of involved atomic centers is illustrated in numerical case studies. It is also discussed how the present reformulation of the CC equations opens new possibilities for reducing the number of involved parameters and thereby the computational cost.     

Calcul de l'Énergie de Corrélation pour l'État Fondamental de la Molécule d'Acétylène

Correlation effects in the acetylene molecule are studied by a method using the technique of Feynman-like diagrams. The introduction of the reduced atomic electronic integrals decreases considerably the value of correlation energy.     

The "JK-only" approximation in density matrix functional and wave function theory

Various energy functionals applying the "JK-only" approximation which leads to two-index two-electron integrals instead of four-index two-electron integrals in the electron-electron interaction term of the electronic energy are presented. Numerical results of multiconfiguration self-consistent field calculations for the best possible "JK-only" wave function are compared to those obtained from the pair excitation multiconfiguration self-consistent (PEMCSCF) method and two versions of density matrix functional theory. One of these is derived making explicit use of some necessary conditions for N representability of the second-order density matrix. It is shown that this method models the energy functional based on the best possible "JK-only" wave function with good accuracy. The calculations also indicate that only a minor fraction of the total correlation energy is incorporated by "JK-only" approaches for larger molecules.     

A systematic sequence of relativistic approximations

An approach to the development of a systematic sequence of relativistic approximations is reviewed. The approach depends on the atomically localized nature of relativistic effects, and is based on the normalized elimination of the small component in the matrix modified Dirac equation. Errors in the approximations are assessed relative to four-component Dirac-Hartree-Fock calculations or other reference points. Projection onto the positive energy states of the isolated atoms provides an approximation in which the energy-dependent parts of the matrices can be evaluated in separate atomic calculations and implemented in terms of two sets of contraction coefficients. The errors in this approximation are extremely small, of the order of 0.001 pm in bond lengths and tens of microhartrees in absolute energies. From this approximation it is possible to partition the atoms into relativistic and nonrelativistic groups and to treat the latter with the standard operators of nonrelativistic quantum mechanics. This partitioning is shared with the relativistic effective core potential approximation. For atoms in the second period, errors in the approximation are of the order of a few hundredths of a picometer in bond lengths and less than 1 kJ mol(-1) in dissociation energies; for atoms in the third period, errors are a few tenths of a picometer and a few kilojoule/mole, respectively. A third approximation for scalar relativistic effects replaces the relativistic two-electron integrals with the nonrelativistic integrals evaluated with the atomic Foldy-Wouthuysen coefficients as contraction coefficients. It is similar to the Douglas-Kroll-Hess approximation, and is accurate to about 0.1 pm and a few tenths of a kilojoule/mole. The integrals in all the approximations are no more complicated than the integrals in the full relativistic methods, and their derivatives are correspondingly easy to formulate and evaluate.     

Singular and nonsingular three-body integrals for exponential wave functions

Integrals which are individually singular, but which may be combined to yield convergent expressions, are needed for computations of relativistic effects and various properties of atomic and quasiatomic systems. As computations become more detailed and precise, more such integrals are required. This paper presents general formulas for the radial parts of the singular and nonsingular (regular) integrals that occur when three-body systems are described using wave functions that include exponentials in all three interparticle coordinates. Our results are compared with those found in the literature for some of the integrals, and are also shown to be consistent with previously reported results for Hylleraas functions (a limiting case in which one of the exponential parameters is set to zero).     

Quantum-chemical investigation of the structures and electronic spectra of the nucleic acid bases at the coupled cluster CC2 level

We study the ground-state structures and singlet- and triplet-excited states of the nucleic acid bases by applying the coupled cluster model CC2 in combination with a resolution-of-the-identity approximation for electron interaction integrals. Both basis set effects and the influence of dynamic electron correlation on the molecular structures are elucidated; the latter by comparing CC2 with Hartree-Fock and M?ller-Plesset perturbation theory to second order. Furthermore, we investigate basis set and electron correlation effects on the vertical excitation energies and compare our highest-level results with experiment and other theoretical approaches. It is shown that small basis sets are insufficient for obtaining accurate results for excited states of these molecules and that the CC2 approach to dynamic electron correlation is a reliable and efficient tool for electronic structure calculations on medium-sized molecules.     

Theory of correlation measurements of resonance light scattering

The boson nature of radiation is shown to give rise to a purely quantum mechanical exchange contribution to the intensity-intensity correlation function for resonant light scattering by an atomic or molecular system. The exchange contribution can be decomposed into three components, one involving the intensity correlation for a pair of coherently scattered photons (“resonant Raman” processes), another for a pair of incoherently scattered ones (“resonance fluorescence”), and the last involving the exchange correlation one of each. The intensity correlation measurements of Kimble et al., on optically pumped atomic beams of sodium atoms are interpreted with the theory, producing values of the decay rate of the excited sodium atoms and of the coherence time of the exciting radiation in good agreement with expectations.     

Ab initio program for treatment of related molecules. II. Integral transformation in extended systems

An integral transformation procedure is reported as a part of an ab initio program (named SYCETY) developed recently for spatially extended molecules. It is shown that all advantages of the program SYCETY can be kept. The procedure is especially suitable to transform the integrals to a localized basis set. By using the MBPT with localized basis the most important part of the correlation can be calculated even for large (spatially extended) systems.     

Electronic correlation studies. III. Self-correlated field method. Application to 2S ground state and 2P excited state of three-electron atomic systems

The self-correlated field method is based on the insertion in the group product wave function of pair functions built upon a set of correlated “local” functions and of “nonlocal” functions. This work is an application to three-electron systems. The effects of the outer electron on the inner pair are studied. The total electronic energy and some intermediary results such as pair energies, Coulomb and exchange “correlated” integrals, are given. The results are always better than those given by conventional SCF computations and reach the same level of accuracy as those given by more laborious methods used in correlation studies.     

Orbital-invariant formulation and second-order gradient evaluation in Møller-Plesset perturbation theory

Based on the Hylleraas functional form, the second and third orders of Møller-Plesset perturbation theory are reformulated in terms of arbitrary (e.g., localized) internal orbitals, and atomic orbitals in the virtual space. The results are strictly equivalent to the canonical formulation if no further approximations are introduced. The new formalism permits the extension of the local correlation method to Møller-Plesset theory. It also facilitates the treatment of weak pairs at a lower (e.g., second order) level of theory in CI and coupled cluster methods. Based on our formalism, an MP2 gradient algorithm is outlined which does not require the storage of derivative integrals, integrals with three external MO indices, and, using the method of Handy and Schaefer, the repeated solution of the coupled-perturbed SCF equations.