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

J. Terasaki
Department of Physics and Astronomy, CB3255, University of North Carolina, Chapel Hill, NC 27599-3255
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996
Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831
Joint Institute for Heavy-Ion Research, Oak Ridge, TN 37831
J. Engel
Department of Physics and Astronomy, CB3255, University of North Carolina, Chapel Hill, NC 27599-3255
M. Bender
Physics Division, Argonne National Laboratory, Argonne, IL 60439
J. Dobaczewski
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996
Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831
Joint Institute for Heavy-Ion Research, Oak Ridge, TN 37831
Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, Poland
W. Nazarewicz
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996
Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831
Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, Poland
M. Stoitsov
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996
Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831
Joint Institute for Heavy-Ion Research, Oak Ridge, TN 37831
Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Science, Sofia 1784, Bulgaria

Abstract
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.