An improved algorithm for the normalized elimination of the small-component method

A new algorithm for the iterative solution of the normalized elimination of the small component (NESC) method is presented that is less costly than previous algorithms and that is based on (1) solving the NESC equations for the uncontracted rather than contracted basis (First-Diagonalize-then-Contract), (2) a new iterative procedure for obtaining the NESC Hamiltonian (iterative TU algorithm), (3) the renormalization scheme connected to the picture change, and (4) a finite nucleus model with a Gaussian charge distribution. The accuracy of NESC energies, which match those of 4-component Dirac calculations, is demonstrated. Test calculations with CCSD(T), DFT, and large basis sets including high angular momentum basis functions (f,g,h,i) are presented to prove the general applicability of the new NESC algorithm. Comparison with other algorithms of solving the NESC equations are shortly discussed and time savings are presented.