Determining the numerical stability of quantum chemistry algorithms

Gerald Knizia, Wenbin Li, Sven Simon, Hans Joachim Werner

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Moller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided.

Original languageEnglish (US)
Pages (from-to)2387-2398
Number of pages12
JournalJournal of Chemical Theory and Computation
Volume7
Issue number8
DOIs
StatePublished - Aug 9 2011

Fingerprint

Quantum chemistry
numerical stability
Convergence of numerical methods
quantum chemistry
Angular momentum
compilers
programming
floating
angular momentum
perturbation theory
injection
evaluation
output
estimates
approximation

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Knizia, Gerald ; Li, Wenbin ; Simon, Sven ; Werner, Hans Joachim. / Determining the numerical stability of quantum chemistry algorithms. In: Journal of Chemical Theory and Computation. 2011 ; Vol. 7, No. 8. pp. 2387-2398.
@article{160da32e8b2c45929025f7a06ddb3cda,
title = "Determining the numerical stability of quantum chemistry algorithms",
abstract = "We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Moller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided.",
author = "Gerald Knizia and Wenbin Li and Sven Simon and Werner, {Hans Joachim}",
year = "2011",
month = "8",
day = "9",
doi = "10.1021/ct200239p",
language = "English (US)",
volume = "7",
pages = "2387--2398",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "8",

}

Determining the numerical stability of quantum chemistry algorithms. / Knizia, Gerald; Li, Wenbin; Simon, Sven; Werner, Hans Joachim.

In: Journal of Chemical Theory and Computation, Vol. 7, No. 8, 09.08.2011, p. 2387-2398.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Determining the numerical stability of quantum chemistry algorithms

AU - Knizia, Gerald

AU - Li, Wenbin

AU - Simon, Sven

AU - Werner, Hans Joachim

PY - 2011/8/9

Y1 - 2011/8/9

N2 - We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Moller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided.

AB - We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Moller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided.

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

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

U2 - 10.1021/ct200239p

DO - 10.1021/ct200239p

M3 - Article

C2 - 26606614

AN - SCOPUS:80051612714

VL - 7

SP - 2387

EP - 2398

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 8

ER -