theories with extra spacetime dimensions have drawn considerable attention ever since the original proposal by kaluza and klein .
there has been renewed interest in such theories since the emergence of string theory .
several new ideas in this context have been proposed and have interesting consequences for particle phenomenology and cosmology @xcite . in these higher - dimensional models
, spacetime is usually taken to be a product of a four - dimensional spacetime and a compact manifold of dimension @xmath0 .
while gravity can propagate freely through the extra dimensions , standard model particles are confined to the four dimensional spacetime .
observers in this three spatial - dimensional wall ( a  3-brane " ) will measure an effective planck scale @xmath1 , where @xmath2 is the volume of the compact space .
if @xmath2 is large enough it could make planck scale of the order of tev , thus removing the hierarchy between the planck and the weak scale .
subsequently , randall _ et al . _
@xcite proposed a higher - dimensional scenario that is based on nonfactorizable geometry and accounts for the hierarchy without introducing large extra dimensions .
however , the braneworld model itself is not stable and it was shown in ref .
@xcite that by introducing a scalar field in the bulk , the modulus - namely the brane separation in the rs model - can be stabilized without any fine - tuning .
assumption of negligibly small scalar backreaction on the metric in the gw approach prompted further work in this direction , where the modification of the rs metric due to backreaction of the bulk fields has been derived ( see @xcite ) .
the stability issues in such cases have been reexamined for time - dependent cases @xcite ; also the effect of gauge fields or higher form fields have been studied in several works ( see @xcite ) .    in an effort to search for the signatures of extra dimensions , the roles of the kaluza - klein modes of different bulk fields on the phenomenology at the standard model brane are of crucial importance .    for the five - dimensional rs model ,
@xcite determined the bulk scalar kk modes and their self - interactions .
it is found that due to the exponential redshift factor in the rs model , kk scalar modes in this spacetime have tev scale mass splitting and inverse tev couplings ( see @xcite ) .
this is in sharp contrast to the kk decomposition in product spacetimes , which for large compactified dimensions , give rise to a large number of light kk modes ( see @xcite ) with a very small coupling with brane fields . due to this very distinct feature , the rs model has interesting consequences @xcite .
motivated by string theory and other extra - dimensional models where one can have several extra dimensions , in this paper we extend the results of the bulk scalar model in five - dimensional rs space - time to arbitrary number of warped dimensions and have obtained the kk decomposition of the scalar kk masses .
we have shown that in these multiply warped models we have much larger number of kk modes than the five - dimensional rs counterpart with effective couplings in the inverse tev range .
we have also discussed possible numerical values for various parameters in our theory and have used them to get possible numerical values of low - lying kk mode masses in our multiply warped model , showing explicit localization in tev brane .
our results also establish a general formula for determining these kk masses and couplings in the presence of any arbitrary number of extra dimensions .
the paper is organized as follows : we give a brief explanation for six - dimensional doubly warped spacetime in sec .
[ bmsix ] , the calculation for bulk scalar field has been done in this six - dimensional spacetime in sec .
[ bmbulk ] .
the same calculations have been finally extended to higher - dimensional spacetime with arbitrary number of extra dimensions in secs .
[ bmseven ] and [ bmbulk7 ] .
the paper ends with a short discussion of our results .
in this section we shall discuss doubly compactified six - dimensional spacetime with @xmath3 orbifolding along each compactified direction . for a detailed discussion
we refer the reader to @xcite .
the manifold under consideration is given by , @xmath4\times s^{1}/z_{2}$ ] @xcite .
we let the compactified dimensions to @xmath5 and @xmath6 , respectively .
the noncompactified dimensions are taken to be , @xmath7 .
the moduli along the compact dimensions are given by @xmath8 and @xmath9 , respectively .
the corresponding metric ansatz is taken as @xmath10+r_{z}^{2}dz^{2},\ ] ] with @xmath11 .
thus we have four orbifold fixed points , which are given by @xmath12 , respectively .
the total bulk - brane action could be given by @xmath13 + \int d^{4}xdydz \left[v_{3}\delta(z)+v_{4}\delta(z-\pi)\right ] \\ & & s_{4}=\sum _
{ i=1}^{2}\sum _ { j=1}^{2 } \int d^{4}xdydz \sqrt{-g_{vis}}\left(l - v\right)\delta ( y - y_{i})\delta ( z - z_{j}).\end{aligned}\ ] ] here the brane potentials in general have the particular functional dependence @xmath14 and @xmath15 .
finally the full six - dimensional einstein s equation is given by , @xmath16 in the above expression @xmath17 , @xmath18 are bulk indices , @xmath19 , @xmath20 run over the usual four spacetime coordinates given by @xmath21 .
the quantities @xmath22 and @xmath23 are the metric in @xmath24 and @xmath25 hypersurfaces , respectively .
the line element derived from the above einstein s equation turns out to be @xcite @xmath26+r_{z}^{2}dz^{2}.\ ] ] in the above line element we have the following identification for the constants @xmath27 and @xmath28 given by @xmath29 the boundary terms lead to the brane tensions and using the einstein s equation across the two boundaries at @xmath30 , @xmath31 , respectively , thus we readily obtain @xmath32 thus the two 4-branes situated at @xmath30 and @xmath31 would have a @xmath6-dependent brane tension .
the fact that the two tensions are equal and opposite is reminiscent of the original rs - form .
similarly we get the brane tensions for other two 4-branes as @xmath33 here @xmath34 were introduced to account orbifolding along the @xmath6-direction and with @xmath35 being a constant , the resulting hypersurface should have only a constant energy density .
the fact that @xmath36 is dependent on the coordinate @xmath6 makes the two hypersurfaces for y orbifolding to have a z - dependent energy density .
the 3-brane located at @xmath37 suffers no warping and can be identified with the planck brane .
the other three can be valid choices for standard model ( visible ) brane .
however , if we assume that there is no brane having lower energy scale than ours , we are forced to identify the sm brane to be located at @xmath38 .
the suppression factor on the tev brane can be given by @xmath39 the desired suppression of @xmath40 on the tev brane can be obtained by choosing different combinations of the parameters @xmath28 and @xmath27 .
however we also have an extra relation as presented in eq .
( [ bm5 ] ) , which shows that in order to avoid large hierarchy between the two moduli @xmath8 and @xmath9 , either of the two parameters @xmath28 and @xmath27 must be large and the other should be small .
for example , we can easily assume @xmath41 and @xmath42 .
however a small hierarchy also exists in the original rs model , where there exists an one order of magnitude hierarchy between @xmath43 and @xmath27 , satisfying planck - to - tev scale warping by @xmath44 .
a natural question that arises with this discussion is whether stabilization of these moduli to the desired values is possible .
for the five - dimensional rs model this has been shown in @xcite by introducing a bulk stabilizing scalar field and tuning the vev of the scalar field at the boundaries .
the modulus in the theory is stabilized near planck length without any fine - tuning .    in this case
as well we can adopt a similar procedure by introducing a bulk stabilizing scalar field . again choosing appropriate vev at the boundary , we can stabilize @xmath8 and @xmath9 to desired values @xcite . in our six - dimensional braneworld scenario with @xmath5 and @xmath6 dependence ,
the action for the scalar field can be expressed such that @xmath45 where the coupling parameters , @xmath46 tend to infinity as the scalar field approach to the following values , @xmath47 , @xmath48 , @xmath49 and @xmath50 .
we take @xmath51 . now following ref .
@xcite we can obtain the equation of motion in the separable form as @xmath52= \left(r_{y}^{2}b^{2}m^{2}-p\right)\chi ( z),\end{aligned}\ ] ] where we have made the following decomposition , @xmath53 and @xmath54 is the separability constant .
also in the above expression prime denote differentiation with respect to @xmath5 and dot denotes differentiation with respect to z. finally , the above equations with appropriate boundary condition @xcite can be solved , which when substituted into the action leads to an effective potential for the moduli as @xmath55,\end{aligned}\ ] ] where we have used the following shorthand notations , @xmath56 , @xmath57 and @xmath58 . then solving the equations , @xmath59 and @xmath60 and then through the second derivatives with respect to @xmath61 and @xmath27 we can arrive at the minimum values of @xmath28 and @xmath27 . as an illustrative example we can start with @xmath62 and @xmath63 , leading to @xmath64 and @xmath65 .
note that these values can resolve the gauge hierarchy problem .
thus along this line any higher - dimensional braneworld models can have their moduli stabilized .
from now on we shall assume that such a stabilization has been performed and all the moduli hence forth will have those stabilized values .
in this analysis , following the stabilizing bulk scalar model , we have assumed that the backreaction of the bulk stabilizing field is negligibly small .
moreover from the action of the bulk stabilizing scalar it may be observed that at the boundaries the stabilizing scalar tends to their vevs @xmath66 when the coupling @xmath67 tends to infinity . this is exactly similar to the five - dimensional counter part of the goldberger - wise calculation of modulus stabilization . as a result at the boundaries ,
the stabilizing scalar is frozen to different values @xmath66 and hence does not contribute to the dynamics of the model .
in this section we carry out the kaluza - klein decomposition of a nongravitational bulk scalar field propagating in the spacetime described by eq .
( @xmath68 ) in the spirit of the work @xcite with bulk scalar field .
we find that in these multiply warped spacetime the sm brane contains larger number of tev scale scalar kk modes than the five - dimensional rs model .
this has significant phenomenological consequences @xcite .
we consider a free scalar field in the bulk for which the action is given by @xmath69 where @xmath70 with @xmath71 is given by eq .
( @xmath68 ) , and m is of order of @xmath72 . after an integration by parts ,
this can be written as @xmath73,\end{aligned}\ ] ] where @xmath74 .
now we make the following substitution for kk decomposition , @xmath75 the following normalization conditions are imposed on the fields @xmath19 and @xmath20 , @xmath76 @xmath77 the differential equation satisfied by the function @xmath78 is @xmath79 where @xmath80 stands for kk mode mass eigenvalue .
the above differential equation can be further simplified and cast to the following form , @xmath81 the above equation can be solved in terms of bessel functions of first and second order as @xmath82,\ ] ] with @xmath83 representing an overall normalization .     with extra - dimension parameter @xmath5 .
the vertical line represents the @xmath31 line showing the fact that the quantity @xmath84 is maximum at @xmath31 , the position of tev brane.,width=384,height=288 ]    now we can proceed much further .
the mass modes determined by @xmath80 must be real .
this reality condition imposed on the mass modes requires the differential operator on the left hand side of eq .
( [ bm14a ] ) to be self - adjoint .
this self - adjointness imply that derivatives of @xmath85 should be continuous at the orbifold fixed points .
these gives two conditions on the parameters @xmath80 and @xmath86 , expressed as , @xmath87 \left[2y_{2}(x_{m}e^{-c\pi})+x_{m}e^{-c\pi}y'_{2}\left(x_{m}e^{-c\pi}\right)\right ] \nonumber \\ & -&\left[2y_{2}\left(x_{m}\right)+x_{m}y'_{2}\left(x_{m}\right)\right ] \left[2j_{2}(x_{m}e^{-c\pi})+x_{m}e^{-c\pi}j'_{2}\left(x_{m}e^{-c\pi}\right)\right]\end{aligned}\ ] ] where @xmath88 .
since to make planck scale down to tev scale we should have @xmath89 .
then the above equation reduces to the following form , @xmath90 .
then for light mode masses we have @xmath91 to be order of unity @xcite .
this keeps @xmath92 also order of unity . then from eq .
( [ bm14a ] ) as well as from fig.[fig1 ] we observe that the modes @xmath85 are larger near the 3-brane at @xmath31 , which makes these low - mass kaluza - klein modes to be found preferentially near the @xmath31 region ( see fig.[fig1 ] ) .
thus , with the tev brane being situated at @xmath31 , we observe that the low - mass kk modes are exponentially suppressed and hence confined to the tev brane .
also for these low - lying kk mass modes the coefficient @xmath86 is of the order of @xmath93 , which shows that we can ignore the @xmath94 part compared to @xmath95 , while performing integrals involving @xmath84 .
similar analysis for @xmath96 yields @xmath97 the solution , apart from an overall normalization , can be expressed as @xmath98 h_{\sqrt{5/2}kz}\left(\frac{-10k^{2}+b_{m}^{2}r_{z}^{2}(1+\cosh(2k\pi))}{10k^{2}}\right ) \nonumber \\ & + & e_{m}exp\left[-\frac{5}{2}k^{2}z^{2}\right]~_{1}f_{1 } \left(-\frac{-10k^{2}+b_{m}^{2}r_{z}^{2}(1+\cosh(2k\pi))}{10k^{2}},\frac{1}{2},\frac{5k^{2}z^{2}}{2}\right),\end{aligned}\ ] ] where @xmath99 is the kummer confluent hypergeometric function and @xmath100 is the hermite polynomial of degree n. then from fig.[fig2 ] we observe that this function is also maximum at @xmath101 . hence the low - lying kk mass modes are confined to the tev brane located at @xmath101 . from the behavior of both @xmath84 and @xmath102 we find that all the low - lying kk mass modes are confined to @xmath103 brane , i.e. , the tev brane .
a possible experimental signature of the bulk scalar kk modes can originate via coupling of the bulk scalar to dihiggs in the form @xmath104 . for @xmath105
the dominant decay channels are @xmath106 and @xmath107 which leads to multijets as final product which though may be difficult to differentiate from the qcd back ground @xcite .
also when the mass of bulk scalar is in the range of @xmath108 to @xmath109 gev then enhanced production of @xmath110 occurs .
moreover for bulk scalar mass in the range @xmath111 to @xmath108 gev we have a relatively larger cross sections for the diphoton channel . in this region due to
small mixing the branching ratio is dominated by @xmath106 and @xmath107 .
the diphoton channel is a very promising search channel as branching ratio remains more or less at constant level even up to @xmath112 threshold @xcite .
this might become possible if the lhc runs extends the diphoton searches for invariant masses above existing @xmath113gev .
+ to determine the parameters of the solution we proceed as follows : we want @xmath114 to be real , as it appears in the mass modes .
thus self - adjointness also applies in this case .
this implies that derivatives of @xmath96 to be continuous around the orbifold fixed points . at @xmath101 ,
this is trivially satisfied irrespective of the quantity @xmath115 .
however at @xmath116 all the terms are suppressed by @xmath117 , thus the self - adjointness there leads to @xmath118 at large values of @xmath6 , confluent hypergeometric series have a large value .
being in the denominator the term can be neglected for practical purposes .     with extra - dimension parameter @xmath6 .
the graph clearly shows the fact that the quantity @xmath102 is maximum at @xmath101 , the position of tev brane.,width=384,height=288 ]    using the above equations we readily obtain the following action for the field @xmath119 as , @xmath120 where we have the following expression for the element @xmath121 , @xmath122 and @xmath123 as , @xmath124 now from the previous discussion we have the solution for these two sets of functions @xmath78 and @xmath125 , which can be used to determine @xmath121 and @xmath126 in order to obtain the masses of the kk modes by evaluating the quantity @xmath127 .
in contrast to the five - dimensional situation ( see @xcite ) where the masses of the bulk fields appear as a diagonalized mass matrix , in this case the bulk field @xmath128 manifests itself to some four - dimensional observer as an infinite kk tower with mass being determined by the the quantity @xmath127 such that a scalar @xmath129 has a mass @xmath130 after an appropriate diagonalization procedure .
the solution for @xmath78 is presented in equation ( [ bm14a ] ) .
a similar solution was obtained by wise et .
al ( see @xcite ) except for the fact that we have bessel functions of second order .
following their discussion we can argue in a similar manner that the lightest kk modes have mass parameter @xmath80 suppressed exponentially with respecanct to the the scale @xmath131 appearing in eq .
( [ bm8 ] ) . since we have taken @xmath131 to be order of planck scale and @xmath28 to be around @xmath132 , by stabilization these mass modes @xmath80 are in the tev range . also from the solution
we could observe that the modes @xmath78 are larger in the region @xmath31 .
this has been explained earlier through graphical presentation of the function @xmath84 .
the solution for @xmath125 has been presented in eq .
( [ bm15a ] ) .
though we have argued following the graphical presentation of the function @xmath102 , we now provide a theoretical motivation for our above - mentioned results .
the solution has an overall factor of @xmath133 $ ] and we see that the solution has maximum value around @xmath101 .
hence the bulk field being a product of these two functions @xmath85 and @xmath134 as shown in eq .
( [ bm9 ] ) has mass parameter in the tev range and has maximum value around @xmath38 .
now from sec .
[ bmsix ] this is precisely the sm brane .
hence , the bulk field has a maximum in the sm brane ; i.e. , the kk modes are most likely to be found in that region where @xmath80 and @xmath114 are in the tev range .
this sets the stage for kk excitations to have tev scale mass splitting on the sm brane .
now we would like to compute some low - lying kk mode masses numerically .
for that we need to fix some parameters , @xmath27 , @xmath28 , @xmath9 and the bulk mass of the scalar field @xmath131 .
we shall take the bulk mass to be in planck scale .
then we can determine the remaining parameters , by making the following demands : ( a ) if we have a gauge boson field in this multiply warped scenario , its lowest massive kk modes should lead to @xmath135 and @xmath136 boson masses , ( b ) the suppression @xmath137 as presented in eq .
( [ branem01 ] ) should be @xmath138 , and finally ( c ) the hierarchy between @xmath8 and @xmath9 should be small .
the kk mode of the gauge boson in this multiply warped spacetime can be obtained from ref .
@xcite .
this desired mass for @xmath135 and @xmath136 boson @xmath139gev can be obtained with @xmath140 and @xmath141gev , about 14 times smaller compared to planck scale . the other parameters @xmath27 and @xmath28
can be determined using small hierarchy between @xmath8 and @xmath9 along with desired warping of @xmath140 .
this finally leads to , the following estimation : @xmath142 , @xmath143 and the ratio between moduli being @xmath144 .
the suppression factor turns out to be @xmath145 .
thus we will calculate the low - lying kk masses for our bulk scalar field with these sets of parameters ( see table.[tab1 ] ) .
we now present the self - interactions of the bulk scalar field . from the four - dimensional point of view
these self - interactions can induce couplings between the kk modes . in this
case self - couplings of the light modes are suppressed by the warp factor and hence if the planck scale set the six - dimensional couplings , the low - lying kk modes have tev range self - interactions .
we present the interaction term in the action with coupling parameter @xmath146 such that @xmath147 where the coupling @xmath146 is of the order of unity
. then we can expand in modes and the self - interactions of light kk states become @xmath148 thus the effective four dimensional coupling constants are @xmath149 which reduces to , @xmath150^{2m}\int _
{ 0}^{\pi}(\beta _ { q})^{2m}dz\ ] ] in the large @xmath151 and @xmath152 limit .
hence we observe that the relevant scale for four - dimensional physics is not the scale set by planck scale , i.e. , @xmath17 , but this is @xmath153 .
hence the kk reduction has lead the couplings from planck scale to the tev scale by the warp factor on the sm brane located at @xmath38 .
.the masses of the kk modes of the scalar field are given in gev units .
we have chosen the following values , @xmath141 gev , @xmath142 , @xmath143 .
some representative masses of low - lying kk modes are given .
[ cols="^,^,^,^ " , ]     from the above discussion we now try to obtain some bounds on the parameters in our model , e.g. , @xmath8 , @xmath9 from the requirement of precision electroweak test .
for that purpose we can use the same setup and put a bulk gauge boson whose kk modes can be detected in precision electroweak tests .
we define a quantity denoted by @xmath154 where @xmath155 is the mass of @xmath135 gauge boson and @xmath156 is the mass of higher kk modes of the bulk gauge boson and @xmath157 is the effective four - dimensional gauge coupling along with @xmath158 to be the gauge couplings for higher kk modes .
then from ref .
@xcite we could argue that for precision electroweak test we should have @xmath159 with @xmath160 percent confidence level . from this result
we can get the following bounds on the parameters of this model , @xmath161gev .
this leads to a bound on @xmath162 as well by assuming a small hierarchy between the two moduli as , @xmath163gev . from ref .
@xcite it can be easily verified that this bound is respected by gauge couplings and kk mode masses .
thus these multiply warped models indeed satisfy precision electroweak tests .
with an aim to arrive at a generic result we shall now try to extend our analysis with one more extra dimension . for that purpose
we start with a seven - dimensional spacetime where the spacelike dimensions are successively warped . in other words
the manifold of interest could be given by @xmath164\right\rbrace \times s^{1}/z_{2}\right]\times s^{1}/z_{2}$ ] .
then the total brane - bulk action can be given by @xmath165 \nonumber \\
& + & \int d^{4}xdydzdw \left[v_{3}\delta(z)+v_{4}\delta(z-\pi)\right ] \nonumber \\ & + & \int d^{4}xdydzdw \left[v_{5}\delta(y)+v_{6}\delta(y-\pi)\right],\end{aligned}\ ] ] with appropriate actions ( @xmath166 ) for 12 possible 4-branes at the edges @xmath167 , @xmath168 and @xmath169 .
we also have eight possible 3-branes at the corners @xmath170 . by natural extension of the method
as illustrated in the previous section we get the line element and other parameters such that @xcite , @xmath171+r_{z}^{2}dz^{2}\right\rbrace + \re _ { w}^{2}dw^{2 } \nonumber \\
\ell ^{2}&=&\frac{-\lambda _ { 7}\re _ { w}^{2}}{15 } \nonumber \\ k&= & \frac{\ell r_{z}}{\re _ { w}\cosh(\ell \pi ) } \nonumber \\
c&=&\frac{\ell r_{y}}{\re _ { w}\cosh(k\pi)cosh(\ell \pi)}=\frac{kr_{y}}{r_{z}\cosh(k\pi)}\end{aligned}\ ] ] it may be of interest that the 5-brane at @xmath172 does not represent a flat metric ( y and z dependencies ) . in order to obtain substantial warping along the @xmath173 direction ( from @xmath172 to @xmath174 ) , one need to make @xmath175 substantial ( same order of magnitude as rs scenario ) .
the seven - dimensional or triply warped model has a structure analogous to that of six - dimensional one , not only in the , form of functional dependence but also on the nature of warping .
this method can easily be extended to even higher dimensions .
also note that the orbifolding requires branes situated at edges of n - dimensional hypercube with 3-branes at the corners .
if one of the direction suffers a large warping any other direction should have small warping so that there is no large hierarchy coming from the moduli . in this case
also we have several candidates for our sm brane .
however applying the fact that no brane should have less energy than ours , leads to @xmath176 to be sm brane .
following the methods of previous sections , we shall carry out the kaluza - klein decomposition of a bulk scalar field propagating in the spacetime given by eq .
( @xmath177 ) . as in the previous section in this case as well we can write the bulk scalar field in terms of product of four functions . by making kk decomposition
we again end up with kk mass modes having tev scale masses and splittings .
the action for the bulk scalar field in this seven - dimensional spacetime can be given as @xmath178.\ ] ] from the line element as given by eq .
( @xmath179 ) , we readily obtain the following form for the action @xmath180,\end{aligned}\ ] ] where @xmath181 .
we make the following substitution for the bulk field : @xmath182 we also impose the following normalization for the functions @xmath183 , @xmath184 and @xmath185 , @xmath186 now applying integration by parts to the integral as presented in eq .
( @xmath187 ) we readily obtain @xmath188 \nonumber \\ & -&\frac{1}{2}\frac{1}{r_{y}^{2}}\left[\sum_{pqrabc } \frac{\cosh^{4}(\ell w)}{\cosh^{4}(\ell \pi)}\frac{\cosh^{3}(kz)}{\cosh^{3}(k\pi ) } \phi _ { pqr}\phi _ { abc}\beta _ { q}\beta _ { b}\gamma _ { r}\gamma _ { c}\alpha _ { p } \partial _ { y}\left(e^{-4\sigma}\partial _ { y}\alpha _ { a}\right)\right ] \nonumber \\ & -&\frac{1}{2}\frac{1}{r_{z}^{2}}\left[\sum_{pqrabc } \frac{\cosh^{4}(\ell w)}{\cosh^{4}(\ell \pi)}e^{-4\sigma } \phi _ { pqr}\phi _ { abc}\alpha _ { p}\alpha _ { a}\gamma _ { r}\gamma _ { c}\beta _ { q } \partial _ { z}\left(\frac{\cosh^{5}(kz)}{\cosh^{5}(k\pi)}\partial _ { z}\beta _ { b}\right)\right ] \nonumber \\ & -&\frac{1}{2}\frac{1}{\re _ { w}^{2}}\left[\sum \frac{\cosh^{5}(kz)}{\cosh^{5}(k\pi)}e^{-4\sigma } \phi _ { pqr}\psi _ { abc}\alpha _ { p}\alpha _ { a}\beta _ { q}\beta _ { b}\gamma _ { r } \partial _ { w}\left(\frac{\cosh^{6}(\ell w)}{\cosh^{6}(\ell \pi ) }    \partial
_ { w}\gamma _ { c}\right)\right ] \nonumber \\ & + & \frac{1}{2}m^{2}\left[\sum \frac{\cosh^{6}(\ell w)}{\cosh^{6}(\ell \pi)}\frac{\cosh^{5}(kz)}{\cosh^{5}(k\pi)}e^{-4\sigma } \phi _ { pqr}\phi _ { abc}\alpha _ { p}\alpha _ { a}\beta _ { q}\beta _ { b}\gamma _ { r}\gamma _ { c}\right ] \big \rbrace.\end{aligned}\ ] ] then we make the following choice for the differential equations satisfied by the functions @xmath189 , @xmath102 and @xmath190 , @xmath191 the first equation as presented in ( [ bmm26a ] ) can be solved and has an identical solution as that obtained in the previous section . however , for convenience we rewrite the solution , @xmath192\ ] ] the second equation as given by eq .
( [ bmm26b ] ) has the same solution as presented in eq .
( [ bm15a ] ) but we rewrite it here , @xmath193 h_{\sqrt{5/2}kz}\left(\frac{-10k^{2}+b_{q}^{2}r_{z}^{2}(1+\cosh(2k\pi))}{10k^{2}}\right ) \nonumber \\ & + & e_{q}exp\left[-\frac{5}{2}k^{2}z^{2}\right]~_{1}f_{1 } \left(-\frac{-10k^{2}+b_{q}^{2}r_{z}^{2}(1+\cosh(2k\pi))}{10k^{2}},\frac{1}{2},\frac{5k^{2}z^{2}}{2}\right)\end{aligned}\ ] ] the third eq .
( [ bmm26c ] ) has the following solution along with an overall normalization , @xmath194 h_{\sqrt{3}\ell w } \left(\frac{-12\ell ^{2}+c_{r}^{2}\re_{w}^{2}(1+\cosh(2\ell \pi))}{12\ell ^{2}}\right ) \nonumber \\ & + & f_{r}exp\left[-3\ell ^{2}w^{2}\right]~_{1}f_{1 } \left(-\frac{-12\ell ^{2}+c_{r}^{2}\re_{w}^{2}(1+\cosh(2\ell \pi))}{24\ell ^{2}},\frac{1}{2},3\ell ^{2}w^{2}\right)\end{aligned}\ ] ] here also we have @xmath195 and @xmath196 to be bessel functions of first and second order respectively .
along with these @xmath100 represents hermite polynomials and @xmath99 is the kummer confluent hypergeometric series .
the arbitrary constants @xmath197 , @xmath198 and @xmath199 can be determined by the self - adjoint criteria and have the following expressions @xmath200 hence our final expression for the action is given by @xmath201 \label{bm27a } \\
m_{pqrabc}&=&\left\lbrace a_{p}^{2}\delta _ { pa}\delta _ { qb}\delta _ { rc}+b_{p}^{2}p_{rc}\delta _ { pa}\delta _ { qb } + c_{p}^{2}p_{qb}q_{rc}\delta _ { pa}+~m^{2}p_{pa}q_{qb}r_{rc}\right\rbrace , \label{bm27b}\end{aligned}\ ] ] where we have defined the following quantities , @xmath202 now from the previous discussion we can find the solution for the three sets of functions @xmath78 , @xmath125 and @xmath203 which in turn determine @xmath123 , @xmath204 and @xmath205 .
therefore from these three functions the explicit expression for kk mass modes can be determined from the quantity @xmath206 as given by eq .
( [ bm27b ] ) . in this case
as well the bulk field @xmath128 manifests itself to some four - dimensional observer as a scalar @xmath207 whose mass is determined by eqs .
( [ bm27a ] ) and ( [ bm27b ] ) .
the solution for @xmath78 as presented in eq .
( [ bm26a ] ) has the same nature as obtained by wise _
( see @xcite ) . in this case as well lightest kk modes have mass modes determined by @xmath80 , suppressed exponentially with respect to the the scaler mass @xmath131 which we have taken to be order of planck scale . thus these mass modes @xmath80 are in the tev range whereas m is of order of @xmath72 .
the solution for @xmath125 and @xmath203 has been presented in eqs .
( [ bm26b ] ) and ( [ bm26c ] ) .
the solution can be seen to include exponential factors such as @xmath133 $ ] , @xmath208 $ ] and we see that when mass parameter @xmath114 is of the order of tev , solutions have maximum value around @xmath101 as obtained earlier in sec .
[ bmbulk ] as well . from the solution of @xmath203 it is evident that the solution has maximum value around @xmath174 .
hence the bulk field being a product of these three functions @xmath209 , @xmath210 and @xmath211 as shown in eq .
( [ bmm4 ] ) , has mass parameter in the tev range and has maximum value to find the modes around @xmath212 which is the location of the sm brane . also the bulk field is maximum in the sm brane ; i.e. , the kk modes are most likely to be found in the tev region as the @xmath80 , @xmath114 and @xmath213 are in the tev range . along with the above line of arguments we could in principle
have plotted all the functions @xmath209 , @xmath210 and @xmath211 and for all of them we have the functions to take maximum value at @xmath31 , @xmath101 and @xmath174 , precisely at the location of the tev brane .
for completeness we present the self - interactions of the bulk scalar field in this seven - dimensional spacetime . from the four - dimensional point of view
these self - interactions induce couplings between the kk modes . in this case
also the effective self - couplings are suppressed by the warp factor and if the planck scale sets the six - dimensional couplings and the low - lying kk modes have tev range self - interactions .
we present the interaction term in the action with coupling parameter @xmath146 such that @xmath214 where the coupling @xmath146 is of the order of unity
. then we could expand in modes and hence the self - interactions of light kk states are given by @xmath215 the effective four - dimensional coupling constants therefore are being given by @xmath216 which in the large @xmath151 , @xmath152 and @xmath217 limit reduces to @xmath218^{2 m } \int _ { 0}^{\pi } dz(\beta _ { q})^{2 m } \int _ { 0}^{\pi } dw(\gamma _ { r})^{2m}.\end{aligned}\ ] ] hence we observe that the relevant scale for four - dimensional physics is not the scale set by planck scale but @xmath219 .
hence the kk reduction lead to the tev scale couplings by the warp factor on the sm brane located at @xmath212 .    now
this result can easily be extended to any higher dimension spacetime . for n extra dimensions we can write the action for the bulk field as , @xmath220\ ] ] where @xmath70 with @xmath221 is given by a generalization of eq .
( @xmath177 ) , and m is of order of @xmath72 .
then the kk splitting for the bulk field can be expressed as the following decomposition , @xmath222 thus among these n extra dimensions , one will have the solution given by eq . ( [ bm26a ] ) , and then the other @xmath223 solutions are being given by generalization of eq .
( [ bm26b ] ) such that the numerical values will be different but form of the solution remains unaltered .
for nth extra dimension @xmath224 the solution for the mode can therefore be expressed as @xmath225 exp\left[-\frac{1}{2}n k^{2}w^{2}\right ] h_{\frac{\sqrt{3+n}}{2}kw } \left(\frac{-6k^{2}-2nk^{2}+m_{r}^{2}r^{2}(1+\cosh(2k \pi))}{2k^{2}(3+n)}\right ) \nonumber \\ & + & f_{r}exp\left[-\frac{3}{2}k^{2}w^{2}\right]\exp \left[-\frac{1}{2}n k^{2}w^{2}\right ]
\nonumber \\ & \times & _ { 1}f_{1 } \left(-\frac{-6k^{2}-2nk^{2}+m_{r}^{2}r^{2}(1+\cosh(2k \pi))}{4k^{2}(3+n)},\frac{1}{2},\frac{1}{2}(3+n)k^{2}w^{2}\right)\end{aligned}\ ] ] hence the bulk field as viewed by a four - dimensional observer leads to a mass matrix whose components can be obtained by solving the eigenvalue problem as presented by each separable functions in the expansion given by eq .
( [ bmm6 ] ) .
also all these eigenvalues have tev scale masses and the bulk field also has maximum value at @xmath226 , which is the sm brane .
hence the standard model particles can be taken as low - lying kaluza - klein modes of a bulk field propagating in any number of extra - dimensional spacetime .
the effective self - coupling is this case turns out to be @xmath227    @xmath228^{2 m } \int _ { 0}^{\pi } dz(\beta _ { q})^{2 m } \int _ { 0}^{\pi } dw(\gamma _ { r})^{2m}\cdots.\ ] ]    thus finally we have obtained the kk mass modes and their self - interactions for n extra dimensions .
we have also observed that in all these cases the kk mass modes and self - interactions are suppressed by the warp factor near the sm brane and hence all are in tev scale .
hence this properties can be used to search for the tev range kk mass modes and self - interactions in next generation colliders .
in this paper we generalize the work presented in ref .
@xcite on the bulk scalar field to determine its kk modes and the effective self - interaction in a multiple warped spacetime . for arbitrary number of extra dimensions ,
we have derived the expressions for the kk mode masses and their self - interactions .
various components determining the masses are in the tev range because of the warp factor suppression .
moreover the bulk scalar field has been shown to have maximum value at the sm brane .
hence the low - lying kk modes for the bulk scalar fields lie in the tev range with inverse tev self - coupling .
thus the appearance of kk mode masses and couplings at tev scale are generic features of warped dimensional models with any number of extra warped dimensions as long as we want to resolve the gauge hierarchy problem without introducing any hierarchical moduli .
we have also introduced the moduli stabilization mechanism in these multiply warped models and have obtained the stabilized values for the moduli .
then we have presented a compact and generic formula to determine all the mass modes and their couplings for models with any arbitrary number of warped extra dimensions .
this work now can be extended to other forms of bulk fields , which in turn may lead to the possibility of identifying various standard model particles as the low - lying kk excitation of various bulk fields , where the small warping in multiple directions can explain mass splitting in standard model particles as discussed in @xcite and @xcite .
the close spacing of the low - lying kk modes along with enhanced coupling makes it likely for them to be seen as a series of close - lying resonances . in order to investigate the role of kk mass modes through collider - based experiments
, we consider the interaction of various modes with thermselves , i.e. , self - interactions , and we have obtained that all of them are suppressed to the tev scale by the warp factor . also from the numerical values of masses for low - lying kk modes
, we readily observe that the masses in the standard rs model get split into infinite number of mass modes , with very close spacings , which is a very interesting feature of these multiply warped models and can be probed in future runs of lhc .
the other things to be noted are that with stronger coupling of the kk modes of the bulk scalar field , one expects the decay widths to be larger , and thus the peaks to be broader as we go to higher and higher dimensions .
the nature of the line shapes , therefore , will be an interesting benchmark to distinguish between higher - dimensional and lower - dimensional kk signals if such excitations appear during the high - luminosity runs of the lhc .
s.c . is funded by a spm fellowship from csir , government of india .
i. antoniadis , phys .
b * 246 * , 377 ( 1990 ) ; + j.d .
lykken , phys .
d * 54 * , r3693 ( 1996 ) ; + r. sundrum , _ ibid .
_ * 59 * , 085009 ( 1999 ) ; + k.r .
dienes , e. dudas and t. gherghetta , phys .
b * 436 * , 55 ( 1998 ) .
o. dewolfe , d.z .
freedman , s.s .
gubser and a. karch , phys .
d * 62 * , 046008 ( 2000 ) ; + g.f .
giudice , r. ratazzi and j.d .
wells , nucl .
b544 * , 3 ( 1999 ) ; + e.a .
mirabelli , m. perelstein and m.e .
peskin , phys .
lett . * 82 * , 2236 ( 1999 ) .
w.d . goldberger and m.b .
wise , phys .
d * 60 * , 107505 ( 1999 ) ( hep - ph/9907218 ) ; + s. kachru , m.b .
schulz and e. silverstein , phys .
d * 62 * , 045021 ( 2000 ) ; + e. dudas and m. quiros , nucl . phys . *
b721 * , 309 ( 2005 ) .
giudice , r. ratazzi and j.d .
wells , nucl .
* b544 * 3 ( 1999 ) ; + e.a .
mirabelli , m. perelstein and m.e .
peskin , phys .
lett . * 82 * 2236 ( 1999 ) ; + t. han , j.d . lykken and r.j.zhang , phys . rev . *
d59 * 105006 ( 1999 ) ; + j.l .
hewett , phys .
82 * 4765 ( 1999 ) .


quasielastic nucleon knockout reactions have been extensively investigated with the aim of obtaining precise information on nuclear structure .
the present work focuses on exclusive proton - induced @xmath8 and @xmath9 processes , whereby the residual @xmath10 nucleus is left in the discrete part of its energy spectrum .
a sketch of the @xmath8 reaction is given in fig .
[ fig : schem_repres ] . for a historical overview of the research into proton - induced nucleon emission off nuclear targets the reader
is referred to refs .
@xcite .     reaction .
the incoming proton undergoes `` soft '' initial - state interactions with the target before knocking out a bound proton through the primary high - momentum - transfer @xmath11 scattering .
both the scattered and the ejected proton suffer final - state interactions while leaving the nucleus . the scattered and the ejected proton are detected in coincidence , while the residual nucleus remains unobserved.,scaledwidth=50.0% ]    in quasielastic @xmath0 scattering the projectile is elastically scattered from a single bound nucleon in the target nucleus , resulting in the struck nucleon being knocked out of the target nucleus .
this relatively simple reaction mechanism of one `` hard '' nucleon - nucleon collision is obscured by the `` soft '' initial- and final - state interactions ( ifsi ) of the incident and two outgoing nucleons with the nuclear medium .
consequently , every model for @xmath8 and @xmath9 reactions has two issues to address  : first , the description of the `` hard '' wide - angle scattering that leads to the ejection of the struck nucleon and second , the distorting mechanisms of the `` soft '' small - angle ifsi .
concerning the treatment of the `` hard '' @xmath12 scattering part , essentially two methods exist .
a so - called `` cross - section factorized '' approximation @xcite can be adopted so that the nucleon - nucleon scattering cross section enters as a multiplicative factor in the differential @xmath0 cross section .
some results of exclusive @xmath8 measurements interpreted with this cross - section factorized form can be found in refs .
the inclusion of spin - dependence in the description of the ifsi , however , breaks this factorization scheme .
in that situation an alternative technique can be used  : the amplitude factorized form of the cross section @xcite . in this approach
, the two - body @xmath12 interaction can be approximated by the interpolation of phase shifts @xcite from free elastic @xmath12 scattering .
various phenomenological forms to fit the amplitudes have also been used in the past .
traditionally , the nucleon - nucleon scattering matrix has been parametrized in terms of five lorentz invariants @xcite , a method usually dubbed as the ia1 model or the spvat ( scalar , pseudoscalar , vector , axial vector , tensor ) form of the @xmath12 scattering matrix .
differential cross section calculations adopting these five - term representations have been reported in refs  @xcite .
it should be noted , however , that although the spvat form gives reasonable predictions of @xmath0 observables , it is , in principle , not correct , as a five - term parametrization of the relativistic @xmath12 scattering matrix is inherently ambiguous @xcite .
tjon and wallace @xcite have shown that a complete expansion of the @xmath12 scattering matrix ( commonly called the ia2 model ) contains @xmath13 independent invariant amplitudes . to date , the only calculations employing this general lorentz invariant representation have been performed in the context of the relativistic plane wave impulse approximation ( rpwia ) , i.e. , a model which ignores all ifsi mechanisms @xcite .
the ifsi effects are typically computed by means of the distorted wave impulse approximation ( dwia ) theoretical framework @xcite .
generally , in a dwia approach the scattering wave functions of the incoming and two outgoing nucleons are generated by solving the schrdinger or dirac equation with complex optical potentials .
parametrizations for these optical potentials are usually not gained from basic grounds , but are obtained by fitting elastic nucleon - nucleus scattering data .
several optical - potential parameter sets @xcite have been used in the description of quasifree proton scattering off nuclei . in the past ,
both nonrelativistic and relativistic dwia versions have proven successful in predicting @xmath0 cross sections over a wide energy range ( @xmath14@xmath15  mev ) and for a whole scope of target nuclei .
most of the available calculations for the exclusive @xmath0 process addressed incident proton kinetic energies of a few hundred mev . in this work we aim at extending the formalism to scattering in the gev energy regime .
as a matter of fact , the majority of dwia frameworks rely on partial - wave expansions of the exact solution to the scattering problem , an approach which becomes increasingly cumbersome at higher energies . in this energy range
the eikonal approximation @xcite , which belongs to the class of high - energy semiclassical methods , offers a valid alternative for describing the ifsi .
nonrelativistic eikonal studies of the @xmath8 reaction used in combination with optical potentials can be found in refs .
@xcite .    in this paper
we propose a relativistic and cross - section factorized formalism based on the eikonal approximation for computing exclusive @xmath0 cross sections at incident proton energies in the few hundred mev to gev range .
the eikonal formalism is implemented relativistically in combination with optical potentials @xcite , as well as with glauber theory @xcite , which is a multiple - scattering extension of the eikonal approximation .
the two frameworks only differ in the way they treat the ifsi
. the main focus will be on the optical - potential approach , as this method turns out to be the most suitable for the description of the ifsi for the kinematical settings discussed in this work .
the paper is organized as follows . in secs .
[ sec : cross_sect_matr_elt ] and [ sec : rpwia ] the factorized cross section is derived in the rpwia formalism .
thereafter , the different methods to deal with the ifsi are developed in sec .
[ sec : ifsi ] .
[ sec : ifsi_factor ] is devoted to a presentation of optical - potential and glauber results for the ifsi factor .
this is a function which accounts for all ifsi effects when computing the @xmath0 observables .
the optical - potential predictions of our model are compared with cross section data that have been collected at pnpi and triumf in sec .
[ sec : cross_sect_results ] .
first , we present our calculations for the @xmath2c , @xmath3o , and @xmath4ca@xmath6 and @xmath16 pnpi data for @xmath1  gev incoming proton energies @xcite .
second , the cross sections for @xmath5he@xmath6 scattering at @xmath7  mev are compared to the triumf data by van oers _ et al .
finally , sec .
[ sec : concl ] states our conclusions .
in this section the formalism for the description of @xmath8 reactions is outlined .
the generalization to @xmath9 reactions is straightforward .
we conform to the conventions of bjorken and drell @xcite for the @xmath17 matrices and the dirac spinors , and take @xmath18 .
the four - momenta of the incident and scattered proton are denoted as @xmath19 and @xmath20 .
the proton momenta @xmath21 and @xmath22 define the scattering plane .
the four - momentum transfer is given by @xmath23 , where @xmath24 , @xmath25 , and @xmath26 are the four - momenta of the target nucleus , residual nucleus , and the ejected proton .
the standard convention @xmath27 is followed for the four - momentum transfer .    in the laboratory frame
, the fivefold differential cross section can be written as @xmath28 here , @xmath29 is the invariant matrix element which reflects the transition between the initial and final states .
the hadronic recoil factor is given by @xmath30 with the energy transfer @xmath31 , the three - momentum transfer @xmath32 , and @xmath33 the angle between @xmath34 and @xmath35 .
the @xmath8 matrix element is given by @xmath36 where @xmath37 is the unknown two - body operator describing the high - momentum transfer `` hard '' @xmath11 scattering , @xmath38 the ground state of the even - even target nucleus and @xmath39 the discrete state in which the residual nucleus is left .
in coordinate space the matrix element takes on the form @xmath40 for the sake of brevity of the notations , only the spatial coordinates are explicitly written .      in this section , the @xmath8 matrix element of eq .
( [ eq : matrix_elt ] ) will be analyzed in the rpwia . in this approach , only one hard collision between the projectile and a bound nucleon
is assumed to occur , knocking the bound nucleon out of the target nucleus .
the modelling of the `` soft '' ifsi processes , which affect both the incoming and outgoing protons , will be considered in sec .
[ sec : ifsi ] .    in evaluating the @xmath8 matrix element of eq .
( [ eq : matrix_elt ] ) , a mean - field approximation for the nuclear wave functions is adopted .
we also assume factorization between the `` hard '' @xmath12 coupling and the nuclear dynamics . for reasons of conciseness ,
the forthcoming derivations are explained for the @xmath41 case .
the generalization to arbitrary mass number @xmath42 is rather straightforward .
the antisymmetrized @xmath43-body wave function in the initial state is of the slater determinant form @xmath44 details on the bound - state single - particle wave functions @xmath45 entering this mean - field @xmath43-body wave function can be found in appendix  [ app : bound ] .
the wave function of the incoming proton is given by a relativistic plane wave @xmath46 e^{i \vec{k } \cdot \vec{r } } \chi_{\frac{1}{2}m_{s } } = e^{i \vec{k } \cdot \vec{r } } \ : u ( \vec{k } , m_{s } ) \ ; .
\label{eq : relat_plane_wave}\end{aligned}\ ] ] the @xmath43-body wave function in the final state reads @xmath47 relative to the target nucleus ground state written in eq .
( [ eq : ini_wave_fctn_slater ] ) , the wave function of eq .
( [ eq : fin_wave_fctn_slater ] ) refers to the situation whereby the struck proton resides in a state `` @xmath48 '' , leaving the residual @xmath10 nucleus as a hole state in that particular single - particle level .
the outgoing protons are represented by relativistic plane waves .
since both the initial and the final wave functions are fully antisymmetrized , one can choose the operator @xmath49 to act on two particular coordinates ( @xmath50 and @xmath51 ) . without any loss of generality , the @xmath8 matrix element of eq .
( [ eq : matrix_elt ] ) can be written as @xmath52 with @xmath53 the levi - civita tensor . in the rpwia , @xmath54 inserting this expression in eq .
( [ eq : matrix_elt_part_coord ] ) one obtains @xmath55 there are @xmath56 possible choices ( permutations ) for the indices @xmath57 , @xmath58 ,  , all giving the same contribution to the matrix element .
accordingly , the above expression can be rewritten as @xmath59 because the bound - state wave functions are normalized to unity ( @xmath60 ) and @xmath61 , the matrix element can be further simplified to @xmath62 including a direct and an exchange term .
substitution of the general form of the scattering operator @xmath63 , with @xmath64 the @xmath12 scattering amplitude in momentum space , in the above expression ( [ eq : rpwia_dir_exch_matrix_elt ] ) leads to @xmath65 for the direct term and a similar expression for the exchange term . here
@xmath66 is the relativistic wave function for the bound nucleon in momentum space ( for details see appendix  [ app : bound ] ) and @xmath67 is the missing momentum . in order to arrive at a cross - section factorized expression for eq .
( [ eq : fivefold_diff_cross_section ] ) , the quasielastic off - shell proton - proton scattering matrix element will be related to the free on - shell proton - proton cross section . for this purpose , we insert the completeness relation @xmath68 = 1 \label{eq : completeness}\ ] ] in @xmath69 and obtain the following expression for the matrix element @xmath70 here , @xmath71 is the matrix element for free @xmath11 scattering @xmath72 factorization breaks down , even when ifsi are disregarded , owing to the negative - energy projection term . to recover factorization ,
the negative - energy projection term is neglected in the remainder of this work @xmath73    using the expression of the relativistic bound - nucleon wave function in momentum space given in appendix  [ app : bound ] , the @xmath74 contraction in eq .
( [ eq : pos_energy_proj_matr_elt ] ) reduces to @xcite @xmath75 where @xmath76 , @xmath77 indicates the spin projection of the spin spherical harmonic @xmath78 on a spin state @xmath79 , and the radial function in momentum space @xmath80 is given by @xmath81 with @xmath82 and @xmath83 the bessel transforms of the standard upper and lower radial functions of the bound - nucleon wave function in coordinate space ( see appendix  [ app : bound ] for details ) and @xmath84 .    upon squaring eq .
( [ eq : pos_energy_proj_matr_elt ] ) , the @xmath11 and nuclear bound - state parts get coupled by the summation over the intermediate spins @xmath85 and @xmath86  : @xmath87 after summation over @xmath57 , the struck nucleon s generalized angular momentum quantum number , the square of @xmath88 yields a @xmath89 , i.e. , becomes diagonal in @xmath85 .
thereby , use is made of the following identity @xmath90 this leads to the decoupling between the @xmath11 scattering and the bound - state part in the matrix element  : @xmath91 with @xmath92    the last factor in eq .
( [ eq : decoupled_squared_matrix_elt ] ) can be related to the free @xmath11 scattering center - of - mass cross section @xmath93 so that the rpwia differential @xmath8 cross section of eq .
( [ eq : fivefold_diff_cross_section ] ) can be written in the cross - section factorized form @xmath94 here , @xmath85 is the mandelstam variable for the @xmath11 scattering , not to be confused with the intermediate spin from the preceding eqs .
( [ eq : completeness])([eq : decoupled_squared_matrix_elt ] ) . in the numerical calculations which will be presented in the forthcoming sections , the free proton - proton cross section @xmath95 is obtained from the said code @xcite .
it is well known that the factorized rpwia result of eq .
( [ eq : factorized_rpwia_cross_sect ] ) adopts an oversimplified description of the reaction mechanism .
the momentum distribution @xmath96 , which represents the probability of finding a proton in the target nucleus with missing momentum @xmath97 , will be modified by the scatterings of the incoming and outgoing protons in the nucleus .
therefore it is necessary to incorporate the effects of these ifsi in the model .
first , in sec .
[ sec : dist_factoriz ] , the differential @xmath8 cross section is written in a factorized form taking ifsi effects into account .
next , the relativistic eikonal methods used for dealing with the ifsi effects in this work , are discussed in depth .
two methods will be used .
the relativistic optical model eikonal approximation ( romea ) is the subject of sec .
[ sec : romea ] , whereas the relativistic multiple - scattering glauber approximation ( rmsga ) is discussed in sec .
[ sec : rmsga ] .      in both versions of the relativistic eikonal framework for @xmath0 reactions presented here ( romea and rmsga ) ,
the antisymmetrized initial- and final - state @xmath43-body wave functions , @xmath98 \label{eq : dist_ini_wave_fctn}\ ] ] and @xmath99 \ ; , \label{eq : dist_fin_wave_fctn}\end{aligned}\ ] ] differ from their respective rpwia expressions of eqs .
( [ eq : ini_wave_fctn_slater ] ) and ( [ eq : fin_wave_fctn_slater ] ) through the presence of the operators @xmath100 , @xmath101 , and @xmath102 . these define the accumulated effect of all interactions that the incoming and emerging protons undergo in their way into and out of the target nucleus .
since the ifsi violate factorization , some additional approximations are in order .
first , only central ifsi are considered , i.e. , spin - orbit contributions are omitted .
further , the zero - range approximation is adopted for the `` hard '' @xmath12 interaction , allowing one to replace the coordinates of the two interacting protons ( @xmath50 and @xmath51 ) by one single collision point in the distorting functions @xmath100 , @xmath101 , and @xmath102 .
this leads to the distorted momentum - space wave function @xmath103 similar to eq .
( [ eq : fourier_bound_wave ] ) , but with the additional ifsi factor @xmath104 accounting for the soft ifsi effects .    now , along the lines of @xcite
, it is natural to define a distorted wave amplitude @xmath105 so that the distorted momentum distribution is given by the square of this amplitude , @xmath106 this distorted momentum distribution has the following properties .
first , it takes into account the distortions for the incoming and outgoing protons .
second , it reduces to the plane wave momentum distribution @xmath96 in the plane wave limit when assuming that @xmath107 satisfies the relation @xmath108 between the upper and lower components .    using the ansatz ( [ eq : dist_mom_dist ] ) for the distorted momentum distribution
, the differential @xmath8 cross section can be cast in the form @xmath109 it differs from the rpwia expression ( [ eq : factorized_rpwia_cross_sect ] ) through the introduction of a `` distorted '' momentum distribution @xmath110 .      as shown for example in refs .
@xcite , in the relativistic eikonal limit the scattering wave function of a nucleon with energy @xmath111 and spin state @xmath112 subject to a scalar ( @xmath113 ) and a vector potential ( @xmath114 ) takes on the form @xmath115 e^{i \vec{k } \cdot \vec{r } } e^{i s(\vec{r } ) } \chi_{\frac{1}{2 } m_{s } } \ ; , \label{eq : eikonal_scatt_wave}\ ] ] where the eikonal phase @xmath116 reads @xmath117 \biggr ] \ ; , \label{eq : eikonal_phase}\ ] ] with @xmath118 and the average momentum @xmath119 pointing along the @xmath120-axis . the central and spin - orbit potentials @xmath121 and @xmath122 in the above expression are determined by @xmath113 and @xmath114 and their derivatives .
in general , a fraction of the strength from the incident beam is removed from the elastic channel into the inelastic ones .
these inelasticities are commonly implemented by means of the imaginary part of the optical potential .    in evaluating the ifsi effects ,
three approximations are introduced .
first , the dynamical enhancement of the lower component of the scattering wave function ( [ eq : eikonal_scatt_wave ] ) , which is due to the combination of the scalar and vector potentials , is neglected .
second , the impulse operator @xmath123 is replaced by the asymptotic momentum @xmath124 of the nucleon . as mentioned before , the spin - orbit potential @xmath122 is also omitted .    as a result ,
the effects of the interactions of the incoming and outgoing protons with the residual nucleus are implemented in the distorted momentum - space wave function of eq .
( [ eq : dist_mom_space_wave_fctn ] ) through the following phase factors    @xmath125    @xmath126    @xmath127    [ eq : eikonal_ifsi_operators ]    with the @xmath120-axes of the different coordinate systems lying along the trajectories of the respective particles ( @xmath120 along the direction of the incoming proton @xmath21 , @xmath128 along the trajectory of the scattered proton @xmath22 , and @xmath129 along the path of the ejected nucleon @xmath34 ) and @xmath130 , @xmath131 , and @xmath132 the coordinates of the collision point @xmath133 in the respective coordinate systems .
the geometry of the scattering process is illustrated in fig .
[ fig : geometry ] .
the integration limits guarantee that the incoming proton only undergoes isi up to the point where the `` hard '' @xmath12 collision occurs and the outgoing protons are only subject to fsi after this `` hard '' collision .    ,
@xmath134 , and @xmath135 are the impact parameters for each of the three paths for a collision occurring at @xmath133 .
@xmath136 , @xmath137 , and @xmath138 are the @xmath120 coordinates of the collision point in the respective coordinate systems . @xmath139 and @xmath140 are the angles of the outgoing nucleons relative to the incoming proton direction.,scaledwidth=50.0% ]    it is worth remarking that the eikonal ifsi operators of eq .
( [ eq : eikonal_ifsi_operators ] ) are one - body operators , i.e. , they do not depend on the coordinates @xmath141 of the residual nucleons .
the normalization of the bound - state wave functions simplifies the ifsi factor ( [ eq : ifsi_factor ] ) considerably to @xmath142 in the romea case .    in the numerical calculations ,
we have employed the global @xmath143 parametrizations of cooper _ et al .
_ @xcite and the optical potential of van oers _ et al .
_ @xcite to describe the pnpi and triumf data , respectively .
hereafter , the @xmath8 calculations which adopt eq .
( [ eq : eikonal_ifsi_operators ] ) as a starting basis are labeled the relativistic optical model eikonal approximation ( romea ) .      in the romea approach
, all the ifsi effects are parametrized in terms of mean - field like optical potentials , i.e. , the ifsi are seen as a scattering of the nucleon with the residual nucleus as a whole . as the energy increases ,
shorter distances are probed and the scattering with the individual nucleons becomes more relevant .
for proton kinetic energies
@xmath144  gev , the highly inelastic and diffractive character of the underlying elementary proton - nucleon scattering cross sections makes the glauber approach @xcite more natural .
this method reestablishes the link between proton - nucleus interactions and the elementary proton - proton and proton - neutron scattering .
it essentially relies on the eikonal , or , equivalently , the small - angle approximation and the assumption of consecutive cumulative scatterings of a fast nucleon on a composite target containing `` frozen '' point scatterers ( nucleons ) .    in relativistic glauber theory ,
the scattering wave function of a nucleon with energy @xmath111 and spin state @xmath112 reads @xcite @xmath145 e^{i \vec{k } \cdot \vec{r } } \chi_{\frac{1}{2}m_{s } } \ ; . \label{eq : glauber_scatt_wave}\ ] ] the operator @xmath146 implements the subsequent elastic or `` mildly inelastic '' collisions of the fast nucleon with the `` frozen '' spectator nucleons @xmath147 \ ; , \label{eq : glauber_operator}\ ] ] where @xmath148 ensures that the nucleon only interacts with other nucleons if they are localized in its forward propagation path .
given the diffractive nature of @xmath149 collisions at gev energies , the profile function @xmath150 for central elastic @xmath149 scattering is parametrized in a functional form of the type @xmath151 at lower energies that part of the profile function proportional to @xmath152 is non - gaussian and makes significant contributions to nuclear scattering . rather than eq .
( [ eq : profile_fctn ] ) a parametrization in terms of the arndt @xmath12 phases @xcite is appropriate at lower energies . for the calculations presented here , which address higher energies ,
the gaussian - like real part of @xmath150 is the dominant contributor and the use of eq .
( [ eq : profile_fctn ] ) is justified .
the parameters in eq .
( [ eq : profile_fctn ] ) can be determined directly from elementary nucleon - nucleon scattering experiments and include the total @xmath149 cross sections @xmath153 , the slope parameters @xmath154 , and the ratios of the real to the imaginary part of the scattering amplitude @xmath152 .
we obtained these glauber parameters through interpolation of the data base available from the particle data group @xcite ( for more details , see ref .
@xcite ) . as in the romea framework
, only the central spin - independent contribution is retained and the impulse operator is replaced by the nucleon momentum .
the glauber operators in eq .
( [ eq : ifsi_factor ] ) take the following form    @xmath155 \ ; , \ ] ]    @xmath156   \ ; , \ ] ]    @xmath157 \ ; , \ ] ]    [ eq : glauber_ifsi_operators ]    where @xmath133 denotes the collision point and @xmath141 are the positions of the frozen spectator protons and neutrons in the target .
the @xmath158 , @xmath159 , and @xmath160 coordinate systems are defined as in the previous section .
the step functions make sure that the incoming proton can only interact with the spectator nucleons which it encounters before the `` hard '' collision and the outgoing protons can only interact with the spectator nucleons which they find in their forward propagation paths .
contrary to romea , the glauber ifsi operators of eq .
( [ eq : glauber_ifsi_operators ] ) are genuine @xmath42-body operators , so the integration over the coordinates of the spectator nucleons in eq .
( [ eq : ifsi_factor ] ) has to be carried out explicitly .
this makes the numerical evaluation of the glauber ifsi factor very challenging .
standard numerical integration techniques were adopted to evaluate the ifsi factor and no additional approximations , such as the commonly used thickness - function approximation , were introduced .
henceforth , we refer to calculations based on eq .
( [ eq : glauber_ifsi_operators ] ) as the relativistic multiple - scattering glauber approximation ( rmsga ) .
in this section , results for the ifsi factor ( [ eq : ifsi_factor ] ) are given for the knockout of nucleons from the fermi level in @xmath2c , @xmath3o , and @xmath4ca , at an incident energy @xmath161  gev and a scattered proton kinetic energy @xmath162  mev .
thereby , we adopt coplanar scattering angles @xmath163 on opposite sides of the incident beam , i.e. , kinematics corresponding with the pnpi experiment of ref .
all ifsi effects are included in the ifsi factor @xmath164 .
note that in the absence of initial- and final - state interactions the real part of the ifsi factor equals one , whereas the imaginary part vanishes identically .
the @xmath8 ifsi factor is a function of three independent variables @xmath165 .
the @xmath120-axis is chosen along the direction of the incoming beam @xmath21 , the @xmath166-axis lies along @xmath167 and the @xmath168-axis lies in the scattering plane defined by the proton momenta @xmath21 and @xmath22 . @xmath169 and @xmath170 denote the polar and azimuthal angles with respect to the @xmath120-axis and the @xmath168-axis , respectively .
the radial coordinate @xmath171 represents the distance relative to the center of the target nucleus .      to gain a better insight into the dependence of the ifsi factor on @xmath171 , @xmath169 , and @xmath170 , we calculated the contribution of the three distorting functions @xmath100 , @xmath101 , and @xmath102 to the ifsi factor . in figs .
[ fig : re.12c.p.edad.phi0.3d ] and [ fig : im.12c.p.edad.phi0.3d ] results are displayed for the computed real and imaginary part of @xmath172 , for proton emission from the fermi level in @xmath2c .
the results were computed within the romea framework , using the edad1 optical - potential parametrization of @xcite .
the @xmath169 dependence can be interpreted as follows . for a given @xmath171 , the distance that the incoming proton travels through the target nucleus before colliding `` hard '' with a target nucleon , decreases with increasing angle @xmath169 . as a consequence , small values of @xmath169
induce the largest isi . for the fsi of the scattered proton
, the opposite holds true and @xmath173 for a large @xmath171 value corresponds to an event whereby the `` hard '' collision transpires at the outskirts of the nucleus and the scattered proton has to travel through the whole nucleus before it becomes asymptotically free , thus giving rise to the smallest ( largest ) values for the real ( imaginary ) part of the ifsi factor . unlike the scattered proton , which moves almost collinear to the @xmath120-axis , the ejected nucleon leaves the nucleus under a large scattering angle @xmath140 .
hence , the fsi are minimal for @xmath169 close to @xmath174 or @xmath175 and maximal for @xmath169 around @xmath176 . finally , the @xmath169 dependence of the complete ifsi factor is the result of the interplay between the three distorting effects , with the strongest scattering and absorption observed at @xmath169 close to @xmath174 , @xmath176 , and @xmath175 .
figure  [ fig : re.40ca.n.edai.phi0.radial ] displays the real part of the ifsi factor as a function of @xmath171 at various values of @xmath169 .
the romea calculations were performed for the same kinematics as in figs .
[ fig : re.12c.p.edad.phi0.3d ] and [ fig : im.12c.p.edad.phi0.3d ] , and use the edai optical - potential fit of @xcite .    turning first to the upper left panel
, it can be inferred from the picture that the isi effects increase with growing @xmath171 for @xmath177 . for @xmath178 and increasing @xmath171 , initially , the growing distance the proton has to travel through the nucleus leads to a decrease of the real part of @xmath100 .
this is followed by an increase for larger @xmath171 up to @xmath179 .
this reduction in isi effects with increasing @xmath171 is brought about by the incoming proton s path through the nucleus moving away from the nuclear interior and closer to the less dense nuclear surface .
the other curves of the upper left figure illustrate a general trend for @xmath180  : as @xmath171 increases , the real part of the ifsi factor grows correspondingly . as can be appreciated from fig .
[ fig : re.40ca.n.edai.phi0.radial ] as well as from the previous figures , the global behavior of the @xmath101 factor describing the scattered proton s fsi can be related to that of the isi factor @xmath100 through the substitution @xmath181 .
this approximate symmetry can be attributed to the small scattering angle @xmath139 , i.e. , the scattered proton leaves the nucleus almost parallel to the incoming proton s direction . in the bottom left panel , the additional curve ( @xmath182 , i.e. , close to @xmath176 ) represents the situation of maximal fsi of the ejected nucleon . for this @xmath169 value , the path of the ejected nucleon passes through the center of the nucleus and the distance travelled through the nucleus increases with @xmath171 .
accordingly , the real part of @xmath102 is a monotonously decreasing function of @xmath171 .
the other extreme is the @xmath183 case where increasing @xmath171 means less fsi . for the other @xmath169 values
, the absorption reaches its maximum for some intermediate @xmath171 value .
again , the combination of @xmath100 , @xmath101 , and @xmath102 determines the total ifsi factor with the strongest attenuation predicted in the nuclear interior .
the dependence of the ifsi factor on the azimuthal angle of the collision point is quite straightforward .
one representative result is displayed in fig .
[ fig : re.16o.p.edad2.r3fm.3d ] . here , @xmath184 ( @xmath185 ) refers to a situation where the `` hard '' @xmath12 collision occurs in the upper ( lower ) hemisphere with respect to the @xmath186-plane ( see fig .  [
fig : geometry ] for a collision point located in the upper hemisphere ) . due to the cylindrical symmetry about the @xmath120-axis ,
the factor describing the isi of the incoming proton is independent of @xmath170 .
regarding the scattered proton , we observe the least fsi in the upper hemisphere , since the proton then avoids passing through the highly absorbing nuclear interior . for the ejected nucleon
, the contrary applies and the strongest fsi effects are found for @xmath187 . as the @xmath188-plane is defined as the scattering plane , the ifsi factor possesses the symmetry @xmath189 for coplanar scattering .
results for the emission from levels other than the fermi level have not been plotted here , as it turns out that the ifsi factors hardly depend on the single - particle level in which the struck nucleon resides .
the peculiar spatial characteristics of the different single - particle orbits have an impact on the observables , though . indeed , the distorted momentum - space wave function @xmath190 of eq .
( [ eq : dist_mom_space_wave_fctn ] ) is determined by the values of the ifsi factor folded with a relativistic bound - state wave function @xmath191 .
as the particles experience less ifsi close to the nuclear surface , one obtains a stronger reduction of the quasifree cross section for nucleon knockout from a level which has a larger fraction of its density in the nuclear interior .
this will become apparent in fig .
[ fig : delta_r ] , but even more so in sec .
[ sec : pnpi ] .
figure  [ fig : delta_r ] shows a function @xmath192 which represents the contribution of the nuclear region with radial coordinate @xmath171 to the differential cross section .
the procedure to calculate this function is similar to the method exposed in ref .
@xcite and is developed in appendix  [ app : radial_contrib ] .
comparison of the upper and lower panels illustrates that ifsi mechanisms make the @xmath8 cross sections reflect surface mechanisms , unlike the @xmath193 reaction where the weakly interacting electron probes the entire nuclear volume and only the outgoing proton interacts with the residual system .
apart from the shift to higher @xmath171 , the ifsi brings about a strong reduction in the magnitude of the cross sections , whereby the fermi level is least affected .
even though @xmath192 is concentrated in the surface region , the average density seen through this reaction still amounts to @xmath194  @xmath195 ( @xmath196  @xmath195 ) or @xmath197 ( @xmath198 ) of the central density in the case of @xmath199 knockout from @xmath2c ( @xmath4ca ) . in the case of emission from the fermi level ,
on the other hand , the average density is only @xmath200 ( @xmath201 ) of the central density for a @xmath2c ( @xmath4ca ) target .
also , the ifsi factors for neutron emission are almost identical to the corresponding ifsi factors for proton knockout and , as expected , the overall effect of ifsi is more pronounced for heavier target nuclei .      in this subsection
, we investigate the sensitivity of the computed ifsi factors to the adopted parametrizations for the optical potentials and compare the romea results with the rmsga predictions . as can be seen in fig .
[ fig : re.12c.n.romea_compar.radial_polar_azimuth ] , the ifsi factor depends on whether @xmath42-dependent ( edad1/edad2 ) or @xmath42-independent ( edai ) fits for the potentials are selected , but the global features are comparable . figure  [ fig : re.12c.p.rmsga.phi0.3d ] , as contrasted to fig .
[ fig : re.12c.p.edad.phi0.3d ] , demonstrates that the rmsga method adequately describes the isi of the incoming proton and the fsi of the scattered proton .
however , the discrepancies between romea and rmsga become significant in the calculation of the fsi of the ejected nucleon ( note the different scales in the bottom left panels of figs .
[ fig : re.12c.p.edad.phi0.3d ] and [ fig : re.12c.p.rmsga.phi0.3d ] ) , and therefore also in the complete ifsi factor .
the noted difference is attributed to the low ejectile kinetic energy ( @xmath202  mev for the specific case of fig .
[ fig : re.12c.p.rmsga.phi0.3d ] , and comparable values for knockout from other levels and other nuclei ) . at such low energies ,
the rmsga predictions are not realistic due to the underlying approximations , mostly the postulation of linear trajectories and frozen spectator nucleons .
so , for the kinematics discussed here , the romea method is to be preferred over the rmsga one , as the latter overestimates the distortion for the low - energetic ejectile nucleon .
c to the adopted choice for the parametrization of the optical potentials .
results of romea calculations with the edad1 ( solid curve ) , edad2 ( dashed curve ) , and edai ( dot - dashed curve ) optical potentials are shown .
kinematics as in fig .
[ fig : re.12c.p.edad.phi0.3d].,title="fig:",scaledwidth=50.0% ] c to the adopted choice for the parametrization of the optical potentials
. results of romea calculations with the edad1 ( solid curve ) , edad2 ( dashed curve ) , and edai ( dot - dashed curve ) optical potentials are shown .
kinematics as in fig .
[ fig : re.12c.p.edad.phi0.3d].,title="fig:",scaledwidth=50.0% ] c to the adopted choice for the parametrization of the optical potentials . results of romea calculations with the edad1 ( solid curve ) , edad2 ( dashed curve ) , and edai ( dot - dashed curve ) optical potentials are shown .
kinematics as in fig .
[ fig : re.12c.p.edad.phi0.3d].,title="fig:",scaledwidth=50.0% ]
the pnpi experiment @xcite was carried out with an incident proton beam of energy @xmath1  gev .
the scattered proton was detected at @xmath203 with a kinetic energy between @xmath204 and @xmath205  mev , while the knocked - out nucleon was observed at @xmath206 having a kinetic energy below @xmath207  mev .
figures  [ fig : diffcross.12c.p.romea.edai][fig : diffcross.40ca.n.romea.edai ] display a selection of differential cross section results as a function of the kinetic energy of the most energetic nucleon in the final state . the edai optical potential @xcite was used for the romea calculations .
the other parametrizations of ref .
@xcite produce similar predictions , whereas the rmsga approach fails to give an adequate description of the data because of the low kinetic energy of the ejected nucleon .
since the experiment of ref .
@xcite only measured relative cross sections , the romea results were normalized to the experimental data .
the romea calculations reproduce the shapes of the measured differential cross sections .
furthermore , comparison of the rpwia and romea calculations shows that the effect of the ifsi is twofold .
first , ifsi result in a reduction of the rpwia cross section that is both level and @xmath42 dependent . from the figures it is clear that ejection of a nucleon from a deeper lying level leads to stronger initial- and final - state distortions .
this reflects the fact that the incoming and outgoing nucleons encounter more obstacles when a deeper lying bound nucleon is probed .
the @xmath42 dependence also conforms with our expectations , i.e. , the ifsi effects are larger for heavier nuclei . besides the attenuation , the ifsi also make the measured missing momentum different from the initial momentum of the struck nucleon . as can be inferred from fig .
[ fig : momdistr.12c.p.1p32.romea.edai ] , this momentum shift leads to an asymmetry between the positive and negative missing - momentum side of the momentum distribution .
note that positive missing momentum corresponds to @xmath208 .
c@xmath6 momentum distribution for the @xmath209 state as a function of the missing momentum .
the solid ( dashed ) curve are romea ( rpwia ) calculations.,scaledwidth=50.0% ]      finally , we present some results for the @xmath5he@xmath6 reaction at an incident proton energy of @xmath7  mev .
figure  [ fig : diffcross.4he.p.vanoers.250 ] compares the data from the triumf experiment @xcite with romea calculations using the optical potential of van oers _ et al .
_ @xcite . the typical shape for knockout of an @xmath85-state proton is reproduced by the romea predictions .
this fair agreement of the romea results with the data demonstrates that our romea model also works satisfactorily at lower incident energies .
a relativistic and cross - section factorized framework to describe the ifsi in quasielastic @xmath0 reactions has been outlined .
the model , which relies on the eikonal approximation , can use either optical potentials or glauber multiple - scattering theory to deal with ifsi .
thanks to the freedom of choice between these two substantially different techniques , our model is expected to be applicable at both intermediate and high incident energies .
the properties of the ifsi factor have been investigated for an incident proton energy of @xmath1  gev .
not surprisingly , the strongest attenuation occurs in the nuclear interior and heavier target nuclei are found to induce larger ifsi effects .
also , the surface - peaked character of the @xmath0 reaction was clearly established to be a consequence of the ifsi . whereas the different types of optical - potential sets contained in ref .
@xcite yield comparable ifsi factors , the rmsga calculations exhibit an unrealistic behavior for these kinematics .
the romea model has been used to calculate cross sections for the kinematics of two different experiments  : quasielastic proton scattering from @xmath2c , @xmath3o , and @xmath4ca at @xmath1  gev , and @xmath5he@xmath6 scattering at @xmath7  mev .
the predictions are shown to reproduce the shape of the data reasonably well at both incoming energies , thereby providing support for the wide applicability range of our model .
although the rmsga approach was deemed unsuitable for the kinematics discussed here , it should prove useful when trying to describe nuclear transparency data .
work in this direction is in progress .
this work was supported by the fund for scientific research , flanders ( fwo ) and the research council of ghent university .
for spherically symmetric potentials , the solutions @xmath210 to a single - particle dirac equation have the form @xcite @xmath211 \ ; , \label{eq : bound_state}\ ] ] where @xmath58 denotes the principal , @xmath212 and @xmath57 the generalized angular momentum quantum numbers .
the @xmath213 are the spin spherical harmonics and determine the angular and spin parts of the wave function , @xmath214 with @xmath215    the fourier transform of the relativistic bound - nucleon wave function is given by @xmath216 \ ; , \label{eq : fourier_bound_wave}\ ] ] with @xmath84 . the radial functions @xmath82 and @xmath83 in momentum space are obtained from their counterparts in coordinate space  :    @xmath217    @xmath218    with @xmath219 the spherical bessel function of the first kind .
the differential @xmath0 cross section ( [ eq : factorized_ifsi_cross_sect ] ) is proportional to the distorted momentum distribution @xmath110 of eq .
( [ eq : dist_mom_dist ] ) . when approximating the completeness relation ( [ eq : completeness ] ) as @xmath220 i.e. , neglecting the negative - energy term as in sec .
[ sec : rpwia ] , this amounts to @xmath221    thus , with @xmath222 defined as @xmath223 the function @xmath224 \label{eq : delta_r } \\ & = & \sum_{m } \bigl ( \overline{d } \left ( r_{1 } \right ) \ : \phi_{\alpha_1}^{d } + \overline{\phi_{\alpha_1}^{d } } \ : d \left ( r_{1 } \right ) \bigr ) \nonumber\end{aligned}\ ] ] represents the contribution of an infinitesimal interval in @xmath171 around @xmath225 to the @xmath0 cross section .
this procedure also enables us to estimate the average density seen through this reaction as @xmath226      g. jacob and th .
j. maris , rev .
mod . phys . * 38 * , 121 ( 1966 ) .
g. jacob and th .
j. maris , rev .
* 45 * , 6 ( 1973 ) .
p. kitching , w. j. mcdonald , th .
j. maris , and c. a. z. vasconcellos , adv .
. phys . * 15 * , 43 ( 1985 ) . n. s. chant and p. g. roos , phys .
c * 15 * , 57 ( 1977 ) .
s. kullander , f. lemeilleur , p. u. renberg , g. landaud , j. yonnet , b. fagerstrm , a. johansson , and g. tibell , nucl . phys . *
a173 * , 357 ( 1971 ) .
r. k. bhowmik , c. c. chang , p. g. roos , and h. d. holmgren , nucl .
a226 * , 365 ( 1974 ) .
r. bengtsson , t. berggren , and ch .
gustafsson , phys .
rep . * 41 * , 191 ( 1978 ) .
p. kitching , c. a. miller , w. c. olsen , d. a. hutcheon , w. j. mcdonald , and a. w. stetz , nucl . phys . * a340 * , 423 ( 1980 ) . w. t. h. van oers , b. t. murdoch , b. k. s. koene , d. k. hasell , r. abegg , d. j. margaziotis , m. b. epstein , g. a. moss , l. g. greeniaus , j. m. greben , j. m. cameron , j. g. rogers , and a. w. stetz , phys . rev .
c * 25 * , 390 ( 1982 ) .
d. s. carman , l. c. bland , n. s. chant , t. gu , g. m. huber , j. huffman , a. klyachko , b. c. markham , p. g. roos , p. schwandt , and k. solberg , phys .
b452 * , 8 ( 1999 ) . n. s. chant and p. g. roos , phys .
c * 27 * , 1060 ( 1983 ) .
r. a. arndt , i. i. strakovsky , and r. l. workman , phys .
c * 62 * , 034005 ( 2000 ) ; scattering analysis interactive dial - in program ( said ) , http://gwdac.phys.gwu.edu/. j. a. mcneil , l. ray , and s. j. wallace , phys .
c * 27 * , 2123 ( 1983 ) ; j. a. mcneil , j. r. shepard , and s. j. wallace , phys .
lett . * 50 * , 1439 ( 1983 ) ; j. r. shepard , j. a. mcneil , and s. j. wallace , _ ibid .
_ * 50 * , 1443 ( 1983 ) . c. j. horowitz , phys .
c * 31 * , 1340 ( 1985 ) ; c. j. horowitz and m. j. iqbal , _ ibid .
_ * 33 * , 2059 ( 1986 ) ; d. p. murdock and c. j. horowitz , _ ibid . _
* 35 * , 1442 ( 1987 ) ; c. j. horowitz and d. p. murdock , _ ibid . _
* 37 * , 2032 ( 1988 ) .
o. v. maxwell , nucl .
a600 * , 509 ( 1996 ) ; * a638 * , 747 ( 1998 ) ; o. v. maxwell and e. d. cooper , _ ibid . _ * a656 * , 231 ( 1999 ) .
e. d. cooper and o. v. maxwell , nucl .
a493 * , 468 ( 1989 ) ; o. v. maxwell and e. d. cooper , _ ibid . _
* a513 * , 584 ( 1990 ) ; * a565 * , 740 ( 1993 ) ; * a574 * , 819 ( 1994 ) ; * a603 * , 441 ( 1996 ) . y. ikebata , phys .
c * 52 * , 890 ( 1995 ) .
j. mano and y. kudo , prog .
theor . phys . *
100 * , 91 ( 1998 ) .
c. a. miller , k. h. hicks , r. abegg , m. ahmad , n. s. chant , d. frekers , p. w. green , l. g. greeniaus , d. a. hutcheon , p. kitching , d. j. mack , w. j. mcdonald , w. c. olsen , r. schubank , p. g. roos , and y. ye , phys .
c * 57 * , 1756 ( 1998 ) .
r. neveling , a. a. cowley , g. f. steyn , s. v. frtsch , g. c. hillhouse , j. mano , and s. m. wyngaardt , phys .
c * 66 * , 034602 ( 2002 ) .
d. l. adams and m. bleszynski , phys .
* b136 * , 10 ( 1984 ) . j. a. tjon and s. j. wallace , phys .
c * 32 * , 1667 ( 1985 ) .
b. i. s. van der ventel and g. c. hillhouse , phys .
c * 69 * , 024618 ( 2004 ) .
d. f. jackson and t. berggren , nucl .
62 * , 353 ( 1965 ) .
a. nadasen , p. schwandt , p. p. singh , w. w. jacobs , a. d. bacher , p. t. debevec , m. d. kaitchuck , and j. t. meek , phys .
c * 23 * , 1023 ( 1981 ) .
p. schwandt , h. o. meyer , w. w. jacobs , a. d. bacher , s. e. vigdor , m. d. kaitchuck , and t. r. donoghue , phys .
c * 26 * , 55 ( 1982 ) .
d. g. madland , los alamos national laboratory report la - ur-87 - 3382 ( unpublished ) .
e. d. cooper , b. c. clark , r. kozack , s. shim , s. hama , j. i. johansson , h. s. sherif , r. l. mercer , and b. d. serot , phys . rev .
c * 36 * , r2170 ( 1987 ) .
e. d. cooper , b. c. clark , s. hama , and r. l. mercer , phys .
b206 * , 588 ( 1988 ) . s. hama , b. c. clark , e. d. cooper , h. s. sherif , and r. l. mercer , phys .
c * 41 * , 2737 ( 1990 ) .
e. d. cooper , s. hama , b. c. clark , and r. l. mercer , phys .
c * 47 * , 297 ( 1993 ) .
l. antonuk , p. kitching , c. a. miller , d. a. hutcheon , w. j. mcdonald , g. c. neilson , and w. c. olsen , nucl . phys . *
a370 * , 389 ( 1981 ) . c. samanta , n. s. chant , p. g. roos , a. nadasen , j. wesick , and a. a. cowley , phys . rev .
c * 34 * , 1610 ( 1986 ) .
a. a. cowley , j. j. lawrie , g. c. hillhouse , d. m. whittal , s. v. frtsch , j. v. pilcher , f. d. smit , and p. g. roos , phys .
c * 44 * , 329 ( 1991 ) .
a. a. cowley , g. j. arendse , j. a. stander , and w. a. richter , phys . lett . * b359 * , 300 ( 1995 ) .
a. a. cowley , g. j. arendse , r. f. visser , g. f. steyn , s. v. frtsch , j. j. lawrie , j. v. pilcher , t. noro , t. baba , k. hatanaka , m. kawabata , n. matsuoka , y. mizuno , m. nomachi , k. takahisa , k. tamura , y. yuasa , h. sakaguchi , t. itoh , h. takeda , and y. watanabe , phys .
c * 57 * , 3185 ( 1998 ) .
g. p. mccauley and g. e. brown , proc .
london * 71 * , 893 ( 1958 ) .
c. j. joachain , _ quantum collision theory _
( elsevier , amsterdam , 1975 ) .
r. j. glauber and g. matthiae , nucl .
b21 * , 135 ( 1970 ) .
s. j. wallace , phys .
c * 12 * , 179 ( 1975 ) .
d. r. yennie , in _ hadronic interactions of electrons and photons _ , edited by j. cummings and d. osborn ( academic press , new york , 1971 ) , p. 321 .
s. l. belostotsky , yu . v. dotsenko , n. p. kuropatkin , o. v. miklukho , v. n. nikulin , o. e. prokofiev , yu .
a. scheglov , v. e. starodubsky , a. yu .
tsaregorodtsev , a. a. vorobyov , and m. b. zhalov , in _ proceedings of the international symposium on modern developments in nuclear physics _
( novosibirsk , 1987 ) , p. 191 . j. d. bjorken and s. d. drell , _ relativistic quantum mechanics _ ( mcgraw - hill , new york , 1964 ) .
j. a. caballero , t. w. donnelly , e. moya de guerra , and j. m. udas , nucl . phys . *
a632 * , 323 ( 1998 ) .
j. r. vignote , m. c. martnez , j. a. caballero , e. moya de guerra , and j. m. udas , phys .
c * 70 * , 044608 ( 2004 ) .
r. d. amado , j. piekarewicz , d. a. sparrow , and j. a. mcneil , phys .
c * 28 * , 1663 ( 1983 ) .
d. debruyne , j. ryckebusch , w. van nespen , and s. janssen , phys .
c * 62 * , 024611 ( 2000 ) .
d. debruyne , j. ryckebusch , s. janssen , and t. van cauteren , phys . lett . *
b527 * , 62 ( 2002 ) .
j. ryckebusch , d. debruyne , p. lava , s. janssen , b. van overmeire , and t. van cauteren , nucl . phys . *
a728 * , 226 ( 2003 ) .
k. hagiwara _
et al . _ ,
d * 66 * , 010001 ( 2002 ) , http://pdg.lbl.gov . k. hatanaka , m. kawabata , n. matsuoka , y. mizuno , s. morinobu , m. nakamura , t. noro , a. okihana , k. sagara , k. takahisa , h. takeda , k. tamura , m. tanaka , s. toyama , h. yamazaki , and y.
yuasa , phys .
78 * , 1014 ( 1997 ) .
j. d. walecka , _ electron scattering for nuclear and nucleon structure _ ( cambridge university press , cambridge , 2001 ) .