The computation of a number of the smallest eigenvalues of large and sparse matrices is crucial in various scientific applications, as the Finite Element solution of PDEs, electronic structure calculations or Laplacian of graphs, to mention a few. We propose in this contribution a parallel algorithm that is based on the spectral low-rank modification of a factorized sparse approximate inverse preconditioner (RFSAI) to accelerate the Newton-based iterative eigensolvers. Numerical results onto matrices arising from various realistic problems with size up to 5 million unknowns and 2.2 x 10^8 nonzero elements account for the efficiency and the scalability of the proposed RFSAI-updated preconditioner.