Quantum-Chemical Calculation of Matrix Elements in a Basis of Functions with Polynomial Tails

A numerical integration method is suggested for calculating Hamiltonian matrix elements in a basis of functions with polynomial tails with allowance for discontinuities of higher-order derivatives of the basis function within the domain of integration. The method is tested by calculating matrix elements for a copper crystal. The results for the overlap matrix elements are presented demonstrating efficiency of the method.     

New analytical expressions,symmetry relations and numerical solutions for the rotational overlap integrals

In this article, extremely simple analytical formulas are obtained for rotational overlap integrals which occur in integrals over two reduced rotation matrix elements. The analytical derivations are based on the properties of the Jacobi polynomials and beta functions. Numerical results and special values for rotational overlap integrals are obtained by using symmetry properties and recurrence relationships for reduced rotation matrix elements. The final results are of surprisingly simple structures and very useful for practical applications. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Energy and energy gradient matrix elements with N-particle explicitly correlated complex Gaussian basis functions with L=1

In this work we consider explicitly correlated complex Gaussian basis functions for expanding the wave function of an N-particle system with the L=1 total orbital angular momentum. We derive analytical expressions for various matrix elements with these basis functions including the overlap, kinetic energy, and potential energy (Coulomb interaction) matrix elements, as well as matrix elements of other quantities. The derivatives of the overlap, kinetic, and potential energy integrals with respect to the Gaussian exponential parameters are also derived and used to calculate the energy gradient. All the derivations are performed using the formalism of the matrix differential calculus that facilitates a way of expressing the integrals in an elegant matrix form, which is convenient for the theoretical analysis and the computer implementation. The new method is tested in calculations of two systems: the lowest P state of the beryllium atom and the bound P state of the positronium molecule (with the negative parity). Both calculations yielded new, lowest-to-date, variational upper bounds, while the number of basis functions used was significantly smaller than in previous studies. It was possible to accomplish this due to the use of the analytic energy gradient in the minimization of the variational energy.     

Calculation of molecular polarizabilities for large systems within the random phase approximation

We examine the static-field molecular polarizability from a sum over uncoupled Hartree—Fock states (SOS), the Tamm—Dancoff approximation (TDA), and the random phase approximation (RPA). An efficient algorithm for the inversion of the TDA or RPA matrix is outlined, which avoids matrix diagonalization and explicit construction of matrix elements over states, allowing for rapid calculation of molecular polarizabilities. The extension of the method is straightforward; third-order hyperpolarizability is developed as an example. Test cases are reported for molecules represented by an intermediate neglect of differential overlap (INDO) wavefunction.     

Adjustment of Born‐Oppenheimer electronic wave functions to simplify close coupling calculations

Technical problems connected with use of the Born‐Oppenheimer clamped‐nuclei approximation to generate electronic wave functions, potential energy surfaces (PES), and associated properties are discussed. A computational procedure for adjusting the phases of the wave functions, as well as their order when potential crossings occur, is presented which is based on the calculation of overlaps between sets of molecular orbitals and configuration interaction eigenfunctions obtained at neighboring nuclear conformations. This approach has significant advantages for theoretical treatments describing atomic collisions and photo‐dissociation processes by means of ab initio PES, electronic transition moments, and nonadiabatic radial and rotational coupling matrix elements. It ensures that the electronic wave functions are continuous over the entire range of nuclear conformations considered, thereby greatly simplifying the process of obtaining the above quantities from the results of single‐point Born‐Oppenheimer calculations. The overlap results are also used to define a diabatic transformation of the wave functions obtained for conical intersections that greatly simplifies the computation of off‐diagonal matrix elements by eliminating the need for complex phase factors. © 2013 Wiley Periodicals, Inc.     

Transformation of Bloch representation into atomic orbital representation and combined matrix element calculation technique for spatially bounded basis functions in crystals

A method is proposed for transforming the Hamiltonian from Bloch to atomic function representation. For spatially bounded functions, this is a rigorous method based on solution of a certain algebraic system of equations. Unlike the conventional procedure based on integration over the Brillouin zone, the new method requires knowledge of the matrix elements of the Bloch representation only at several points of the Brillouin zone. The number of these points is determined by the trimming radius for the spatially bounded functions and by the lattice constant. The method can be used for calculating matrix elements in a basis of atomic functions and for reducing computations in matrix element calculations of the Bloch representation for procedures using numerical integration.     

Integrals and derivatives for correlated Gaussian functions using matrix differential calculus

The matrix differential calculus is applied for the first time to a quantum chemical problem via new matrix derivations of integral formulas and gradients for Hamiltonian matrix elements in a basis of correlated Gaussian functions. Requisite mathematical background material on Kronecker products, Hadamard products, the vec and vech operators, linear structures, and matrix differential calculus is presented. New matrix forms for the kinetic and potential energy operators are presented. Integrals for overlap, kinetic energy, and potential energy matrix elements are derived in matrix form using matrix calculus. The gradient of the energy functional with respect to the correlated Gaussian exponent matrices is derived. Burdensome summation notation is entirely replaced with a compact matrix notation that is both theoretically and computationally insightful. © 1996 John Wiley & Sons, Inc.     

Limitations of the diatomics-in-molecules (DIM) method with neglect of overlap

The approximations made in the diatomics-in-molecules (DIM) method are examined by comparing the matrix elements of the DIM calculation with those of the ab initio calculation. It was found that if the overlap between atomic orbitals is small, the ab initio matrix elements can be expressed as the sums of DIM matrix elements and three-center terms ignored in the DIM method. The zero overlap of the atomic orbital approximation used in the DIM method makes the three-center terms simpler and apparently smaller which justifies the applicability of this approximation to the DIM method. Results of calculations show that deviations of the DIM values from accurate values are particularly large in the linear HXH configuration where X is a many-electron atom. It is suggested that a supplementary term be added to the DIM energy (DIM-SC) to take into account the neglected three-center terms.     

Matrix elements of multi-electron operator in restricted Hartree-Fock theory

A theorem has been established for the matrix elements of a general t-electron operator between N-dimensional determinantal wave functions arising in the solution of the atomic and molecular multi-electron problem by the restricted Hartree-Fock (RHF) method (1 ≤ t ≤ N). The required matrix elements of this operator are sums of matrix elements over t-dimensional basic determinantal wave functions. The final results are especially useful in the determination of multi-electron properties for atoms and molecules when the Hylleraas approach in RHF theory is employed.     

Improved method for calculating vibronic spectra of complex molecules

The previously suggested approximate method for calculating the overlap integrals of vibrational wave functions is considerably improved for the purpose of maximally accurate calculation of excitation-induced mixing of normal coordinates. A general formula is obtained for all types of overlap integrals as a finite power series of the potential surface shift parameters; the coefficients are derivatives of the corresponding generating Junctions represented as polynomials of the shift vector elements of the normal coordinates and the mixing matrix. The spectra of model molecules of decatetraene and tetra- and hexadecaheptaene were calculated using the expressions derived in this work and a semiempirical parametric method for determination of excitation-induced changes in the potential surface of molecules. The calculations confirmed the high efficiency of both the parametric method and the new technique. Translated fromZhumal Strukturnoi Khimii, Vol. 38, No. 2, pp. 240–247, March–April, 1997.     

Spin-free approach for evaluation of electronic matrix elements using character operators of ℒN

A spin-free method is presented for evaluating electronic matrix elements over a spin-independent many-electron Hamiltonian. The spin-adapted basis of configuration state functions is obtained using a nonorthogonal spin basis consisting of projected spin eigenfunctions. The general expressions for the matrix elements are given explicitly, and it is demonstrated how the matrix elements may be obtained simply from the knowledge of the irreducible characters of the permutation group ℒ_N. The presented formulas are very general and may be applied in connection with both spin-coupled valence bond studies and in conventional configuration interaction (CI) methods based on an orthonormal orbital basis. © 1996 John Wiley & Sons, Inc.     

Nonadiabatic theory of electron–vibrational coupling: New basis for microscopic interpretation of superconductivity

A systematic procedure has been developed to construct an electronic energy matrix for diatomics in the basis of antisymmetrized products of atomic wave functions represented as linear combinations of coupled momenta eigenfunctions. The exchange matrix element is expanded in powers of electronic interchange between atoms. General expressions of many-electron angular coefficients have been obtained for all types of products of one- and two-electron and overlap integrals in energy matrix elements. © 1992 John Wiley & Sons, Inc.     

Fock transforms in reciprocal-space quantum theory

In Fock's reciprocal-space treatment of the hydrogen atom,k-space is mapped onto the surface of a 4-dimensional hypersphere, and the solutions (apart from an invariant factor) are 4-dimensional hyperspherical harmonics. Fock's method can be generalized to provide solutions for the Schrödinger equation of a charged particle moving in a many-center Coulomb potential, and in this case the solutions are found by diagonalizing an overlap matrix involving products of hyperspherical harmonics. The present paper discusses a transform which can conveniently be used to evaluate the elements of the overlap matrix.     

Errors in molecular integrals: An influence on RHF energy values

The influence of errors in molecular integrals on the calculated RHF energy values is considered for two models which correspond to round-off and shift errors. The energy variations induced by errors in elements of one- and two-electron matrices and the overlap matrix are represented in a quadratic approximation and the same degree of accuracy is maintained for mean values and standard deviations. The formulas given point out that mean values are less than zero when some of the errors are non-shifted. Special care is required when Gaussian lobe functions are used.     

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.     

Paired-permanent approach to valence bond theory

A new function called paired-permanent is defined and widely discussed, and a practicable procedure for evaluations of paired-permanents is proposed, which is similar to the Laplace method for determinants. Using the concept of paired-permanents, an efficient algorithm is presented for evaluating the Hamiltonian and overlap matrix elements in the spin-free form of valence bond (VB) theory. With the new algorithm, a spin-free wavefunction is simply written as a paired-permanent, and an overlap matrix element may be obtained by evaluating a corresponding paired-permanent. Meanwhile, the Hamiltonian matrix element is expressed in terms of the summation of the products of electronic integrals and the corresponding sub-paired-permanents     

Transforming Gaussians into Wannier functions using inversion of cyclic matrix

One-dimensional Gaussians are transformed into orthonormalized atomic orbitals (OAOS ) using the inverse of the cyclic overlap matrix. During that process it was found that these OAOS are nothing but the Wannier functions and have properties similar to the Wannier functions obtained by Wannier himself from Gaussians using a different method. The present method gives rise to a better Wannier function with less distortion.     

A diagrammatic valence bond method for configuration interaction calculations in atoms and molecules

A diagrammatic valence bond method based on Rumer-Pauling rules for configuration interaction calculations is described. The advantages of this method are that it is simple and flexible and is expected to be computationally efficient as the basis functions can be coded as increasing integers. Evaluation of Hamiltonian matrix elements involves simple bit manipulations and binary searches. The basis, being represented pictorially, should also help in utilizing spatial symmetries for further block-diagonalizing the Hamiltonian matrix. The eigenfunctions of the Hamiltonian can also be used to compute matrix elements between different electronic states.     

Efficient algorithm for the spin-free valence bond theory. I. New strategy and primary expressions

A new algorithm for nonorthogonal valence bond (VB) method is presented by using symmetric group approach. In the present algorithm, a new function, called paired-permanent-determinant (PPD), is defined, which is an algebrant and has the same symmetry of a corresponding VB structure. The evaluation of a PPD is carried out by using a recursion formula similar to the Laplace expansion method for determinants. An overlap matrix element in the spin-free VB method may be obtained by evaluating a corresponding PPD, while the Hamiltonian matrix element is expressed in terms of the products of electronic integrals and sub-PPDs. In the present work, some important properties of PPDs are discussed, and the primary procedure for the evaluation of PPD is deduced. Furthermore, the expressions for evaluating both the overlap and Hamiltonian matrix elements are also given in details, which are essential to develop an efficient algorithm for nonorthogonal VB calculations. In the present study, some further effective technical considerations will be adopted, and a new ab initio VB program will be introduced. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 287–297, 1998     

Berry phase approach to longitudinal dipole moments of infinite chains in electronic-structure methods with local basis sets

The authors provide a reformulation of the modern theory of polarization for one-dimensional stereoregular polymers, at the level of the single determinant Hartree-Fock and Kohn-Sham methods within a basis set of local orbitals. By starting with localization of one-electron orbitals, their approach naturally arrives to the Berry phases of Bloch orbitals. Then they describe a novel numerical algorithm for evaluation of longitudinal dipole moments, computationally more convenient than those presently implemented within the local basis periodic codes. This method is based on the straightforward evaluation of the usual direct space dipole matrix elements between local orbitals, as well as overlap matrices between wave functions at two neighboring k points of the reciprocal space mesh. The practical behavior of the algorithm and its convergence properties with respect to the k-point mesh density are illustrated in benchmark calculations for water chains and fluorinated trans-polyacetylene.