**J. Toivanen, B.G. Carlsson, J. Dobaczewski, K. Mizuyama,
R.R. Rodriguez-Guzman, P. Toivanen, P. Vesely**

**January 30, 2010**

We report on an implementation of a new method to calculate RPA strength functions with iterative non-hermitian Arnoldi diagonalization method, which does not explicitly calculate and store the RPA matrix. We discuss the treatment of spurious modes, numerical stability, and how the method scales as the used model space is enlarged. We perform the particle-hole RPA benchmark calculations for double magic nucleus and compare the resulting electromagnetic strength functions against those obtained within the standard RPA.

- Introduction
- RPA from linearized TDHF
- Iterative solution of the RPA equations
- Treatment of spurious RPA modes
- Convergence of strength functions
- Convergence as a function of the number of Arnoldi iterations
- Convergence in function of the number of HO shells

- Scaling of iterative solution method
- Summary and conclusions
- Bibliography
- About this document ...

Jacek Dobaczewski 2010-01-30