the possible modification of nucleon properties in the nuclear medium is currently a much discussed topic in low energy hadron physics .
one way to consider such problems is to describe the nucleon as a topological soliton and then study the influence of baryonic matter on the properties of such solitons , making use of the independent particle picture which has been so successful in describing many properties of nuclei .
there exist already some works where nucleons described as skyrme  type solitons embedded in infinite nuclear matter have been considered  @xcite .
the results of these studies were in qualitative agreement with experimental indications and with results of other authors using different approaches .
on the quantitative level , however , there is a too large renormalization of the nucleons effective mass in nuclear matter .
the quantitative value of this renormalization is about @xmath2 for the normal nuclear matter density , @xmath3  @xcite .
it is therefore difficult to relate this modification of the nucleon self  energy in the medium to the nucleon mass in free space because such calculations consider only the special case when a nucleon is placed in the center of the heavy nucleus .
the infinite nuclear matter approach allows one only to consider properties of nucleons placed near the center of heavy nuclei where the density is constant . one can expect that taking into account non - spherical effects , i.e. deformation of the skyrmion in the finite nucleus , would improve the results also on the quantitative level
density changes play an important role when the nucleon is placed at sufficiently large distances from the center of nucleus .
the role of deformation effects ( for an early work see  @xcite ) has already been investigated in connection with the nucleon
nucleon interaction in ref .
the results presented in this work demonstrated the necessity of non - spherical deformed solutions of the field equations to obtain a reliable description of the @xmath4 potential . a similar approach to the deformation for two interacting skyrmions has been used in the work  @xcite , where changes in the nucleon shape are investigated by letting the nucleon deform under the strong interactions with another nucleon . on the other hand , rotational and vibrational excitations of deformed skyrmions
have also been considered for baryon number equal to one ( @xmath5 ) system  @xcite . in a series of papers  @xcite axially symmetric few  skyrmion systems with baryon number
@xmath6 larger than one were investigated , because spherically symmetric configurations in the skyrme model do not give bound skyrmionic systems  @xcite even for the solutions of the model within the finite nuclei as considered in ref .
@xcite ( for similar or related work , see also ref .
@xcite ) . in the works  @xcite the deformation
is introduced by a modification of the unit spherical vector and additionally through the modification of the profile function  @xcite in the expression of the ansatz for the chiral field ( for more precise definitions , see below ) .
the results of all investigations for @xmath7 systems lead to the conclusion that the system has to be deformed to obtain the ground state .
some other papers that deal with deformation effects and non - spherical configurations are collected in refs .
@xcite .    in the present work we consider properties of the deformed nucleon embedded in light , medium  heavy and heavy nuclei .
deformation effects are introduced by the distortion of the profile function of chiral field under the action of the external field ( which parameterizes the baryonic density within a given nucleus ) .
in contrast to most previous investigations , we do not consider a constant nuclear density as appropriate for nuclear matter but rather the distance - dependent density as given for finite nuclei .
we calculate the modifications of the nucleon properties in finite nuclei , and their dependence on the distance between the topological center of the skyrmion inside the nucleus and the geometrical center of the nucleus under consideration will be considered .
the paper is organized in the following way .
after the formulation of the problem in section  [ sec : form ] , we discuss the classical soliton mass and distribution of baryon charge for a skyrmion inside finite nuclei in section  [ sec : mass ] . based upon the results for the mass and baryon charge distribution of the deformed classical soliton , we present in section  [ sec : quant ] the quantization of axially symmetric systems with baryon number @xmath5 .
section  [ sec : prop ] is devoted to the quantum properties of the solitons , in particular the effective mass of nucleons in finite nuclei and the consequences for @xmath0 mass splittings .
we also consider the density dependence of the nucleons electromagnetic radii and calculate the intrinsic quadrupole moment of the deformed nucleon .
finally , in section  [ sec : concl ] the conclusions and an outlook for further investigations are given .
our starting point is the medium
modified skyrme lagrangian  @xcite @xmath8 ^ 2 + { { \displaystyle}\frac}{f_{\pi}^2m_{\pi}^2}{16}\alpha_s(\vec r)\mbox{tr}[u+u^{\dagger}-2]\,\ , , \label{lag } \end{array}\ ] ] where @xmath9 , @xmath10mev is the pion decay constant , @xmath11 is a dimensionless parameter and @xmath12mev is the pion mass .
the medium functionals @xmath13 and @xmath14 are expressed via the functionals @xmath15 and @xmath16 of @xmath17-wave and @xmath18-wave pion
nucleon scattering lengths and the nuclear density @xmath19 : @xmath20 here , @xmath16 has the meaning of the pion dipole susceptibility of the medium .
the influence of the medium enters through these functionals ( for explicit expressions , see section  [ sec : sus ] ) .
note that we employ here the same approximations as in  @xcite .
since the fourth order term can be obtained from the infinite mass limit of @xmath21 meson exchange , a possible density dependence of that term could enter via in  medium modified @xmath21 meson properties as it was done e.g. in  @xcite .
similarly , one could add a sixth order term @xmath22 arising from @xmath23 meson exchange or an symmetric fourth order term due to scalar mesons .
we do not consider such terms here .
also , since the empirical information on possible changes of vector meson properties in the medium has not yet lead to an established picture , we treat the coupling @xmath11 as a density independent number .
this could be released in more sophisticated model but goes beyond the scope of the present investigation . for a more detailed discussion about how the specific form of the lagrangian in eq .
( [ lag ] ) arises , we refer to ref .
@xcite .
the ansatz for the chiral field @xmath24 has the hedgehog form @xmath25 \label{anzsph}\ ] ] for the skyrmion placed at the center of the nucleus .    when the skyrmion is placed at a distance @xmath26 from the center of the spherical nucleus ( see fig .
[ fig1 ] ) , the chiral field @xmath24 can be written in the form @xmath27\ , , \label{ansatz}\ ] ] where @xmath28 is the unit isotopic vector correlated with the unit space vector @xmath29 . in the more general case , @xmath30 is a functional of some angular functions @xmath31 , @xmath32  @xcite . the dependence @xmath33 $ ] and the dependence of the profile function @xmath34 on the direction of the radial co - ordinate space vector @xmath35 characterizes the deformation of the skyrmion . as it was mentioned above , in the works  @xcite such deformation
was taken into account when considering properties of multi - skyrmionic systems .
these modifications are dictated by the topological properties of the model , in particular the @xmath36 dependence is dictated by the constraint that the baryon number is larger than one , @xmath7 .
here we consider the @xmath5 case and simply set @xmath37 .
consequently , we can represent the vector @xmath28 in the following form @xmath38 however , this is not the only deformation effect .
we also modify the profile function @xmath34 by allowing for a @xmath39 dependence , @xmath40 note that any azimuthal @xmath41 dependence of @xmath34 is removed due to the axial symmetry of the system ( see the next subsection ) .
obviously these deformation effects also exist in the case when a skyrmion is placed at the center of nucleus , @xmath42 , but in this case the deformation has only a radial dependence . in the previous work  @xcite this dependence
was determined by solving the equations for the profile function of the ansatz  ( [ anzsph ] ) in nuclear matter , @xmath43 .
it is clear that in the present case the deformation has in general both a radial and an angular dependence .
the angular dependence will come into play when the skyrmion is located at some distance from the center of the nucleus .
the radial part of the deformation represents the breathing mode of the skyrmion and the angular part represents its deviations from the spherical form . finally , taking into account our choice of @xmath28 and @xmath34 ( eqs .
( [ greattheta ] ) and  ( [ profform ] ) ) , the following ansatz for the chiral field @xmath24 arises @xmath44~. \label{ansatzlast}\ ] ] this form will be used in the present work . as a result some modifications beyond the one induced by the breathing mode will appear .
for example , one expects that an intrinsic quadrupole moment of the skyrmion emerges .
this will be discussed in more detail below .    before discussing the in - medium skyrmion properties ,
we point out that throughout this paper we neglect the effects due to fermi motion of the solitons ( nucleons ) . in a strict large-@xmath45 counting ,
such effects are suppressed . however , since we do not expect deformation effects to be large , one can doubt the validity of this formal suppression . at present , it is simply not known how to include such effects systematically because one not only would we have to boost the solitons ( where the boost velocity could be obtained in the local density approximation ) but also the coupling to the rotational modes would considerably complicate the quantization procedure to obtain physical states . while the neglect of fermi motion might to some extent alter the quantitative results of our study , we believe , however , that the qualitative ones would remain unaltered .
the mass functional in the static case can be obtained from the lagrangian  ( [ lag ] ) by using the ansatz  ( [ ansatzlast ] ) @xmath46d^3\vec r\ , .
\end{array}\ ] ] to simplify the calculations one introduces the variable @xmath47 putting the origin of the coordinate system @xmath48 at the topological center of the skyrmion @xmath49 ( see fig .
[ fig1 ] ) .
then the expression for the mass functional can be written as @xmath50d^3\vec r\,'\ , .
\label{massr } \end{array}\ ] ] we note that the distance @xmath51 between the skyrmion and the nucleus is a matter of convention .
it depends on how one determines the center of the soliton .
if the skyrmion is deformed its mass and topological centers may not coincide . in this paper
we determine @xmath51 as the distance between the center or the nucleus under consideration and topological center of the skyrmion under consideration .    from figure  [ fig1 ]
it is easy to see that the system has an axial symmetry .
one can choose the direction of the @xmath52 axis along the direction of the vector @xmath26 so that the integrand in  eq .
( [ massr ] ) has no azimuthal @xmath53 dependence .
consequently the profile function @xmath34 only depends on the variables @xmath54 and @xmath55 .
therefore , the lagrangian  ( [ lag ] ) in the static case can be written in the form @xmath56\alpha_p(r\,',\theta',r)\\
\quad \\\qquad\qquad -{{\displaystyle}\frac}{1}{4e^2 } [ ( \partial_i\vec n\cdot\partial_i \vec n)^2- ( \partial_i\vec n\cdot\partial_j \vec n)^2 ] \sin^4f\\ \quad \\\qquad\qquad -{{\displaystyle}\frac}{1}{2e^2 } [ ( \partial_i\vec n\cdot\partial_i \vec n)\partial_j f\partial_jf- ( \partial_i\vec n\cdot\partial_j \vec n ) \partial_i f\partial_jf ] \sin^2f\\ \quad \\\qquad\qquad -{{\displaystyle}\frac}{m_\pi^2f_\pi^2}{4}(1-\cos f)\alpha_s(r\,',\theta',r)~ , \label{lstatic } \end{array}\ ] ] which ( after some algebraic manipulations ) leads to following mass functional for convenience . ] @xmath57 f_\theta^2+{{\displaystyle}\frac}{m_\pi^2}{e^2f_\pi^2}\left(\sin{{\displaystyle}\frac}{f}{2}\right)^2 \tilde r^2\alpha_s(\tilde r,\theta,\tilde r)\\ \quad \\\qquad \qquad \left
. + \left[\sin^2f\left({{\displaystyle}\frac}{\sin^2\theta}{\sin^2\theta}+\theta_\theta^2\right ) + { { \displaystyle}\frac}{\tilde r^2\alpha_p(\tilde r,\theta,\tilde r)}{4}\right ] f_{\tilde r}^2\right\ } \sin \theta d\theta d\tilde r\ , .
\label{masstat } \end{array}\ ] ] here @xmath58 , @xmath59 , @xmath60 and we have introduced the dimensionless variables @xmath61 , @xmath62 .
the corresponding baryon number of the skyrmion is given by @xmath63    the static soliton solution is found by minimizing the functional  ( [ masstat ] ) .
for doing that , one has to solve a system of coupled differential equations @xmath64 for the functions @xmath65 and @xmath31 subject to the boundary conditions determined from the baryon number condition @xmath5  ( [ baryonnumber ] ) . to simplify the numerical computations we perform various approximations .
first , consider the profile function @xmath66 for an undeformed ( spherical ) skyrmion .
it is well known that @xmath66 can be parameterized to a very good approximation as @xmath67 for the free skyrmions with baryon number @xmath5 .
the value of the parameter @xmath68 is obtained from the constraint of having the minimum mass functional for the free skyrmion .
this parameterization in principle can be used for a bound skyrmion and we apply this ansatz in constructing our profile function ( despite some short - coming , see section  [ sec : mass ] ) . also ,
from previous calculations we know that the dependence of the profile function on the density is weak , see e.g.  @xcite .
renormalization of the mass functional  ( [ masstat ] ) mainly occurs via the medium functionals @xmath69 and @xmath70 .
furthermore , from fig .
( [ fig1 ] ) we see that changes of @xmath55 mainly affect the value of medium functionals and @xmath70 are functionals of the medium density @xmath71 . ] and following the results from ref .
@xcite , we assume changes of profile function to be small . using these facts we can represent our profile function as @xmath72\right\}~ , \label{proflast}\ ] ] where the @xmath73 functions are chosen to maintain periodicity in @xmath39 . for the function @xmath74 we choose the parameterization @xmath75 in order to avoid singularities from factors like @xmath76 in the expression of the mass functional  ( [ masstat])@xcite . in eqs .
( [ proflast ] ) and ( [ thetalast ] ) , @xmath68 , @xmath77 , @xmath78 , @xmath79 , @xmath80 and @xmath81 , @xmath82 , @xmath83 , @xmath80 are variational ( or in other words deformation ) parameters .    our aim is the minimization of mass the functional  ( [ masstat ] ) using the expressions  ( [ proflast ] ) and  ( [ thetalast ] ) for the profile function . before doing that
, we have to specify the medium dependence that enters through the @xmath17- and @xmath18-wave pion - nucleon functionals @xmath84 and @xmath85 , respectively .
the medium functionals @xmath84 and @xmath85 introduced in eqs .
( [ lag ] ) ,  ( [ medfunc ] ) have the form  @xcite . for asymmetric nuclei
, additional isovector effects will be present but are not considered here . ]
@xmath86 where @xmath87 is a kinematical factor , @xmath88mev the mass of the nucleon and @xmath89 is the migdal parameter which takes into account the short
range pair correlations of the dipole centers .
the empirical parameters @xmath90 , @xmath91  @xcite can be taken from analyses of low energy pion - nucleus scattering data .
we choose the following parameterization for the density of the finite spherical nuclei considered here  @xcite : @xmath92 \exp\left\{-{{\displaystyle}\frac}{r^2}{r_{0}^{2}}\right\}\,,\quad & \mbox{for $ 4 < a < 20$}\,,\\ \quad \\
\rho(r )   = \left({{\displaystyle}\frac}{a-1}{a}\right){{\displaystyle}\frac}{\rho_0 }   { 1+\exp\left\{{{\displaystyle}\frac}{r - r'}{a}\right\}}\ , , & \mbox{for $ a\ge 20$}\ , .
\label{density } \end{array}\ ] ] in the actual calculations ( minimization procedure ) we have considered the nuclei @xmath93c , @xmath94o , @xmath95ca , @xmath96fe , @xmath97au and @xmath98pb . the first two represent _ light _ , the next two _ medium - heavy _ and the last two _ heavy _ nuclei .
the parameter @xmath99 depends on the type of nucleus considered .
its value is @xmath100  fm for @xmath93c and @xmath101  fm for @xmath94o .
@xmath102 is the charge of the nucleus and @xmath3 is the normal nuclear matter density .
furthermore , @xmath103  fm and @xmath104  fm .
@xmath105 is the baryon number of the nucleus and the factor @xmath106 allows one to consider properties of a given nucleon which is one of the nucleons of the nucleus . for the accuracy of our calculations
, such a simple rescaling of the nuclear density should be sufficient .
the density is pameterized as @xmath107 as dictated by the axial symmetry of the system .
for the input parameters of the skyrme model we use @xmath108  mev and @xmath109 in order to reproduce the free space masses of the nucleon and the @xmath0 as in the seminal work  @xcite .
one could also work with the empirical value of @xmath110 and adjust @xmath11 to the @xmath111 mass splitting . in that case , however , the nucleon mass would be too large , but relative changes would be affected considerably less  @xcite . for definiteness , we prefer here to work with the physical nucleon mass .
our calculations have been performed using different sets of the variational parameters @xmath68 , @xmath77 , @xmath78 , @xmath112 , @xmath82 , @xmath113 .
we represent the skyrmion size parameter @xmath68 as @xmath114 where @xmath115 is the value for the skyrmion in free space ( @xmath116  fm ) .
the minimization procedure of the mass functional  ( [ masstat ] ) showed that considering the set of parameters @xmath117 is sufficiently accurate .
the strength of the parameter @xmath118 decreases with increasing index @xmath119 .
for example , @xmath79 is one order of magnitude smaller than @xmath77 at the distances where deformation effects are maximal .
the parameter @xmath120 decrease very fast with increasing index @xmath121 .
for example , @xmath82 is one order of magnitude smaller than @xmath81 and @xmath83 is almost zero . for controlling our results and keeping a high accuracy ,
we have performed calculations with the following set of parameters : @xmath122 for the reasons just mentioned , we only show results of @xmath77 , @xmath78 and @xmath81 .    the dependence of the skyrmion mass from the distance between the centers of the nucleus considered and the soliton ( which represents one of the nucleons within the nucleus ) is shown in figs .
[ co ] ,  [ cafe ] ,  [ aupb ] , for the light , medium - heavy and heavy nuclei , in order .
the bands labeled `` 1 '' in these figures represent the soliton mass normalized to its free space value and for convenience we display by the bands denoted `` 2 '' the nuclear density normalized to its value at the center of the nucleus .
these bands are shown in the top panels of these figures .
the behavior of the soliton mass is similar for all cases when the soliton is placed within various nuclei , i.e. the mass increases with falling medium density and approaches its free space value at the border of the nucleus .
it has its minimum value at the center of the nuclei .
this is even the case for the light nuclei , which do not have their maximum density exactly in their center .
the deformation parameters which represent the deviations of the nucleon shape from the spherical one are presented in the bottom panels of figs .
[ co ] ,  [ cafe ] ,  [ aupb ] .
the bands denoted `` 3 '' represent the dependence of the relative value of the radial deformation parameter @xmath123 from soliton location inside nucleus and the bands `` 4 '' and `` 5 '' give the dependence for @xmath77 and @xmath78 , respectively .
the band `` 6 '' represents the parameter @xmath81 multiplied by factor 10 .
the other deformation parameters are very small and we do not show them here .
the leading parameter @xmath77 has its maximum value at the distances where the density of the nuclei falls off by @xmath12420% for @xmath93c and @xmath94o , @xmath124 40% for @xmath95ca , @xmath96fe , @xmath97au and @xmath98pb , whereas the radial deformation parameter @xmath125 has its maximal value at the center of the nuclei .
the parameters @xmath78 and @xmath81 have maximum values where the density of the nucleus falls off @xmath12480% . in these regions , the density gradient is maximal leading to the largest angular deformation . from this one
concludes that the skyrmions inside nuclei is in a deformed state .
the topological baryon charge distribution for a skyrmion embedded in the nuclei @xmath93c and @xmath96fe is shown in figs .
[ c12bc ] ,  [ fe56bc ] , respectively , where figures * a * ) represent the baryon charge distribution in the @xmath126 plane and figures * b * ) give its projection onto this plane .
each area at the bottom figure represents the projection of surfaces of equal baryon charge density .
the point @xmath127 indicates the topological center of the skyrmion and the center of nucleus is situated at the negative side of the @xmath128 axis .
one can see that opposite directions in the @xmath128 direction are not equivalent in this plane
. stated differently , going for the same distance in opposite directions in @xmath128 gives different values of the baryon charge density .
we remark that the holes at the origin of the system ( @xmath127 ) are caused by the parameterization of the profile function ( see eq .
( [ proflast ] ) ) . in this approximation , @xmath129 as @xmath130 .
this derivative takes a finite value if one directly solves the equation coming from minimization of the mass functional  ( [ masstat ] ) not using the ansatz eq .
( [ profpar ] ) , thus no such holes would appear . however , the effects discussed so far and in the following are genuine and independent of the ( dis)appearance of these computational artifacts . the light area at positive values of the axis @xmath128 is the area of maximal density .
the mass distribution , i.e. integrand in the expression  ( [ masstat ] ) shows a behavior similar to the one of baryon charge distribution .
an exotic behavior of the mass and baryon charge distributions is observed when the skyrmion is placed at sufficiently large distances from the center of the nucleus , where density changes are significant ( see figs .
[ c12bc ] ,  [ fe56bc ] ) .
one can see that the maximum of the baryon density is shifted away from the center of the nuclei .
when the skyrmion is located at the center of nucleus , the baryon density distribution is symmetric and thus the shift becomes more effective with increasing of @xmath131 .
the shift of the density maximum from the topological center of the skyrmion in opposite direction to the direction of center of the nucleus can be explained in the following way .
the mass of skyrmion decreases in the nuclear medium and if the density is decreasing , the skyrmion mass increases . while we are going from the center of the nucleus to its surface the density fall - off and
mass concentration of the skyrmion occurs at the side of the smaller nuclear density as seen from its topological center .
a similar behavior of the baryon charge density for other nuclei is observed , so we do not present it here .
finally , at the end of this section we remark that rotations of the figures * b ) * around the @xmath128axis give three dimensional surfaces with an equal baryon charge distribution . from this
we conclude that the soliton has two equal components of the moment of inertia which are not equal to the third one , @xmath132 .
we will use this fact in the next section .
as it is well known , baryons emerge in the skyrme model as quantized solitons  @xcite .
this quantization is performed in terms of an adiabatic rotation of the soliton and gives rise to the appropriate quantum numbers for the emerging physical states .
in such an approach , the grand spin @xmath133 , i.e. the sum of angular momentum @xmath134 and isospin @xmath135 , is conserved .
this leads to a tower of states with @xmath136 .
it is well known that in this case only the nucleon ( @xmath137 ) and the @xmath0 ( @xmath138 have physical meaning .
all other particles ( higher rotational states ) should be considered artifacts of the model and will be discarded .    to quantize the solitons of our model we use the standard canonical quantization procedure performing time
dependent rotations in isospace and in co - ordinate space @xmath139 where the matrix @xmath140 specifies the rotation in isospace and @xmath141 the one in co - ordinate space  @xcite .
the parameters of the rotation matrices @xmath140 and @xmath142 serve as collective variables describing the rotational degrees of freedom with respect to a body - fixed system of reference .
the latter has its origin at the topological center @xmath49 of the skyrmion and points in the direction of the distance vector @xmath26 from the center of the nucleus to the center of the skyrmion ( see fig .
[ fig1 ] ) .
the quantization procedure used here is similar to the standard one  @xcite but now the axial symmetry of the system should be taken into account  @xcite .
the angular velocities of isorotations @xmath143 and space rotations @xmath144 can be determined from the respective equations @xmath145 where dot means differentiation with respect to the time component , and @xmath146 is the totally antisymmetric tensor in three dimensions .
using now the anstze  ( [ ansatzlast ] ) ,  ( [ rotansatz ] ) in the lagrange density  ( [ lag ] ) and integrating over the whole space one gets following expression @xmath147 @xmath148
\sin\theta d\theta d\tilde r\,,\\ \label{angmom2 } i_{\omega\omega}^{(12)}(r ) & = & { { \displaystyle}\frac}{\pi}{2e^3f_\pi}{\displaystyle}\int\limits_0^\infty\int\limits_{0}^\pi \left[{{\displaystyle}\frac}{\tilde r^2}{4}{{\displaystyle}\frac}{\sin^2\theta \cos^2\theta}{\sin^2\theta}\sin^2f+ \left({{\displaystyle}\frac}{\tilde r^2}{4}+{{\displaystyle}\frac}{(3+\cos2\theta)}{2 } { { \displaystyle}\frac}{\sin^2\theta}{\sin^2\theta}\sin^2f\right)f_\theta^2 \right.{\nonumber}\\ & & \qquad\qquad\left .
+ { { \displaystyle}\frac}{\tilde r}{2}\left ( \tilde r^2 + 4{{\displaystyle}\frac}{\sin^2\theta}{\sin^2\theta}\sin^2f\right ) f_\theta f_{\tilde r}\sin(\theta-\theta ) \right.{\nonumber}\\ & & \qquad\qquad\left .
+ { { \displaystyle}\frac}{\tilde r^2}{4}\left(2\,[\,2-\cos2(\theta\!-\!\theta)+\cos2\theta\ , ] { { \displaystyle}\frac}{\sin^2\theta}{\sin^2\theta}\sin^2f+\tilde r^2\sin^2(\theta\!-\!\theta ) \right ) f_{\tilde r}^2 \right.{\nonumber}\\ & & \qquad\qquad\left . + { { \displaystyle}\frac}{1}{4}\sin^2f\left(\tilde r^2 + 2(3+\cos2\theta){{\displaystyle}\frac}{\sin^2\theta}{\sin^2\theta } \sin^2f \right)\theta_\theta^2 \right .
{ \nonumber}\\ & & \qquad\qquad\left .
-{{\displaystyle}\frac}{\tilde r^2}{2}\sin^2f[-3+\cos2(\theta-\theta)\ , ] f_{\tilde r}^2\theta_\theta^2\right ] \sin\theta d\theta d\tilde r\,,\\ \label{angmom3 } i_{\omega\omega}^{(12)}(r ) & = &   { { \displaystyle}\frac}{\pi}{2e^3f_\pi}{\displaystyle}\int\limits_0^\infty\int\limits_{0}^\pi \sin^2f\left [ \left({{\displaystyle}\frac}{f_\theta^2}{2}+{{\displaystyle}\frac}{\tilde r^2}{8 } ( 1 + 4f_{\tilde r}^2)\right){{\displaystyle}\frac}{\cos\theta}{\sin\theta}\sin2\theta \right .
{ \nonumber}\\ & & \qquad\qquad\left .
+ \left({{\displaystyle}\frac}{\tilde r^2}{4}+ { { \displaystyle}\frac}{\sin^2\theta}{\sin^2\theta}\sin^2f- \tilde rf_{\tilde r } [ f_\theta\sin(\theta-\theta)-\tilde rf_{\tilde r}]\right)\theta_\theta \right .
{ \nonumber}\\ & & \qquad\qquad\left .
+ { { \displaystyle}\frac}{\cos\theta}{2\sin\theta}\theta_\theta^2\sin^2f\sin2\theta\right ] \sin\theta d\theta d\tilde r\,,\\ i_{\omega\omega}^{(33)}(r ) & = &   { { \displaystyle}\frac}{\pi}{2e^3f_\pi}{\displaystyle}\int\limits_0^\infty\int\limits_{0}^\pi \sin^2f\left ( { { \displaystyle}\frac}{\tilde r^2}{2}+ 2f_\theta^2 + 2\tilde r^2f_{\tilde r}^2 + 2\theta_\theta^2\sin^2f\right)\sin^2\theta \sin\theta d\theta d\tilde r\ , , \label{angmom4}\end{aligned}\ ] ] where @xmath149 is static mass functional ( eq .  ( [ masstat ] ) ) and @xmath150 , @xmath151 , @xmath152 , @xmath153 are four moments of inertia of the skyrmion . in deriving equations  ( [ angmom])-([angmom4 ] )
we have made use of the axial symmetry of the system .
defining now canonical conjugate variables in the body - fixed reference system @xmath154 one obtains the hamiltonian of the system @xmath155}+{{\displaystyle}\frac}{\hat t_3 ^ 2}{4i_{\omega\omega}^{(33 ) } } \nonumber\\ & = & m+{{\displaystyle}\frac}{\hat { \bf t}^2i_{\omega\omega}^{(12 ) } + \hat { \bf j}^2i_{\omega\omega}^{(12 ) } + 2\hat { \bf t } \cdot \hat{\bf
j}i_{\omega\omega}^{(12 ) } -\left(i_{\omega\omega}^{(12 ) } + i_{\omega\omega}^{(12 ) } -2i_{\omega\omega}^{(12)}\right)\hat t_3 ^ 2 } { 4\left[i_{\omega\omega}^{(12)}i_{\omega\omega}^{(12)}- \left(i_{\omega\omega}^{(12)}\right)^2\right]}+{{\displaystyle}\frac}{\hat t_3 ^ 2}{4i_{\omega\omega}^{(33)}}\ ,
. \label{hamiltonian}\end{aligned}\ ] ] here @xmath156 is the isospin operator , @xmath157 the spin operator . if in addition one considers states with @xmath158 , where @xmath159 is the grand spin , and sandwiches the hamiltonian between wave functions @xmath160 in the body - fixed system of reference , the energy of the system takes the form @xmath161
\left(t(t+1)-t_3 ^ 2\right ) } { 4\left[i_{\omega\omega}^{(12)}i_{\omega\omega}^{(12)}- \left(i_{\omega\omega}^{(12)}\right)^2\right ] } + { { \displaystyle}\frac}{t_3 ^ 2}{4i_{\omega\omega}^{(33)}}\ , , \end{array } \label{energy}\ ] ] where we used the constraint @xmath162 which is due to the axial symmetry of ansatz  ( [ ansatzlast ] ) .
the degeneracy on the quantum number @xmath163 is partially lifted , i.e.  the energy of the skyrmion only coincides for states with @xmath164 . for example , the following pairs of baryons have the same masses : @xmath165 but the baryons in the states @xmath166 and @xmath167 do _ not _ have the same masses .
the above specified quantum numbers refer to the body - fixed system of reference .
in contrast to a deformed skyrmion in an isotropic background ( e.g.  the vacuum or isotropic and isosymmetric nuclear matter ) @xcite , our present case corresponds to a deformed nucleon or delta in an unisotropic background given by the nuclear density profile .
therefore , the projection onto good quantum numbers in the laboratory system of reference does not follow here from a simple rotation of the intrinsically deformed object .
in fact , only the third components of angular momentum and isospin in the body - fixed system , @xmath163 and @xmath168 , are good quantum numbers .    as a check ,
let us consider whether we can obtain the canonical spherically symmetric solutions from our deformed ones .
first , it is easy to see that @xmath169 where @xmath170 is the moment of inertia of the spherically symmetric hedgehog ansatz , @xmath171 consequently @xmath172 and one obtains the well known energy expression for a baryon in the spherically symmetric case @xmath173 more generally , the following relation for all observables can be shown to hold : @xmath174 this shows that our construction has the proper limit when all deformation parameters are set to zero .
predictions for the masses of the quantized solitons are presented in fig .
[ nmass ] , where we show the nucleon mass normalized to its free space value .
note that the distance @xmath131 between the centers of the nucleon and the nucleus is normalized to the corresponding radius @xmath175 of the nucleus for given baryon number @xmath105 .
we show the behavior of the nucleon mass only for one representative nucleus of the three classes , i.e. for @xmath94o , @xmath96fe and @xmath97au .
the other nuclei from the same class exhibit a similar behavior .
for the radii of the nuclei we use formula @xmath176 the behavior of the nucleon mass in the nucleus is similar for all light , medium  heavy and heavy nuclei
. when the density decreases , the mass of the nucleon increases approaching its free space value at the edge of the nucleus .
we conclude that the mean effective mass of the nucleon in real nuclei will not decrease as much as  @xmath2 in contrast to the case of homogeneous nuclear matter with normal density @xmath177 .
the finite nuclei approach for the medium essentially improves the results , especially for light nuclei , i.e. taking into account realistic mass distributions one obtains more reasonable values for the mass renormalization . for example , this renormalization is about  @xmath178 for a nucleon placed in the center of @xmath93c or @xmath94o . in principle , one could now calculate an average medium modification in a given nucleus by simply folding the calculated effective mass with the nuclear density .
we do not show the results of such calculations here because the calculation could be further refined .
if one uses e.g. the shell structure of the nuclei , then all nucleons will be situated at some distances ( within the shells ) from their mass center and the mean effective value of the nucleon mass in such systems will be further reduced for nucleons in heavy nuclei .
further improvements of the results can be made by including a dilaton field  @xcite .
we also note that a similar renormalization of the nucleon mass was found in ref.@xcite .
as it was noted in the previous section , the mass of the @xmath0 states depends on the absolute value of the third component of isospin .
we found small deviation from degeneracy for delta states with @xmath179 and @xmath180 .
the difference of the masses @xmath181 of @xmath0 in the states @xmath166 and @xmath167 ( see eq .
( [ states ] ) ) has a negative value near the center of the nucleus and after reaching its minimum begins to increase towards positive values reaching at some distance its maximal value and then smoothly drops to zero .
for example , it has its minimal value @xmath182mev at @xmath183  fm and maximal value @xmath184mev at @xmath185  fm for deltas in @xmath93c . for the case of @xmath96fe ,
the corresponding values are @xmath186(@xmath187  @xmath188mev and @xmath186(@xmath189  @xmath190mev .
we remind the reader that all these numbers refer to states defined in the body - fixed system of reference as discussed before .    for the relative mass difference @xmath191 we obtained a behavior similar to the one of nucleon mass ( cf .
fig.([nmass ] ) ) and do not graphically present it here . the difference @xmath191 decreases when the density increases and has its minimal value at the center of the nucleus .
more specifically , @xmath192 has the value 0.867 for @xmath93c , 0.615 for @xmath96fe and 0.609 for @xmath97au , in the respective center of the nucleus .
note that in free space @xmath193mev in the current approach .
as to the ratio of @xmath0 and @xmath194 masses , @xmath195 , it is almost constant inside the nucleus and approximately equal to that ratio in free space .
the calculation of the electromagnetic currents reduces to the calculation of the isoscalar and the isovector currents . in topological models and in particular in the skyrme model ,
the isoscalar current is equal to the model  independent baryon number current ( for a discussion on the role of vector mesons , see ref.@xcite ) @xmath196 with @xmath197 the totally antisymmetric tensor in four dimensions . similarly , the isovector current , which is the noether current associated to the global transformations @xmath198 ( @xmath199 are @xmath200 matrixes ) , has the form in the equation  ( [ veccur ] ) while there is one over the index @xmath201 . ]
@xmath202 + [ r_{\nu},[r_{\mu},r_{\nu}]]\}\ , , \label{veccur}\ ] ] where @xmath203 then , the operator of the electric charge of the nucleon is determined by the zeroth components of the isoscalar and the isovector currents , @xmath204 where @xmath205 is
the third component of the isovector current and @xmath206 is the third component of the isotopic spin operator . sandwiching this current between corresponding nucleon states , we obtain the proton and neutron charge distributions .
the normalized isoscalar root - mean - square radius along each axis , @xmath207 , ( normalized with the respect to the isoscalar charge ) can be determined from the zeroth component of the baryon current and has the form @xmath208 where @xmath209 ( @xmath210 ) indicates the projection of the radius vector on the corresponding axis .
similarly , the isovector root - mean - square radii ( normalized to the isovector charge ) , @xmath211 , are determined by the zeroth component of the vector current @xmath205 @xmath212 where @xmath213\,.\ ] ] one can see that the influence of the medium enters indirectly via the modified profile functions .
results for the in  medium radii are collected in table  [ table ] .
we give the isoscalar @xmath214 and isovector @xmath215 root - mean - square radii of the nucleons along the @xmath216 and @xmath128 axes inside the nuclei @xmath94o , @xmath96fe and @xmath97au , evaluated in the center of each nucleus and at those distances where the intrinsic quadrupole moment of the nucleons has its first and second extremum , correspondingly ( see the discussion below ) .
we note that due to the axial symmetry @xmath217 .
clearly , these radii swell , as first pointed out in ref.@xcite in the analysis of data on quasi - elastic electron - nucleus scattering . however , this swelling is not isotropic and depends very sensitively on the distance from the center of the nucleus .
the concept of a common swelling by one scale factor is certainly much too simple .
our results are therefore not in contradiction with the y - scaling analysis of ref.@xcite ( see also the discussion in @xcite on this issue ) .
as noted before due to the deformation the nucleon can have a nonvanishing _ intrinsic _ quadrupole moment  @xcite @xmath218 that is characterized by the third component @xmath163 of the angular momentum in the body - fixed frame of reference .
the corresponding isoscalar and isovector intrinsic quadrupole moments follow from ( [ isosr ] ) and ( [ isovr ] ) by appropriate projection @xmath219 consequently , the proton and neutron intrinsic quadrupole moments follow from @xmath220 the dependence of the isoscalar and isovector intrinsic quadrupole moments as a function of @xmath221 is shown in fig .
again we present results only for the nuclei @xmath94o , @xmath96fe and @xmath97ag . for the other nuclei the dependences are similar .
one can see that intrinsic quadrupole moments are small .
this smallness of the intrinsic quadrupole moments means that nucleons in the nucleus are weakly deformed .
the intrinsic quadrupole moments has two extrema .
near the center of nucleus they have negative values and a corresponding minimum at that distances where density changes are @xmath12420% for @xmath94o and @xmath1247% for @xmath96fe and @xmath97au , respectively .
then with increasing density the intrinsic quadrupole moments become positive and have a maximum at the distance where the density of the nucleus has fallen by @xmath222% .
note also that the isoscalar intrinsic quadrupole moment is about two times smaller than the isovector one . consequently near the center of the nucleus protons
are flattened along the axis @xmath128 ( oblate ) and and near the border of the nucleus they are extended ( prolate ) .
in contrast neutrons near the center of nucleus slightly prolate and near the border of nucleus oblate . consequently , at some intermediate distance the nucleons must have spherical form .
finally , let us recall that the observable static quadrupole moment in the laboratory system of an intrinsically deformed baryon in an isotropic background is linked to its intrinsic quadrupole moment by the wigner - eckhardt theorem   as follows  @xcite @xmath223 where the relations are valid for maximally projected states with @xmath224 .
thus for a deformed baryon of total angular momentum @xmath225 the quadrupole moment in the _ laboratory _ system is zero for its ground state , i.e.  for the nucleon .
it would be interesting to investigate the transition quadrupole moments to excited states , extending the work of @xcite . as mentioned before
, our case corresponds to a deformed baryon in an unisotropic background .
a simple projection to the laboratory system is therefore not possible .
in summary , we have considered properties of a nucleon embedded in finite nuclei making use of the skyrme model .
deformation effects are taken into account by the distortion of chiral profile function under the action of the nuclear medium , parameterized in terms of the functionals eq.([funcrho ] ) .
nucleons and deltas emerge as quantized solitons after performing an adiabatic rotation in co - ordinate and isospace which takes into account the axial symmetry of the system ( leading to two different moments of inertia , @xmath226 .
the main results of this study can be summarized as follows :    * the density dependence of the nucleon mass shows a more realistic behavior than in the case of a uniform density as e.g. for homogeneous nuclear matter .
the effective nucleon mass has its minimum at the center of the nucleus and approaches its free space value at the surface of the nucleus . *
the nucleons in finite nuclei are weakly deformed .
they acquire a small intrinsic quadrupole moment which , however , is strongly dependent on the distance from the center of the nucleus , e.g. the proton ( neutron ) deformation changes from a oblate ( prolate ) to a prolate ( oblate ) form as it is moved toward the surface of the nucleus . *
similarly , there is a direction - dependent swelling of the isoscalar and isovector root - mean - square radii .
we have stressed that the concept of a uniform swelling factor is too simple a concept to apply to real nuclei . * as consequence of the axial symmetry of the system , the @xmath227 @xmath0 states ( @xmath228 ) and the @xmath229 ones ( @xmath230 ) have slightly different masses in finite nuclei .
let us briefly give an outlook for further studies based on these results . as it was mentioned above this approach can be used for model investigations of the shell structure of nuclei and to study the properties of a remote nucleon in such a system ( halo nuclei ) . in the latter case the nucleon in the finite nucleus
should not be considered as one within the nucleus but weakly coupled to the nucleus at large distances .
also , the model naturally allows to study nucleon properties in the lightest nuclei , especially the electromagnetic structure of nucleons in @xmath231he , which is of current experimental interest .
work along these lines is under way .
the work of u. yakhshiev has been supported by intas ys fellowship @xmath23200 - 51 .
we thank andy jackson for a useful correspondence .    99 a.m. rakhimov , m.m .
musakhanov , f.c . khanna and u.t .
yakhshiev , phys . rev . *
c58 * , 1738 ( 1998 ) .
a.m. rakhimov , f.c .
khanna , u.t .
yakhshiev , m.m .
musakhanov , nucl . phys . *
a643 * , 383 ( 1998 ) ; + m. musakhanov , a. rakhimov , u. yakhshiev , z. kanokov , phys . atom .
j. nucl . phys . )
* 62 * , 1988 ( 1999 ) .
meiner , phys .
lett . * 62 * , 1013 ( 1989 ) ; phys . lett . *
b220 * , ( 1989 ) ; nucl .
* a503 * , 801 ( 1989 ) . i. klebanov , nucl . phys . *
b262 * , 133 ( 1985 ) ; + n.k .
glendenning , phys . rev . *
c34 * , 1072 ( 1986 ) ; + e. wst , g.e .
brown and a.d .
jackson , nucl .
phys . * a 468 * , 450 ( 1987 ) ; + a.s . goldhaber and n.s .
manton , phys .
b198 * , 231 ( 1987 ) ; + t.s .
walhout , nucl . phys . *
a484 * , 397 ( 1988 ) ; nucl . phys . * a 519 * , 816 ( 1990 ) ; + a.d .
jackson , a. wirzba and l.s .
castillejo , phys . lett . *
b198 * , 315 ( 1987 ) ; nucl . phys . *
a486 * , 634 ( 1988 ) ; + m. kugler and s. shtrikman , phys . lett . *
b208 * , 491 ( 1988 ) ; phys . rev . *
d40 * , 3421 ( 1989 ) ; + l.s .
et al . _ ,
a501 * , 801 ( 1989 ) ; + h. forkel _
et al . _ ,
. phys . * a504 * , 818 ( 1989 ) ; + a.d .
jackson , c. weiss and a. wirzba , nucl .
a529 * , 741 ( 1991 ) ; + a. wirzba , in `` baryons as skyrme solitons '' , g. holzwarth ( ed . ) , world scientific , singapore , 1993 ; + w.k .
baskerville , nucl . phys . *
a596 * , 611 ( 1996 ) ; + p. amore , nuovo cim . * a111 * , 493 ( 1998 ) .
v. vento , g. baym and a.d .
jackson , phys .
b102 * , 97 ( 1981 ) .
u.b . kaulfu and ulf - g .
meiner , phys . rev . *
d31 * , 3024 ( 1985 ) .
a. rakhimov , t. okazaki , m.m .
musakhanov , f.c .
khanna , phys .
b378 * , 12 ( 1996 ) . c. hajduk , b. schwesinger , phys
. lett . * b140 * , 172 ( 1984 ) ; + c. hajduk , b. schwesinger , nucl . phys . *
a453 * , 620 ( 1986 ) .
kopeliovich , b.e .
stern , sov .
j. jetp lett .
* 45 * , 165 ( 1987 ) ; + v.b .
kopeliovich , sov .
j. nucl.phys .
* 47 * , 1495 ( 1988 ) ; + b.e .
stern , sov .
j. nucl.phys .
* 49 * , 530 ( 1989 ) ; + v.b .
kopeliovich , russ .
j. nucl.phys .
* 56 * , 160 ( 1993 ) .
nikolaev , sov .
* 20 * , 401 ( 1989 ) ; + v.a .
nikolaev , o.g .
tkachev , sov .
* 21 * , 1499 ( 1990 ) ; + r.m .
nikolaeva , v.a .
nikolaev , o.g .
tkachev , sov .
* 23 * , 542 ( 1992 ) .
t. kurihara , h. kanada , t. otofuji , s. saito , prog .
. phys . * 81 * , 858 ( 1989 ) .
bogomolniy , v.a .
fateev , sov .
. phys . * 37 * , 228 ( 1983 ) .
yakhshiev , n.a.taylanov , uzb . j. phys . * 2 * , 114 ( 2000 ) .
e. braaten and l. carson , phys .
* 56 * , 1897 ( 1986 ) ; + e. braaten and l. carson , phys .
rev . * d38
* , 3525 ( 1988 ) ; + e. braaten , s. townsend and l. carson , phys . lett . *
b235 * , 147 ( 1990 ) ; + n.s .
manton , phys .
b192 * , 177 ( 1987 ) ; + m.f .
atiyah and n.s .
manton , commun .
phys . * 153 * , 391 ( 1993 ) ; + r.a .
leese and n.s .
manton , nucl .
a572 * , 575 ( 1994 ) ; + c.j .
houghton , n.s .
manton and p.m. sutcliffe , nucl .
* b510 * , 507 ( 1998 ) .
i. zahed _ et al .
_ , phys . rev . *
d31 * , 1114 ( 1985 ) ; + c .- k .
lin , nucl .
* a449 * , 673 ( 1986 ) ; + j. wambach , h.w .
wyld and h.m .
sommermann , phys . lett . * b186 * , 272 ( 1987 ) ; + t. otofuji _
. phys . * 78 * , 527 ( 1987 ) ; + i. mishustin , sov .
jetp * 71 * , 21 ( 1990 ) ; + f. leblond and l. marleau , phys . rev . *
d58 * , 054002 ( 1998 ) .
a. wirzba and v. thorsson , nucl .
* a589 * , 633 ( 1985 ) . t. ericson and w. weise , _ pions and nuclei _
( claredon , oxford , 1988 ) .
a. akhiezer , a.g .
sitenko , v.k .
tartakovski , _ nuclear electrodynamics _ ( springer verlag , 1994 ) .
adkins , c.r .
nappi and e. witten , nucl . phys . *
b228 * , 552 ( 1983 ) ; + g.s .
adkins and c.r .
nappi , _ ibid . _ * b233 * , 109 ( 1984 ) .
j. schechter , phys rev .
* d21 * , 3393 ( 1980 ) ; h. gomm _ et al .
_ , phys . rev . *
d33 * , 3476 ( 1986 ) .
meissner , a. rakhimov and u.t .
yakhshiev , phys . lett . *
b473 * , 200 ( 2000 ) .
meiner , phys . rep .
* 161 * , 213 ( 1988 ) .
noble , phys .
* 46 * , 412 ( 1981 ) . i. sick , phys
b157 * , 13 ( 1985 ) .
a. bohr and b.r .
mottelson , _ nuclear structure , vol .
ii , nuclear deformations _ ( w.a . benjamin , reading , massachusetts , 1975 ) .    .
[ table ] the isoscalar @xmath214 and isovector @xmath215 nucleon root - mean - square radii along the @xmath216 and @xmath128 axes inside various nuclei .
the values of these quantities are given for the center of the nuclei and at those distances where the nucleon intrinsic quadrupole moments have their first and second extremum values , correspondingly .
@xmath131 is the distance between centers of the nucleus and skyrmion topological center and the last column is the ratio of the density of nucleus @xmath233 to its value at the center @xmath234 .
we note that @xmath235 . for comparison
we also give the free space values of the root - mean - square radii . [ cols="^,^,^,^,^,^,^",options="header " , ]