𝔖 Bobbio Scriptorium
✦   LIBER   ✦

Numerical instabilities in the computation of pseudopotential matrix elements

✍ Scribed by Christoph van Wüllen

Publisher
John Wiley and Sons
Year
2005
Tongue
English
Weight
99 KB
Volume
27
Category
Article
ISSN
0192-8651

No coin nor oath required. For personal study only.

✦ Synopsis

Abstract

Steep high angular momentum Gaussian basis functions in the vicinity of a nucleus whose inner electrons are replaced by an effective core potential may lead to numerical instabilities when calculating matrix elements of the core potential. Numerical roundoff errors may be amplified to an extent that spoils any result obtained in such a calculation. Effective core potential matrix elements for a model problem are computed with high numerical accuracy using the standard algorithm used in quantum chemical codes and compared to results of the MOLPRO program. Thus, it is demonstrated how the relative and absolute errors depend an basis function angular momenta, basis function exponents and the distance between the off‐center basis function and the center carrying the effective core potential. Then, the problem is analyzed and closed expressions are derived for the expected numerical error in the limit of large basis function exponents. It is briefly discussed how other algorithms would behave in the critical case, and they are found to have problems as well. The numerical stability could be increased a little bit if the type 1 matrix elements were computed without making use of a partial wave expansion. © 2005 Wiley Periodicals, Inc., J Comput Chem 27: 135–141 2006

📜 SIMILAR VOLUMES

On the computation of matrix elements be
On the computation of matrix elements between numerical wave functions: The canonical functions method
✍ H. Kobeissi; M. Dagher; A. El-Hajj; M. Kobeissi 📂 Article 📅 1989 🏛 John Wiley and Sons 🌐 English ⚖ 727 KB

The problem of the computation of the matrix elements Z ( v , v ' ; k ) = [ q U ( r ) ( r -rJkqUr(r)dr, is considered when q u ( r ) and q J r ) are eigenfunctions related to a diatomic potential of the RKR type (defined by the coordinates of its turning points Pi with polynomial interpolations). Th

Numerical aspects of finite element simu
Numerical aspects of finite element simulations of residual stresses in metal matrix composites
✍ F. Teixeira-Dias; L. F. Menezes 📂 Article 📅 2001 🏛 John Wiley and Sons 🌐 English ⚖ 352 KB 👁 2 views

When a metal matrix composite (MMC) is cooled down from the fabrication or annealing temperature to room temperature, residual stresses are induced in the composite due to the mismatch of the thermal expansion coe cients of the matrix and reinforcement. A thermomechanical model describing these proc

A Numerical Study Of The Coefficient Mat
A Numerical Study Of The Coefficient Matrix Of The Boundary Element Method Near Characteristic Frequencies
✍ P. Juhl 📂 Article 📅 1994 🏛 Elsevier Science 🌐 English ⚖ 446 KB

In this paper the well-known failure in numerical computations of the exterior surface Helmholtz integral equation at certain characteristic frequencies is investigated. The problem of characteristic frequencies is numerically reflected in an ill-conditioning of the coefficient matrix arising from t