the study of nuclei far from stability is an increasingly important part of nuclear physics @xcite .
as radioactive beams allow more experiments on these nuclei , theoretical modeling is changing in significant ways .
new ideas and progress in computer technology have allowed nuclear theorists to understand bits and pieces of nuclear structure quantitatively @xcite .
short - lived exotic nuclei offer unique tests of those aspects of our developing many - body theories that depend on neutron excess  @xcite .
the major challenge is to predict or describe in detail exotic new properties of nuclei far from the stability valley , and to understand the origins of these properties .    for medium - mass and heavy nuclei
, an important goal is obtaining a universal energy - density functional , which will be able to describe static and dynamic properties of finite nuclei and extended nucleonic matter with arbitrary neutron - to - proton ratio .
self - consistent methods based on density - functional theory are already sophisticated enough to allow precise analysis of ground - state properties ( e.g.  binding energies ) in heavy nuclei @xcite .
they can also help describe nuclear decays and excited states .
their predictions for collective excitations as we approach the neutron drip line are especially interesting .
what happens to low- and high - frequency multipole modes when the neutron excess is unusually large ?    to address these questions we use the quasiparticle random - phase approximation ( qrpa ) , a powerful tool for understanding both low - lying vibrational states and giant resonances @xcite .
the qrpa is a microscopic approach that is nevertheless simple enough to allow `` no - core '' calculations .
the approximation , which should be good for collective vibrations as long as their amplitudes are small , is especially effective in conjunction with skyrme energy functionals .
our work is a part of a broad program to test and improve these functionals , which thus far have been fitted mainly to ground - state observables , by applying them to collective excitations , particularly near the drip line .
this paper lays out our approach and evaluates its accuracy . for these purposes
we restrict ourselves to a single skyrme functional , skm@xmath0 .
a forthcoming study will examine the performance of skyrme functionals more generally .
the qrpa is a standard method for describing collective excitations in open - shell superconducting nuclei with stable mean - field solutions , either spherical or deformed .
what is not standard , and at the same time is extremely important for weakly bound nuclei , is the treatment of the particle continuum .
continuum extensions of the random phase approximation ( rpa ) or qrpa are usually carried out in coordinate space , facilitating treatment of decay channels and guaranteeing correct asymptotics .
surprisingly , as we discuss below , the rich literature on the rpa and qrpa , which includes many coordinate - space calculations , contains few treatments of the continuum that exploit the entire skyrme functional in a fully self - consistent way .    to avoid confusion
, we state what we mean by a fully self - consistent rpa or qrpa calculation . first
, the underlying mean - field calculation must be self - consistent in the usual sense .
next , the residual interaction used in the rpa or qrpa must be derived from the same force or energy functional that determines the mean field .
an important consequence of this condition , and of other more detailed technical conditions discussed below , is that spurious excitations arising from symmetry breaking by the mean field have zero or nearly zero energy , leaving the physical intrinsic excitations completely uncontaminated by spurious motion .
finally , energy - weighted sum rules must be satisfied to high accuracy .
we elaborate on these requirements below ; refs .
@xcite discuss ways in which rpa calculations commonly violate them .
the literature applying rpa or qrpa to nuclear structure is huge , and a complete review is beyond the scope of our paper .
we do , however , present an overview of the studies that are related in one way or another to nuclear density functionals , self consistency , pairing , and the key issue of the particle continuum .
the standard version of qrpa , the so - called matrix formulation , is carried out in the configuration space @xcite of single - quasiparticle states .
a number of papers treat collective states in spherical nuclei in the skyrme - rpa and qrpa matrix formulation ( see refs . @xcite and references cited therein ) , in which the positive - energy continuum is discretized , e.g.  by solving the hartree - fock - bogoliubov ( hfb ) and qrpa equations in a harmonic - oscillator single - particle basis . within this group , the first fully self - consistent calculations that properly account for continuum effects are those of refs .
@xcite , in which the localized canonical basis of coordinate - space hfb is used to calculate beta - decay rates of neutron - rich r - process nuclei and gamow - teller strength distributions .
recently , fully self - consistent hfb+qrpa calculations have also been carried out with the finite - range gogny force @xcite . unlike many previous gogny+hfb studies that employed a harmonic oscillator basis , ref .
@xcite solves the hfb equations in the eigenbasis of a woods - saxon potential , the particle continuum of which is discretized by enclosing the system in a box .
coordinate - space green s functions as a method of implementing the rpa through linear response were first used in ref .
@xcite and subsequently applied to the description of low- and high - energy nuclear modes ( see , e.g. , refs .
many of those calculations are not realistic enough , however , because they ignore the spin - orbit and coulomb residual interactions in the rpa @xcite .
coordinate - space greens - function qrpa was studied in ref .
@xcite , in the bcs approximation , with a phenomenological woods - saxon average potential .
coordinate - space hfb+qrpa for spherical nuclei was formulated in refs .
@xcite and applied to excitations of neutron - rich nuclei .
as in @xcite , the hartree - fock ( hf ) field in refs .
@xcite was approximated by a woods - saxon potential .
while the calculations of refs .
@xcite are based on skyrme - hfb fields , they violate full self consistency by replacing the residual velocity - dependent terms of the skyrme force by the landau - migdal force in the qrpa , and neglecting spin - spin , spin - orbit , and coulomb residual interactions entirely . within this approach ,
extensive skyrme - hf+bcs qrpa calculations of e1-strength in neutron - rich nuclei were carried out in refs .
@xcite .
an alternative coordinate - representation approach , also based on green s functions , was formulated in refs .
@xcite within migdal s finite - fermi - systems theory .
most of practical applications of this method , however , involve approximations that break self consistency in one way or another , including the use of highly truncated pairing spaces , different interactions in hfb and qrpa , and the so - called diagonal pairing approximation @xcite .
properties of excited states and strength functions have also been investigated within the relativistic rpa @xcite or qrpa @xcite .
the qrpa work employs the matrix formulation and is fully self - consistent , since it uses the same lagrangian in the relativistic hartree - bogoliubov calculation of the ground state and in the qrpa matrix equations , which are solved in the canonical basis .    at the present ,
no fully self - consistent continuum hfb+qrpa calculations exist in deformed nuclei .
@xcite studied giant resonances in deformed nuclei within time - dependent hf theory , formulated in coordinate space with a complex absorbing boundary condition imposed .
symmetry - unrestricted rpa calculations , with no pairing , were carried out in ref .
@xcite in a `` mixed representation '' @xcite on a cartesian mesh in a box , while ref .
@xcite contains examples of bcs+qrpa calculations in the single - particle basis of a deformed woods - saxon potential .
the work described in this paper is fully self - consistent : among other things we use precisely the same interaction in the hfb and qrpa calculations so as to preserve the small - amplitude limit of time - dependent hfb .
we formulate the qrpa in the canonical eigenbasis of the one - body particle - density matrix @xcite which is calculated in the coordinate representation in a large spherical box .
as mentioned above , the canonical basis has been used previously to study @xmath1 decay and gamow - teller strength @xcite ; its use in charge - conserving modes near the drip line is more challenging , however , because of the existence of spurious states in the monopole and dipole channels .
these zero - energy modes can mix with physical states unless the qrpa equations are solved with high accuracy .
a less precise implementation of our approach was used to calculate neutrino - nucleus cross sections in @xmath2pb in ref .
@xcite .
this paper is organized as follows .
section  [ sec : method ] below presents our approach . in sec .
[ sec : accuracy ] we check the qrpa solutions carefully , focusing on spurious modes .
section [ sec : conclusion ] contains the main conclusions of our work .
mathematical details are in two appendices , the first of which is on the qrpa equations and the second on calculating the derivatives of the skyrme functionals that enter the formalism .
our first step in the self - consistent treatment of excitations is to solve the spherical hfb equations in coordinate space ( without mixing neutron and proton quasiparticle wave functions @xcite ) , with the method developed in ref .
@xcite ( see also refs .
we can use arbitrary skyrme functionals in the particle - hole and pairing ( particle - particle ) channels .
we modify the code used in refs .
@xcite so that it solves the hfb equations with higher accuracy , which we need because the qrpa uses all the single - quasiparticle states produced by the hfb equations , even those that are essentially unoccupied .
our modifications are : ( i ) the use of quadruple precision ( though in solving the qrpa equations we use double precision ) ; ( ii ) a smaller discretization length ( 0.05 fm ) ; and ( iii ) a high quasiparticle - energy cutoff ( 200 mev ) and a maximum angular momentum @xmath3=15/2 ( @xmath4 ) or 21/2 ( @xmath5 ) . in a 20 fm box , this cutoff corresponds to 200300 quasiparticle states for each kind of nucleon .
we include all these quasiparticle states in the hfb calculation because a very large energy cutoff is essential for the accuracy of self - consistent qrpa calculations @xcite .
hence , the effective pairing window in our hfb calculations is also very large , with the pairing functional fitted to experimental pairing gaps extracted as in ref .
@xcite from the measured odd - even mass differences in several sn , ni , and ca isotopes .
next , we construct the canonical basis , the eigenstates of single - particle density matrix @xmath6 . to avoid poor accuracy
( see ref .
@xcite ) in the wave functions of the nearly empty canonical particle states , we do not diagonalize @xmath6 directly in coordinate space .
instead we construct an intermediate basis by orthonormalizing a set of functions @xmath7 , where @xmath8 and @xmath9 are the upper and lower components of the quasiparticle wave function with energy @xmath10 @xcite .
we use the density matrix in coordinate space to calculate the matrix in this basis , which we then diagonalize to obtain the canonical states .
the reason for using the sum of @xmath8 and @xmath11 is that solutions of the hfb equations expressed in the canonical basis ( eqs .
( 4.14 ) of ref .
@xcite ) are , in the new basis , _ guaranteed _ to be numerically consistent with those of the original hfb problem .
this is because the configuration space is the same in both cases , independent of the pairing cutoff ( see ref .
@xcite for a discussion relevant to this point ) . without pairing , when either @xmath12 or @xmath9 is equal to zero , our method is equivalent to taking a certain number of hf states , including many unoccupied states .    in the canonical basis ,
the hfb+qrpa equations have a form almost identical to that of the bcs+qrpa approximation , the only difference being the presence of off - diagonal terms in the single - quasiparticle energies .
the qrpa+hfb formalism employs more pairing matrix elements than the qrpa+bcs , however .    as noted already , full self consistency requires the use of the same interaction in the qrpa as in the hfb approximation . more specifically , this means that the matrix elements that enter the qrpa equation are related to second derivatives of a mean - field energy functional .
we describe the densities and the form of the functional carefully in the appendices .
but we must meet other conditions as well for qrpa calculation to be self - consistent .
essentially all the single - particle or quasiparticle states produced by the hfb calculation must be used in the space of two - quasiparticle qrpa excitations .
this requirement is rather stringent , so we truncate the two - quasiparticle space at several levels and check for convergence of the qrpa solution
. first we omit canonical - basis wave functions that have occupation probabilities @xmath13 less than some small @xmath14 , ( or hf energies greater than some @xmath15 if there is no pairing ) .
then we exclude from the qrpa pairs of canonical states for which the occupation probabilities are both larger than @xmath16 .
this second cut is based on the assumption that two - particle transfer modes are not strongly coupled to particle - hole excitations .
in addition , if the factors containing @xmath17 and @xmath18 in the qrpa equation  see eqs .
( [ a_spherical ] ) and ( [ b_spherical ] )  are very small , in practice smaller than 0.01 , then we set the corresponding matrix elements equal to zero .
this does not affect the size of the qrpa space , but significantly speeds up the calculations . for good performance
we diagonalize qrpa - hamiltonian matrices of order @xmath19 in neutron - rich sn isotopes .
having solved the qrpa equations , we can then calculate the strength function @xmath20 for the multipole operator @xmath21 .
the smoothing width @xmath22 is supposed to be large enough to remove spurious oscillations in @xmath23 associated with a finite box radius @xmath24 @xcite .
a reasonable form , based on a single - particle estimate , for the smoothing width ( app .
b of ref .
@xcite ) , is @xmath25 where @xmath26 is the neutron fermi level and @xmath27 is the nucleon mass . in deriving eq .
( [ smwidth ] ) we assumed that the single - proton continuum is effectively shifted up several mev by the coulomb barrier .
in other words , we associate the threshold energy with the neutron fermi level .    in all the tests below
, we use the skyrme functional skm@xmath28 @xcite and a volume pairing functional @xcite ( @xmath29 a constant in eq .
( [ pairing_energy_functional ] ) ) .
the pairing parameter in eq .
( [ volpair ] ) is @xmath30 mev @xmath31 .
usually we work in a box of radius 20 fm , though we vary this radius below to see its effects . in several tests
we examine the weakly bound nucleus @xmath32sn , which is very close to the two - neutron drip line . in this system ,
the protons are unpaired and the neutrons paired ( with @xmath33=1.016mev ) in the hfb ground state .
neutron quasiparticle energies @xmath34 for @xmath35 states in @xmath32sn , calculated by diagonalizing the hfb hamiltonian in the canonical basis , versus the quasiparticle energies @xmath36 obtained by directly solving the hfb equations in coordinate space .
standard ( solid line , diamonds ) and improved ( dotted line , dots ) methods are used to obtain the canonical states .
see text for details . , width=453 ]    benchmark tests of the hfb part of our calculations are reported in ref .
since the accuracy of the canonical wave functions , in which the qrpa calculations are carried out , strongly affects the quality of results ( in particular qrpa self consistency ) , we take special care to compute them precisely . as discussed in sec .
[ sec : method ] , we obtain canonical states by diagonalizing the single - particle density matrix @xmath6 represented in the orthonormalized set of functions @xmath37 .
the accuracy of this method is illustrated in fig .
[ fig : ehfb_eigh11 ] , which plots the quasiparticle energies @xmath34 , obtained by diagonalizing the hfb hamiltonian in the canonical basis ( eq . ( 4.20 ) of ref .
@xcite ) , versus the quasiparticle energies @xmath36 obtained by solving the hfb differential equations directly in coordinate space ( eq . ( 4.10 ) of ref .
two sets of canonical states are used : ( i ) those obtained through the procedure outlined above ( dotted line ) and ( ii ) those obtained in the standard way by diagonalizing the density matrix @xmath38 in discretized coordinate space ( eq .
( 3.24a ) of ref .
@xcite ; solid line ) .
if the canonical basis is precisely determined , @xmath34=@xmath36 and the two sets of @xmath34 coincide . within the standard approach , however , the high canonical energies deviate visibly from their hfb counterparts , i.e. , the accuracy of the underlying canonical wave functions is poor . on the other hand ,
the quasiparticle energies and canonical wave functions calculated within the modified approach introduced above are as accurate as the original solutions to the hfb equations , even for high - lying nearly - empty states .
( see also sec .
vi.d of ref .
@xcite for a discussion relevant to this point . )    having examined the canonical basis , we turn to the accuracy of the qrpa part of the calculation . to test it , we first consider solutions related to symmetries .
if a hamiltonian is invariant under a symmetry operator @xmath39 and the hfb state @xmath40 spontaneously breaks the symmetry , then @xmath41 , with @xmath42 an arbitrary c - number , is degenerate with the state @xmath40 .
the qrpa equations have a spurious solution at zero energy associated with the symmetry breaking @xcite , while all other solutions are free of the spurious motion .
this property is important for strength functions , and gives us a way of testing the calculations . since our qrpa equations ,
which assume spherical symmetry , are based on mean fields that include pairing and are localized in space , there appear spurious states associated with particle - number nonconservation ( proton and/or neutron ; 0@xmath43 channel ) and center - of - mass motion ( 1@xmath44 channel ) .
these two cases are discussed below in sec .
[ zeroplus ] and sec .
[ oneminus ] .
in addition to the spurious state associated with nonconservation of particle - number by the hfb , the @xmath45 channel contains the important `` breathing mode '' . in table  [ tab : strength_number ] we display results from a run with @xmath46=@xmath47 for neutrons and @xmath48 mev for protons , resulting in the inclusion of 310 proton quasiparticle states and the same number of neutron states , with angular momentum up to @xmath49 .
the table shows the qrpa energies and transition matrix elements of the particle - number operator .
the spurious state is below 200 kev , well separated from the other states , all of which have negligible `` number - strength '' .
the nonzero number strength in the spurious state , like the nonzero energy of that state , is a measure of numerical error . if the space of two - quasiparticle states is smaller , with @xmath50 mev and @xmath51 , the energy of the spurious state and the number strength barely change .
isoscalar @xmath45 strength function in @xmath32sn for ( i ) the single - proton energy cutoff @xmath15=100 mev and the neutron - quasiparticle occupation cutoff @xmath52 ( thin solid line ) ; ( ii ) @xmath15=150 mev and @xmath53 ( dotted line ) ; and ( iii ) @xmath15=200 mev and @xmath54 ( thick solid line ) .
results corresponding to ( ii ) and ( iii ) practically coincide.,width=453 ]    figure [ fig:0+_sn174_ecrit ] shows the strength function @xmath23 for the isoscalar @xmath45 transition operator , cf .
@xcite , @xmath55 we have plotted three curves with successively more quasiparticle levels ( from 246 proton levels and 203 neutron levels to 341 proton levels and 374 neutron levels ) , with cutoff parameters given in the figure caption .
the major structures in the strength function are stable .
the error remaining after to the gentlest truncation is extremely small .
the dependence of the strength function on the box size and quasiparticle cutoff is shown in fig .
[ fig:0+_sn174_box ] .
the upper part of the figure ( panels a - c ) corresponds to a constant smoothing width of @xmath22=0.5 mev .
this relatively small value is not sufficient to eliminate the finite - box effects but it allows us to assess the stability of the qrpa solutions as a function of @xmath24 .
the large structure corresponding to the giant monopole resonance ( gmr ) is independent of box size no matter what the cutoff , but increasing the number of configurations magnifies the dependence on box size of local fluctuations in @xmath23 . the lower part of the figure ( panels d - f ) are smoothed more realistically , as in eq .
( [ smwidth ] ) .
it is gratifying to see that the resulting strength functions are practically identical , i.e. , the remaining dependence on @xmath24 and the cutoff is very weak .
isoscalar @xmath45 strength function in @xmath32sn for the box radii : @xmath24=20fm ( solid line ) and @xmath24=25fm ( dotted line ) . in ( a ) , ( b ) , and ( c ) the smoothing - width parameter @xmath22 is 0.5 mev for all energies , while in ( d ) , ( e ) , and ( f ) @xmath56 is given by eq .
( [ smwidth ] ) .
we use the same three sets of cutoff conditions as in fig .
[ fig:0+_sn174_ecrit ] , namely ( i ) in parts ( a ) and ( d ) , ( ii ) in parts ( b ) and ( e ) , and ( iii ) in parts ( c ) and ( f ) . , width=566 ]     isoscalar @xmath57 strength function in @xmath58sn for the corrected dipole operator in eq .
( [ dipolecorr ] ) ( solid line ) and the uncorrected operator in eq .
( [ dipole ] ) ( dotted line ) . the cutoff @xmath15 is 140 mev and @xmath59 is @xmath60 .
the self consistency of our calculations makes the solid and dotted curves coincide nearly exactly over the whole energy range.,width=491 ]     isoscalar @xmath57 strength function in @xmath61sn ( left ) and @xmath32sn ( right ) for two box radii : @xmath24=20fm ( solid line ) and @xmath24=25fm ( dotted line ) . in ( a ) and ( b )
the smoothing - width parameter is constant ( @xmath22=0.5 mev ) , while in ( c ) and ( d ) @xmath56 is given by eq .
( [ smwidth ] ) .
, width=566 ]    isoscalar and isovector strength functions for ( a ) the 0@xmath43 channel of @xmath61sn , ( b ) the 0@xmath43 channel of @xmath32sn , ( c ) the 2@xmath43 channel of @xmath61sn , and ( d ) the 2@xmath43 channel of @xmath32sn .
the cutoff @xmath62 is 150 mev and @xmath59 is @xmath47 .
, width=529 ]    the energy - weighted sum rule ( ewsr ) for the isoscalar @xmath45 mode @xcite is given by @xmath63 where the expectation value is evaluated in the hfb ground state .
this sum rule provides a stringent test of self consistency in the qrpa . in @xmath32sn , the right - hand side of eq.([isoscalar_0+_sum_rule ] ) is 35215 @xmath64  mev  @xmath65 and the left - hand side 34985@xmath66 @xmath64  mev  @xmath65 for all of the calculations of fig .
[ fig:0+_sn174_box ] ; the qrpa strength essentially exhausts the sum rule .
( the qrpa values of the ewsr in this paper are obtained by summing up to @xmath67  mev . )
the @xmath57 channel , home of the giant dipole resonance , the isoscalar squeezing resonance , and as yet incompletely understood low - energy peaks in neutron - rich nuclei ( sometimes associated with skin excitations ) , has a spurious isoscalar mode associated with center - of - mass motion that can seriously compromise the low - energy spectrum if not handled with extreme care .
we test the ability of our qrpa to do so in @xmath68sn , @xmath61sn , @xmath32sn , and @xmath69sn .
( the nuclei @xmath68sn and @xmath69sn are the two - proton and two - neutron drip - line systems predicted by the hfb calculation with skm@xmath70 .
neither nucleus has any static pairing , i.e. , @xmath33=@xmath71=0 . ) in the following calculations , we take @xmath72 mev for the protons and @xmath73 for the neutrons . as discussed above ,
smoothed strength functions are practically independent of small changes in the cutoff .
they are also independent of the cutoff in quasiparticle angular momentum provided we include all states with @xmath74@xmath7515/2 .
figure [ fig:1-_sn100_sn176 ] shows the predicted isoscalar dipole strength function for @xmath58sn . for the transition operator
, we use @xmath76 and the corrected operator @xmath77 to remove as completely as possible residual pieces of the spurious state from the physical states @xcite .
the fact that the strength functions produced by these two operators  displayed in fig .
[ fig:1-_sn100_sn176 ]  coincide so closely shows the extreme accuracy of our qrpa solutions ; they are uncontaminated by spurious motion even without the operator correction .
the spurious - state energies @xmath78 are 0.964 mev for @xmath68sn and 0.713 mev for @xmath61sn , and the energies of the first physical excited states are 7.958 mev for @xmath68sn and 7.729 mev for @xmath61sn . in @xmath32sn ( @xmath69sn ) , @xmath78 is 0.319 mev ( 0.349 mev ) and the first physical state is at 3.485 mev ( 2.710 mev ) , lower than in the more stable isotopes .
pairing correlations do not affect accuracy ; the neutrons in @xmath61sn and @xmath32sn are paired , while those in @xmath68sn and @xmath69sn are not .
we display the fine structure of the isoscalar @xmath57 strength functions in @xmath61sn and @xmath32sn in fig .
[ fig:0-_sn174_box ] , which also illustrates the dependence of the results on @xmath24 .
the dependence is consistent with that of fig .
[ fig:0+_sn174_box ] for the isoscalar 0@xmath43 strength ; the low - amplitude fluctuations in @xmath23 that are unstable as a function of @xmath24 disappear , and the smoothed strength function depends only weakly on @xmath24 . in @xmath61sn ,
the two sharp peaks below 10 mev correspond to discrete states while the broad maxima centered around 15 mev and 27 mev are in the continuum , well above neutron - emission threshold .
a similar three - peaked structure emerges in @xmath32sn , though most of the strength there is concentrated in the low - energy peak at @xmath79 mev .
[ fig:1-_sn100_sn176 ] shows ( as we will discuss in our forthcoming paper @xcite ) that the appearance of the low - energy isoscalar dipole strength is a real and dramatic feature of neutron - rich dripline nuclei @xcite .
the ewsr for the isoscalar @xmath57 mode @xcite is @xmath80 in @xmath32sn , the right - hand side is 403310 @xmath64  mev  @xmath81 , while the left - hand side is 400200 @xmath64  mev  @xmath81 . for @xmath69sn ,
the corresponding numbers are 406576 @xmath64  mev  @xmath81 and 407100 @xmath64  mev  @xmath81 .
this level of agreement is very good .
figure  [ fig : is_iv_0+_2 + ] displays strength functions for the 0@xmath43 and 2@xmath43 channels in @xmath61sn and @xmath32sn .
( we discussed the isoscalar @xmath45 mode above to illustrate the accuracy of our solutions , but include it here as well for completeness . )
the calculations show the appearance of low - energy @xmath45 strength  both isovector and isoscalar  and low - energy isovector 2@xmath43 strength in @xmath32sn , though in none of these instances is the phenomenon quite as dramatic as in the isoscalar @xmath57 channel .
the ewsr for the isoscalar 2@xmath43 transition operator , @xmath82 can be written as @xcite @xmath83 the sum rule is obeyed as well in the 2@xmath43 isoscalar channel as in the @xmath45 and @xmath57 channels , the only difference being that one needs to include quasiparticle states with @xmath84 for @xmath32sn . for @xmath61sn ( @xmath32sn ) from fig .  [
fig : is_iv_0+_2 + ] , the ewsr is 37222 ( 34971 ) @xmath64mev@xmath65 while the qrpa value is 37030 ( 35010 ) @xmath64mev@xmath65 .    while on the topic of the sum rule , we display in table  [ tab : jmax_dependence_of_ewsr ] the @xmath85-dependence of the ewsr for several channels in @xmath86sn , with @xmath87 fm . by taking @xmath88
we appear to obtain essentially the entire strength in all three cases .
in this work we have reported on the development and detailed testing of a fully self - consistent skyrme - qrpa framework that employs the canonical hfb basis . the method can be used to calculate strength distributions in any spin - isospin channel and in any spherical even - even nucleus .
a good calculation requires a large single - quasiparticle space .
our results show that our space is large enough in nuclei as heavy as the sn isotopes .
we are currently investigating the predictions of a range of skyrme functionals across the ca , ni , and sn isotope chains . the initial results presented here point to increases in low - lying strength at the neutron drip line , particularly in the isoscalar - dipole channel . in a forthcoming paper @xcite
we will report on the robustness of these effects , on the physics underlying them , on their variation with atomic mass and number , and on their implications for the future of skyrme functionals .
we gratefully acknowledge useful discussions with shalom shlomo , nils paar , dario vretenar , and the japanese members of the us - japan cooperative project on  mean - field approach to collective excitations in unstable medium - mass and heavy nuclei " .
mb is grateful for the warm hospitality at the physics division at ornl and the theory division at gsi darmstadt , where most of his contribution to this work was made .
this work was supported in part by the u.s .
department of energy under contracts nos .
de - fg02 - 97er41019 ( university of north carolina at chapel hill ) , de - fg02 - 96er40963 ( university of tennessee ) , de - ac05 - 00or22725 with ut - battelle , llc ( oak ridge national laboratory ) , de - fg05 - 87er40361 ( joint institute for heavy ion research ) , and w-31 - 109-eng-38 ( argonne national laboratory ) ; by the national science foundation contract no .
0124053 ( u.s .- japan cooperative science award ) ; by the polish committee for scientific research ( kbn ) ; and by the foundation for polish science ( fnp ) .
we used the parallel computers of the center for computational sciences at oak ridge national laboratory and information technology services at the university of north carolina at chapel hill .
the qrpa equations are the small - oscillations limit of the time - dependent hartree - fock - bogoliubov approximation , see , e.g. ,  @xcite . in the canonical basis
the most general equations take the form @xmath89
@xmath90 @xmath91 @xmath92 } { \delta\rho_{k^\prime k}\delta\rho_{l^\prime l } } , \label{ph_matrix_element_uncoupled}\ ] ] @xmath93 } { \delta\kappa^\ast_{k^\prime k }   \delta\kappa_{l^\prime l } } , \label{pp_matrix_element_uncoupled}\ ] ] @xmath94 } { \delta\kappa^\ast_{k^\prime k }   \delta\rho_{l l^\prime } } = { \bar{v}^{\rm 1p3h\ \ast}}_{l l^\prime k^\prime k } , \label{3p1h_matrix_element_uncoupled}\ ] ] where @xmath95 and @xmath96 are single - particle indices for the canonical basis , and the states are assumed to be ordered .
the symbol @xmath97 refers to the conjugate partner of @xmath95 , @xmath98 and @xmath99 come from the bcs transformation associated with the canonical basis , and the @xmath100 are the one - quasiparticle matrix elements of the hfb hamiltonian ( cf .
( 4.14b ) of ref .
@xcite ) . @xmath101 and @xmath102 are the forward and backward amplitudes of the qrpa solution @xmath103 , and @xmath104 is the corresponding excitation energy .
@xmath105 $ ] is the energy functional ( see app.[app : matrix_elements ] for an explicit definition ) and @xmath6 and @xmath106 are the density matrix and pairing tensor . after taking the functional derivatives , we replace @xmath6 and @xmath106 by their hfb solutions , in complete analogy with an ordinary taylor - series expansion .
to write the equations in coupled form , we introduce the notation @xmath107 where @xmath108 denote spherical quantum numbers . using ( i )
rotational , time - reversal , and parity symmetries of the hfb state , ( ii ) the conjugate single - particle state where @xmath109 and @xmath110 are real radial and spin wave functions .
thus , we multiply all hfb @xmath111 by @xmath112 in the qrpa calculations . ]
@xmath113    and ( iii ) the relations @xmath114j_k } \times\left\{\begin{array}{cc}\sqrt{2 } , & \mu=\mu^\prime,\\ 1 , & \mbox{otherwise},\end{array}\right .
\\ & & y^k_{kk^\prime } = ( -)^{j_\mu - m_\mu}(-)^{j_{\mu^\prime}-m_{\mu^\prime } } \langle j_\mu -m_\mu j_{\mu^\prime } -m_{\mu^\prime}|j_k m_k\rangle \bar{y}^k_{[\mu\mu^\prime]j_k } \times\left\{\begin{array}{cc}\sqrt{2 } , & \mu=\mu^\prime,\\ 1 , & \mbox{otherwise},\end{array}\right.\end{aligned}\ ] ] with @xmath115 the angular momentum of the state @xmath103 and the factor @xmath116 for convenience @xcite , one can rewrite the qrpa equation as    @xmath117j_k,[\nu\nu^\prime]j_k } & b_{[\mu\mu^\prime]j_k,[\bar{\nu}\bar{\nu}^\prime]j_k } \\
-b^\ast_{[\mu\mu^\prime]j_k,[\bar{\nu}\bar{\nu}^\prime]j_k } & -a^\ast_{[\mu\mu^\prime]j_k,[\nu\nu^\prime]j_k } \end{array } \right ) \left ( \begin{array}{c } \bar{x}^k_{[\nu\nu^\prime]j_k}\\ \bar{y}^k_{[\nu\nu^\prime]j_k } \end{array } \right ) = e_k \left ( \begin{array}{c } \bar{x}^k_{[\mu\mu^\prime]j_k}\\ \bar{y}^k_{[\mu\mu^\prime]j_k } \end{array } \right),\ \ \mu \leq
\mu^\prime , \label{qrpa_eq}\ ] ]    @xmath118j_k,[\nu\nu^\prime]j_k } & = & \frac{1}{\sqrt{1+\delta_{\mu\mu^\prime } } } \frac{1}{\sqrt{1+\delta_{\nu\nu^\prime } } } \left\ { e_{\mu\nu}\:\delta_{\mu^\prime\nu^\prime } -e_{\mu^\prime\nu}\:\delta_{\mu\nu^\prime}(-)^{j_\mu+j_{\mu^\prime}-j_k } \right
. \nonumber\\ & & -e_{\mu\nu^\prime}\:\delta_{\mu^\prime\nu}(-)^{j_\mu+j_{\mu^\prime}-j_k } + e_{\mu^\prime\nu^\prime}\:\delta_{\mu\nu}\nonumber\\ & & + g(\mu\mu^\prime\nu\nu^\prime;j_k ) ( u_{\mu^\prime } u_\mu u_\nu u_{\nu^\prime } + v_{\nu } v_{\nu^\prime } v_{\mu^\prime } v_{\mu } ) \nonumber\\ & & + f(\mu\mu^\prime\nu\nu^\prime;j_k ) ( u_{\mu } v_{\nu^\prime } u_\nu v_{\mu^\prime } + u_{\mu^\prime } v_{\nu } u_{\nu^\prime } v_{\mu } ) \nonumber\\ & &   -(-)^ { j_{\nu^\prime}+j_\nu - j_k } f(\mu\mu^\prime\nu^\prime\nu;j_k ) ( u_{\mu } v_{\nu } u_{\nu^\prime } v_{\mu^\prime } + u_{\mu^\prime } v_{\nu^\prime } u_{\nu } v_{\mu } ) \nonumber\\ & & -h(\mu\mu'\nu\nu';j_k ) ( v_\nu v_{\nu ' } u_\mu v_{\mu ' } + u_{\mu ' } v_\mu u_\nu u_{\nu ' } ) \nonumber \\ & & + ( -)^{j_\mu+j_{\mu'}-j_k } h(\mu'\mu\nu\nu';j_k ) ( v_\nu v_{\nu ' } u_{\mu ' } v_\mu + u_\mu v_{\mu ' } u_\nu u_{\nu ' } ) \nonumber \\ & & -h^\ast(\nu\nu'\mu\mu';j_k ) ( u_\mu u_{\mu ' } u_{\nu ' } v_\nu + u_\nu v_{\nu ' } v_{\mu ' } v_\mu ) \nonumber \\ & & + ( -)^{j_\nu+j_{\nu'}-j_k}h^\ast(\nu'\nu\mu\mu';j_k ) ( u_\mu u_{\mu ' } u_\nu v_{\nu ' } + u_{\nu ' } v_\nu v_{\mu ' } v_\mu ) \left.\right\ } , \label{a_spherical}\end{aligned}\ ] ]    @xmath119j_k,[\bar{\nu}\bar{\nu}^\prime]j_k } & = & \frac{1}{\sqrt{1+\delta_{\mu\mu^\prime } } } \frac{1}{\sqrt{1+\delta_{\nu\nu^\prime } } } \left\ {   -g(\mu\mu^\prime\nu\nu^\prime;j_k ) ( u_{\mu^\prime } u_\mu v_\nu v_{\nu^\prime } + u_{\nu } u_{\nu^\prime } v_{\mu^\prime } v_{\mu})\right .
\nonumber\\ & & -(-)^ { j_\nu+j_{\nu^\prime}-j_k } f(\mu\mu^\prime\nu^\prime\nu;j_k ) ( u_{\mu } u_{\nu } v_{\mu^\prime } v_{\nu^\prime } + u_{\mu^\prime } u_{\nu^\prime } v_{\nu } v_{\mu } ) \nonumber\\ & & + ( -)^ { j_\nu+j_{\nu^\prime}+j_\mu+j_{\mu^\prime } } f(\mu^\prime\mu\nu^\prime\nu;j_k ) ( u_{\mu^\prime } u_{\nu } v_{\nu^\prime } v_{\mu } + u_{\mu } u_{\nu^\prime } v_{\mu^\prime } v_{\nu } ) \nonumber\\ & & + h(\mu\mu'\nu\nu';j_k)(v_\nu v_{\nu ' } u_{\mu ' } v_\mu + u_\mu v_{\mu ' } u_\nu u_{\nu ' } ) \nonumber \\ & & -(-)^{j_\mu+j_{\mu'}-j_k } h(\mu'\mu\nu\nu';j_k ) ( v_\nu v_{\nu ' } u_\mu v_{\mu ' } + u_{\mu ' } v_\mu u_\nu u_{\nu ' } ) \nonumber\\ & & -h^\ast(\nu\nu'\mu\mu';j_k ) ( v_\mu v_{\mu ' } u_{\nu ' } v_\nu + u_\nu v_{\nu ' } u_\mu u_{\mu ' } ) \nonumber\\ & & + ( -)^{j_\nu + j_{\nu'}-j_k}h^\ast(\nu'\nu\mu\mu';j_k ) ( v_\mu v_{\mu ' } u_\nu v_{\nu ' } + u_{\nu ' } v_\nu u_\mu u_{\mu ' } ) \left.\right\ } , \label{b_spherical}\end{aligned}\ ] ]    @xmath120j_k | \bar{v}^{\rm pp }    f(\mu\mu^\prime\nu\nu^\prime;j_k)&=&\sum _ { m_\mu m_{\mu^\prime}m_\nu m_{\nu^\prime } } \langle j_\mu m_\mu j_{\mu^\prime } m_{\mu^\prime } | j_k m_k\rangle \langle j_\nu m_\nu j_{\nu^\prime } m_{\nu^\prime } | j_k m_k\rangle \bar{v}^{\rm ph}_{k\bar{l^\prime}\bar{k^\prime}l } \nonumber\\ & = & \sum_{j^\prime } ( -)^{j _ { \mu^\prime}+j_\nu + j^\prime }   \left\ {   \begin{array}{ccc }   j_\mu & j_{\mu^\prime } & j_k \\
j_\nu & j_{\nu^\prime } & j^\prime   \end{array }   \right\ }   ( 2j^\prime + 1 ) \langle [ \mu\nu^\prime]j_k | \bar{v}^{\rm ph }
| [ \mu^\prime \nu]j_k \rangle , \label{f_matrix_element }
\\ h(\mu\mu'\nu\nu';j_k ) & = & \sum_{m_\mu m_{\mu ' } m_\nu m_{\nu ' } } \langle j_\mu m_\mu j_{\mu ' } m_{\mu ' } | jm \rangle \langle j_\nu m_\nu j_{\nu ' } m_{\nu ' } | jm \rangle \bar{v}^{\rm 3p1h}_{\bar{l}\bar{l } ' k \bar{k } ' }   \nonumber\\ & = & \sum_{j^\prime } ( -)^{j _ { \mu}+j_\nu + 1 + l_{\nu^\prime } -j_k - j^\prime }   \left\ {   \begin{array}{ccc }   j_\mu & j_{\mu^\prime } & j_k \\   j_{\nu^\prime } & j_{\nu } & j^\prime   \end{array }   \right\ }   ( 2j^\prime + 1 ) \langle [ \mu\nu]j_k | \bar{v}^{\rm 3p1h } | [ \mu^\prime \nu^\prime]j_k \rangle . \nonumber\\
\label{h_matrix_element}\end{aligned}\ ] ]    we have represented the second derivatives of the energy functional @xmath105 $ ] as unsymmetrized matrix elements of effective interactions @xmath121 , @xmath122 , and @xmath123 .
these effective interactions are given in app .
[ app : matrix_elements ] .
although the `` matrix elements '' are unsymmetrized , the underlying two - quasiparticle states are of course antisymmetric . as a consequence , @xmath124j_k,[\nu\nu^\prime]j_k}= b_{[\mu\mu^\prime]j_k , [ \bar{\nu } \bar{\nu}^\prime ]
j_k}=0 $ ] if @xmath115 is odd and either @xmath125 or @xmath126 .    the nuclear energy functional @xmath105 $ ] is usually separated into particle - hole ( ph ) and pairing pieces ( again , see app .
b for explicit expressions ) .
if the pairing functional , which we will call @xmath127 $ ] , depends on @xmath6 then the derivatives of @xmath127 $ ] with respect to @xmath128 are called pairing - rearrangement terms @xcite . in the qrpa ,
two kinds of pairing - rearrangement terms can arise in general .
one has particle - hole character and is included in @xmath129 ; the other affects 3-particle-1-hole ( 3p1h ) and 1-particle-3-hole ( 1p3h ) configurations and is represented by @xmath130 and @xmath131 .
if the @xmath6-dependence of @xmath127 $ ] is linear , then the ph - type pairing - rearrangement term does not appear .
furthermore , the 3p1h and 1p3h pairing - rearrangement terms arise only for @xmath132 modes if the hfb state has @xmath133 .
most existing work uses a pairing functional that is linear in @xmath6 , and so needs no pairing - rearrangement terms in @xmath134 channels .
in this appendix , we discuss interaction matrix elements coming from @xmath105 $ ] , which we take to contain separate skyrme ( i.e.  strong - force , @xmath106-independent ) , coulomb , and pairing energy functionals : @xmath135 = e_{\rm skyrme}[\rho ] + e_{\rm coul}[\rho_{\rm p } ] + e_{\rm pair}[\rho,\kappa,\kappa^\ast],\ ] ] where @xmath136 is the proton density matrix . the most general skyrme energy functional in common use
is given by @xmath137 & = & \sum_{t=0,1 } \int d^3r      \bigg\ {   c_t^{\rho}[\rho_{00}]\ , \rho^2_{t0 } ( { \bm{r } } )             + c_t^{\delta\rho } \rho_{t0}({\bm{r } } ) \delta \rho_{t0}({\bm{r } } )            + c_t^{\tau } \left [   \rho_{t0}({\bm{r } } ) \
, \tau_{t0}({\bm{r } } )                                -{\bm{j}}^2_{t0}({\bm{r } } )                          \right ]      \nonumber \\ & &         + c_t^{s}[\rho_{00}]\ , { \bm{s}}^2_{t0}({\bm{r } } )             + c_t^{\delta s}{\bm{s}}_{t0}({\bm{r}})\cdot\delta{\bm{s}}_{t0}({\bm{r } } )             + c_t^{t } \left [ { \bm{s}}_{00}({\bm{r } } ) \cdot { \bm{t}}_{t0}({\bm{r } } )                            -\tensor{j}{}^2_{t0}({\bm{r } } )                       \right ]      \nonumber \\ & &         + c_t^{\nabla j }                       \left [ \rho_{t0}({\bm{r}}){\bm{\nabla}}\cdot{\bm{j}}_{t0}({\bm{r } } )                            + { \bm{s}}_{t0 } ( { \bm{r } } ) \cdot { \bm{\nabla } }                             \times { \bm{j}}_{t0}({\bm{r } } )                       \right ]      \bigg\ } . \ ] ] ( see , e.g. , @xcite and references therein for a general discussion . )
all densities are labeled by isospin indices @xmath138 , where @xmath139 takes values zero and one and @xmath140 is always equal to 0 .
a more general theory could violate isospin at the single - quasiparticle level , leading to additional densities @xmath141 @xcite .
we do not consider such densities here .
the @xmath142 are the coupling constants for the effective interaction . as usual , two of them are chosen to be density dependent : @xmath143 & = & a_t^\rho + b_t^\rho \ , \rho_{00}^\alpha({\bm{r } } ) ,        \nonumber \\ c_t^s [ \rho_{00 } ] & = & a_t^s + b_t^s \ , \rho_{00}^\alpha({\bm{r } } ) .\end{aligned}\ ] ] here @xmath144 , @xmath145 , @xmath146 , @xmath147 , @xmath148 , @xmath149 , and @xmath150 are local densities and currents , which are derived from the general density matrices for protons and neutrons @xmath151 where @xmath152 and @xmath153 labels the spin components so that , e.g. , @xmath154 is a spin component of the single - particle wave function associated with the state @xmath95 . defining @xmath155 where @xmath156 is a matrix element of the vector of pauli spin matrices , we write the local densities and currents as @xmath157 the coulomb energy functional is given by @xmath158 =   \frac{e^2}{2 } \iint \ ! d^3r \ , d^3 r '     \frac{\rho_{\rm p}({\bm{r } } ) \rho_{\rm p}({\bm{r}}^\prime ) }          { |{\bm{r}}-{\bm{r}}'| }    -\frac{3}{4 } e^2 \left(\frac{3}{\pi}\right)^{\frac{1}{3 } }     \int \ !
d^3 r \ ; \rho_{\text{p}}^{4/3 } ( { \bm{r } } ) , \ ] ] where we make the usual slater approximation @xcite for the exchange term .    for the pairing functional
we take the quite general form @xmath159 = e_{\text{pair}}[\rho,\tilde{\rho},\tilde{\rho}^\ast ] =   \int \ !
d^3 r \ ;    c^{\tilde\rho}[\rho_{00}({\bm{r } } ) ]    \sum_{\tau={\rm p , n } } |\tilde{\rho}_\tau({\bm{r}})|^2 , \label{pairing_energy_functional}\ ] ] where the density - dependent pairing coupling constant @xmath160 $ ] is an arbitrary function of @xmath161 .
the quantity @xmath162 is defined as @xcite @xmath163 with @xmath164 being the standard pairing tensor in the coordinate representation .    the second derivatives of the energy functional in eq.([ph_matrix_element_uncoupled ] ) , as the equation indicates and we ve already noted , can be written as _ unsymmetrized _
matrix elements @xmath165 of an effective interaction between uncoupled pairs of single - particle states .
the particle - hole matrix elements take the form @xmath166 the last term contains the pairing rearrangement discussed at the end of the previous appendix .
the effective skyrme interaction in eq .
( [ decomp ] ) is given by @xmath167 \
, \rho_{00}^{\alpha-1}({\bm{r } } ) \ , \delta({\bm{r}}-{\bm{r } } ' )        \nonumber\\ &    & + \big [ f_3 \ , \rho^2_{10}({\bm{r } } )               + h_3 \ , { \bm{s}}^2_{00}({\bm{r } } )               + n_3 \ , { \bm{s}}^2_{10}({\bm{r } } )          \big ] \ , \rho_{00}^{\alpha-2}({\bm{r } } ) \ , \delta({\bm{r}}-{\bm{r } } ' )        \nonumber\\ &    & + (    a_4            + c_4 \ , { \vec{\tau } } \cdot { \vec{\tau } } ' ) \ ,                    ( { \bm{\sigma } } + { \bm{\sigma } } ' ) \cdot                    { { \bm{k}}^\dagger}\times \delta({\bm{r}}-{\bm{r } } ' ) \ , { \bm{k } } .
\label{operator_v_skyrme}\end{aligned}\ ] ] @xmath168 is the vector of pauli matrices in isospin space and @xmath169 the coefficients in eq .
( [ skyrme_operator ] ) are defined in tables  [ tab : coefficients_abc ] and  [ tab : coefficients_efg ] .
equation ( [ operator_v_skyrme ] ) contains the usual skyrme - interaction operators , but , the energy functional ( [ e : skyrme ] ) does not necessarily correspond to a real ( density - dependent ) two - body skyrme interaction because the matrix elements are not antisymmetrized . compared to the case usually discussed in the literature , the more general functional relaxes relations that would otherwise restrict the spin - isospin structure of the effective interaction in eq .
( [ skyrme_operator ] ) ; see , e.g. , @xcite for a discussion of the increased freedom .
the densities and currents that appear in eq .
( [ skyrme_operator ] ) come mostly from rearrangement terms and take the values given by the hfb ground state .
the isoscalar and isovector spin densities @xmath170 vanish when the hfb ground state is time - reversal invariant or spherical as assumed here .
the terms containing them will therefore not appear in the expressions for the matrix elements of the effective interaction for such states given below .       the effective coulomb interaction in eq .  ( [ decomp ] ) ,
acting between protons , is given by @xmath175 finally , the ph - type pairing - rearrangement terms in eq .
( [ decomp ] ) come from an effective interaction @xmath176}{d \rho_{00}^2({\bm{r } } ) } \sum_{\tau={\rm p , n}}|\tilde{\rho}_\tau({\bm{r}})|^2\ , \delta({\bm{r}}'-{\bm{r}}).\ ]
]    the second derivatives with respect to @xmath177 also can be written as unsymmetrized matrix elements of effective interactions , this time in the particle - particle channel .
the particle - particle effective interaction entering the matrix elements @xmath178 is obtained from eq .
( [ pairing_energy_functional ] ) through eq .
( [ pp_matrix_element_uncoupled ] ) as @xmath179    \left (   3           - { \bm{\sigma } } \cdot { \bm{\sigma } } '           - { \vec{\tau } }    \cdot { \vec{\tau } } '           - { \bm{\sigma } } \cdot { \bm{\sigma } } ' \ , { \vec{\tau } }    \cdot { \vec{\tau } } '    \right ) \ , \delta({\bm{r}}-{\bm{r } } ' ) .\ ] ] in the numerical calculations of this paper , we use a volume pairing - energy functional , i.e. , @xmath180    last of all are the mixed functional derivatives involving both @xmath6 and @xmath106 ( or @xmath181 ) in eq .
( [ 3p1h_matrix_element_uncoupled ] ) .
they also can be written as the unsymmetrized matrix elements of an effective interaction : @xmath182 where @xmath183 denotes the time - reversed state of @xmath95 , and the 3p1h effective interaction itself is @xmath184}{d\rho_{00}({\bm{r } } ) } \left [ \tilde{\rho}_{\rm p } ( { \bm{r}})\left(1-\tau_z^\prime\right ) + \tilde{\rho}_{\rm n } ( { \bm{r}})\left(1+\tau_z^\prime\right)\right ] \delta({\bm{r}}'-{\bm{r } } ) , \ ] ] where @xmath185 acts on the single - particle states @xmath186 and @xmath183 in eq .
( [ v_3p1h_matrix_element ] ) , and the eigenvalues 1 and @xmath187 are assigned to the neutron and proton , respectively .      to calculate the coupled matrix elements in eqs .
( [ g_matrix_element])([h_matrix_element ] ) , we use an intermediate ls scheme : @xmath188j_k|\bar{v}|[\nu\nu^\prime]j_k\rangle }   \nonumber\\ & = & \sum_{ll^\prime s } \hat{j}_\mu \hat{j}_{\mu^\prime } \hat{j}_\nu \hat{j}_{\nu^\prime }   \hat{l } \hat{l^\prime } \hat{s}^2 \left\ { \begin{array}{ccc } l_\mu & l_{\mu^\prime } & l \\ 1/2 & 1/2 & s \\
j_\mu & j_{\mu^\prime } & j_k \end{array } \right\ } \left\ { \begin{array}{ccc } l_\nu & l_{\nu^\prime } & l^\prime \\ 1/2 & 1/2 & s \\
j_\nu & j_{\nu^\prime } & j_k \end{array } \right\}\nonumber \\ & & \times \langle ( l_\mu l_{\mu^\prime})ls;j_k|\bar{v}|(l_\nu l_{\nu^\prime})l^\prime s;j_k\rangle,\end{aligned}\ ] ] @xmath189 eq .
( [ skyrme_operator ] ) gives    \i ) proton - proton or neutron - neutron matrix elements : @xmath190 ii ) proton - neutron matrix elements : @xmath191 we use the canonical ( and real ) radial wave functions @xmath192 , the angular wave functions @xmath193 , and the spin wave functions to write the nontrivial matrix elements included in eqs.([v_skyrme_proton_proton ] ) and ( [ v_skyrme_proton_neutron ] ) as @xmath194^\ast_{l0 } \left [ y_{l_\nu}(\omega ) y_{l_{\nu^\prime}}(\omega)\right]_{l0 } , \nonumber\\ \label{matrix_element_delta}\end{aligned}\ ] ] @xmath195^\ast_{l0 }   \left [ y_{l_\nu}(\omega ) y_{l_{\nu^\prime}}(\omega)\right]_{l0 } = \frac{1}{4\pi}\hat{l_\mu } \hat{l_{\mu^\prime } } \hat{l_\nu } \hat{l_{\nu^\prime } } \left(\begin{array}{ccc}l_\mu & l_{\mu^\prime } & l \\ 0 & 0 & 0\end{array}\right ) \left(\begin{array}{ccc}l_\nu & l_{\nu^\prime } & l \\ 0 & 0 & 0\end{array}\right ) , \ ] ] @xmath196^\ast_{l0 }   \left [ y_{l_\nu}(\omega ) y_{l_{\nu^\prime}}(\omega)\right]_{l0 }   \nonumber\\ & & \times \int dr\:r^2 r_\mu(r ) r_{\mu^\prime}(r ) \left\ { \left [ \left ( \frac{d^2}{dr^2}+\frac{2}{r}\frac{d}{dr}-\frac{l_\nu(l_\nu+1)}{r^2 } \right)r_\nu(r)\right]r_{\nu^\prime}(r)\right.\nonumber\\ & & \left.+r_\nu(r)\left [ \left ( \frac{d^2}{dr^2}+\frac{2}{r}\frac{d}{dr }   -\frac{l_{\nu^\prime}(l_{\nu^\prime}+1)}{r^2}\right)r_{\nu^\prime}(r )   \right]\right\}\nonumber\\ & & -\sum_{\delta l_\nu=0,1 } \sum_{\delta l_{\nu^\prime}=0,1 }   \frac{1}{2}\sqrt{(l_\nu+\delta l_\nu)(l_{\nu^\prime}+\delta l_{\nu^\prime } ) }   \left\{\begin{array}{ccc}l_{\nu^\prime}-1 + 2\delta l_{\nu^\prime } &   l_{\nu}-1 + 2\delta l_{\nu } & l \\
l_\nu   & l_{\nu^\prime } & 1 \end{array}\right\ } \nonumber\\ & & \times\int dr\:r^2 r_\mu(r ) r_{\mu^\prime}(r ) \left\ { ( l_\nu+1-\delta l_\nu)\frac{r_\nu(r)}{r }   + ( -)^{\delta l_\nu}\frac{dr_\nu(r)}{dr}\right\}\nonumber\\ & & \times \left\ { ( l_{\nu^\prime}+1-\delta l_{\nu^\prime})\frac{r_{\nu^\prime}(r)}{r }   + ( -)^{\delta l_{\nu^\prime}}\frac{dr_{\nu^\prime}(r)}{dr}\right\}\nonumber\\ & & \times   \delta_{ll^\prime }   \int d\omega \left [ y_{l_\mu}(\omega ) y_{l_{\mu^\prime}}(\omega)\right]^\ast_{l0 }   \left [ y _ { l_{\nu^\prime}-1 + 2\delta l_{\nu^\prime } } ( \omega )   y_{l_\nu-1 + 2\delta l_\nu } ( \omega)\right]_{l0 } , \end{aligned}\ ] ] @xmath197 @xmath198 the square brackets around products of spherical harmonics and the parentheses surrounding products of operators indicate angular - momentum coupling .    to evaluate eqs .
( [ t2 ] ) and ( [ spin_orbit ] ) , one can use @xmath199   -(-)^{l_\nu+l_{\nu^\prime}+l^\prime}[\nu \leftrightarrow \nu^\prime ]   -(-)^{l_\mu+l_{\mu^\prime}+l}[\mu \leftrightarrow \mu^\prime ] \nonumber\\ & &   + [ \mu \leftrightarrow \mu^\prime { \rm \ and\ } \nu \leftrightarrow \nu^\prime ] , \end{aligned}\ ] ] where for reduced matrix elements we have used the convention @xmath200 and made the abbreviation @xmath201
-(-)^{l_\nu+l_{\nu^\prime}+l^\prime}[\nu \leftrightarrow \nu^\prime ]   -(-)^{l_\mu+l_{\mu^\prime}+l}[\mu \leftrightarrow \mu^\prime ]   + [ \mu \leftrightarrow \mu^\prime { \rm \ and\ }      \nu \leftrightarrow \nu^\prime ] \nonumber\\ & \equiv & a_{\mu\mu^\prime\nu\nu^\prime } - ( -)^{l_\nu+l_{\nu^\prime}+l^\prime } a_{\mu\mu^\prime\nu^\prime\nu } - ( -)^{l_\mu+l_{\mu^\prime}+l}a_{\mu^\prime\mu\nu\nu^\prime } + a_{\mu^\prime\mu\nu^\prime\nu } .\end{aligned}\ ] ]    eq .
( [ matrix_element_delta ] ) , modified to include additional factors in the radial integral , can also be used ( together with the subsequent equations ) to evaluate the matrix elements of the terms involving @xmath202 in @xmath203 , the coulomb - exchange interaction , and the contributions of the pairing functional to the effective ph , pp , and 3p1h interactions .
the coulomb - direct term can be evaluated in a similar but slightly more complicated way , via a multipole expansion .    in the main part of this paper we used the skyrme functional skm@xmath28 , which is usually parameterized as in interaction in terms of coefficients @xmath204 and @xmath205 .
the relations between these coefficients and those used here , if no terms are neglected , are @xcite @xmath206 in the hf fits that originally determined the skm@xmath28 parameters , the effects of @xmath207 ( the `` @xmath208 terms '' ) were neglected because of technical difficulties .
these terms have often been included in subsequent rpa calculations . to maintain self consistency here ,
we have set them to zero both in the hfb calculation and in the qrpa .