Spline representation. I. Linear spline bases for atomic calculations

摘要

The use of a basis of cardinal splines for atomic calculations with the expansion method is proposed. The basis has the property that the coefficientsCpin the orbital expansionjgr; lpar;rrpar; equals; Sgr;pequals;minus;1Nplus;1thinsp;Cpkhgr;plpar;rrpar;are given byCminus;1equals;jgr;prime;(r0),Cpequals;jgr;(rp),pequals;0,1, middot;middot;middot;,N, andCNplus;1equals;jgr;prime;(rN), wherer0,r1, middot;middot;middot;,rNis the mesh of knots on which the spline is defined. The basis functions khgr;p(r) are continuous with continuous first and second derivatives, and may readily be calculated from the expansionkhgr;plpar;rrpar; equals; Sgr;kequals;1Nthinsp;Sgr;sequals;03thinsp;rsdklpar;rrpar;Sksp, wheredk(r)equals;1 ifrkminus;1le;rrkand 0 otherwise andSkspis a matrix of coefficients which are uniquely determined by the mesh alone and may be constructed by simple algebraic operations, the most complex being the inversion of a single (Nplus;1)times;(Nplus;1) matrix. Atomic Hartreehyphen;Fock matrix elements are given by simple polynomial expressions. The number of twohyphen;electron integrals increases only as the square of the number of basis functions and can be stored in factored form, so that the storage requirements for the twohyphen;electron integrals are trivial. The largest two arrays (equal in size) are the density matrix and the matrix elements of the exchange operator. The results of numerical tests for the hydrogen atom using several different distributions of mesh points are presented. It is found that the number of mesh points required is smaller than for standard numerical methods by a factor of about 4. The accuracy of a linear spline representation withNplus; 1 knots is comparable to that of a Gaussian expansion withNbasis functions.