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.     

Matrix elements of N-particle explicitly correlated Gaussian basis functions with complex exponential parameters

In this work we present analytical expressions for Hamiltonian matrix elements with spherically symmetric, explicitly correlated Gaussian basis functions with complex exponential parameters for an arbitrary number of particles. The expressions are derived using the formalism of matrix differential calculus. In addition, we present expressions for the energy gradient that includes derivatives of the Hamiltonian integrals with respect to the exponential parameters. The gradient is used in the variational optimization of the parameters. All the expressions are presented in the matrix form suitable for both numerical implementation and theoretical analysis. The energy and gradient formulas have been programmed and used to calculate ground and excited states of the He atom using an approach that does not involve the Born-Oppenheimer approximation.     

The evaluation of matrix elements in the analysis of anharmonic molecular vibrations: Functional taylor series expansions

This article presents methods for computing Cartesian Gaussian matrix elements using a Taylor series for general potential energy operators that admit well-behaved radial derivatives. These operators arise in the analyses of anharmonic vibrations in molecules. Application to the evaluation of matrix elements for hydrogen associated two wells illustrates the method. © 1995 John Wiley & Sons, Inc.     

Second quantization for systems with a constant number of particles in the Dirac notation

The relation between the completeness condition for an appropriate one-particle basis set and the occupation number representation (second quantization) is shown for the time-independent case. The explicit expressions for the basic symmetric operators are derived in the Dirac bra–ket notation. The physical meaning of these operators, the algebra as well as the connections with the one-electron density matrix and with the projector on the Fermi sea in the one-electron approximation, follow directly from these expressions. The generalization for a nonorthogonal basis and the algebra for corresponding basic operators are formulated. The connection with the notion of the molecular diagrams of different kinds for the nonorthogonal atomic orbitals is shown. The Mulliken populations and the Chirgwin–Coulson bond orders are equal to the diagonal and offdiagonal elements of the molecular diagram 1, respectively. The matrix elements of the projector on the Fermi sea in the one-electron approximation in the representation of nonorthogonal atomic orbitals are elements of the molecular diagram 2.     

Analytical treatment of Coriolis coupling for three-body systems

In a previous article [J. Chem. Phys. 108 (1998) 5216], an efficient method was presented for performing “exact” quantum calculations for the three-body rovibrational Hamiltonian, within the helicity-conserving approximation. This approach makes use of a certain three-body “effective potential,” enabling the same bend angle basis set to be employed for all values of the rotational quantum numbers, J, K and M. In the present work, the method is extended to incorporate Coriolis coupling, for which the relevant matrix elements are derived exactly. These can be used to solve the full three-body rovibrational problem, in the standard Jacobi coordinate vector embedding. Generalization of the method for coupled kinetic energy operators arising from other coordinate systems, embeddings, and/or system sizes, is also discussed.     

Hierarchical expansion of the kinetic energy operator in curvilinear coordinates for the vibrational self-consistent field method

A new hierarchical expansion of the kinetic energy operator in curvilinear coordinates is presented and modified vibrational self-consistent field (VSCF) equations are derived including all kinematic effects within the mean field approximation. The new concept for the kinetic energy operator is based on many-body expansions for all G matrix elements and its determinant. As a test application VSCF computations were performed on the H(2)O(2) molecule using an analytic potential (PCPSDE) and different hierarchical approximations for the kinetic energy operator. The results indicate that coordinate-dependent reduced masses account for the largest part of the kinetic energy. Neither kinematic couplings nor derivatives of the G matrix nor its determinant had significant effects on the VSCF energies. Only the zero-point value of the pseudopotential yields an offset to absolute energies which, however, is irrelevant for spectroscopic problems.     

Tensor formalism in anharmonic calculations

A new method is presented to compute cartesian tensors in the expansion of curvilinear internal coordinates. Second- and higher-order coefficients are related to the metrics of the space of displacements. Components of the metric tensor are taken from existing tables of inverse kinetic energy matrix elements or, when rotations are involved, derived from general invariance conditions of scalars within a molecule. This leads to a tensor formalism particularly convenient in dealing with curvilinear coordinates in anharmonic calculations of vibrational frequencies. Formulae are given for elements of the potential energy matrix, related to quadratic and cubic force constants in terms of Christoffel symbols. The latter quantities are also used in the expansion of redundancy relations, with explicit coefficients given up to the third order.     

Second-order correlation potential in the Kohn–Sham approximation for atoms

A new type of correlation functional derived from the second-order expression for the correlation energy of an atom is proposed. The derived correlation potential contains one free parameter, which is determined by fitting the known pair correlation energy. The calculations with this potential in the Kohn–Sham approximation give rather accurate values for the matrix elements of different operators.     

Real‐space grid representation of momentum and kinetic energy operators for electronic structure calculations

We show that the central finite difference formula for the first and the second derivative of a function can be derived, in the context of quantum mechanics, as matrix elements of the momentum and kinetic energy operators on discrete coordinate eigenkets defined on a uniform grid. Starting from the discretization of integrals involving canonical commutations, simple closed‐form expressions of the matrix elements are obtained. A detailed analysis of the convergence toward the continuum limit with respect to both the grid spacing and the derivative approximation order is presented. It is shown that the convergence from below of the eigenvalues in electronic structure calculations is an intrinsic feature of the finite difference method. © 2018 Wiley Periodicals, Inc.     

Functional derivative of the kinetic energy functional for spherically symmetric systems

Ensemble non-interacting kinetic energy functional is constructed for spherically symmetric systems. The differential virial theorem is derived for the ensemble. A first-order differential equation for the functional derivative of the ensemble non-interacting kinetic energy functional and the ensemble Pauli potential is presented. This equation can be solved and a special case of the solution provides the original non-interacting kinetic energy of the density functional theory.     

Kratzer potential algebraic representation and matrix elements recurrence formulae

Generalized recurrence relations for the calculation of multipole matrix elements for Kratzer potential wave functions are obtained operationally. These formulas have been determined by using a non-analytical procedure based on the algebraic representation of the Kratzer eigenfunctions along with the usual ladder properties and commutation relations. For that, the creation and annihilation operators are adequately derived by means of an alternative approach to the factorization method and the exact expressions for matrix elements are achieved with the aid of a relationship between the ladder operators associated with the bra and theket. The proposed algebraic approach as well as the formulas for the calculation of matrix elements thus derived are quite simple and direct when compared with other alternative expressions already obtained analytically or pseudo-algebraically by means of the hypervirial theorem commutator algebra.     

Matrix elements for Ĵ2 and Ĵz operators over explicitly correlated Cartesian Gaussian functions

General formalism for evaluation of multiparticle integrals involving J?² and J?_z operators over explicitly correlated Cartesian Gaussian functions is presented. The integrals are expressed in terms of the general overlap integrals. An explicitly correlated Cartesian Gaussian function is a product of spherical orbital Gaussian functions, powers of the Cartesian coordinates of the particle, and exponential Gaussian factors, which depend on interparticular distances. This development is relevant to both adiabatic and nonadiabatic calculations of energy and properties of multiparticle systems. © 1995 John Wiley & Sons, Inc.     

Introduction

The formalism for a configuration interaction approach is presented, in which explicit introduction of interelectronic coordinates into individual configurations is accomplished through the use of spherical Gaussian correlation factors. Formulas for all required matrix elements over these configurations are derived, and a computationally convenient algorithm to evaluate the matrix elements is presented. In addition, analysis of the form of the correlation factor shows that both angular and radial correlation can be obtained using spherical Gaussian correlation factors.     

CALCULATION OF THE TRANSITION MATRIX ELEMENTS FOR AN ATOM-DIATOMIC MOLECULE SCATTERING

In this article, the IOS and the GIOS approximate transition matrix elements and theirapplicable conditions, the potential energy condition, mass conditions and kinetic energycondition, have been derived by using the energy shift operator.     

Implementation of gradient formulas for correlated gaussians: He, ∞He,Ps2, 9Be,and ∞Be test results

New formulas in the basis of explicitly correlated Gaussian basis functions, derived in a previous article using powerful matrix calculus, are implemented and applied to find variational upper bounds for nonrelativistic ground states of ⁴He, ^∞He, Ps₂, ⁹Be, and ^∞Be. Analytic gradients of the energy are included to speed optimization of the exponential variational parameters. Five different nonlinear optimization subroutines (algorithms) are compared: TN, truncated Newton; DUMING, quasi-Newton; DUMIDH, modified Newton; DUMCGG, conjugate gradient; and POWELL, direction set (nongradient). The new analytic gradient formulas are found to significantly accelerate optimizations that require gradients. We found that the truncated Newton algorithm out-performs the other optimizers for the selected test cases. Computer timings and energy bounds are reported. © 1997 John Wiley & Sons, Inc.     

Analytical evaluation of pseudopotential matrix elements with Gaussian‐type solid harmonics of arbitrary angular momentum

An efficient formalism for evaluating pseudopotential matrix elements with Gaussian‐type solid harmonics of arbitrary angular momentum is presented. It is based on the tensor coupling technique, which is especially well suited for treating Gaussian‐type solid harmonics of arbitrary angular momentum. Closed analytical expressions are derived for the matrix elements as well as for their nuclear displacement derivatives. The efficiency of the implementation into our new parallel density functional program PARA GAUSS and the quality of the pseudopotential approach is tested for a set of representative molecules and cluster models. To this end the results of pseudopotential calculations are compared to those of nonrelativistic and scalar‐relativistic all‐electron calculations. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 209–221, 2000     

Harmonic oscillator tensors. II. angular momentum expressions of matrix elements of vibrational operators

The spatial symmetries of the harmonic oscillator and the recently found irreducible tensors constructed from the associated annihilation and creation operators are exploited to obtain new expressions for the elements of the matrix representatives of several examples of vibrational operators. Since all vibrational operators are expressible in terms of the irreducible tensors, their matrix elements reflect the angular momentum symmetry inherent in them, for the results derived here are in terms of the Clebsch–Gordan coefficients and the isoscalar factors that arise from the couplinig rule of the irreducible tensors. Familiarity with the mathematical properties of these quantities derived from the elementary theory of angular momentum facilitates the evaluation of many vibrational operators that may be of importance in the study of potentials in this basis. In particular, it is shown that the nonvanishing of matrix elements is governed by a law of conservation of angular momentum along the axis of quantization of the nondegenerate harmonic oscillator. © 1993 John Wiley & Sons, Inc.     

The algebraic approach for the derivation of ladder operators and coherent states for the Goldman and Krivchenkov oscillator by the use of supersymmetric quantum mechanics

The ladder operators for the Goldman and Krivchenkov anharmonic potential have been derived within the algebraic approach. The method is extended to include the rotating oscillator. The coherent states for the Goldman and Krivchenkov oscillator, which are the eigenstates of the annihilation operator and minimize the generalized position-momentum uncertainty relation, are constructed within the framework of supersymmetric quantum mechanics. The constructed ladder operators can be a useful tool in quantum chemistry computations of non-trivial matrix elements. In particular, they can be employed in molecular vibrational–rotational spectroscopy of diatomic molecules to compute transition energies and dipole matrix elements.     

Unitary group approach to the many-electron problem. II. Adjoint tensor operators for U(n)

In this paper we present a derivation of the U(n) adjoint coupling coefficients for the representations appropriate to many-electron systems. Since the states of a many-fermion system are to comprise the totally antisymmetric Nth rank tensor representation of U(2n), the work of this paper enables the matrix elements of the U(2n) generators to be evaluated directly in the U(n) × U(2) (i.e., spin orbit) basis using their transformation properties as adjoint tensor operators. A connection between the adjoint coupling coefficients, as derived in this paper, and the matrix elements of certain (spin independent) two-body operators is also presented. This indicates that in CI calculations, one may obtain the matrix elements of spin-dependent operators from the known matrix elements of certain spin-independent two-body operators. In particular this implies a segment-level formula for the matrix elements of the U(2n) generators in the spin-orbit basis.     

Correlated wave-function theory for many-center many-electron problems

The correlated electronic wave-function theory developed by S. Obara and K. Hirao [Bull. Chem. Soc. Jpn. 66 , 3300 (1993)], as applied to two-electron molecular systems, is generalized to many-center many-electron problems. The exact formulas for effective Hamiltonian operators are given. The rules for the calculation of matrix elements with three-electron operators over Slater determinants are formulated. From the energy-minimum principle, the system of master equations is derived for variational coefficients of a trial wave function for the molecules with closed electronic shells. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 639–648, 1998