TY - JOUR

T1 - Molecular integrals by numerical quadrature. I. Radial integration

AU - Lindh, Roland

AU - Malmqvist, Per Åke

AU - Gagliardi, Laura

PY - 2001/7

Y1 - 2001/7

N2 - This article presents a numerical quadrature intended primarily for evaluating integrals in quantum chemistry programs based on molecular orbital theory, in particular density functional methods. Typically, many integrals must be computed. They are divided up into different classes, on the basis of the required accuracy and spatial extent. Ideally, each batch should be integrated using the minimal set of integration points that at the same time guarantees the required precision. Currently used quadrature schemes are far from optimal in this sense, and we are now developing new algorithms. They are designed to be flexible, such that given the range of functions to be integrated, and the required precision, the integration is performed as economically as possible with error bounds within specification. A standard approach is to partition space into a set of regions, where each region is integrated using a spherically polar grid. This article presents a radial quadrature which allows error control, uniform error distribution and uniform error reduction with increased number of radial grid points. A relative error less than 10-14 for all s-type Gaussian integrands with an exponent range of 14 orders of magnitude is achieved with about 200 grid points. Higher angular l quantum numbers, lower precision or narrower exponent ranges require fewer points. The quadrature also allows controlled pruning of the angular grid in the vicinity of the nuclei.

AB - This article presents a numerical quadrature intended primarily for evaluating integrals in quantum chemistry programs based on molecular orbital theory, in particular density functional methods. Typically, many integrals must be computed. They are divided up into different classes, on the basis of the required accuracy and spatial extent. Ideally, each batch should be integrated using the minimal set of integration points that at the same time guarantees the required precision. Currently used quadrature schemes are far from optimal in this sense, and we are now developing new algorithms. They are designed to be flexible, such that given the range of functions to be integrated, and the required precision, the integration is performed as economically as possible with error bounds within specification. A standard approach is to partition space into a set of regions, where each region is integrated using a spherically polar grid. This article presents a radial quadrature which allows error control, uniform error distribution and uniform error reduction with increased number of radial grid points. A relative error less than 10-14 for all s-type Gaussian integrands with an exponent range of 14 orders of magnitude is achieved with about 200 grid points. Higher angular l quantum numbers, lower precision or narrower exponent ranges require fewer points. The quadrature also allows controlled pruning of the angular grid in the vicinity of the nuclei.

KW - Density functional theory

KW - Molecular orbital theory

KW - Numerical quadrature

UR - http://www.scopus.com/inward/record.url?scp=0035590955&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0035590955&partnerID=8YFLogxK

U2 - 10.1007/s002140100263

DO - 10.1007/s002140100263

M3 - Article

AN - SCOPUS:0035590955

VL - 106

SP - 178

EP - 187

JO - Theoretical Chemistry Accounts

JF - Theoretical Chemistry Accounts

SN - 1432-881X

IS - 3

ER -