Gaussian and finite-element Coulomb method for the fast evaluation of Coulomb integrals

The authors propose a new linear-scaling method for the fast evaluation of Coulomb integrals with Gaussian basis functions called the Gaussian and finite-element Coulomb (GFC) method. In this method, the Coulomb potential is expanded in a basis of mixed Gaussian and finite-element auxiliary functions that express the core and smooth Coulomb potentials, respectively. Coulomb integrals can be evaluated by three-center one-electron overlap integrals among two Gaussian basis functions and one mixed auxiliary function. Thus, the computational cost and scaling for large molecules are drastically reduced. Several applications to molecular systems show that the GFC method is more efficient than the analytical integration approach that requires four-center two-electron repulsion integrals. The GFC method realizes a near linear scaling for both one-dimensional alanine alpha-helix chains and three-dimensional diamond pieces.     

New recurrence relations for the rapid evaluation of electron repulsion integrals based on the accompanying coordinate expansion formula

We present an algorithm for the rapid computation of electron repulsion integrals (ERIs) over Gaussian basis functions based on the accompanying coordinate expansion (ACE) formula. The present algorithm uses equations termed angular momentum reduced expressions and introduces two types of recurrence relations to ACE formulas. Numerical efficiencies are assessed for (p pmid R:p p) and (sp spmid R:sp sp) ERIs by using the floating-point operation count. The algorithm is suitable for calculating ERIs for the same exponents but different angular momentum functions, such as L shells and derivatives of ERIs. The present algorithm is also capable of calculating ERIs with highly contracted Gaussian basis functions.     

Methods for efficient evaluation of integrals for Gaussian type basis sets

Methods are described that allow for an efficient evaluation of two-electron integrals over contracted Gaussian lobe functions. The improvement in computational speed is achieved by avoiding the computation of integrals that are: 1. sufficiently small on numerical reasons, 2. zero by symmetry, 3. identical to other integrals by symmetry. Examples of the effectiveness of these techniques are included. We also report the timings for a further processing of two-electron integrals in a Hartree Fock and correlation energy computation.     

The prism algorithm for two-electron integrals

A new approach to the evaluation of two-electron repulsion integrals over contracted Gaussian basis functions is developed. The new scheme encompasses 20 distinct, but interrelated, paths from simple shell-quartet parameters to the target integrals, and, for any given integral class, the path requiring the fewest floating-point operations (FLOPS ) is that used. Both theoretical (FLOP counting) and practical (CPU timing) measures indicate that the method represents a substantial improvement over the HGP algorithm.     

Auxiliary functions for molecular integrals with Slater‐type orbitals. I. Translation methods

Many types of molecular integrals involving Slater functions can be expressed, with the ζ‐function method in terms of sets of one‐dimensional auxiliary integrals whose integrands contain two‐range functions. After reviewing the properties of these functions (including recurrence relations, derivatives, integral representations, and series expansions), we carry out a detailed study of the auxiliary integrals aimed to facilitate both the formal and computational applications of the ζ‐function method. The usefulness of this study in formal applications is illustrated with an example. The high performance in numerical applications is proved by the development of a very efficient program for the calculation of two‐center integrals with Slater functions corresponding to electrostatic potential, electric field, and electric field gradient. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A general formulation for the efficient evaluation of n-electron integrals over products of Gaussian charge distributions with Gaussian geminal functions

In this work, we present a general formulation for the evaluation of many-electron integrals which arise when traditional one particle expansions are augmented with explicitly correlated Gaussian geminal functions. The integrand is expressed as a product of charge distributions, one for each electron, multiplied by one or more Gaussian geminal factors. Our formulation begins by focusing on the quadratic form that arises in the general n-electron integral. Using the Rys polynomial method for the evaluation of potential energy integrals, we derive a general formula for the evaluation of any n-electron integral. This general expression contains four parameters ω, θ, v, and h, which can be evaluated by an examination of the general quadratic form. Our analysis contains general expressions for any n-electron integral over s-type functions as well as the recursion needed to build up arbitrary angular momentum. The general recursion relation requires at most n + 1 terms for any n-electron integral. To illustrate the general method, we develop explicit expressions for the evaluation of two, three, and four particle electron repulsion integrals as well as two and three particle overlap and nuclear attraction integrals. We conclude our exposition with a discussion of a preliminary computational implementation as well as general computational requirements. Implementation on parallel computers is briefly discussed.     

Recursion formula for electron repulsion integrals over Hermite polynomials

A method for computing electron repulsion integrals over contracted Gaussian functions is described in which intermediate integrals over Hermite polynomials are generated by a “pre‐Hermite” recursion (PHR) step before the conversion to regular integrals. This greatly reduces the floating‐point operation counts inside the contraction loops, where only simple “scaling”‐type operations are required, making the method efficient for contracted Gaussians, particularly of high angular momentum. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Improved partition–expansion of two‐center distributions involving slater functions

The calculation of the electronic structure of large systems is facilitated by the substitution of the two‐center distributions by their projections on auxiliary basis sets of one‐center functions. An alternative is the partition–expansion method in which one first decides what part of the distribution is assigned to each center, and next expands each part in spherical harmonics times radial factors. The method is exact, requires neither auxiliary basis sets nor projections, and can be applied to Gaussian and Slater basis sets. Two improvements in the partition–expansion method for Slater functions are reported: general expressions valid for arbitrary quantum numbers are derived and the efficiency of the procedure is increased giving analytical solutions to integrals previously computed by numerical quadrature. The efficiency of the new version is assessed in several molecules and the advantages over the projection methods are pointed out. © 2013 Wiley Periodicals, Inc.     

Analytical GIAO and hybrid-basis integral derivatives: application to geometry optimization of molecules in strong magnetic fields

Analytical integral evaluation is a central task of modern quantum chemistry. Here we present a general method for evaluating differentiated integrals over standard Gaussian and mixed Gaussian/plane-wave hybrid orbitals. The main idea is to have a representation of basis sets that is flexible enough to enable differentiated integrals to be reinterpreted as standard integrals over modified basis functions. As an illustration of the method, we report a very simple implementation of Hartree-Fock level geometrical derivatives in finite magnetic fields for gauge-origin independent atomic orbitals, within the London program. As a quantum-chemical application, we optimize the structure of helium clusters and some well-known covalently bound molecules (water, ammonia and benzene) subject to strong magnetic fields.     

A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians

Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals. For two- and three-center two-electron integrals, the solid-harmonic Hermite scheme is particularly efficient, significantly reducing the cost relative to the Cartesian scheme.     

Libcint: An efficient general integral library for Gaussian basis functions

An efficient integral library Libcint was designed to automatically implement general integrals for Gaussian‐type scalar and spinor basis functions. The library is able to evaluate arbitrary integral expressions on top of p, r and σ operators with one‐electron overlap and nuclear attraction, two‐electron Coulomb and Gaunt operators for segmented contracted and/or generated contracted basis in Cartesian, spherical or spinor form. Using a symbolic algebra tool, new integrals are derived and translated to C code programmatically. The generated integrals can be used in various types of molecular properties. To demonstrate the capability of the integral library, we computed the analytical gradients and NMR shielding constants at both nonrelativistic and 4‐component relativistic Hartree–Fock level in this work. Due to the use of kinetically balanced basis and gauge including atomic orbitals, the relativistic analytical gradients and shielding constants requires the integral library to handle the fifth‐order electron repulsion integral derivatives. The generality of the integral library is achieved without losing efficiency. On the modern multi‐CPU platform, Libcint can easily reach the overall throughput being many times of the I/O bandwidth. On a 20‐core node, we are able to achieve an average output 8.3 GB/s for C₆₀ molecule with cc‐pVTZ basis. © 2015 Wiley Periodicals, Inc.     

Molecular integrals for Gaussian and exponential‐type functions: Shift operators

Basis functions with arbitrary quantum numbers can be attained from those with the lowest numbers by applying shift operators. We derive the general expressions and the recurrence relations of these operators for Cartesian basis sets with Gaussian and exponential radial factors. In correspondence, the expressions of molecular integrals involving functions with arbitrary quantum numbers can be obtained by applying these operators on the integrals with the lowest quantum numbers. Since the original form of the shift operators is not appropriate to deal with integrals, we give their representation in terms of derivatives with respect to the parameters on which these integrals explicitly depend. Moreover, we translate the recurrence relations to the new representation and, finally, we analyze the general expressions ot the molecular integrals. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 137–145, 2000     

Auxiliary functions for molecular integrals with Slater‐type orbitals. II. Gauss transform methods

The Gauss transform of Slater‐type orbitals is used to express several types of molecular integrals involving these functions in terms of simple auxiliary functions. After reviewing this transform and the way it can be combined with the shift operator technique, a master formula for overlap integrals is derived and used to obtain multipolar moments associated to fragments of two‐center distributions and overlaps of derivatives of Slater functions. Moreover, it is proved that integrals involving two‐center distributions and irregular harmonics placed at arbitrary points (which determine the electrostatic potential, field and field gradient, as well as higher order derivatives of the potential) can be expressed in terms of auxiliary functions of the same type as those appearing in the overlap. The recurrence relations and series expansions of these functions are thoroughly studied, and algorithms for their calculation are presented. The usefulness and efficiency of this procedure are tested by developing two independent codes: one for the derivatives of the overlap integrals with respect to the centers of the functions, and another for derivatives of the potential (electrostatic field, field gradient, and so forth) at arbitrary points. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Multicenter and multiparticle integrals for explicitly correlated cartesian gaussian-type functions

An analytical derivation of multicenter and multiparticle integrals for explicitly correlated Cartesian Gaussian-type cluster functions is demonstrated. The evaluation method is based on the application of raising operators that transform spherical cluster Gaussian functions into Cartesian Gaussian functions.     

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.     

Strategies for an efficient implementation of the Gauss-Bessel quadrature for the evaluation of multicenter integral over STFs

In a previous work, a new Gauss quadrature was introduced with a view to evaluate multicenter integrals over Slater-type functions efficiently. The complexity analysis of the new approach, carried out using the three-center nuclear integral as a case study, has shown that for low-order polynomials its efficiency is comparable to the SD. The latter was developed in connection with multi-center integrals evaluated by means of the Fourier transform of B functions. In this work we investigate the numerical properties of the Gauss-Bessel quadrature and devise strategies for an efficient implementation of the numerical algorithms for the evaluation of multi-center integrals in the framework of the Gaussian transform/Gauss-Bessel approach. The success of these strategies are essential to elaborate a fast and reliable algorithm for the evaluation of multi-center integrals over STFs.     

Basis set construction for molecular electronic structure theory: natural orbital and Gauss-Slater basis for smooth pseudopotentials

A simple yet general method for constructing basis sets for molecular electronic structure calculations is presented. These basis sets consist of atomic natural orbitals from a multiconfigurational self-consistent field calculation supplemented with primitive functions, chosen such that the asymptotics are appropriate for the potential of the system. Primitives are optimized for the homonuclear diatomic molecule to produce a balanced basis set. Two general features that facilitate this basis construction are demonstrated. First, weak coupling exists between the optimal exponents of primitives with different angular momenta. Second, the optimal primitive exponents for a chosen system depend weakly on the particular level of theory employed for optimization. The explicit case considered here is a basis set appropriate for the Burkatzki-Filippi-Dolg pseudopotentials. Since these pseudopotentials are finite at nuclei and have a Coulomb tail, the recently proposed Gauss-Slater functions are the appropriate primitives. Double- and triple-zeta bases are developed for elements hydrogen through argon. These new bases offer significant gains over the corresponding Burkatzki-Filippi-Dolg bases at various levels of theory. Using a Gaussian expansion of the basis functions, these bases can be employed in any electronic structure method. Quantum Monte Carlo provides an added benefit: expansions are unnecessary since the integrals are evaluated numerically.     

Adaptive density partitioning technique in the auxiliary plane wave method

We have developed the adaptive density partitioning technique (ADPT) in the auxiliary plane wave method, in which a part of the density is expanded to plane waves, for the fast evaluation of Coulomb matrix. Our partitioning is based on the error estimations and allows us to control the accuracy and efficiency. Moreover, we can drastically reduce the core Gaussian products that are left in Gaussian representation (its analytical integrals is the bottleneck in this method). For the taxol molecule with 6-31G** basis, the core Gaussian products accounted only for 5% in submicrohartree error.     

Gaussian matrix elements of the operator X: Reduction to confluent hypergeometric functions

A new method has been proposed for calculating the Gaussian matrix elements of the interaction operator represented by an arbitrary degree of interelectron distance X. The method is based on the expansion of two-electron integrals as the sum of one-electron integrals which in turn admit compact operator representation in terms of confluent hypergeometric functions. The generating differential operator has been shown to be related to the modified Hermitian polynomials. The standard structure of the special functions encountered in this approach is useful in studying the analytical behavior of the integrals and makes it possible to obtain for these integrals recurrence relations, direct algebraic expressions in the forms of finite sum of confluent hypergeometric functions, integral representations, and asymptotic properties. Unlike the usual methods based on integral transformation of the interaction operator, the proposed approach has a wider field of application, and in addition, leads to compact and convenient analytical expressions. The idea of using differential properties of integrals to simplify the integrand structure gives the proposed approach a certain resemblance to that suggested by Boys but not developed in detail in his pioneer work.     

Many-Electron,multicenter integrals in superposition of correlated configurations method. I. one- and two-electron integrals

The calculations by means of the superposition of correlated configurations method (Hylleraas-CI ), that is, the combination of configuration interaction with the Hylleraas-type correlation factors, needs the effective evaluation of some nontrivial integrals. This series of papers gives the formulas for all types of integrals needed for molecular calculations when Gaussian lobe functions are used as a basis set. The formulas for two-electron integrals are given in the present paper. The preliminary results for two-electron systems are presented.