Summary: The fully consistent relativistic continuum random phase approximation (RCRPA) has been constructed in the momentum representation in the first part of this paper. In this part we describe the numerical details for solving the Bethe-Salpeter equation. The numerical results are checked by the inverse energy weighted sum rules in the isoscalar giant monopole resonance, which are obtained from the constraint relativistic mean field theory and also calculated with the integration of the RCRPA strengths. Good agreement between them is achieved. We study the effects of the self-consistency violation, particularly the currents and Coulomb interaction to various collective multipole excitations. Using the fully consistent RCRPA method, we investigate the properties of isoscalar and isovector collective multipole excitations for some stable and exotic from light to heavy nuclei. The properties of the resonances, such as the centroid energies and strength distributions are compared with the experimental data as well as with results calculated in other models.
For part I see [the authors, ibid. 53, No. 4, 716–722 (2010; Zbl 1228.81283)].