the determination of the binding energy of nuclei starting from realistic @xmath1 potentials is one of the subjects which has received permanent attention from the early days when the brueckner - bethe - goldstone ( bbg ) equation introduced methods to overcome the strong repulsion of the nuclear forces at short distances .    at present , several many body techniques compete to accurately determine the binding energy of nuclear matter starting from the realistic nn potentials .
one of them is the correlated basis functions ( cbf ) @xcite , which follows the line of the hypernetted chain approach ( fhnc ) @xcite .
another one follows the traditional bbg approach @xcite , and , although costly numerically , the variational monte carlo method ( mc ) has allowed to make , in principle , exact calculations , although limited to nuclei with small or medium value of a @xcite .
methods like the selfconsistent treatment of the nucleon selfenergy have also introduced new advances in the field @xcite .
the need for three body forces has also been emphasized and the present status is that it is difficult to be quantitative on the strength of these forces , and usually they are parameterized in order to adjust the precise value of the binding energy @xcite .
it has also been noted in @xcite that short range correlations play an important role when considering these three body forces .
a common feature of these approaches is that they start from the realistic nucleon nucleon interaction , obtained from fits to nn data and deuteron data .
they use hence the free nn interaction as input .
one of the important ingredients of this interaction is the one pion exchange ( ope ) .
however , from detailed studies of the pion nuclear interaction it is well known that the pion properties in the nuclear medium are sizably renormalized @xcite .
there is also the question of the intermediate range attraction , which is basic in the binding of nuclei .
models for this interaction would contain @xmath0 exchange , uncorrelated two pion exchange and omega exchange @xcite . in as much as the pion properties are changed in the medium , so should the two pion exchange be modified .
medium effects in the two - pion exchange have been investigated in early works like in ref .
@xcite restricting themselves in this case to a subset of two pion exchange diagrams with no @xmath3-isobar intermediate states , by including pauli blocking in the intermediate nucleons
.    the medium modification of the two pion exchange got a new impulse after the models of the @xmath4 interaction in the medium showed large modifications @xcite , later on softened by the introduction of chiral constraints  @xcite .
the implementation of this medium modified @xmath4 interaction in the correlated two pion exchange @xmath1 potential increased appreciably the @xmath1 attraction in nuclear matter .
this was partially reduced by the consideration of the chiral constraints in refs .
@xcite . in these latter references
the importance of short range correlations which modify the @xmath5 nucleus selfenergy was already discussed .
the further use of medium modified vector meson masses led to improvements in the nuclear matter saturation curve .
a new perspective into this problem has been made possible by studies of meson meson interaction within chiral unitary approaches  @xcite which allow to improve the description of the correlated two pion exchange @xmath1 interaction  @xcite , as an alternative to the conventional @xmath0-meson exchange interaction .
this picture of the @xmath0 exchange was mandatory after extensive studies showing that the @xmath0 is not a genuine resonance , made up of @xmath6 but just the manifestation of a pion - pion resonance state created by the interaction of the pions , what is called a dynamically generated resonance .
this shows up naturally within the context of chiral unitary approaches which use the input of the chiral lagrangians for the meson meson interaction and extends chiral perturbation theory to implement exactly unitarity in coupled channels @xcite .
this means the @xmath0 exchange inside a nuclear medium will also be modified as direct consequence of the change of the pion properties .
our aim in the present paper is to start from this new picture for the @xmath0 exchange , use also the standard approach for the uncorrelated two pion exchange and modify these in the nuclear medium to see what changes one finds from these sources .
further improvements come from the consideration of short range correlations not only in the pion selfenergy but also in the vertex functions appearing in the model .
the changes obtained are moderate , thanks to the simultaneous consideration of these nuclear short range effects in the calculation . in the absence of these
, the renormalization of the @xmath1 interaction is huge . yet , even the moderate results obtained are large enough to motivate further calculations of the nuclear binding and other properties of matter .    the paper is organized as follows . in sec .
ii we provide those elements of the chiral lagrangian which are relevant for the present calculations and briefly discuss peculiarities of the finite baryon density . in sec .
iii we consider the modification of the one pion exchange @xmath1 force in the nuclear medium and in sec .
iv we discuss the propagation of two pions in the nuclear matter .
v and  vi are devoted to the in - medium two pion exchange in the scalar - isoscalar channel , both , correlated and uncorrelated .
the technical details are relegated to the appendix .
in this section we will briefly specify those elements of the effective chiral lagrangian in the meson - baryon sector which are relevant for the subsequent calculations .
the effective chiral lagrangian is written as the sum of a purely mesonic lagrangian @xmath7 and the baryonic lagrangian @xmath8 @xmath9 both are organized in a derivative and quark mass expansion .
the lowest order mesonic lagrangian @xmath10 is given by @xmath11 and contains the most general low - energy interactions of the pseudo - scalar meson octet .
( [ lpipi ] ) the symbol @xmath12 indicates the trace in flavor space , the goldstone fields are collected in a unitary matrix @xmath13 , @xmath14  mev is the pseudoscalar decay constant and the leading symmetry - breaking term @xmath15 is linear in the quark masses . for @xmath16 and in the isospin limit @xmath17 .
the lowest order baryon octet - meson octet lagrangian reads @xmath18 \rangle\ ] ] where the brackets @xmath19 $ ] and @xmath20 denote commutators and anti - commutators , respectively .
the covariant derivative of the @xmath21 baryon matrix @xmath22 is defined as @xmath23\ ] ] in the absence of external field eqs .
( [ barlagr ] ) and ( [ covarderiv ] ) involve other quantities @xmath24 the @xmath21 axial vector coupling constants are determined by neutron and hyperon @xmath25-decay .
one finds @xmath26 , @xmath27 and the axial coupling constant is @xmath28 . in the @xmath16
limit the lagrangian simplifies to @xmath29 where @xmath30 is a two component dirac field @xmath31 .
the @xmath32 vertex ( a ) , the contact @xmath33 ( b ) and pion pole ( c ) terms . ]       the opep diagrams .
the first graph is the vacuum contribution .
the second and third diagrams of the upper and first two diagrams of the lower sets of diagrams correspond to the @xmath44 and @xmath45 rpa series .
the last two diagrams account for the excitation of @xmath46 states and the contribution of the @xmath47-wave optical potential , respectively .
the latter diagrams play a minor role and will not be further discussed . ]
we start with the one - pion exchange @xmath1 potential ( opep ) .
the typical diagrams modifying it are shown in fig .
[ opepmedimgraphs ] where the propagation of exchanged pions is distorted by interactions with nucleons forming the fermi sea .
the fermionic bubbles describe the decay of the pion in @xmath44 and @xmath45 states and take into account the conventional nuclear matter polarization effects .
all these diagrams are responsible for the interaction of two nucleons in the particle - particle ladder channel with the in - medium virtual excitations . in fig .
[ opepmedimgraphs ] the first graph is the vacuum contribution .
the second and third diagrams correspond to the @xmath44 rpa series and the last diagrams accounts for the excitation of @xmath46 states and the contribution of the @xmath47-wave optical potential , respectively .
the @xmath48-wave pion self energy is given by @xmath49^{-1}\ ] ] where @xmath50 is the landau - migdal parameter  @xcite , @xmath51 is the nuclear matter density and @xmath52 is the lindhard function accounting for the direct and crossed contributions of @xmath44 and @xmath45 excitations with the normalization of the appendix of ref .
@xcite .     left panel : the opep in the vacuum ( solid curve ) and at normal @xmath53  @xmath54 nuclear matter density ( dashed curve ) .
right panel : opep at several densities .
, title="fig : " ]   left panel : the opep in the vacuum ( solid curve ) and at normal @xmath53  @xmath54 nuclear matter density ( dashed curve ) .
right panel : opep at several densities .
, title="fig : " ]    the opep in the momentum space takes the form @xmath55 where we have defined @xmath56 where @xmath57 is the pion propagator in the medium @xmath58 and @xmath59 stands for a monopole form factor @xmath60 with the cut off scale @xmath61  gev , @xmath62 is a momentum transfer and @xmath63 .
note , that opep depends on the real part of the polarization operator only , since for @xmath64 one has @xmath65 . the well - known vacuum @xmath1 amplitude is recovered in eq .
( [ wopep ] ) by setting @xmath66 or in the limit @xmath67 .
our results for @xmath68 are presented in fig .
[ opepmedium ] ( left ) .
the standard vacuum behavior is shown by the solid curve .
the dashed curve represents the modified opep at normal nuclear matter density were we observe an additional strong increase of the strength associated with the attractive excitation of @xmath44 and @xmath69 collective states , with @xmath44 playing a dominant role .
their individual contributions are shown by dot - dot - dashed and dot - dashed curves , respectively .
the opep with short range correlations ( solid curves ) as a function of the nuclear matter density .
the left and right panels are the central and tensor parts , respectively.,title="fig : " ]   the opep with short range correlations ( solid curves ) as a function of the nuclear matter density .
the left and right panels are the central and tensor parts , respectively.,title="fig : " ]    the analytic properties of the in - medium pion propagator which we use here may be verified by using the dispersion representation for the green function @xmath70 in this case the opep of eq .
( [ wopep ] ) can be written in terms of the absorptive part of the pion propagator only @xmath71\end{aligned}\ ] ] the causality requires that both equations must produce the same result .
indeed as one can see in fig .
[ opepmedium ] the curves calculated with dispersive and absorptive parts of the pion propagator are practically indistinguishable .
we would like to note that a strong modification of the opep at finite baryonic density observed here is not new and was predicted long time ago by migdal  @xcite .
there it was also shown that in - medium modified opep helps to explain the unnatural parity states in finite nuclei , for instance , the shift of @xmath72 state in closed shell nuclei .    in the right panel of fig .
[ opepmedium ] we show our combined plot for a few densities @xmath73 and @xmath74 . here
, we also show the sensitivity of our results to the value of the landau - migdal parameter @xmath75 .
we find that the increase of @xmath75 from 0.6 to 0.7 makes the opep softer .
this fact suggests that the proper treatment of the @xmath1 short range correlations is important for understanding the in - medium properties of the opep .    at this point
we would like to mention that in a realistic calculation one will have to add strong repulsive forces at short distances .
this can be done in a straightforward way using any of many body schemes discussed in the introduction .
the correlations of this part of the interaction would effectively modulate the @xmath5 exchange interaction , introducing the correlation parameter @xmath75  @xcite .
the denominator in eq .
( [ inmediumgf ] ) takes into account this effect between @xmath48-wave bubbles in the diagrams of fig .
[ opepmedimgraphs ] , but not between the external nucleon and the contiguous bubble . to account for this
we make the separation between the longitudinal and transverse parts of the pion effective interaction  @xcite @xmath76 where @xmath77 is the cartesian component of the unit vector @xmath63 and @xmath78f({\mbox{\boldmath $ q$}})^2\ ] ] @xmath79 when we perform the sum of diagrams in fig .
[ opepmedimgraphs ] then we get @xmath80\sigma_1^i \sigma^j_2    { \mbox{\boldmath $ \tau$}}_1 { \mbox{\boldmath $ \tau$}}_2   \nonumber \\ & = & w_0(q,\rho ) { \mbox{\boldmath $ \sigma$}}_1 { \mbox{\boldmath $ \sigma$}}_2   { \mbox{\boldmath $ \tau$}}_1 { \mbox{\boldmath $ \tau$}}_2
+ w_2(q,\rho ) ( { \mbox{\boldmath $ \sigma$}}_1 \hat{{\mbox{\boldmath $ q$ } } }   { \mbox{\boldmath $ \sigma$}}_2 \hat{{\mbox{\boldmath $ q$ } } } -   \frac{1}{3}{\mbox{\boldmath $ \sigma$}}_1 { \mbox{\boldmath $ \sigma$}}_2 )   { \mbox{\boldmath $ \tau$}}_1 { \mbox{\boldmath $ \tau$}}_2 \end{aligned}\ ] ] @xmath81 \\ w_2(q,\rho)&= & \frac{\mathcal{v}_l(q)}{1-\mathcal{u}(q,\rho)\mathcal{v}_l(q ) } -   \frac{\mathcal{v}_t(q)}{1-\mathcal{u}(q,\rho)\mathcal{v}_t(q)}\end{aligned}\ ] ] where we have explicitely separated the central and the tensor parts of the interaction . in fig .  [ opepmediummodif ]
we show the results for the central @xmath82 ( left panel ) and tensor @xmath83 ( right panel ) parts ( omitting the spin - isospin operators ) as a function of the baryonic density .
note that in the limit @xmath84 , @xmath85 where @xmath86 is given by eq .
[ wopep ] .
in this section we turn to the dynamics of the two pion system in the nuclear medium .
our main interest here is the propagation of two @xmath47-wave pion pairs .
as it is well known in the @xmath47-wave scattering the use of the proper unitarization schemes lead to the generation of the @xmath0-meson . later on
we will use this result for the in - medium @xmath1 interaction mediated by exchange of two correlated pions in the scalar - isoscalar @xmath0-meson channel .
the unitary series representing @xmath4 scattering amplitude . ]    for the @xmath87 scattering process , defined by the cartesian isospin indices @xmath88 , the use of the standard @xmath15pt procedure in expanding the @xmath89 of eq .
( [ lpipi ] ) to order @xmath90 results in the tree level contact interaction @xmath91 where @xmath92 and @xmath93 is the off - shell part of the invariant @xmath4 amplitude . at this order of the pion field expansion the isoscalar @xmath47-wave @xmath4 partial amplitude ( @xmath94 ) is obtained from the standard decomposition @xmath95 where @xmath96 are the legendre polynomials and @xmath97 accounts for the statistical factor occurring in states with identical particles : @xmath98 for @xmath43 in the unitary normalization of the states  @xcite .
the tree level scalar - isoscalar @xmath4 scattering amplitude is @xmath99 in eq .
( [ vpipi00 ] ) the off shell part depends on choice of @xmath13 and is equal to zero for on mass shell pions .    following ref .
@xcite and using the bethe - salpeter equation we unitarize the @xmath47-wave @xmath4 scattering amplitude ( see fig .  [ unitseries ] ) @xmath100^{-1},\ ] ] where @xmath101 is a scalar two - pion loop function @xmath102   ( k^2-m_{\pi}^2+i0^+)}\ ] ] where @xmath103 and @xmath104 .
the @xmath105 function is analytic with a cut along the positive real axis starting at the @xmath4 threshold .
note that , eq .  ( [ tpipiunitary ] ) contains a pole in the second riemann sheet corresponding to the @xmath0-meson with mass and width @xmath106  mev .
renormalization of the in - medium @xmath4 amplitude including @xmath48-hole and @xmath3-hole excitations . ]    in the nuclear medium the @xmath47-wave @xmath4 scattering amplitude and therefore the @xmath0-meson get renormalized and explicit calculations were done in refs .
the diagrams at one loop level are shown in fig .
[ renormpipiinmedium ] .
for instance , the amplitudes corresponding to the insertion of fermion bubbles in the upper meson line are given by @xmath107 ^ 2 \\
v_{\mbox{\tiny ii}}^{up}(s ) & = &   \frac{-i2}{3f_{\pi}^2 } \left ( \frac{d+f}{2 f_{\pi } } \right)^2 \int \frac{d^4 k}{(2 \pi)^4 } { \mbox{\boldmath $ k$}}^2 \mathcal{u}(k ) d_{\pi}(k ) d_{\pi}(p - k ) \big [ v_{on}(s ) + \frac{1}{3 f_{\pi}^2 } \sum_i \lambda_i \big ] \\ v _ { \mbox{\tiny iii}}^{up}(s ) & = & \frac{i}{9f_{\pi}^4 } \left ( \frac{d+f}{2 f_{\pi } } \right)^2 \int \frac{d^4 k}{(2 \pi)^4 }   { \mbox{\boldmath $ k$}}^2 \mathcal{u}(k ) d_{\pi}(p - k)\end{aligned}\ ] ] as was shown in refs .
@xcite , in the center of mass frame of the two pions the off shell part ( @xmath108 terms ) in @xmath109 cancels exactly the @xmath110 and @xmath111 terms . and one is left only with the diagrams of type ( i ) but with the on shell @xmath4 amplitude @xmath112 it is straightforward to iterate the @xmath44 and @xmath45 excitations in fig .
[ renormpipiinmedium ] ( i ) and the loop function , @xmath113 , is given by @xmath114 where @xmath115 is in - medium modified scalar loop integral @xmath116 using the spectral representation for the in - medium pion propagators , eq .
( [ inmediumgf ] ) , we get for the loop function @xmath117     \int   \limits_{x}^{\infty } d y \big[\mbox{im } \tilde{d}_{\pi}(y - x,{\mbox{\boldmath $ p$}}-{\mbox{\boldmath $ k$}})\big ] \nonumber \\ \times
\frac{y}{(p_0+y)(p_0-y+i0^+ ) } \end{aligned}\ ] ]    we refer to ref .
@xcite where different aspects of @xmath47-wave @xmath4 scattering in the nuclear medium are discussed and also the behavior of the @xmath0-meson mass and width at finite baryonic density is addressed . but here , we would like to illustrate the impact of the nuclear medium on the @xmath4 system .
for that consider the imaginary part of the loop function in the @xmath4 center of mass frame @xmath118 with @xmath119 .
this situation is relevant for the in - medium @xmath4 scattering and contains the proper information about the dynamics of the pole position of @xmath0 at finite density .
our results for the imaginary part of the scalar loop function for several densities are shown in fig .
[ imgpipimedium](left panel ) . the solid curve correspond to the vacuum loop function which can be obtained from eq .
( [ gpipigen ] ) by substituting the imaginary part of in - medium pion propagators by their vacuum expressions @xmath120 which is proportional to the two - body phase space of two particles ( pions ) . as one can see in fig .
[ imgpipimedium ]  ( left ) the effect of the medium is remarkable due to the increase of available for pions phase space because of additional pion decay branches like @xmath44 , @xmath45 and @xmath46 .     the imaginary part of the in - medium scalar loop function @xmath121 in the @xmath122 cm frame ( left panel ) with the normal dressing of the pion propagators including @xmath44 , @xmath45 and @xmath46 excitations .
the right panel shows the loop function in the space like region @xmath123.,title="fig : " ]   the imaginary part of the in - medium scalar loop function @xmath121 in the @xmath122 cm frame ( left panel ) with the normal dressing of the pion propagators including @xmath44 , @xmath45 and @xmath46 excitations .
the right panel shows the loop function in the space like region @xmath123.,title="fig : " ]    the kinematics relevant for the @xmath1 force , where two nucleons interact by exchange of mesons is defined by the moving reference frame where @xmath123 . in this case
the @xmath4 interaction in the medium has to be generalized from the results in  @xcite since there @xmath124 but @xmath125 .
we repeat all the steps that led to cancellations in  @xcite and find again the same cancellations as before but with some remnant @xmath126 dependent terms vanishing in the limit @xmath127 . to evaluate these terms
we simplify the calculation assuming @xmath62 relatively small ( this is fine for momenta below the fermi momentum ) .
concretely , we assume @xmath128 where @xmath129 is the cut off in the three momentum ( of the order of @xmath130  gev in  @xcite ) that one uses to regularize the @xmath121 function .
we also assume @xmath131 which implies @xmath132 as it is the case .
and 3 ) we expand @xmath133 in terms of @xmath134 and vice - versa to relate different terms .
after all this is done we find that the corrections can be taken into account by means of the change in the lindhard function @xmath135\ ] ] the expression in brackets in eq .
( [ umodif ] ) is @xmath136 for very small @xmath137 and also @xmath136 for @xmath138 .
hence with very good approximation we can take the bracket equal to unity and thus there are no other corrections to be done to the result of  @xcite except the obvious one of changing @xmath139 .
for these value of @xmath140 the @xmath121 is only real , contrary to the case studied in  @xcite .
results for @xmath141 can be seen in fig .
[ imgpipimedium ] ( right panel ) for different nuclear densities @xmath142 , @xmath143 and @xmath74 .
we can see that the corrections are sizable particularly at small values of @xmath144 .
the correlated two pion exchange ( crtpe ) in the scalar - isoscalar channel , the equivalent to a @xmath0 exchange in meson exchange models  @xcite , or correlated two pion exchange in the dispersion relations in  @xcite was studied in  @xcite within the context of chiral lagrangians .
one starts from the diagrams of fig .
[ ctpevacuum.eps ] , where the @xmath145 scattering shows the off shell ambiguities . to avoid these ambiguities with isoscalar exchange
it was stated in ref .
@xcite that one must include the subset of diagrams of fig .  [ setcans ] to find cancellations of the off shell @xmath4 isoscalar amplitude .
this statement was rigorously verified in ref .
@xcite , where it was shown that the consideration of these subset of chiral diagrams , fig .
[ setcans ] , including the contact @xmath41 interactions , results in the cancellation of the off - shell part of the @xmath4 amplitude , and the on - shell part of the @xmath4 amplitude can be factorized out from the loop integrals .
one step forward was given in  @xcite , where iteration of the @xmath4 interaction , through the bethe - salpeter equation , was done by means of which a simple analytical expression was obtained for the correlated two pion exchange in the scalar - isoscalar channel @xmath146^{-1}\ ] ] where @xmath147 in the @xmath1 c.m .
the vertex function @xmath148 for the triangle loop with two mesons and one baryon propagator , including @xmath149 and @xmath3 intermediate states , is evaluated in  @xcite using a cut off in @xmath150 of about 1  gev .
note that , the bracket in eq .
( [ vcrtpesigma ] ) contains a pole in the @xmath140-channel corresponding to the @xmath0-meson .
this restores the relation to the @xmath0-meson exchange , which now enters the formalism as a dynamical resonance in the @xmath4 system  @xcite ( see also related discussions in refs .
@xcite ) .
diagrams representing the correlated two pion exchange with @xmath1 , @xmath151 and @xmath152 intermediate states . ]
the set of diagrams which cancel the off - mass shell part of the scalar - isoscalar correlated two - pion exchange . ]
the in - medium diagrams involving the @xmath48-wave pion dressing and vertex corrections . the diagrams with the intermediate @xmath3 states are not shown . ]
the diagrams responsible for the renormalization of crtpe in the nuclear medium are shown in fig .
[ corrtpecrmedium ] . there , in analogy to what was done for the @xmath4 interaction we include @xmath5 selfenergy corrections as well as vertex corrections .
we find that the cancellation of the off shell part of the vacuum @xmath0-exchange discussed in  @xcite is also exact at finite baryonic density but for zero momentum transfer only ( see appendix b ) .
the results of the derivation for @xmath153 can be summarized as follows :    \1 ) the expressions for the triangle vertex functions @xmath154 and @xmath155 of ref .
@xcite are obtained in the same way replacing the two free pion propagators by the renormalized ones .
\2 ) the lindhard function entering the pion selfenergy is changed to @xmath156\ ] ] to account for the corrections obtained at @xmath153 .
these corrections are negligible for small @xmath144 and for @xmath157 of the order of @xmath158 , hence moderate in all cases .
\3 ) the final expression for the potential @xmath159 is given by eq .
( [ vcrtpesigma ] ) , which accounts for the @xmath4 rescattering , by substituting @xmath160 by @xmath161 of eq .
( [ gpipigen ] ) and taking the expression for the in - medium vertex function @xmath162 where @xmath163 \nonumber   \\
\nonumber   \\ & & \times \int \limits_0^{\infty } dy \frac{x+y+e({\mbox{\boldmath $ k$}})-m_n } { ( x+y)(x+e({\mbox{\boldmath $ k$}})-m_n)(y+e({\mbox{\boldmath $ k$}})-m_n ) } \big[\mbox{im } \ , \tilde{d}_{\pi}(y,{\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$}})\big]\end{aligned}\ ] ] @xmath164   \nonumber \\ \nonumber \\ & & \times \int \limits_0^{\infty } dy \frac{x+y+e_{\delta}({\mbox{\boldmath $ k$}})-m_n }       { ( x+y)(x+e_{\delta}({\mbox{\boldmath $ k$}})-m_n )        ( y+e_{\delta}({\mbox{\boldmath $ k$}})-m_n ) } \big[\mbox{im } \ , \tilde{d}_{\pi}(y,{\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$}})\big ] \end{aligned}\ ] ] the coupling constants @xmath165 are defined by @xmath166     the momentum space unitary @xmath0 meson exchange @xmath1 potential at finite baryon density .
left panel without short range correlations .
right panel with short range correlations.,title="fig : " ]   the momentum space unitary @xmath0 meson exchange @xmath1 potential at finite baryon density . left panel without short range correlations .
right panel with short range correlations.,title="fig : " ]    there is still one more correction to be done in order to account for short range correlations .
so far we have replaced the free pion propagator by the renormalized one of eq .
( [ inmediumgf ] ) .
however , since the interaction between the @xmath44 bubble and the external nucleons is affected by correlations one should take @xmath167 ( see eq .
( [ vl ] ) ) instead of @xmath168 ^ 2 d_{\pi}$ ] .
thus we would get the series @xmath169 where @xmath170 and @xmath171 are the free and dressed pion propagators , respectively .
hence , the expressions for @xmath154 and @xmath155 get modified by including inside the @xmath172 integral the factors @xmath173    our results for @xmath0-meson exchange in the momentum space are shown if fig .  [ sigmaexchmediummod ] for both cases , no short range correlations ( left ) and with short range correlations ( right )
in this section we consider another sort of intermediate distance contributions to the @xmath1 force generated by the uncorrelated two pion exchange . the material presented here for the vacuum
@xmath1 scattering is standard and we merely generalize it to the nuclear medium .     the planar box diagrams involving the nucleon and @xmath3 intermediate states . ]    in the perturbative expansion of the @xmath1 force we must to take into account the planar and crossed box diagrams shown in figs .  [ planarmedium ] and  [ crossedmedium ] , respectively
. we will discuss the scalar - isoscalar part of this contributions only .
the contribution of the isovector exchange is small and can be found , for instance , in ref .
@xcite    in vacuum , the expression for the planar box diagrams , fig .
[ planarmedium ] , including the nucleon pole and @xmath151 and @xmath152 intermediate states reads @xmath174   \big[{\mbox{\boldmath $ \sigma$}}^{(1)}_n \cdot { \mbox{\boldmath $ k$}}\big]^{\dagger } \big[{\mbox{\boldmath $ \sigma$}}^{(2)}_m \cdot ( { \mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$}})\big ]   \big[{\mbox{\boldmath $ \sigma$}}^{(2)}_m \cdot { \mbox{\boldmath $ k$}}\big]^{\dagger }    \ , \big\{\mathcal{r}_{nm}^{(p)}({\mbox{\boldmath $ k$ } } , { \mbox{\boldmath $ q$ } } ) \
, i_{nm}^{(p)}\big\}\ ] ] in eq .
( [ planar ] ) the sum over @xmath175 is assumed and the spin transition operators are @xmath176 , @xmath177 .
the coupling constants @xmath165 are defined in eq .
( [ kappaconstants ] ) .
the isospin factors are given by @xmath178 for which we find @xmath179 the function @xmath180 in eq .
( [ planar ] ) contains the integration over the time - like component of the four vector @xmath181 .
@xmath182 where @xmath170 and @xmath183 is the pion and nonrelativistic baryon propagators , respectively .
@xmath183 is given by @xmath184 this integration can be carried out explicitely @xmath185[e_m - e]\big)(e_n+e_m-2 e )    \nonumber\\ & & +   \big(\omega_1 + \omega_2\big ) \big(2\omega_1\omega_2+[e_n+e_m-2 e]^2\big)\big ]   \nonumber\\ & & \times \frac{1}{\omega_1\omega_2 ( \omega_1 + \omega_2 ) }   \left\{\frac{1 } { e_n+e_m- 2 e - i 0^+}\right\ }     \nonumber\\ & & \times   \frac{m_n}{e_n [ \omega_1 + e_n - e][\omega_2 + e_n - e ] } \nonumber\\ & & \times \frac{m_m}{e_m [ \omega_1 + e_m - e][\omega_2 + e_m - e ] } \end{aligned}\ ] ] recall that @xmath64 , @xmath186 , @xmath187 in @xmath1 c.m .
frame and in eq .
( [ rnmfull1 ] ) @xmath188 are the on - shell energies of intermediate pions and baryons with c.m .
energy of the initial nucleons @xmath189 .    in the nuclear matter @xmath190 , and the expression for @xmath191 is obtained by using the dispersion representation for the in - medium pion propagator @xmath192 where @xmath193 is given by eq .
( [ rnmfull1 ] ) with substitution @xmath194 and @xmath195 .     the crossed box diagrams with nucleon and @xmath3 intermediate states .
]    for the crossed box diagrams , fig .
[ crossedmedium ] , we have @xmath196   \big[{\mbox{\boldmath $ \sigma$}}^{(1)}_n \cdot { \mbox{\boldmath $ k$}}\big]^{\dagger } \big[{\mbox{\boldmath $ \sigma$}}^{(2)}_m \cdot { \mbox{\boldmath $ k$}}\big ]   \big[{\mbox{\boldmath $ \sigma$}}^{(2)}_m \cdot({\mbox{\boldmath
$ k$}}+{\mbox{\boldmath $ q$}})\big]^{\dagger }   \ , \big\{\mathcal{r}_{nm}^{(c)}({\mbox{\boldmath $ k$ } } , { \mbox{\boldmath $ q$ } } ) \ , i_{nm}^{(c)}\big\}\ ] ] the isospin factors are the following @xmath197 and the vacuum expression for @xmath198 is given by @xmath199 the integration can be done analytically and our result reads @xmath200   \nonumber\\ & & \times \frac{1}{\omega_1\omega_2 ( \omega_1 + \omega_2 ) }     \nonumber\\ & & \times   \frac{m_n}{\tilde{e}_n [ \omega_1 + \tilde{e}_n - e][\omega_2 + \tilde{e}_n - e ] } \nonumber\\ & & \times \frac{m_m}{e_m [ \omega_1 + e_m - e][\omega_2 + e_m - e ] } \end{aligned}\ ] ] here @xmath201 , @xmath202 and @xmath203 are defined in eq .
( [ kin ] ) and @xmath204 the corresponding expression for the in - medium crossed box diagrams takes the form @xmath205     the momentum space uncorrelated two - pion exchange at finite baryonic density .
the left panel correspond to normal pion dressing without sr correlations and right panel with sr correlations.,title="fig : " ]   the momentum space uncorrelated two - pion exchange at finite baryonic density .
the left panel correspond to normal pion dressing without sr correlations and right panel with sr correlations.,title="fig : " ]    for the generic case of non - vanishing initial momenta the analytical structure of eq .
( [ rnmfull1 ] ) is driven by the term in figure brackets @xmath206 here we have to pay a special attention to the case where two intermediate @xmath1 states @xmath207 appear , because in time ordered perturbation theory these diagrams are generated by iterations of the opep in a lippmann - schwinger equation ( lse ) . considering the @xmath1 intermediate state only
, the iterated tpe can be easily identified and comes from the nucleon pole in eq .
( [ rnmfull ] ) ( in the lower half of the complex plane ) corresponding to @xmath208 with @xmath209 @xmath210\ ] ] expanding eq .
( [ nucleonpole ] ) in powers of @xmath211 @xmath212 and inserting the leading order result in eq .
( [ planar ] ) one can get the second order term in the non - relativistic lippmann - schwinger equation .
it is instructive to derive the contributions of the two remaining poles ( in lower half - plane ) from the pion propagators @xmath213 after the expansion of this result in powers of @xmath211 we get @xmath214 in this limit eq .
( [ pipoles ] ) cancels exactly the corresponding crossed box diagram with two intermediate nucleons in the isoscalar channel .
indeed , considering the crossed box diagram with two intermediate nucleons , the leading term of @xmath215 in the @xmath211 expansion is given by @xmath216 note that for the isovector @xmath4 exchange because of the different sign of the @xmath217 term in @xmath218 , @xmath219 in eqs .
( [ ipl ] ) and  ( [ icr ] ) these two contributions would add .
the cancellation discussed above hold also in the nuclear medium .
the result of the vacuum scalar - isoscalar @xmath1 force generated by the planar and crossed box diagrams , with the nucleon pole diagrams excluded , is shown in the left panel of fig .
[ uncormainmedium ] by the solid curve .
it is in agreement with ref .
@xcite where it was shown that the consistent use of the cut off regularization in both the correlated two pion exchange , and uncorrelated two pion exchange , together with the contribution of a repulsive @xmath220-exchange lead to a scalar - isoscalar potential in good agreement with the argonne  @xcite potential in the whole range of relevant distances .
the corresponding results for the nuclear matter are shown in fig .
[ uncormainmedium ] ( left ) for three densities @xmath142 , @xmath143 and @xmath74 .
qualitatively the behavior is similar to the vacuum case but with a strong enhancement toward the small momentum transfer .
as we have seen this feature is generic in present calculations .
uncorrelated @xmath4 exchange becomes extremely attractive at intermediate distances .
the spin sum over intermediate baryon states of the @xmath221 operator in eqs .
( [ planar ] ) and  ( [ crossed ] ) gives in the scalar channel @xmath222 ^ 2\ ] ] where @xmath223 @xmath224 , @xmath225 @xmath226 and @xmath227 @xmath228 .
now we again wish to take into account the short range correlations @xmath229 ^ 2   \nonumber \\
\longrightarrow \tilde{\mathcal{w}}_l(k+q ) \tilde{\mathcal{w}}_l(k ) \left [ ( \widehat{{\mbox{\boldmath $ k$}}+ { \mbox{\boldmath $ q$ } } } ) \cdot \hat{{\mbox{\boldmath $ k$ } } } \right]^2   \nonumber \\
+ \left [   \tilde{\mathcal{w}}_l(k+q)\tilde{\mathcal{w}}_t(k )   +   \tilde{\mathcal{w}}_l(k ) \tilde{\mathcal{w}}_t(k+q)\right ] \left\ {   1 - \left [ ( \widehat{{\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } } )   \cdot \hat{{\mbox{\boldmath $ k$}}}\right]^2 \right\ } \nonumber \\
+ \tilde{\mathcal{w}}_t(k+q ) \tilde{\mathcal{w}}_t(k ) \left\ { 1 + \left [ ( \widehat{{\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } } ) \cdot \hat{{\mbox{\boldmath $ k$ } } } \right]^2 \right\ } \nonumber \\\end{aligned}\ ] ] where @xmath230 where @xmath231 is the lindhard function and @xmath167 and @xmath232 are given by eq .
( [ vl ] ) and  ( [ vt ] ) , respectively .
by comparing the results in fig .  [ opepmedium ] and fig .  [ opepmediummodif ]
we observe that the effect of correlations has been essential and reduces drastically the medium effects found in fig .
[ opepmedium ] without corrections .
we observe that the medium corrections weaken the strength of the central part of the opep . on the other hand the effect of the medium corrections on the tensor part
are more moderate .
the @xmath0 exchange presents similar features as we can see in fig .
[ sigmaexchmediummod ] .
there we see that the medium effects in the absence of short range correlations are rather large and increase the strength of the interaction in about a factor two at @xmath233 .
however , as soon as the correlations are taken into account the medium modifications are reduced drastically and at @xmath233 just reduce the strength in about 25 percent .
the medium effects in the uncorrelated two pion exchange are also relatively small , with respect to the size of the interaction in the vacuum , see fig .
14 ( right ) .
once again the consideration of the correlations has been essential and reducing the size of the medium effects
. however , the fact that the strength of this interaction is bigger than that of the correlated two pion exchange makes the absolute correction relevant and of the same size as the corrections discussed before .
the modification of the @xmath1 interaction at normal nuclear matter density in momentum space ( left panel ) and in configuration space ( right panel ) . shown
are the central part of the opep ( solid curve ) , tensor part of the opep ( dashed curve ) , correlated @xmath0-meson exchange ( dot - dashed curve ) and uncorrelated @xmath4 exchange ( dot - dot - dashed curve).,title="fig : " ]   the modification of the @xmath1 interaction at normal nuclear matter density in momentum space ( left panel ) and in configuration space ( right panel ) .
shown are the central part of the opep ( solid curve ) , tensor part of the opep ( dashed curve ) , correlated @xmath0-meson exchange ( dot - dashed curve ) and uncorrelated @xmath4 exchange ( dot - dot - dashed curve).,title="fig : " ]    since we have taken also the effect of correlations in the free part of the interaction and any realistic calculation of the binding energy of matter will explicitly account for these correlations , the use of our modified in - medium potentials would lead to double - counting if any of these methods is used . for this purpose the relevant results from this work should be the differences between the in - medium potential and the free one . in this case
we start already with one bubble and the correlations account for the repulsion between the external leg and the one in the bubble .
the two external legs still have to be correlated and this will be done with the use , for instance , of the ordinary bethe golstone equation .
after this discussion of the medium effects in the different terms , we show the differences in momentum space between the medium and free parts of the different terms in fig .
[ uncormainmediummm ] ( left ) .
the results are shown at @xmath233 and we see that these corrections are moderate , but they could have a relevant role in the binding of matter . in order to have a qualitative idea of the relevance of these corrections to the potentials we rewrite them in coordinate space and show these results in fig .
[ uncormainmediummm ] ( right ) .
the effects seem sizeable at short distances , but this will be irrelevant in any realistic calculation of binding energies since the consideration of the short range correlations between the external legs will make this interaction inoperative .
more interesting is the strength of the corrections around 1 fm , and there we see that all them are relatively small , of the order of 20 mev or less , the biggest one being the central part of the ope .
however , given the size of the empirical scalar isoscalar attraction @xcite which is of the order of 20 mev at intermediate densities , the corrections found here are not negligible .
it would be thus interesting to see the effects of the results obtained here in observables like the binding of nuclear matter and other properties , which we hope to stimulate with the results obtained here .
in this paper we studied the modification of the one pion exchange , as well as two pion exchange potential inside a nuclear medium . for this purpose
we separated the two pion exchange into an uncorrelated two pion exchange and the correlated two pion exchange .
we study both in the scalar isoscalar channel , which is by large the most important one generated by this interaction .
the correlated two pion exchange gave rise to the equivalent of the @xmath0 exchange in other models and we studied the medium modification to it . on the other hand , for the uncorrelated two pion exchange we have followed a traditional approach in which only terms with at least one intermediate @xmath3 state are considered .
one of the important findings here was the effect played by the nn short range correlations which drastically moderated the medium corrections to the potential . in the absence of these , the corrections where unrealistically large . yet , even if relatively moderate , the medium corrections found in this paper are sizeable enough to have relevant repercussions in the binding and other properties of nuclear matter and we would like to encourage calculations in this direction .
we would like to acknowledge a.  ramos for a critical reading of the manuscript and useful suggestions .
this work is partly supported by dgicyt contract number bfm2003 - 00856 , and the e.u .
euridice network contract no .
hprn - ct-2002 - 00311 .
this research is part of the eu integrated infrastructure initiative hadron physics project under contract number rii3-ct-2004 - 506078 .
the set of feynman diagrams shown in fig .
[ 3pinnvertcharge ] appear in the construction of the @xmath40 transition amplitude .
the corresponding @xmath41 vertex functions are given by     contact terms appearing in the construction of the @xmath40 transition amplitude . ]
i l_b & = & \frac{d+f}{2 }   \frac{2 \sqrt{2}}{12 f_{\pi}^3 } \big [ { \mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ q$}}_2 + 2 { \mbox{\boldmath $ k$}}_2 - { \mbox{\boldmath $ k$}}_1 ) \big ] \\
i l_c & = & \pm \frac{d+f}{2 }   \frac{2}{12 f_{\pi}^3 } \big [ { \mbox{\boldmath $ \sigma$}}(2 { \mbox{\boldmath $ q$}}_1 + { \mbox{\boldmath $ k$}}_1 + { \mbox{\boldmath $ k$}}_2 ) \big ] \\
i l_d & = & \frac{d+f}{2 }   \frac{2 \sqrt{2}}{12 f_{\pi}^3 } \big [ { \mbox{\boldmath $ \sigma$}}(2{\mbox{\boldmath $ q$}}_1 + { \mbox{\boldmath $ k$}}_1 + { \mbox{\boldmath $ k$}}_2 ) \big ] \\
i l_e & = & \frac{d+f}{2 }   \frac{2 \sqrt{2}}{12 f_{\pi}^3 } \big [ { \mbox{\boldmath $ \sigma$}}(2{\mbox{\boldmath $ q$}}_2 + { \mbox{\boldmath $ k$}}_1 + { \mbox{\boldmath $ k$}}_2 ) \big ] \\
i l_f & = & 0\end{aligned}\ ] ]
in this appendix we demonstrate the cancellation of the off shell part of the correlated two pion exchange in the scalar isoscalar channel . for simplicity
we consider the nucleon intermediate states only .
the generalization to @xmath151 and @xmath152 is straightforward .
the diagrams are shown in fig .
[ corrtpecrmedium ] and corresponding amplitudes are given by @xmath235 \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big[{\mbox{\boldmath $ k$}}^2   \mathcal{u}(k)\big ] s_f(p_1-k ) d^2(k ) d(k+q ) \\
-i t_a^{(1 ) } & = & 2 \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }   \big[v_{n\pi\pi}(t ) \big ] \nonumber\\ & & \times    \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big[{\mbox{\boldmath $ k$}}^2   \mathcal{u}(k)\big ] s_f(p_1-k ) d(k ) d(k+q )   \\
-i t_a^{(2 ) } & = & 2 \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }    \big[v_{n\pi\pi}(t ) \big ] \nonumber\\ & & \times \int \frac{d^4 k}{(2 \pi)^4 } \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[{\mbox{\boldmath $ k$}}^2   \mathcal{u}(k)\big ] s_f(p_1-k ) d^2(k ) \\ \nonumber\end{aligned}\ ] ] @xmath236 \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big[({\mbox{\boldmath $ k$}}+ { \mbox{\boldmath $ q$}})^2   \mathcal{u}(k+q)\big ] s_f(p_1-k ) d(k ) d^2(k+q )   \\
-it_b^{(1 ) } & = & 2 \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }   \big[v_{n\pi\pi}(t ) \big ] \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big[({\mbox{\boldmath $ k$}}+ { \mbox{\boldmath $ q$}})^2   \mathcal{u}(k+q)\big ] s_f(p_1-k ) d(k ) d(k+q ) \\
-it_b^{(2 ) } & = & 2 \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 } \big[v_{n\pi\pi}(t ) \big ] \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[({\mbox{\boldmath $ k$}}+ { \mbox{\boldmath $ q$}})^2   \mathcal{u}(k+q)\big ] s_f(p_1-k ) d^2(k+q ) \nonumber \\ -i
t_c & = & -i t_a(k \to l , q \to - q ,
p_1 \to p_2 ) \nonumber \\ -i t_d & = & -i t_b(k \to l , q \to - q , p_1 \to p_2)\end{aligned}\ ] ] where @xmath237 is the triangle loop integral @xmath238   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ l$}}\big ]   s_f(p_2-l ) d(l - q ) d(l)\ ] ] and the on - mass - shell @xmath4 scattering amplitude is @xmath239    @xmath240 \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big [ ( 3 { \mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}+ 2 { \mbox{\boldmath $ k$}}^2 )    \mathcal{u}(k)\big ] s_f(p_1-k ) d(k ) d(k+q ) \\
-i t_f & = & \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }
\big[v_{n\pi\pi}(t ) \big ] \\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big [ ( { \mbox{\boldmath $ q$}}^2 - 2 { \mbox{\boldmath $ k$}}^2 - { \mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$ } } )   \mathcal{u}(k+q)\big ] s_f(p_1-k ) d(k )   d(k+q ) \nonumber \\
t_g & = & -i t_e(k \to l , q \to - q ,
p_1 \to p_2 ) \\
-i t_h & = & -i t_f(k \to l , q \to - q ,
p_1 \to p_2 ) \\
-i t_i & = &   \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }   \big[v_{n\pi\pi}(t ) \big ]   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ q$}}-2{\mbox{\boldmath $ k$ } } ) \big ] \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big[{\mbox{\boldmath $ k$}}^2 \mathcal{u}(k)\big ] s_f(p_1-k ) d^2(k ) \\
-i t_j & = & - \left(\frac{d+f}{2 f_{\pi}}\right)^4   \frac{1}{f_{\pi}^2 }   \big[v_{n\pi\pi}(t ) \big ] \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ q$}}+2{\mbox{\boldmath $ k$ } } ) \big ]    \big[{\mbox{\boldmath $ k$}}^2 \mathcal{u}(k)\big ] s_f(p_1+q - k ) d^2(k ) \\ -i t_k & = & -i t_i(k \to l , q \to - q ,
p_1 \to p_2 ) \\
-i t_l & = & -i t_j(k \to l , q \to - q ,
p_1 \to p_2 ) \end{aligned}\ ] ]    in the present case @xmath64 and for @xmath241 the amplitude @xmath242 cancels exactly @xmath243 .
the same cancelation is found for @xmath244 and @xmath245 .
the sum of @xmath246 and @xmath247 is given by @xmath248 \nonumber\\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[{\mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}\ ,   \mathcal{u}(k)\big ] s_f(p_1-k ) d(k ) d(k+q)\end{aligned}\ ] ] the sum of @xmath249 and @xmath250 takes the form @xmath251 \\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) { \mbox{\boldmath $ q$}}\mathcal{u}(k+q)\big ] s_f(p_1-k ) d(k ) d(k+q ) \nonumber\end{aligned}\ ] ] from this the sum of four diagrams is given by @xmath252 \\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$ } } ) \big ]   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]    \big [ { \mbox{\boldmath $ q$}}^2 \mathcal{u}(k+q )   + { \mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}\big\{\mathcal{u}(k+q)-\mathcal{u}(k)\big\ } \big ]
s_f(p_1-k ) d(k ) d(k+q ) \nonumber\end{aligned}\ ] ] alternatively , the spin flip parts of eqs .
( [ mm ] ) and  ( [ mn ] ) can be canceled if we change in eq .
( [ mn ] ) @xmath253 and using again the fact that @xmath64 and the properties of the lindhard function @xmath254 we get @xmath255 \\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[{\mbox{\boldmath $ \sigma$}}{\mbox{\boldmath $ k$}}\big ]   \big[{\mbox{\boldmath $ \sigma$}}({\mbox{\boldmath $ k$}}+{\mbox{\boldmath $ q$}})\big ]   \big
[ { \mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}\mathcal{u}(k)\big ] s_f(p_1-k ) d(k ) d(k+q ) \nonumber\end{aligned}\ ] ] and @xmath256 \\ & & \times   \int \frac{d^4 k}{(2 \pi)^4 }   \big[({\mbox{\boldmath $ k$}}^2 + { \mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}){\mbox{\boldmath $ k$}}{\mbox{\boldmath $ q$}}\big ] \mathcal{u}(k ) s_f(p_1-k ) d(k ) d(k+q ) \nonumber\end{aligned}\ ] ] finally @xmath257 \\ & & \times   \int \frac{d^4 l}{(2 \pi)^4 }   \big[({\mbox{\boldmath $ l$}}{\mbox{\boldmath $ q$}}-{\mbox{\boldmath $ l$}}^2){\mbox{\boldmath $ l$}}{\mbox{\boldmath $ q$}}\big ] \mathcal{u}(l ) s_f(p_2-l ) d(l ) d(l - q ) \nonumber\end{aligned}\ ] ] as one can see in the limit @xmath241 we get an exact cancellation of the off mass shell part of the @xmath4 amplitude and only @xmath258 and @xmath259 are left .
a. fabrocini , and s. fantoni , in _ first course on condensed matter _ , acif series , vol.viii , d. prosperi , s. rosati , and g. violini eds .
, world scientific , singapore ( 1986)87 . s. fantoni and a. fabrocini , in _ microscopic quantum many - body theories and their applications _ edited by j. navarro and a. polls , lecture notes in physics , vol .
510 ( springer - verlag , berlin 1998 ) .
e.  oset , h.  toki and w.  weise , phys .
* 83 * ( 1982 ) 281 .
w.  r.  gibbs and b.  f.  gibson , ann .
nucl .  part .
* 37 * ( 1987 ) 411 .
t.  e.  o.  ericson and w.  weise , `` pions and nuclei '' , oxford university press , 1988 . c.  m.  chen , d.  j.  ernst and m.  b.  johnson , phys .
c * 47 * ( 1993 ) 9 .
j.  nieves , e.  oset and c.  garcia - recio , nucl .
phys .  a * 554 * , 554 ( 1993 ) .
m.  nuseirat , m.  a.  k.  lodhi and w.  r.  gibbs , phys .
c * 58 * ( 1998 ) 314 .
t.  s.  h.  lee and r.  p.  redwine , ann .
nucl .  part .
sci .   * 52 * ( 2002 ) 23 .
r.  machleidt , k.  holinde and c.  elster , phys .
* 149 * ( 1987 ) 1 .
h.  muther , a.  faessler , m.  r.  anastasio , k.  holinde and r.  machleidt , phys .
c * 22 * , 1744 ( 1980 ) .
t.  s.  jensen and a.  f.  miranda , phys .
c * 55 * , 1039 ( 1997 ) .
o.  jaekel , h.  w.  ortner , m.  dillig and c.  a.  z.  vasconcellos , nucl .
a * 511 * , 733 ( 1990 ) .
s.  kamefuchi , l.  oraifeartaigh and a.  salam , nucl .
* 28 * , 529 ( 1961 )                m.  m.  kaskulov and h.  clement , phys .
c * 70 * , 014002 ( 2004 ) .
m.  m.  kaskulov and h.  clement , phys .
c * 70 * , 057001 ( 2004 ) ; h.  a.  clement , m.  m.  kaskulov and e.  a.  doroshkevich , int .
j.  mod .
a * 20 * , 674 ( 2005 ) .
d.  jido , e.  oset and j.  e.  palomar , nucl .
a * 694 * , 525 ( 2001 ) .