in 1974 , t hooft described the behavior of the large @xmath1 limit of qcd in @xmath0 dimensions @xcite .
he showed that the model confined quarks and was exactly soluable .
his work was followed quickly by several papers @xcite discussing various aspects of the model .
the work by einhorn @xcite discussed the t hooft model predictions for deep inelastic scattering ( dis ) .
einhorn found that the dis cross section was proportional to the square of the undressed quark masses , and hence would approach zero as the quark masses approach zero .
this pathlogical result made it difficult to apply the t hooft model to the phenomenology of dis ; clearly someting major was missing .
of course , limiting the physics to @xmath0 dimensions also excludes transverse spin degrees of freedom , so @xmath3 was zero and the callan - gross relation @xcite did not hold .
since then qcd in @xmath0 diminsions has been studied from a variety of points of view @xcite and applied recently to the study of duality in heavy meson decays @xcite .
however , to the best of our knowledge , no one has followed up on einhorn s study @xcite of dis .
a few years ago we found that a very simple toy model of scalar and spinor fields in @xmath0 dimension gave nice phenomenological results for dis @xcite . by extending the spinor degrees of freedom to 1 + 3 dimensions we were able to recover the callan - gross relation .
we decided to extend our approach to the t hooft model , and learned of the pathology in its description of dis .
it turns out that this pathology can be removed by extending the spinor degrees of freedom to 1 + 3 dimensions , just as we did before , and this paper grew out of that study .
this paper is divided into five sections , with seven appendices that include many of the details . following this brief introduction , sec .
ii reviews the t hooft model results for the strong sector .
our discussion focuses on ( i ) a method of treatment of the confining interaction that removes all singularities from the theory and gives finite dressed quark masses , and ( ii ) the consequence of including transverse degrees of freedom for the spin of the quarks .
we calculate the dressed quark mass , and the properties and spectrum of @xmath4 bound states , and also study the consequences of the completness of the bound states .
we show how to construct an off - shell @xmath4 scattering matrix , and show that the on - shell scattering matrix must be zero , as required for confined particles .
sections iii and iv study electromagnetic interactions and deep inelastic scatterng . here
the presence of transverse spin degrees of freedom completly alters the einhorn discussion , giving us phenomonologically useful results .
we conclude with a brief discussion .    in a subsequent paper @xcite we plan to present a numerical study of the approach to scaling in the dis limit .
we propose that the action have the following form @xmath5 where the qcd part of the action ( gluons and quark momenta or coordinates ) live in @xmath0 dimensions only , but the quark spin and the qed part ( leptons and the photon ) live in the full @xmath2 dimensions @xmath6 in his original paper , t hooft @xcite discussed the strong interactions only , and in this section we review the 1 + 1 dimensional model of qcd that he presented .
the qed part of the action ( [ a2 ] ) , and reduced qcd ( rqcd ) , will be discussed in sec .
[ rqcd ] .
the qcd lagrangian density is @xmath7 + \sum_i \bar{q}_i\left(i d_{\mu } \gamma^{\mu}-m_{0i } \right)q_i\ , , \label{1eq1}\ ] ] where @xmath8 is the quark field with flavor @xmath9 and bare mass @xmath10 .
the gluon field quantities are @xmath11 \\
\nonumber d_{\mu}&=&\partial_{\mu}+ig_0 a_{\mu } , \label{1eq2}\end{aligned}\ ] ] where @xmath12 is the gluon field tensor and @xmath13 are the gluon fields with the lorentz index @xmath14 and the color index @xmath15 .
the gluon fields have components in the 0 and 3 direction only , so the sum over the index @xmath14 is restricted to 0 and 3 .
the qcd coupling constant is @xmath16 and @xmath17 are the generators of the @xmath18 color group , normalized to @xmath19 = 2\,\delta^{ab}\nonumber\\ & & \sum_{a}\lambda^a\lambda^a= \frac{2(n_c^2 - 1)}{n_c}\;\openone \ , , \label{1eq2aa}\end{aligned}\ ] ] where matrix multiplication of the @xmath20 matrices is implied . as @xmath21 , finite results are obtained if @xmath22 as @xmath23 where the effective coupling @xmath24 is a constant .
while the space - time coordinates , momenta , and lorentz vector sums in eq .
( [ 1eq1 ] ) are restricted to 1 + 1 dimension , we assume that the gamma matrices have the usual @xmath25 dirac structure .
the quark fields are therefore a direct product of four dirac dimensions and @xmath1 color dimensions .
note that the action ( [ a2 ] ) is lorentz invariant under the subgroup @xmath26 of lorentz transformations that leave the @xmath27 plane invariant .
specifically , the group is generated by the hamiltoninian @xmath28 , and the boost @xmath29 , momentum @xmath30 , and angular momentum @xmath31 , operators that generate boosts and translations along the @xmath32 axis , and rotations about the @xmath32 axis .
four of the six commutation relations between these generators are zero , and the other two close @xmath33=-i{\cal h } \qquad   [ { \cal k}_3,{\cal h}]=-i{\cal p}_3   \ , , \label{lgroup}\end{aligned}\ ] ] guaranteeing that @xmath26 is a group .
following t hooft @xcite we introduce light cone variables : @xmath34 for any arbitrary vector @xmath35 ( in the qcd sector , the only perpendicular components come from matrix elements of @xmath36 ) .
note that the scalar product of any two vectors @xmath15 and @xmath35 is @xmath37 the derivatives are defined @xmath38 so that the divergence of a two - vector is @xmath39 in the same way , we can define the @xmath40 , @xmath41 , and @xmath42 components of the @xmath43 matrices .
the anticommutation relations are all zero except for @xmath44    since the gluon fields are confined to 1 + 1 dimensions , there is only one nonvanishing component of the gluon field strength tensor @xmath45\ , .\ ] ] the qcd part of the theory is simplified if we choose the light cone gauge , where @xmath46 , so that the commutator contained in the field tensor @xmath47 disappears and there is only one nonzero component ( @xmath48 ) of the gluon field . the lagrangian density ( [ 1eq1 ] ) then reduces to @xmath49 + \sum_i \bar{q}_i\left(i
\partial_+ \gamma_- + i \partial_- \gamma_+ -g_0 \gamma_{-}a_{+ }   -m_{0i }   \right)q_i\ , .
\label{1eq4}\ ] ] the equation of motion for this field is then @xmath50 the solution of ( [ 1eq5 ] ) is @xmath51 where the green s function @xmath52 is @xmath53 the coefficients @xmath54 and @xmath55 can not be determined without knowing the boundary conditions , so they are free parameters .
the gauge condition did not eliminate all superfluous degrees of freedom , just as the coulomb gauge , or the lorentz gauge , do not determine uniquely the photon propagator in qed ( gribov ambiguity ) .
we can therefore set the coefficients @xmath54 and @xmath55 equal to zero ( a specific choice of gauge ) in order to simplify our calculations .
einhorn discussed the gauge issues related to these two parameters and showed that the eigenvalues of the two body bound state equation are independent of the choice of @xmath54 and @xmath55 .
however , the dressed quark mass does depend on the choice of @xmath55 , and this in turn implies that the location of the quark mass pole is gauge dependent .
since any physically meaningful quantity is gauge invariant , we can conclude that the location of the mass pole is not physically meaningful , and this can happen only if the quarks are confined , so that free quark states do not exist .
this is an indirect consequence of confinement .
direct consequences of confinement will be discussed in the next subsection when the two - body bound state equation is discussed .
equation  ( [ 1eq6 ] ) shows that the gluonic field is no longer a dynamical variable , and that there are therefore no ghosts .
the fourier transform of the green s function ( [ 1eq7 ] ) gives us the gluon `` propagator '' , or more precisely the momentum dependence of the effective quark - quark interaction . with @xmath54 and @xmath55 equal to zero ,
it can be written @xmath56 this is a singular operator that is well defined only when it is part of an integral over @xmath57 . the second term in eq . ( [ 1eq8 ] )
was introduced by gross and milana @xcite in a different context .
its purpose is to preserve the condition @xmath58 and it also insures that the potential is finite for any finite value of @xmath59 .
the only variables remaining are the quark fields .
the feynman rule for the ( undressed ) quark propagator is @xmath60 and the quark-``gluon '' coupling is @xmath61 leading to the following result for the exchange of gluons with momentum @xmath57 between quarks with dirac indicies 1 and 2 @xmath62\nonumber\\ & & =    -ig_0 ^ 2\ , \frac{2(n_c^2 - 1)}{n_c } \gamma^1_{-}{\cal o}_{12}\gamma^2_-\ , \tilde{\cal g}(k_- ) % & & { \displaystyle\lim_{n_c\to\infty } } \to -ig^2 \gamma^1_{-}{\cal o}_{12}\gamma^2_-\ , \tilde{\cal g}(k_- ) \
, , \label{1eq11}\end{aligned}\ ] ] where we used eq .
( [ 1eq2aa ] ) ( with matrix multiplication of the @xmath20 matrices implied ) , and took the limit as @xmath21 using the definition ( [ 1eq2a ] ) .
the quark self - energy and the dressed quark mass are calculated in the next subsection .
we determine the dressed single quark propagator , @xmath63 , using the ( one body ) dyson schwinger equation ( dse ) : @xmath64s_{0}(p)\ , , \label{dse}\end{aligned}\ ] ] shown graphically in fig .
[ ds ] . the gluon interaction does not mix quark flavors , and the calculation is identical for each flavor of quark , so the flavor index is suppressed . the rainbow approximation ( undressed vertices and the absence of the quark loops from the gluon propagator ) is justified in the large @xmath1 limit @xcite . since for every internal loop
there is a factor of @xmath65 , and a multiplicative factor of @xmath66 , the color dependence disappears .
the vertex corrections and the quark - gluon vertices do not have a multiplicative factor of @xmath1 , and are therefore supressed in the large @xmath1 limit .    in eq .
( [ dse ] ) , @xmath67 , and since @xmath68 does not depend on @xmath69 , it follows immediately that the self - energy integral does not depend on @xmath70 either , and must have the form @xmath71 .
hence the dressed propagator is of the form
@xmath72 \gamma_- -i\epsilon } \ , , \label{1eq11a}\ ] ] where , using @xmath73 , the self - energy contribution is then @xmath74 performing the @xmath69 integral gives @xmath75 and substituting this back into eq .
( [ 1eq12 ] ) gives @xmath76 using ( [ 1eq8 ] ) for @xmath77 gives @xmath78 substituting this back into eq .
( [ 1eq11 ] ) gives @xmath79 \gamma_- -i\epsilon } = \frac{m_0+p_-\gamma_+ + \left(p_+ + { \displaystyle\frac{g^2}{2\pi p_- } } \right ) \gamma_- }   { m_0 ^ 2-{\displaystyle\frac{g^2}{\pi}}-2 p_+p_-   -i \epsilon}\ , .
\label{1eq16}\end{aligned}\ ] ] note that the quark mass pole has been shifted to a _ smaller _ value @xmath80 where @xmath81 is the dressed quark mass .
it turns out that choosing @xmath82 will not effect the dressing of the mass , but choosing @xmath83 would give @xmath84 this gives a mass that is momentum dependent and not covariant .
hence the dressed mass is gauge dependent and unphysical . with these cautionary remarks ,
we choose @xmath85 because the mass is covariant .
one undesirable feature of this choice is that @xmath86 in the chiral limit .
this could be corrected by choosing @xmath55 large and positive .
having obtained the dressed propagator , we are able to proceed with the two body bound state calculation .        consider a bound state of a @xmath4 pair .
the quark has dressed mass @xmath87 and electric charge @xmath88 , and the antiquark ( which might be of a different flavor ) has dressed mass @xmath89 and charge @xmath90 .
the momentum of the bound state is @xmath91 , the momentum of the quark is @xmath92 and the momentum of the outgoing antiquark is @xmath93 , as shown in fig .
[ bs ] . in much of the following discussion
, we will treat the outgoing antiquark as an incoming quark with momentum @xmath94 , and following this convention label the bound state vertex function @xmath95 . for a color neutral state we may carry out the color sums using ( [ 1eq11 ] ) , giving the bethe - salpeter equation for the bound state vertex function : @xmath96 where @xmath97 and @xmath98 are dressed quark propagators , and
the kernel is the singular operator @xmath99 with the substitution @xmath100 @xcite , eq .
( [ bse ] ) becomes : @xmath101 \left[m_2 ^ 2-(k - r)^2 -i\epsilon\right]}\ , .
\label{bse2}\ ] ] the formula shows that @xmath102 does not depend on @xmath70 , and hence the @xmath69 integration can be carried out immediately .
integrations of the variable @xmath69 over two quark propagators occurr frequently , and are worked out in detail in appendix [ appen : a ] .
it is convenient to introduce the momentum fractions @xmath103 if the fraction @xmath32 lies in the interval [ 0,1 ] , eq .
( [ bse2 ] ) reduces to @xmath104\right ) \
, , \label{phiin}\end{aligned}\ ] ] where @xmath105 is the principal value integral , and @xmath106 , @xmath107 , and @xmath108 are dimensionless parameters @xmath109 with @xmath110 the anticipated mass eigenvalues that will emerge from the solution of the equation , and @xmath111 the _ dressed _ quark masses . the two - body `` wave function '' @xmath112 , is _ defined _ by @xmath113 where @xmath114 is the anticipated eigenfunction solution for the vertex corresponding to the @xmath115th bound state , and @xmath116 is a normalization constant .
it is chosen so that @xmath117 and will be calculated in subsection [ elasticff ] .
note that eq .
( [ phiin ] ) _ defines _ the two - quark hamiltonian , @xmath118 , on the interval @xmath119 .
the _ bare _
quark masses enter into the hamiltonian , while the _ dressed _ quark masses enter into the relation ( [ 2eq3d ] ) between the vertex function and the wave function .    if the momentum fraction @xmath32 lies _ outside _ of the interval @xmath119 , the vertex function defined by eq .
( [ bse2 ] ) is not zero , but can be obtained by quadrature from the vertex function defined _ inside _ of the interval [ 0,1 ] . using ( [ 2eqac ] ) and
the definition ( [ 2eq3d ] ) gives @xmath120\right ) \
, .\qquad \label{phiout}\end{aligned}\ ] ] these definitions and equations will be used in the following sections when the hadronic currents are computed .
an alternative form of eq .
( [ phiin ] ) follows if we integrate the last term by parts : @xmath121 where @xmath122 .
if @xmath123 , which is true in all but the chiral limit , this form of the equation follows directly from t hooft s principle value prescription @xcite , so that our two - body equation is identical with t hooft s , even though we have finite dressed quark masses .    to avoid confusion later , we call attention to the fact that the wave function defined by eqs .
( [ 2eq3d ] ) and ( [ phiout ] ) is not identical to the bethe - salpeter wave function .
the latter can be defined @xmath124 where the factor @xmath125 in the numerator comes from the @xmath126 factors in the quark propagators @xmath127 .
this wave function carries the @xmath70 dependence of the quark propagators . if we integrate ( [ bswave ] ) over @xmath70 ( using the fact that @xmath128 does not depend on @xmath70 ) we obtain @xmath129 $ \cr 0 & if $ z\notin [ 0,1]$ }   \ , .
\label{equalpp}\ ] ] this is the `` equal @xmath130 '' wave function , and equals the wave function we are using _ only in the region _
@xmath131 $ ] . as we have seen , the wave function @xmath132 that solves the t hooft equation is non zero _ for all _ @xmath32 .
we conclude this discussion of the two body wave equation by returning to eq .
( [ bse ] ) . in place of the replacement @xmath133 , it could equally well have been assumed that @xmath134 where @xmath135 is either @xmath136 or @xmath137 ( recall that we are assuming that the dirac matrices live in 1 + 3 dimensional space ) .
note that the bethe - salpeter equations with the substitutions ( [ 2eqeee ] ) reduce to the same equation ( [ bse2 ] ) for @xmath138 , and hence _ each two - body state is triplely degererate , with spin structure given by @xmath126 , @xmath139 , or @xmath140_. in each of these three cases the wave functions and bound state mass are identical .
the additional states with the @xmath141 or @xmath142 structure are outside the scope of the original t hooft model , but must be considered in reduced qcd .
they are states with quark spins in the transverse direction , and will play a major role in the discussion of dis below .
the chiral limit is exactly soluable and of some interest . if @xmath143 , then eq .  ( [ phiin ] ) has the normalized ground state solution @xmath144 outside of the interval @xmath145 $ ] , eq .
( [ phiout ] ) gives the same result , so that @xmath146 for all @xmath32 in @xmath147 .
the t hooft model illustrates how confinement can be realized in two ( apparently different ) ways .
if the singularities of the confining interaction are left unregularized , as they were in t hooft s original paper @xcite , the quark masses are infinite . even if the singularity is regularized as in this paper , whenever @xmath148 the dressed quark mass will be unphysical , and there is no real quark mass pole .
this insures that the quarks are confined .    in the second case (
@xmath149 and @xmath150 ) the dressed quark mass is real and positive and the there is a quark mass pole .
what keeps the quarks confined in this case ?
the answer is that the quark vertex function @xmath151 has a _ zero _ at precisely the value of @xmath32 where all of the quarks in the bound state could be on mass shell .
this requirement is illustrated diagramatically in fig .
[ confinement ] , and seems miraculous .
however , as discussed in ref .
@xcite , it is only a consequence of fact that the _ wave function _ is finite at the pole ( because the wave equation permits no singularities ) , and therefore the vertex function @xmath152 \frac{\phi(z , r)}{{\cal n } r_-}\end{aligned}\ ] ] must be zero .
precisely the same mechanism works in the schrdinger description of two nonrelativistic particles bound by a linear ( or any other ) confining potential .    in summary , quark mass poles
are always unphysical because the dressed mass is gauge dependent , and bound state poles are physical because they are gauge invariant . in this respect
the undressed quark mass is the only physically significant parameter .
the two ways of viewing confinement are therefore not exclusive , but equivalent , since they are related through gauge transformations .
the bound state equation was solved numerically using modified cubic splines .
the splines and the technique are fully described in appendix [ appen : splines ] .
each spline is a smooth function ( continuous up to its third derivative ) and has support over a small region of the @xmath32 axis .
the splines overlap , so that by independently adjusting the height of each spline it is possible to smoothly approximate any function defined on the interval @xmath145 $ ] .
.numerical values of @xmath153 for @xmath154 and various @xmath155 ( number of splines ) and @xmath156 ( number of guass points for each spline segment ) .
[ cols=">,>,^,^,^,^,^,^,^,^,^,^,^,^ " , ]     [ tab:1 ]    we have modified the standard spline technique to insure that the numerical solutions satisfy the boundary conditions exactly .
these boundary conditions are determined by examing eq .
( [ phiin ] ) in the vicinity of @xmath157 and @xmath158 @xcite , and require the wave function go like @xmath159 note that the @xmath160 in the chiral limit ( when @xmath161 ) , consistent with the exact solution ( [ chiralsol ] ) . to fit these boundary conditions the splines closest to @xmath162 and @xmath163 have these fractional powers built into their functional form , as described in appendix [ appen : splines ] .
the wave function is expanded into @xmath155 splines , each with four segments ( except for those at the boundary , which have only three ) .
the splines overlap , and taking into account that there must be at least three splines ( one regular spline in the center and one modifed spline at each end ) , the total number of segments is @xmath164 , dividing the interval into segments of length @xmath165 . the equation is turned into a matrix equation by integrating over each segment using @xmath156 gaussian points .
tables [ tab:2 ] and [ tab:1 ] show the numerical stability .
if the quark mass is larger than unity , better than 1% accuracy is achieved with only two gauss points per interval , provided @xmath166 . as the level number @xmath115 increases above @xmath167 toward @xmath155 , the reliability of the calculation decreases . to study highly excited states , it is sufficient to use @xmath168 , and push @xmath155 as high as possible .
the situation is somewhat different if the bare quark masses are less than unity ( table [ tab:1 ] ) . in this case
the wave function is very steep at the boundaries and 4 or 6 gaussian points per interval are needed to get 1% accuracy for the low lying states . in a subsequent paper @xcite , we will limit our numerical discussion of duality to cases with @xmath169 where the solutions are very stable , and there are no poles in the region outside of [ 0,1 ] [ see the brief remarks following eq.([rlimit ] ) ] .
figure [ mass ] shows the linearly rising regge trajectory characteristic of the t hooft model .
figures [ 1st]-[10th ] show the ground , first excited and 10th level for three cases with equal quark masses of different values .
note the `` tails '' of the wave functions outside of the region [ 0,1 ] , given by eq .
( [ phiout ] ) .
these tails are large for small quark masses , but shrink quickly as the quark mass increases , and are also small for highly excited states .
figure [ 20th ] shows the 20th level for spline numbers 20 , 40 , and 60 .
note that the state is not fully described even with 40 splines , supporting the observation that an accurate description of states of order @xmath115 requires a spline number of approximately @xmath170 .
the hamiltonian defined in eq .
( [ phiin ] ) is hermitian on the interval @xmath145 $ ] and it is straightforward to show that the different solutions are orthogonal on this interval . with the normalization given in eq .
( [ norm ] ) @xmath171 the completeness of the eigenfunctions implies that @xmath172\ , .
\label{complete}\ ] ] the greens function is therefore @xmath173 where @xmath174 is any external momentum .
the completness relation shows that @xmath128 has the property of a greens function @xmath175 \ , .\end{aligned}\ ] ] the same argument also gives @xmath176 \ , .\end{aligned}\ ] ]      in preparation for the study of electromagnetic interactions , it is useful to introduce the @xmath4 scattering matrix , defined by the following infinite series @xmath177_1 \otimes\left[\gamma_{-}s_2(k - q)\gamma_-\right]_2 \nonumber\\ & & \qquad+i^2\int\frac{d^2k'}{(2\pi)^2}\int\frac{d^2k}{(2\pi)^2 } v(p',k')\,v(k',k)\,v(k , p)\,\nonumber\\ & & \qquad\qquad\times \left[\gamma_-s_1(k')\gamma_{-}s_1(k)\gamma_-\right]_1 \otimes\left[\gamma_{-}s_2(k - q)\gamma_-s_2(k'-q)\gamma_-\right]_2   \nonumber\\ & & \qquad + \cdots \ , .\quad \label{eq248}\end{aligned}\ ] ] the momentum are defined as for the two - body bound state , except that here @xmath174 is the unconstrained momentum of the @xmath4 pair and @xmath178 is not an eigenvalue .
the kernel @xmath179 is the singular operator ( [ kernel2 ] ) .
the singularities of this operator will never be exposed because we will limit application of this scattering equation to cases where we integrate over either the initial or final momentum .
note that the structure of the equation anticipates the reduction of the dressed propagators @xmath180_i=\frac{k_-\ , [ \gamma_-\gamma_+\gamma_-]_i } { m_i^2 - 2k_-k_+}=\frac{2k_-}{m_i^2 - 2k_-k_+}\,\gamma_{i-}\ , , \end{aligned}\ ] ] so that the factor of @xmath181 is common , and the series ( [ eq248 ] ) reduces to @xmath182 where @xmath183 was defined in eq .
( [ d2 ] ) . from this result
we can conclude that @xmath184 does not depend on @xmath185 and @xmath70 .    in applications below , we will encounter this scattering series in the following form @xmath186 where @xmath187 is an operator that , by assumption , does not depend on @xmath70 . in appendix [ appen : b ] we show how to use the greens function to write this series as @xmath188 where we replaced @xmath189 and @xmath190 .
this equation has a nice physical interpretation .
it displays the `` scattering amplitude '' as the sum over the propagation of the bound states , which is all confinement will allow ( without meson decay mechanisms , which are ignored in this paper , no cuts are possible ) .
the only singularities in @xmath178 that can occur are poles at the bound state masses @xmath153 .    with these tools in place we are able to study the electromagnetic interaction of hadrons .
we return to the discussion of the lagrangian density introduced in eq .
( [ a2 ] ) and look at the electromagnetic sector . in rqcd the electromagnetic fields ( and the electron currents ) are extended to 1 + 3 dimensions .
the lagrangian density is therefore @xmath191 where @xmath192 is the electromagnetic field ( with the subscript 0 distinguishing it from the gluon field ) , the vector sums are now over 4 dimensions with @xmath193 , @xmath194 is the electromagnetic charge of quarks with flavor @xmath9 , @xmath195 is the electron current that produces the electromagnetic field , and @xmath196 is the electromagnetic field tensor . because the electron current exists in four dimensions , all four components of the electromagnetic field will be non - zero , in general , and all four components of the quark current will be excited by electromagnetic scattering .
we now apply this discussion to electron scattering , where @xmath197 , so that @xmath198 .
the initial particle ( quark or hadron ) has momentum @xmath199 and the final particle momentum @xmath200 , with momentum fractions @xmath201 if the final particle is on shell , then @xmath202 and @xmath203 , and @xmath204 as @xmath205 .
we will use the notation @xmath206 to denote the _ bare _ quark current _ operator _ , which is @xmath207    while the dressing of the strong quark - gluon vertex is supressed by the large @xmath1 limit @xcite , the quark electromagnetic vertex can be dressed by gluon exchanges .
the longitudinal part of the dressed electromagnetic vertex ( @xmath208 ) can be computed directly from the ward - takahashi identity @xmath209 where @xmath210 . using the dressed propagator ( [ 1eq16 ] ) , we find that @xmath211 consequently the @xmath212 component is unmodified .
it is instructive to obtain the full quark current directly from the dressing of the quark - photon vertex .
the vertex is dressed by successive gluon exchanges , as shown diagramatically in fig .
[ quarkcurrent ] . when coupling to a photon , there is no flavor change , and hence , for each flavor @xmath9 , this series is @xmath213 this series is evaluated in appendix [ appen : c ] using the methods developed in appendix [ appen : b ] .
the final results are conveniently expressed in term of a quark form factor , defined to be @xmath214 where @xmath215 and @xmath216 in anticipatation of applications to electron scattering .
note that the form factor is expressed as a sum of bound state poles at @xmath217 , arising from the propagation of the bound states that couple to the photon .
vector dominance is a rigorous consequence of this model .    using this form factor
, the currents reduce to @xmath218 note that the terms independent of the quark form factor @xmath219 are identical to the expected result ( [ current1 ] ) , and that the contribution from the form factor is _ purely transverse _ , in that @xmath220\big|_{f_i\;{\rm term } } & = & \gamma_- \,f_i(z , q^2)\,\left\ { q_+-q_-\frac{q_+}{q_-}\right\}=0\ , .\qquad \label{ff5}\end{aligned}\ ] ] this result mimics the method for insuring current conservation developed in ref .
@xcite , and gives some evidence that that method is dynamically sound .
these results are used in the next subsection to calculate the hadronic transition current .      in the presence of confinement , and without any meson decay mechanisms , the quarks in a @xmath4 bound state must remain bound , even after they absorb an energetic photon . in this section
we calculate the transition current for the process @xmath221 , where @xmath222 is the mass of a ( possibly excited ) final bound state and @xmath223 is the mass of the initial ( ground ) state .
the transition current consists of the sum of two contributions : one in which the photon is absorbed by the quark with dressed mass @xmath87 , and a second in which the photon is absorbed by the antiquark with mass @xmath89 .
the momentum used to label each of these processes are shown in fig .
[ hadronicamplitude ] ; the initial bound state has momentum @xmath224 and the final state has momentum @xmath225 .
the transition current coming from the first term is @xmath226\nonumber\\ & & \qquad\quad = i \int\frac{d^2k}{(2 \pi)^2}\frac { e_1n_c\,g_f\left(-p_f - k ,- k\right ) g_i\left(p+k , k \right)\,n^\mu } { [ m_1 ^ 2-(p_f+k)^2-i\epsilon ] [ m_1 ^ 2-(p+k)^2-i\epsilon ] [ m_2 ^ 2-k^2-i\epsilon ] } \ , , \qquad \label{3eq1-}\end{aligned}\ ] ] where @xmath227 is the quark current operator for quark 1 , as worked out in subsection  [ quarkem ] , @xmath228 is the spin dependent numerator discussed below , and we have allowed for the possibility that the final state may have _ either _ the @xmath126 _ or _
the @xmath139 structure discussed above ( we will not consider the @xmath140 states ; they are needed only if we consider the @xmath137 component of the transverse current and give results identical to the results for the @xmath136 components ) .
note the factor of @xmath1 coming from the sum over all possible quark colors .
only the term in the numerator of the @xmath98 propagator equal to @xmath229 can give a nonzero result , but the terms that contribute from the @xmath97 propagators depend on the matrix form of the final state , and of @xmath227 . because the trace of a product of gamma matrices with only one factor of @xmath141 or @xmath142 is zero , the transverse currents ( with @xmath141 or @xmath142 ) will be zero unless the final state also includes @xmath141 ( or @xmath142 for the @xmath137 component of the current ) .
hence we have the rules : ( i ) @xmath136-type transverse currents couple @xmath126 states to @xmath139 states only , and ( ii ) longitudinal currents couple @xmath126 states to each other .
because of the exact degeneracy of these states , the results in both of these cases use the same momentum - space wave functions .
the integral ( [ 3eq1- ] ) is evaluated in appendix [ appen : d ] .
the result for the @xmath41 component of the current is @xmath230 \ , , \qquad\qquad\qquad\qquad \label{3eq1b}\end{aligned}\ ] ] where @xmath231 [ for definitions of all the momentum fractions see eqs .
( [ frac2 ] ) and ( [ frac3 ] ) ] , the state @xmath232 has the same structure as the initial state ( assumed to be @xmath126 ) , @xmath233 is the quark form factor defined in eq .
( [ ff1 ] ) and @xmath234 is given in eq .
( [ rlimit ] ) .
the plus and transverse components of the current are @xmath235 } { m_{01}\,\eta\,(1-\eta)\,q^2 } \nonumber\\ \left < f_-|{\cal j}_+(p_f , p)|i_-\right > \big|_{e_1\,{\rm
term}}&=&\frac{16 e_1n_c\,p_{-}}{{(2 \pi)\cal n}^2 } \left\ { \int_{y}^1   +   \int_{0}^{y } { \cal r}\right\ } d\xi \;\phi_f(\xi',p_f)\phi_i(\xi , p)\nonumber\\ & & \qquad\qquad\times \frac{q_+\,\left[m_1 ^ 2- \eta(1-\eta)\,q^2 \,f_1(\eta , q^2)\right ] } { q_-\,\eta\,(1-\eta)\,q^2 } \ , , \qquad\qquad \label{currentd}\end{aligned}\ ] ] where the state @xmath236 has the structure @xmath139 .
these results are used in the next subsection , and in the discussion of dis .    using the definitions of the momenta shown in fig.[hadronicamplitude ]
, it is easy to see that the second term in the transition current can be obtained from the first by letting @xmath237 , as discussed in appendix [ appen : d ] . the exact expression for the @xmath136 component of the transition current , which is the sum of the @xmath88 and @xmath90 contributions ,
can therefore be written @xmath238 where the reduced current @xmath239 is @xmath240 with @xmath241 and @xmath242 and @xmath243 \ , .
\label{qffeff}\end{aligned}\ ] ] finally , the function @xmath244 , was defined in eq .
( [ rlimit ] ) ; the substitutions that convert the @xmath88 term into the @xmath90 term therefore give @xmath245 .
if the scattering is elastic , @xmath246 and the momentum fraction @xmath247 must be computed from the exact expressions for @xmath174 and @xmath224 in the c.m .
frame @xcite .
the result , for elastic scattering from the state @xmath115 , is @xmath248 as @xmath205 , @xmath249 and @xmath250 [ defined in eq .
( [ fraccur ] ) ] approaches @xmath251 , @xmath252 , and the _ exact _ result for the elastic curent becomes @xmath253 ^ 2=2\,e_1\,p_-\ , { \cal f}_n(0)\ , , \label{3eq1c}\end{aligned}\ ] ] where charge normalization requires that the form factor @xmath254 be unity at @xmath255 .
hence @xmath256 gives the normalization condition ( [ norm ] ) .
furthermore , elastic scattering requires that , as @xmath205 , @xmath257 and hence @xmath258 which follows from the wave function identity ( [ idet3 ] ) .
note that _ the charge is properly normalized and conserved for _ both _ components because we used the correct dressed quark current _ ( [ current1 ] ) .
finally , we point out that the elastic transverse current is zero for all @xmath259 , as required by the summetries of the states @xmath260 we now turn to a discussion of deep inelastic scattering .
using the standard conventions ( also defined in ref .
@xcite ) , the dis cross section is @xmath261\nonumber\\ & = & \sigma_m\left[w_2 + 2w_1\tan^2\left(\frac{\theta}{2}\right ) \right]\ , , \end{aligned}\ ] ] where the structure functions @xmath3 and @xmath262 are part of the hadronic tensor , defined by @xmath263 where @xmath264 is the target , or ground state mass .
the structure function @xmath3 is immediately extracted from the @xmath265 ( or @xmath266 ) component of the tensor @xmath267 here we anticipate that @xmath268 scales in the deep inelastic limit to the function @xmath269 of the variable @xmath270
. note that @xmath3 must be identically zero if the @xmath136 ( and @xmath137 ) components of the current are zero .
the structure function @xmath262 can be extracted from the @xmath40 and @xmath41 components of the current using @xmath271 in the deep inelastic limit this gives @xmath272 the contribution from @xmath273 will be shown to go like @xmath274 in the following sections .
hence , _ if @xmath233 is non - zero _ , we obtain the callan - gross relation .
however , if the transverse components of the current are omitted , @xmath275 and @xmath276 does not scale . in this case
it is the quanity @xmath277 that scales , and we recover the results of einhorn @xcite .
we now turn to a calculation of the structure functions .
we calculate them first ( i ) in the _ partonic picture _
, where all final state interations are ignored and the quarks in the final state are assumed to be free , and then ( ii ) in the _ hadronic _ picture , where confinement ensures that the only possible final states are the bound @xmath4 states we have already discussed .
we emphasize that the partonic picture can never actually occur because the confining interaction , which also acts in the final state , can never be ignored .
consider deep inelastic scattering ( dis ) in the partonic picture , which assumes a final state composed of free quarks with no interaction .
this picture is illustrated in fig .
[ partonicamplitude ] .
the full result is the sum of a contributions in which the photon is absorbed by the quark with charge @xmath88 , and one in which it is absorbed by the antiquark with charge @xmath90 . in the center of mass of the ejectiles ,
either particle can go forwards or backwards and the unpolarized cross section is the sum over both helicities of the particles .
begin by calculating the transverse component ( chosen be be @xmath136 ) of the current .
the outgoing quark has a mass of @xmath87 and momentum @xmath278 in the center of mass frame of the ejectiles , while the antiquark has a mass of @xmath89 and momentum of @xmath279 ( in this section @xmath279 will be the physical momentum of the outgoing antiquark , opposite in sign to that used in the earlier sections of this paper ) .
the current is the sum of two diagrams , only one of which is shown in fig .
[ partonicamplitude ] . since the @xmath4 vertex contains the factor of @xmath126 , only the factor @xmath280 in the numerator of the dressed propagator ( [ 1eq16 ] ) survives , but we will retain only the factor of @xmath281 which dominates the result at large @xmath259 .
the following result is obtained for the @xmath136 component of the current for the disintegration of the ground state ( the subscript @xmath282 usually attached to the ground state will be omitted in this section ) : @xmath283 \frac{(p_1-q)_-\,g(p_1-q ,-
p_2)}{m_1 ^ 2-(p_1-q)^2 } \nonumber\\ & & + e_2\,[{\bar u}(p_1)\,\gamma_-\gamma_+ \gamma_x\ , v(p_2 ) ] \frac{(p_2-q)_-\ , g(p_1,-p_2+q )   } { m_2 ^ 2-(p_2-q)^2 } \ , .
\label{5eq1}\end{aligned}\ ] ] introducing the momentum of the bound state , @xmath284 , and defining @xmath285 [ as in eq .
( [ defz ] ) ] , the first denominator can be written @xmath286(p - p_2)_- \nonumber\\ & = & m_1 ^ 2-z\,\left(m^2-\frac{m_2 ^ 2}{1-z}\right ) = z\left[\frac{m_1 ^ 2}{z}+\frac{m_2 ^ 2}{1-z}-m^2\right ] \ , .
\label{5eq1a}\end{aligned}\ ] ] here we have made use of the constraints @xmath287 and @xmath288 .
a similar result follows for the second denominator if we define a new momentum fraction @xmath289 @xmath290 with these definitions we obtain @xmath291 \;\frac { p_- g(z , p)}{\delta(z , m^2 ) } \nonumber\\ & & + e_2 [ { \bar u}(p_1)\ , \gamma_-\gamma_+ \gamma_x\ , v(p_2 ) ] \ ;       \frac { p_- g(1-z',p ) } { % { \displaystyle\frac{m_2 ^ 2}{z'}+\frac{m_1 ^ 2}{1-z'}- m^2 } \delta(1-z ' , m^2 ) }   \nonumber\\ \nonumber\\ & = & -\frac{e_1}{{\cal n } } \left[{\bar u}(p_1)\,\gamma_x\,\gamma_+\gamma_-\ , v(p_2 ) \right ] \phi(z , p ) \nonumber\\ & & + \frac{e_2}{\cal n}\left[{\bar u}(p_1)\ , \gamma_-\gamma_+ \gamma_x\ , v(p_2 ) \right]\phi(1-z',p ) \ , , \label{5eq1b}\end{aligned}\ ] ] where the answer has been expressed in terms of the wave function ( [ 2eq3d ] ) .
the current and the dis cross section are evaluated in the c.m .
frame , using the kinematics defined in appendix [ appen : e ] .
the unpolarized cross section is obtained by squaring the current and summing over spins using @xmath292\nonumber\\ & = & 32 \,p_{1+}p_{2- } \
, , \end{aligned}\ ] ] and neglecting the interference term ( which vanishes in the large @xmath259 limit ) .
observing that @xmath293 in the dis limit , and summing over all colors of the outgoing quarks ( which do not interfere ) gives @xmath294 where the momentum are evaluated in appendix [ appen : e ] .
the hadronic tensor defined in eq .
( [ wmn ] ) is related to the square of the currents by @xmath295 the @xmath296 component of the tensor is @xcite is therefore @xmath297 this gives the familiar parton model result for the structure function @xmath3 @xmath298    the structure function @xmath262 , given in eq .
( [ w2 ] ) , is then @xmath299 where we anticipate the result @xmath300 as @xmath301 .
this is the callan - gross relation .
now look at the plus and minus components of the current .
if the dirac space were restricted to two dimensions , only these components of the current would exist .
the minus component can be obtained immediately from ( [ 5eq1b ] ) by replacing @xmath141 with @xmath126 , giving @xmath302\,\frac{2}{\cal n}\,\left\{e_1 \phi(z , p)-e_2 \phi(1-z',p)\right\ } \
, , \label{5eq1d}\end{aligned}\ ] ] the plus component could include the two contributions from the quark current , eq .  ( [ current1 ] ) , but in the partonic picture the outgoing quark is not dressed , so the second part of the current proportional to @xmath303 should be ignored .
however , in order to study the contribution of this term we will retain it for now . the leading contribution from both terms in the @xmath304 current
is @xmath305 + \frac{g^2}{\pi}\frac{\phi(z , p ) } { p_{1-}^{(1)}(p_{1-}^{(1)}-q_-)}\,\left[{\bar u}(p_1)\,\gamma_-\ , v(p_2)\right ]   \right\}\nonumber\\   & & + \frac{e_2}{\cal n}\left\{m_{02 } \frac{\phi(1-z',p)}{p_{2-}^{(2)}-q_-}\,\left[{\bar u}(p_1)\ , \gamma_-\,\gamma_+\ , v(p_2 ) \right ] + \frac{g^2}{\pi } \frac{\phi(z , p ) } { p_{2-}^{(2)}(p_{2-}^{(2)}-q_-)}\ , \left[{\bar u}(p_1)\,\gamma_-\ , v(p_2)\right ] \right\ }   .\qquad\quad \label{5eq1e}\end{aligned}\ ] ] the @xmath273 component of the hadron tensor is then @xmath306 using ( [ 5eq1d ] ) and ( [ 5eq1e ] ) , and summing over colors , the @xmath307 contribution to the current sum becomes , at large @xmath259 , @xmath308\left(m_{01}^2-\frac{g^2}{\pi}\right ) \nonumber\\ & = & 8\pi\,e_1 ^ 2\,m_{1}^2 \left(\frac{1-x}{x}\right ) \;\phi^2(x , p)\ , .
\label{5eq1g}\end{aligned}\ ] ] hence @xmath309 and this term is subleading at high @xmath259 .
note that the effect of neglecting the term proportional to @xmath303 in the @xmath40 component of the quark current ( as we are instructed to do in the partonic picture ) is to replace @xmath310 by the bare quark mass @xmath311 .
if @xmath312 , then @xmath313 scales as @xmath314    up to a factor , this is the result originally obtained by einhorn @xcite . in that paper ,
the dirac space was restricted to 2 dimensions , @xmath3 was necessarily zero , and @xmath277 scaled as ( [ einhorn ] ) .
this distribution amplitude does not have the correct physical interpretation .
in particular , it predicts that the dis cross section depends on the square of bare quark mass , a clearly unphysical result . only by considering the full four dimensional dirac space
are we able to obtain scaling for @xmath276 , and the familiar physical result ( [ 5eq2c ] ) for @xmath315 .
another consequence of the fact that @xmath273 is subleading is that we can avoid coming to grips with the fact that the current components @xmath316 and @xmath304 , as defined above , do not conserve current .
since these components do not contribute to the final result , we need not discuss how they can be redefined in order to conserve current .
the _ exact _ hadronic tensor ( [ wmn ] ) for electroproduction of a _ single _ final state @xmath317 with mass @xmath318 ( denoted by @xmath115 ) from the ground state @xmath9 with mass @xmath319 ( denoted by 0 ) is @xmath320 this is a delta function . in any physical measurement
, the detectors will accept a finite range of values of the final electron energy , @xmath321 , requiring that this theoretical cross section be averaged over @xmath321 .
this averaging process will be discussed in subsection [ sec : duality ] below , where the deta function will be written as a delta function in the bjorken variable @xmath136 , which spikes at values of @xmath322 corresponding to the excitation of the final state @xmath115 . in the c.m .
system , this leads to the following _ exact _ transformation @xmath323 with @xmath322 the value of @xmath136 at the bound state mass @xmath318 , given in eq .
( [ xf ] ) below .    in the dis limit ,
the transverse currents dominate the scattering .
the exact result for the transverse tensor , @xmath296 , that describes the scattering from the ground state with @xmath126 structure , to a _ single _
transverse final state @xmath115 with @xmath139 structure , is @xmath324 where the reduced transverse current was defined in eq.([tranff ] ) .
( here the @xmath325 subscript is replaced by @xmath115 , and denotes a transition between the ground and @xmath115th state . )
next we observe that the value of @xmath136 at which the final state is excited is @xmath326 and that the spacing between neighboring levels is @xmath327 with @xmath328 .
hence the _ exact _ result for the total inelastic tensor @xmath296 is @xmath329 where the sum is over all possible final states @xmath115 .
the dis limit is defined by the requirement that @xmath301 with @xmath136 held constant . in this limit
both @xmath259 and @xmath153 become large , with @xmath330 study of the dis limit therefore requires estimating the transition form factors for large final state number @xmath115 , but for all @xmath136 ( since @xmath322 depends on @xmath259 ) . in the dis limit ,
@xmath331 [ easily seen from the expansions ( [ 5eq5 ] ) ] , so @xmath137 will be replaced by @xmath136 in the following discussion .
when @xmath115 is large the `` tails '' of the wave functions are very small , and the exact transition current ( [ tranff ] ) can be approximated by neglecting the contribution from the region @xmath332 $ ] ( where @xmath333 ) , giving @xmath334,p ) } { \xi'[\xi'(1-x)+x]}\ ,   f_{1}^{\rm eff}(\eta , q^2 ) \ , , \qquad\label{translim}\end{aligned}\ ] ] where we continue to suppress the label @xmath282 on the ground ( initial ) state wave function , and @xmath335 in the region of integration .
when @xmath336 and @xmath301 , the relation ( [ ff3 ] ) may be used to approximate the quark form factor by @xmath337 giving @xmath338 equation  ( [ translim ] ) is now further reduced by expanding @xmath339 around @xmath340 @xmath341 where the higher order terms are proportional to powers of @xmath342 @xmath343 this replacement gives the following series for @xmath344 @xmath345 where @xmath346 using the identity ( [ idet4 ] ) and the other estimates given in appendix [ appen : a2 ] it follows that @xmath347 hence the large @xmath115 approximation to ( [ translimseries ] ) is @xmath348    the plus component of the transition current is smaller than the transverse current by a factor of @xmath349 [ compare the two results in eq .
( [ currentd ] ) ] , and hence is obtained immediately from the result ( [ translim2 ] ) . using eq .
( [ 5eq5 ] ) for @xmath350 , and eq .  ( [ xf2 ] ) to replace @xmath351 , gives @xmath352 the minus component can also be obtained from eq .
( [ translim2 ] ) .
reviewing the derivation , we see that the minus component , in the dis limit , is obtained from the transverse component by replacing @xmath353 hence , multiplying ( [ translim2 ] ) by @xmath354 @xmath355 where the @xmath356 integral was evaluated using identity ( [ b10 ] ) .
the transition currents can be shown to be gauge invariant by useing the ward - takahashi identity on the quark - photon vertex , and reducing the result using the bethe - salpeter equation .
these currents will now be used to determine the structure functions @xmath3 and @xmath262 in the hadron picture .      the cross section in the dis limit
can now be obtained from the exact formula ( [ tensorxx2 ] ) . as long as @xmath136 is not too close to unity
each term in the sum must correspond to some state with large @xmath115 ( the states with small @xmath115 `` pile up '' near @xmath246 where their total contribution is very small ) , and the approximate result for the transverse current , ( [ translim2 ] ) , may be used .
furthermore , for large @xmath115 , @xmath357 may be approximated by using ( [ b11 ] ) , and @xmath358 .
hence the final hadronic result [ denoted @xmath359 to distinguish it from the partonic result @xmath317 of eq .
( [ 5eq2c ] ) ] in the dis limit is @xmath360 note that @xmath361 is a series of separate spikes , zero for all @xmath136 except at particular values @xmath322 , where it is infinite .
the partonic function , @xmath269 , is a smooth function of @xmath136 . to compare @xmath361 to @xmath269
, we must average over @xmath136 .
this will lead to the concept of duality and will be discussed in the next subsection below .
now look at the @xmath273 component of the tensor . replacing the transverse components of the current in eq .
( [ tensorxx2 ] ) with plus and minus components , and taking the dis limit , gives @xmath362 substituting eqs .
( [ translimplus ] ) and ( [ translimminus ] ) , and droping the term proportional to @xmath363 ( which averages to zero ) gives @xmath364 this is a nonleading term and can be neglected in the dis limit .
hence the hadronic picture also gives the callan - gross relation , and it is sufficient to compare the functions @xmath233 and @xmath365 only .    if @xmath3 were zero , eq .
( [ wpm2 ] ) would give the following scaling relation for @xmath313 @xmath366 this is to be compared with ( [ einhorn ] ) .
hence , the duality of @xmath3 and the duality of the nonleading terms in @xmath313 depends on comparison of the same functions , @xmath315 and @xmath367 .
we discuss this comparison now .
in any physical measurement , the detectors which define the final state will accept a range of final electron energies @xmath368 .
for fixed @xmath259 , this can be converted into the acceptance of a range of values of @xmath136 centered at @xmath369 with width @xmath370 so that @xmath136 varies over the interval bounded by @xmath371 .
then the _ experimentally measured _
hadronic structure function in the dis limit can be computed from ( [ f1 ] ) @xmath372 note that this is the sum of smooth terms proportional to @xmath307 and @xmath373 , and a rapidly oscillating interference term proportional to @xmath374 . in the dis limit ,
the separation @xmath375 between states approaches zero , and therefore a large number of states are necessarily included in any interval @xmath370 .
the interference term , which changes sign as @xmath363 , therefore averages to zero .
this cancellation of the interference term was pointed out originally by einhorn @xcite , and empahsized recently by close and isgur in the context of the nonrelativistic quark model @xcite .
the smooth terms can be approximated by their value at the center of the interval .
using the fact that the number of states in the interval times the @xmath136 spacing between them must necessarily equal the width @xmath370 , so that @xmath376 we can reduce ( [ f1hadronb ] ) to @xmath377 where @xmath315 was defined in eq .
( [ 5eq2c ] ) .
this result should be compared with the partonic result , averaged in the same way . since
the partonic result is smooth ( and the interference term already neglected ) the partonic average is trival @xmath378 hence we see that @xmath379 proving duality .
note that a similar identity holds for the @xmath115th moments of @xmath317 . in the high @xmath259 limit when the final states are numerous and close enough together so that @xmath380 , the @xmath81th moment of the hadronic structure function is @xmath381 \nonumber\\ & \to & \int_0 ^ 1 dx\ ; x^{m}\;\left[e^2_1\,\phi^2(x , p ) + e^2_2\,\phi^2(1-x\ , , p ) \right]\nonumber\\ & = & \int_0 ^ 1dx\;x^m f(x ) \ , .
\label{f1a}\end{aligned}\ ] ] quark - hadron duality means that properly averaged hadronic observables in a certain kinematic regime ( high @xmath259 ) can be described by perturbative qcd , as schematically represented in fig .
[ hadronic ] and demonstrated in eq.([f1hadron ] ) .
recent data from jefferson laboratory exhibit duality , and this has focused attention on this old subject .
recent studies show that the weak decays within the t hooft model exhibit duality @xcite , and in the framework of dis , close and isgur @xcite studied its emergence in the context of realistic nonrelativistic quark models . a recent paper by isgur , jeschonnek , melnitchouk and van orden @xcite studied duality using a relativistic klein - gordon equation with a confining interaction .
this work is complementary to these other studies .
while we have proved that duality emerges in the dis limit , and have corrected the pathologies of earlier work , we have not yet studied how duality emerges as a function of @xmath259 .
a future paper will present numerical studies of the onset of duality for rqcd @xcite .
we have shown that a very satisfactory description of deep inelastic scattering can be obtained from large @xmath1 qcd in 1 + 1 dimensions , _ provided _ the quark spin degrees of freedom are extended to the full 1 + 3 dimensions ( along with the leptons and the electromagnetic field ) .
we refer to this new model as reduced qcd , and suggest that rqcd may be used to provide insight into the behavior of many other physical processes . in the applications to dis discussed in this paper , we are able to confirm the callan - gross relation , and that duality emerges in the deep inelastic limit .
our theoretical discussion teaches several lessons .
electromagnetic gauge invariance and duality emerge only if the quark current is fully and consistently dressed .
the way in which a quark form factor emerges from this dynamical model mimics the phenomenlogical approach of gross and riska @xcite , showing that their prescription has a dynamical justification .
our treatment of confinement removes all singularities , giving a finite mass for the dressed quark .
this mass is gauge dependent , showing that it is not a physical quantity , which can be taken as a demonstration of confinement .
but the fact that we have finte dressed quark masses in the presence of confinement provides justification for the constituent quark model .
the simultaneous presence of confinement and finite quark mass is possible only if the `` on - shell '' quark scattering amplitude is identicallyt zero , and we have shown that this is indeed true , as modeled previously in ref .
@xcite .
the treatment of reduced qcd presented here lays a foundation for much further study .
this work has been supported in part by the doe through grant no .
de - fg05 - 88er40435 . the southeastrn universities research association ( sura ) operates the thomas jefferson national accelerator facility under doe contract de - aco5 - 84er40150 .
one of us ( z. b. ) has been supported by a predoctoral fellowship at the college of william & mary and by the funda~ ao para a ci^ encia e a tecnologia under the post doctoral grant sfrh / bpd/5661/2001 the project grant cfif - plurianual .
integrals over the variable @xmath69 like the one encountered in eq .
( [ bse2 ] ) appear several places in this paper . here
we evaluate it in full generality .
consider the integral @xmath382 \left[m_2 ^ 2-(k - r)^2 -i\epsilon\right]}\ , .
\label{int1}\end{aligned}\ ] ] substituting for @xmath179 gives @xmath383 where the denominator @xmath183 is @xmath384 \left[\frac{m_2 ^ 2}{(k - r)_-}-2(k - r)_+ -i\,\frac{\epsilon}{(k - r)_-}\right]\nonumber\\ & = & \left[\frac{m_1 ^ 2}{y\,r_-}-2k_+ -i\,\frac{\epsilon}{y\,r_-}\right ] \left[-\frac{m_2 ^ 2}{(1-y)\,r_-}-2(k - r)_+ + i\,\frac{\epsilon}{(1-y)\,r_-}\right]\ , , \qquad \label{d2}\end{aligned}\ ] ] and we have introduced the momentum fractions ( [ defz ] ) .
this shows that the denominator has only two poles in @xmath69 , and that they will both be in the same half of the complex plane ( giving zero for the integral ) unless @xmath385 doing the integral over @xmath69 , and expressing the answer in terms of the momentum fractions @xmath137 and @xmath32 gives @xmath386 where the new denominator @xmath387 is @xmath388 \ , .\label{d1}\ ] ] [ note that this result holds for _ both _ signs of @xmath389 .
this is because the sign of the countor integral , which depends on the sign of @xmath389 , is cancelled by the sign of the integral over @xmath390 , which must be changed from @xmath391 when @xmath392 . ]
it is important to note that the integral over @xmath137 is not restricted in the term in ( [ 2eqaa ] ) that is proportional to @xmath393 because for that term the restriction ( [ r1 ] ) applies to @xmath32 and not @xmath137 .
this term not only regulates the singularity in the first term ( which arises only if @xmath394 ) , but also contibutes an additional contribution .
furthermore , if @xmath32 lies outside of the interval @xmath394 , the first term is non singular ( except at the end points of the interval ) and is nonzero !
hence , in the region @xmath394 the integral ( [ 2eqaa ] ) becomes @xmath395\right)\ , , \qquad\quad \label{2eqab}\end{aligned}\ ] ] and outside of this region it is @xmath396\right)\ , , \qquad \label{2eqac}\end{aligned}\ ] ] where we have defined @xmath397$\cr { \rm v}_0(z , y)\equiv { \displaystyle\frac{g^2}{\pi } \frac{1}{(z - y)^2 } } & if $ z\notin[0,1]$ } \ , .\qquad \label{eq251a}\end{aligned}\ ] ] the multiplication by @xmath398 in ( [ eq251a ] ) is needed to scale the momenta from @xmath399 , etc .
we emphasize that the identities ( [ 2eqab ] ) and ( [ 2eqac ] ) hold for both @xmath400 and @xmath392 .
the first case is needed for the bound state equation and the second for electron scattering .
in this appendix we derive a number of identities that hold when @xmath123 ( i.e. in all cases but the chiral limit ) . equation ( [ 2eqee ] )
leads immediately to the following identity @xmath401 integrating by parts gives @xmath402 an additional identity is derived by multiplying ( [ 2eqee ] ) by @xmath403 and integrating .
first observe that @xmath404 hence , integrating by parts and using the normalization condition ( [ norm ] ) for the wave function , we obtain @xmath405 similarly we can prove that @xmath406 additional identities are needed for the study of duality .
first , integrate eq .
( [ phiin ] ) over @xmath32 to obtain @xmath407 next , the commutator of the operator @xmath408 with the hamiltonian defined in ( [ phiin ] ) is @xmath409\,\phi_n(x , r)=\int_0 ^ 1dy \left [ \frac{m^2_{01}}{xy } -\frac{m^2_{02}}{(1-x)(1-y)}\right]\phi_n(y , r ) \ , , \label{com1}\end{aligned}\ ] ] as discussed in refs .
@xcite . multiplying this by @xmath410 and integtating over @xmath136 gives zero on the l.h.s . , and
hence the relation @xmath411 taking the square root of both sides of this equation gives @xmath412 where the sign is positive if the phase of @xmath413 is the same as @xmath414 and negative if it is opposite .
since the states are non - degenerate , the phase of the @xmath115th eigenstate is @xmath415 . combining this with eq .
( [ b6 ] ) gives @xmath416 c_n\ , , \label{b10}\end{aligned}\ ] ] where the constant @xmath417 is yet to be determined . in this paper
we have chosen the sign of @xmath112 so that it is always positive as @xmath418 , so that @xmath419 for all @xmath115 .
our demonstration of duality depends upon the relations @xmath420 which can be derived approximately using arguments given by t hooft @xcite , and callan , coote , and gross @xcite , given here for completeness .
first , for large @xmath115 , the boundary conditions require that the normalized wave function go like @xmath421 one test of this approximation is shown in fig .
[ fig : idet1 ] .
if @xmath115 is very large , so that @xmath422 , the behavior of the wave function at the endpoints can be ignored , and near @xmath423 eq .  ( [ 2eqee ] ) can be approximated @xmath424 substituting the ansatz ( [ ansatz ] ) into this equation gives @xmath425 hence @xmath426 as @xmath427 and identity ( [ b11 ] ) is proved . to prove ( [ identb ] ) we consider the matrix element of ( [ com1 ] ) @xmath428\,\phi_n(x , r)&= & 2c_{n+1}c_n   \nonumber\\ & = & \left(m_{n+1}^2-m_n^2\right)\dashint_0 ^ 1 dxdy \frac{\phi_{n+1}(x , r)\phi_n(y , r)}{y - x}\ , .\end{aligned}\ ] ] where we used the fact that the phases of @xmath429 and @xmath430 are opposite .
the integral is evaluated in the large @xmath115 approximation @xmath431\pi x)\sin(n\pi y)}{y - x - i\epsilon } \nonumber\\ & = & 2 % \frac{2\pi}{n_c } \int_0 ^ 1 dx \sin([n+1]\pi x)\cos(n\pi x)\nonumber\\ & = & 2 % \frac{2}{n_c } \left(1 + \frac{1}{2n+1}\right )    \to 2 % \frac{2}{n_c }    \ , .\end{aligned}\ ] ] hence , @xmath432 and ( [ identb ] ) is proved .
a numerical demonstration of the identities ( [ idet4 ] ) , ( [ identb ] ) , and ( [ b11 ] ) is given in figs.[fig : idet1 ] .
the identities are satisfied to better than a few percent for state numbers between 40 and 80 , provided the number of splines is at least 160 ( twice the number of the maximum state of interest ) .    two additional relations that follow from the completeness relation ( [ complete ] ) and the bound state equation are needed . if @xmath433 $ ] , then : @xmath434 the proof of the second relation ( [ id2 ] ) follows from application of the bound state equation @xmath435\int_0 ^ 1dz'\ , \phi_n(z',r_n)\nonumber\\ & & \qquad=\sum_n\left\{\int_0^{z-\epsilon } + \int_{z-\epsilon}^{z+\epsilon}+\int_{z+\epsilon}^1\right\ } dy\ , \left[\frac{\phi_n(y , r_n)-\phi_n(z , r_n ) } { ( y - z)^2 } \right]\int_0 ^ 1dz'\ , \phi_n(z',r_n)\nonumber\\ & & \qquad\to\sum_n\frac{1}{2}\int_{z-\epsilon}^{z+\epsilon}dy\ ; \phi_n''(z , r_n)\int_0 ^ 1dz'\,\phi_n(z',r_n ) = \nonumber\\ & & \qquad= \epsilon\,\sum_n \phi_n''(z , r_n)\int_0 ^ 1dz'\,\phi_n(z',r_n)\to0 \ , , \end{aligned}\ ] ] where the identity ( [ id1 ] ) is used to eliminate the nonsingular integrals over the intervals @xmath436 $ ] and @xmath437 $ ] , and the singular part of the integral is proportional to @xmath438 which vanishes as @xmath439 .    if eq .
( [ phiout ] ) is used , modifications of the completeness relation ( [ complete ] ) and the identity ( [ id1 ] ) , which hold in the region @xmath440 $ ] , can be derived .
the results are @xmath441 \label{complete2}\\ & & % n_c
\sum_n\left(\frac{m_1 ^ 2}{z}+\frac{m^2_2}{1-z}- m_n^2\right)\phi_n(z , r_n)\int_0 ^ 1dz'\ , \phi_n(z',r_n)= -\frac{g^2/\pi}{z(1-z ) } \label{id3 }   \ , .\end{aligned}\ ] ] note that the completeness relation and all of these identities hold for _ any _ two quark hamiltonian . in particular , the hamiltonian may describe bound states with flavor ( where , in general , @xmath442 ) as well as flavorless bound states where @xmath443 . in some applications we will study the coupling of a photon to a @xmath4 pair .
such a coupling necessarily involves _ flavorless _ states only , and we will use the identities for the flavorless states .
here we show how to sum this series ( [ eq249 ] ) and obtain a useful form for @xmath184 .
since @xmath179 is independent of the plus components of momentum , the integrations over the plus components of momentum can be carried out . using the identity ( [ 2eqab ] ) [ with @xmath198 in this case ]
, the first term on the rhs of ( [ eq249 ] ) becomes @xmath444\ , , \quad\quad\quad\label{first}\end{aligned}\ ] ] or @xmath445\ , .\quad\quad\quad\label{firsta}\end{aligned}\ ] ] here the replacement ( [ eq251a ] ) has been used to convert @xmath446 ( or v@xmath447 ) .
similarly , the second term becomes @xmath448 \ , , \quad\quad\quad \label{second}\end{aligned}\ ] ] or @xmath449 \ , , \quad\quad\quad \label{seconda}\end{aligned}\ ] ] and the third term follows the same pattern . replacing @xmath189 and @xmath190 , the series ( [ eq249 ] ) , in the interval @xmath450 $ ] ,
is clearly summed by the following equation @xmath451 \
, .\qquad\qquad \label{eq253}\end{aligned}\ ] ] introducing the shorthand notation @xmath452 eq .  ( [ eq253 ] ) can be rearranged into the following form @xmath453 which requires that @xmath454 satisfy the equation @xmath455 this equation is easily solved with the use of the greens function : @xmath456 putting it all together gives eq .
( [ eq267 ] ) .
if @xmath457 $ ] , the series reduces to the same formal result @xmath458 where the completeness relation ( [ complete ] ) and eq .
( [ phiout ] ) were used in the last step .
this shows that ( [ eq267 ] ) can be applied for _ any _ @xmath459 $ ] .
the series ( [ current1a ] ) for the dressed quark current is evaluated in this appendix .
the `` first '' integral over @xmath69 in each term in the series is evaluated as in eq .
( [ first ] ) or ( [ firsta ] ) @xmath460$\cr   { \displaystyle\int_0 ^ 1dy\,{\rm v}_0(z',y)\,\frac{n_i^\mu(k_-,q_-)}{\delta_i(y , q^2 ) } } & if $ z'\notin[0,1]\ , , \qquad\quad$ } \label{current1b}\end{aligned}\ ] ] where @xmath461\,\gamma^\mu\,\left[m_{0i } + ( k_--q_-)\gamma_+\right]\,\gamma_-}{4\,k_-(k_--q_-)}\ , , \label{current1bb}\end{aligned}\ ] ] and @xmath462 and @xmath463 are ( [ d2 ] ) and ( [ d1 ] ) with equal masses @xmath464 .
note that the only dependence of the integrand on @xmath69 was in the denominator @xmath462 , so the @xmath69 integral could be evaluated using the methods of appendix [ appen : b ] .
the second integral is similarily evaluated .
if @xmath450 $ ] , the following result is obtained for the series eq .
( [ current1a ] ) @xmath465 where @xmath466 , and the series is summed using eq .
( [ eq253 ] ) for the scattering matrix m. using the solution ( [ eq267 ] ) for the scattering matrix the result , for all values of @xmath215 , is @xmath467 where @xmath468 , and @xmath463 is the @xmath469 of eq .
( [ d1 ] ) with equal masses @xmath464 .
the same expression also holds for @xmath457 $ ] , but the derivation differs . as before , the full result can be deduced from the form of the first two terms of the series ( [ current1a ] ) . guided by ( [ current1b ] ) , and using the steps already displayed in eq .
( [ eq267a ] ) , we obtain @xmath470 we see , as we did with eq .
( [ eq267 ] ) , that ( [ current2a ] ) holds for _ all _ @xmath289 .    equation ( [ current2a ] ) expresses the dressed current as a sum over contributions from the bound states that couple to the photon , showing that vector dominance is a rigorous consequence of this model .
the structure of the numerator n@xmath471 depends on the bare current : @xmath472 note the interesting structure of the correction term to the bare current @xmath141 , which will be of central importance in our discussion below . substituting these results into ( [ current2a ] ) gives the following    @xmath473    we organize these expressions by introducing the quark form factor , defined in eq .
( [ ff1 ] ) .
this immediately gives eq.([ff3aa ] ) for the @xmath212 current .
the other currents in ( [ current5 ] ) can also be written in terms of @xmath219 . to this end note
that the bound state equation for equal mass quarks implies that @xmath474 hence the quark mass can be removed from the second term in @xmath475 , giving @xmath476 where @xmath477 defined by this equation should not be confused with the two body greens function or the bound state vertex function .
if @xmath450 $ ] , the second term can be reduced using the identities ( [ id1 ] ) and ( [ id2 ] ) @xmath478 \phi_n(z',r_n)\ ; \int_0 ^ 1dy\,\phi_n(y , r_n ) = -\frac{g^2}{\pi}\frac{1}{z'(1-z ' ) } \ , .\qquad\quad \label{ff3}\end{aligned}\ ] ] the same result is obtained directly from identity ( [ id3 ] ) if @xmath457 $ ] .
hence , for all @xmath289 , @xmath479    finally , the @xmath480 component of the current can be similarly reduced .
the second term in @xmath480 is proportional to @xmath481 the part of the second term in @xmath482 is the same constant derived above , and combining the results gives eq .
( [ ff6 ] ) .
kinematics for dis in the c.m .
frame are given in this appendix .
when @xmath259 is very large , the components of @xmath224 and @xmath174 in this frame are ( see ref .
@xcite ) @xmath483      \nonumber\\ q_-=&&-q\sqrt{\frac{x}{2(1-x ) } } \left[1+\frac{m_0 ^ 2 x ( 1 - 2x)}{2q^2(1-x)}+ ... \right ]   \nonumber\\   p_+=&&\frac{m_0 ^ 2}{q}\sqrt{\frac{x(1-x)}{2 } } \left[1-\frac{m_0 ^ 2x(1 - 2x^2)}{2q^2(1-x)}+ ... \right ]
\nonumber\\   p_-=&&\frac{q}{\sqrt{2x(1-x ) } } \left[1+\frac{m_0 ^ 2x(1 - 2x^2)}{2q^2(1-x)}+ ... \right]\ , .
\label{5eq5}\end{aligned}\ ] ] the energy and momentum of particles 1 and 2 in the final state are @xmath484 \nonumber\\ p_{10}&=&\sqrt{m_1 ^ 2+p_z^2}= \frac{q}{2}\sqrt{\frac{1-x}{x } } \left[1+\frac{x(m_0 ^ 2 + 2m_1 ^ 2 - 2m_2 ^ 2)}{2q^2(1-x)}+\cdots\right ] \nonumber\\
p_{20}&=&\sqrt{m_2 ^ 2+p_z^2}= \frac{q}{2}\sqrt{\frac{1-x}{x } } \left[1+\frac{x(m_0 ^ 2 - 2m_1 ^ 2 + 2m_2 ^ 2)}{2q^2(1-x)}+\cdots\right ] \ , .\end{aligned}\ ] ] there are two possibilities corresponding to the two terms in the current ( [ 5eq1b ] )
. the first term in the current will be large only if the momentum of particle 1 is in the direction of @xmath485 ( i.e.  @xmath486 ) .
[ choosing it in the opposite direction requires a very large momentum flow through the wave function , and is suppressed . ] in this case @xmath487
\nonumber\\ p^{(1)}_{1-}&=&\sqrt{\frac{x}{2(1-x)}}\ ; \frac{m_1 ^ 2}{q } \left[1-\frac{x(m_0 ^ 2 - 2m_2 ^ 2)}{2q^2(1-x)}+\cdots\right ] \nonumber\\ p^{(1)}_{2+}&=&\sqrt{\frac{x}{2(1-x)}}\ ; \frac{m_2 ^ 2}{q }   \left[1-\frac{x(m_0 ^ 2 - 2m_1 ^ 2)}{2q^2(1-x)}+\cdots\right ] \nonumber\\ p^{(1)}_{2-}&=&q\sqrt{\frac{1-x}{2x } } \left[1+\frac{x(m_0 ^ 2 - 2m_1 ^ 2)}{2q^2(1-x)}+\cdots\right ] \ , , \label{5eq5a}\end{aligned}\ ] ] where the superscript @xmath488 is a reminder that these relations hold only for the _ first _ term in the current .
the momentum relations for the _ second _ term in the current ( [ 5eq1b ] ) follow by interchanging 1 and 2 on both sides of these equations .
we see that the first term in the current gives @xmath489 and @xmath490 , while the second term gives @xmath491 and @xmath492 .
the two terms describe kinematically distinct regions of phase space , and their interference can safely be neglected when the total cross section is computed . using the definitions ( [ defza ] ) and ( [ defzp ] ) , @xmath32 and @xmath289 can be easily related to the bjorken variable @xmath136 .
we have @xmath493
this appendix includes details of the evaluation of the transition currents shown in fig .
[ hadronicamplitude ] . the @xmath88 term will be evaluated first .
the @xmath90 term can then be obtained by a simple substitution , as discussed at the end of this appendix .
the @xmath88 term , given in eq .
( [ 3eq1- ] ) , is reduced by introducing the momentum fractions @xmath494 so that @xmath495 in this notation the @xmath136-type transverse quark current , evaluated in eq .
( [ ff6 ] ) , is @xmath496 first evaluate the numerators , @xmath228 , of the traces in ( [ 3eq1- ] ) .
the numerator @xmath497 is @xmath498\nonumber\\ & = & 16\,k_- \,\left[m_{01}\,q_- + 2(p_f+k)_-(p+k)_- h_1(\eta , q)\right]\nonumber\\ & = & 16\,k_-q_- \frac{m_1 ^ 2}{m_{01}}-32(p_f+k)_-(p+k)_-\frac{q_+}{m_{01}}\ , f_1(\eta , q^2)\nonumber\\ & = & \frac{16k_-q_-}{m_{01}}\left(m_1 ^ 2- \eta(1-\eta)\,q^2 \,f_1(\eta , q^2)\right)\ , , \label{nx}\end{aligned}\ ] ] where the trace has been evaluated in a dirac space of four dimensions .
the numerator of the minus component of the current has only one term @xmath499\,\left[1+f_1(\eta , q^2)\right]\nonumber\\ & = & 32\,k_- ( p_f+k)_-(p+k)_-\,\left[1+f_1(\eta , q^2)\right ] \nonumber\\ & = & 16k_-\,\frac{q_-}{q_+}\,\eta(1-\eta)\,q^2\ , \left[1+f_1(\eta , q^2)\right]\ ,   .
\label{minus}\end{aligned}\ ] ] the trace for the plus component is @xmath500\nonumber\\ & = & 16k_-m_1 ^ 2 - 32k_-(p_f+k)_-(p+k)_-\,\frac{q_+}{q_-}\ , f_1(\eta , q^2 ) \nonumber\\ & = & 16k_-\left(m_1 ^ 2-\eta(1-\eta)\,q^2\ , f_1(\eta , q^2)\right)\ , .\label{plusc}\end{aligned}\ ] ]    inserting the numerator ( [ minus ] ) into the general result ( [ 3eq1- ] ) gives the following result for the minus component of the current @xmath501[m_1 ^ 2-(p+k)^2-i\epsilon ] [ m_2 ^ 2-k^2-i\epsilon]}\,\left[1+f_1(\eta , q^2)\right ]    \nonumber\\   & & \!\!\!\!=32i e_1\int \frac{d^2k}{(2 \pi)^2 }   \;\frac{1}{d_3 } g_f\left(-p_f - k ,-
k\right)\,g_i\left(p+k , k\right ) \,\left[1+f_1(\eta , q^2)\right ] \ , , \label{3eq1}\end{aligned}\ ] ] where the subscript - on @xmath232 reminds us that the final state must have a pure @xmath126 structure (
if the ground state does ) , @xmath502 is the product of the denominators of the three quark propagators , with three poles in @xmath69 @xmath503 \left[\frac{m_1 ^ 2}{\xi\,p_-}-\frac{m_i^2}{p_-}-2k_+-i \epsilon_b\right ] % \nonumber\\&&\times \left[\frac{m_2 ^ 2}{(\xi-1)\,p_-}-2k_+-i\epsilon_c\right]\ , , \qquad \label{3poles}\end{aligned}\ ] ] and the @xmath438 s change sign according to @xmath504\nonumber\\ \epsilon_b&=&\epsilon/(p+k)_-=\epsilon/(\xi\,p_- ) \nonumber\\ \epsilon_c&=&\epsilon / k_-=\epsilon/[(\xi-1)\,p_- ] \ , .
\label{epsilon}\end{aligned}\ ] ] since the vertex functions @xmath128 do not depend on @xmath69 , we can evaluate the @xmath69 integral .
it will be nonzero only when the three poles of ( [ 3poles ] ) do not all lie in the same half plane . since @xmath198 for electron scattering , @xmath505 , which implies that @xmath506 hence @xmath507 and all three poles will be in the same half of the complex plane ( giving zero for the integral ) unless @xmath508 there are two terms , depending on the sign of @xmath509 . closing the @xmath69 contour in the upper half plane , and using @xmath510 and @xmath511 gives @xmath512 \nonumber\\      & & -16 e_1\,p_-^2p_{f- } \int_{0}^{y}\frac{d\xi}{(2\pi ) } \;\frac{g_f\left(\xi',p_f\right)\,g_i\left(\xi , p\right ) \left[1+f_1(\eta , q^2)\right ] } { \delta(\xi',p_f^2)\,\left(\delta(\xi , p^2)-\displaystyle { \frac{\delta(\xi',p_f^2)}{(1-y ) } } \right ) } \ , , \qquad\quad \label{3eq1a}\end{aligned}\ ] ] where the first term is the contribution from the pole due to the zero in the third term in ( [ 3poles ] ) above , and the second from the pole due to the zero in the first term .
using the definition ( [ 2eq3d ] ) of the wave function and combining the two contributions in the region @xmath513 $ ] we get eq .
( [ 3eq1b ] ) with @xmath514 } { ( 1-\xi')\left\{\xi\,\xi'\left(m_f^2-(1-y)m_i^2 \right)-y\ , m_1 ^ 2\right\ } } \quad\to 1\quad{\rm if}\quad m_f^2\to \infty\ , .
\label{rlimit}\end{aligned}\ ] ] this factor is needed for positive values of @xmath515 in the interval @xmath513 $ ] and _ negative _ values of @xmath356 in the interval @xmath516 $ ] . in this region
the denominator has a zero only if @xmath517 . in order to avoid the discussion of such cases we limit numerical applications to cases with @xmath518 .
results for the other components of the transition current , eq.([currentd ] ) , are obtained by a similar arguement using ( [ nx ] ) and ( [ plusc ] ) in place of ( [ minus ] ) .
now consider the modifications required in order to evaluate the @xmath90 term .
using the momenta defined in fig.[hadronicamplitude ] , and the same definitions of momentum fractions ( [ frac2 ] ) and ( [ frac3 ] ) , the @xmath90 term is obtained from the @xmath88 term simply by substituting @xmath519 and @xmath520 . however , the momentum fractions in the wave functions are , by convention , the fraction of the momentum carried by the quark @xmath87 , and the momentum fraction in the quark form factor is that of the the _ outgoing _ quark , and hence these functions must be written in terms of @xmath521 these observations lead immediately to the final result ( [ exacttran ] ) .
the two body equations were solved using a modification of the standard cubic splines employed previously in many problems .
the standard splines are defined on 4 segments of length @xmath522 , bounded by the 5 points @xmath523 : @xmath524 &     if $ b < x < c$\cr   1 + 3\left[\frac{(d - x)}{h}+\frac{(d - x)^2}{h^2}-\frac{(d - x)^3}{h^3}\right ]   &      if $ c < x < d$\cr \frac{(e - x)^3}{h^3 }    &    if $ d < x
< e$\ , .}\end{aligned}\ ] ] if the interval [ 0,1 ] is spanned by @xmath525 standard splines , there must be @xmath164 segments , of length @xmath526 .
the first spline is numbered @xmath527 beginning at @xmath528 and the last is numbered @xmath529 ending at @xmath246 .
the modified splines used in this paper consist of the standard splines plus one additional , non - standard spline inserted at the beginning and end of the interval [ 0,1 ] .
these non - standard splines are defined on only three segments , and will be numbered 1 and @xmath155 . they are @xmath530   &      if $ h < x<2h$\cr \frac{(3h - x)^3}{4h^3 }    &    if $ 2h < x<3h$ } \nonumber\\ s_{n_s}(x)&=&\cases{\frac{(x-1 + 3h)^3}{4h^3 }     &    if $ 1 - 3h < x<1 - 2h$\cr   \frac{1}{4 }   + \frac{3}{4}\left[\frac{(1-h - x)}{h}+\frac{(1-h - x)^2}{h^2 }   -\frac{(1-h - x)^3}{h^3}\right ]   &      if $ 1 - 2h < x<1-h$\cr d_1\frac{(1-x)^{\beta_2}}{h^{\beta_2}}+d_2\frac{(1-x)^2}{h^2}+ d_3\frac{(1-x)^3}{h^3 }    &    if
$ 1-h < x<1$\ , . } \end{aligned}\ ] ] where the @xmath531 are the fractional exponents determined by the boundary conditions [ given in eq .
( [ betas ] ) ] , and the coefficients @xmath532 depend on the exponents @xmath533 @xmath534 where @xmath535 for @xmath97 or @xmath536 for @xmath537 .
these coefficients were fixed by the requirement that the spline and its first two derivatives be continuous .
figure [ fig : spline ] shows two examples of the modified splines for the minimal number @xmath538 . in this case the central spline spans the full interval [ 0,1 ] and has 4 segments .
99 g.  t hooft , nucl .
* b75 * , 461 ( 1974 ) .
c.  g.  callan , n.  coote and d.  j.  gross , phys .
d * 13 * , 1649 ( 1976 ) .
m.  einhorn , phys .
d * 14 * , 3451 ( 1976 ) .
m. i bars , m. b. green , phys .
rev * d * 15 , ( 1978 ) . c.  g.  callan and d.  j.  gross , phys .
* 22 * , 156 ( 1969 ) .
m.  burkardt , phys .
d * 62 * , 094003 ( 2000 ) [ arxiv : hep - ph/0005209 ] . m.  burkardt and n.  uraltsev , phys .
d * 63 * , 014004 ( 2001 ) [ arxiv : hep - ph/0005278 ] .
v.  schon and m.  thies , contribution to the festschrift in honor of boris ioffe , m.  shifman , ( ed ) . in _ at the frontier of particle physics , vol . 3
_ p.  1945
( 2000 ) [ arxiv : hep - th/0008175 ] .
m.  burkardt , f.  lenz and m.  thies , phys .  rev .
d * 65 * , 125002 ( 2002 ) [ arxiv : hep - th/0201235 ] .
y.  s.  kalashnikova and a.  v.  nefediev , phys .
* 45 * , 347 ( 2002 ) [ usp .
nauk * 172 * , 378 ( 2002 ) ] [ arxiv : hep - ph/0111225 ] . b. grinstein , r. f. lebed , phys .  rev .
d * 57*:1366 ( 1998 ) [ hep - ph/9708396 ] .
r. f. lebed , talk given at 3rd international conference in quark confinement and hadron spectrum ( confinement iii ) , newport news , va , 7 - 12 jun 1998 ; published in _ newport news 1998 , quark confinement and the hadron spectrum iii _ , p.161 ( 1998 ) [ hep - ph/9808452 ] .
r. f. lebed , n. g. uraltsev , phys .
d * 62*:094011 ( 2000 ) [ hep - ph/0006346 ] .
z.  batiz and f.  gross , phys .
c * 58 * , 2963 ( 1998 ) .
+ note that the square bracket in the second of eqs .
( 5.8 ) should be @xmath539 $ ] instead of @xmath540 $ ] .
+ z.  batiz and f.  gross , in preparation .
f. gross and j. milana , phys .
d * 43*(1991 ) 2401 .
+ f. gross and j. milana , phys .
d * 45 * ( 1992 ) 969 . + f. gross and j. milana , phys .
d * 50 * ( 1994 ) 3332 . c.  savkli and f.  gross , phys .
c * 63 * , 035208 ( 2001 ) . f.  gross and d.  o.  riska , phys .
c * 36 * , 1928 ( 1987 ) .
f. e. close , n. isgur , phys
.  lett .
b * 509 * , 81 ( 2001 ) [ hep - ph/0102067 ] .
n. isgur , s. jeschonnek , w. melnitchouk and j. w. van orden , phys .
d * 64*:054005 ( 2001 ) .