one of the spectacular predictions of quantum electrodynamics is the emergence of forces that arise purely from quantum fluctuations of the electromagnetic vacuum  @xcite .
also known as london - van der waals  @xcite or casimir forces  @xcite in different regimes , these forces are often dominant at short distances and can give rise to undesirable effects in nanoscale systems , such as nanomechanical stiction  @xcite .
these forces have also attracted increasing attention in the fields of atomic physics and quantum optics . in particular ,
significant efforts have been made in recent years to interface cold atoms with the evanescent fields of dielectric micro- and nano - photonic systems  @xcite .
these systems are expected to facilitate strong , tunable interactions between individual atoms and photons for applications such as quantum information processing  @xcite and the investigation of quantum many - body physics  @xcite . in practice
, efficient atomic coupling to the evanescent fields of the nanophotonic systems requires that atoms be trapped within sub - wavelength distances of these structures .    at these scales
, quantum vacuum forces can overwhelm the forces associated with conventional optical dipole traps , typically resulting in a loss of trap stability at distances @xmath0 nm from dielectric surfaces .
given the ability to engineer the properties of nanophotonic structures , an interesting question arises as to whether such systems could be used to significantly modify vacuum forces , perhaps changing their sign from being attractive to repulsive , or even creating local potential minima .
the strength of nanoscale vacuum forces should lead to unprecedented energy and length scales for atomic traps , which would find use beyond nanophotonic interfaces , such as in quantum simulation protocols based upon ultracold atoms  @xcite and control of inter - atomic interactions  @xcite .    here
we propose a novel mechanism in which engineered vacuum forces can enable the formation of a nanoscale atomic trap . while a recent no - go theorem  @xcite forbids a vacuum trap for atoms in their electronic ground states , tailored nanophotonic systems can yield strong repulsive potentials for atomic excited states .
we show that a weak external optical field can give rise to an overall trapping potential for a dressed state .
remarkably , absent fundamental limits on the losses of the surrounding dielectric structure , the fraction of excited - state population in the dressed state can become infinitesimal , which greatly enhances the trapping lifetime and stability .
we identify and carefully analyze the actual limiting mechanisms . while we present calculations on a simple model where analytical results can be obtained , we also discuss the generality of our protocol to realistic systems such as photonic crystal structures .    the no - go theorem for non - magnetic media states that a dielectric object in vacuum can not be stably trapped with vacuum forces for any surrounding configuration of dielectric objects , provided that the system is in thermal equilibrium  @xcite . in analogy with earnshaw s theorem , which prohibits trapping of charged objects with static electric potentials , at best
one can create a saddle - point potential  ( see , _
e.g. _ , ref .
@xcite involving metal particles and ref .
@xcite for a hybrid vacuum / optical trap for atoms ) . for dielectric objects ,
potential loopholes involve embedding the system in a high - index fluid  @xcite , or using blackbody radiation pressure associated with strong temperature gradients  @xcite .
applied to atoms , this theorem essentially forbids stable vacuum trapping near a dielectric structure when the atom is in its electronic ground state .
no such constraint exists for atoms in their excited states  @xcite , although the robustness of any possible trap would be limited by the excited state lifetime . here
, we show that an atom only weakly dressed by its excited state can be trapped by properly tailoring the dispersion of the underlying structure .
we begin by briefly reviewing how quantum vacuum fluctuations give rise to forces on the ground and excited states . within an effective two - level approximation of an isotropic atom with ground and excited states @xmath1  ( see appendix for description of how isotropy is incorporated into the model ) , the interaction between the electromagnetic field and an atom at position @xmath2 is given by the dipole hamiltonian @xmath3 . following a complete decomposition of the electric field into its normal modes @xmath4 , one can write @xmath5 , where @xmath6 . the energy non - conserving terms  ( @xmath7 and @xmath8 ) enable an atom in its ground state @xmath9 to couple virtually to the excited state and create a photon , @xmath10 , which is subsequently re - absorbed .
the corresponding frequency shift for the ground state within second - order perturbation theory is @xmath11 , where @xmath12 is the unperturbed atomic transition frequency .
this shift results in a mechanical `` vacuum '' potential when translational symmetry is broken due to dielectric surfaces . using a quantization technique for electrodynamics in the presence of dispersive dielectric media ,
the shift can be expressed in terms of the scattered component of the dyadic electromagnetic green s function evaluated at imaginary frequencies @xmath13  @xcite  ( see appendix ) , _
g()=_0^dug_(,,iu ) , where @xmath14 is the free - space spontaneous emission rate of the atom .    a similar process can occur for an atom in its excited state @xmath15 , in which the atom virtually emits and re - absorbs an off - resonant photon  ( @xmath16 ) .
however , the excited state is unique in that it can also emit a resonant photon .
the total excited state shift is given by  @xcite _ e()=--g_(,,_0).[eq : omegae ] the first term on the right arises from off - resonant processes , while the second term can be interpreted as the interaction between the atom and its own resonantly emitted photon .
naturally , the dielectric environment can modify the spontaneous emission rate as well , ( ) = _ 0+g_(,,_0).[eq : gamma ] note that the modified emission and resonant shifts are complementary in that they emerge from different quadratures of the green s function .
the versatility of micro- and nanophotonic systems has already been exploited to greatly modify emission rates of atoms and molecules in a variety of contexts  @xcite .
we propose that these same systems can be used to tailor excited - state potentials in order to facilitate vacuum trapping .
of course , the excited - state shift can be affected by coupling to even higher electronic levels . however , since our goal is to engineer giant resonant shifts of @xmath17 , the off - resonant couplings to higher - lying states yield only minor corrections and consequently are ignored here .
in what follows , we present a simple model system that enables one - dimensional trapping in a plane parallel to a semi - infinite dielectric slab , as illustrated in fig .
[ fig : schematic ] . the use of a homogeneous dielectric  ( as opposed to , _
e.g. _ , a photonic crystal ) greatly simplifies the calculation and enables one to understand the relevant trap properties analytically , although we later argue why the qualitative features should still hold in more complex geometries . in the dielectric slab model , dispersion engineering will be realized via the frequency - dependent electric permittivity @xmath18 .
while this particular system achieves a one - dimensional trap along @xmath19 , our electrodynamic calculations are performed in three dimensions .
enables large repulsion in the excited state @xmath20 , in contrast to the attractive vacuum potential experienced by an atom in its ground state .
an overall trapping potential normal to the surface ( along @xmath19 ) can be generated by driving the atom with an external laser of frequency @xmath21 , which is greater than the natural resonance frequency @xmath12 of the atom .
the laser comes into local resonance with the atom at @xmath22 .
an atom far from @xmath22 is unaffected by the largely detuned field and feels an attractive force @xmath23 toward the dielectric interface .
b ) the atom comes closer to local resonance with the laser as @xmath24 .
this creates an atomic dressed state , whose excited state component yields a net repulsive force on the atom .
c ) a slowly moving atom experiences the adiabatic dressed state potential @xmath25 shown in red , which characterizes the position - dependent mixing between ground and excited states .
the minimum of the dressed potential is located at @xmath26.[fig : schematic],width=491 ]
a short - distance expansion of the green s function , valid for sub - wavelength scales , shows that the excited - state emission rate @xmath27 and resonant contribution @xmath28 to the shift are given by & & ( ) , [ eq : gammageneral ] + & & -(),[eq : omegaegeneral ] where the permittivity @xmath29 is evaluated at the atomic resonance frequency , @xmath19 is the distance from the surface , and @xmath30 is the resonant free - space wavevector .
a large repulsive potential is generated when @xmath31 , which is analogous to the interaction between a classical dipole at position @xmath19 and its large induced image dipole in the dielectric .
the modified spontaneous emission rate originates from the absorption or quenching of the atomic emission due to material losses .
we choose a drude model , @xmath32 , as a simple dielectric function .
note that this function satisfies kramers - kronig relations  ( causality ) , as a choice of a non - causal function could result in an apparent violation of the no - go theorem  @xcite as well . in the limit of vanishing material loss parameter @xmath33
, the system passes through @xmath34 at the plasmon resonance frequency @xmath35
. the system is conveniently parameterized by the quality factor @xmath36 and a dimensionless detuning @xmath37 . then , for large @xmath38 and @xmath39 , @xmath40 and @xmath41 .
the dispersion and dissipation scale like @xmath42 and @xmath43 , respectively .
thus , for high @xmath38 , it is possible to choose detunings where the atom still sees significant repulsive forces , but where spontaneous emission into the material is not significantly enhanced  ( in addition to material - induced emission of eq .
( [ eq : gammageneral ] ) , there still is emission into free - space at a rate @xmath44 ) .
this scaling behavior is in close analogy with conventional optical trapping of atoms , where simultaneously using large trapping intensity and detuning maintains reasonable trap depths but suppresses unwanted photon scattering .
the ground - state shift for atomic frequencies @xmath45 and in the near - field is given by @xmath46 .
we now describe how a trap can form for an atom subject to these potentials , and in the presence of a weak driving laser of frequency @xmath21 and rabi frequency @xmath47  ( see fig .
[ fig : schematic ] ) . for conceptual simplicity ,
we choose the rabi frequency to be spatially uniform , such that all atomic forces are attributable to the vacuum potentials alone .
qualitatively , the strong position dependence of the ground and excited state frequencies causes the laser to come into resonance with the atom at a single , tunable point @xmath48 . under certain conditions  ( described in detail later )
, an atom starting at position @xmath49 will be far detuned , such that it is essentially in the ground state and is attracted by the pure ground - state potential .
however , moving closer to @xmath22 brings the atom closer to resonance .
the dressing of the atom with a small fraction of excited state population causes a repulsive barrier to form near @xmath50 , yielding a metastable trap .
significantly , our approach does not attempt to directly counteract the attractive ground - state potential with large external optical potentials , in sharp contrast with other trapping schemes near dielectric surfaces  @xcite .
the order of our calculation of the trap properties is as follows .
we first find the adiabatic potential experienced by a slowly - moving atom .
we derive relevant properties around the location of the trap minimum @xmath51  ( fig .
[ fig : schematic ] ) .
in particular , we show that absent material constraints  ( arbitrarily high @xmath38 ) , the dressed state can in principle have infinitesimal excited state population and scattering .
we then quantize the motion to find the motional eigenstates , binding energy , and position uncertainty @xmath52 .
finally , we analyze the mechanisms that limit the trap lifetime .    first treating the motion classically
, the force experienced by the atom is given by @xmath53 , while the internal atomic dynamics satisfies the usual bloch equations , _
e.g. _ , = ( i_a(z)-(z)/2)_ge+i(_ee-_gg)[eq : bloch ] . here
@xmath54 is the detuning between the laser and the local atomic resonance frequency .
when the atom moves slowly on the time scales of the internal dynamics , the atomic coherence can be adiabatically eliminated , @xmath55 , such that the atomic populations are functions of position alone , @xmath56 .
we are primarily interested in the regime where the atom is weakly driven , @xmath57 .
a more careful calculation reveals that the green s function @xmath58 for the excited state resonant shift and emission rate in eqs .
( [ eq : omegae ] ) and  ( [ eq : gamma ] ) must in fact be evaluated at the laser frequency in this regime , which reflects that the atom primarily acts as a rayleigh scatterer  ( see appendix ) .
consequently , @xmath59 in eqs .
( [ eq : gammageneral ] ) and  ( [ eq : omegaegeneral ] ) is replaced with @xmath60 .
the adiabatic potential seen by a slowly - moving atom is @xmath61 .
three independent frequencies characterize the system  ( @xmath62 , and @xmath63 ) , and two relative frequencies must be specified .
one is determined by selecting the trap barrier position , determined by the relation @xmath64 .
a second is the dimensionless detuning between the plasmon resonance and laser , @xmath65 , which we treat as a free parameter that determines the relative strengths of the excited state shift and dissipation . in the near - field ,
the excited state repulsion exceeds the ground state attraction by a factor @xmath66 .
the population at @xmath22 must then exceed the inverse of this quantity , @xmath67 , in order to provide a classical barrier in the adiabatic potential , which translates to a minimum rabi frequency @xmath68 .
since @xmath69 , this scheme potentially enables atom trapping near surfaces with greatly reduced intensities as compared to conventional trapping . in the latter case
, an optical dipole potential requires a minimum driving field of @xmath70 to overcome the attractive ground state potential .
the excited - state population at the trap minimum is fixed by the ratio of ground to excited state forces , @xmath71 , and is surprisingly independent of @xmath47  ( provided that @xmath72 ) .
an increasing rabi frequency has the effect of increasing the local detuning @xmath73 to preserve the same population .
this manifests itself as a `` self - selection '' of the trapping position and an increasing distance from the barrier , ( ^2 - 1)^1/2 , where we have defined the dimensionless rabi frequency @xmath74 .
the photon scattering rate at the trap position , @xmath75 , leads to lifetime limitations of our trapping scheme .
fixing other parameters , the detuning @xmath76 from the plasmon resonance can be optimized to yield the minimum scattering rate , @xmath77 .
thus , absent fundamental limitations to the magnitude of @xmath38 , the photon scattering rate can in principle be suppressed to an arbitrary degree .
this implies that an `` infinitesimal '' violation of the assumptions underlying the no - go theorem can in fact allow for trapping .    in a conventional optical dipole trap ,
the potential @xmath78 seen by an atom depends on the local field intensity and is proportional to a spatially - constant atomic polarizability evaluated at the laser frequency .
in contrast , here we exploit a novel mechanism arising from strong spatial variations in the polarizability itself , induced through the vacuum level shifts .
this creates a back - action effect analogous to that producing a dynamical `` optical spring '' in opto - mechanical systems  @xcite , and enables extremely steep trap barriers around @xmath22  ( as in fig .
[ fig : schematic]c ) . in particular , the repulsive force is given by @xmath79 , where @xmath80 is engineered to be large .
however , the excited state population @xmath81 itself varies strongly with position through the spatially - dependent detuning  ( for simplicity , here we omit the varying linewidth @xmath27 ) .
this provides an overall repulsive force that scales like @xmath82 , f~-()^2 .     and an infinite barrier at @xmath48 .
the maximum depth of the potential is denoted by @xmath83 .
the quantum wave - function amplitudes @xmath84  ( in arbitrary units ) for the two lowest eigenstates are illustrated in red , while the binding energy @xmath85 and position uncertainty @xmath52 of the ground state are depicted as well .
b ) properties of the ideal trapping potential versus barrier position @xmath22 , for the case of a cesium atom .
the blue solid and dashed curves give the potential depth at the barrier position , @xmath86 , and the ground - state binding energy @xmath85 , respectively .
the green curve depicts the ground - state uncertainty @xmath52 .
c ) for the non - ideal trapping potential , the depth of the classical potential @xmath87 is numerically calculated . here
we plot it normalized by the maximum possible value @xmath83 , versus detuning parameter @xmath76 and normalized rabi freuqency @xmath88 .
d ) the trapping potential , ground - state binding energy @xmath85 , and wave - function amplitude @xmath84 for the parameter sets a , b , as illustrated in c ) .
the color coding is the same as in a).[fig : properties],width=566 ]    given the possible steepness of this barrier , an `` ideal '' limit for the overall potential is shown in fig .
[ fig : properties]a .
the attractive part consists of the pure ground - state potential and arises from quantum fluctuations , while an infinite `` hard wall '' is created at @xmath48 due to the combination of strong resonant excited - state shifts  ( originating from the interaction of the atom with a `` classical '' image dipole ) and opto - mechanical back - action .
this idealized model enables one to understand the best possible scaling of the trap properties .
in particular , the trap depth @xmath87  ( the energy required for a classical particle to become unbound ) can not exceed the value of the ground - state potential at @xmath22 , @xmath89  ( _ e.g. _ , @xmath90  mk when @xmath91  nm for a cesium atom ) .
this trapping potential can be easily quantized numerically . in fig .
[ fig : properties]b , we plot the quantum binding energy @xmath85  ( where @xmath92 ) and position uncertainty @xmath52 for the motional ground - state wave function . for the numerical results in figs .
[ fig : properties ] and  [ fig : lifetime ] , we use atomic properties corresponding to cesium  ( @xmath93  nm , @xmath94  mhz , recoil frequency @xmath95 ) . the strength of the vacuum potentials is reflected in the strong confinement of the wave - functions  ( _ e.g. _ , @xmath96  nm for distances @xmath97  nm ) , which are an order of magnitude smaller than in conventional optical traps .
to contrast our results for cesium with those of a less massive atom , we have considered lithium as well .
results for li similar to those for cs in fig .
[ fig : properties]a , b show that at @xmath91  nm , the potential depth is smaller by roughly a factor of two , while the confinement @xmath52 sees roughly a five - fold increase .
generally , within the wkb approximation and to lowest order in @xmath52 , the ground - state localization scales like @xmath98 , while the binding energy @xmath99 can be nearly the entire classical potential depth .
this ideal model describes qualitatively well the actual potentials over a large parameter regime . in fig .
[ fig : properties]c , we plot the ratio of the numerically obtained trap depth to the theoretical maximum , @xmath100 , as functions of @xmath76 and dimensionless driving amplitude @xmath101 . here , we have chosen a quality factor of @xmath102 and a barrier position of @xmath91  nm . although we have provided near - field approximations for the energy shifts above , the plots in figs .  [
fig : properties]c , d and fig .
[ fig : lifetime ] are calculated using the full green s functions .
a trap depth comparable to @xmath86 is achievable over a large range of values . for low values of @xmath101 , the trap depth
is weakened due to escape over the potential barrier , while for large values , the trap position becomes shallow as @xmath51 is pulled away from @xmath22  ( as discussed previously ) .
two representative trap potentials , along with the ground - state wave functions and confinement @xmath52 are illustrated in fig .
[ fig : properties]d .
we have identified four sources of motional heating  recoil heating , non - adiabatic motion of the atom within the trap , fluctuations of the atomic coherence , and tunneling over the finite - height potential barrier  that limit the trap lifetime in absence of external cooling mechanisms .
tunneling is exponentially suppressed with barrier height and thus is easily suppressed with increasing rabi frequency .
this calculation can be found in the appendix , while we discuss the more important mechanisms here .    as in conventional optical trapping ,
the random momentum recoil imparted on the atom by scattered photons yields an increase in motional energy at the rate @xmath103 .
two qualitative differences emerge relative to free space , however .
first , the photon scattering rate @xmath104 is modified in the presence of dielectric surfaces .
second , the effective momentum @xmath105 imparted by a photon is enhanced compared to the free - space momentum .
this momentum scales at close distances like @xmath106  ( see appendix ) , due to emission into high - wavevector photons in the dielectric . starting from the motional ground state ,
this heating causes the atom to become unbound over a time scale @xmath107 .
it should be noted that the heating rate @xmath108 does not depend on the atom confinement  ( _ i.e. _ , whether we operate in the lamb - dicke limit ) , as the suppression of motional jumps in tight traps is compensated by the increase in energy gained per jump  @xcite .
our derivation of the atomic potential assumes that the atomic coherence adjusts rapidly to the local trap properties , @xmath109 , which results in a purely conservative potential . as we now show , non - adiabatic corrections result in a momentum - dependent force , @xmath110  ( anti - damping ) .
we can solve eq .
( [ eq : bloch ] ) for the atomic motion perturbatively by writing @xmath111 , where @xmath112 is the adiabatic solution .
it depends only on the instantaneous position and is given by @xmath113 in the weak saturation limit . here
, @xmath114 is the complex detuning . substituting this solution into eq .
( [ eq : bloch ] ) yields a velocity - dependent correction to the coherence , @xmath115 .
this in turn yields a new contribution to the force , @xmath116 , where the anti - damping rate is given by @xmath117 $ ] evaluated at the trap position @xmath51 . intuitively , the atom is anti - damped  ( @xmath118 ) because the laser frequency is blue - detuned relative to the atomic frequency at @xmath51 , giving rise to preferential stokes scattering over anti - stokes . aside from the term proportional to @xmath119 ,
the heating rate is exactly analogous to that occurring in an opto - mechanical system in the blue - detuned regime  @xcite .
we define a characteristic time @xmath120 needed for the energy increase produced by anti - damping to equal the ground - state binding energy @xmath85 .
thus far , we have focused on the motion of the atomic dressed state , in which a trap emerges due to weak mixing with the excited state potential . following a spontaneous emission event
, however , the atom returns to the ground state and sees the pure ground - state potential for a transient time @xmath121 until the atomic coherence is restored , during which the atomic wave - packet accelerates toward the surface and gains energy . while a simulation of the full dynamics of the atom  ( including the spatial and internal degrees of freedom ) is quite challenging , here we estimate the heating rate from these transient processes . in particular , we first find the spatial propagator under the pure ground - state potential , @xmath122  ( see appendix ) .
then , assuming that a spontaneous emission event occurs at @xmath123 with the atom initially in the motional ground - state @xmath124 , the new wave - function that emerges once the coherence equilibrates is given by @xmath125 .
the probability of atom loss @xmath126 per spontaneous emission event is estimated from the overlap between the initial and final states , @xmath127 , with a corresponding trap lifetime of @xmath128 .    in fig .
[ fig : lifetime ] , we plot the overall trap lifetime , which is obtained by adding the individual lifetimes in parallel , @xmath129 .
the overall trap lifetime exhibits a complicated dependence on rabi frequency @xmath47 and detuning @xmath76 due to the different scalings of the constituent heating mechanisms . as an example , however , for external parameters @xmath130 and @xmath131 for atomic cs , one achieves a lifetime of @xmath132  ms , ground - state confinement @xmath133  nm , trap depth @xmath134  mk , and trap distance @xmath135  nm .     versus dimensionless rabi frequency @xmath74 and detuning @xmath76 .
here we have assumed a plamson resonance quality factor of @xmath102 and atomic properties corresponding to cesium .
[ fig : lifetime],width=302 ]
we have described a protocol for an atomic trap based upon engineered vacuum forces , and have analyzed in detail a model case of one - dimensional trapping near a drude material , where analytical results are possible to obtain .
the drude response approximates well a number of metals , such as silver and gold , in the optical domain  @xcite . in practice , however , these metals are limited by their low quality factors  ( @xmath136 ) , and the thermal fluctuations of such conducting materials can give rise to strong noise - induced trap heating and decoherence  @xcite . on the other hand , photonic crystal structures enable engineering of resonances through geometry , and quality factors approaching @xmath137 have been observed  @xcite . while the field profiles of such structures are more complex , we anticipate that our trapping protocol could be quite generally applied , and the resulting trap behavior would qualitatively remain the same . in particular , the scaling of the engineered excited - state shift  ( @xmath138 ) and dissipation  ( @xmath139 ) are generic to coupling between an atom and resonator , and only the spatial dependence is determined by the field profiles .
photonic crystals should thus offer great flexibility in tailoring the properties of vacuum traps  ( such as dimensionality ) , which we plan to investigate in detail in future work .
we note that related efforts have already been made to design realistic hybrid atom traps in photonic crystal structures , where a combination of vacuum forces and optical forces are required to stabilize the trap in all dimensions  @xcite .
in addition , while we have focused on trapping of atoms weakly dressed by an external laser , it should also be possible to achieve laser - free traps by using atoms pumped into metastable electronic levels and engineering resonant shifts .    the ability to create atomic traps with parameter sets  ( such as depth , confinement , and proximity to surfaces ) that are not possible in conventional traps  ( whether optical , electrical , or magnetic ) should lead to a number of intriguing applications .
for example , it has been proposed that the trapping of atoms near dielectric surfaces with nano - scale features  ( such as a sub - wavelength lattice ) would facilitate large interaction strengths and long - range interactions for ultracold atoms  @xcite . however , the major limitation of optical trapping is the divergent intensity needed to overcome the ground - state vacuum forces of these structures .
the use of resonantly - enhanced excited state repulsion to create a trap could significantly reduce the power requirements and enable smaller lattice constants .
furthermore , traps with localization on the nanometer scale can induce inter - atomic forces that are of comparable strength to molecular van der waals forces , which enables novel opportunities to control ultracold atomic collisions  @xcite and realize exotic interactions such as @xmath140-wave scattering of fermions  @xcite .
atoms trapped near surfaces could also act as exquisite nanoscale probes of surface physics , including the precision measurement of vacuum forces themselves .
the authors thank n. stern and o. painter for helpful discussions .
dec acknowledges support from fundaci privada cellex barcelona .
ks was funded by the nsf physics frontier center at the jqi .
jmt acknowledges funding from the nsf physics frontier center at the jqi and the us army research office muri award w911nf0910406 .
hjk acknowledges funding from the iqim , an nsf physics frontier center with support of the moore foundation , by the afosr qumpass muri , by the dod nsseff program , and by nsf phy-1205729 .
in this section we briefly describe the quantization scheme of the electromagnetic field in the presence of arbitrary linear dielectric media , which follows that of ref .
the free field is characterized by a set of bosonic modes with frequencies @xmath141 and annihilation operators @xmath142  ( @xmath143 ) , with a corresponding hamiltonian h_f = _ jd _ j^(,)_j ( , ) .
the bosonic operators satisfy canonical commutation relations @xmath144=\delta_{jk}\delta(\bfr-\bfr')\delta(\omega-\omega')$ ] and are physically associated with noise polarization sources in the material .
the electric field operator at frequency @xmath141 is driven by the sources at the same frequency and is given in the coulomb gauge by  @xcite e_i(,)= i_0 ^2 d g_ij(,,)_j(,).[eq : e ] here and in the following , it is assumed that all repeated vector or tensor indices are summed over , _
e.g. _ , @xmath145 .
@xmath146 is the frequency - dependent permittivity at position @xmath2 , and @xmath147 is the dyadic green s function which is the ( gauge - invariant ) solution to g(,,)=(-)i.[eq : greeneq ] the total electric field is @xmath148 similarly , the magnetic field is @xmath149 . it can be verified that these expressions preserve the same field commutation relations as in free - space field quantization .    in the following ,
we consider a two - level atom at position @xmath150 and with ground and excited states @xmath151,@xmath17 interacting with the electromagnetic field within the electric dipole approximation .
the interaction hamiltonian is given by h _
& = & -*d*(_a ) + & = & -_i(_ege_i(_a)+h.c.),[eq : haf ] where @xmath152 is the atomic dipole matrix element .
our system of interest consists of an atom in a region of vacuum with nearby dielectric objects . in general , for @xmath153 in the vacuum region , one can express the green s function as a sum @xmath154 , where the free component is the solution in uniform vacuum  ( @xmath155 ) and the scattered component physically describes waves scattered from the dielectric objects .    in a planar geometry such as that considered in the main text ,
the green s function can be evaluated exactly using a plane wave expansion .
specifically , suppose that we have an interface with vacuum in the region @xmath156 and a dielectric with permittivity @xmath18 in the region @xmath157 .
the free and scattered components of the dyadic green s function in the region @xmath158 are then given by g _ & = & -d _ e^i_(-),[eq : gfree ] + g _ & = & d_e^i_(- ) e^i(z+z ) .[eq : gsc ] here @xmath159 and @xmath160 are the wavevector and position along the direction parallel to the interface , @xmath161 , @xmath162 , and @xmath163 and @xmath164 are the fresnel reflection coefficients for @xmath165 and @xmath140 polarized waves , respectively .
the interaction hamiltonian @xmath166 couples the ground state of the bare system  @xmath9  ( consisting of the atomic ground state and vacuum ) to states @xmath167 . within second - order perturbation theory ,
this interaction induces a shift of the bare ground state energy by an amount _
g(_a)=-_jdd , where @xmath12 is the bare resonance frequency of the atom . substituting eqs .
( [ eq : e ] ) and  ( [ eq : haf ] ) into the expression above , and using analyticity of the integrand to rotate the frequency integral onto the positive imaginary axis @xmath13 , one finds  @xcite _ g(_a ) = _ 0^ du _
i g_,ij(_a,_a , iu ) _ j.[eq : omegag ] here we have only included the scattered component of the green s function , as the free component gives a contribution that is independent of position . while we have considered a two - level atom for simplicity , a more realistic model of an isotropic atom consists of multiple excited states and equal polarizabilities in each direction  ( _ i.e. _ , the ground state can emit virtual photons of any polarization with equal strength ) . since the contribution to the ground - state shift from each excited state transition is additive , eq .
( [ eq : omegag ] ) can be modified to account for this by setting @xmath168 for all @xmath169 and associating @xmath170 with the free - space spontaneous emission rate of any excited state .
this yields eq .
( 1 ) in the main text .
in this section , because we want to focus exclusively on the effect of vacuum interactions with the excited state , it is sufficient to include only energy - conserving terms in the atom - field interaction hamiltonian @xmath166  ( _ i.e. _ , allowing a transition from the excited to ground state accompanied by creation of a photon from vacuum ) . we also want to consider the specific case where the atom is weakly driven by a classical field of frequency @xmath21 , which motivates working in a rotating frame where the free atomic evolution is given by @xmath171 and @xmath172 is the detuning between the laser frequency and atomic resonance frequency @xmath12 . as will be seen , the inclusion of the laser frequency in the free hamiltonian causes the effective excited state shift and decay rates to depend on the green s function evaluated at @xmath21 once the vacuum modes are eliminated .
physically , this describes the effect of the weakly driven atom  ( which primarily acts as a rayleigh scatterer ) interacting with its own rayleigh - scattered field , which can be altered due to the presence of nearby dielectric surfaces .
the standard technique to derive the excited - state shift and decay rate follows that of appendix  [ sec : ground ] for the ground state and uses time - independent second - order perturbation theory  @xcite . in this approach
, however , one can not find the jump operator associated with photon emission , which is needed in order to calculate the effect of emission on atomic motion  ( _ e.g. _ , recoil heating ) . to rectify this
, we employ an alternative approach based on deriving an atomic master equation within the born - markov approximation  @xcite . in particular , the equation of motion for the reduced atomic density matrix @xmath173
is given by _
a = _ f_0^,[eq : rhodot ] where @xmath174 is the interaction hamiltonian in the interaction picture  ( with respect to the free hamiltonian @xmath175 ) and @xmath176 is the electromagnetic vacuum state . any exponential of time appearing in the integral can be evaluated using the relation @xmath177 . using the field expansion of eq .
( [ eq : e ] ) in @xmath174 and following some manipulation , eq .
( [ eq : rhodot ] ) can be written in the form _
a=-(i/)[h_,_a]+[_a ] .
the effective hamiltonian describes coherent interactions between the atom and vacuum field and takes the form h_=_e(_a)_ee , which can be interpreted as a position - dependent energy shift of the excited state .
a realistic model of an atom allows for the possibility for an excited state to emit photons of different polarizations and transition into different states in a ground - state manifold .
the different clebsch - gordan coefficients of these transitions , along with the symmetry breaking introduced by the nearby dielectric structure , imply that different excited states can experience different shifts absent some special engineering .
to arrive at a simpler two - level atomic model that captures the essential physics of our trapping mechanism , we assume that the excited state is equally coupled to all polarizations  ( but still with total emission rate @xmath14 ) , such that one should average over all dipole orientations . in this case , the excited - state shift is given by _
e(_a)=-_0^dug_(_a,_a , iu)-g_(_a,_a,_l),[eq : omegaeappendix ] as stated in the main text .
the liouvillian term @xmath178 $ ] can be written in a lindblad form [ _ a]=_jd _ j_a _ j^-(_j^_j_a+_a_j^_j ) .
the jump operators take the form _
i g^_ij(_a,,_l)_ge.[eq : ojump ] the excited - state spontaneous emission rate is in turn given through the relation @xmath179 , or ( _ a ) = _ 0+g_(_a,_a,_l).[eq : gammaappendix ] eq .
( [ eq : gammaappendix ] ) again averages over dipole orientations and reproduces the expression for the emission rate provided in the main text . to derive this
, we have used the identity d  ( (  , ) ) g_ij(,,)g_jk^(,,)=g_ik(, , ) .
the expressions for the excited - state emission rate @xmath180 and level shift @xmath181 agree with previous derivations  @xcite , if the laser frequency @xmath21 is replaced by the atomic resonance frequency @xmath12 .
the appearance of the laser frequency here is associated with the restriction to the weak - driving limit , which was not considered in earlier derivations . the main new result obtained here
, however , is an explicit expression for the jump operator associated with photon emission , which we will use to calculate photon recoil heating in appendix  [ sec : recoil ] .
here we derive in more detail the equations relating to heating of the atomic motion around the minimum of the vacuum force trap .
recoil heating can be calculated by treating the atomic position @xmath150 in eq .
( [ eq : ojump ] ) as an operator . for notational simplicity
, we shall consider recoil heating only along one direction , say @xmath19 , although our results can be easily generalized . in the main text
, @xmath19 corresponds to the direction in which the atom is trapped  ( normal to the dielectric surface ) .
the increase in momentum uncertainty due to photon scattering is given by = ( p_z^2 _ a ) .
here we will focus solely on the liouvillian term in the density matrix evolution  ( _ i.e. _ , momentum incurred from spontaneous emission ) .
it can be shown that the liouvillian term produces no net force , _
i.e. _ , @xmath182 . using the fact that @xmath183=-i\hbar\frac{\partial f}{\partial z_a}$ ] for any function @xmath184 , one can derive the following general expression for the recoil heating rate in terms of green s functions , = _
0_l^2 _ i _ j ( _ 1 ^ 2 + 2_1_2-_2 ^ 2).[eq : recoil ] here , the derivatives @xmath185 and @xmath186 act on the first and second spatial arguments of the green s functions , respectively , and we have suppressed the atomic spatial variables in the directions that are not of interest .
we further assume that the internal and spatial degrees of freedom can be de - correlated , @xmath187 .
formally , we can write the above equation in the convenient form = ( k_)^2 ( _ a ) p_e , where @xmath188 is the excited - state population .
this form is intuitive as it describes the momentum increase arising from a random walk process , where the atom experiences random momentum kicks of size @xmath189 at a rate @xmath190  ( _ i.e. _ , the photon scattering rate ) .
@xmath191 is a derived quantity that thus characterizes the effective momentum associated with scattered photons in the vicinity of a dielectric surface .    as a simple example
, we can consider an atom in free space polarized parallel to the surface  ( say along @xmath192 ) .
evaluating eq .
( [ eq : recoil ] ) with the free - space green s function @xmath193 from eq .
( [ eq : gfree ] ) yields @xmath194  ( and @xmath195 ) , recovering the known result  @xcite . this recoil momentum is smaller than the full momentum @xmath196 of the scattered photon because the atom emits in a dipole pattern , but only the projection of the photon momentum along @xmath19 contributes to motional heating in that direction .    as a more complicated example
, we can consider the recoil heating near a surface in the trapping scheme considered in the main text .
for the drude model and at small distances between the atom and surface , we find that the scattered component of the green s function dominates the recoil heating , leading to the asymptotic expansion of @xmath197 given in the main text . in all of our numerical calculations ,
the full green s function is used rather than asymptotic results .      as described in the main text , following a spontaneous emission event ,
the spatial wave - function of the atom is temporarily untrapped and evolves under the pure ground - state potential for a characteristic time @xmath198 before the internal dynamics equilibrates .
here we discuss how to calculate an approximate propagator @xmath122 for the wave - function evolving under the ground - state potential .
the problem significantly simplifies if the ground - state potential @xmath199 is linearized around the atom trapping position @xmath200 , h+u_g(z_t)+mg(z - z_t ) , where @xmath201 .
this linearization is justified in our regime of interest where the atom is tightly trapped and the wave - function undergoes little evolution in the time @xmath198 . for simplicity , here we will ignore the constant energy and position offsets and consider the simplified hamiltonian @xmath202 , which corresponds to a mass accelerating under a constant gravitational force .
this observation enables one to solve for the propagator by utilizing the equivalence principle , that the system of interest is equivalent to that of a free particle viewed in an accelerating frame .
one finds that the propagator is given by u(z_f , z_i , t,0)=u_(z_f , z_i , t)(--),[eq : prop ] where @xmath203 is the well - known propagator for a free particle in an inertial frame  @xcite .    in our simulations ,
the ground state wave - function for an arbitrary trapping potential is obtained numerically on a lattice , and so the propagator must be spatially discretized as well . in particular , while the exponential term in eq .
( [ eq : prop ] ) is well - behaved on a lattice , the free propagation @xmath204 is not at short times . we instead replace @xmath204 with the propagator @xmath205 for the free discrete schrodinger equation , = i(_n+1(t)-2_n(t)+_n-1(t ) ) , where @xmath206 , @xmath207 is the lattice site index , and @xmath208 is the lattice constant .
the propagator for this equation ( defined as the matrix satisfying @xmath209 ) is given by d_mn(t)=i^m - nj_m - n(2t)e^-2it , where @xmath210 is the @xmath207-th order bessel function .
the escape rate of an atom from the metastable potential formed by the vacuum trapping scheme due to quantum tunneling is calculated in the wkb approximation in two steps .
first , we find an approximation to the real part of the bound state energy @xmath211 .
this is implemented by replacing the actual metastable dressed - state potential @xmath25 with a stable one , @xmath212 , according to the formula  ( also see fig .
[ fig : wkb ] ) u_s(z ) & = & u_d(z ) ( z > z _ ) , + & = & u_d(z _ ) ( zz _ ) . here , @xmath213 is the position corresponding to the maximum of the metastable barrier .
the ground - state energy of @xmath212 is solved numerically .    once the approximate ground - state binding energy is obtained , the tunneling probability @xmath214 for a particle hitting the metastable barrier is calculated using the wkb approximation , p_t=4(-2_z_1^z_2 dz ) . here
@xmath215 are the solutions to @xmath216 and are thus the classical turning points of a particle with energy @xmath211 . the rate that the trapped atom collides with the barrier is approximated by @xmath217 , where the position and momentum uncertainties are obtained from the numerical ground - state solution .
thus we estimate the inverse of the tunneling - limited lifetime to be @xmath218 .
( in arbitrary units ) seen by the atom .
the maximum of the metastable barrier is located at @xmath213 .
red curve : stable potential @xmath212 used to approximately calculate the binding energy @xmath211  ( blue line ) of the atom in the trap .
the classical turning points @xmath215 of a particle with this energy , which are used to evaluate the tunneling rate in the wkb approximation , are labeled as well.[fig : wkb],width=264 ]


very unstable nuclei with unusual ratios of neutron number ( @xmath2 ) to proton number ( @xmath5 ) are currently the subject of intense theoretical and experimental work . a particularly important issue is the effect on nuclear properties of a large neutron excess and the accompanying low - lying continuum of excited states . in ground
states these conditions sometimes lead to halos , thick skins , and changes in magic numbers .
less is known about excited states , which we study here by calculating transition strengths and densities .
most of the work on excited states is in light neutron - rich nuclei ; a low - energy strength - function peak has been studied for years .
recent work @xcite has focused on enhanced @xmath4 strength in nuclei like @xmath6li or in neutron - rich oxygen isotopes @xcite .
but in heavier nuclei , calculations are more spotty .
several groups have calculated the isovector ( iv ) @xmath4 strength functions , or photoabsorption cross sections , in the ca isotopes @xcite .
references @xcite predict a small peak at around @xmath7 mev in @xmath8ca .
the first paper also calculates the strength function in @xmath9ca , finding no low - energy peak , but ref .
@xcite , a calculation based on density - functional theory , finds a very broad bump at @xmath1010 mev .
it is not clear , therefore , at what value of @xmath2 the low - energy strength is first discernible .
the method we use here , the quasiparticle random phase approximation ( qrpa ) , has been applied mainly to selected @xmath4 excitations in particular isotopes of o , ca , and ni ; one example is ref.@xcite , which asserts that the low - energy excitations reflect strong particle - particle (  dineutron " ) correlations .
references @xcite show extensive and systematic qrpa calculations of photoabsorption cross sections , but focus primarily on radiative neutron - capture cross sections .
the relativistic qrpa @xcite has been applied more with an eye to nuclear structure issues ( for a recent review , see ref .
@xcite ) , asserting , for example , that low - lying isoscalar ( is ) @xmath4 states exhibit toroidal flow . references @xcite go beyond rpa , altering details of predicted low - lying strength .
other channels have been discussed less than the iv @xmath4 .
the authors of ref .
@xcite calculate @xmath3 strength functions in ca isotopes up to the neutron drip line , finding an enhancement at low energy in @xmath8ca .
reference @xcite examines transition strengths in the @xmath3 , @xmath4 , and @xmath1 channels near the neutron drip line , finding , for example , a sharp @xmath1 peak near @xmath11 in @xmath12ni .
the above represents a considerable amount of work on excitations in exotic nuclei .
its range and results can be summarized , roughly , as follows : 1 ) most of the calculations are for spherical nuclei , and the iv @xmath4 channel has been studied the most .
2 ) a low - energy peak grows as @xmath2 increases .
3 ) the low - energy excitations do not seem collective in the @xmath4 channel ( e.g.  ref .
@xcite ) , i.e.  they do not contain more than a few particle - hole or two - quasiparticle configurations .
4 ) the corresponding transition densities have long neutron tails extending beyond the bulk nuclear radius .
5 ) the neutron skin thickness is correlated with the amount of low - energy strength @xcite .
6 ) calculations reproduce measurements of the amount of low - energy strength fairly well if they include many - particle many - hole correlations @xcite .
most of this work has focused on nuclei right at the drip line and , as already noted , has emphasized the iv @xmath4 channel .
it is still not known whether theory makes similar predictions for other neutron - rich nuclei , or to what extent the properties of the low - energy peaks vary from one multipole to another .
the purpose of this paper , therefore , is to calculate strength functions in as many nuclei as possible , in all multipoles up to @xmath1 , and investigate how the properties of the strength functions vary as @xmath2 and the multipolarity change .
we use the qrpa with skyrme density functionals ( which we refer to as interactions from here on ) and volume - type delta pairing interactions .
currently , the qrpa is the most sophisticated method that can be used systematically for resonances of all energies in heavy nuclei without assuming that part of the nucleus forms an inert core . here
we will display results of qrpa calculations for ca , ni , and sn isotopes from the proton drip line to the neutron drip line , describing the corresponding strength functions and transition densities in detail .
we will also examine whether results depend on which skyrme interaction we use , so that we can identify important measurements / nuclei that will help distinguish among them .
section [ sec : calculation ] below describes our method .
section [ sec : strength function ] shows strength functions in the three isotopic chains and looks for @xmath13-dependent trends in their behavior .
section [ sec : transition density ] discusses transition densities to states in low - lying peaks .
section [ sec : conclusion ] is a conclusion .
our approach to the qrpa is described in detail in ref .
@xcite , and we discuss it here only briefly . the starting point is a coordinate - space hartree - fock - bogoliubov ( hfb ) calculation , carried out with the code used in refs .
the code outputs quasiparticle wave functions , with spherical symmetry assumed .
it takes the continuum into account by placing the nucleus in a 20-fm box .
we use as large a single - particle space as possible and solve the hfb equations accurately , so as to ease the removal of spurious states from the qrpa and satisfy sum rules .
we diagonalize the hfb density matrix to obtain a set of canonical states  as many as there are quasiparticle states  that make up our qrpa basis .
the qrpa - hamiltonian matrix then contains off - diagonal elements from the quasiparticle one - body hamiltonian as well as the matrix elements of the residual two - body interaction .
the two - body matrix elements are functional derivatives of the energy with respect to the density and pairing tensor , after diagonalizing the qrpa hamiltonian , we use the eigenvectors to obtain strength functions @xcite and transition densities ( see appendix ) .
we can not handle an infinite number of canonical quasiparticle states , and so truncate our basis in the following way : for @xmath14 the maximum single - particle angular momentum @xmath15 is 15/2 , and for @xmath16 @xmath15 is 21/2 .
if the ground state has a finite pairing gap in the hfb calculation we truncate the canonical quasiparticle basis in the qrpa by omitting states with very small occupation probability @xmath17 . in the @xmath18 and @xmath4 channel
we keep states with @xmath19 , and in the @xmath20 channel we keep those with @xmath21 .
if the mean field has no pairing we use the hartree - fock energy to truncate ; in the @xmath18 and @xmath4 channels we keep single - particle states with energy @xmath22 mev , and in the @xmath20 channel we keep those with @xmath23 mev . in the hfb calculations we always cut off the quasiparticle spectrum at 200 mev .
the largest qrpa - hamiltonian matrix that we encounter here ( in the @xmath4 channel for sn near the neutron drip line ) is about 14,000 , and the minimum is about 2,000 ( in the @xmath3 channel for ca near the proton drip line ) . and
@xmath24 mev were used for the @xmath1 channel , and then the maximum dimension of the hamiltonian matrix was about 20,000 . ]
the separation of spurious excitations from physical states is as complete as in ref .
@xcite in most of cases .
we use 3 skyrme interactions : skm@xmath25 , sly4 , and skp .
our pairing interaction has zero range with no density dependence ( this is `` volume pairing '' ) . for reasons we do not understand , we can not obtain solutions of the hfb equation with skp for ca and ni isotopes .
the other interactions give solutions that are unstable against qrpa quadrupole vibrations for isotopes with : @xmath26 in ni ( skm@xmath25 ) ; @xmath27 ( skm@xmath25 ) and @xmath28 , @xmath29 ( sly4 ) in sn , and so we do not discuss these nuclei .
strength functions are important because they capture information about all excited states and can often be measured .
we define the strength functions through the transition operators    @xmath30 : @xmath31    @xmath32 : @xmath33    @xmath34 : @xmath35 where @xmath36 is total - nucleon number , and @xmath37 is a single - particle index .
the upper limit for @xmath5 ( @xmath2 ) on a summation sign means that the sum is over protons ( neutrons ) .    in the next few sections we show the strength functions , with spurious states removed , produced by skm@xmath25 , and discuss those produced by the other skyrme interactions when appropriate .
we then display energy - weighted sums ( up to several different energies ) produced by all our interactions , in all isotopes .
figures [ fig : str_ca_0+][fig : str_sn_0 + ] show the is and iv strength functions in the @xmath18 channels . in ca ( fig.[fig : str_ca_0 + ] ) , most of the functions are smooth except for small low - energy spikes in several nuclei and the large low - energy peaks in @xmath39ca .
the is giant resonances peak at @xmath40 mev for @xmath41ca and gradually move to 15 mev by @xmath42ca
. as @xmath2 approaches the neutron drip line , a peak near @xmath11 grows , becoming quite large at the drip line .
the iv giant resonance is broad in the light isotopes and develops a low - energy component as @xmath2 increases , apparently increasing the summed strength .
the larger peak varies only slightly in energy over the same range of @xmath2 . near the drip line the low - energy iv strength develops a peak that exactly mirrors the is peak , indicating that the strength in both channels is produced solely by neutrons
this distinct low - energy bump seems to appear first at 10 mev in @xmath43ca .
this nucleus is not as short - lived as the isotopes near the drip line ; perhaps the appearance of the low - energy mode is experimentally testable .
the same nucleus marks the beginning of the low energy peak when we use sly4 .
the @xmath2 dependence of the strength function of ni ( fig .  [ fig : str_ni_0 + ] ) is similar in many ways to that of ca ; the energy of the is giant - resonance energy , for example , falls in a similar way past @xmath44 .
a large low - energy peak again develops near the neutron - drip line , though not as dramatically as in ca , and with the complication that the strength shrinks after @xmath45ni .
this peak first is noticeable in both the is and iv channel at 6 mev in @xmath46ni , here and in calculations with sly4 .
this example and that in ca suggest more generally that the nucleus with 2 neutrons outside a closed shell will mark the appearance of the low - energy peak .    in sn ( fig .
[ fig : str_sn_0 + ] ) , the peak near zero energy around the drip line is much smaller than in ca and ni .
the bump first appears in @xmath47sn at @xmath48 mev ( violating by two nucleons the conjecture made a few lines up ) in both the is and iv channels .
this threshold nucleus is the same with other skyrme interactions .
although we can make these kinds of comparisons in words , we would like to display the results of both ( or all three where possible ) skyrme interactions in a more quantitative and graphic way . to limit the number of our figures , we do so by showing predictions for the sum of energy weighted strength up to 4 successively higher energies .
we define the `` partial '' sum , up to an energy @xmath49 as @xmath50 where @xmath51 is a strength function , eq .
( [ eq : part_esum ] ) overestimates the energy - weighted sum , because the tail in the negative energy region is not included .
this problem occurs near the neutron drip line , and the energy - weighted sum is then overestimated by 10 % at most .
we did not calculate @xmath52 from negative @xmath53 , because then @xmath54 is sometimes negative . ]
( see eq .
( 1 ) of ref .
figure [ fig : part_esum_is_ca0 + ] shows @xmath52 with @xmath55 , 20 , 30 , and 100 mev for the is @xmath3 channel in ca .
all skyrme interactions show similar @xmath2-dependence and @xmath49-dependence ; the difference between the @xmath56 and 100 mev curves is comparable with or smaller than the differences among the other curves ( separated by 10 mev ) , reflecting the eventual vanishing of the strength as @xmath49 increases .
the @xmath40 mev curve increases more rapidly with @xmath2 than the others in the region @xmath57 , indicating that the strength around @xmath40 mev is shifted down at large @xmath2 .
the differences among the interactions are most apparent in the location of the neutron drip line , a static property .
the differences in strength functions , at least in this channel , are for the most part minor ( though large enough to make significant differences in , e.g , the compression modulus ) .
figure [ fig : part_esum_iv_ca0 + ] shows the partial sums in the iv channel . here
the two interactions do yield somewhat different results . though the @xmath2-dependence is similar
, skm@xmath25 produces more integrated strength than sly4 . to see how the total energy - weighted strength can differ from one interaction to another , we examine the iv energy - weighted sum rule ( ewsr ) @xmath58
@xmath59 @xmath60 @xmath61 @xmath62 @xmath63 here @xmath64 is the energy of the excited state @xmath65 , @xmath66 is the ground state , @xmath67 is the nucleon mass , and @xmath68 @xmath69 is the mean square proton ( neutron ) ground - state radius .
the sum over states @xmath70 does not include spurious states .
the second term @xmath71 in eq .
( [ eq : ewsr_iv ] ) , with @xmath72 and @xmath73 of eq .
( [ eq:9 ] ) the proton and neutron ground - state densities , is a correction arising from the momentum - dependent terms of the skyrme interaction .
the skyrme parameters @xmath74 are defined in refs .
@xcite ( skm@xmath25 ) , @xcite ( sly4 ) , and @xcite ( skp ) .
another term in the sum rule depending on the proton - neutron mixed density is always zero in our work .    in tab .
[ tab : ewsr_iv_ca0 + ] we show @xmath75 , @xmath76 , @xmath77 , and ewsr ( the full energy - weighted sum ) for @xmath8ca .
the interactions reproduce ground - state observables about equally well , so the difference in the ewsr , as the table shows , comes from the term involving @xmath78 .
the degree by which the ewsr exceeds @xmath79 is in line with expectations ( see , e.g. , chap .  10 of @xcite and refs .
@xcite ) .
the properties of ni are similar to those of ca .
@xmath80 in the is channel of ni increases more rapidly than the other curves ( fig .
[ fig : part_esum_is_ni0 + ] ) .
skm@xmath25 again has the larger iv sum - rule value .
[ fig : part_esum_iv_ni0 + ] ) .
@xmath80 in the is channel of sn is much higher than @xmath81 ( fig .
[ fig : part_esum_is_sn0 + ] ) , reflecting the lower energy of the giant resonance .
figure [ fig : part_esum_iv_sn0 + ] shows that the iv sum rules follow the same order as in the other isotopes , with skp slightly larger than sly4 .
.[tab : ewsr_iv_ca0 + ] components of the ewsr of iv @xmath3 channel of @xmath8ca . [ cols="^,^,^,^,^ " , ]      is ( solid ) and iv ( dashed ) 0@xmath82 strength functions for even ca isotopes ( skm@xmath25).,scaledwidth=100.0% ]     the same as fig .
[ fig : str_ca_0 + ] but for ni isotopes .
, scaledwidth=100.0% ]     the same as fig .
[ fig : str_ca_0 + ] but for sn isotopes . , scaledwidth=100.0% ]     @xmath2-dependence of the partial energy - weighted sums in the is @xmath3 channel of ca .
the upper limits of the sums are 10 mev ( plus symbol ) , 20 mev ( x ) , 30 mev ( asterisk ) , and 100 mev ( open square ) . , scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_ca0 + ] but for the iv channel .
the open squares now correspond to an upper limit of 150 mev .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_ca0 + ] but for ni .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_iv_ca0 + ] but for ni .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_ca0 + ] but for sn .
, scaledwidth=100.0% ]     the same as fig .
[ fig : part_esum_iv_ca0 + ] but for sn .
, scaledwidth=100.0% ]      having analyzed the @xmath3 strength functions in detail , we will , for the most part , let the higher - multipole figures speak for themselves .
just a few remarks about the @xmath4 channel : the is strength functions in ca appear in fig .
[ fig : str_ca_1- ] . here , unlike in the @xmath3 channel , there is a clear low - energy peak in all nuclei , increasing in sizeca are so sharp because we smooth the strength functions with a smaller width at bound states than in the continuum ( see ref .
@xcite ) . ] after @xmath9ca .
( the is strength functions of @xmath84ca are multiplied by 0.1 in the figure . ) in @xmath42ca 45 % of the energy - weighted strength is below 10 mev .
the iv strength functions also have increasing low - energy strength at large @xmath2 .
the small low - energy bump in @xmath43ca is produced by sly4 as well as skm@xmath85 .
we see similar @xmath2-dependence in ni ( see fig .  [
fig : str_ni_1- ] ) , though after @xmath45ni the low - energy is peak shrinks .
the threshold for rapid growth with both skyrme interactions is the closed - shell nucleus @xmath86ni . in sn , once again , the low - energy is peak grows with @xmath2 , though more smoothly than in the other isotopes ( see fig .  [
fig : str_sn_1- ] ) .
in particular , there is apparently no sudden increase at the closed - shell - plus - two nucleus @xmath87sn .
recently , ref .
@xcite reported a clear iv bump around 10 mev in that nucleus .
the bumps in our calculated spectra do not seem to really set in until higher @xmath2 .
the @xmath4 @xmath52 tell a slightly different story ( figs .
[ fig : part_esum_is_ca1-][fig : part_esum_iv_sn1- ] ) .
the partial sums in the is channel of sn ( fig .
[ fig : part_esum_is_sn1- ] ) show a noticeable kink at @xmath88 , due clearly to an increase in low - energy strength .
as remarked above , the increase is not obvious in the strength - function figures .
the choice of interaction makes very little difference in the @xmath52 except in the blank regions ( indicating quadrupole deformation ) and near the neutron drip line .    according to the analytical argument in the previous subsection , if the iv ewsr depends on the interaction in one multipole , then it should in all multipoles . as figs .
[ fig : part_esum_iv_ca1- ] , [ fig : part_esum_iv_ni1- ] , and [ fig : part_esum_iv_sn1- ] show , this is indeed the case .
we have not commented on the question of whether protons are involved in the transitions to low - lying states in this channel because the is and iv operators do not have the same radial dependence .
we will address the issue in later when we display transition densities .
is ( solid , scale on left ) and iv ( dashed , scale on right ) 1@xmath89 strength functions for even ca isotopes ( skm@xmath25 ) .
the is strength functions in @xmath84ca are reduced by a factor of 10.,scaledwidth=100.0% ]     is ( solid , scale on left ) and iv ( dashed , scale on right ) 1@xmath89 strength functions for even ni isotopes ( skm@xmath25 ) .
, scaledwidth=100.0% ]     the same as fig .
[ fig : str_ni_1- ] but for sn isotopes .
, scaledwidth=100.0% ]     the same as fig .
[ fig : part_esum_is_ca0 + ] but for the @xmath4 channel . ,
scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_iv_ca0 + ] but for the @xmath4 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_ni0 + ] but for the @xmath4 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_iv_ni0 + ] but for the @xmath4 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_sn0 + ] but for the @xmath4 channel . , scaledwidth=100.0% ]     the same as fig .
[ fig : part_esum_iv_sn0 + ] but for @xmath4 channel .
, scaledwidth=100.0% ]      the strength functions in the @xmath1 channels are qualitatively different from those in the lower - multipole channels . in ca , for example ( fig .  [
fig : str_ca_2 + ] ) , all isotopes strength functions of @xmath8ca without low - energy peaks ; we do not know the reasons for the difference between that calculation and ours . ]
except @xmath91 have considerable low - energy strength , in many instances because of sharp quadrupole vibrations .
the low - energy strength grows with @xmath2 , but unlike in the @xmath3 channel , both neutrons and protons are involved in the transitions to these states , even near the drip line , because is and iv distributions do not coincide .
this is so for all skyrme interactions .
in fact the is strength at low energies is much larger than the iv strength in all isotopes except @xmath42ca .
even there , the is strength to the low - energy peak with skm@xmath85 is 256 @xmath92@xmath93 and the iv strength only 136 @xmath92@xmath93 .    the @xmath1 strength functions in ni ( fig .  [
fig : str_ni_2 + ] ) are similar , though with two differences : all the nickel isotopes have low - energy peaks , and for @xmath94ni the iv peak is at a different energy than the is peak .
the lowest - lying peak seems to be nearly pure is , like the familiar surface - quadrupole vibrations in stable nuclei . in the next section
we will discuss the extent to which the states near the neutron - drip line display the same vibrational features we observe near stability .
sly4 puts the drip line at @xmath95ni , much closer to stability , and in isotopes close to @xmath95ni puts the is and iv peaks at about the same location . in the lighter ( @xmath96 ) ni isotopes
, the lowest peak always has a little bit of iv strength mixed with the is strength , as with skm@xmath25 . in @xmath46ni ( skm@xmath25 ) ,
the iv strength in the peak at 400 kev is 5 % of the is strength .
things are quite similar in sn ( fig .
[ fig : str_sn_2 + ] ) ; for @xmath97sn the lowest - energy peak has both is and the iv components but in @xmath98sn it has no iv component .
the other interactions give the same results near the neutron drip line , indicating again that proton excitations are important in these modes everywhere .
the predictions for @xmath52 in ca ( fig .
[ fig : part_esum_is_ca2 + ] ) show that sly4 evinces slightly more of a closed shell ( a small kink in the 10-mev curve ) at @xmath91 than does skm@xmath85 , while showing less of a shell effect at @xmath99 ( no small kink in the 20-mev curve ) . as was the case in the @xmath3 and @xmath4 channels , the iv @xmath52 ( fig .
[ fig : part_esum_iv_ca2 + ] ) depends on the interaction .
similar differences in shell - strength among the skyrme interactions are apparent in the ni and sn @xmath52 curves for the is channel ( figs .
[ fig : part_esum_is_ni2 + ] and [ fig : part_esum_is_sn2 + ] ) , while the difference is not obvious in the iv channels ( figs .
[ fig : part_esum_iv_ni2 + ] and [ fig : part_esum_iv_sn2 + ] ) .
the is @xmath52 for sn ( fig .
[ fig : part_esum_is_sn2 + ] ) indicate that skp has more low - energy strength @xmath100 than the other parameter sets .
the iv curves in sn show interaction differences similar to those in the @xmath4 channel .     is ( solid ) and iv ( dashed ) 2@xmath82 strength functions in even ca isotopes ( skm@xmath25 ) .
, scaledwidth=100.0% ]     the same as fig .
[ fig : str_ca_2 + ] but for ni isotopes .
, scaledwidth=100.0% ]     the same as fig .
[ fig : str_ca_2 + ] but for sn isotopes .
, scaledwidth=100.0% ]     the same as fig .
[ fig : part_esum_is_ca0 + ] but for the @xmath1 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_iv_ca0 + ] but for the @xmath1 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_ni0 + ] but for the @xmath1 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_iv_ni0 + ] but for the @xmath1 channel .
, scaledwidth=80.0% ]     the same as fig .
[ fig : part_esum_is_sn0 + ] but for the @xmath1 channel .
, scaledwidth=100.0% ]     the same as fig .
[ fig : part_esum_iv_sn0 + ] but for @xmath1 channel .
, scaledwidth=100.0% ]
giant resonances lie in the continuum and their @xmath2-dependence is therefore smooth .
the low - energy peaks , on the other hand , can vary dramatically with @xmath2 , especially near the drip line . in this section , to learn more about their collectivity and @xmath2 dependence , we examine transition densities to representative low - lying states .    the transition density @xmath101 is defined in the appendix . here
we focus on the radial transition density @xmath102 the transition amplitude depends on this quantity through the relation @xmath103 where @xmath104 is the radial part of the multipole operator @xmath105 .
the radial transition density therefore specifies the radius at which the strength is large or small . the transition density , as eq .
( [ eq : tr_den ] ) shows , is proportional to the product of particle and hole wave functions and so is localized as long as all hole states are localized , i.e.  as long as the ground state is bound . in this section , as before , we will be using skm@xmath25 unless we state otherwise .
figure [ fig : tr_den_ca0+n16 ] shows the transition densities to 4 states below the giant resonances in @xmath106ca with large strengths ( the energies label the results of the qrpa calculation before smoothing ) .
proton excitations are important , particularly because of the factor @xmath107 in the @xmath3 transition operator .
the rms radii of the protons and neutrons in the ground state are 3.40 and 3.19 fm .
the important part for the strength comes from outside the rms radius .
this is not true of the giant resonances .
figure [ fig : tr_den_gr_ca0+n16 ] shows transition densities in the peaks of the giant resonances of @xmath106ca .
the upper two panels represent the is giant resonance , and as is expected the protons and neutrons have similar transition densities . in the lower panels , corresponding to the iv giant resonance , the two densities are out of phase .
the peaks of the proton transition densities are closer to the center of the nucleus than those of the lower - energy states , particularly in the is channel .
we saw that a small bump appeared after magic numbers in neutron - rich nuclei . in ca ,
the relevant isotope is @xmath43ca .
[ fig : tr_den_ca0+n30 ] ) shows that the excitation is created entirely by neutrons .
but the state is not at all collective .
the largest two - quasiparticle component in the qrpa wave function is @xmath108(@xmath109 mev , @xmath110(@xmath111 mev , @xmath112 , where @xmath113 is the quasiparticle energy in the hfb quasiparticle basis and @xmath114 is the diagonal ground - state occupation probability ( norm of the lower component of the quasiparticle wave function ) .
the radial quantum numbers are assigned in order of @xmath113 ( @xmath115 for those having @xmath116 ) .
@xmath117(@xmath118 ) is the forward ( backward ) amplitude of the qrpa solution .
the squared amplitude @xmath119 of the @xmath120 component is @xmath121 .
the component corresponds to the excitation of a neutron from a @xmath122 orbit near the fermi surface to another @xmath122 orbit about 9 mev above .
the first orbit is strongly affected by pairing correlations .
the two - quasiparticle energy is 10.055 mev , and the energy of the excited state itself is 10.005 mev ; the excitation is clearly almost unaffected by the residual interaction .
yet this nearly pure two - quasiparticle state , because the neutron transition density has such a long tail , carries 5 % of the total is transition strength .
furthermore , protons contribute essentially nothing to the transition ; the maximum of the @xmath119 for protons is @xmath123 .
this remarkable situation  dominance by a single two - neutron - quasiparticle configuration with very large strength , occurs in the @xmath3 channel repeatedly .
take , for example , the drip - line nucleus @xmath42ca .
two low - lying states , at @xmath124 and 2.995 mev , together account for 25.8 % of the total is strength and 10.2 % of the ewsr .
transition densities , shown in fig .  [ fig : tr_den_ca0+n56 ] , are again dominated by neutrons , with outer peaks at 911 fm .
the rms ground - state radii are 3.738 fm for protons and 4.615 fm for neutrons , so the highest ( inner ) peaks of @xmath125 are close to the neutron rms radius .
the main component of the excited state with @xmath124 mev , with a squared amplitude of 0.807 , is @xmath126 mev , @xmath127 mev , @xmath128 ) .
the main component of the state at @xmath129 mev , with a squared amplitude of 0.930 , is @xmath130 mev , @xmath131 mev , @xmath132 .
like before , these are nearly pure two - neutron - quasiparticle states with long tails in their product .
turning now to ni near the closed @xmath133 shell , we see that @xmath46ni has a small bump around 5 mev , while @xmath86ni does not ( fig .
[ fig : str_ni_0 + ] ) . the transition density in fig .
[ fig : tr_den_ni0+n52 ] again shows neutron dominance .
this excited state has 2.8 % of the total is strength , and its main component is @xmath134 mev , @xmath135 mev , @xmath136 , with @xmath137 . at the drip - line nucleus @xmath138ni
, the main component of the excited state at @xmath139 mev , which is one of the states forming a low - energy bump , is @xmath140 mev , @xmath141 ; it has @xmath142 .    in sn the nucleus just past the closed neutron shell
is @xmath87sn . here , as we noted earlier , a new bump is not visible ( fig .
[ fig : str_sn_0 + ] ) , until @xmath47sn .
we show the transition density to the corresponding state in fig .
[ fig : tr_den_sn0+n86 ] .
the excitation s main component is @xmath143 mev , @xmath144 mev , @xmath145 , with @xmath146 .
@xmath147 mev is the smallest neutron quasiparticle energy of the hfb calculation .
figure  [ fig : tr_den_sn0+n126 ] shows the transition densities to three states in the low - energy bump in @xmath148sn .
the behavior is similar to that in the drip - line ca and ni nuclei .
we turn finally to the dependence of the transition densities on the skyrme interaction .
figure [ fig : tr_den_ca0+n16_sly4_skoprime ] shows , unsurprisingly by now , that skm@xmath25 and sly4 give similar transition densities to states below the giant resonances in @xmath106ca . we show transition densities for some giant - resonance states in fig .  [
fig : tr_den_ca0+n16_gr_sly4_skoprime ] .
the is / iv character of those states is clear .
next we compare predictions for @xmath43ca .
figure [ fig : tr_den_ca0+n30_sly4_skoprime ] shows the sly4 transition densities to the states in the small bump in the low - energy continuum .
the bump emerges there with both interactions .
@xmath149ca is the neutron - drip - line nucleus obtained with sly4 ; figure [ fig : str_ca0+n46_sly4 ] displays the strength functions in that nucleus .
apparently they are similar not to those of the skm@xmath85 drip - line nucleus @xmath42ca but rather to those of the skm@xmath85 @xmath149ca ( see fig .  [
fig : str_ca_0 + ] ) . clearly the quantum numbers of the quasiparticle orbits involved in the transition affect the transition density much more than the chemical potential .
figures [ fig : tr_den_ca0+n46_sly4 ] ( sly4 ) and [ fig : tr_den_ca0+n46 ] ( skm@xmath25 ) show the transition densities to the excited states in the shoulder at about 9 mev in @xmath149ca .
the positions of the outer neutron peak are nearly the same in both , but the neutron transition densities for @xmath150 fm have some differences .
the proton transition density is small except in the right panel of fig .
[ fig : tr_den_ca0+n46 ] .
the tails of the transition densities are quite different from those in fig .
[ fig : tr_den_ca0+n56 ] , for the nucleus with 10 more neutrons .
the two - quasiparticle components with the largest amplitude involve either neutron @xmath151 s or @xmath122 s , depending on the state and interaction , and the largest amplitude itself ranges from 0.5 to 0.99 .
proton and neutron transition densities to the four states in @xmath106ca in the peaks of the strength function ( @xmath152 mev ) in fig .
[ fig : str_ca_0 + ] . ]     the same as fig .
[ fig : tr_den_ca0+n16 ] but for the giant resonances .
the states in the upper panels correspond to the peaks of the is giant resonance , and in the lower to the iv giant resonance . ]
transition densities to the excited state in the small low - energy bump in the @xmath3 strength function of @xmath43ca . ]     transition densities to the two excited states in the low - energy peak in the @xmath3 strength functions of @xmath42ca . ]
transition densities to the excited state in the small low - energy bump in the @xmath3 strength function of @xmath46ni . ]     transition densities to the excited state in a small bump in the @xmath3 strength function of @xmath47sn . ]
transition densities to the excited state in a small bump in the @xmath3 strength function of @xmath148sn . ]     transition densities to the @xmath3 excited states below the giant resonances in @xmath106ca with relatively large strength ( sly4 ) . ]
transition densities to the @xmath3 excited states of @xmath106ca in the giant resonances ( sly4 ) . ]     transition densities to the @xmath3 excited states of @xmath43ca in a small bump at low - energy ( sly4 ) . ]
0@xmath82 strength functions of @xmath149ca ( sly4 ) . ]
transition densities to the excited states in the low - energy shoulder ( is ) or peak ( iv ) around @xmath7 mev of the strength function of fig .
[ fig : str_ca0+n46_sly4 ] ( sly4 ) . ]
transition densities to the excited states in the low - energy shoulders around @xmath7 mev and 11 mev in the @xmath3 strength function of @xmath149ca ( skm@xmath25 ) . ]      as we have just seen , the large neutron tails in the transition density are responsible for enhanced low - lying strength near the drip line .
the robust enhancement of low - lying is @xmath4 strength is due in large part to the factor @xmath153 , which emphasizes large radii , in the transition operator .
one stark difference between these transition densities and those of the @xmath3 channel is the role played by protons , which is non - negligible here .
figure [ fig : tr_den_ca1-n30 ] shows the transition densities of two states in the low - lying is peak in @xmath43ca .
protons and neutrons contribute coherently .
the component of the state at @xmath154 mev with the largest amplitude is @xmath155 mev , @xmath156 mev , @xmath157 .
the squared amplitude @xmath119 is only 0.36 .
two other configurations , @xmath158 and @xmath159 have @xmath160 and 0.12 , not much smaller . and since the two - quasiparticle energy of the largest component is 9.099 mev , the excited state gains more correlation energy than those in the @xmath3 channel .
the proton contributions remain even at the neutron drip line .
figure [ fig : tr_den_ca1-n56 ] shows the transition densities to 3 excited states in the large low - energy is peak in @xmath42ca ( fig .  [ fig : str_ca_1- ] ) .
these states have is strengths 9227 @xmath92@xmath161 mev ) , 3219 @xmath92@xmath162 mev ) , and 4856 @xmath92@xmath163 mev ) , representing 42 % , 15 % , and 22 % of the total strength .
while protons play no role for @xmath164 fm , they are active around @xmath165 fm .
the main difference between these states and those of @xmath43ca is that here the neutrons both extend much further out and contribute more inside the nucleus .
it is the tail that is mainly responsible for the increased strength .
interestingly , the states in @xmath42ca are mainly of two - quasiparticle character even though there is a coherent proton contribution .
the component in the state at @xmath166 mev with the largest amplitude ( @xmath167 ) is @xmath168 mev , @xmath169 mev , @xmath170 . in the state at @xmath171 mev
the component with the largest amplitude ( @xmath172 ) is @xmath173 mev , @xmath174 , and at @xmath175 mev it is @xmath176 mev , @xmath177 , with @xmath178 .
the proton two - quasiparticle configuration with the largest amplitude has @xmath179 .    in ni ,
a low - energy peak rises suddenly at @xmath46ni ( fig .
[ fig : str_ni_1- ] ) .
the transition densities of the states in the peak , shown in fig .
[ fig : tr_den_ni1-n52 ] , are similar to those of the neutron - rich ca isotopes .
the density in @xmath138ni , at the drip line , is similar though larger . in sn
there is no sudden increase in the low - energy peak around @xmath88 ( see fig .
[ fig : str_sn_1- ] ) . fig .
[ fig : tr_den_sn1-n80_n84 ] shows the transition densities to low - energy states in three isotopes around @xmath88 . here
the proton contributions stretch to larger radii than in ca or ni , even though they are dwarfed by the neutron tails there .
the other skyrme interactions give very similar results , both in the is and iv channels and for all three isotope chains .
transition densities to the excited states in the low - energy peak of the is @xmath4 strength function of @xmath43ca . ]
transition densities to the excited states in the low - energy peak of the is @xmath4 strength function of @xmath42ca . ]
transition densities to the excited states in the low - energy peak of the is @xmath4 strength function of @xmath46ni . ]
transition densities to the excited states in the low - energy peak of the is @xmath4 strength function of @xmath180sn ( upper panels ) , @xmath181sn ( middle ) , and @xmath87sn ( lower ) . ]
figures [ fig : str_ca_2+][fig : str_sn_2 + ] showed that at the neutron drip line the is and iv strength functions in ni and sn have distinct low - energy peaks ( though the is strength is always larger , even at the low energy iv peak ) while in ca the peaks coincided .
this is true , except in ni with sly4ni has a small iv peak at the is peak .
] , no matter which skyrme interaction we use . with skm@xmath85 ,
the is strength of the large low - energy peak at @xmath182 mev in @xmath138ni is 92 % of the total strength ( 8.9 % of the ewsr ) , and that of the peak at @xmath183 mev in @xmath148sn is 23 % of the total strength ( 2.2 % of the ewsr ) . in fig .
[ fig : tr_den_sn2+n126_ni2+n70_ca2+n56 ] we compare the transition densities of the states in these low - energy peaks with those in @xmath42ca . as fig .
[ fig : str_ca_2 + ] suggests , though the protons play a role , the neutrons are proportionately more important in ca and do not seem to correlate with the protons . in @xmath138ni and @xmath148sn , by contrast , the neutrons and protons are completely in phase ( as they must be if the iv strength is negligible ) .
the other skyrme interactions bear out these conclusions , though in the sly4 calculation of the drip - line ca isotope the ratio of proton transition density to neutron transition density is larger than that with skm@xmath25 ( shown in the left panel of fig .
[ fig : tr_den_sn2+n126_ni2+n70_ca2+n56 ] ) . in all nuclei near the neutron drip line , the states in low - energy peaks are mostly above neutron - emission threshold , yet the neutron tails cut off at much smaller radii than to those in the @xmath3 and @xmath4 channels .
in addition , the transition densities in those other channels generally have a node at about @xmath165 fm , while the @xmath1 transition densities show no such feature .
we put off for just a moment the obvious question of whether these low - lying is peaks are familiar quadrupole vibrations and turn briefly to the iv channel .
figure [ fig : tr_den_2_sn2+n126_ni2+n70 ] shows the transition densities to states in the low - energy iv strength - function peak in @xmath148sn and @xmath138ni .
the states in @xmath148sn display remarkably is character , despite being at the peak of the iv distribution .
the fact that the iv distribution peaks there is a little misleading because the is strength is still considerably larger , even though it is larger still at lower energies .
we get some iv strength because the ratio of the transition density of protons to that of neutrons is not exactly equal to @xmath184 ( see eq .
( [ eq : transition_operator_2 + ] ) ) .    returning to the
is strength , we note , in accord with fig.[fig : tr_den_sn2+n126_ni2+n70_ca2+n56 ] , that the transitions appear more collective in the ni and sn drip - line isotopes than in that of ca .
the component with the largest amplitude ( @xmath185 ) in the state at @xmath183 mev of @xmath42ca is @xmath186 mev , @xmath187 mev , @xmath188 , and that with the second largest ( @xmath189 ) is @xmath190 .
this is still a predominantly single two - quasiparticle excitation .
by contrast , the components of the state with at @xmath182 mev in @xmath138ni with the largest amplitudes are @xmath191 while still not very collective , this state is a little more so than its counterpart in ca .
( @xmath138ni may be in the spherical - deformed transition region . ) in @xmath148sn , the component of the state at @xmath192 mev with the largest amplitude is ( fig .
[ fig : tr_den_sn2+n126_ni2+n70_ca2+n56 ] ) @xmath193 mev , @xmath194 mev , @xmath195 , ( @xmath196 ) and that with the next largest is @xmath197 mev , @xmath198 mev , @xmath199 , ( @xmath200 ) .
this excitation is mainly a single neutron particle and neutron hole .
still , the proton contribution is much larger than that in other channels at the neutron drip line .
this excited state has 23.7 % of the total is strength and 8.6 % of the ewsr .
we see both in this channel and others that the large transition strength near the neutron drip line does not always indicate the coherent contributions of many two - quasiparticle excitations ( i.e. , collectivity ) . yet
the transition densities make these strong low - lying states , particularly in ni and sn , look very much like the surface vibrations that characterize nuclei closer to stability . to get a better handle on the degree to which the states we examine here can be called vibrations , we show in fig .
[ fig : n - e - be2-largest_sn_ni_ca2 + ] the energy , transition probability @xmath201 , and largest @xmath119 for neutrons in the lowest - energy @xmath1 states of ca , ni , and sn .
the ca panel shows , first of all , that we do not reproduce experimental data very well except in a few nuclei , and second , that the states are dominated by a single two - quasiparticle excitation for @xmath202 . for @xmath203
something peculiar happens ; when the energy increases , so does the @xmath201 , and vice versa , with the largest @xmath119 , a measure of collectivity , varying irregularly . in nuclei with collective vibrations , transition strength is usually correlated inversely with energy @xcite ; we see the opposite here . near the neutron drip line of ni ,
the classic signs of increasing quadrupole collectivity  decreasing energy , rising @xmath201 , and decreasing maximum @xmath119  are all present .
the sn isotopes look much more familiar , exhibiting the usual trends , except at the magic numbers . and the states close to the drip line with @xmath204 are as classically collective as ( or more so than ) those in more stable isotopes with @xmath205 .
finally we show the interaction dependence of the collectivity measures for the lowest @xmath1 states in sn ( the easiest chain to understand ) in fig .
[ fig : n - e - be2-largest_sly4_skp_skoprime ] . in the mid - shell region , the energies are low , the @xmath201 s are high , and the largest amplitudes low , as expected from systematics in stable nuclei in regions of large dynamical deformation .
the values of these quantities depend significantly on the interaction , however .
ca and ni show similar behavior but with more irregularities .
transition densities to the excited states in the low - energy peak of the is @xmath1 strength functions in the neutron - drip - line nuclei ( wth skm@xmath25 ) @xmath42ca , @xmath138ni , and @xmath148sn .
both proton and neutron curves are reduced by a factor of 5 in @xmath138ni to fit on the figure .. ]     transition densities to the excited states in the low - energy peak of the iv @xmath1 strength functions ( skm@xmath25 ) in @xmath148sn and @xmath138ni . ]
@xmath2-dependence of the energy ( upper ) , @xmath201 ( middle ) , and largest @xmath119 ( lower ) of the lowest @xmath1 states in ca , ni , and sn , with skm@xmath25 .
experimental data are also shown where available .
, scaledwidth=100.0% ]     the same as fig .
[ fig : n - e - be2-largest_sn_ni_ca2 + ] but for sn only , sly4 and skp . ]
the feature that leaps out of our calculations is the enhancement of low - energy strength as the neutron drip line is approached .
it is present in all multipoles ( and is most pronounced in the is @xmath4 channel , where the transition operator contains a factor @xmath153 ) , with all skyrme interactions , and in every isotopic chain , beginning just after a closed neutron shell .
interestingly , the recent measurement of an iv @xmath4 strength function at gsi shows a low - energy peak before a closed shell , in @xmath180sn .
this may be due to correlations not included in the qrpa .
the transition densities to the low - lying states reveal some interesting features . in the @xmath3 and @xmath4 channels ,
the low - energy peaks are usually created essentially entirely by a single two - neutron - quasiparticle configuration . in the @xmath3 channel , where the total number of two - quasiparticle configurations is restricted , even small admixtures of two - proton - quasiparticle states
are absent ; the proton transition density is nearly zero .
protons play a noticeably larger role in the transition density of the @xmath4 channel despite the lack of any two - proton - quasiparticle configurations with sizable amplitude . in both these channels ,
the transition densities have a peak around the rms radius and very long tail outside the radius , which allows single quasiparticle configuration to have such large strength .
the transition densities to @xmath1 states are different ; they have no real nodes and a proton component that is of the same order as the neutron component .
the transition densities near the neutron - drip line resemble those of surface vibrations , and in many cases the corresponding states clearly are vibrational .
but the usual systematics of collective vibrations is not always present ; energies and @xmath492 strengths are sometimes correlated in the wrong way .
some of the most significant differences in the predictions of the three skyrme interactions we used are in static properties : the locations of the neutron drip line and quadrupole - deformed regions .
but some dynamical properties also depend significantly on the interaction , most notably ( and easiest to measure ) the energies of and transition - strengths to the lowest @xmath1 states near regions of deformation and the total energy - weighted iv strength .
this paper , in a way , is a catalog of predicted dynamic properties of nuclei near the drip line . for the most part we have not offered simple explanations for the anomalies that appear .
that , we hope , will be the subject of a future paper .
this work was supported in part by the u.s .
department of energy under grant de - fg02 - 97er41019 .
we used computers belonging to the national center for computational sciences at oak ridge national laboratory and information technology services at university of north carolina at chapel hill .
parts of this research were done when one of us ( j.t . ) was at the university of north carolina at chapel hill and at riken .
the transition density is defined by @xmath206 where @xmath66 and @xmath65 denote the ground and excited states of the system , and @xmath207 in the right - hand side is an operator acting on the nucleon @xmath37 .
the transition density is evaluated in the qrpa ( we omit @xmath208 ) as @xmath209 where @xmath210 and @xmath211 are the forward and backward amplitudes of the qrpa solution @xmath70 ( see appendix a of ref .
@xcite ) , and @xmath212 stands for a single - particle canonical - basis wave function .
the @xmath213 factors associated with the basis are represented by @xmath214 and @xmath215 .
assuming spherical symmetry for the system and using @xmath216 with @xmath217 a set of quantum numbers , @xmath218 a radial wave function , and @xmath219 a spin wave function , one can derive @xmath220 where @xmath221 , and @xmath222 denotes the angular momentum of the state @xmath70 . @xmath223 and @xmath224 are the forward and backward amplitudes in the spherical representation .