Self-Consistent Description of Multipole Strength in Exotic Nuclei I: Method

We use the canonical Hartree-Fock-Bogoliubov basis to implement a completely self-consistent quasiparticle-random-phase approximation with arbitrary Skyrme energy density functionals and density-dependent pairing functionals. The point of the approach is to accurately describe multipole strength functions in spherical even-even nuclei, including weakly-bound drip-line systems. We describe the method and carefully test its accuracy, particularly in handling spurious modes. To illustrate our approach, we calculate isoscalar and isovector monopole, dipole, and quadrupole strength functions in several Sn isotopes, both in the stable region and at the drip lines.