Chebyshev kernel polynomial method for efficient calculation of the quasiparticle random phase approximation response function

Abstract

Efficient and accurate algorithms for the calculation of the multipole response of deformed atomic nuclei are very important tools in nuclear structure, especially for large scale calculations. In this paper we present an implementation of the algorithm based on the expansion of the response function in terms of the Chebyshev polynomials in conjunction with the kernel polynomial method for a very efficient calculation of the quasiparticle random phase approximation response function. Several test calculations are performed in order to asses the applicability and feasibility of this algorithm, already used successfully in the field of computational solid state physics, in various nuclear structure calculations.